Relation between the work function and Young's modulus of RhSi and estimate of Schottky-barrier height at RhSi/Si interface: An ab-initio study

Density-functional calculations are performed to explore the relationship between the work function and Young's modulus of RhSi, and to estimate the p-Schottky-barrier height (SBH) at the Si/RhSi(010) interface. It is shown that the Young's modulus and the workfunction of RhSi satisfy the generic sextic relation, proposed recently for elemental metals. The calculated p-SBH at the Si/RhSi interface is found to differ only by 0.04 eV in opposite limits, viz., no-pinning and strong pinning. We find that the p-SBH is reduced as much as by 0.28 eV due to vacancies at the interface. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4761994]

Dissociation Constants of Weak Acids from ab Initio Molecular Dynamics Using Metadynamics: Influence of the Inductive Effect and Hydrogen Bonding on pK(a) Values

The theoretical estimation of the dissociation constant, or pK(a), of weak acids continues to be a challenging field. Here, we show that ab initio CarParrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free-energy profile of the dissociation reaction provide reasonable estimates of the pK(a) value. Water molecules, sufficient to complete the three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. The free-energy profiles exhibit two distinct minima corresponding to the dissociated and neutral states of the acid, and the difference in their values provides the estimate for pK(a). We show for a series of organic acids that CPMD simulations in conjunction with metadynamics can provide reasonable estimates of pK(a) values. The acids investigated were aliphatic carboxylic acids, chlorine-substituted carboxylic acids, cis- and trans-butenedioic acid, and the isomers of hydroxybenzoic acid. These systems were chosen to highlight that the procedure could correctly account for the influence of the inductive effect as well as hydrogen bonding on pK(a) values of weak organic acids. In both situations, the CPMD metadynamics procedure faithfully reproduces the experimentally observed trend and the magnitudes of the pK(a) values.

Pt-Poisoning-Free Efficient CO Oxidation on Pt3Co Supported on MgO(100): An Ab Initio Study

Late-transition-metal-doped Pt clusters are prevalent in CO oxidation catalysis, as they exhibit better catalytic activity than pure Pt, while reducing the effective cost and poisoning However, completely eliminating the critical problem of Pt poisoning still poses a big challenge. Here, we report for the first time that, among the bimetallic clusters ((Pt3M where M = Co, Ni, and Cu)/MgO(100)), the CO adsorption site inverts for Pt3Co/MgO(100) from Pt to Co, due to the complete uptake of Pt d-states by lattice oxygen. While this resolves the problem of Pt poisoning, good reaction kinetics are predicted through low barriers for Langmuir-Hinshelwood and Mars van Krevelen (MvK) mechanisms of CO oxidation for Pt3Co/MgO(100) and Li-doped MgO(100), respectively. Li doping in MgO(100) compensates for the charge imbalance caused by a spontaneous oxygen vacancy formation. Pt-3 Co/Li-doped MgO(100) stands out as an exceptional CO oxidation catalyst, giving an MvK reaction barrier as low as 0.11 eV. We thereby propose a novel design strategy of d-band center inversion for CO oxidation catalysts with no Pt poisoning and excellent reaction kinetics.

Estimating successive pK(a) values of polyprotic acids from ab initio molecular dynamics using metadynamics: the dissociation of phthalic acid and its isomers

Estimation of the dissociation constant, or pK(a), of weak acids continues to be a central goal in theoretical chemistry. Here we show that ab initio Car-Parrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free energy profile of the dissociation reaction can provide reasonable estimates of the successive pK(a) values of polyprotic acids. We use the distance-dependent coordination number of the protons bound to the hydroxyl oxygen of the carboxylic group as the collective variable to explore the free energy profile of the dissociation process. Water molecules, sufficient to complete three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. Two distinct minima corresponding to the dissociated and un-dissociated states of the acid are observed and the difference in their free energy values provides the estimate for pK(a), the acid dissociation constant. We show that the method predicts the pK(a) value of benzoic acid in good agreement with experiment and then show using phthalic acid (benzene dicarboxylic acid) as a test system that both the first and second pK(a) values as well, as the subtle difference in their values for different isomers can be predicted in reasonable agreement with experimental data.

Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pK(a) Values from Metadynamics Sampling

Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation deprotonation reaction of the 20 canonical alpha amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metad-ynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pK(a) values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pK(a) values with a mean relative error, with respect to experimental results, of 0.2 pK(a) units.

H2O-CH4 and H2S-CH4 complexes: a direct comparison through molecular beam experiments and ab initio calculations

New molecular beam scattering experiments have been performed to measure the total ( elastic plus inelastic) cross sections as a function of the velocity in collisions between water and hydrogen sulfide projectile molecules and the methane target. Measured data have been exploited to characterize the range and strength of the intermolecular interaction in such systems, which are of relevance as they drive the gas phase molecular dynamics and the clathrate formation. Complementary information has been obtained by rotational spectra, recorded for the hydrogen sulfide-methane complex, with a pulsed nozzle Fourier transform microwave spectrometer. Extensive ab initio calculations have been performed to rationalize all the experimental findings. The combination of experimental and theoretical information has established the ground for the understanding of the nature of the interaction and allows for its basic components to be modelled, including charge transfer, in these weakly bound systems. The intermolecular potential for H2S-CH4 is significantly less anisotropic than for H2O-CH4, although both of them have potential minima that can be characterized as `hydrogen bonded'.

Non-covalent C-Cl center dot center dot center dot pi interaction in acetylene-carbon tetrachloride adducts: Matrix isolation infrared and ab initio computational studies

Non-covalent halogen-bonding interactions between n cloud of acetylene (C2H2) and chlorine atom of carbon tetrachloride (CCl4) have been investigated using matrix isolation infrared spectroscopy and quantum chemical computations. The structure and the energies of the 1:1 C2H2-CCl4 adducts were computed at the B3LYP, MP2 and M05-2X levels of theory using 6-311++G(d,p) basis set. The computations indicated two minima for the 1:1 C2H2-CCl4 adducts; with the C-Cl center dot center dot center dot pi adduct being the global minimum, where pi cloud of C2H2 is the electron donor. The second minimum corresponded to a C-H...Cl adduct, in which C2H2 is the proton donor. The interaction energies for the adducts A and B were found to be nearly identical. Experimentally, both C-Cl center dot center dot center dot pi and C-H center dot center dot center dot Cl adducts were generated in Ar and N2 matrixes and characterized using infrared spectroscopy. This is the first report on halogen bonded adduct, stabilized through C-Cl center dot center dot center dot pi interaction being identified at low temperatures using matrix isolation infrared spectroscopy. Atoms in Molecules (AIM) and Natural Bond Orbital (NBO) analyses were performed to support the experimental results. The structures of 2:1 ((C2H2)(2)-CCl4) and 1:2 (C2H2-(CCl4)(2)) multimers and their identification in the low temperature matrixes were also discussed. (C) 2015 Elsevier B.V. All rights reserved.

Ab initio pair potentials at metal-ceramic interfaces

A systematic approach is proposed to obtain the interfacial interatomic potentials. By inverting ab initio adhesive energy curves for the metal-MgO ceramic interfaces, We derive interfacial potentials between Ag and O2-, Ag and Mg2+, Al and O2-, Al and Mg2+. The interfacial potentials, obtained from this method, demonstrate general features of bondings between metal atoms and ceramic ions.

Ab initio investigation of the elasticity and stability of aluminium

On the basis of the pseudopotential plane-wave (PP-PW) method in combination with the local density functional theory (LDFT), complete stress-strain curves for the uniaxial loading and uniaxial deformation along the [001] and [111] directions, and the biaxial proportional extension along [010] and [001] for aluminium are obtained. During the uniaxial loading, certain general behaviours of the energy versus the stretch and the load versus the stretch are confirmed; in each case, there exist three special unstressed structures: f.c.c., b.c.c., and f.c.t. for [001]; f.c.c., s.c., and b.c.c. for [111]. Using stability criteria, we find that all of these states are unstable, and always occur together with shear instability, except the natural f.c.c. structure. A Pain transformation from the stable f.c.c. structure to the stable b.c.c. configuration cannot be obtained by uniaxial compression along any equivalent [001] and [111] direction. The tensile strengths are similar for the two directions. For the higher energy barrier of the [111] direction, the compressive strength is greater than that for the [001] direction. With increase in the ratio of the biaxial proportional extension, the stress and tensile strength increase; however, the critical strain does not change significantly. Our results add to the existing ab initio database for use in fitting and testing interatomic potentials.

MoS2 eta WS2 monogeruzen egitura elektroniko eta magnetikoaren ab-initio kalkuloak

Bi dimentsiotako materialetan presente diren propietate elektroniko bereziek betidanik piztu izan dute komunitate zientifikoaren interesa. Idealki atomo bakarreko lodierako materialak diren hauek hasiera batean joku teoriko huts zirela uste bazen ere, A.K. Geim eta K.S. Novoselov-ek kontrakoa frogatu zuten lehenengo aldiz grafenoa sintetizatuz[1]. Grafitoa osatzen duen geruzetako bakoitza den grafenoak guztiz anomaloak diren pro- pietate elektronikoak dauzka, Dirac-en motako sei puntuz besterik ez osatutako Fermi gainazala duelarik. Honen ondorioz, eroapen elektroiak masa gabekoak balira bezala higitzen dira mobilitate elektronikoa areagotuz. Propietate berezi hauetaz baliatuko liratekeen aplikazio teknologiko posibleek[2] material honekiko interesa egun arlo zienti- fikotik at ere hedatzea eragin du. Grafenoaren sintesiaren errekonozimendu gisa Geim eta Novoselov-ek 2010ean fisikaren Nobel saria lortu zuten. Hala ere, grafenoa ez da sintetiza daitekeen material bidimentsional bakarra. Grafenoa lortzeko teknika bera erabiliz (banantze mikromekanikoa), Geim eta Novoselov-ek zu- zendutako taldeak M oS2 eta N bSe2 sintetizatzea lortu zuen[3]. Konkretuki, M oS2 mo- nogeruza erdieroalea izanik transistoreak minimizatzeko prozesuan silizioaren ordezkari gisa jarduteko hautagaia da. Hala ere, hau egin ahal izateko bere propietate elektro- nikoak sakonkiago aztertzea komeni da. Gradu amaierako lan honetan material honen egitura elektronikoaren eta magnetikoaren karakterizazio teorikoan aurrerapauso txiki bat egitea izan dugu helburu. Horrez gain, W S2 materiala ere era berean landu da, tungsteno atomoa pisutsuagoa izatean, spin-orbita elkarrekintzaren eragina nabariagoa izatea espero baita. Modu honetan, lan hau hiru atal nagusitan banatzen da. Lehenengoa teoriari dago- kio, DF T (Dentsitatearen Funtzionalaren Teoria) inplementatzeko oinarri teorikoa lan- du delarik. Magnetizazioa aztertzeko ezinbestekoa den espina inplementatzeko modua ere aztertu da, eta baita egin beharreko hurbilketen eta pseudopotentzialen metodoaren azalpen bat eman ere. Bigarren atalean QuantumEspresso kodea erabiliz burututako ab-initio kalkuluen deskripzio eta emaitzak aurkeztu dira, azkenei dagokien interpreta- zioa eginez. Bertan M oS2 -n bolumenetiketik monogeruzara pasatzeak egitura elektroni- koan duen eragina aztertu da, ondoren M oS2 eta W S2 monogeruzen banda egitura eta magnetizazioan analisi sakonagoa eginez. Azkenengo atalean ateratako ondorioak idatzi dira, etorkizunerako lanetarako ateak zabalduz.

Ab Initio Study Of Zno-Based Gas-Sensing Mechanisms: Surface Reconstruction And Charge Transfer

Density functional theory (DFT) calculations were employed to explore the gas-sensing mechanisms of zinc oxide (ZnO) with surface reconstruction taken into consideration. Mix-terminated (10 (1) over bar0) ZnO surfaces were examined. By simulating the adsorption process of various gases, i.e., H-2, NH3, CO, and ethanol (C2H5OH) gases, on the ZnO (10 (1) over bar0) surface, the changes of configuration and electronic structure were compared. Based on these calculations, two gas-sensing mechanisms were proposed and revealed that both surface reconstruction and charge transfer result in a change of electronic conductance of ZnO. Also, the calculations were compared with existing experiments.

Band Engineering in Thermoelectric Materials Using Optical, Electronic, and Ab-Initio Computed Properties

Thermoelectric materials have demanded a significant amount of attention for their ability to convert waste heat directly to electricity with no moving parts. A resurgence in thermoelectrics research has led to significant enhancements in the thermoelectric figure of merit, zT, even for materials that were already well studied. This thesis approaches thermoelectric zT optimization by developing a detailed understanding of the electronic structure using a combination of electronic/thermoelectric properties, optical properties, and ab-initio computed electronic band structures. This is accomplished by applying these techniques to three important classes of thermoelectric materials: IV-VI materials (the lead chalcogenides), Half-Heusler’s (XNiSn where X=Zr, Ti, Hf), and CoSb₃ skutterudites.

In the IV-VI materials (PbTe, PbSe, PbS) I present a shifting temperature-dependent optical absorption edge which correlates well to the computed ab-initio molecular dynamics result. Contrary to prior literature that suggests convergence of the primary and secondary bands at 400 K, I suggest a higher convergence temperature of 700, 900, and 1000 K for PbTe, PbSe, and PbS, respectively. This finding can help guide electronic properties modelling by providing a concrete value for the band gap and valence band offset as a function of temperature.

Another important thermoelectric material, ZrNiSn (half-Heusler), is analyzed for both its optical and electronic properties; transport properties indicate a largely different band gap depending on whether the material is doped n-type or p-type. By measuring and reporting the optical band gap value of 0.13 eV, I resolve the discrepancy in the gap calculated from electronic properties (maximum Seebeck and resistivity) by correlating these estimates to the electron-to-hole weighted mobility ratio, A, in narrow gap materials (A is found to be approximately 5.0 in ZrNiSn).

I also show that CoSb₃ contains multiple conduction bands that contribute to the thermoelectric properties. These bands are also observed to shift towards each other with temperature, eventually reaching effective convergence for T>500 K. This implies that the electronic structure in CoSb₃ is critically important (and possibly engineerable) with regards to its high thermoelectric figure of merit.

Ab initio study of anharmonicity in high pressure Cmca-4 metallic hydrogen

Hydrogen is the only atom for which the Schr odinger equation is solvable. Consisting only of a proton and an electron, hydrogen is the lightest element and, nevertheless, is far from being simple. Under ambient conditions, it forms diatomic molecules H2 in gas phase, but di erent temperature and pressures lead to a complex phase diagram, which is not completely known yet. Solid hydrogen was rst documented in 1899 [1] and was found to be isolating. At higher pressures, however, hydrogen can be metallized. In 1935 Wigner and Huntington predicted that the metallization pressure would be 25 GPa [2], where molecules would disociate to form a monoatomic metal, as alkali metals that lie below hydrogen in the periodic table. The prediction of the metallization pressure turned out to be wrong: metallic hydrogen has not been found yet, even under a pressure as high as 320 GPa. Nevertheless, extrapolations based on optical measurements suggest that a metallic phase may be attained at 450 GPa [3]. The interest of material scientist in metallic hydrogen can be attributed, at least to a great extent, to Ashcroft, who in 1968 suggested that such a system could be a hightemperature superconductor [4]. The temperature at which this material would exhibit a transition from a superconducting to a non-superconducting state (Tc) was estimated to be around room temperature. The implications of such a statement are very interesting in the eld of astrophysics: in planets that contain a big quantity of hydrogen and whose temperature is below Tc, superconducting hydrogen may be found, specially at the center, where the gravitational pressure is high. This might be the case of Jupiter, whose proportion of hydrogen is about 90%. There are also speculations suggesting that the high magnetic eld of Jupiter is due to persistent currents related to the superconducting phase [5]. Metallization and superconductivity of hydrogen has puzzled scientists for decades, and the community is trying to answer several questions. For instance, what is the structure of hydrogen at very high pressures? Or a more general one: what is the maximum Tc a phonon-mediated superconductor can have [6]? A great experimental e ort has been carried out pursuing metallic hydrogen and trying to answer the questions above; however, the characterization of solid phases of hydrogen is a hard task. Achieving the high pressures needed to get the sought phases requires advanced technologies. Diamond anvil cells (DAC) are commonly used devices. These devices consist of two diamonds with a tip of small area; for this reason, when a force is applied, the pressure exerted is very big. This pressure is uniaxial, but it can be turned into hydrostatic pressure using transmitting media. Nowadays, this method makes it possible to reach pressures higher than 300 GPa, but even at this pressure hydrogen does not show metallic properties. A recently developed technique that is an improvement of DAC can reach pressures as high as 600 GPa [7], so it is a promising step forward in high pressure physics. Another drawback is that the electronic density of the structures is so low that X-ray di raction patterns have low resolution. For these reasons, ab initio studies are an important source of knowledge in this eld, within their limitations. When treating hydrogen, there are many subtleties in the calculations: as the atoms are so light, the ions forming the crystalline lattice have signi cant displacements even when temperatures are very low, and even at T=0 K, due to Heisenberg's uncertainty principle. Thus, the energy corresponding to this zero-point (ZP) motion is signi cant and has to be included in an accurate determination of the most stable phase. This has been done including ZP vibrational energies within the harmonic approximation for a range of pressures and at T=0 K, giving rise to a series of structures that are stable in their respective pressure ranges [8]. Very recently, a treatment of the phases of hydrogen that includes anharmonicity in ZP energies has suggested that relative stability of the phases may change with respect to the calculations within the harmonic approximation [9]. Many of the proposed structures for solid hydrogen have been investigated. Particularly, the Cmca-4 structure, which was found to be the stable one from 385-490 GPa [8], is metallic. Calculations for this structure, within the harmonic approximation for the ionic motion, predict a Tc up to 242 K at 450 GPa [10]. Nonetheless, due to the big ionic displacements, the harmonic approximation may not su ce to describe correctly the system. The aim of this work is to apply a recently developed method to treat anharmonicity, the stochastic self-consistent harmonic approximation (SSCHA) [11], to Cmca-4 metallic hydrogen. This way, we will be able to study the e ects of anharmonicity in the phonon spectrum and to try to understand the changes it may provoque in the value of Tc. The work is structured as follows. First we present the theoretical basis of the calculations: Density Functional Theory (DFT) for the electronic calculations, phonons in the harmonic approximation and the SSCHA. Then we apply these methods to Cmca-4 hydrogen and we discuss the results obtained. In the last chapter we draw some conclusions and propose possible future work.