WorldWideScience

Sample records for ab initio studies

  1. Ab Initio Studies of Stratospheric Ozone Depletion Chemistry

    Lee, Timothy J.; Head-Gordon, Martin; Langhoff, Stephen R. (Technical Monitor)

    1995-01-01

    An overview of the current understanding of ozone depletion chemistry, particularly with regards the formation of the so-called Antarctic ozone hole, will be presented together with an outline as to how ab initio quantum chemistry can be used to further our understanding of stratospheric chemistry. The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of fluorine, chlorine, bromine and nitrogen oxide species will be demonstrated by presentation of some example studies. The ab initio results will be shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the theoretical results are shown to fill in the gaps and to resolve experimental controversies. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of halogen oxide species will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of experimental studies.

  2. Ab-initio study of thermoelectricity of layered tellurium compounds

    Ibarra Hernández, Wilfredo

    2015-01-01

    In this thesis, we explore the electronic, dynamic and thermoelectric properties of different tellurium-based compounds. We perform ab-initio calculations within the Vienna Ab-initio Simulation Package (VASP) that works in the framework of Density Functional Theory (DFT). For the thermoelectric properties, we use the Boltztrap code that solves the Boltzmann Transport Equations (BTE) for electrons within the Constant Relaxation Time Approximation (CRTA). This computational pa...

  3. Ab-initio study of napthelene based conducting polymer

    Ruhela, Ankur [Advanced Materials Research Group, Computational Nanoscience and Technology Lab (CNTL), ABV-Indian Institute of Information Technology and Management, Gwalior -474010, India and Amity Institute of Nanotechnology, Amity University, Noida-201303 (India); Kanchan, Reena, E-mail: reena.kanchan1977@gmail.com [Department of Chemistry, Jiwaji University, Gwalior-474001 (India); Srivastava, Anurag [Advanced Materials Research Group, Computational Nanoscience and Technology Lab (CNTL), ABV-Indian Institute of Information Technology and Management, Gwalior -474010 (India); Sinha, O. P. [Amity Institute of Nanotechnology, Amity University, Noida-201303 (India)

    2014-04-24

    In this paper, we have identified structural and electronic properties of conducting polymers by using DFT based ATK-VNL ab-initio tool. Naphthalene derivative structures were stabilized by varying the bond length between two atoms of the molecule C-N and C-C. We have also studied the molecular energy spectrum of naphthalene derivatives and found the HOMOLUMO for the same. A comparison of structural and electronic properties of naphthalene derivatives by attaching the functional group of amine, have been performed and found that they show good semi conducting properties.

  4. Ab-initio study of napthelene based conducting polymer

    In this paper, we have identified structural and electronic properties of conducting polymers by using DFT based ATK-VNL ab-initio tool. Naphthalene derivative structures were stabilized by varying the bond length between two atoms of the molecule C-N and C-C. We have also studied the molecular energy spectrum of naphthalene derivatives and found the HOMOLUMO for the same. A comparison of structural and electronic properties of naphthalene derivatives by attaching the functional group of amine, have been performed and found that they show good semi conducting properties

  5. Ab initio study of phase equilibria in TiCx

    Korzhavyi, P.A.; Pourovskii, L.V.; Hugosson, H.W.;

    2002-01-01

    The phase diagram for the vacancy-ordered structures in the substoichiometric TiCx (x = 0.5-1.0) has been established from Monte Carlo simulations with the long-range pair and multisite effective interactions obtained from ab initio calculations. Three ordered superstructures of vacancies (Ti2C, ...

  6. Relaxation of Small Molecules: an ab initio Study

    CAO Yi-Gang; JIAO Zheng-Kuan; A. Antons; K. Schroeder; S. Blügel2

    2002-01-01

    Using an ab initio total energy and force method, we have relaxed several group IV and group V elementalclusters, in detail the arsenic and antimony dimers, silicon, phosphorus, arsenic and antimony tetramers. The obtainedbond lengths and cohesive energies are more accurate than other calculating methods, and in excellent agreement withthe experimental results.

  7. Molecular ion LiHe+: ab initio study

    Highlights: ► Excited electronic states of LiHe+ are studied. ► Potential energy curves of thirteen states are calculated. ► Dipole moment and transition dipole moment functions are determined. ► Basic spectroscopic properties of the electronic states are derived. - Abstract: High level ab initio calculations are performed on the molecular ion LiHe+. Potential energy curves for the low-lying singlet and triplet electronic states are calculated using the multi-reference configuration interaction and single-reference coupled cluster methods with large basis sets. The corresponding dipole moments and transition dipole moments functions are also determined. The basic spectroscopic properties and excitation energies of the electronic states are derived from rovibrational bound state calculations.

  8. Serine Proteases an Ab Initio Molecular Dynamics Study

    De Santis, L

    1999-01-01

    In serine proteases (SP's), the H-bond between His-57 and Asp-102, and that between Gly-193 and the transition state intermediate play a crucial role for enzymatic function. To shed light on the nature of these interactions, we have carried out ab initio molecular dynamics simulations on complexes representing adducts between the reaction intermediate and elastase (one protein belonging to the SP family). Our calculations indicate the presence of a low--barrier H-bond between His-57 and Asp-102, in complete agreement with NMR experiments on enzyme--transition state analog complexes. Comparison with an ab initio molecular dynamics simulation on a model of the substrate--enzyme adduct indicates that the Gly-193--induced strong stabilization of the intermediate is accomplished by charge/dipole interactions and not by H-bonding as previously suggested. Inclusion of the protein electric field in the calculations does not affect significantly the charge distribution.

  9. Experimental, ab initio and density functional theory studies on sulfadiazine

    Ogruc-Ildiz, Gulce; Akyuz, Sevim; Ozel, Aysen E.

    2009-04-01

    In the present study, combined experimental and computational study on molecular vibrations of free sulfadiazine has been reported. The theoretically possible stable conformers of free sulfadiazine molecule in electronically ground state were searched by means of torsion potential energy surfaces scan studies through C1 sbnd C7 sbnd S8 sbnd N8, C7 sbnd S8 sbnd N9 sbnd C10 and S8 sbnd N9 sbnd C10 sbnd N11 dihedral angles, at both semi-empirical PM3 and B3LYP/3-21G levels of theory. The final geometrical parameters for the obtained stable conformers were determined by means of geometry optimization carried out at ab initio HF/6-31G++(d,p) and DFT/B3LYP/6-31G++(d,p) theory levels. The harmonic and anharmonic vibrational wavenumbers and IR intensities were calculated at the same theory levels used in geometry optimization. The modes of the fundamental vibrations were characterized depending on their the total energy distribution (TED%). In order to fit the calculated harmonic wavenumbers to experimental ones, dual scale factors were used. The experimental infrared and Raman spectra of sulfadiazine in solid phase have been measured and compared with the calculated vibrational spectra of each conformer.

  10. Ab initio study of optical excitations in VO2

    Coulter, John; Gali, Adam; Manousakis, Efstratios

    2014-03-01

    Motivated by recent experimental efforts to fabricate p-n junctions from transition metal oxides (TMOs) and a recent theoretical study claiming TMOs to be good absorbers and promising materials for efficient carrier multiplication, we study the optical properties of a prototypical TMO, the insulator M1 phase of vanadium dioxide (VO2), by ab initio methods. We applied the Bethe-Salpeter equations (BSE) to calculate the optical properties, starting from self-consistent GW quasi-particle energy levels and states. In contrast to expectations, the exciton binding energy obtained by BSE is in good agreement with the experiment. We find that the electron-electron interaction is very strong which makes this material promising for efficient carrier multiplication that might lead to an enhanced efficiency in photo-voltaics applications. To illustrate this more quantitatively, we calculated the impact ionization rate within the independent quasiparticle approximation, and find that the rate is significantly higher than silicon in the region of highest solar intensity, due to the strong multiple carrier excitations.

  11. Ab initio study of neutron drops with chiral Hamiltonians

    We report ab initio calculations for neutron drops in a 10 MeV external harmonic-oscillator trap using chiral nucleon–nucleon plus three-nucleon interactions. We present total binding energies, internal energies, radii and odd–even energy differences for neutron numbers N=2–18 using the no-core shell model with and without importance truncation. Furthermore, we present total binding energies for N=8,16,20,28,40,50 obtained in a coupled-cluster approach. Comparisons with quantum Monte Carlo results, where available, using Argonne v8′ with three-nucleon interactions reveal important dependences on the chosen Hamiltonian

  12. Tailoring magnetoresistance at the atomic level: An ab initio study

    Tao, Kun

    2012-01-05

    The possibility of manipulating the tunneling magnetoresistance (TMR) of antiferromagnetic nanostructures is predicted in the framework of ab initio calculations. By the example of a junction composed of an antiferromagnetic dimer and a spin-polarized scanning tunneling microscopy tip we show that the TMR can be tuned and even reversed in sign by lateral and vertical movements of the tip. Moreover, our finite-bias calculations demonstrate that the magnitude and the sign of the TMR can also be tuned by an external voltage. © 2012 American Physical Society.

  13. Ab initio study of neutron drops with chiral Hamiltonians

    H.D. Potter

    2014-12-01

    Full Text Available We report ab initio calculations for neutron drops in a 10 MeV external harmonic-oscillator trap using chiral nucleon–nucleon plus three-nucleon interactions. We present total binding energies, internal energies, radii and odd–even energy differences for neutron numbers N=2–18 using the no-core shell model with and without importance truncation. Furthermore, we present total binding energies for N=8,16,20,28,40,50 obtained in a coupled-cluster approach. Comparisons with quantum Monte Carlo results, where available, using Argonne v8′ with three-nucleon interactions reveal important dependences on the chosen Hamiltonian.

  14. Ab initio and kinetic modeling studies of formic acid oxidation

    Marshall, Paul; Glarborg, Peter

    2015-01-01

    A detailed chemical kinetic model for oxidation of formic acid (HOCHO) in flames has been developed, based on theoretical work and data from literature. Ab initio calculations were used to obtain rate coefficients for reactions of HOCHO with H, O, and HO2. Modeling predictions with the mechanism...... on calculations with the kinetic model. Formic acid is consumed mainly by reaction with OH, yielding OCHO, which dissociates rapidly to CO2 + H, and HOCO, which may dissociate to CO + OH or CO2 + H, or react with H, OH, or O2 to form more stable products. The branching fraction of the HOCHO + OH...

  15. Ab Initio Studies of Halogen and Nitrogen Oxide Species of Interest in Stratospheric Chemistry

    Lee, Timothy J.; Langhoff, Stephen R. (Technical Monitor)

    1995-01-01

    The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of fluorine, chlorine, bromine and nitrogen oxide species will be demonstrated by presentation of some example studies. The ab initio results are shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the theoretical results are shown to fill in the gaps and to resolve experimental controversies. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of halogen oxide species will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of experimental studies.

  16. Lithium Insertion In Silicon Nanowires: An ab Initio Study

    Zhang, Qianfan

    2010-09-08

    The ultrahigh specific lithium ion storage capacity of Si nanowires (SiNWs) has been demonstrated recently and has opened up exciting opportunities for energy storage. However, a systematic theoretical study on lithium insertion in SiNWs remains a challenge, and as a result, understanding of the fundamental interaction and microscopic dynamics during lithium insertion is still lacking. This paper focuses on the study of single Li atom insertion into SiNWs with different sizes and axis orientations by using full ab initio calculations. We show that the binding energy of interstitial Li increases as the SiNW diameter grows. The binding energies at different insertion sites, which can be classified as surface, intermediate, and core sites, are quite different. We find that surface sites are energetically the most favorable insertion positions and that intermediate sites are the most unfavorable insertion positions. Compared with the other growth directions, the [110] SiNWs with different diameters always present the highest binding energies on various insertion locations, which indicates that [110] SiNWs are more favorable by Li doping. Furthermore, we study Li diffusion inside SiNWs. The results show that the Li surface diffusion has a much higher chance to occur than the surface to core diffusion, which is consistent with the experimental observation that the Li insertion in SiNWs is layer by layer from surface to inner region. After overcoming a large barrier crossing surface-to-intermediate region, the diffusion toward center has a higher possibility to occur than the inverse process. © 2010 American Chemical Society.

  17. An ab initio study of plutonium oxides surfaces

    By means of first-principles calculations, we have studied the atomic structure as well as the thermodynamic stability of various plutonium dioxide surfaces in function of their environment (in terms of oxygen partial pressure and temperature). All these simulations have been performed with the ABINIT code. It is well known that DFT fails to describe correctly plutonium-based materials since 5f electrons in such systems are strongly correlated. In order to go beyond DFT, we have treated PuO2 and β-Pu2O3 in a DFT+U framework. We show that the couple of parameters (U,J) that works well for pure Pu is also well designed for describing ground state (GS) properties of these two oxides. The major improvement with respect with DFT is that we are able to predict an insulating GS in agreement with experiments. The presence of a gap in the DOS (Density of States) of plutonium oxides should play a significant role in the predicted surface reactivity. However, performing DFT+U calculations on surfaces of plutonium oxide from scratch was too ambitious. That is why we decided, as a first step, to study the stability of the (100), (110) and (111) surfaces of PuO2 in a DFT-GGA framework. For each of these orientations, we considered various terminations. These ab initio results have been introduced in a thermodynamic model which allows us to predict the relative stability of the different terminations as a function of temperature and oxygen partial pressure (pO2). We conclude that at room temperature and for pO2∼10 atm., the polar O2-(100) termination is favoured. The stabilization of such a polar stoichiometric surface is surprising and should be confirmed by DFT+U calculations before any final conclusion. (authors)

  18. Highly anisotropic thermal conductivity of arsenene: An ab initio study

    Zeraati, Majid; Vaez Allaei, S. Mehdi; Abdolhosseini Sarsari, I.; Pourfath, Mahdi; Donadio, Davide

    2016-02-01

    Elemental two-dimensional (2D) materials exhibit intriguing heat transport and phononic properties. Here we have investigated the lattice thermal conductivity of newly proposed arsenene, the 2D honeycomb structure of arsenic, using ab initio calculations. Solving the Boltzmann transport equation for phonons, we predict a highly anisotropic thermal conductivity of 30.4 and 7.8 W/mK along the zigzag and armchair directions, respectively, at room temperature. Our calculations reveal that phonons with mean free paths between 20 nm and 1 μ m provide the main contribution to the large thermal conductivity in the zigzag direction; mean free paths of phonons contributing to heat transport in the armchair directions range between 20 and 100 nm. The obtained anisotropic thermal conductivity and feasibility of synthesis, in addition to high electron mobility reported elsewhere, make arsenene a promising material for nanoelectronic applications and thermal management.

  19. Ab initio study of II-(VI){sub 2} dichalcogenides

    Olsson, P; Vidal, J; Lincot, D, E-mail: polsson@kth.se [Institut de R and D sur l' energie photovoltaique (IRDEP), UMR 7174-EDF-CNRS-ENSCP, 6 quai Watier, 78401 Chatou Cedex (France)

    2011-10-12

    The structural stabilities of the (Zn,Cd)(S,Se,Te){sub 2} dichalcogenides have been determined ab initio. These compounds are shown to be stable in the pyrite phase, in agreement with available experiments. Structural parameters for the ZnTe{sub 2} pyrite semiconductor compound proposed here are presented. The opto-electronic properties of these dichalcogenide compounds have been calculated using quasiparticle GW theory. Bandgaps, band structures and effective masses are proposed as well as absorption coefficients and refraction indices. The compounds are all indirect semiconductors with very flat conduction band dispersion and high absorption coefficients. The work functions and surface properties are predicted. The Te and Se based compounds could be of interest as absorber materials in photovoltaic applications. (paper)

  20. Ab initio study of the transition-metal carbene cations

    李吉海; 冯大诚; 冯圣玉

    1999-01-01

    The geometries and bonding characteristics of the first-row transition-metal carbene cations MCH2+ were investigated by ab initio molecular orbital theory (HF/LANL2DZ). All of MCH2+ are coplanar. In the closed shell structures the C bonds to M with double bonds; while in the open shell structures the partial double bonds are formed, because one of the σ and π orbitals is singly occupied. It is mainly the π-type overlap between the 2px orbital of C and 4px, 3dxz, orbitals of M+ that forms the π orbitals. The dissociation energies of C—M bond appear in periodic trend from Sc to Cu. Most of the calculated bond dissociation energies are close to the experimental ones.

  1. Ab initio studies of phoshorene island single electron transistor.

    Ray, S J; Venkata Kamalakar, M; Chowdhury, R

    2016-05-18

    Phosphorene is a newly unveiled two-dimensional crystal with immense potential for nanoelectronic and optoelectronic applications. Its unique electronic structure and two dimensionality also present opportunities for single electron devices. Here we report the behaviour of a single electron transistor (SET) made of a phosphorene island, explored for the first time using ab initio calculations. We find that the band gap and the charging energy decrease monotonically with increasing layer numbers due to weak quantum confinement. When compared to two other novel 2D crystals such as graphene and MoS2, our investigation reveals larger adsorption energies of gas molecules on phosphorene, which indicates better a sensing ability. The calculated charge stability diagrams show distinct changes in the presence of an individual molecule which can be applied to detect the presence of different molecules with sensitivity at a single molecular level. The higher charging energies of the molecules within the SET display operational viability at room temperature, which is promising for possible ultra sensitive detection applications. PMID:27093536

  2. Ab initio studies of phosphorene island single electron transistor

    Ray, S. J.; Venkata Kamalakar, M.; Chowdhury, R.

    2016-05-01

    Phosphorene is a newly unveiled two-dimensional crystal with immense potential for nanoelectronic and optoelectronic applications. Its unique electronic structure and two dimensionality also present opportunities for single electron devices. Here we report the behaviour of a single electron transistor (SET) made of a phosphorene island, explored for the first time using ab initio calculations. We find that the band gap and the charging energy decrease monotonically with increasing layer numbers due to weak quantum confinement. When compared to two other novel 2D crystals such as graphene and MoS2, our investigation reveals larger adsorption energies of gas molecules on phosphorene, which indicates better a sensing ability. The calculated charge stability diagrams show distinct changes in the presence of an individual molecule which can be applied to detect the presence of different molecules with sensitivity at a single molecular level. The higher charging energies of the molecules within the SET display operational viability at room temperature, which is promising for possible ultra sensitive detection applications.

  3. Ab initio study of MoS2 nanotube bundles

    Verstraete, Matthieu; Charlier, Jean-Christophe

    2003-07-01

    Recently, the synthesis of a new phase of MoS2I1/3 stoichiometry was reported [M. Remskar, A. Mrzel, Z. Skraba, A. Jesih, M. Ceh, J. Demšar, P. Stadelmann, F. Lévy, and D. Mihailovic, Science 292, 479 (2001)]. Electron microscope images and diffraction data were interpreted to indicate bundles of sub-nanometer-diameter single-wall MoS2 nanotubes. After experimental characterization, the structure was attributed to an assembly of “armchair” nanotubes with interstitial iodine. Using first-principles total-energy calculations, bundles of MoS2 nanotubes with different topologies and stoichiometries are investigated. All of the systems are strongly metallic. Configurations with “zigzag” structures are found to be more stable energetically than the “armchair” ones, though all of the structures have similar stabilities. After relaxation, there remain several candidates which give a lattice parameter in relative agreement with experiment. Further, spin-polarized calculations indicate that a structure with armchair tubes iodine atoms in their center acquires a very large spontaneous magnetic moment of 12μB, while the other structures are nonmagnetic. Our ab initio calculations show that in most of the other structures, the tubes are very strongly bound together, and that the compounds should be considered as a crystal, rather than as a bundle of tubes in the habitual sense.

  4. An ab initio study of plutonium oxides surfaces; Etude ab initio des surfaces d'oxydes de Pu

    Jomard, G.; Bottin, F.; Amadon, B

    2007-07-01

    By means of first-principles calculations, we have studied the atomic structure as well as the thermodynamic stability of various plutonium dioxide surfaces in function of their environment (in terms of oxygen partial pressure and temperature). All these simulations have been performed with the ABINIT code. It is well known that DFT fails to describe correctly plutonium-based materials since 5f electrons in such systems are strongly correlated. In order to go beyond DFT, we have treated PuO{sub 2} and {beta}-Pu{sub 2}O{sub 3} in a DFT+U framework. We show that the couple of parameters (U,J) that works well for pure Pu is also well designed for describing ground state (GS) properties of these two oxides. The major improvement with respect with DFT is that we are able to predict an insulating GS in agreement with experiments. The presence of a gap in the DOS (Density of States) of plutonium oxides should play a significant role in the predicted surface reactivity. However, performing DFT+U calculations on surfaces of plutonium oxide from scratch was too ambitious. That is why we decided, as a first step, to study the stability of the (100), (110) and (111) surfaces of PuO{sub 2} in a DFT-GGA framework. For each of these orientations, we considered various terminations. These ab initio results have been introduced in a thermodynamic model which allows us to predict the relative stability of the different terminations as a function of temperature and oxygen partial pressure (p{sub O{sub 2}}). We conclude that at room temperature and for p{sub O{sub 2}}{approx}10 atm., the polar O{sub 2}-(100) termination is favoured. The stabilization of such a polar stoichiometric surface is surprising and should be confirmed by DFT+U calculations before any final conclusion. (authors)

  5. Mechanisms for catalytic carbon nanofiber growth studied by ab initio density functional theory calculations

    Abild-Pedersen, Frank; Nørskov, Jens Kehlet; Rostrup-Nielsen, Jens;

    2006-01-01

    Mechanisms and energetics of graphene growth catalyzed by nickel nanoclusters were studied using ab initio density functional theory calculations. It is demonstrated that nickel step-edge sites act as the preferential growth centers for graphene layers on the nickel surface. Carbon is transported...

  6. Ab Initio Studies of Chlorine Oxide and Nitrogen Oxide Species of Interest in Stratospheric Chemistry

    Lee, Timothy J.; Langhoff, Stephen R. (Technical Monitor)

    1995-01-01

    The ability of modern state-of-the art ab initio quantum chemical techniques to characterize reliably the gas-phase molecular structure, vibrational spectrum, electronic spectrum, and thermal stability of chlorine oxide and nitrogen oxide species will be demonstrated by presentation of some example studies. In particular the geometrical structures, vibrational spectra, and heats of formation Of ClNO2, CisClONO, and trans-ClONO are shown to be in excellent agreement with the available experimental data, and where the experimental data are either not known or are inconclusive, the ab initio results are shown to fill in the gaps and to resolve the experimental controversy. In addition, ab initio studies in which the electronic spectra and the characterization of excited electronic states of ClONO2, HONO2, ClOOC17 ClOOH, and HOOH will also be presented. Again where available, the ab initio results are compared to experimental observations, and are used to aid in the interpretation of the experimental studies.

  7. Relaxation of Small Molecules:an ab initio Study

    CAOYi-Gang; A.Antons; 等

    2002-01-01

    Using an ab inito total energy and force method,we have relaxed several group IV and group V elemental clusters,in detail the arsenic and antimony dimers,silicon,phosphorus,arsenic and antimony tetraners,The obtained bond lengths and cohesive energies are more accurate than other calculating methods,and in excellent agreement with the experimental results.

  8. Ab initio studies of equations of state and chemical reactions of reactive structural materials

    Zaharieva, Roussislava

    subject of studies of the shock or thermally induced chemical reactions of the two solids comprising these reactive materials, from first principles, is a relatively new field of study. The published literature on ab initio techniques or quantum mechanics based approaches consists of the ab initio or ab initio-molecular dynamics studies in related fields that contain a solid and a gas. One such study in the literature involves a gas and a solid. This is an investigation of the adsorption of gasses such as carbon monoxide (CO) on Tungsten. The motivation for these studies is to synthesize alternate or synthetic fuel technology by Fischer-Tropsch process. In this thesis these studies are first to establish the procedure for solid-solid reaction and then to extend that to consider the effects of mechanical strain and temperature on the binding energy and chemisorptions of CO on tungsten. Then in this thesis, similar studies are also conducted on the effect of mechanical strain and temperature on the binding energies of Titanium and hydrogen. The motivations are again to understand the method and extend the method to such solid-solid reactions. A second motivation is to seek strained conditions that favor hydrogen storage and strain conditions that release hydrogen easily when needed. Following the establishment of ab initio and ab initio studies of chemical reactions between a solid and a gas, the next step of research is to study thermally induced chemical reaction between two solids (Ni+Al). Thus, specific new studies of the thesis are as follows: (1) Ab initio Studies of Binding energies associated with chemisorption of (a) CO on W surfaces (111, and 100) at elevated temperatures and strains and (b) adsorption of hydrogen in titanium base. (2) Equations of state of mixtures of reactive material structures from ab initio methods. (3) Ab initio studies of the reaction initiation, transition states and reaction products of intermetallic mixtures of (Ni+Al) at elevated

  9. Ab Initio Study on the Anti-HIV Activities of Hydroxyflavones

    ZHANG Yu

    2005-01-01

    Flavone and 95 hydroxyflavones have been studied with ab initio method, and their total energies, atomic charges, dipole moments, multipole moments, molecular orbital compositions, orbital energies etc. were obtained. Among them the relationship between total atomic charges and activities against HIV is basically in accordance with the experimental results. The beneficial references are provided for the extraction and synthesis of strong active anti-HIV medicines.

  10. Ab initio study on the magnetostructural properties of MnAs

    Sanvito, Stefano; RUNGGER, IVAN

    2006-01-01

    The magnetic and structural properties of MnAs are studied with ab initio methods, and by mapping total energies onto a Heisenberg model. The stability of the different phases is found to depend mainly on the volume and on the amount of magnetic order, confirming previous experimental findings and phenomenological models. It is generally found that for large lattice constants the ferromagnetic state is favored, whereas for small lattice constants different antiferromagnetic states can be stab...

  11. Ab Initio Studies on Hematite Surface and the Adsorption of Phosphate

    Swati Chaudhury; Chandrika Varadachari; Kunal Ghosh

    2014-01-01

    This investigation explores the ab initio DFT method for understanding surface structure of hematite and the nature and energetics of phosphate adsorption. Using the full potential linearized plane wave method (FP-LAPW), we derived the structure and energies of various magnetic forms of hematite. The antiferromagnetic (AFM) form was observed to be the most stable. Hematite surfaces with Fe-termination, O-termination, or OH-termination were studied. The OH-terminated surface was the most stabl...

  12. Hydrogen adsorption in ZIF-7: A DFT and ab-initio molecular dynamics study

    Dixit, Mudit; Major, Dan Thomas; Pal, Sourav

    2016-05-01

    Primary H2 adsorption sites in a zeolitic imidazolate framework, ZIF-7, are identified using ab-initio density functional theory (DFT) based molecular dynamics annealing simulations. The simulations suggest several low energy adsorption sites. The effect of light transition metal decoration on hydrogen storage properties was studied. Our ab-intio DFT calculations illustrate that decorating the ZIF with Sc increases both the number of H2 molecules, as well as the H2 binding energy. The binding energy (∼25 kJ/mol per H2) at 8H2 loading in the pore, suggests that Sc-ZIFs can be potential candidates for hydrogen storage.

  13. Ab-initio approach to study hydrogen diffusion in 9Cr steels

    We calculate the equilibrium energies and migration barriers of Fe, Cr and H interstitial defects in α-FeX(X=Cr). We use the ab-initio electronic structure code, SIESTA, coupled to the monomer method to find activated states (or migration barriers), in order to study atomic defects diffusion. Ab-initio calculations reveal that in the presence of Cr the H migration barriers are higher than in pure α-Fe. On the other hand, our permeation tests on 9%Cr-91%Fe alloys reveal a permeation coefficient 10 times lower and a diffusion coefficient 200 times lower than in pure, annealed iron. Focusing on our experimental results, we explore very simple model of new H trapping sites and possible migration paths that can explain the experimental observations.

  14. Study of atomic structure of liquid Hg-In alloys using ab-initio molecular dynamics

    Sharma, Nalini; Ahluwalia, P. K. [Department of Physics, Himachal Pradesh University, Shimla(HP)-171005 (India); Thakur, Anil [Department of Physics, Govt. P. G. College Solan (HP)-173212 (India)

    2015-05-15

    Ab-initio molecular dynamics simulations are performed to study the structural properties of liquid Hg-In alloys. The interatomic interactions are described by ab-initio pseudopotentials given by Troullier and Martins. Five liquid Hg-In mixtures (Hg{sub 10}In{sub 90}, Hg{sub 30}In{sub 70}, Hg{sub 50}In{sub 50}, Hg{sub 70}In{sub 30} and Hg{sub 90}In{sub 10}) at 299K are considered. The radial distribution function g(r) and structure factor S(q) of considered alloys are compared with respective experimental results for liquid Hg (l-Hg) and (l-In). The radial distribution function g(r) shows the presence of short range order in the systems considered. Smooth curves of Bhatia-Thornton partial structure factors factor shows the presence of liquid state in the considered alloys.

  15. Study of atomic structure of liquid Hg-In alloys using ab-initio molecular dynamics

    Ab-initio molecular dynamics simulations are performed to study the structural properties of liquid Hg-In alloys. The interatomic interactions are described by ab-initio pseudopotentials given by Troullier and Martins. Five liquid Hg-In mixtures (Hg10In90, Hg30In70, Hg50In50, Hg70In30 and Hg90In10) at 299K are considered. The radial distribution function g(r) and structure factor S(q) of considered alloys are compared with respective experimental results for liquid Hg (l-Hg) and (l-In). The radial distribution function g(r) shows the presence of short range order in the systems considered. Smooth curves of Bhatia-Thornton partial structure factors factor shows the presence of liquid state in the considered alloys

  16. Liquid Be, Ca and Ba. An orbital-free ab-initio molecular dynamics study

    Several static and dynamic properties of liquid beryllium (l-Be), liquid calcium (l-Ca) and liquid barium (l-Ba) near their triple point have been evaluated by the orbital-free ab initio molecular dynamics method (OF-AIMD), where the interaction between valence electrons and ions is described by means of local pseudopotentials. These local pseudopotentials used were constructed through a force-matching process with those obtained from a Kohn-Sham ab initio molecular dynamics study (KS-AIMD) of a reduced system with non-local pseudopotentials. The calculated static structures show good agreement with the available experimental data, including an asymmetric second peak in the structure factor which has been linked to the existence of a marked icosahedral short-range order in the liquid. As for the dynamic properties, we obtain collective density excitations whose associated dispersion relations exhibit a positive dispersion

  17. Electronic properties of liquid Hg-In alloys : Ab-initio molecular dynamics study

    Sharma, Nalini; Thakur, Anil; Ahluwalia, P. K.

    2016-05-01

    Ab-initio molecular dynamics simulations are performed to study the structural properties of liquid Hg-In alloys. The interatomic interactions are described by ab-initio pseudopotentials given by Troullier and Martins. Three liquid Hg-In alloys (Hg10In90, Hg30In70,. Hg50In50, Hg70In30, and Hg90Pb10) at 299 K are considered. The calculated results for liquid Hg (l-Hg) and lead (l-In) are also drawn. Along with the calculated results of considered five liquid alloys of Hg-In alloy. The results obtained from electronic properties namely total density of state and partial density of states help to find the local arrangement of Hg and In atoms and the presence of liquid state in the considered five alloys.

  18. Theoretical study of silicon carbide under irradiation at the nano scale: classical and ab initio modelling

    The behaviour of silicon carbide under irradiation has been studied using classical and ab initio simulations, focusing on the nano scale elementary processes. First, we have been interested in the calculation of threshold displacement energies, which are difficult to determine both experimentally and theoretically, and also the associated Frenkel pairs. In the framework of this thesis, we have carried out simulations in classical and ab initio molecular dynamics. For the classical approach, two types of potentials have been used: the Tersoff potential, which led to non satisfactory results, and a new one which has been developed during this thesis. This potential allows a better modelling of SiC under irradiation than most of the empirical potentials available for SiC. It is based on the EDIP potential, initially developed to describe defects in silicon, that we have generalized to SiC. For the ab initio approach, the feasibility of the calculations has been validated and average energies of 19 eV for the C and 38 eV for the Si sublattices have been determined, close to the values empirically used in the fusion community. The results obtained with the new potential EDIP are globally in agreement with those values. Finally, the elementary processes involved in the crystal recovery have been studied by calculating the stability of the created Frenkel pairs and determining possible recombination mechanisms with the nudged elastic band method. (author)

  19. Ab Initio Studies of Shock-Induced Chemical Reactions of Inter-Metallics

    Zaharieva, Roussislava; Hanagud, Sathya

    2009-06-01

    Shock-induced and shock assisted chemical reactions of intermetallic mixtures are studied by many researchers, using both experimental and theoretical techniques. The theoretical studies are primarily at continuum scales. The model frameworks include mixture theories and meso-scale models of grains of porous mixtures. The reaction models vary from equilibrium thermodynamic model to several non-equilibrium thermodynamic models. The shock-effects are primarily studied using appropriate conservation equations and numerical techniques to integrate the equations. All these models require material constants from experiments and estimates of transition states. Thus, the objective of this paper is to present studies based on ab initio techniques. The ab inito studies, to date, use ab inito molecular dynamics. This paper presents a study that uses shock pressures, and associated temperatures as starting variables. Then intermetallic mixtures are modeled as slabs. The required shock stresses are created by straining the lattice. Then, ab initio binding energy calculations are used to examine the stability of the reactions. Binding energies are obtained for different strain components super imposed on uniform compression and finite temperatures. Then, vibrational frequencies and nudge elastic band techniques are used to study reactivity and transition states. Examples include Ni and Al.

  20. Ab initio studies of magnetism in the iron chalcogenides FeTe and FeSe

    The iron chalcogenides FeTe and FeSe belong to the family of iron-based superconductors. We study the magnetism in these compounds in the normal state using the ab initio downfolding scheme developed for strongly correlated electron systems. In deriving ab initio low-energy effective models, we employ the constrained GW method to eliminate the double counting of electron correlations originating from the exchange correlations already taken into account in the density functional theory. By solving the derived ab initio effective models, we reveal that the elimination of the double counting is important in reproducing the bicollinear antiferromagnetic order in FeTe, as is observed in experiments. We also show that the elimination of the double counting induces a unique degeneracy of several magnetic orders in FeSe, which may explain the absence of the magnetic ordering. We discuss the relationship between the degeneracy and the recently found puzzling phenomena in FeSe as well as the magnetic ordering found under pressure. (author)

  1. Temperature dependent surface relaxation for Al(110) and Mg(10-10) studied by orbital free ab initio molecular dynamics

    Gonzalez, Luis E.; Gonzalez, David J

    2006-01-01

    We have performed orbital free ab initio molecular dynamics simulations in order to study the thermal behaviour of two open surfaces of solid metallic systems, namely the (110) face of fcc Al and the (10-10) face of hcp Mg. Our results reproduce qualitatively both the experimental measurements and previous ab initio calculations performed with the more costly Kohn-Sham approach of Density Functional Theory. These calculations can be viewed as a validation test of the orbital free method for s...

  2. Silicene on metal and metallized surfaces: ab initio studies

    The deposition of silicene on several metals is investigated. For fcc crystals the (111) surfaces while for hexagonal ones the (0001) surfaces are used. The Ca(111)1 × 1 substrate is found to be the most promising candidate. The silicene adsorption on Ca-functionalized Si(111)1 × 1 and 2 × 1 surfaces is also studied. The 1 × 1 substrates lead to overlayer silicene with hexagonal symmetry and Dirac cones. However, the Dirac points are below the Fermi level, and small energy gaps are opened. In the case of 2 × 1 surfaces, strong lattice relaxation occurs. Only rudiments of conical linear bands remain visible. (paper)

  3. Ab initio studies of complexation of anions to neutral species

    Johansson, Patrik [Department of Applied Physics, Chalmers University of Technology, SE-412 96 Goeteborg (Sweden)]. E-mail: patrikj@fy.chalmers.se; Jacobsson, Per [Department of Applied Physics, Chalmers University of Technology, SE-412 96 Goeteborg (Sweden)

    2005-06-30

    The complexation of simple anions (F{sup -} and Cl{sup -}) to different neutral species, anion-coordinating agents, has been studied using electronic structure calculations. The obtained changes in the equilibrium constants for salt dissolution reactions in different typical electrolyte systems are reported. In addition the lithium ion affinities of the obtained anionic complexes have been calculated. Using the present results we discuss strategies for future usage of anion complexing agents and make recommendations of salt and agent combinations for better lithium battery electrolyte performance.

  4. Ab initio molecular dynamics study of liquid methanol

    Handgraaf, J W; Meijer, E J; Handgraaf, Jan-Willem; Erp, Titus S. van; Meijer, Evert Jan

    2003-01-01

    We present a density-functional theory based molecular-dynamics study of the structural, dynamical, and electronic properties of liquid methanol under ambient conditions. The calculated radial distribution functions involving the oxygen and hydroxyl hydrogen show a pronounced hydrogen bonding and compare well with recent neutron diffraction data, except for an underestimate of the oxygen-oxygen correlation. We observe that, in line with infrared spectroscopic data, the hydroxyl stretching mode is significantly red-shifted in the liquid. A substantial enhancement of the dipole moment is accompanied by significant fluctuations due to thermal motion. Our results provide valuable data for improvement of empirical potentials.

  5. Ab initio studies of isolated hydrogen vacancies in graphane

    Mapasha, R. E.; Molepo, M. P.; Chetty, N.

    2016-05-01

    We present a density functional study of various hydrogen vacancies located on a single hexagonal ring of graphane (fully hydrogenated graphene) considering the effects of charge states and the position of the Fermi level. We find that uncharged vacancies that lead to a carbon sublattice balance are energetically favorable and are wide band gap systems just like pristine graphane. Vacancies that do create a sublattice imbalance introduce spin polarized states into the band gap, and exhibit a half-metallic behavior with a magnetic moment of 1.00 μB per vacancy. The results show the possibility of using vacancies in graphane for novel spin-based applications. When charging such vacancy configurations, the deep donor (+1/0) and deep acceptor (0/-1) transition levels within the band gap are noted. We also note a half-metallic to metallic transition and a significant reduction of the induced magnetic moment due to both negative and positive charge doping.

  6. Hydrogen adsorption on boron doped graphene: an {\\it ab initio} study

    Miwa, R. H.; Martins, T B; Fazzio, A.

    2007-01-01

    The electronic and structural properties of (i) boron doped graphene sheets, and (ii) the chemisorption processes of hydrogen adatoms on the boron doped graphene sheets have been examined by {\\it ab initio} total energy calculations.

  7. Time-resolved photoelectron spectroscopy and ab initio multiple spawning studies of hexamethylcyclopentadiene

    Wolf, T. J. A.; Kuhlman, Thomas Scheby; Schalk, O.;

    2014-01-01

    Time-resolved photoelectron spectroscopy and ab initio multiple spawning were applied to the ultrafast non-adiabatic dynamics of hexamethylcyclopentadiene. The high level of agreement between experiment and theory associates wavepacket motion with a distinct degree of freedom....

  8. Ab initio study of the buckling on silicene and germanene

    Martinez-Guerra, Edgar; Hernández, Karla; Cifuentes-Quintal, Eduardo; de Coss, Romeo

    2013-03-01

    Recently, a new graphene-like silicon structure was discovered: silicene. Since its discovery, silicene has been more exciting than graphene because this is a semiconductor and it should be compatible with silicon-based electronic. Silicon and germanium atoms have similar electronic configurations as those of carbon and this the reason that the bandstructure of silicene and germacene exhibits the Dirac cones at K point, with a very similar linear dispersion around it, like in graphene. The disvintage is that sp2 bonded Si is much less stable than for carbon resulting that to be stable in the planar layer their atoms must buckle. In this work, we calculated the sp character on silicene and germacene to correlate its hibridization with the velocity of electrons and holes at Dirac cones. The calculations were performed using the pseudopotential LCAO method with GGA for the exchange-correlation energy functional. The buckling of silicene and germacene layer was 0.50 and 0.69 Å, respectively. In addition, the sp- character of silicene and germacene buckled was 2.33 and 2.64, respectively. Thus, a detailed analysis on the electronic band structure of these system show that as sp character goes from sp2 to sp3 it is correlated with a decrease of velocity of electrons and holes at Dirac cones. This study is primarly important and it could address a new future to modulate carrier velocities on bidimensional systems. This research was supported by Conacyt under Grant No. 133022.

  9. Ab initio study of boron nitride lines on graphene

    Mata-Carrizal, Berenice; Sanginés-Mendoza, Raúl; Martinez, Edgar

    2013-03-01

    Graphene has unusual electronic properties which make it a promising material for electronic devices. Neverthless, the absence of a band gap sets limitations on its practical applications. Thus, it is crucial to find methods to create and tune the band gap of systems based on graphene. In this way, we explore the modulation of the electronic properties of graphene through doping with boron nitride lines. In particular, we studied the electronic structure of graphene sheets doped with boron nitride lines armchair and zigzag type. The calculations were performed using the pseudopotential LCAO method with a Generalized Gradient Approximation (GGA) for the exchange-correlation energy functional. We found that both doping lines type induce a bandgap and that the energy gap increases as the length of doping lines increases. Accordingly to our DFT calculations, we found that the energy gap on graphene doped with armchair and zigzag lines is due to a two different mechanisms to drain charge from pi- to sigma- orbitals. Thus, we found that doping graphene with boron nitride lines is a useful way to induce and modulate the bandgap on graphene. This research was supported by Consejo Nacional de Ciencia y Tecnología (Conacyt) under Grant No. 133022.

  10. Ab-initio study of thermal expansion in pure graphene

    Mann, Sarita; Rani, Pooja; Kumar, Ranjan; Jindal, V. K.

    2016-05-01

    Graphene is a zero band gap semiconductor with exceptionally high thermal conductivity. The electronic properties having been studied, therole of phonon in contributing to thermal expansion, thermal conductivity and other thermodynamic properties, is required to be investigated. This paper focuses more on thermal expansion. Some others results like phonon dispersion, Grüneisenparameters and bulk modulus,which are essential to estimation of thermal expansion, are also presented. The dynamical matrix was calculated using VASP code using both DFT and DFPTand the phonon frequencies were calculated using phonopy code under harmonic approximation. The linear thermal expansion coefficient of graphene is found to be strongly dependent on temperature but remains negative upto 470 K and positive thereafter, with a room temperature value of -1.44×10-6. The negative expansion coefficient is very interesting and is found to be in conformity with experimental as well as with recent theoretical estimates. There is only qualitative agreement of our results with experimental data and motivates further investigation, primarily on the high negative values of Grüneisen parameters.

  11. Ab initio study of tungsten defects near the surface

    Guerrero, C. L.; Gordillo, N.; Iglesias, R.; Perlado, J. M.; Gonzalez, C.

    2016-05-01

    A first principles analysis of the behaviour of point defects, namely, self-interstitial atoms, a single vacancy and light impurity atoms such as H and He in tungsten is reported. These defects can be produced in the first wall of the future nuclear fusion reactors due to the high radiation fluxes present. The evolution of defects that appear in the bulk and end up reaching the surface has been followed. An energetic study has been combined with a detailed charge density analysis of the system by means of the SIESTA code. The resulting data have been validated by confronting them with those obtained with a more precise plane wave code, namely VASP. Meanwhile, the structural and the mechanical properties of the system have been positively compared with experimental measurements. Such comparisons have led us to present a new SIESTA basis for tungsten. This complete analysis establishes a nanoscopic view of the phenomena involving the presence of light atoms at native defects in tungsten, paying special attention to the vicinity of surfaces.

  12. Ab initio Study of Naptho-Homologated DNA Bases

    Sumpter, Bobby G [ORNL; Vazquez-Mayagoitia, Alvaro [ORNL; Huertas, Oscar [Universitat de Barcelona; Fuentes-Cabrera, Miguel A [ORNL; Orozco, Modesto [Institut de Recerca Biomedica, Parc Cientific de Barcelona, Barcelona, Spain; Luque, Javier [Universitat de Barcelona

    2008-01-01

    Naptho-homologated DNA bases have been recently used to build a new type of size expanded DNA known as yyDNA. We have used theoretical techniques to investigate the structure, tautomeric preferences, base-pairing ability, stacking interactions, and HOMO-LUMO gaps of the naptho-bases. The structure of these bases is found to be similar to that of the benzo-fused predecessors (y-bases) with respect to the planarity of the aromatic rings and amino groups. Tautomeric studies reveal that the canonical-like form of naptho-thymine (yyT) and naptho-adenine (yyA) are the most stable tautomers, leading to hydrogen-bonded dimers with the corresponding natural nucleobases that mimic the Watson-Crick pairing. However, the canonical-like species of naptho-guanine (yyG) and naptho-cytosine (yyC) are not the most stable tautomers, and the most favorable hydrogen-bonded dimers involve wobble-like pairings. The expanded size of the naphto-bases leads to stacking interactions notably larger than those found for the natural bases, and they should presumably play a dominant contribution in modulating the structure of yyDNA duplexes. Finally, the HOMO-LUMO gap of the naptho-bases is smaller than that of their benzo-base counterparts, indicating that size-expansion of DNA bases is an efficient way of reducing their HOMO-LUMO gap. These results are examined in light of the available experimental evidence reported for yyT and yyC.

  13. Multifunctional properties of zinc oxide: An ab-initio study

    Gopal, Priya

    In this dissertation, we present our results of first-principles calculations of the electronic and magnetic properties of a range of ZnO-based systems. ZnO is a II-VI wide band gap semiconductor, is a strong piezoelectric, is transparent and can be easily fabricated in a variety of microstructurures from bulk single crystals to thin films to nanostructures. Recent reports also suggest the possibility of magnetic ZnO as a potential spintronic material. Due to its range of properties, ZnO has numerous applications including varistors, piezoelectric transducers, sonars, acoustic devices, solar cells, UV blocking, chemical catalysts and transparent conducting oxides. In this work, we explore the following "multifunctional aspects" of ZnO-based systems. First we calculate the structural, polarization and piezoelectric properties of wurtzite-structure ZnO, MgO and CdO (used in band gap engineering of ZnO) to explore the relationship between optical and piezoelectric functionalities. Second, we investigate wurtzite-structure MnO to understand coexisting magnetic and piezoelectric functionalities. Finally, we study the nature of magnetic interactions in ZnO doped with 3d transition metals to explore the simultaneous magnetic, semiconducting and piezoelectric properties. Our main findings are: (i) Polarization differences between MgO and ZnO, and CdO and ZnO are large suggesting that strong polarization gradients will accompany the changes in band gap induced by Cd or Mg doping. (ii) wurtzite-structure MnO is antiferromagnetic, with a strong magnetostructural coupling and piezoelectric response; this indicates the possibility of coexistence of magnetism and piezoelectricity. (iii) In transition-metal doped ZnO we do not observe a room-temperature ferromagnetic interaction. We find that the magnetic ground state is highly sensitive to details of the calculations such as k-point mesh, details of the geometry optimizations and the choice of exchange correlation functional.

  14. An ab initio molecular dynamics study of the roaming mechanism of the H2+HOC+ reaction

    Yu, Hua-Gen

    2011-08-01

    We report here a direct ab initio molecular dynamics study of the p-/o-H2+HOC+ reaction on the basis of the accurate SAC-MP2 potential energy surface. The quasi-classical trajectory method was employed. This work largely focuses on the study of reaction mechanisms. A roaming mechanism was identified for this molecular ion-molecule reaction. The driving forces behind the roaming mechanism were thoroughly investigated by using a trajectory dynamics approach. In addition, the thermal rate coefficients of the H2+HOC+ reaction were calculated in the temperature range [25, 300] K and are in good agreement with experiments.

  15. Ab initio study of the EFG at the N sites in imidazole

    We study the nuclear quadrupole interaction at the nitrogen sites in the molecular and crystalline phases of the imidazole compound. We use PAW which is a state-of-the-art method to calculate the electronic structure and electric field gradient at the nucleus in the framework of the density functional theory. The quadrupole frequencies at both imino and amino N sites are in excellent agreement with measurements. This is the first time that the electric field gradient at crystalline imidazole is correctly treated by an ab initio theoretical approach.

  16. Ab-initio density functional theory study of a WO3 NH3-sensing mechanism

    WO3 bulk and various surfaces are studied by an ab-initio density functional theory technique. The band structures and electronic density states of WO3 bulk are investigated. The surface energies of different WO3 surfaces are compared and then the (002) surface with minimum energy is computed for its NH3 sensing mechanism which explains the results in the experiments. Three adsorption sites are considered. According to the comparisons of the energy and the charge change between before and after adsorption in the optimal adsorption site O1c, the NH3 sensing mechanism is obtained. (nuclear physics)

  17. Ab initio molecular dynamics simulation study of successive hydrogenation reactions of carbon monoxide producing methanol

    Pham, Thi Nu; Ono, Shota; Ohno, Kaoru

    2016-04-01

    Doing ab initio molecular dynamics simulations, we demonstrate a possibility of hydrogenation of carbon monoxide producing methanol step by step. At first, the hydrogen atom reacts with the carbon monoxide molecule at the excited state forming the formyl radical. Formaldehyde was formed after adding one more hydrogen atom to the system. Finally, absorption of two hydrogen atoms to formaldehyde produces methanol molecule. This study is performed by using the all-electron mixed basis approach based on the time dependent density functional theory within the adiabatic local density approximation for an electronic ground-state configuration and the one-shot GW approximation for an electronic excited state configuration.

  18. Ab initio simulations and neutron scattering studies of structure and dynamics in PdH

    The work presented in this PhD thesis is concerned with the interpretation of the neutron scattering measurements from the palladium hydrogen system by means of ab initio electronic structure calculations. The motivation of performing such calculations was due to recent neutron scattering studies on this system that showed a strong directional dependence to the dynamical structure factor together with a complex dependence on energy. Here we attempt to describe the origin of these features by ab initio simulations of the dynamical structure factor. The method assumes an adiabatic separation of the motion of the proton and palladium atoms. The proton wave functions are calculated by a direct solution of the associated single-particle Schroedinger equation using a plane wave basis set method and a mapping of the adiabatic surface. The Fourier components of the adiabatic potential are obtained from LDA pseudopotential calculations. Using Fermi's golden rule within the Born approximation we were then able to calculate the dynamical structure factor, S(Q,ω), for exciting the proton from its ground state to various excited states as a function of the magnitude and direction of the scattering wave vector. The results are in agreement with the inelastic neutron scattering spectra and allow us to identify the origin of previous inexplicable features, in particular the strong directional dependence to the experimental data. The method was extended to investigate the expansion of the equilibrium lattice constant as a function of the H isotope when the zero-point energy of the proton/deuterium is explicitly taken into account in the relaxation process. The results we obtained predicted a bigger lattice constant for the hydride, as expected. Furthermore, other complex ab initio calculations were carried out in order to describe the origin of the large optic dispersion, seen previously in the coherent neutron scattering data. Our calculated dispersion proved to be in good

  19. Ab Initio Study of Hot Carriers in the First Picosecond after Sunlight Absorption in Silicon

    Bernardi, Marco; Vigil-Fowler, Derek; Lischner, Johannes; Neaton, Jeffrey B.; Louie, Steven G.

    2014-06-01

    Hot carrier thermalization is a major source of efficiency loss in solar cells. Because of the subpicosecond time scale and complex physics involved, a microscopic characterization of hot carriers is challenging even for the simplest materials. We develop and apply an ab initio approach based on density functional theory and many-body perturbation theory to investigate hot carriers in semiconductors. Our calculations include electron-electron and electron-phonon interactions, and require no experimental input other than the structure of the material. We apply our approach to study the relaxation time and mean free path of hot carriers in Si, and map the band and k dependence of these quantities. We demonstrate that a hot carrier distribution characteristic of Si under solar illumination thermalizes within 350 fs, in excellent agreement with pump-probe experiments. Our work sheds light on the subpicosecond time scale after sunlight absorption in Si, and constitutes a first step towards ab initio quantification of hot carrier dynamics in materials.

  20. Ab initio, density functional theory and structural studies of 4-amino-2-methylquinoline

    Arjunan, V.; Saravanan, I.; Ravindran, P.; Mohan, S.

    2009-10-01

    The Fourier transform infrared (FTIR) and FT-Raman spectra of 4-amino-2-methylquinoline (AMQ) have been recorded in the range 4000-400 and 4000-100 cm -1, respectively. The experimental vibrational frequency was compared with the wavenumbers obtained theoretically by ab initio HF and DFT-B3LYP gradient calculations employing the standard 6-31G** and high level 6-311++G** basis sets for optimised geometry of the compound. The complete vibrational assignment and analysis of the fundamental modes of the compounds were carried out using the experimental FTIR and FT-Raman data, and quantum mechanical studies. The geometry and normal modes of vibration obtained from the HF and DFT methods are in good agreement with the experimental data. The potential energy distribution of the fundamental modes was calculated with ab initio force fields utilising Wilson's FG matrix method. The NH -π interactions and the influence of amino and methyl groups on the skeletal modes are investigated.

  1. Ab initio study of Ni2MnGa under shear deformation

    Zelený Martin

    2015-01-01

    Full Text Available The effect of shear deformation on Ni2MnGa magnetic shape memory alloy has been investigated using ab initio electronic structure calculations. We used the projector-augmented wave method for the calculations of total energies and stresses as functions of applied affine shear deformation. The studied nonmodulated martensite (NM phase exhibits a tetragonally distorted L21 structure with c/a > 1. A large strain corresponding to simple shears in {001}, {100} and {100} systems was applied to describe a full path between two equivalent NM lattices. We also studied {101} shear which is related to twining of NM phase. Twin reorientation in this system is possible, because applied positive shear results in path with significantly smaller energetic barrier than for negative shear and for shears in other studied systems. When the full relaxation of lattice parameters is allowed, the barriers further strongly decrease and the structures along the twinning path can be considered as orthorhombic.

  2. Comparative studies of the spectroscopy of CuCl2: DFT versus standard ab initio approaches.

    Ramírez-Solís, A; Poteau, R; Vela, A; Daudey, J P

    2005-04-22

    The X2Pi g-2Sigma g+, X2Pi g-2Delta g, X2Pi g-2Sigma u+, X2Pi g-2Pi u transitions on CuCl2 have been studied using several exchange-correlation functionals from the various types of density functional theory (DFT) approaches like local density approximation (LDA), generalized gradient approximation (GGA), hybrid and meta-GGA. The results are compared with the experience and with those coming from the most sophisticated nondynamic and dynamic electronic correlation treatments using the same relativistic effective core potentials and especially developed basis sets to study the electronic structure of the five lowest states and the corresponding vertical and adiabatic transition energies. The calculated transition energies for three of the hybrid functionals (B3LYP, B97-2, and PBE0) are in very good agreement with the benchmark ab initio results and experimental figures. All of the other functionals largely overestimate the X2Pi g-2Sigma g+ and X2Pi g-2Delta g transition energies, many of them even placing the 2Delta g ligand field state above the charge transfer 2Pi u and 2Sigma u+ states. The relative weight of the Hartree-Fock exchange in the definition of the functional used appears to play a key role in the accurate description of the LambdaSSigma density defined by the orientation of the 3d hole (sigma, pi, or delta) on Cu in the field of both chlorine atoms, but no simple connection of this weight with the quality of the spectra has been found. Mulliken charges and spin densities are carefully analyzed; a possible link between the extent of spin density on the metal for the X2Pi g state and the performance of the various functionals was observed, suggesting that those that lead to the largest values (close to 0.65) are the ones that best reproduce these four transitions. Most functionals lead to a remarkably low ionicity for the three ligand field states even for the best performing functionals, compared to the complete active space (SCF) (21, 14) ab initio

  3. Ab-initio study of magnetism behavior in TiO2 semiconductor with structural defects

    Zarhri, Z.; Houmad, M.; Ziat, Y.; El Rhazouani, O.; Slassi, A.; Benyoussef, A.; El Kenz, A.

    2016-05-01

    Magnetic, electronic and structural properties of titanium dioxide material with different structural defects are studied using the first-principles ab-initio calculations and the Korringa-Kohn-Rostoker method (KKR) combined with the coherent potential approximation (CPA) method in connection with the local density approximation (LDA). We investigated all structural defects in rutile TiO2 such as Titanium interstitial (Tii), Titanium anti-sites (Tio), Titanium vacancies (VTi), Oxygen interstitial (Oi), Oxygen anti-sites (OTi) and oxygen vacancies (Vo). Mechanisms of hybridization and interaction between magnetic atoms are investigated. The transition temperature is computed using the Mean Field Approximation (MFA).Magnetic stability energy of ferromagnetic and disordered local moment states is calculated to determine the most stable state. Titanium anti-sites have a half-metallic aspect. We also studied the change type caused by structural defects in this material.

  4. Tunneling of electrons via rotor-stator molecular interfaces: combined ab initio and model study

    Petreska, Irina; Pejov, Ljupco; Kocarev, Ljupco

    2015-01-01

    Tunneling of electrons through rotor-stator anthracene aldehyde molecular interfaces is studied with a combined ab initio and model approach. Molecular electronic structure calculated from first principles is utilized to model different shapes of tunneling barriers. Together with a rectangular barrier, we also consider a sinusoidal shape that captures the effects of the molecular internal structure more realistically. Quasiclassical approach with the Simmons' formula for current density is implemented. Special attention is paid on conformational dependence of the tunneling current. Our results confirm that the presence of the side aldehyde group enhances the interesting electronic properties of the pure anthracene molecule, making it a bistable system with geometry dependent transport properties. We also investigate the transition voltage and we show that confirmation dependent field emission could be observed in these molecular interfaces at realistically low voltages. The present study accompanies our previ...

  5. Hydrogen-Water Mixtures in Giant Planet Interiors Studied with Ab Initio Simulations

    Soubiran, Francois

    2015-01-01

    We study water-hydrogen mixtures under planetary interior conditions using ab initio molecular dynamics simulations. We determine the thermodynamic properties of various water-hydrogen mixing ratios at temperatures of 2000 and 6000 K for pressures of a few tens of GPa. These conditions are relevant for ice giant planets and for the outer envelope of the gas giants. We find that at 2000 K the mixture is in a molecular regime, while at 6000 K the dissociation of hydrogen and water is important and affects the thermodynamic properties. We study the structure of the liquid and analyze the radial distribution function. We provide estimates for the transport properties, diffusion and viscosity, based on autocorrelation functions. We obtained viscosity estimates of the order of a few tenths of mPa.s for the conditions under consideration. These results are relevant for dynamo simulations of ice giant planets.

  6. Atomic carbon chains as spin-transmitters: An ab initio transport study

    Fürst, Joachim Alexander; Brandbyge, Mads; Jauho, Antti-Pekka

    2010-01-01

    An atomic carbon chain joining two graphene flakes was recently realized in a ground-breaking experiment by Jin et al. (Phys. Rev. Lett., 102 (2009) 205501). We present ab initio results for the electron transport properties of such chains and demonstrate complete spin-polarization of the transmi...

  7. Ab initio study of long-period superstructures in close-packed A3B compounds

    Rosengaard, N. M.; Skriver, Hans Lomholt

    1994-01-01

    We have performed ab initio calculations of the stability of one-dimensional long-period superstructures in Cu3Pd, Cu3Al, and Ag3Mg by means of an interface Green's function technique based on the linear-muffin-tin-orbitals method within the tight-binding and atomic-sphere approximations. The...

  8. Study on the surface hydroxyl group on solid breeding materials by ab-initio calculations

    Tanaka, Satoru; Taniguchi, Masaki [Tokyo Univ. (Japan). Faculty of Engineering

    1996-10-01

    The nature of -OH on the surface of Li{sub 2}O was analyzed with the ab-initio quantum chemical calculation technique. Calculation results showed that the stretching vibration of O-H is affected by the chemical species around the -OH. (author)

  9. Raman and ab initio studies of simple and binary 1-alkyl-3-methylimidazolium ionic liquids

    Berg, R.W.; Deetlefs, M.; Seddon, K.R.;

    2005-01-01

    Raman spectra of the ionic liquids, 1-butyl-3-methylimidazolium hexafluorophosphate ([C(4)mim][PF6]), 1-hexyl-3-methylimidazolium chloride ([C(6)mim]Cl), and 1-hexyl-3-methylimidazolium hexafluorophosphate ([C(6)mim][PF6]), and binary mixtures thereof, have been assigned using ab initio MP2...

  10. Structural and electronic properties of lead nanowires: Ab-initio study

    Highlights: → In the present revised manuscript entitled 'Structural and Electronic Properties of Lead Nanowires: Ab-initio study', we have analyzed the stability, electronic properties as well as ground state properties of various atomic configurations of Lead nanowires. → The two-atom zigzag shaped lead nanowire with highest binding energy and lowest total energy has been confirmed as the most stable structure out of the six atomic configurations taken into consideration. → The electronic band structure and density of states have been described in detail with a remarkable observation in case of three-atom triangular lead nanowire having a very small band gap while other atomic configurations are found to be metallic. → The bulk modulus and pressure derivatives for all the stable geometries have also been computed and discussed in the manuscript. The mechanical strength of nanowires has also been discussed in terms of its bulk modulus. → The two-atom ladder shaped nanowire with highest bulk modulus, defends this structure as mechanically stronger than the other tested structure. - Abstract: Ab-initio self-consistent study has been performed to analyze the stability of lead nanowires in its six stable configurations like linear, zigzag, triangular, ladder, square and dumbbell. In the present study, the lowest energy structures have been analyzed under the revised Perdew-Burke-Ernzerhof (revPBE) parameterization of generalized gradient approximation (GGA) potential. The two-atom zigzag shaped atomic configuration with highest binding energy and lowest total energy has been confirmed as the most stable structure out of the six atomic configurations. The electronic band structure and density of states have been discussed in detail with a remarkable observation in case of three-atom triangular lead nanowire having a very small band gap while other configurations are found to be metallic. Bulk modulus, pressure derivatives and lattice parameters for different

  11. Structural properties of iron nitride on Cu(100): An ab-initio molecular dynamics study

    Heryadi, Dodi

    2011-01-01

    Due to their potential applications in magnetic storage devices, iron nitrides have been a subject of numerous experimental and theoretical investigations. Thin films of iron nitride have been successfully grown on different substrates. To study the structural properties of a single monolayer film of FeN we have performed an ab-initio molecular dynamics simulation of its formation on a Cu(100) substrate. The iron nitride layer formed in our simulation shows a p4gm(2x2) reconstructed surface, in agreement with experimental results. In addition to its structural properties, we are also able to determine the magnetization of this thin film. Our results show that one monolayer of iron nitride on Cu(100) is ferromagnetic with a magnetic moment of 1.67 μ B. © 2011 Materials Research Society.

  12. Ab Initio Study of BiFeO3: Thermodynamic Stability Conditions.

    Heifets, E; Kotomin, E A; Bagaturyants, A A; Maier, J

    2015-07-16

    BiFeO3 is investigated intensively, mainly as a multiferroic material. In this paper, the state-of-the-art ab initio hybrid functional approach with atomic basis sets was employed for a study of the stability range of BiFeO3 with respect to its decomposition into binary oxides and elementary metals, as a function of temperature and oxygen partial pressure. The calculated atomic and electronic structure of BiFeO3 was compared with previous LDA+U calculations using plane-wave basis sets. Based on performed calculations, the phase diagram was constructed, which allows us to predict the stability region of stoichiometric BiFeO3. PMID:26266871

  13. Ab Initio Studies on Hematite Surface and the Adsorption of Phosphate

    Swati Chaudhury

    2014-01-01

    Full Text Available This investigation explores the ab initio DFT method for understanding surface structure of hematite and the nature and energetics of phosphate adsorption. Using the full potential linearized plane wave method (FP-LAPW, we derived the structure and energies of various magnetic forms of hematite. The antiferromagnetic (AFM form was observed to be the most stable. Hematite surfaces with Fe-termination, O-termination, or OH-termination were studied. The OH-terminated surface was the most stable. Stability of hematite surfaces follows the order OH-termination > Fe-termination > O-termination. Thus, surface reaction with hematite would occur with the OH at the surface and not with Fe atoms. The structure of phosphate adsorbed on hematite was derived. Bonding is through the H atom of the OH at the surface. An alternative mechanism of phosphate adsorption on hematite has been derived. Adsorption energy is high and suggests chemisorption rather than physisorption of phosphate on hematite.

  14. An ab initio cluster model study of the magnetic coupling in KNiF3

    Casanovas, Jordi; Illas, Francesc

    1994-06-01

    Cluster models of increasing complexity have been used to model magnetic interactions in KNiF3. These clusters contain two or four magnetic centers plus the bridge F- anions and different representations of the remaining of the crystal. The magnetic coupling constant has been obtained by computing ab initio wave functions for different spin states. These wave functions explicitly include internal and external correlation effects. Several sets of Gaussian functions have been tested and many sets of molecular orbitals have been considered in order to study the physical origin of magnetism in KNiF3. The calculated magnetic coupling constant differs from model to model but shows a fairly good convergence to the experimental result. The use of different cluster models permits to separate the magnetic coupling constant in several contributions. These are the delocalization of magnetic orbitals, the external correlation, and the collective effects normally hidden in the two body operator of the Heisenberg Hamiltonian.

  15. Ab initio study of one-dimensional disorder on III-V semiconductor surfaces

    Romanyuk, O. [Institute of Physics, Academy of Sciences of the Czech Republic, Cukrovarnicka 10, 162 00 Prague 6 (Czech Republic); Grosse, F.; Braun, W. [Paul-Drude-Institut fuer Festkoerperelektronik, Hausvogteiplatz 5-7, 10117 Berlin (Germany)

    2010-02-15

    Atomic disorder on GaSb(001) and GaAs(001) surfaces is studied by ab initio calculations within density functional theory (DFT). Surface energies are computed for GaSb(001) {beta} (4 x 3), {beta} (4 x 3)A{sub 1}, (4 x 6), and GaAs(001) {beta} 2(2 x 4), c (4 x 4) reconstructions. The computed energy differences do not exceed 1 meV / (1 x 1) which is in agreement with observed one-dimensional disorder on GaSb(001) and GaAs(001) surfaces at elevated temperatures for {beta} and {beta} 2 stoichiometries, respectively. Deviations in bond lengths due to disorder with respect to the ordered ground state phases are calculated (copyright 2010 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim) (orig.)

  16. An ab initio study of three (ethane-1,2 diol/water) complexes

    Manivet, Philippe; Masella, Michel

    1998-05-01

    Three (ethane-1,2 diol/water) complexes have been studied using ab initio calculations at the MP2 level. In two complexes, the ethane-1,2 diol structure is close to its gas phase experimental structure (presence of an intramolecular hydrogen bond HB and the O-C-C-O dihedral angle is gauche) while the intramolecular HB is disrupted by the presence of a water molecule in the third ( tGg'a). Computations have shown that most of the experimental observations regarding the solvation of ethane-1,2 diol in water may be reproduced only by considering the tGg'a complex (absence of intramolecular HB, O-C-C-O dihedral angle of 72-74°), which is also more stable than the other two by 2 kcal mol -1.

  17. Study of the atomic rare gas behaviour by ab initio calculations

    The atomic behaviour of helium and krypton in uranium dioxide has been studied using an ab initio simulation technique. Incorporation energies and solution energies of these two rare gases have been calculated. Krypton atoms are found to be insoluble in this nuclear fuel whatever the trap considered and their presence in the lattice induces swelling when they are located in interstitial sites or in oxygen vacancies. Due to its small atomic size, the predicted helium behaviour is very different. Indeed, helium is found to be soluble in stoichiometric and hyper-stoichiometric uranium dioxide in the presence of uranium vacancies or divacancies constituted by one uranium and one oxygen vacancy. Moreover helium atoms induce a lattice parameter contraction except in interstitial sites where a slight expansion is calculated. (author)

  18. High multiplicity states in disordered carbon systems: Ab initio and semiempirical study

    Stability of non-zero spin projection states for disordered carbon clusters of low symmetry were examined using semiempirical and ab initio methods. The study proves previous results of V.D. Khavryuchenko, Y.A. Tarasenko, V.V. Strelko, O.V. Khavryuchenko, V.V. Lisnyak, Int. J. Mod. Phys. B 21 (2007) 4507, obtained for the large polyaromatic hydrocarbons clusters and shows that the phenomenon is intrinsic for carbon-rich systems and independent of their symmetries. The electronic properties of the carbon clusters may alter from insulating to semiconducting upon change of C/H ratio and stabilization of non-zero spin projection states. A partial collectivization of the electrons is observed in deeply carbonized carbon clusters in higher Sz states.

  19. High-pressure study of ScVO4 by Raman scattering and ab initio calculations

    Panchal, V.; Manjón, F. J.; Errandonea, D.; Rodriguez-Hernandez, P.; López-Solano, J.; Muñoz, A.; Achary, S. N.; Tyagi, A. K.

    2011-02-01

    We report results of experimental and theoretical lattice-dynamics studies on scandium orthovanadate up to 35 GPa. Raman-active modes of the low-pressure zircon phase are measured up to 8.2 GPa, where the onset of an irreversible zircon-to-scheelite phase transition is detected. Raman-active modes in the scheelite structure are observed up to 16.5 GPa. Beyond 18.2 GPa we detected a gradual splitting of the Eg modes of the scheelite phase, indicating the onset of a second phase transition. Raman symmetries, frequencies, and pressure coefficients in the three phases of ScVO4 are discussed in the light of ab initio lattice-dynamics calculations that support the experimental results. The results on all the three phases of ScVO4 are compared with those previously reported for related orthovanadates.

  20. Femtosecond Laser Processing of Germanium: An Ab Initio Molecular Dynamics Study

    Ji, Pengfei

    2016-01-01

    An ab initio molecular dynamics study of femtosecond laser processing of germanium is presented in this paper. The method based on the finite temperature density functional theory is adopted to probe the structural change, thermal motion of the atoms, dynamic property of the velocity autocorrelation, and the vibrational density of states. Starting from a cubic system at room temperature (300 K) containing 64 germanium atoms with an ordered arrangement of 1.132 nm in each dimension, the femtosecond laser processing is simulated by imposing the Nose Hoover thermostat to the electronic subsystem lasting for ~100 fs and continuing with microcanonical ensemble simulation of ~200 fs. The simulation results show solid, liquid and gas phases of germanium under adjusted intensities of the femtosecond laser irradiation. We find the irradiated germanium distinguishes from the usual germanium crystal by analyzing their melting and dynamic properties.

  1. Manipulating magnetism and conductance of an adatom-molecule junction on a metal surface: An ab initio study

    Tao, Kun; Stepanyuk, V.S.; Bruno, P.;

    2008-01-01

    The state of the art ab initio calculations reveal the effect of a scanning tunneling microscopy tip on magnetic properties and conductance of a benzene-adatom sandwich on Cu(001). We concentrate on a benzene-Co system interacting with a Cr tip. Our studies give a clear evidence that magnetism an...

  2. A Photoelectron Spectroscopy and ab initio Study of B3- and B4- Anions and Their Neutrals

    Zhai, Hua-Jin; Wang, Lai S.; Alexandrova, A N.; Boldyrev, Alexander I.; Zakrzewski, V G.

    2003-11-06

    The two smallest boron clusters (B3 and B4) in their neutral and anionic forms were studied by photoelectron spectroscopy and ab initio calculations. Vibrationally resolved photoelectron spectra were observed for B3- at three photon energies (355, 266, and 193 nm) and the electron affinity of B3 was measured to be+0.02 eV.

  3. Ab-initio study of magnetic properties and phase transitions in Ga (Mn) N with Monte Carlo approach

    Sbai, Y.; Ait Raiss, A.; Salmani, E. [LMPHE (URAC 12), Faculty of Science, Mohammed V University, Av. Ibn Batouta, Rabat (Morocco); Bahmad, L., E-mail: Bahmad@fsr.ac.ma [LMPHE (URAC 12), Faculty of Science, Mohammed V University, Av. Ibn Batouta, Rabat (Morocco); Benyoussef, A. [LMPHE (URAC 12), Faculty of Science, Mohammed V University, Av. Ibn Batouta, Rabat (Morocco); Institute of Nanomaterials and Nanotechnologies, MAScIR, Rabat (Morocco); Hassan II Academy of Science and Technology, Rabat (Morocco)

    2015-12-15

    On the basis of ab-initio calculations and Monte Carlo simulations the magnetic and electronic properties of Gallium nitride (GaN) doped with the transition metal Manganese (Mn) were studied. The ab initio calculations were done using the AKAI–KKR–CPA method within the Local Density Approximation (LDA) approximation. We doped our Diluted Magnetic Semiconductor (DMS), with different concentrations of magnetic impurities Mn and plotted the density of state (DOS) for each one. Showing a half-metallic behavior and ferromagnetic state especially for Ga{sub 0.95}Mn{sub 0.05}N making this DMS a strong candidate for spintronic applications. Moreover, the magnetization and susceptibility of our system as a function of the temperature has been calculated and give for various system size L to study the size effect. In addition, the transition temperature was deduced from the peak of the susceptibility. The Ab initio results are in good agreement with literature especially for (x=0.05) of Mn which gives the most interesting results. - Highlights: • The AKAI–KKR–CPA method has been applied to study the doped compound GaN:Mn. • The local density approximation (LDA) has been applied. • The ab-initio calculations have been performed. • The density of states (DOS) have been plotted for differents doping concentrations, using Monte Carlo simulations.

  4. Conformational space of clindamycin studied by ab initio and full-atom molecular dynamics.

    Kulczycka-Mierzejewska, Katarzyna; Trylska, Joanna; Sadlej, Joanna

    2016-01-01

    Molecular dynamics (MD) simulations allow determining internal flexibility of molecules at atomic level. Using ab initio Born-Oppenheimer molecular dynamics (BOMD), one can simulate in a reasonable time frame small systems with hundreds of atoms, usually in vacuum. With quantum mechanics/molecular mechanics (QM/MM) or full-atom molecular dynamics (FAMD), the influence of the environment can also be simulated. Here, we compare three types of MD calculations: ab initio BOMD, hybrid QM/MM, and classical FAMD. As a model system, we use a small antibiotic molecule, clindamycin, which is one of the lincosamide antibiotics. Clindamycin acquires two energetically stable forms and we investigated the transition between these two experimentally known conformers. We performed 60-ps BOMD simulations in vacuum, 50-ps QM/MM, and 100-ns FAMD in explicit water. The transition between two antibiotic conformers was observed using both BOMD and FAMD methods but was not noted in the QM/MM simulations. PMID:26733483

  5. Melting of sodium under high pressure. An ab-initio study

    González, D. J.; González, L. E. [Departamento de Física Teórica, Atómica y Optica, Universidad de Valladolid, 47011 Valladolid (Spain)

    2015-08-17

    We report ab-initio molecular dynamics simulations of dense liquid/solid sodium for a pressure range from 0 to 100 GPa. The simulations have been performed with the orbital free ab-initio molecular dynamics method which, by using the electron density as the basic variable, allows to perform simulations with large samples and for long runs. The calculated melting curve shows a maximum at a pressure ≈ 30 GPa and it is followed by a long, steep decrease. These features are in good agreement with the experimental data. For various pressures along the melting curve, we have calculated several liquid static properties (pair distribution functions, static structure factors and short-range order parameters) in order to analyze the structural effects of pressure.

  6. Experimental studies and ab initio calculations on characteristics of the C state of SF2 radical

    SF2 radicals were generated by a pulsed dc discharge in the mixture gas beam of SF2 and Ar. The (2+1) resonance-enhanced multiphoton ionization (REMPI) excitation spectroscopy of SF2 radical was obtained between 325 and 365 nm. The SF+ ion signals were also observed in the same wavelength range. The analysis shows that the spectrum can be assigned as the two-photon resonant excitation of SF2 radical (B-tilde1 B1 and (C-tilde1 A1 states). And also, ab initio calculations suggest that the C-tilde state is a bonding state with Rydberg characteristic. The potential energy surfaces (PESs) of SF2 and SF2+ by ab initio calculations suggest that SF+ ions originate from dissociation processes of excited SF2+ ions. (author)

  7. Ab InitioStudy of Hot Carriers in the First Picosecond after Sunlight Absorption in Silicon

    Bernardi, Marco; Vigil-Fowler, Derek; Lischner, Johannes; Neaton, Jeffrey B.; Louie, Steven G.

    2014-01-01

    Hot carrier thermalization is a major source of efficiency loss in solar cells. Because of the subpicosecond time scale and complex physics involved, a microscopic characterization of hot carriers is challenging even for the simplest materials. We develop and apply an ab initio approach based on density functional theory and many-body perturbation theory to investigate hot carriers in semiconductors. Our calculations include electron-electron and electron-phonon interactions, and require no e...

  8. How Is Acetylcholinesterase Phosphonylated by Soman? An Ab Initio QM/MM Molecular Dynamics Study

    Sirin, Gulseher Sarah; Zhang, Yingkai

    2014-01-01

    Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nerve system that hydrolyzes acetylcholine (ACh) and terminates synaptic signals by reducing the effective concentration of ACh in the synaptic clefts. Organophosphate compounds irreversibly inhibit AChEs, leading to irreparable damage to nerve cells. By employing Born–Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characteriz...

  9. Ab initio study of one-dimensional disorder on III-V semiconductor surfaces

    Romanyuk, Olexandr; Grosse, F.; Braun, W.

    2010-01-01

    Roč. 7, č. 2 (2010), s. 330-333. ISSN 1862-6351 R&D Projects: GA AV ČR KAN300100802 Grant ostatní: German Resarch Fondation(DE) GZ:436 TSE 113/62/0-1 Institutional research plan: CEZ:AV0Z10100521 Keywords : atomic disorder * ab initio * semiconductor * reconstruction Subject RIV: BM - Solid Matter Physics ; Magnetism

  10. Ab initio study of antiphase boundaries and stacking faults in L12 and DO22 compounds

    Rosengaard, N. M.; Skriver, Hans Lomholt

    1994-01-01

    We have performed ab initio calculations of the energies of antiphase boundaries as well as complex and superlattice intrinsic stacking faults in nine intermetallic compounds observed in the face-centered-cubic L1(2) and DO22 structures. The calculations were performed by means of a Green's funct......'s function technique based on the linear-muffin-tin-orbitals method within the tight-binding and atomic-sphere approximations....

  11. Ab initio study of phase transformations in transition-metal disilicides

    Káňa, Tomáš; Šob, Mojmír; Vitek, V.

    2011-01-01

    Roč. 19, č. 7 (2011), s. 919-926. ISSN 0966-9795 R&D Projects: GA AV ČR IAA100100920; GA MŠk(CZ) OC10008 Institutional research plan: CEZ:AV0Z20410507 Keywords : silicides various * phase transformation * plastic deformation mechanisms * defects * dislocation geometry and arrangement * ab-initio calculations * aero- engine components Subject RIV: BM - Solid Matter Physics ; Magnetism Impact factor: 1.649, year: 2011

  12. Ab initio study of ideal tensile strength and mechanical stability of transition-metal disilicides

    Friák, Martin; Šob, Mojmír; Vitek, V.

    2003-01-01

    Roč. 68, č. 18 (2003), s. 184101-1 - 181101-10. ISSN 0163-1829 R&D Projects: GA AV ČR IAA1041302; GA ČR GA202/03/1351; GA MŠk OC 523.90 Institutional research plan: CEZ:AV0Z2041904 Keywords : ab initio calculations * electronic structure * theoretical tensile strength Subject RIV: BM - Solid Matter Physics ; Magnetism Impact factor: 2.962, year: 2003

  13. Capillary affinity electrophoresis and ab initio calculation studies on complexation of valinomycin with Na+ ion

    Ehala, Sille; Dybal, Jiří; Makrlík, E.; Kašička, Václav

    2009-01-01

    Roč. 32, č. 4 (2009), s. 597-604. ISSN 1615-9306 R&D Projects: GA ČR(CZ) GA203/06/1044; GA ČR(CZ) GA203/08/1428; GA AV ČR 1ET400500402 Institutional research plan: CEZ:AV0Z40550506; CEZ:AV0Z40500505 Keywords : capillary affinity electrophoresis * valinomycin * ab initio calculation Subject RIV: CB - Analytical Chemistry, Separation Impact factor: 2.551, year: 2009

  14. Ab initio study of thermodynamic, electronic, magnetic, structural, and elastic properties of Ni4N allotropes

    Hemzalová, P.; Friák, Martin; Šob, Mojmír; Ma, D.; Udyansky, A.; Raabe, D.; Neugebauer, J.

    2013-01-01

    Roč. 88, č. 17 (2013), Art. no. 174103. ISSN 1098-0121 R&D Projects: GA ČR(CZ) GAP108/12/0311; GA ČR GD106/09/H035; GA AV ČR IAA100100920 Institutional support: RVO:68081723 Keywords : nitrides * ab initio * thermodynamics * elasticity Subject RIV: BM - Solid Matter Physics ; Magnetism Impact factor: 3.664, year: 2013

  15. Ab initio studies of the kinetic isotope effect of the CH4 + OH atmospheric reaction

    Lasaga, Antonio C.; Gibbs, G. V.

    1991-07-01

    High level ab initio calculations have been carried out on the C-13 - C-12 kinetic isotope effect of the CH4 + OH reaction in the atmosphere. The results agree quite well with both the absolute value of the isotope effect and the temperature dependence of the effect, based on new experimental data. The calculated kinetic isotope effect supports a bigger effect of biomass burning on the methane global budget.

  16. Ab initio studies on the lattice thermal conductivity of silicon clathrate frameworks II and VIII

    Härkönen, Ville J.; Karttunen, Antti J.

    2015-01-01

    The lattice thermal conductivities of silicon clathrate frameworks II and VIII are investigated by using ab initio lattice dynamics and iterative solution of the linearized Boltzmann transport equation(BTE) for phonons. Within the temperature range 100-350 K, the clathrate structures II and VIII were found to have lower lattice thermal conductivity values than silicon diamond structure (d-Si) by factors of 1/2 and 1/5, respectively. The main reason for the lower lattice thermal conductivity o...

  17. Ab initio molecular dynamics study of liquid sodium and cesium up to critical point

    Yuryev, Anatoly A. [Institute of Metallurgy of Ural Branch of the Russian Academy of Sciences, Amundsen st. 101,620016, Yekaterinburg (Russian Federation); Ural Federal University, Vira st. 19, 620002, Yekaterinburg (Russian Federation); Gelchinski, Boris R. [Institute of Metallurgy of Ural Branch of the Russian Academy of Sciences, Amundsen st. 101,620016, Yekaterinburg (Russian Federation)

    2015-08-17

    Ab initio modeling of liquid metals Na and K is carried out using the program SIESTA. We have determined the parameters of the model (the optimal step, the number of particles, the initial state etc) and calculated a wide range of properties: the total energy, pair correlation function, coefficient of self-diffusion, heat capacity, statistics of Voronoi polyhedra, the density of electronic states up to the critical temperature.

  18. Discovering chemistry with an ab initio nanoreactor

    Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.

    2014-01-01

    Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provides detailed physical insight. While theoretical and computational studies have generally focused on specific processes or mechanistic hypotheses, recent methodological and computational advances harken the advent of their principal role in discovery. Here we report the development and application of the ab initio nanoreactor – a highly accelerat...

  19. On Possible Pitfalls in ab initio QM/MM Minimization Approaches For Studies of Enzymatic Reactions

    Klähn, Marco; Braun-Sand, Sonja; Rosta, Edina; Warshel, Arieh

    2006-01-01

    Reliable studies of enzymatic reactions by combined quantum mechanics /molecular mechanics (QM/MM) approaches, with an ab initio description of the quantum region, presents a major challenge to computational chemists. The main problem is the need for a very large computer time for the evaluation of the QM energy, which in turn makes it extremely challenging to perform proper configurational sampling. A seemingly reasonable alternative is to perform energy minimization studies of the type used in gas phase ab initio studies. However, it is hard to see why such an approach should give reliable results in protein active sites. In order to examine the problems with energy minimization QM/MM approaches we chose the hypothetical reaction of a metaphosphate ion with water in the Ras•GAP complex. This hypothetical reaction served as a simple benchmark reaction. The possible problems with the QM/MM minimization were explored by generating several protein configurations from long MD simulations and using energy minimization and scanning of the reaction coordinates to evaluate the corresponding potential energy surfaces of the reaction for each of these different protein configurations. Comparing these potential energy surfaces, we found major variations of the minima of the different total potential energy surfaces. Furthermore, the reaction energies and activation energies also varied significantly even for similar protein configurations. The specific coordination of a magnesium ion, present in the active center of the protein complex, turned out to influence the energetics of the reaction in a major way and a direct coordination to the reactant leads to an increase of the activation energy by 17 kcal/mol. This study demonstrates that energy minimizations starting from a single protein structure could lead to major errors in calculations of activation free energies and binding free energies. Thus we believe that extensive samplings of the configurational space of the

  20. Ab-Initio Molecular Dynamics

    Kühne, Thomas D

    2012-01-01

    Computer simulations and molecular dynamics in particular, is a very powerful method to provide detailed and essentially exact informations of classical many-body problems. With the advent of \\textit{ab-initio} molecular dynamics, where the forces are computed on-the-fly by accurate electronic structure calculations, the scope of either method has been greatly extended. This new approach, which unifies Newton's and Schr\\"odinger's equations, allows for complex simulations without relying on any adjustable parameter. This review is intended to outline the basic principles as well as a survey of the field. Beginning with the derivation of Born-Oppenheimer molecular dynamics, the Car-Parrinello method as well as novel hybrid scheme that unifies best of either approach are discussed. The predictive power is demonstrated by a series of applications ranging from insulators to semiconductors and even metals in condensed phases.

  1. The crystallization process of liquid vanadium studied by ab initio molecular dynamics

    We present a study of the crystallization process in liquid vanadium over a temperature range from 3000 K down to 1500 K using ab initio molecular dynamics simulations. Short-range order evolution during solidification is studied using various structural analysis methods. We show that the icosahedral-like short-range order is detected in the stable liquid phase and grows upon supercooling. The system undergoes a first-order phase transition (from a liquid to a solid state) at a temperature of about 1600 K. The crystal nucleation process is further studied using the time–temperature transformation mechanism by annealing the system at 1650 K. The nucleation is examined using bond-orientational order and density fluctuation analysis. Our finding is that various precursors appear in the region of high bond-orientational order with the majority having body-centered cubic (bcc)-like symmetry. This bcc-like region grows on annealing via thermal fluctuations. Our results reveal that the bond-orientational order precedes the density fluctuation, and is the main driving factor for nucleation. (papers)

  2. Ab initio study of basal slip in Nb2AlC

    Using ab initio calculations, we have studied shearing in Nb2AlC, where NbC and Al layers are interleaved. The stress-strain analysis of this deformation mode reveals Nb-Al bond breaking, while the Nb-C bond length decreases by 4.1%. Furthermore, there is no evidence for phase transformation during deformation. This is consistent with basal slip and may be understood on the basis of the electronic structure: bands below the Fermi level are responsible for the dd bonding between NbC basal planes and only a single band with a weak dd interaction is not resistant to shearing, while all other bands are unaffected. The Al-Nb bonding character can be described as mainly metallic with weak covalent-ionic contributions. Our study demonstrates that Al layers move with relative ease under shear strain. Phase conservation upon shearing is unusual for carbides and may be due to the layered nature of the phase studied. Here, we describe the electronic origin of basal slip in Nb2AlC, the atomic mechanism which enables reversible plasticity in this class of materials

  3. Modulation of electronic properties of silicon carbide nanotubes via sulphur-doping: An ab initio study

    Singh, Ram Sevak; Solanki, Ankit

    2016-03-01

    Silicon carbide nanotubes (SiCNTs) have received a great deal of scientific and commercial interest due to their intriguing properties that include high temperature stability and electronic properties. For their efficient and widespread applications, tuning of electronic properties of SiCNTs is an attractive study. In this article, electronic properties of sulphur doped (S-doped) zigzag (9 , 0) SiCNT is investigated by ab initio calculations based on density functional theory (DFT). Energy band structures and density of states of fully optimized undoped and doped structures with varying dopant concentration are calculated. S-doped on C-site of the nanotube exhibits a monotonic reduction of energy gap with increase in dopant concentration, and the nanotube transforms from semiconductor to metal at high dopant concentration. In case of S-doped on Si-site doping has less influence on modulating electronic structures, which results in reduction of energy gap up to a moderate doping concentration. Importantly, S preferential substitutes of Si-sites and the nanotube with S-doped on Si-site are energetically more stable as compared to the nanotube with S-doped on C-site. The study of tunable electronic properties in S-doped SiCNT may have potential in fabricating nanoelectronic devices, hydrogen storage and gas sensing applications.

  4. Tunneling of electrons via rotor-stator molecular interfaces: Combined ab initio and model study

    Petreska, Irina; Ohanesjan, Vladimir; Pejov, Ljupčo; Kocarev, Ljupčo

    2016-07-01

    Tunneling of electrons through rotor-stator anthracene aldehyde molecular interfaces is studied with a combined ab initio and model approach. Molecular electronic structure calculated from first principles is utilized to model different shapes of tunneling barriers. Together with a rectangular barrier, we also consider a sinusoidal shape that captures the effects of the molecular internal structure more realistically. Quasiclassical approach with the Simmons' formula for current density is implemented. Special attention is paid on conformational dependence of the tunneling current. Our results confirm that the presence of the side aldehyde group enhances the interesting electronic properties of the pure anthracene molecule, making it a bistable system with geometry dependent transport properties. We also investigate the transition voltage and we show that conformation-dependent field emission could be observed in these molecular interfaces at realistically low voltages. The present study accompanies our previous work where we investigated the coherent transport via strongly coupled delocalized orbital by application of Non-equilibrium Green's Function Formalism.

  5. Elastic properties of amorphous boron suboxide based solids studied using ab initio molecular dynamics

    We have studied the correlation between chemical composition, structure, chemical bonding and elastic properties of amorphous B6O based solids using ab initio molecular dynamics. These solids are of different chemical compositions, but the elasticity data appear to be a function of density. This is in agreement with previous experimental observations. As the density increases from 1.64 to 2.38 g cm-3, the elastic modulus increases from 74 to 253 GPa. This may be understood by analyzing the cohesive energy and the chemical bonding of these compounds. The cohesive energy decreases from -7.051 to -7.584 eV/atom in the elastic modulus range studied. On the basis of the electron density distributions, Mulliken analysis and radial distribution functions, icosahedral bonding is the dominating bonding type. C and N promote cross-linking of icosahedra and thus increase the density, while H hinders the cross-linking by forming OH groups. The presence of icosahedral bonding is independent of the density

  6. Numerical study of two-photon ionization of helium using an ab initio numerical framework

    Few-photon-induced breakup of helium is studied using a newly developed ab initio numerical framework for solving the six-dimensional time-dependent Schroedinger equation. We present details of the method and calculate (generalized) cross sections for the process of two-photon nonsequential (direct) double ionization at photon energies ranging from 39.4 to 54.4 eV, a process that has been very much debated in recent years and is not yet fully understood. In particular, we have studied the convergence property of the total cross section in the vicinity of the upper threshold (∼ 54.4 eV) versus the pulse duration of the applied laser field. We find that the cross section exhibits an increasing trend near the threshold, as has also been observed by others, and show that this rise cannot solely be attributed to an unintended inclusion of the sequential two-photon double ionization process caused by the bandwidth of the applied field.

  7. Ab initio theoretical comparative study of magnetic coupling in KNiF3 and K2NiF4s

    Moreira, Ibério de Pinho Ribeiro; Illas i Riera, Francesc

    1997-01-01

    The origin of magnetic coupling in KNiF3 and K2 NiF4 is studied by means of an ab initio cluster model approach. By a detailed study of the mapping between eigenstates of the exact nonrelativistic and spin model Hamiltonians it is possible to obtain the magnetic coupling constant J and to compare ab initio cluster-model values with those resulting from ab initio periodic Hartree-Fock calculations. This comparison shows that J is strongly determined by two-body interactions; this is a surprisi...

  8. Yttrium aluminium garnet under pressure: Structural, elastic, and vibrational properties from ab initio studies

    Monteseguro, V. [Departamento de Física and MALTA Consolider Team, Universidad de La Laguna, 38200 San Cristóbal de La Laguna, Santa Cruz de Tenerife (Spain); Rodríguez-Hernández, P.; Muñoz, A., E-mail: amunoz@ull.es [Departamento de Física and MALTA Consolider Team, Universidad de La Laguna, 38200 San Cristóbal de La Laguna, Santa Cruz de Tenerife (Spain); Instituto de Materiales y Nanotecnología. Universidad de La Laguna, 38200 San Cristóbal de La Laguna, Santa Cruz de Tenerife (Spain)

    2015-12-28

    The structural, elastic, and vibrational properties of yttrium aluminum garnet Y{sub 3}Al{sub 5}O{sub 12} are studied under high pressure by ab initio calculations in the framework of the density functional theory. The calculated ground state properties are in good agreement with the available experimental data. Pressure dependences of bond length and bulk moduli of the constituent polyhedra are reported. The evolution of the elastic constants and the major elastic properties, Young and shear modulus, Poisson's ratios, and Zener anisotropy ratio, are described. The mechanical stability is analyzed, on the light of “Born generalized stability criteria,” showing that the garnet is mechanically unstable above 116 GPa. Symmetries, frequencies, and pressure coefficients of the Raman-active modes are discussed on the basis of the calculated total and partial phonon density of states, which reflect the dynamical contribution of each atom. The relations between the phonon modes of Y{sub 3}Al{sub 5}O{sub 12} and the internal and external molecular modes of the different polyhedra are discussed. Infrared-active modes, as well as the silent modes, and their pressure dependence are also investigated. No dynamical instabilities were found below 116 GPa.

  9. Ab initio study of the magnetostructural properties of MnAs

    Rungger, Ivan; Sanvito, Stefano

    2006-07-01

    The magnetic and structural properties of MnAs are studied with ab initio methods and by mapping total energies onto a Heisenberg model. The stability of the different phases is found to depend mainly on the volume and on the amount of magnetic order, confirming previous experimental findings and phenomenological models. It is generally found that for large lattice constants the ferromagnetic state is favored, whereas for small lattice constants different antiferromagnetic states can be stabilized. In the ferromagnetic state the structure with minimal energy is always hexagonal, whereas it becomes orthorhombically distorted if there is an antiferromagnetic alignment of the magnetic moments in the hexagonal plane. For the paramagnetic state the stable cell is found to be orthorhombic up to a critical lattice constant of about 3.7Å , above which it remains hexagonal. This leads to the second-order structural phase transition between paramagnetic states at about 400K , where the lattice parameter increases above this critical value with rising temperature due to the thermal expansion. We also evaluate the magnetic susceptibility as a function of temperature, from which a semiquantitative description of the MnAs phase diagram emerges.

  10. Ab initio study of a Y-doped Σ31 grain boundary in alumina

    2008-01-01

    The atomic structures and energetics of clean and Y-doped general grain boundary (GB) Σ31/(0001) models in α-Al2O3 are studied by a series of high precision ab initio calculations. A large supercell with 700 atoms and periodic boundary conditions is adopted for undoped and Y-doped GB with different substitution sites and con-centrations. It is shown that Y atoms preferably segregate to the central column of the 7-member Al ring. This is explained as more favorable bond formation for Y in this position and lower GB energy. The calculated GB formation energy for the clean and Y-doped cases is respectively 3.99 and 3.67 J/m2. On the average, the GB region in Σ31 has a slightly lower charge density than the bulk crystalline region. In addtition, the GB induces a long ranged asymmetric electrostatic potential distri-bution on each side of the grain boundary.

  11. Ab initio kinetics studies of hydrogen atom abstraction from methyl propanoate.

    Tan, Ting; Yang, Xueliang; Ju, Yiguang; Carter, Emily A

    2016-02-14

    The kinetics of hydrogen abstraction by five radicals (H, CH3, O((3)P), OH, and HO2) from a biodiesel surrogate, methyl propanoate (MP), is theoretically investigated. We employ high-level ab initio quantum chemistry methods, coupled-cluster singles and doubles with perturbative triples correction (CCSD(T)) and multi-reference singles and doubles configuration interaction (MRSDCI) with the Davidson-Silver (DS) correction, and obtain chemically accurate reaction energetics. Overall, MRSDCI + DS predicts comparable energetics to CCSD(T) for MP + H/CH3/O/OH. The rate constants are computed using transition state theory (TST-Rice-Ramsperger-Kassel-Marcus theory) in conjunction with the separable-hindered-rotor approximation, variable reaction coordinate TST, and the multi-structure all-structure (MS-AS) approach. A simplified method, semi-multi-structure, is also employed for MP + OH/HO2, and the rate coefficients with this less expensive method are in good agreement with the results obtained with the MS-AS method. The fitted modified Arrhenius expressions are provided over a temperature range of 250 to 2000 K. The predicted rate coefficients for MP + OH agree remarkably well with experimental data over a wide temperature range. Branching ratio analysis of all the studied reactions shows that abstractions of the secondary H atoms within MP are expected to dominate the consumption of fuel at low temperatures, and the contributions of abstractions from the two methyl groups increase with temperature for all abstracting radicals. PMID:26796249

  12. Why does bromine square palladium off? An ab initio study of brominated palladium and its nanomorphology.

    Yoo, Su-Hyun; Lee, Ji-Hwan; Delley, Bernard; Soon, Aloysius

    2014-09-14

    A first-principles description and prediction of brominated nanocrystals of Pd is presented. In particular, we conducted an extensive study of the adsorption behaviour of Br on various Pd surfaces (including both low and high Miller-index surfaces) as a function of its surface coverage. By coupling our calculated surface energies with ab initio (electrochemical) thermodynamics and the Gibbs-Wulff shape model, we find that the relative stability of the Pd surfaces is strongly modified by Br, allowing high Miller-index surfaces of Pd (namely the (210) surface) to become competitively favourable at moderate concentrations of Br. We also show that Pd nanoparticles assume a cube-like crystal shape at high concentrations of Br, exposing mainly the (100) facets with a Br surface coverage of 0.5 ML. This not only confirms and explains recent solution synthesis results, but also provides a quantitative atomic picture of the exposed surface facets, which is crucial in understanding the local surface chemistry of shape-controlled nanoparticles for better nanocatalyst design. PMID:25075669

  13. Ab initio computational study of vincristine as a biological active compound: NMR and NBO analyses

    Shiva Joohari

    2015-06-01

    Full Text Available Vincristine is a biological active alkaloid that has been used clinically against a variety of neoplasms. In the current study we have theoretically investigated the magnetic properties of titled compound to predict physical and chemical properties of vincristine as a biological inhibitor. Ab initio computation using HF and B3LYP with 3-21G(d and 6-31G(d level of theory have been performed and then magnetic shielding tensor (, ppm, shielding asymmetry (, magnetic shielding anisotropy (aniso, ppm, the skew of a tensor (K, chemical shift anisotropy ( and chemical shift ( were calculated to indicate the details of the interaction mechanism between microtubules and vincristine. Moreover, EHOMO, ELUMO and Ebg were evaluated. The maximum and minimum values of Ebg were found in HF/3-21g and B3LYP/3-21g respectively. It was also uggested that O24, O37, O49 and O55 with minimum values of iso, are active sites of titled compound. Furthermore the calculated chemical shifts were compared with experimental data in DMSO and CDCl3 solvents.

  14. ab initio Studies on Molecular Conductor (BEDSe-TTF)2[Fe(CN)5NO

    YAO Kai-Lun; TU Hai-Bo; WANG Wei-Zhong

    2001-01-01

    In this paper the ab initio study using pseudopotential plane wave method with the local spin density functional approximation is presented for the molecular conductor (BEDSe-TTF)2[Fe(CN)5NO]. The mean electronic density distributions are obtained, and we find that the extended π orbital of the selenium does not affect the properties of material as assumed in other papers and the "side-by-side" type S...S interaction is the primary interaction between donors. From band structure calculations we analyze the influence of the NO groups on the electronic structure and magnetic properties of molecule. It is shown that the itinerant electrons important to electronic properties in these types of hybrids are delocalized electrons contributed by NO groups, instead of by the 3d electrons of Fe. Additionally, we have found that the localized magnetic moment is also contributed by the NO groups in this molecular conductor. From total energy calculations the molecular structure with the lowest energy is found due to the interaction between split spins, and the particular positions of the NO groups are obtained.

  15. Influence of the "second gap" on the transparency of transparent conducting oxides: An ab initio study

    Ha, Viet-Anh; Waroquiers, David; Rignanese, Gian-Marco; Hautier, Geoffroy

    2016-05-01

    Transparent conducting oxides (TCOs) are essential to many technologies. These materials are doped (n- or p-type) oxides with a large enough band gap (ideally >3 eV) to ensure transparency. However, the high carrier concentration present in TCOs leads additionally to the possibility for optical transitions from the occupied conduction bands to higher states for n-type materials and from lower states to the unoccupied valence bands for p-type TCOs. The "second gap" formed by these transitions might limit transparency, and a large second gap has been sometimes proposed as a design criteria for high performance TCOs. Here, we study the influence of this second gap on optical absorption using ab initio computations for several well-known n- and p-type TCOs. Our work demonstrates that most known n-type TCOs do not suffer from second gap absorption in the visible even at very high carrier concentrations. On the contrary, p-type oxides show lowering of their optical transmission for high carrier concentrations due to second gap effects. We link this dissimilarity to the different chemistries involved in n- versus typical p-type TCOs. Quantitatively, we show that second gap effects lead to only moderate loss of transmission (even in p-type TCOs) and suggest that a wide second gap, while beneficial, should not be considered as a needed criteria for a working TCO.

  16. A theoretical study of benzaldehyde derivatives as tyrosinase inhibitors using Ab initio calculated NQCC parameters

    Marjan Rafiee

    2015-09-01

    Full Text Available Tyrosinase is a multifunctional copper-containing enzyme. It can catalyze two distinct reactions of melanin synthesis and benzaldehyde derivatives, which are potential tyrosinase inhibitors. To find the relationships between charge distributions of benzaldehyde and their pharmaceutical behavior, the present study aimed at investigating nuclear quadrupole coupling constants of quadrupolare nuclei in the functional benzaldehyde group and calculating some its derivatives. In addition, the differences between the electronic structures of various derivatives of this depigmenting drug were examined. All ab initio calculations were carried out using Gaussian 03. The results predicted benzaldehyde derivatives to be bicentral inhibitors; nevertheless, the oxygen or hydrogen contents of the aldehyde group were not found to be the only active sites. Furthermore with the presence of the aldehyde group, the terminal methoxy group in C4 was found to contribute to tyrosinase inhibitory activities. In addition, an oxygen atom with high charge density in the side chain was found to play an important role in its inhibitory effect.

  17. Ab initio studies on the lattice thermal conductivity of silicon clathrate frameworks II and VIII

    Härkönen, Ville J.; Karttunen, Antti J.

    2016-01-01

    The lattice thermal conductivities of silicon clathrate frameworks II and VIII are investigated by using ab initio lattice dynamics and an iterative solution of the linearized Boltzmann transport equation (BTE) for phonons. Within the temperature range 100-350 K, the clathrate structures II and VIII were found to have lower lattice thermal conductivity values than the silicon diamond structure (d -Si) by factors of 1/2 and 1/3, respectively. The main reason for the lower lattice thermal conductivity of the clathrate structure II in comparison to d -Si was found to be the harmonic phonon spectra, while in the case of the clathrate structure VIII, the difference is mainly due to the harmonic phonon spectra and partly due to the shorter relaxation times of phonons. In the studied clathrate frameworks, the anharmonic effects have larger impact on the lattice thermal conductivity than the size of the unit cell. For the structure II, the predicted lattice thermal conductivity differs approximately by a factor of 20 from the previous experimental results obtained for a polycrystalline sample at room temperature.

  18. Ab initio computational studies on molecular conformation of N-methyl-glyphosate

    Kaliannan, P.; Naseer Ali, M. Mohamed; Venuvanalingam, P.

    Conformational analysis of N-methyl-glyphosate has been carried out using an ab initio molecular orbital (MO) method at the HF/3-21G* levels of theory and the results are compared with the results of a previously studied compound, namely glyphosate. The potential energy surface of the molecule obtained by varying the central torsion angles (Φ, Ψ) was investigated in detail. Fourteen conformers with 5 kcal mol-1 energy cut-off have been selected from the potential energy surface for geometry optimization to locate the true minimum on the conformational space. The minimum has been found to be at (-62°, 110°) for the central torsion angles. This conformation is stabilized by hydrogen bond interactions (O-H···O and C-H···O) and the interactions due to protons nearer to each other. This cationic field leads to the formation of a hydrophobic patch in this structure, as well as in the structures closer to the global minimum. This patch may destabilize the favourable interaction of N-methyl-glyphosate with the surrounding amino acid residues in the binding cavity as they form the cationic field throughout the glyphosate binding region.

  19. Energetics, diffusion, and magnetic properties of cobalt atom in a monolayer graphene: An ab initio study

    We use ab initio methods to study the binding, diffusion, and magnetic properties of cobalt atom embedded in graphene vacancies. We investigate the diffusion of Co-monovacancy (Co-MV) and Co-divacancy (Co-DV) defect complexes, and determine the minimum energy path (MEP), as well as the activation energy barrier of migration. We obtained similar activation energy barriers, of ∼5.8 eV, for Co-MV and Co-DV diffusion, respectively. Our calculations also suggest that, at electron–irradiation energy of 200 keV as used in a related experiment, the maximum energy transfer to the Co atom, of approximately 9.0 eV is sufficiently high to break metal-carbon bonding. The incident electron energy is also high enough to displace graphene's carbon atoms from their lattice positions. The breaking of metal-carbon bonding and the displacement of graphene atoms may act to facilitate the migration of Co. We conclude therefore that the detrapping and diffusion of cobalt as observed experimentally is likely to be radiation-induced, similar to what has been observed for Au and Fe in electron-irradiated graphene. Furthermore, we show that Co migration in graphene is such that its magnetic moment varies along the diffusion path. The magnetic moment of Co is consistently higher in Co-DV diffusion when compared to that of Co-MV diffusion

  20. Ab initio pseudopotential studies of cubic BC2N under high pressure

    We present the results of a systematic study of the structural, electronic, and vibrational properties of various cubic BC2N phases under high pressure. Ab initio pseudopotential total-energy and phonon calculations have been carried out to examine the changes in the structural parameters, bonding behaviours, band structures, and dynamic instabilities caused by phonon softening in these phases. We find that an experimentally synthesized high-density phase of cubic BC2N exhibits outstanding stability in the structural and electronic properties up to very high pressures. On the other hand, another experimentally identified phase with lower density and lower symmetry undergoes a dramatic structural transformation with a volume and bond-length collapse and a concomitant semi-metal to semiconductor transition. A third phase is predicted to be favourable over the above-mentioned lower-density phase by the enthalpy calculations. However, the dynamic phonon calculations reveal that it develops imaginary phonon modes and, therefore, is unstable in the experimental pressure range. The calculations indicate that its synthesis may be achieved at reduced pressures. These results provide a comprehensive understanding for the high-pressure behaviour of the cubic BC2N phases and reveal their interesting properties that can be verified by experiments

  1. AB INITIO HF AND DFT STUDIES ON MOLECULAR STRUCTURE AND VIBRATIONAL ANALYSIS OF 2,5-DIBROMOPYRIDINE

    ÇIRAK, Çağrı; KOÇ, Nurettin

    2014-01-01

    Theoretical study on molecular structure and vibrational spectra of 2,5-dibromopyridine (2,5-DBP) have been investigated. The optimized geometry, theoretical vibration frequencies and intensities were calculated by using ab initio Hartree-Fock and density functional B3LYP method with 6-31G(d,p) basis sets. The vibrational analysis of title molecule was done and its optimized geometry parameters (bond lengths and bond angles) were given. Scaled theoretical frequencies have been compared with e...

  2. Structural, dynamical, electronic, and bonding properties of laser-heated silicon: An ab initio molecular-dynamics study

    Silvestrelli, P.-L.; Alavi, A; Parrinello, M.; Frenkel, D

    1997-01-01

    The method of ab initio molecular dynamics, based on finite-temperature density-functional theory, is used to simulate laser heating of crystalline silicon. We found that a high concentration of excited electrons dramatically weakens the covalent bonding. As a result the system undergoes a melting transition to a metallic state. We studied several structural, dynamical, electronic, and bonding properties of this phase of silicon. In contrast to ordinary liquid silicon, this liquid is characte...

  3. On the hierarchical parallelization of ab initio simulations

    Ruiz-Barragan, Sergi; Shiga, Motoyuki

    2016-01-01

    A hierarchical parallelization has been implemented in a new unified code PIMD-SMASH for ab initio simulation where the replicas and the Born-Oppenheimer forces are parallelized. It is demonstrated that ab initio path integral molecular dynamics simulations can be carried out very efficiently for systems up to a few tens of water molecules. The code was then used to study a Diels-Alder reaction of cyclopentadiene and butenone by ab initio string method. A reduction in the reaction energy barrier is found in the presence of hydrogen-bonded water, in accordance with experiment.

  4. Experimental and ab initio study of the photofragmentation of DNA and RNA sugars

    The photoelectron-photoion-photoion coincidence method is used to measure the photodissociation of doubly charged D-ribose (C5H10O5), the RNA sugar molecules, and 2-deoxy-D-ribose (C5H10O4), the DNA sugar molecules, following normal Auger decay after initial C 1s and O 1s core ionizations. The fragment identification is facilitated by measuring isotopically labeled D-ribose, such as D-ribose deuterated at C(1), and with 13C at the C(5) position. Ab initio quantum chemistry calculations are used to gain further insight into the abundant appearance of the CHO+ fragment.

  5. Ab initio study of beryllium-decorated fullerenes for hydrogen storage

    Lee, Hoonkyung; Huang, Bing; Duan, Wenhui; Ihm, Jisoon

    2010-01-01

    We have found that a beryllium (Be) atom on nanostructured materials with H2 molecules generates a Kubas-like dihydrogen complex [H. Lee et al. arXiv:1002.2247v1 (2010)]. Here, we investigate the feasibility of Be-decorated fullerenes for hydrogen storage using ab initio calculations. We find that the aggregation of Be atoms on pristine fullerenes is energetically preferred, resulting in the dissociation of the dihydrogen. In contrast, for boron (B)-doped fullerenes, Be atoms prefer to be ind...

  6. The electronic structure of wurtzite and zincblende AlN: an ab initio comparative study

    This work deals with the electronic properties, in different crystal phases, of AlN (wurtzite and zincblende) compounds computed using an all electron ab initio linearized augmented plane wave method. Results include band structure, total and partial density of states, charge density and the ionicity factor. Most of the calculated band parameters, of direct bandgap, total- and upper-valence bandwidths and antisymmetric gap for wurtzite-AlN are close to those of c-AlN to within 1%. The charge distributions have similar features, meaning that AlN has the same ionicity factor in both structures. (author)

  7. Low-lying electronic states of nitrosyl cyanide (NCNO): An ab initio MCHF study

    Ab initio multiconfiguration Hartree--Fock calculations of the electronic structure of several low-lying electronic states of nitrosyl cyanide (NCNO) are reported. The essential features of the electronic structure of these states were analyzed to provide a qualitative correlation diagram for the dissociation process NCNO → NC+NO. It is found that the four lowest-lying states 1/sup ,/3A' and 1/sup ,/3A'' are connected directly to ground state products. Excitation energies and geometrical parameters for these states and for the 2- 1A'' and 2- 1A' states are presented. Implications of these findings for NCNO photodissociation processes are discussed

  8. Ab initio studies of ionization potentials of hydrated hydroxide and hydronium

    Swartz, Charles W

    2013-01-01

    The ionization potential distributions of hydrated hydroxide and hydronium are computed with many-body approach for electron excitations with configurations generated by {\\it ab initio} molecular dynamics. The experimental features are well reproduced and found to be closely related to the molecular excitations. In the stable configurations, the ionization potential is mainly perturbed by water molecules within the first solvation shell. On the other hand, electron excitation is delocalized on both proton receiving and donating complex during proton transfer, which shifts the excitation energies and broadens the spectra for both hydrated ions.

  9. Ab-initio study of magnetic properties and phase transitions in Ga (Mn) N with Monte Carlo approach

    Sbai, Y.; Ait Raiss, A.; Salmani, E.; Bahmad, L.; Benyoussef, A.

    2015-12-01

    On the basis of ab-initio calculations and Monte Carlo simulations the magnetic and electronic properties of Gallium nitride (GaN) doped with the transition metal Manganese (Mn) were studied. The ab initio calculations were done using the AKAI-KKR-CPA method within the Local Density Approximation (LDA) approximation. We doped our Diluted Magnetic Semiconductor (DMS), with different concentrations of magnetic impurities Mn and plotted the density of state (DOS) for each one. Showing a half-metallic behavior and ferromagnetic state especially for Ga0.95Mn0.05N making this DMS a strong candidate for spintronic applications. Moreover, the magnetization and susceptibility of our system as a function of the temperature has been calculated and give for various system size L to study the size effect. In addition, the transition temperature was deduced from the peak of the susceptibility. The Ab initio results are in good agreement with literature especially for (x=0.05) of Mn which gives the most interesting results.

  10. Revealing halogen bonding interactions with anomeric systems: an ab initio quantum chemical studies.

    Lo, Rabindranath; Ganguly, Bishwajit

    2015-02-01

    A computational study has been performed using MP2 and CCSD(T) methods on a series of O⋯X (X=Br, Cl and I) halogen bonds to evaluate the strength and characteristic of such highly directional noncovalent interactions. The study has been carried out on a series of dimeric complexes formed between interhalogen compounds (such as BrF, BrCl and BrI) and oxygen containing electron donor molecule. The existence and consequences of the anomeric effect of the electron donor molecule has also been investigated through an exploration of halogen bonding interactions in this halogen bonded complexes. The ab initio quantum chemical calculations have been employed to study both the nature and directionality of the halogen molecules toward the sp(3) oxygen atom in anomeric systems. The presence of anomeric nO→σ*CN interaction involves a dominant role for the availability of the axial and equatorial lone pairs of donor O atom to participate with interhalogen compounds in the halogen-bonded complexes. The energy difference between the axial and equatorial conformers with interhalogen compounds reaches up to 4.60 kJ/mol, which however depends upon the interacting halogen atoms and its attaching atoms. The energy decomposition analysis further suggests that the total halogen bond interaction energies are mainly contributed by the attractive electrostatic and dispersion components. The role of substituents attached with the halogen atoms has also been evaluated in this study. With the increase of halogen atom size and the positive nature of σ-hole, the halogen atom interacted more with the electron donor atom and the electrostatic contribution to the total interaction energy enhances appreciably. Further, noncovalent interaction (NCI) studies have been carried out to locate the noncovalent halogen bonding interactions in real space. PMID:25522359

  11. Origin of the Hadži ABC structure: An ab initio study

    Medium and strong hydrogen bonds are well known to give rise to broad features in the vibrational spectrum often spanning several hundred wavenumbers. In some cases, these features can span over 1000 cm−1 and even contain multiple broad peaks. One class of strongly hydrogen-bonded dimers that includes many different phosphinic, phosphoric, sulfinic, and selenic acid homodimers exhibits a three-peaked structure over 1500 cm−1 broad. This unusual feature is often referred to as the Hadži ABC structure. The origin of this feature has been debated since its discovery in the 1950s. Only a couple of theoretical studies have attempted to interpret the origin of this feature; however, no previous study has been able to reproduce this feature from first principles. Here, we present the first ab initio calculation of the Hadži ABC structure. Using a reduced dimensionality calculation that includes four vibrational modes, we are able to reproduce the three-peak structure and much of the broadness of the feature. Our results indicate that Fermi resonances of the in-plane bend, out-of-plane bend, and combination of these bends play significant roles in explaining this feature. Much of the broadness of the feature and the ability of the OH stretch mode to couple with many overtone bending modes are captured by including an adiabatically separated dimer stretch mode in the model. This mode modulates the distance between the monomer units and accordingly the strength of the hydrogen-bonds causing the OH stretch frequency to shift from 2000 to 3000 cm−1. Using this model, we were also able to reproduce the vibrational spectrum of the deuterated isotopologue which consists of a single 500 cm−1 broad feature. Whereas previous empirical studies have asserted that Fermi resonances contribute very little to this feature, our study indicates that while not appearing as a separate peak, a Fermi resonance of the in-plane bend contributes substantially to the feature

  12. Vibrational spectroscopic investigations, ab initio and DFT studies on 7-bromo-5-chloro-8-hydroxyquinoline

    Arjunan, V.; Mohan, S.; Ravindran, P.; Mythili, C. V.

    2009-05-01

    The Fourier transform infrared (FTIR) and FT-Raman spectra of 7-bromo-5-chloro-8-hydroxyquinoline (BCHQ) have been measured in the range 4000-400 and 4000-100 cm -1, respectively. Complete vibrational assignment and analysis of the fundamental modes of the compound were carried out using the observed FTIR and FT-Raman data. The geometry was optimised without any symmetry constrains using the DFT/B3LYP and HF methods with 6-31G** basis set. The vibrational frequencies which were determined experimentally are compared with those obtained theoretically from ab initio HF and density functional theory (DFT) gradient calculations employing the HF/6-31G** and B3LYP/6-31G** methods for the optimised geometry of the compound. The structural parameters and normal modes of vibration obtained from HF and DFT methods are in good agreement with the experimental data. Normal coordinate analysis was also carried out with ab initio force fields utilising Wilson's FG matrix method.

  13. Ab initio study on electron excitation and electron transfer in tryptophan-tyrosine system

    In this article, ab initio calculation has been performed to evaluate the transition energy of electronic excitation in tryptophan and tyrosine by using semiempirical molecular orbital method AM1 and complete active space self-consistent field method. The solvent effect has been considered by means of the conductor-like screening model. After geometric optimizations of isolated tryptophan and tyrosine, and their corresponding radicals and cations, reaction heat of these electron transfer reactions have been obtained by the means of complete active space self-consistent field method. The transition energies from the ground state, respectively, to the lowest excited state and to the lowest triplet state of these two amino acids are also calculated and compared with the experimentally observed values. The ionization potential and electron affinity are also calculated for tryptophan and tyrosine employing Koopmans' theorem and ab initio calculation. Compared with the experimental measurements, the theoretical results are found satisfactory. Theoretical results give good explanations on the experimental phenomena that N3· can preferably oxide the side chain of tryptophan residue and then the electron transfer from tyrosine residue to tryptophan residue follows in peptides involving tryptophan and tyrosine

  14. An ab initio cluster study of the structure of the Si(001) surface

    Ab initio calculations, employing double zeta plus polarization (DZP) basis sets and generalized valence bond (GVB) wave functions, have been performed on clusters of varying size, to investigate the utility of such clusters as prototypes for the study of silicon surfaces, and to investigate the effect of the level of theory used on predicted results. This work builds on landmark papers by Goddard in 1982 and Paulus in 1998 that demonstrate that a single reference wave function description of the silicon dimer bond is incorrect, and that a multireference description results in a symmetric dimer in a silicon cluster containing one dimer. In this work, it is shown that the imposition of arbitrary geometrical constraints (fixing subsurface atoms at lattice positions) on cluster models of the Si(100) surface can also lead to nonphysical results. Calculations on the largest clusters, without geometrical constraints, reveal that surface rearrangement due to dimer bond formation is ''felt'' several layers into the bulk. The predicted subsurface displacements compare favorably to experiment. Thus, small clusters, such as Si9H12, cannot adequately represent bulk behavior. Vibrational analysis shows that dimer buckling modes require minimal excitation energy, so the experimental observation of buckled dimers on silicon surfaces may reflect the ease with which a symmetric dimer can be perturbed from its minimum energy structure. In the study of surface reconstruction and relaxation, and the associated issue of the buckling of dimer surfaces, it is critical to use adequate wave functions. As shown in this work and previously by Goddard and Paulus, this generally means that multireference treatments are needed to correctly treat the dangling bonds. (c) 2000 American Institute of Physics

  15. Reaction mechanisms and kinetics of the iminovinylidene radical with NO: Ab initio study

    Hsiao, Ming-Kai; Chung, Yi-Hua; Hung, Yu-Ming; Chen, Hui-Lung, E-mail: chl3@faculty.pccu.edu.tw [Department of Chemistry and Institute of Applied Chemistry, Chinese Culture University, Taipei 111, Taiwan (China)

    2014-05-28

    The nitric oxide (NO) is a notorious compound for polluting environment. Recent year, removing nitric oxide from the atmosphere becomes a focus of the investigation. In our work, we study the iminovinylidene (HNCC) radical reacted with NO molecule. The mechanism and kinetic for reaction of the HNCC radical with the NO molecule is investigated via considering the possible channels of the N and O atoms of NO attacking the N and C atoms of the HNCC based on the high level ab initio molecular orbital calculations in conjunction with variational TST and RRKM calculations. The species involved have been optimized at the B3LYP/6-311++G(3df,2p) level and their single-point energies are refined by the CCSD(T)/aug-cc-PVQZ//B3LYP/6-311++G(3df,2p) method. The calculated potential energy surfaces indicated that energetically the most favorable channel for the HNCC + NO reaction was predicted to be the formation of HNC+CNO (P8) product via the addition reaction of the C atom of HNCC radical and the N atom of NO with the head to head orientation. To rationalize the scenario of the calculated results, we also employ the Fukui functions and HSAB theory to seek for a possible explanation. In addition, the reaction rate constants were calculated using VariFlex code, and the results show that the total rate coefficient, k{sub total}, at Ar pressure 760 Torr can be represented with an equation: k{sub total} = 6.433 × 10{sup −11} T {sup 0.100} exp(0.275 kcal mol{sup −1}/RT) at T = 298–3000 K, in units of cm{sup 3} molecule{sup −1} s{sup −1}.

  16. Diffusion within α-CuI studied using ab initio molecular dynamics simulations

    The structure and dynamics of superionic α-CuI are studied in detail by means of ab initio Born-Oppenheimer molecular dynamics simulations. The extreme cation disorder and a soft immobile face centred cubic sublattice are evident from the highly diffuse atomic density profiles. The Cu-Cu pair distribution function and distribution of Cu-I-Cu bond angles possess distinct peaks at 2.6 A and 60 deg. respectively, which are markedly lower than the values expected from the average cationic density, pointing to the presence of pronounced short-range copper-copper correlations. Comparison with lattice static calculations shows that these correlations and the marked shift in the cationic density profile in the (111) directions are associated with a locally distorted cation sublattice, and that the movements within the tetrahedral cavities involve rapid jumps into and out of shallow basins on the system potential energy surface. On average, the iodines are surrounded by three coppers within their first coordination shell, with the fourth copper being located in a transition zone between two neighbouring iodine cavities. However, time-resolved analysis reveals that the local structure actually involves a mixture of threefold-, fourfold- and fivefold-coordinated iodines. Examination of the ionic trajectories shows that the copper ions jump rapidly to nearest neighbouring tetrahedral cavities (aligned in the (100) directions) following a markedly curved trajectory and often involving short-lived (∼1 ps) interstitial positions. The nature of the correlated diffusion underlying the unusually high fraction of coppers with short residence time can be attributed to the presence of a large number of 'unsuccessful' jumps and the likelihood of cooperative motion of pairs of coppers. The calculated diffusion coefficient at 750 K, DCu = 2.8 x 10-5 cm2 s-1, is in excellent agreement with that found experimentally.

  17. Diffusion within α-CuI studied using ab initio molecular dynamics simulations

    Mohn, Chris E.; Stølen, Svein; Hull, Stephen

    2009-08-01

    The structure and dynamics of superionic α-CuI are studied in detail by means of ab initio Born-Oppenheimer molecular dynamics simulations. The extreme cation disorder and a soft immobile face centred cubic sublattice are evident from the highly diffuse atomic density profiles. The Cu-Cu pair distribution function and distribution of Cu-I-Cu bond angles possess distinct peaks at 2.6 Å and 60° respectively, which are markedly lower than the values expected from the average cationic density, pointing to the presence of pronounced short-range copper-copper correlations. Comparison with lattice static calculations shows that these correlations and the marked shift in the cationic density profile in the lang111rang directions are associated with a locally distorted cation sublattice, and that the movements within the tetrahedral cavities involve rapid jumps into and out of shallow basins on the system potential energy surface. On average, the iodines are surrounded by three coppers within their first coordination shell, with the fourth copper being located in a transition zone between two neighbouring iodine cavities. However, time-resolved analysis reveals that the local structure actually involves a mixture of threefold-, fourfold- and fivefold-coordinated iodines. Examination of the ionic trajectories shows that the copper ions jump rapidly to nearest neighbouring tetrahedral cavities (aligned in the lang100rang directions) following a markedly curved trajectory and often involving short-lived (~1 ps) interstitial positions. The nature of the correlated diffusion underlying the unusually high fraction of coppers with short residence time can be attributed to the presence of a large number of 'unsuccessful' jumps and the likelihood of cooperative motion of pairs of coppers. The calculated diffusion coefficient at 750 K, DCu = 2.8 × 10-5 cm2 s-1, is in excellent agreement with that found experimentally.

  18. Ab initio study of different structures of CaC: Magnetism, bonding, and lattice dynamics

    Highlights: •Electronic structure and phonon calculations are performed on seven CaC structures. •The more covalent structures have lower energies and nonmagnetic ground state. •Ferromagnetism in the ionic phases comes from the sharp C p band at the Fermi level. •Tendency of C atoms for dimerization may lead to structural instabilities in CaC. •Nonmagnetic B33 CaC is stable in a wide range of temperatures and pressures. -- Abstract: On the basis of ab initio pseudopotential calculations, we study structural, magnetic, dynamical, and mechanical properties of the hypothetical CaC ionic compound in the rock-salt (RS), B2, zinc-blende (ZB), wurtzite (WZ), NiAs (NA), anti-NiAs (NA∗), and CrB (B33) structures. It is argued that the ZB, WZ, NA, and RS structures are more ionic while the NA∗, B2, and B33 structures are more covalent systems. As a result of that, the nonmagnetic B33–CaC is the energetically preferred system, while the more ionic structures prefer a ferromagnetic ground state with high Fermi level spin polarization. The observed ferromagnetism in the more ionic systems is attributed to the sharp partially filled p states of carbon atom in the system. In the framework of density functional perturbation theory, the phonon spectra of these systems are computed and the observed dynamical instabilities of the NA∗ and B2 structures are explained in terms of the covalent bonds between carbon atoms. The calculated Helmholtz and Enthalpy free energies indicate the highest stability of the B33 structure in a wide range of temperatures and pressures. Among the ferromagnetic structures, RS–CaC and ZB–CaC are reported, respectively, to be the most and the least metastable systems in various thermodynamics conditions. Several mechanical properties of the dynamically stable structures of CaC are determined from their phonon spectra

  19. Ab initio molecular dynamics study of the interlayer and micropore structure of aqueous montmorillonite clays

    Suter, James L.; Kabalan, Lara; Khader, Mahmoud; Coveney, Peter V.

    2015-11-01

    Ab initio molecular dynamics simulations have been performed to gain an understanding of the interfacial microscopic structure and reactivity of fully hydrated clay edges. The models studied include both micropore and interlayer water. We identify acidic sites through dissociation mechanisms; the resulting ions can be stabilized by both micropore and interlayer water. We find clay edges possess a complex amphoteric behavior, which depends on the face under consideration and the location of isomorphic substitution. For the neutral (1 1 0) surface, we do not observe any dissociation on the timescale accessible. The edge terminating hydroxyl groups participate in a hydrogen bonded network of water molecules that spans the interlayer between periodic images of the clay framework. With isomorphic substitutions in the tetrahedral layer of the (1 1 0) clay edge, we find the adjacent exposed apical oxygen behaves as a Brönsted base and abstracts a proton from a nearby water molecule, which in turn removes a proton from an AlOH2 group. With isomorphic substitutions in the octahedral layer of the (1 1 0) clay edge the adjacent exposed apical oxygen atom does not abstract a proton from the water molecules, but increases the number of hydrogen bonded water molecules (from one to two). Acid treated clays are likely to have both sites protonated. The (0 1 0) surface does not have the same interfacial hydrogen bonding structure; it is much less stable and we observe dissociation of half the terminal SiOH groups (tbnd Sisbnd Osbnd H → tbnd Sisbnd O- + H+) in our models. The resulting anions are stabilized by solvation from both micropore and interlayer water molecules. This suggests that, when fully hydrated, the (0 1 0) surface can act as a Brönsted acid, even at neutral pH.

  20. Ab initio studies on the adsorption and implantation of Al and Fe to nitride materials

    The formation of transfer material products on coated cutting and forming tools is a major failure mechanism leading to various sorts of wear. To describe the atomistic processes behind the formation of transfer materials, we use ab initio to study the adsorption energy as well as the implantation barrier of Al and Fe atoms for (001)-oriented surfaces of TiN, Ti0.50Al0.50N, Ti0.90Si0.10N, CrN, and Cr0.90Si0.10N. The interactions between additional atoms and nitride-surfaces are described for pure adhesion, considering no additional stresses, and for the implantation barrier. The latter, we simplified to the stress required to implant Al and Fe into sub-surface regions of the nitride material. The adsorption energies exhibit pronounced extrema at high-symmetry positions and are generally highest at nitrogen sites. Here, the binary nitrides are comparable to their ternary counterparts and the average adhesive energy is higher (more negative) on CrN than TiN based systems. Contrary, the implantation barrier for Al and Fe atoms is higher for the ternary systems Ti0.50Al0.50N, Ti0.90Si0.10N, and Cr0.90Si0.10N than for their binary counterparts TiN and CrN. Based on our results, we can conclude that TiN based systems outperform CrN based systems with respect to pure adhesion, while the Si-containing ternaries exhibit higher implantation barriers for Al and Fe atoms. The data obtained are important to understand the atomistic interaction of metal atoms with nitride-based materials, which is valid not just for machining operations but also for any combination such as interfaces between coatings and substrates or multilayer and phase arrangements themselves

  1. Study of the behaviour of cesium fission product in uranium dioxide by the ab initio method

    The knowledge of the behaviour of fission products in the nuclear fuel is very important for safety considerations and for understanding the evolution of the fuel properties under irradiation. In this work, we focussed mainly on the behaviour of caesium in UO2 through ab initio studies of its solubility at point defects in the matrix, its diffusion and its contribution to the formation of solid phases in the fuel. The role of electronic correlation effects of the f electrons of uranium on these properties and on the description of the defect free crystal, is assessed. The formation energies of the main point defects are calculated and their concentration as a function of fuel stoichiometry and temperature is estimated. The migration barriers and migration paths for the self-diffusion of oxygen and uranium vacancies and oxygen interstitials in UO2 are discussed. The solubility of Cs is found to be very low in UO2 in agreement with experimental findings. The most favourable trapping sites are determined as a function of oxygen concentration in the fuel. Our results show that in the hyper-stoichiometric regime, the diffusion of Cs from its most favourable trapping site is limited by the uranium vacancy diffusion mechanism. We also considered the formation of the main solid phases of caesium resulting from its oxidation (Cs2O, Cs2O2, CsO2) and from its interaction with the fuel (Cs2UO4), with molybdenum (Cs2MoO4) and with the zirconium of the clad (Cs2ZrO3), since the formation of such phases, their solubility and their interdependence will affect the release of caesium. (author)

  2. Crystal structure and magnetism in α -RuCl3 : An ab initio study

    Kim, Heung-Sik; Kee, Hae-Young

    2016-04-01

    α -RuCl3 has been proposed recently as an excellent playground for exploring Kitaev physics on a two-dimensional (2D) honeycomb lattice. However, structural clarification of the compound has not been completed, which is crucial in understanding the physics of this system. Here, using ab initio electronic structure calculations, we study a full three-dimensional (3D) structure of α -RuCl3 , including the effects of spin-orbit coupling (SOC) and electronic correlations. The three major results are as follows: (i) SOC suppresses dimerization of Ru atoms, which exists in other Ru compounds such as isostructural Li2RuO3 , and makes the honeycomb closer to an ideal one. (ii) The nearest-neighbor Kitaev exchange interaction between the jeff=1 /2 pseudospin strongly depends on the Ru-Ru distance and the Cl position, originating from the nature of the edge-sharing geometry. (iii) The optimized 3D structure without electronic correlations has P 3 ¯1 m space-group symmetry independent of SOC, but including electronic correlation changes the optimized 3D structure to either C 2 /m or C m c 21 within 0.1 meV per formula unit (f.u.) energy difference. The reported P 3112 structure is also close in energy. The interlayer spin-exchange coupling is a few percent of the in-plane spin-exchange terms, confirming that α -RuCl3 is close to a 2D system. We further suggest how to increase the Kitaev term via tensile strain, which sheds light in realizing the Kitaev spin-liquid phase in this system.

  3. On the room-temperature phase diagram of high pressure hydrogen: an ab initio molecular dynamics perspective and a diffusion Monte Carlo study.

    Chen, J; Ren, X.; Li, X Z; Alfè, D.; Wang, E

    2014-01-01

    The finite-temperature phase diagram of hydrogen in the region of phase IV and its neighborhood was studied using the ab initio molecular dynamics (MD) and the ab initio path-integral molecular dynamics (PIMD). The electronic structures were analyzed using the density-functional theory (DFT), the random-phase approximation, and the diffusion Monte Carlo (DMC) methods. Taking the state-of-the-art DMC results as benchmark, comparisons of the energy differences between structures generated from ...

  4. Ab initio study of the low-lying electronic states of the CaO molecule

    Graphical abstract: Highly correlated ab initio calculations have been performed for an accurate determination of the electronic structure of the low-lying electronic states of the CaO molecule. The computations are done using the aug-cc-pV5Z basis set for O and the cc-pCV5Z for Ca. The potential energy curves for the molecular states correlating to the lowest three asymptotes are calculated at the CASSCF level. The potential curves of the lowest five molecular states, X1Σ+, a3Π, A'1Π, b3Σ+ and A1Σ+, and the corresponding dipole moment functions have been determined using internally contracted multi-reference configuration interaction approaches. The spectroscopic constants associated with these electronic states are compared to experimental values. The corresponding electronic wavefunctions have also been analyzed using the dipole moment functions. Display Omitted Highlights: → The five lowest electronic states of Cao have been determined ab initio at a high level of accuracy. → Large active space, core-valence correlation and configuration interaction are required. → The multi-configurational nature of the electronic ground state is confirmed as well as its monovalent and divalent ionic nature using dipole moment analysis. → These interacting potentials will serve for future obtention of spin-rovibronic levels. - Abstract: Highly correlated ab initio calculations have been performed for an accurate determination of the electronic structure of the low-lying electronic states of the CaO molecule. The computations are done using the aug-cc-pV5Z basis set for O and the cc-pCV5Z for Ca. The potential energy curves for the molecular states correlating to the lowest three asymptotes are calculated at the CASSCF level. The potential curves of the lowest five molecular states, X1Σ+, a3Π, A'1Π, b3Σ+ and A1Σ+, and the corresponding dipole moment functions have been determined using internally contracted multi-reference configuration interaction approaches

  5. Theoretical study of silicon carbide under irradiation at the nano scale: classical and ab initio modelling; Etude theorique a l'echelle nanometrique du carbure de silicium sous irradiation: modelisation classique et ab initio

    Lucas, G

    2006-10-15

    The behaviour of silicon carbide under irradiation has been studied using classical and ab initio simulations, focusing on the nano scale elementary processes. First, we have been interested in the calculation of threshold displacement energies, which are difficult to determine both experimentally and theoretically, and also the associated Frenkel pairs. In the framework of this thesis, we have carried out simulations in classical and ab initio molecular dynamics. For the classical approach, two types of potentials have been used: the Tersoff potential, which led to non satisfactory results, and a new one which has been developed during this thesis. This potential allows a better modelling of SiC under irradiation than most of the empirical potentials available for SiC. It is based on the EDIP potential, initially developed to describe defects in silicon, that we have generalized to SiC. For the ab initio approach, the feasibility of the calculations has been validated and average energies of 19 eV for the C and 38 eV for the Si sublattices have been determined, close to the values empirically used in the fusion community. The results obtained with the new potential EDIP are globally in agreement with those values. Finally, the elementary processes involved in the crystal recovery have been studied by calculating the stability of the created Frenkel pairs and determining possible recombination mechanisms with the nudged elastic band method. (author)

  6. Ab initio study of pressure induced structural and electronic properties in TmPo

    Makode, Chandrabhan, E-mail: cbmakode@gmail.com; Pataiya, Jagdish; Sanyal, Sankar P. [Department of Physics, Barkatullah University, Bhopal-462026 (India); Panwar, Y. S.; Aynyas, Mahendra [Department of Physics, C.S.A. Govt. P.G. College, Sehore-466001 (India)

    2015-06-24

    We report an ab initio calculation of pressure induced structural phase transition and electronic properties of Thulium Polonide (TmPo).The total energy as a function of volume is obtained by means of self-consistent tight binding linear muffin-tin-orbital (TB-LMTO) method within the local density approximation (LDA). It is found that TmPo is stable in NaCl-type (B{sub 1}-phase) structure to CsCl-type (B{sub 2}-phase) structure of this compound in the pressure range of 7.0 GPa. We also calculate the lattice parameter (a{sub 0}), bulk modulus (B{sub 0}), band structure and density of states. From energy diagram it is observed that TmPo exhibit metallic behavior. The calculated values of equilibrium lattice parameter and bulk modulus are in general good agreement.

  7. Optoelectronic Properties, Elastic Moduli and Thermoelectricity of SrAlGa: An Ab Initio Study

    Roshan, Ali; Murtaza, G.; Takagiwa, Y.; Khenata, R.; Haleem, Uddin; Ullah, H.; A. Khan, S.

    2014-04-01

    Half-Heusler compounds are an impressive class of materials with a huge potential for different applications such as in future energy, especially in the fields of thermoelectrics and solar cells. We present ab initio total energy calculations within the modified Becke—Johnson generalized gradient approximation (mBJ-GGA) to obtain the physical properties of SrAlGa compounds. The structural, elastic, acoustic, electronic, chemical bonding, optical, and thermoelectric properties are calculated and compared with the available calculation data. The SrAlGa is found to be a small-band-gap (0.125-0.175 eV) material, suitable for thermoelectric applications with a relatively high Seebeck coefficient. Also, SrAlGa has the potential in the optoelectronic applications due to high optical conductivity and reflectivity in the infrared and visible region of electromagnetic spectra.

  8. Isomerism of OBe3F3+ cation: an ab initio study

    Ab initio MP2/6-31G*/HF/6-31G*+ZPE(HF/6-31G*) calculations of the potential energy surface in the vicinity of stationary points and the pathways of intramolecular rearrangements between low-lying structures of the OBe3F3+ cation detected in the mass spectra of μ4-Be4O(CF3COO)6 were carried out. Ten stable isomers with di- and tricoordinate oxygen atoms were localized. The relative energies of six structures lie in the range 0-8 kcal mol-1 and those of the rest four structures lie in the range 20-40 kcal mol-1. two most favorable isomers are a planar C2, isomer and a pyramidal C3 isomer

  9. Relaxation of the excited -(2-hydroxy benzylidene) aniline molecule: An ab initio and TD DFT study

    Biswajit Chowdhury; Rina De; Pinaky Sett; Joydeep Chowdhury

    2010-11-01

    The photophysical behaviour of N-(2-hydroxy benzylidene) aniline or most commonly known as salicylideneaniline (SA) has been investigated using the ab initio and DFT levels of theory. The quantum chemical calculations show that the optimized non planar enol (1) form of the SA molecule is the most stable conformer in the ground state and is marked by the twisting of the phenolic and anilino rings of the molecule. The geometry optimizations and the subsequent frequency calculations of the excited singlet electronic states of the various tautomeric forms of SA molecule were performed with the CIS level of theory. A detail theoretical investigation on the relaxation dynamics of the SA molecule has been presented. Possible explanation on the excitation wavelength dependence of the photochromic yield of the molecule is also reported.

  10. Ab initio and phenomenological studies of the static response of neutron matter

    Buraczynski, Mateusz

    2016-01-01

    We investigate the problem of periodically modulated strongly interacting neutron matter. We carry out ab initio non-perturbative auxiliary-field diffusion Monte Carlo calculations using an external sinusoidal potential in addition to phenomenological two- and three-nucleon interactions. Several choices for the wave function ansatz are explored and special care is taken to extrapolate finite-sized results to the thermodynamic limit. We perform calculations at various densities as well as at different strengths and periodicities of the one-body potential. Our microscopic results are then used to constrain the isovector term from energy-density functional theories of nuclei at many different densities, while making sure to separate isovector contributions from bulk properties. Lastly, we use our results to extract the static density-density linear response function of neutron matter at different densities. Our findings provide insights into inhomogeneous neutron matter and are related to the physics of neutron-...

  11. Ab initio study of the electronic properties of the planar Ga5N5 cluster

    Zheng Hao-Ping; Hao Jing-An

    2005-01-01

    The first-principles, all electron, ab initio calculations have been performed for an the amazing stable planar structure of Ga5N5 cluster based on the density functional theory. Electronic structure, Electron affinity, ionization potential, and binding energy are obtained. No spin magnetic moment is found. The results show that the planar structure of the Ga5N5 cluster is stable. It is found that for the planar structure of Ga5N5 cluster, three nitrogen atoms in the N3 subunit bind together with large electon transfer although no free N3 can exist. This may be important to the stability of the planar structure of the Ga5N5 cluster which has the lowest ground-state energy.

  12. Ab initio study on the dynamics of furfural at the liquid-solid interfaces

    Dang, Hongli; Xue, Wenhua; Shields, Darwin; Liu, Yingdi; Jentoft, Friederike; Resasco, Daniel; Wang, Sanwu

    2013-03-01

    Catalytic biomass conversion sometimes occurs at the liquid-solid interfaces. We report ab initio molecular dynamics simulations at finite temperatures for the catalytic reactions involving furfural at the water-Pd and water-Cu interfaces. We found that, during the dynamic process, the furan ring of furfural prefers to be parallel to the Pd surface and the aldehyde group tends to be away from the Pd surface. On the other hand, at the water-Cu(111) interface, furfural prefers to be tilted to the Cu surface while the aldehyde group is bonded to the surface. In both cases, interaction of liquid water and furfural is identified. The difference of dynamic process of furfural at the two interfaces suggests different catalytic reaction mechanisms for the conversion of furfural, consistent with the experimental investigations. Supported by DOE (DE-SC0004600). Simulations and calculations were performed on XSED's and NERSC's supercomputers

  13. AB INITIO STUDY ON VALENCE INDICES AND REACTIVITIES OF SOME BORAENS

    曹阳; 王友良

    1991-01-01

    In the pressnt paper, 3-21G ab initio molecular orbital calculations arc performed on diborane B2H6 and the substituted bridged-atom species H4B2X2(X=F, C1, OH, NH2, CH3),and these geometries are optimized with the energy gradient technique. According to the quantum chemical definition of atomic valence, the valences of the bridged-atoms are calculated to studtd the characteristion of the bridging bond B-X-B. Some larger boranes B4H10, B5H9, and B5H11 are also calculated to discuss the valence indices and analyze the reactiveities of the bridged-atoms.

  14. Ab-initio Hartree-Fock study of tritium desorption from Li{sub 2}O

    Taniguchi, Masaki; Tanaka, Satoru [Tokyo Univ. (Japan). Faculty of Engineering

    1998-03-01

    Dissociative adsorption of hydrogen on Li{sub 2}O (110) surface has been investigated with ab-initio Hartree-Fock quantum chemical calculation technique. Heat of adsorption and potential energy surface for H{sub 2} dissociative adsorption was evaluated by calculating the total energy of the system. Calculation results on adsorption heat indicated that H{sub 2} adsorption is endothermic. However, when oxygen vacancy exists adjacent to the adsorption sites, heat of adsorption energy became less endothermic and the activation energy required to dissociate the H-H bonding was smaller than that for the terrace site. This is considered to be caused by the excess charge localized near the defect. (author)

  15. Ab initio study of beryllium-decorated fullerenes for hydrogen storage

    Lee, Hoonkyung; Huang, Bing; Duan, Wenhui; Ihm, Jisoon

    2010-04-01

    We have found that a beryllium (Be) atom on nanostructured materials with H2 molecules generates a Kubas-like dihydrogen complex [Lee, Huang, Duan, and Ihm, Appl. Phys. Lett. 96, 143120 (2010)]. Here, we investigate the feasibility of Be-decorated fullerenes for hydrogen storage using ab initio calculations. We find that the aggregation of Be atoms on pristine fullerenes is energetically preferred, resulting in the dissociation of the dihydrogen. In contrast, for boron (B)-doped fullerenes, Be atoms prefer to be individually attached to B sites of the fullerenes, and a maximum of one H2 molecule binds to each Be atom in a form of dihydrogen with a binding energy of ˜0.3 eV. Our results show that individual dispersed Be-decorated B-doped fullerenes can serve as a room-temperature hydrogen storage medium.

  16. Communication: Ab initio study of O4H+: A tracer molecule in the interstellar medium?

    The structure and energetics of the protonated molecular oxygen dimer calculated via ab initio methods is reported. We find structures that share analogies with the eigen and zundel forms for the protonated water dimer although the symmetrical sharing of the proton is more prevalent. Analysis of different fragmentation channels show charge transfer processes which indicate the presence of conical intersections for various states including the ground state. An accurate estimate for the proton affinity of O4 leads to a significantly larger value (5.6 eV) than for O2 (4.4 eV), implying that the reaction H3+ + O4 → O4H+ + H2 is exothermic by 28 Kcal/mol as opposed to the case of O2 which is nearly thermoneutral. This opens up the possibility of using O4H+ as a tracer molecule for oxygen in the interstellar medium

  17. Charge Transfer in FeO: A combined Molecular-Dynamics and Ab Initio Study

    Molecular dynamics simulations and ab initio electronic structure calculations were carried out to determine the rate of charge transfer in stoichiometric w-stite (FeO). The charge transfer of interest occurs by II/III valence interchange between nearest-neighbor Fe atoms, with the Fe(III) constituting a ''hole'' electronic defect. There are two possible nearest-neighbor charge transfers in the FeO lattice, which occur between edge-sharing or corner-sharing FeO6 octahedra. Molecular dynamics simulations predict charge transfer rates of 3.7x1011 and 1.9x109 s-1 for the edge and corner transfers, respectively, in good agreement with those calculated using an ab initio cluster approach (1.6x1011 and 8.0x108 s-1, respectively). The calculated rates are also similar to those along basal and c-axis directions in hematite (?-Fe2O3) determined previously. Therefore, as is the case for hematite, w-stite is predicted to show anisotropic electrical conductivity. Our findings indicate that a rigid ion model does not give acceptable results, thus showing the need to account for the change in polarizability of the system upon charge transfer. Our model achieves this by using a simple mechanical shell model. By calculating the electronic coupling matrix elements for many transition state configurations obtained from the molecular dynamics simulations, we found evidence that the position of the bridging oxygen atoms can greatly affect the amount electronic coupling between the donor and acceptor states. Finally, we address the effect of oxygen vacancies on the charge transfer. It was found that an oxygen vacancy not only creates a driving force for holes to transport away from the vacancy (or equivalently for electrons to diffuse toward the vacancy) but also lowers the free energy barriers for charge transfer. In addition, the reorganization energy significantly differed from the non-defective case in a small radius around the defect

  18. Ionization dynamics of the water trimer: A direct ab initio MD study

    Tachikawa, Hiroto, E-mail: hiroto@eng.hokudai.ac.jp [Division of Materials Chemistry, Graduate School of Engineering, Hokkaido University, Kita-ku, Sapporo 060-8628 (Japan); Takada, Tomoya [Department of Material Chemistry, Asahikawa National College of Technology, Syunkodai, Asahikawa 071-8142 (Japan)

    2013-03-29

    Highlights: ► We calculated ionization dynamics of water trimer. ► Direct ab initio molecular dynamics (MD) method is applied. ► Proton transfer dynamics were discussed. ► The proton transfer process calculated are well reproduced in recent experiments. - Abstract: Ionization dynamics of the cyclic water trimer (H{sub 2}O){sub 3} have been investigated by means of direct ab initio molecular dynamics (AIMD) method. Two reaction channels, complex formation and OH dissociation, were found following the ionization of (H{sub 2}O){sub 3}. In both channels, first, a proton was rapidly transferred from H{sub 2}O{sup +} to H{sub 2}O (time scale is ∼15 fs after the ionization). In complex channel, an ion–radical contact pair (H{sub 3}O{sup +}–OH) solvated by the third water molecule was formed as a long-lived H{sub 3}O{sup +}(OH)H{sub 2}O complex. In OH dissociation channel, the second proton transfer further takes place from H{sub 3}O{sup +}(OH) to H{sub 2}O (time scale is 50–100 fs) and the OH radical is separated from the H{sub 3}O{sup +}. At the same time, the OH dissociation takes place when the excess energy is efficiently transferred into the kinetic energy of OH radical. The OH dissociation channel is significantly minor, and almost all product channels were the complex formation. The reaction mechanism was discussed on the basis of theoretical results.

  19. Modelling of interactions between volatile anaesthetics (halothane, enflurane) and aromatic compounds, ab initio study

    Graphical abstract: Halothane and enflurane form the blue-shifting hydrogen bonded C-H...complexes with benzene. The CCSD(T)/CBS calculated stabilization energies are about -10 kcal mol-1. - Abstract: For many years halothane and enflurane have been used clinically as volatile anaesthetics, however, their mechanism of action is still not fully understood. Recently, it has been suggested that they can act by a direct bonding to neuroreceptors containing the aromatic groups. In this work, the halothane...benzene and enflurane...benzene complexes were studied by the ab initio MP2 and CCSD(T) methods. All possible structures of the complexes were calculated by means of the counterpoise CP-corrected gradient optimization technique. It has been found that among these species, the C-H...π hydrogen bonded complexes are the most stable. The CCSD(T)/CBS calculated stabilization energies for halothane and enflurane complexes are: -10.56 and -9.72 kcal mol-1, respectively. The interaction energy is mainly dominated by the dispersion attraction. In the case of enflurane, the C-H bond shows a very small contraction (by -0.0008 A) upon complexation. This change is accompanied by the blue-shift (20 cm-1) of the C-H stretching frequency and an increase of the infrared intensity of the corresponding mode by 7 km mol-1. Similar results were obtained for the halothane complex: a small contraction of the C-H bond; an increase of the C-H stretching frequency by 11 cm-1 (blue-shift); and an increase of the infrared intensity by 37 km mol-1. In order to explain the nature of these effects, the halothane and enflurane molecules were studied in the electric field generated by benzene atoms, and Natural Bond Orbital (NBO) analyses were performed. The molecular dipole moments of these molecules were calculated with respect to the C-H bond changes. The positive dipole moment derivative obtained for halothane is in agreement with the literature data, while, in the case of enflurane, an unusual

  20. Stacking stability of MoS2 bilayer: An ab initio study

    Tao, Peng; Guo, Huai-Hong; Yang, Teng; Zhang, Zhi-Dong

    2014-10-01

    The study of the stacking stability of bilayer MoS2 is essential since a bilayer has exhibited advantages over single layer MoS2 in many aspects for nanoelectronic applications. We explored the relative stability, optimal sliding path between different stacking orders of bilayer MoS2, and (especially) the effect of inter-layer stress, by combining first-principles density functional total energy calculations and the climbing-image nudge-elastic-band (CI-NEB) method. Among five typical stacking orders, which can be categorized into two kinds (I: AA, AB and II: AA', AB', A'B), we found that stacking orders with Mo and S superposing from both layers, such as AA' and AB, is more stable than the others. With smaller computational efforts than potential energy profile searching, we can study the effect of inter-layer stress on the stacking stability. Under isobaric condition, the sliding barrier increases by a few eV/(ucGPa) from AA' to AB', compared to 0.1 eV/(ucGPa) from AB to [AB]. Moreover, we found that interlayer compressive stress can help enhance the transport properties of AA'. This study can help understand why inter-layer stress by dielectric gating materials can be an effective means to improving MoS2 on nanoelectronic applications.

  1. Ab initio studies on the adsorption and implantation of Al and Fe to nitride materials

    Riedl, H., E-mail: helmut.riedl@tuwien.ac.at [Christian Doppler Laboratory for Application Oriented Coating Development at the Institute of Materials Science and Technology, TU Wien, A-1040 Wien (Austria); Zálešák, J. [Erich Schmid Institute for Materials Science, Austria Academy of Science, A-8700 Leoben (Austria); Department of Physical Metallurgy and Materials Testing, Montanuniversität Leoben, A-8700 Leoben (Austria); Arndt, M. [Oerlikon Balzers, Oerlikon Surface Solutions AG, LI-9496 Balzers (Liechtenstein); Polcik, P. [Plansee Composite Materials GmbH, D-86983 Lechbruck am See (Germany); Holec, D. [Christian Doppler Laboratory for Application Oriented Coating Development at the Institute of Materials Science and Technology, TU Wien, A-1040 Wien (Austria); Department of Physical Metallurgy and Materials Testing, Montanuniversität Leoben, A-8700 Leoben (Austria); Mayrhofer, P. H. [Christian Doppler Laboratory for Application Oriented Coating Development at the Institute of Materials Science and Technology, TU Wien, A-1040 Wien (Austria); Institute of Materials Science and Technology, TU Wien, A-1040 Vienna (Austria)

    2015-09-28

    The formation of transfer material products on coated cutting and forming tools is a major failure mechanism leading to various sorts of wear. To describe the atomistic processes behind the formation of transfer materials, we use ab initio to study the adsorption energy as well as the implantation barrier of Al and Fe atoms for (001)-oriented surfaces of TiN, Ti{sub 0.50}Al{sub 0.50}N, Ti{sub 0.90}Si{sub 0.10}N, CrN, and Cr{sub 0.90}Si{sub 0.10}N. The interactions between additional atoms and nitride-surfaces are described for pure adhesion, considering no additional stresses, and for the implantation barrier. The latter, we simplified to the stress required to implant Al and Fe into sub-surface regions of the nitride material. The adsorption energies exhibit pronounced extrema at high-symmetry positions and are generally highest at nitrogen sites. Here, the binary nitrides are comparable to their ternary counterparts and the average adhesive energy is higher (more negative) on CrN than TiN based systems. Contrary, the implantation barrier for Al and Fe atoms is higher for the ternary systems Ti{sub 0.50}Al{sub 0.50}N, Ti{sub 0.90}Si{sub 0.10}N, and Cr{sub 0.90}Si{sub 0.10}N than for their binary counterparts TiN and CrN. Based on our results, we can conclude that TiN based systems outperform CrN based systems with respect to pure adhesion, while the Si-containing ternaries exhibit higher implantation barriers for Al and Fe atoms. The data obtained are important to understand the atomistic interaction of metal atoms with nitride-based materials, which is valid not just for machining operations but also for any combination such as interfaces between coatings and substrates or multilayer and phase arrangements themselves.

  2. Ab initio study of different structures of CaC: Magnetism, bonding, and lattice dynamics

    Nourbakhsh, Zahra, E-mail: z.nourbakhsh@ph.iut.ac.ir; Hashemifar, S. Javad, E-mail: hashemifar@cc.iut.ac.ir; Akbarzadeh, Hadi

    2013-12-05

    Highlights: •Electronic structure and phonon calculations are performed on seven CaC structures. •The more covalent structures have lower energies and nonmagnetic ground state. •Ferromagnetism in the ionic phases comes from the sharp C p band at the Fermi level. •Tendency of C atoms for dimerization may lead to structural instabilities in CaC. •Nonmagnetic B33 CaC is stable in a wide range of temperatures and pressures. -- Abstract: On the basis of ab initio pseudopotential calculations, we study structural, magnetic, dynamical, and mechanical properties of the hypothetical CaC ionic compound in the rock-salt (RS), B2, zinc-blende (ZB), wurtzite (WZ), NiAs (NA), anti-NiAs (NA{sup ∗}), and CrB (B33) structures. It is argued that the ZB, WZ, NA, and RS structures are more ionic while the NA{sup ∗}, B2, and B33 structures are more covalent systems. As a result of that, the nonmagnetic B33–CaC is the energetically preferred system, while the more ionic structures prefer a ferromagnetic ground state with high Fermi level spin polarization. The observed ferromagnetism in the more ionic systems is attributed to the sharp partially filled p states of carbon atom in the system. In the framework of density functional perturbation theory, the phonon spectra of these systems are computed and the observed dynamical instabilities of the NA{sup ∗} and B2 structures are explained in terms of the covalent bonds between carbon atoms. The calculated Helmholtz and Enthalpy free energies indicate the highest stability of the B33 structure in a wide range of temperatures and pressures. Among the ferromagnetic structures, RS–CaC and ZB–CaC are reported, respectively, to be the most and the least metastable systems in various thermodynamics conditions. Several mechanical properties of the dynamically stable structures of CaC are determined from their phonon spectra.

  3. Intermolecular interactions of trifluorohalomethanes with Lewis bases in the gas phase: An ab initio study

    Wang, Yi-Siang; Yin, Chih-Chien; Chao, Sheng D., E-mail: sdchao@spring.iam.ntu.edu.tw [Institute of Applied Mechanics, National Taiwan University, Taipei 106, Taiwan (China)

    2014-10-07

    We perform an ab initio computational study of molecular complexes with the general formula CF{sub 3}X—B that involve one trifluorohalomethane CF{sub 3}X (X = Cl or Br) and one of a series of Lewis bases B in the gas phase. The Lewis bases are so chosen that they provide a range of electron-donating abilities for comparison. Based on the characteristics of their electron pairs, we consider the Lewis bases with a single n-pair (NH{sub 3} and PH{sub 3}), two n-pairs (H{sub 2}O and H{sub 2}S), two n-pairs with an unsaturated bond (H{sub 2}CO and H{sub 2}CS), and a single π-pair (C{sub 2}H{sub 4}) and two π-pairs (C{sub 2}H{sub 2}). The aim is to systematically investigate the influence of the electron pair characteristics and the central atom substitution effects on the geometries and energetics of the formed complexes. The counterpoise-corrected supermolecule MP2 and coupled-cluster single double with perturbative triple [CCSD(T)] levels of theory have been employed, together with a series of basis sets up to aug-cc-pVTZ. The angular and radial configurations, the binding energies, and the electrostatic potentials of the stable complexes have been compared and discussed as the Lewis base varies. For those complexes where halogen bonding plays a significant role, the calculated geometries and energetics are consistent with the σ-hole model. Upon formation of stable complexes, the C–X bond lengths shorten, while the C–X vibrational frequencies increase, thus rendering blueshifting halogen bonds. The central atom substitution usually enlarges the intermolecular bond distances while it reduces the net charge transfers, thus weakening the bond strengths. The analysis based on the σ-hole model is grossly reliable but requires suitable modifications incorporating the central atom substitution effects, in particular, when interaction components other than electrostatic contributions are involved.

  4. An ab initio molecular dynamics study of the roaming mechanism of the H{sub 2}+HOC{sup +} reaction

    Yu Huagen, E-mail: hgy@bnl.gov [Department of Chemistry, Brookhaven National Laboratory, Upton, NY 11973 (United States)

    2011-08-01

    We report here a direct ab initio molecular dynamics study of the p-/o-H{sub 2}+HOC{sup +} reaction on the basis of the accurate SAC-MP2 potential energy surface. The quasi-classical trajectory method was employed. This work largely focuses on the study of reaction mechanisms. A roaming mechanism was identified for this molecular ion-molecule reaction. The driving forces behind the roaming mechanism were thoroughly investigated by using a trajectory dynamics approach. In addition, the thermal rate coefficients of the H{sub 2}+HOC{sup +} reaction were calculated in the temperature range [25, 300] K and are in good agreement with experiments.

  5. NO sub 3 , the study of molecular properties and photodissociation by ab initio method, spectroscopy, and translational spectroscopy

    Kim, B.

    1990-10-01

    This report discusses the following topics: molecular structure of NO{sub 3} radical studied by laser induced fluorescence; photodissociation and fluorescence spectroscopy of NO{sub 3} in molecular beam; vertical electronic spectrum of NO{sub 3}:{sup 2}A{prime}{sub 2}, {sup 2}E{double prime}({sup 2}A{sub 2}{sup 2}B{sub 1}), and {sup 2}E{prime} states; and Ab initio study of the vibrational spectra of NO{sub 3}.

  6. Theoretical studies of H2--H2 collisions. V. Ab initio calculations of relaxation phenomena in parahydrogen gas

    The temperature dependence of effective Waldmann--Snider cross sections determining relaxation and line broadening phenomena has been studied for p-H2 between 20 and 200 K. In particular, the rotational relaxation cross section and the relaxation cross sections of the rotational angular momentum vector and tensor polarizations and their respective fluxes have been calculated in an entirely ab initio treatment and close coupling formalism. As far as experimental results were available, quantitative agreement has been obtained. The rotational relaxation cross section has been proven to be quite sensitive to the potential anisotropy. Furthermore, the validity of some approximate relations between effective cross sections has been tested

  7. Ab initio study of the optical phonons in 1D antiferromagnet Ca2CuO3

    Hoang, Nam Nhat; Nguyen, Thu Hang; Nguyen, Chau

    2008-01-01

    In the spin 1/2 antiferromagnetic Ca2CuO3 the optical excitation along 1D a Cu-O chain showed the enrichment of forbidden peaks that could not be explained by the group theory. We present the cluster-model ab initio study of these optical phonons based on the Hartree-Fock SCF calculation with the 3-21G basic set. The obtained results showed very good agreement with the observed data. The Cu-O resonances generally showed the lower shifts in Ca2CuO3 than in pure CuO and were primarily composed ...

  8. Chemisorption of group-III metals on the Si(111) and Ge(111) surfaces: An ab initio study

    Ricart, J M; Rubio Martínez, Jaime; Illas i Riera, Francesc

    1990-01-01

    Chemisorption of group-III metal adatoms on Si(111) and Ge(111) has been studied through the ab initio Hartree-Fock method including nonempirical pseudopotentials and using cluster models to simulate the surface. Three different high-symmetry sites (atop, eclipsed, and open) have been considered by using X4H9, X4H7, and X6H9 (X=Si,Ge) cluster models. In a first step, ideal surface geometries have been used. Metal-induced reconstruction upon chemisorption has also been taken into account. Equi...

  9. Ab initio study of vibronic transitions between x2π and 12Σ+ electronic states of HCP+ ion

    Stojanović Ljiljana

    2013-01-01

    Full Text Available The ground and low-lying excited doublet electronic states of the HCP+ ion were studied by means of multireference configuration interaction method. Vibronic energy levels of the X2Π state of Σ, Π, Δ, and Φ symmetry, up to the 2500 cm-1, have been calculated variationally, employing previously developed ab initio methods which take into account vibronic and spin-orbit interactions. Obtained vibronic wave functions were used to estimate transition moments between vibronic energy levels of the X2Π and 12Σ+ electronic states. Results were compared to available experimental and theoretical data. [Projekat Ministarstva nauke Republike Srbije, br. 172040

  10. Mechanism of prooxidant reaction of vitamin E. Kinetic, spectroscopic, and ab initio study of proton-transfer reaction

    Nagoaka, Shin-ichi; Sawada, Kouhei; Fukumoto, Youji [Ehime Univ., Matsuyama (Japan)] [and others

    1992-08-06

    This paper discusses using kinetic, spectroscopic, and ab initio studies of the prooxidant reaction of vitamin E derivatives to determine second-order rate constants for the reaction of six tocopheroxyl radicals with five alkyl hydroperoxides in benzene along with the first adiabatic ionization potentials of the alkyl hydroperoxides. These results suggest that charge transfer and proton tunneling play important roles in this prooxidant reaction, with proton tunneling taking place below the transition state and allowing the proton to cut a corner on the potential energy surface. 31 refs., 10 figs., 2 tabs.

  11. Ab initio study of ferromagnetic La0.5Ba0.5CoO3

    Umesh V Waghmare

    2003-10-01

    We study structure and magnetic properties of La0.5Ba0.5CoO3 (LBCO) using ab initio density functional theory (DFT) method based on pseudopotentials and a plane-wave basis. We find the cubic structure of LBCO is ferromagnetic and lower in energy than the nonmagnetic rhombohedral structure. Through the calculation of -point phonons of LBCO in the cubic structure, we determine its structural instabilities and find that they correspond to the structural phase transition accompanying a para-ferromagnetic transition observed recently.

  12. Stacking stability of MoS2 bilayer: An ab initio study

    The study of the stacking stability of bilayer MoS2 is essential since a bilayer has exhibited advantages over single layer MoS2 in many aspects for nanoelectronic applications. We explored the relative stability, optimal sliding path between different stacking orders of bilayer MoS2, and (especially) the effect of inter-layer stress, by combining first-principles density functional total energy calculations and the climbing-image nudge-elastic-band (CI-NEB) method. Among five typical stacking orders, which can be categorized into two kinds (I: AA, AB and II: AA', AB', A'B), we found that stacking orders with Mo and S superposing from both layers, such as AA' and AB, is more stable than the others. With smaller computational efforts than potential energy profile searching, we can study the effect of inter-layer stress on the stacking stability. Under isobaric condition, the sliding barrier increases by a few eV/(ucGPa) from AA' to AB', compared to 0.1 eV/(ucGPa) from AB to [AB]. Moreover, we found that interlayer compressive stress can help enhance the transport properties of AA'. This study can help understand why inter-layer stress by dielectric gating materials can be an effective means to improving MoS2 on nanoelectronic applications. (condensed matter: structural, mechanical, and thermal properties)

  13. High adsorption capacity of heavy metals on two-dimensional MXenes: an ab initio study with molecular dynamics simulation.

    Guo, Xun; Zhang, Xitong; Zhao, Shijun; Huang, Qing; Xue, Jianming

    2016-01-01

    Density functional theory (DFT) calculation is employed to study the adsorption properties of Pb and Cu on recently synthesized two-dimensional materials MXenes, including Ti3C2, V2C1 and Ti2C1. The influence of surface decoration with functional groups such as H, OH and F have also been investigated. Most of these studied MXenes exhibit excellent capability to adsorb Pb and Cu, especially the adsorption capacity of Pb on Ti2C1 is as high as 2560 mg g(-1). Both the binding energies and the adsorption capacities are sensitive to the functional groups attached to the MXenes' surface. Ab initio molecular dynamics (ab-init MD) simulation confirms that Ti2C1 remains stable at room temperature after adsorbing Pb atoms. Our calculations imply that these newly emerging two-dimensional MXenes are promising candidates for wastewater treatment and ion separation. PMID:26602974

  14. Carbonate-Bridged Lanthanoid Triangles: Single-Molecule Magnet Behavior, Inelastic Neutron Scattering, and Ab Initio Studies.

    Giansiracusa, Marcus J; Vonci, Michele; Van den Heuvel, Willem; Gable, Robert W; Moubaraki, Boujemaa; Murray, Keith S; Yu, Dehong; Mole, Richard A; Soncini, Alessandro; Boskovic, Colette

    2016-06-01

    Optimization of literature synthetic procedures has afforded, in moderate yield, homogeneous and crystalline samples of the five analogues Na11[{RE(OH2)}3CO3(PW9O34)2] (1-RE; RE = Y, Tb, Dy, Ho, and Er). Phase-transfer methods have allowed isolation of the mixed salts (Et4N)9Na2[{RE(OH2)}3CO3(PW9O34)2] (2-RE; RE = Y and Er). The isostructural polyanions in these compounds are comprised of a triangular arrangement of trivalent rare-earth ions bridged by a μ3-carbonate ligand and sandwiched between two trilacunary Keggin {PW9O34} polyoxometalate ligands. Alternating-current (ac) magnetic susceptibility studies of 1-Dy, 1-Er, and 2-Er reveal the onset of frequency dependence for the out-of-phase susceptibility in the presence of an applied magnetic field at the lowest measured temperatures. Inelastic neutron scattering (INS) spectra of 1-Ho and 1-Er exhibit transitions between the lowest-lying crystal-field (CF) split states of the respective J = 8 and (15)/2 ground-state spin-orbit multiplets of the Ho(III) and Er(III) ions. Complementary ab initio calculations performed for these two analogues allow excellent reproduction of the experimental magnetic susceptibility and low-temperature magnetization data and are in reasonable agreement with the experimental INS data. The ab initio calculations reveal that the slight difference in coordination environments of the three Ln(III) ions in each complex gives rise to differences in the CF splitting that are not insignificant. This theoretical result is consistent with the observation of multiple relaxation processes by ac magnetic susceptibility and the broadness of the measured INS peaks. The ab initio calculations also indicate substantial mixing of the MJ contributions to the CF split energy levels of each Ln(III) ion. Calculations indicate that the CF ground states of the Ho(III) centers in 1-Ho are predominantly comprised of contributions from small MJ, while those of the Er(III) centers in 1-Er are predominantly

  15. Ab initio valence calculations in chemistry

    Cook, D B

    1974-01-01

    Ab Initio Valence Calculations in Chemistry describes the theory and practice of ab initio valence calculations in chemistry and applies the ideas to a specific example, linear BeH2. Topics covered include the Schrödinger equation and the orbital approximation to atomic orbitals; molecular orbital and valence bond methods; practical molecular wave functions; and molecular integrals. Open shell systems, molecular symmetry, and localized descriptions of electronic structure are also discussed. This book is comprised of 13 chapters and begins by introducing the reader to the use of the Schrödinge

  16. Ab initio study of intrinsic defects and diffusion behaviors in solid molecular hydrogens

    Sun, Qingqiang; Ye, Yingting; Yang, Tianle; Yang, Li; Peng, Shuming; Long, Xinggui; Zhou, Xiaosong; Zu, Xiaotao; Du, Jincheng

    2015-12-01

    The behaviors of intrinsic defects in solid molecular hydrogens (H2) were investigated using ab initio calculations based on density functional theory. The results show that the formation energy of a vacancy is dependent on molecule orientation in disordered hexagonal close-packed (hcp) H2 crystals, but independent of molecular orientation in face-centered cubic-Pa3 H2. Furthermore, H2 molecules generally prefer to occupy the basal octahedral sites in hcp and octahedral sites in Pa3 structures. The formation energies of an interstitial H2 depend sensitively on the volume of interstitial sites, and also on near spatial distributions of molecular axes of a H2 in hcp H2 crystals. The strong force field introduced by an interstitial H2 might induce the rotation of molecular axes and reduce the formation energy of the interstitial H2. The migration barrier energy of a vacancy in Pa3 is larger than the average in hcp structures. However, the H2 molecules prefer to jump though vacancies in the two structures of solid H2 than the interstitial migrations.

  17. Anomalous structural dynamics in liquid Al80Cu20: An ab initio molecular dynamics study

    In this work, the structural dynamics of liquid Al80Cu20 is systematically investigated in terms of the evolution of its atomic structure, diffusivity, viscosity and fragility through ab initio molecular dynamics simulations. In addition to using pair correlation functions and coordination numbers, the various local ordered clusters are characterized comprehensively by Honeycutt-Anderson bond pairs and Voronoi polyhedra. Compared to the self diffusivity of pure liquid Cu, the tracer diffusion coefficients of Cu in liquid Al80Cu20 are increased, in agreement with the results measured by quasielastic neutron-scattering (QENS). Although the interdiffusion coefficients predicted by Darken’s equation match well to those obtained from the viscosity measurements via the Stokes-Einstein relation, they are smaller than those measured by QENS or X-ray radiography, indicative of an anomalous nature of the structural dynamics, dominated by the local ordered clusters in liquid Al80Cu20. Furthermore, Vogel–Fulcher–Tammann fitting results indicate that the liquid Al80Cu20 can be classified into a strong liquid. The deformation electron density shows that the intrinsic tetrahedral-type bonds in FCC Al and Cu are transformed into an amorphous type in liquid Al80Cu20. The present work provides insights into the understanding of structural dynamics and the kinetic properties of such metallic melts

  18. Ab initio studies of transition metal complexes and related electron transfer properties

    Ab initio electronic structure calculations have been carried out for various aquo and ammine complexes of Fe, Co, and Ru in their 2+ and 3+ oxidation states. The results of these calculations are used as assessing a variety of factors controlling electron transfer kinetics including the charge-state dependence of various inner shell geometrical parameters and associated vibrational frequencies (metal-ligand and intra-ligand and intra-ligand stretching modes, and librational modes), and the dependence of electron transfer matrix elements on the nature of the ligand and the metal orbital type (/sup t/2/sub g/ vs. e/sub g/). The charge-state dependence of the OH bond lengths in hexa-aquo ions is predicted to yield H/D isotope effects (i.e., fractionation between bulk and first shell water), which should be detectable by neutron scattering experiments on appropriate aqueous solutions. The relationship between these thermodynamic isotope effects, which depend critically on strong OH triple bond O hydrogen bonding between first and second-shell water molecules, and kinetic H/D isotope effects in electron exchange involving hexa-aquo ions will be discussed

  19. Ab-initio study of the dielectric response of high-permittivity perovskites for energy storage

    Many of materials based on transition metals have a wide range of applications, such as the storage of energy, due to their peculiar properties (high-dielectric constants, ferro-electricity,...). The knowledge of their bulk properties is essential in designing targeted devices with high performance. For instance, ABO3 perovskites are peculiarly interesting for their atomic structural flexibility, allowing high number of atoms substitution and giving them specific chemical and electrical properties compared to the pure compounds. In this context, first principles calculations can be useful to understand the structural and electronic properties of these materials. The pressure-induced giant dielectric anomaly of ABO3 perovskites has been investigated at the ab initio level. Its mechanism has been analyzed in terms of thermodynamic phase stability, structural and phonon contributions and Born effective charges. It is shown that the IR-active soft phonon is responsible for the anomaly. This mode always involves a displacement and a deformation of the oxygen octahedra, while the roles of A and B ions vary among the materials and between high- and low-pressure phase transitions. A sharp increase in the phonon amplitude near the phase transition gives rise to the dielectric anomaly. The use of hybrid functionals is required for agreement with experimental data. The calculations show that the dielectric anomaly in the pressure-induced phase transitions of these perovskites is a property of the bulk material. (author)

  20. Ab initio Mechanism Study on the Reaction of Chlorine Atom with Formic Acid

    于海涛; 付宏刚; 等

    2003-01-01

    The potential energy surface(PES) for the reaction of Cl atom with HCOOH is predicted using ab initio molecular orbital calculation methods at UQCIDS(T,full)6-311++G(3df,2p)//UMP2(full)/6-311+G(d,P) level of theory with zero-point vibrational energy (ZPVE) correction.The calculated results show that the reaction mechanism of Cl atom with formic acid is a C-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom with a 3.73kJ/mol reaction barrier height,leading to the formation of cis-HOCO radical which will reacts with Cl atom or other molecules in such a reaction system.Because the reaction barrier height of O-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom which leads to the formation of HCO2 radical is 67.95kJ/mol,it is a secondary reaction channel in experiment,This is in good agreement with the prediction based on the previous experiments.

  1. An ab initio study of magneto-electric coupling of YMnO3

    This paper proposes the direct calculation of the microscopic contributions to the magneto-electric coupling, using ab initio methods. The electrostrictive and the Dzyaloshinskii–Moriya contributions were evaluated individually. For this purpose a specific method was designed, combining density functional theory calculations and embedded fragment, explicitly correlated, quantum chemical calculations. This method allowed us to calculate the evolution of the magnetic couplings as a function of an applied electric field. We found that in YMnO3 the Dzyaloshinskii–Moriya contribution to the magneto-electric effect is three orders of magnitude weaker than the electrostrictive contribution. Strictive effects are thus dominant in the magnetic exchange evolution under an applied electric field, and by extension in the magneto-electric effect. These effects however, remain quite small, and the modifications of the magnetic excitations under an applied electric field will be difficult to observe experimentally. Another important conclusion is that it can be shown that the linear magneto-electric tensor is null due to the inter-layer symmetry operations. (paper)

  2. Ab initio energetic study of oxide ceramics with rare-earth elements

    WU Bo; Matvei Zinkevich; WANG Chong; Fritz Aldinger

    2006-01-01

    Ab initio energetic calculations based on the density functional theory (DFT) and the projector augmented wave method (PAW) for determining the polymorphisms of lanthanide sesquioxides Ln2O3 (where Ln = rare-earth element. Y,and Sc), LnMO3 perovskites (where M = Al and Ga), and Ln2B2O7 pyrochlores (where B = Ti, Zr, and Hf) were reported. The relative lattice stabilities agreed well with the critically assessed results or the experimental results except the C-type Ln2O3 with a cubic structure, for which the calculated total energies were considerably more negative. With the increase of the Ln3+-cation radius, the polymorphic structures showed a degenerative tendency. The tendencies and quantities of the enthalpies of formation of the ternary oxide ceramics synthesized from their constituent binary oxides reasonably agreed with the available experimental results, and valuable thermodynamic properties were afforded to the compound, for which no experimental data is available. The enthalpies of formation of both perovskites and pyrochlores tend to become more negative with the increase of the Ln3+-cation radius.

  3. Switchable magnetic moment in cobalt-doped graphene bilayer on Cu(111): An ab initio study

    Souza, Everson S.; Scopel, Wanderlã L.; Miwa, R. H.

    2016-06-01

    In this work, we have performed an ab initio theoretical investigation of substitutional cobalt atoms in the graphene bilayer supported on the Cu(111) surface (Co/GBL/Cu). Initially, we examined the separated systems, namely, graphene bilayer adsorbed on Cu(111) (GBL/Cu) and a free standing Co-doped GBL (Co/GBL). In the former system, the GBL becomes n -type doped, where we map the net electronic charge density distribution along the GBL-Cu(111) interface. The substitutional Co atom in Co/GBL lies between the graphene layers, and present a net magnetic moment mostly due to the unpaired Co-3 dz2 electrons. In Co/GBL/Cu, we found that the Cu(111) substrate rules (i) the energetic stability, and (ii) the magnetic properties of substitutional Co atoms in the graphene bilayer. In (i), the substitutional Co atom becomes energetically more stable lying on the GBL surface, and in (ii), the magnetic moment of Co/GBL has been quenched due to the Cu(111) → Co/GBL electronic charge transfer. We verify that such a charge transfer can be tuned upon the application of an external electric field, and thus mediated by a suitable change on the electronic occupation of the Co-dz2 orbitals, we found a way to switch-on and -off the magnetization of the Co-doped GBL adsorbed on the Cu(111) surface.

  4. A Comparative Study of Ab-Initio Thermal Conductivity Approaches: The Case of Cubic Boron Nitride

    Mukhopadhyay, Saikat; Lindsay, Lucas; Broido, David; Stewart, Derek

    2013-03-01

    Given its high strength and large thermal conductivity, cubic boron nitride (cBN) provides an important complement to diamond films for heat spreading applications. However, cBN, in contrast to diamond, is a polar material with significant LO-TO splitting in the phonon dispersion. In this talk, we examine the lattice thermal conductivity of cBN using several approaches based on first principles calculations. These approaches include: (1) an analytic modified Callaway-Debye model that relies on parameters from ab-initio harmonic force constants, (2) a fully self-consistent calculation of the thermal conductivity that links an iterative solution of the phonon Boltzmann transport equation (BTE) with harmonic and anharmonic interatomic force constants. The force constants for the BTE are calculated using two approaches: density functional perturbation theory and a real-space supercell approach. We will compare the results from these approaches, highlight the role of normal phonon-phonon scattering, and also examine the impact of optical modes and LO-TO splitting. In addition, we will discuss how isotope scattering affects thermal conductivity and compare this to other boron nitride structures (hexagonal BN, BN sheets and BN nanotubes).

  5. Ab initio study of the unusual thermal transport properties of boron arsenide and related materials

    Broido, D. A.; Lindsay, L.; Reinecke, T. L.

    2013-12-01

    Recently, using a first principles approach, we predicted that zinc blende boron arsenide (BAs) will have an ultrahigh lattice thermal conductivity, κ, of over 2000 Wm-1K-1 at room temperature (RT), comparable to that of diamond. Here, we provide a detailed ab initio examination of phonon thermal transport in boron arsenide, contrasting its unconventional behavior with that of other related materials, including the zinc blende crystals boron nitride (BN), boron phosphide, boron antimonide, and gallium nitride (GaN). The unusual vibrational properties of BAs contribute to its weak phonon-phonon scattering and phonon-isotope scattering, which are responsible for its exceptionally high κ. The thermal conductivity of BAs has contributions from phonons with anomalously large mean free paths (˜2 μm), two to three times those of diamond and BN. This makes κ in BAs sensitive to phonon scattering from crystal boundaries. An order of magnitude smaller RT thermal conductivity in a similar material, zinc blende GaN, is connected to more separated acoustic phonon branches, larger anharmonic force constants, and a large isotope mixture on the heavy rather than the light constituent atom. The striking difference in κ for BAs and GaN demonstrates the importance of using a microscopic first principles thermal transport approach for calculating κ. BAs also has an advantageous RT coefficient of thermal expansion, which, combined with the high κ value, suggests that it is a promising material for use in thermal management applications.

  6. Ab initio atomistic thermodynamics study on the oxidation mechanism of binary and ternary alloy surfaces

    Utilizing a combination of ab initio density-functional theory and thermodynamics formalism, we have established the microscopic mechanisms for oxidation of the binary and ternary alloy surfaces and provided a clear explanation for the experimental results of the oxidation. We construct three-dimensional surface phase diagrams (SPDs) for oxygen adsorption on three different Nb-X(110) (X = Ti, Al or Si) binary alloy surfaces. On the basis of the obtained SPDs, we conclude a general microscopic mechanism for the thermodynamic oxidation, that is, under O-rich conditions, a uniform single-phase SPD (type I) and a nonuniform double-phase SPD (type II) correspond to the sustained complete selective oxidation and the non-sustained partial selective oxidation by adding the X element, respectively. Furthermore, by revealing the framework of thermodynamics for the oxidation mechanism of ternary alloys through the comparison of the surface energies of two separated binary alloys, we provide an understanding for the selective oxidation behavior of the Nb ternary alloy surfaces. Using these general microscopic mechanisms, one could predict the oxidation behavior of any binary and multi-component alloy surfaces based on thermodynamics considerations

  7. Ab initio study of Cr interactions with point defects in bcc Fe

    Full text of publication follows. Ferritic martensitic steels are candidate structural materials for fast neutron reactors, and in particular high-Cr reduced-activation steels. In Fe-Cr alloys, Cr plays a major role in the radiation-induced evolution of the mechanical properties. Using ab initio calculations based on density functional theory, the properties of Cr in α-Fe have been investigated. The intrinsic point defect formation energies were found to be larger in model bcc Cr as compared to those in ferromagnetic bcc Fe. The interactions of Cr with point defects (vacancy and self interstitials) have been characterised. Single Cr atoms interact weakly with vacancies but significantly with self-interstitial atoms. Mixed interstitials of any interstitial symmetry are bound. Configurations where two Cr atoms are in nearest neighbour position are generally unfavourable in bcc Fe except when they are a part of a interstitial complex. Mixed interstitials do not have as strong directional stability as pure Fe interstitials have. The effects on the results using the atom description scheme of either the ultrasoft pseudo-potential (USPP) or the projector augmented wave (PAW) formalisms are connected to the differences in local magnetic moments that the two methods predict. As expected for the Fe-Cr system, the results obtained using the PAW method are more reliable than the ones obtained with USPP. (authors)

  8. Laser spectroscopy and ab initio studies of metal-containing free radicals

    Greetham, G M

    2000-01-01

    strontium-containing free radical is reported, that of SrCCH. This new excited electronic state is accessed by the orbitally-forbidden B-tilde' sup 2 DELTA-X-tilde sup 2 SIGMA sup + transition. Spin-orbit and vibrational structure have been seen in spectra of SrCCH and SrCCD and confirmed the assignment. Finally, observation of a new transition in an unidentified gallium-containing molecule is reported. Two progressions corresponding to two different vibrational modes of the molecule are seen in the spectrum. Potential spectral carriers, including Ga sub x clusters and other gallium-containing molecules formed by reaction with impurities, are discussed in an attempt to explain the observed spectrum. This work describes the use of laser spectroscopy and ab initio calculations in the investigation of several new electronic transitions in metal-containing free radicals. These free radicals were prepared in a supersonic jet by laser ablation of solid metal samples in the presence of appropriate precursor molecule...

  9. Ab initio Study on Structures and Isomerization of Magnesium Fluorosilylenoid H2SiFMgF

    Yi-jian Zhang; Mei-jiang Li; Guo-qiao Lai; Da-cheng Feng; Sheng-yu Feng

    2008-01-01

    The structures and isomerization of magnesium fluorosilylenoid H2SiFMgF were investigated by ab initio molecular orbital theory. Four equilibrium structures and three isomeric transition states were located and fully optimized at the B3LYP/6-31G(d,p) and G3MP2B3 levels, respectively. Based on the B3LYP/6-31G(d,p) optimized geometries, harmonic frequencies of various structures were obtained and 29Si chemical shifts were calculated. The solvent effects were investigated by means of the polarizable continuum model using THF as a solvent at B3LYP/6-31G(d,p) level. Isomerization paths for isomers were confirmed by in-trinsic reaction coordinate calculations. The calculated results show that tetrahedral structure has the lowest energy and is the most stable; tetrahedral, three-membered ring, and p-complex structures are suggested to be the experimentally detectable ones; and σ-complex structure has the highest energy and will not exist.

  10. Ferromagnetic ordering of Cr and Fe doped p-type diamond: An ab initio study

    Ferromagnetic ordering of transition metal dopants in semiconductors holds the prospect of combining the capabilities of semiconductors and magnetic systems in single hybrid devices for spintronic applications. Various semiconductors have so far been considered for spintronic applications, but low Curie temperatures have hindered room temperature applications. We report ab initio DFT calculations on the stability and magnetic properties of Fe and Cr impurities in diamond, and show that their ground state magnetic ordering and stabilization energies depend strongly on the charge state and type of co-doping. We predict that divacancy Cr+2 and substitutional Fe+1 order ferromagnetically in p-type diamond, with magnetic stabilization energies (and magnetic moment per impurity ion) of 16.9 meV (2.5 μB) and 33.3 meV (1.0 μB), respectively. These magnetic stabilization energies are much larger than what has been achieved in other semiconductors at comparable impurity concentrations, including the archetypal dilute magnetic semiconductor GaAs:Mn. In addition, substitutional Fe+1 exhibits a strong half-metallic character, with the Fermi level crossing bands in only the spin down channel. These results, combined with diamond’s extreme properties, demonstrate that Cr or Fe dopedp-type diamond may successfully be considered in the search for room temperature spintronic materials

  11. Ab initio atomistic thermodynamics study on the oxidation mechanism of binary and ternary alloy surfaces

    Liu, Shi-Yu, E-mail: buaasyliu@gmail.com [College of Physics and Electronic Information Science, Tianjin Normal University, Tianjin 300387 (China); Department of Mechanical and Biomedical Engineering, City University of Hong Kong, Kowloon, Hong Kong (China); Liu, Shiyang [Institute of Information Optics, Zhejiang Normal University, Jinhua, Zhejiang 321004 (China); Li, De-Jun [College of Physics and Electronic Information Science, Tianjin Normal University, Tianjin 300387 (China); Wang, Sanwu, E-mail: sanwu-wang@utulsa.edu [Department of Physics and Engineering Physics, The University of Tulsa, Tulsa, Oklahoma 74104 (United States); Guo, Jing; Shen, Yaogen, E-mail: meshen@cityu.edu.hk [Department of Mechanical and Biomedical Engineering, City University of Hong Kong, Kowloon, Hong Kong (China)

    2015-02-14

    Utilizing a combination of ab initio density-functional theory and thermodynamics formalism, we have established the microscopic mechanisms for oxidation of the binary and ternary alloy surfaces and provided a clear explanation for the experimental results of the oxidation. We construct three-dimensional surface phase diagrams (SPDs) for oxygen adsorption on three different Nb-X(110) (X = Ti, Al or Si) binary alloy surfaces. On the basis of the obtained SPDs, we conclude a general microscopic mechanism for the thermodynamic oxidation, that is, under O-rich conditions, a uniform single-phase SPD (type I) and a nonuniform double-phase SPD (type II) correspond to the sustained complete selective oxidation and the non-sustained partial selective oxidation by adding the X element, respectively. Furthermore, by revealing the framework of thermodynamics for the oxidation mechanism of ternary alloys through the comparison of the surface energies of two separated binary alloys, we provide an understanding for the selective oxidation behavior of the Nb ternary alloy surfaces. Using these general microscopic mechanisms, one could predict the oxidation behavior of any binary and multi-component alloy surfaces based on thermodynamics considerations.

  12. Interaction of minor alloying elements of high-Cr ferritic steels with lattice defects: An ab initio study

    Bakaev, A. [SCK-CEN, Nuclear Materials Science Institute, Boeretang 200, Mol B2400 (Belgium); Center for Molecular Modeling, Department of Physics and Astronomy, Ghent University, Technologiepark 903, 9052 Zwijnaarde (Belgium); Department of Experimental Nuclear Physics K-89, Faculty of Physics and Mechanics, St. Petersburg State Polytechnical University, 29 Polytekhnicheskaya str., 195251 St. Petersburg (Russian Federation); Terentyev, D., E-mail: dterenty@sckcen.be [SCK-CEN, Nuclear Materials Science Institute, Boeretang 200, Mol B2400 (Belgium); Bonny, G. [SCK-CEN, Nuclear Materials Science Institute, Boeretang 200, Mol B2400 (Belgium); Klaver, T.P.C. [Department of Materials Science and Engineering, Faculty of 3mE, TU Delft, Mekelweg 2, 2628 CD Delft (Netherlands); Olsson, P. [Department of Neutron Research, Angstrom Laboratory, Uppsala University, Box 525, SE-75120 Uppsala (Sweden); Van Neck, D. [Center for Molecular Modeling, Department of Physics and Astronomy, Ghent University, Technologiepark 903, 9052 Zwijnaarde (Belgium)

    2014-01-15

    Highlights: • The interaction of Mo, W, Nb, Ta, V, Mn, Si with point and extended defects is characterized. • Mn and Si exhibit peculiar interaction with both point and extended lattice defects. • The results for substitutional atoms of the refractory metals well follow one specific trend. -- Abstract: Basic properties of minor alloying elements, namely Mo, W, Nb, Ta, V, Mn, Si entering the conventional and reduced-activation structural Fe–(9–12)Cr steels have been analyzed using ab initio calculations. The electronic structure calculations were applied to study the interaction of minor alloying elements with a number of important and well defined lattice structures, such as point defects, the 1/2〈1 1 1〉 screw dislocation core, high angle symmetric grain boundaries and free surfaces. The studied elements were classified according to their similarities and discrepancies regarding the interaction with the above mentioned defects. The refractory alloying elements are found to follow the same trend whereas Mn and Si exhibit peculiar behavior with respect to the interaction with both point and extended lattice defects. The obtained results are discussed and compared with previously published ab initio and available experimental data.

  13. Ab initio theoretical comparative study of magnetic coupling in KNiF3sand K2 NiF4s

    de P. R. Moreira, Ibério; Illas, Francesc

    1997-02-01

    The origin of magnetic coupling in KNiF3 and K2 NiF4 is studied by means of an ab initio cluster model approach. By a detailed study of the mapping between eigenstates of the exact nonrelativistic and spin model Hamiltonians it is possible to obtain the magnetic coupling constant J and to compare ab initio cluster-model values with those resulting from ab initio periodic Hartree-Fock calculations. This comparison shows that J is strongly determined by two-body interactions; this is a surprising and unexpected result. The importance of the ligands surrounding the basic metal-ligand-metal interacting unit is reexamined by using two different partitions and the constrained space orbital variation method of analysis. This decomposition enables us to show that this effect is basically environmental. Finally, dynamical electronic correlation effects have found to be critical in determining the final value of the magnetic coupling constant.

  14. A theoretical-spectroscopy, ab-initio-based study of the electronic ground state of 121SbH3

    Yurchenko, Sergei N.; Carvajal Zaera, Miguel; Yachmenev, Andrey; Thiel, Walter; Jensen, Per

    2010-01-01

    For the stibine isotopologue (SbH3)-Sb-121, we report improved theoretical calculations of the vibrational energies below 8000 cm- and simulations of the rovibrational spectrum in the 0-8000 cm(-1) region. The calculations are based on a refined ab initio potential energy surface and on a new dipole moment surface obtained at the coupled cluster CCSD(T) level. The theoretical results are compared with the available experimental data in order to validate the ab initio surfaces and the TROVE co...

  15. Ab initio studies on the spin-forbidden cooling transitions of the LiRb molecule.

    You, Yang; Yang, Chuan-Lu; Zhang, Qing-Qing; Wang, Mei-Shan; Ma, Xiao-Guang; Liu, Wen-Wang

    2016-07-20

    The spin-forbidden cooling of the LiRb molecule is investigated based on ab initio quantum chemistry calculations. The multireference configuration interaction method is used to generate the potential energy curves (PECs) of the ground state X(1)Σ(+) and the low-lying excited states a(3)Σ(+), B(1)Π, and b(3)Π. The spin-orbit coupling effects for the PECs and the transition dipole moments (TDMs) between the X(1)Σ(+), b(3)Π and a(3)Σ(+) states are also calculated. The analytical functions for the PECs are deduced. The rovibrational energy levels, the spectroscopic parameters and the Franck-Condon factors (FCF) are determined by solving the Schrödinger equation of nuclear movement with the obtained analytical functions. The b(3)Π0 ↔ X(1)Σ(+) and b(3)Π1 ↔ X(1)Σ(+) transitions have highly diagonal distributed FCFs and non-zero TDMs, demonstrating that the LiRb molecule could be a very promising candidate for laser cooling. Therefore, a three-cycle laser cooling scheme for the molecule has been proposed based on these two spin-forbidden transitions. Using the radiative lifetime and linewidth calculated from the obtained TDM functions, we present further analysis of the cooling of LiRb and the corresponding KRb molecule. The transition b(3)Π0 ↔ X(1)Σ(+) is found to be a practical transition to cool the LiRb molecule, and a sub-microkelvin cool temperature could be reached for the KRb molecule using a similar laser cooling scheme. PMID:27388722

  16. Rotational study of the CH4-CO complex: Millimeter-wave measurements and ab initio calculations.

    Surin, L A; Tarabukin, I V; Panfilov, V A; Schlemmer, S; Kalugina, Y N; Faure, A; Rist, C; van der Avoird, A

    2015-10-21

    The rotational spectrum of the van der Waals complex CH4-CO has been measured with the intracavity OROTRON jet spectrometer in the frequency range of 110-145 GHz. Newly observed and assigned transitions belong to the K = 2-1 subband correlating with the rotationless jCH4 = 0 ground state and the K = 2-1 and K = 0-1 subbands correlating with the jCH4 = 2 excited state of free methane. The (approximate) quantum number K is the projection of the total angular momentum J on the intermolecular axis. The new data were analyzed together with the known millimeter-wave and microwave transitions in order to determine the molecular parameters of the CH4-CO complex. Accompanying ab initio calculations of the intermolecular potential energy surface (PES) of CH4-CO have been carried out at the explicitly correlated coupled cluster level of theory with single, double, and perturbative triple excitations [CCSD(T)-F12a] and an augmented correlation-consistent triple zeta (aVTZ) basis set. The global minimum of the five-dimensional PES corresponds to an approximately T-shaped structure with the CH4 face closest to the CO subunit and binding energy De = 177.82 cm(-1). The bound rovibrational levels of the CH4-CO complex were calculated for total angular momentum J = 0-6 on this intermolecular potential surface and compared with the experimental results. The calculated dissociation energies D0 are 91.32, 94.46, and 104.21 cm(-1) for A (jCH4 = 0), F (jCH4 = 1), and E (jCH4 = 2) nuclear spin modifications of CH4-CO, respectively. PMID:26493903

  17. Ab initio study on stacking sequences, free energy, dynamical stability and potential energy surfaces of graphite structures

    Ab initio simulations have been performed to study the structure, energetics and stability of several plausible stacking sequences in graphite. These calculations suggest that in addition to the standard structures, graphite can also exist in AA-simple hexagonal, AB-orthorhombic and ABC-hexagonal type stacking. The free energy difference between these structures is very small (∼1 meV/atom), and hence all the structures can coexist from purely energetic considerations. Calculated x-ray diffraction patterns are similar to those of the standard structures for 2θ ⩽ 70°. Shear elastic constant C44 is negative in AA-simple hexagonal, AB-orthorhombic and ABC-hexagonal structures, suggesting that these structures are mechanically unstable. Phonon dispersions show that the frequencies of some modes along the Γ–A direction in the Brillouin zone are imaginary in all of the new structures, implying that these structures are dynamically unstable. Incorporation of zero point vibrational energy via the quasi-harmonic approximation does not result in the restoration of dynamical stability. Potential energy surfaces for the unstable normal modes are seen to have the topography of a potential hill for all the new structures, confirming that all of the new structures are inherently unstable. The fact that the potential energy surface is not in the form of a double well implies that the structures are linearly as well as globally unstable. (paper)

  18. Magnetic Transitions in Iron Porphyrin Halides by Inelastic Neutron Scattering and Ab-initio Studies of Zero-Field Splittings

    Zero-field splitting (ZFS) parameters of nondeuterated metalloporphyrins [Fe(TPP)X] (X = F, Br, I; H2TPP = tetraphenylporphyrin) are determined by inelastic neutron scattering (INS). The ZFS values are D = 4.49(9) cm-1 for tetragonal polycrystalline [Fe(TPP)F], and D = 8.8(2) cm-1, E = 0.1(2) cm-1 and D = 13.4(6) cm-1, E = 0.3(6) cm-1 for monoclinic polycrystalline [Fe(TPP)Br] and [Fe(TPP)I], respectively. Along with our recent report of the ZFS value of D = 6.33(8) cm-1 for tetragonal polycrystalline [Fe(TPP)Cl], these data provide a rare, complete determination of ZFS parameters in a metalloporphyrin halide series. The electronic structure of [Fe(TPP)X] (X = F, Cl, Br, I) has been studied by multireference ab initio methods: the complete active space self-consistent field (CASSCF) and the N-electron valence perturbation theory (NEVPT2) with the aim of exploring the origin of the large and positive zero-field splitting D of the 6A1 ground state. D was calculated from wave functions of the electronic multiplets spanned by the d5 configuration of Fe(III) along with spin-orbit coupling accounted for by quasi degenerate perturbation theory. Results reproduce trends of D from inelastic neutron scattering data increasing in the order from F, Cl, Br, to I. A mapping of energy eigenvalues and eigenfunctions of the S = 3/2 excited states on ligand field theory was used to characterize the σ- and Π-antibonding effects decreasing from F to I. This is in agreement with similar results deduced from ab initio calculations on CrX63- complexes and also with the spectrochemical series showing a decrease of the ligand field in the same directions. A correlation is found between the increase of D and decrease of the Π- and σ-antibonding energies e(λ)(X) (λ = σ, Π) in the series from X = F to I. Analysis of this correlation using second-order perturbation theory expressions in terms of angular overlap parameters rationalizes the experimentally deduced trend. Furthermore, D

  19. Magnetic Transitions in Iron Porphyrin Halides by Inelastic Neutron Scattering and Ab Initio Studies of Zero-Field Splittings.

    Stavretis, Shelby E; Atanasov, Mihail; Podlesnyak, Andrey A; Hunter, Seth C; Neese, Frank; Xue, Zi-Ling

    2015-10-19

    Zero-field splitting (ZFS) parameters of nondeuterated metalloporphyrins [Fe(TPP)X] (X = F, Br, I; H₂TPP = tetraphenylporphyrin) have been directly determined by inelastic neutron scattering (INS). The ZFS values are D = 4.49(9) cm⁻¹ for tetragonal polycrystalline [Fe(TPP)F], and D = 8.8(2) cm⁻¹, E = 0.1(2) cm⁻¹ and D = 13.4(6) cm⁻¹, E = 0.3(6) cm⁻¹ for monoclinic polycrystalline [Fe(TPP)Br] and [Fe(TPP)I], respectively. Along with our recent report of the ZFS value of D = 6.33(8) cm⁻¹ for tetragonal polycrystalline [Fe(TPP)Cl], these data provide a rare, complete determination of ZFS parameters in a metalloporphyrin halide series. The electronic structure of [Fe(TPP)X] (X = F, Cl, Br, I) has been studied by multireference ab initio methods: the complete active space self-consistent field (CASSCF) and the N-electron valence perturbation theory (NEVPT2) with the aim of exploring the origin of the large and positive zero-field splitting D of the ⁶A₁ ground state. D was calculated from wave functions of the electronic multiplets spanned by the d⁵ configuration of Fe(III) along with spin–orbit coupling accounted for by quasi degenerate perturbation theory. Results reproduce trends of D from inelastic neutron scattering data increasing in the order from F, Cl, Br, to I. A mapping of energy eigenvalues and eigenfunctions of the S = 3/2 excited states on ligand field theory was used to characterize the σ- and π-antibonding effects decreasing from F to I. This is in agreement with similar results deduced from ab initio calculations on CrX₆³⁻ complexes and also with the spectrochemical series showing a decrease of the ligand field in the same directions. A correlation is found between the increase of D and decrease of the π- and σ-antibonding energies e(λ)(X) (λ = σ, π) in the series from X = F to I. Analysis of this correlation using second-order perturbation theory expressions in terms of angular overlap parameters rationalizes the

  20. Structural, dynamical, electronic, and bonding properties of laser-heated silicon: An ab initio molecular-dynamics study

    Silvestrelli, P.-L.; Alavi, A.; Parrinello, M.; Frenkel, D.

    1997-01-01

    The method of ab initio molecular dynamics, based on finite-temperature density-functional theory, is used to simulate laser heating of crystalline silicon. We found that a high concentration of excited electrons dramatically weakens the covalent bonding. As a result the system undergoes a melting t

  1. An ab initio study on the ground and low-lying doublet electronic states of linear C2As

    Graphical abstract: Here we report the results of ab initio calculations of the geometry, adiabatic excitation energies and electric dipole moments of the low-lying doublet electronic states of linear C2As. We also discuss some peculiar details of the vibronic and spin-orbit structure in the X∼2Π and 12Δ states, being the subject of our previous paper. The results for transition moments between the vibronic energy levels within the ground electronic state are presented. Research highlights: → The present study represents a continuation of our theoretical work on spectroscopy of C2As. → We determine the equilibrium geometry and adiabatic transition energies of the low-lying doublet excited states. → We presented the vibronic transitions within the ground X∼2Π electronic state. - Abstract: The present study is closely related to a recently published paper on the X∼2Π-12Δ spectral system of C2As [S. Jerosimic, Lj. Stojanovic, M. Peric, J. Chem. Phys. 133 (2010) 024307]. Here we report the results of ab initio calculations of the geometry, adiabatic excitation energies and electric dipole moments of the low-lying doublet electronic states of linear C2As. We also discuss some peculiar details of the vibronic and spin-orbit structure in the X∼2Π and 12Δ states, being the subject of our previous paper. The results for transition moments between the vibronic energy levels within the ground electronic state are presented. The aim of the present study is to motivate new experimental studies on the title molecule and to facilitate the interpretation of the results of these investigations.

  2. Ab initio study of the lattice stability of β-SiC under intense laser irradiation

    Shen, Yanhong; Gao, Tao, E-mail: gaotao@scu.edu.cn

    2015-10-05

    Highlights: • The band gap of β-SiC vanishes and its metallic character is presented when T{sub e} > 6 eV. • The TA modes of β-SiC are found to be negative T{sub e} = 3.39 eV. • The LO–TO splitting degree of β-SiC at Γ point begin to decline as T{sub e} > 4.5 eV. • The ionic strength of β-SiC is related to laser radiation intensity. - Abstract: We have performed ab initio calculation of electronic properties, lattice-dynamical properties, charge density difference and charge density of β-SiC at different electronic temperatures (T{sub e}) using local density approximation (LDA) pseudopotential method within the density functional perturbation theory (DFPT). The results of electronic density of state display that β-SiC is still semiconductor with band gap of 1.51 eV at T{sub e} = 0 eV. But, beyond a temperature of 6 eV, the band gap of β-SiC vanishes and its metallic character is presented. The calculated phonon frequencies of β-SiC at T{sub e} = 0 eV show a good agreement with the experimental values and other calculations. However, when β-SiC undergoes a sharp increase of its electronic temperature, the phonon frequencies of β-SiC have a significant softening. The transverse acoustic modes of β-SiC are found to be negative T{sub e} = 3.39 eV which lead to the lattice instability. Moreover, the LO–TO splitting degree of β-SiC at Γ point increases at first and then reduces as T{sub e} is raised, the turning point is at T{sub e} = 4.5 eV. By using CUT3D, the results of the charge density difference and charge density of β-SiC indicate that when radiation intensity is only strong enough (e.g. T{sub e} > 4.5 eV), it will make the ionic strength of β-SiC weaken. Otherwise, when radiation intensity is not very high (e.g. 0–4.5 eV), the ionic strength of β-SiC will increase with the rise of T{sub e}.

  3. Ab initio molecular dynamics study of the static, dynamic, and electronic properties of liquid mercury at room temperature.

    Calderín, L; González, L E; González, D J

    2009-05-21

    We report a study on several static, dynamic, and electronic properties of liquid Hg at room temperature. We have performed ab initio molecular dynamics simulations using Kohn-Sham density functional theory combined with a nonlocal ultrasoft pseudopotential. The calculated static structure shows good agreement with the available experimental data. We present results for the single-particle dynamics, and recent experimental data are analyzed. The calculated dynamic structure factors S(q,omega) fairly agree with their experimental counterparts as measured by inelastic x-ray (and neutron) scattering experiments. The dispersion relation exhibits a positive dispersion, which however is not so marked as suggested by the experiment; moreover, its slope at the long-wavelength limit provides a good estimate of the experimental sound velocity. We have also analyzed the dynamical processes behind the S(q,omega) in terms of a model including a relaxation mechanism with both fast and slow characteristic time scales. PMID:19466841

  4. Ab initio molecular dynamics study of the properties of cerium in liquid sodium at 1000 K temperature

    For liquid-sodium-cooled fast nuclear reactor systems, it is crucial to understand the behavior of lanthanides and other potential fission products in liquid sodium or other liquid metal solutions such as liquid cesium-sodium. In this study, we focus on lanthanide behavior in liquid sodium. Using ab initio molecular dynamics, we found that the solubility of cerium in liquid sodium at 1000 K was less than 0.78 at. %, and the diffusion coefficient of cerium in liquid sodium was calculated to be 5.57 × 10−9 m2/s. Furthermore, it was found that cerium in small amounts may significantly alter the heat capacity of the liquid sodium system. Our results are consistent with the experimental results for similar materials under similar conditions

  5. Ab initio molecular dynamics study of the properties of cerium in liquid sodium at 1000 K temperature

    Samin, Adib; Li, Xiang; Zhang, Jinsuo [Nuclear Engineering Program, Department of Mechanical and Aerospace Engineering, The Ohio State University, 201 W 19th Avenue, Columbus, Ohio 43210 (United States); Mariani, R. D. [Idaho National Laboratory, Materials and Fuels Complex, Idaho Falls, Idaho 83415 (United States); Unal, Cetin [Los Alamos National Laboratory, P.O. Box 1663, Los Alamos, New Mexico 87545 (United States)

    2015-12-21

    For liquid-sodium-cooled fast nuclear reactor systems, it is crucial to understand the behavior of lanthanides and other potential fission products in liquid sodium or other liquid metal solutions such as liquid cesium-sodium. In this study, we focus on lanthanide behavior in liquid sodium. Using ab initio molecular dynamics, we found that the solubility of cerium in liquid sodium at 1000 K was less than 0.78 at. %, and the diffusion coefficient of cerium in liquid sodium was calculated to be 5.57 × 10{sup −9} m{sup 2}/s. Furthermore, it was found that cerium in small amounts may significantly alter the heat capacity of the liquid sodium system. Our results are consistent with the experimental results for similar materials under similar conditions.

  6. Ab initio molecular dynamics study of the properties of cerium in liquid sodium at 1000 K temperature

    Samin, Adib; Li, Xiang; Zhang, Jinsuo; Mariani, R. D.; Unal, Cetin

    2015-12-01

    For liquid-sodium-cooled fast nuclear reactor systems, it is crucial to understand the behavior of lanthanides and other potential fission products in liquid sodium or other liquid metal solutions such as liquid cesium-sodium. In this study, we focus on lanthanide behavior in liquid sodium. Using ab initio molecular dynamics, we found that the solubility of cerium in liquid sodium at 1000 K was less than 0.78 at. %, and the diffusion coefficient of cerium in liquid sodium was calculated to be 5.57 × 10-9 m2/s. Furthermore, it was found that cerium in small amounts may significantly alter the heat capacity of the liquid sodium system. Our results are consistent with the experimental results for similar materials under similar conditions.

  7. Fingerprints of antiaromaticity in the negative ion (Li_3Al_4)^- via an ab initio quantum-chemical study

    Grassi, A; Angilella, G G N; March, N H; Pucci, R

    2012-01-01

    Fingerprints of antiaromaticity in the negative ion (Li_3Al_4)^-, this species being realizable via a laser vaporization technique, are revealed by means of an ab initio quantum-chemical investigation. First, the ground-state equilibrium geometry of this ion is predicted. Also, the characteristics of the HOMO are studied, both for the square and the rectangular Al_4 geometry in two low-lying isomers of the negative ion. There is no particular sensitivity to the change in geometry of the Al_4 configuration. Therefore, we have calculated theoretically chemical shifts, which contain remarkable fingerprints of antiaromaticity. As to future directions, some comments are added in relation to the Shannon entropy.

  8. Ab initio mass tensor molecular dynamics

    Tsuchida, Eiji

    2010-01-01

    Mass tensor molecular dynamics was first introduced by Bennett [J. Comput. Phys. 19, 267 (1975)] for efficient sampling of phase space through the use of generalized atomic masses. Here, we show how to apply this method to ab initio molecular dynamics simulations with minimal computational overhead. Test calculations on liquid water show a threefold reduction in computational effort without making the fixed geometry approximation. We also present a simple recipe for estimating the optimal ato...

  9. Ab initio vel ex eventu. II

    Thiessen, P. A.; Treder, H.-J.

    Jedes initium wird durch experimenta crucis zum eventus. Jedes theoretisch interpretierbare ex-eventu-Resultat führt auf ein neues Initium. Gerade dies ist die gemeinsame Aussage von Atomistik, Quantenmechanik und Relativitätstheorie.Translated AbstractAb initio vel ex eventu. IIEvery initium becomes an eventus by experimenta crucis. Every theoretically interpretable ex-eventu result leads to a new initium. Right this is the joint assertion of atomism, quantum mechanics, and relativity.

  10. Highly scalable Ab initio genomic motif identification

    Marchand, Benoît

    2011-01-01

    We present results of scaling an ab initio motif family identification system, Dragon Motif Finder (DMF), to 65,536 processor cores of IBM Blue Gene/P. DMF seeks groups of mutually similar polynucleotide patterns within a set of genomic sequences and builds various motif families from them. Such information is of relevance to many problems in life sciences. Prior attempts to scale such ab initio motif-finding algorithms achieved limited success. We solve the scalability issues using a combination of mixed-mode MPI-OpenMP parallel programming, master-slave work assignment, multi-level workload distribution, multi-level MPI collectives, and serial optimizations. While the scalability of our algorithm was excellent (94% parallel efficiency on 65,536 cores relative to 256 cores on a modest-size problem), the final speedup with respect to the original serial code exceeded 250,000 when serial optimizations are included. This enabled us to carry out many large-scale ab initio motiffinding simulations in a few hours while the original serial code would have needed decades of execution time. Copyright 2011 ACM.

  11. Ab initio and DFT study of the geometric structures and static dipole (hyper)polarizabilities of aromatic anions.

    Castellano, O; Bermúdez, Y; Giffard, M; Mabon, G; Cubillan, N; Sylla, M; Nguyen-Phu, X; Hinchliffe, A; Soscún, H

    2005-11-17

    The geometries and the static dipole (hyper)polarizabilities (alpha, beta, gamma) of a series of aromatic anions were investigated at the ab initio (HF, MP2, and MP4) and density functional theory DFT (B3LYP) levels of theory. The anions chosen for the present study are the benzenethiolate (Ph-S-), benzenecarboxylate (Ph-CO2-), benzenesulfinate (Ph-SO2-), benzenesulfonate (Ph-SO3-), and 1,3-benzenedicarboxylate (1,3-Ph-(CO2)2(2-)). For benzenethiolate anion, additional alpha, beta, and gamma calculations were performed at the coupled cluster CCSD level with MP2 optimized geometries. The standard diffuse and polarized 6-31+G(d,p) basis set was employed in conjunction to the ab initio and DFT methods. Additional HF calculations were performed with the 6-311++G(3d,3p) basis set for all the anions. The correlated electric properties were evaluated numerically within the formalism of finite field. The optimized geometries were analyzed in terms of the few reports about the phenolate and sulfonate ions. The results show that electron correlation effects on the polarizabilities are very important in all the anion series. Was found that Ph-SO2- is highly polarizable in terms of alpha and beta, and the Ph-S- is the highest second hyperpolarizable in the series. The results of alpha were rationalized in terms of the analysis of the polarization of charge based in Mulliken atomic population and the structural features of the optimized geometries of anions, whereas the large differences in the beta and gamma values in the series were respectively interpreted in terms of the bond length alternation BLA and the separation of charge in the aromatic ring by effects of the substitution. These results allowed us to suggest the benzenesulfinate and benzenethiolate anions as promising candidates that should be incorporated in ionic materials for second and third-order nonlinear optical devices. PMID:16833334

  12. An experimental and ab initio study of the electronic spectrum of the jet-cooled F2BO free radical

    We have studied the B~2A1–X~2B2 laser-induced fluorescence (LIF) spectrum of the jet-cooled F2BO radical for the first time. The transition consists of a strong 000 band at 446.5 nm and eight weak sequence bands to shorter wavelengths. Single vibronic level emission spectra obtained by laser excitation of individual levels of the B~ state exhibit two electronic transitions: a very weak, sparse B~–X~ band system in the 450–500 nm region and a stronger, more extensive set of B~2A1–A~2B1 bands in the 580–650 nm region. We have also performed a series of high level ab initio calculations to predict the electronic energies, molecular structures, vibrational frequencies, and rotational and spin-rotation constants in the X~2B2, A~2B1 and B~2A1 electronic states as an aid to the analysis of the experimental data. The theoretical results have been used as input for simulations of the rotationally resolved B~2A1–X~2B2 000 LIF band and Franck-Condon profiles of the LIF and single vibronic level emission spectra. The agreement between the simulations obtained with purely ab initio parameters and the experimental spectra validates the geometries calculated for the ground and excited states and the conclusion that the radical has C2v symmetry in the X~, A~, and B~ states. The spectra provide considerable new information about the vibrational energy levels of the X~ and A~ states, but very little for the B~ state, due to the very restrictive Franck-Condon factors in the LIF spectra

  13. Ab initio studies of O2-(H2O)n and O3-(H2O)n anionic molecular clusters, n≤12

    Bork, Nicolai Christian; Kurtén, T.; Enghoff, Martin Andreas Bødker;

    2011-01-01

    An ab initio study of gaseous clusters of O2− and O2− with water is presented. Based on thorough scans of configurational space, we determine the thermodynamics of cluster growth. The results are in good agreement with benchmark computational methods and existing experimental data. We find...

  14. Ab initio studies of O-2(-) (H2O)(n) and O-3(-) (H2O)(n) anionic molecular clusters, n

    Bork, Nicolai Christian; Kurten, T.; Enghoff, Martin Andreas Bødker;

    2011-01-01

    An ab initio study of gaseous clusters of O-2(-) and O-3(-) with water is presented. Based on thorough scans of configurational space, we determine the thermodynamics of cluster growth. The results are in good agreement with benchmark computational methods and existing experimental data. We find...

  15. On the room-temperature phase diagram of high pressure hydrogen: an ab initio molecular dynamics perspective and a diffusion Monte Carlo study.

    Chen, Ji; Ren, Xinguo; Li, Xin-Zheng; Alfè, Dario; Wang, Enge

    2014-07-14

    The finite-temperature phase diagram of hydrogen in the region of phase IV and its neighborhood was studied using the ab initio molecular dynamics (MD) and the ab initio path-integral molecular dynamics (PIMD). The electronic structures were analyzed using the density-functional theory (DFT), the random-phase approximation, and the diffusion Monte Carlo (DMC) methods. Taking the state-of-the-art DMC results as benchmark, comparisons of the energy differences between structures generated from the MD and PIMD simulations, with molecular and dissociated hydrogens, respectively, in the weak molecular layers of phase IV, indicate that standard functionals in DFT tend to underestimate the dissociation barrier of the weak molecular layers in this mixed phase. Because of this underestimation, inclusion of the quantum nuclear effects (QNEs) in PIMD using electronic structures generated with these functionals leads to artificially dissociated hydrogen layers in phase IV and an error compensation between the neglect of QNEs and the deficiencies of these functionals in standard ab initio MD simulations exists. This analysis partly rationalizes why earlier ab initio MD simulations complement so well the experimental observations. The temperature and pressure dependencies for the stability of phase IV were also studied in the end and compared with earlier results. PMID:25028021

  16. On the room-temperature phase diagram of high pressure hydrogen: An ab initio molecular dynamics perspective and a diffusion Monte Carlo study

    The finite-temperature phase diagram of hydrogen in the region of phase IV and its neighborhood was studied using the ab initio molecular dynamics (MD) and the ab initio path-integral molecular dynamics (PIMD). The electronic structures were analyzed using the density-functional theory (DFT), the random-phase approximation, and the diffusion Monte Carlo (DMC) methods. Taking the state-of-the-art DMC results as benchmark, comparisons of the energy differences between structures generated from the MD and PIMD simulations, with molecular and dissociated hydrogens, respectively, in the weak molecular layers of phase IV, indicate that standard functionals in DFT tend to underestimate the dissociation barrier of the weak molecular layers in this mixed phase. Because of this underestimation, inclusion of the quantum nuclear effects (QNEs) in PIMD using electronic structures generated with these functionals leads to artificially dissociated hydrogen layers in phase IV and an error compensation between the neglect of QNEs and the deficiencies of these functionals in standard ab initio MD simulations exists. This analysis partly rationalizes why earlier ab initio MD simulations complement so well the experimental observations. The temperature and pressure dependencies for the stability of phase IV were also studied in the end and compared with earlier results

  17. Ab initio and experimental study of the interaction of nitrous oxide with the isolated hydroxyl of silica

    Garrone, E.; Ugliengo, P.; Ghiotti, G.; Borello, E.; Saunders, V. R.

    1993-08-01

    High-quality ab initio computational results for the interaction of N 2O with the silanol molecule H 3SiOH, taken as representative of the isolated hydroxyl at the silica surface, are compared with IR data concerning the real system measured at 195 K. Besides a non-specific interaction, features typical of hydrogen bonding are observed: two types of complexes are formed simultaneously, differing in their 2 v1 mode, probably one nitrogen-bonded, the other oxygen-bonded. Ab initio results are in excellent agreement as far as binding energies are concerned: the two complexes show comparable stabilities, once electron correlation has been taken into account; however, SCF-calculated IR features do not allow a definite vibrational assignment of the two complexes.

  18. Ab initio no core full configuration approach for light nuclei

    Kim, Youngman; Shin, Ik Jae; Maris, Pieter; Vary, James P.; Forssén, Christian; Rotureau, Jimmy

    2015-10-01

    Comprehensive understanding of the structure and reactions of light nuclei poses theoretical and computational challenges. Still, a number of ab initio approaches have been developed to calculate the properties of atomic nuclei using fundamental interactions among nucleons. Among them, we work with the ab initio no core full configuration (NCFC) method and ab initio no core Gamow Shell Model (GSM). We first review these approaches and present some recent results.

  19. Ab Initio Calculations of Oxosulfatovanadates

    Frøberg, Torben; Johansen, Helge

    1996-01-01

    Restricted Hartree-Fock and multi-configurational self-consistent-field calculations together with secondorder perturbation theory have been used to study the geometry, the electron density, and the electronicspectrum of (VO2SO4)-. A bidentate sulphate attachment to vanadium was found to be stable...... with anO-V-O angle of 72.5 degrees . The calculated spectrum shows bands in reasonable agreement with anexperimental spectrum which has been attributed to (VO2SO4)-. The geometry and the electron density fortwo binuclear vanadium complexes proposed as intermediates in the vanadium catalyzed SO2...

  20. Ab initio calculations of material strength

    Šob, Mojmír; Friák, Martin; Vitek, V.

    Tokyo : The Japan Society of Mechanical Engineers, 2003, s. 467-475. [International Symposium on Micro-Mechanical Engineering - Heat Transfer, Fluid Dynamics, Reliability and Mechanotronics.. Tsuchiura and Tsukuba (JP), 01.12.2003-03.12.2003] R&D Projects: GA AV ČR IAA1041302; GA ČR GA202/03/1351; GA MŠk OC 523.90 Institutional research plan: CEZ:AV0Z2041904 Keywords : ab initio calculations * electronic structure * theoretical tensile strength Subject RIV: BM - Solid Matter Physics ; Magnetism

  1. Ab Initio Molecular Dynamics: A Virtual Laboratory

    Hobbi Mobarhan, Milad

    2014-01-01

    In this thesis, we perform ab initio molecular dynamics (MD) simulations at the Hartree-Fock level, where the forces are computed on-the-fly using the Born-Oppenheimer approximation. The theory behind the Hartree-Fock method is discussed in detail and an implementation of this method based on Gaussian basis functions is explained. We also demonstrate how to calculate the analytic energy derivatives needed for obtaining the forces acting on the nuclei. Hartree-Fock calculations on the ground s...

  2. Germacrene D Cyclization: An Ab Initio Investigation

    William N. Setzer

    2008-01-01

    Full Text Available Essential oils that contain large concentrations of germacrene D are typically accompanied by cadinane sesquiterpenoids. The acid-catalyzed cyclization of germacrene D to give cadinane and selinane sesquiterpenes has been computationally investigated using both density functional (B3LYP/6-31G* and post Hartree-Fock (MP2/6-31G** ab initio methods. The calculated energies are in general agreement with experimentally observed product distributions, both from acid-catalyzed cyclizations as well as distribution of the compounds in essential oils.

  3. Ab initio non-relativistic spin dynamics

    Ding, Feizhi; Goings, Joshua J.; Li, Xiaosong, E-mail: xsli@uw.edu [Department of Chemistry, University of Washington, Seattle, Washington 98195 (United States); Frisch, Michael J. [Gaussian, Inc., 340 Quinnipiac St, Bldg 40, Wallingford, Connecticut 06492 (United States)

    2014-12-07

    Many magnetic materials do not conform to the (anti-)ferromagnetic paradigm where all electronic spins are aligned to a global magnetization axis. Unfortunately, most electronic structure methods cannot describe such materials with noncollinear electron spin on account of formally requiring spin alignment. To overcome this limitation, it is necessary to generalize electronic structure methods and allow each electron spin to rotate freely. Here, we report the development of an ab initio time-dependent non-relativistic two-component spinor (TDN2C), which is a generalization of the time-dependent Hartree-Fock equations. Propagating the TDN2C equations in the time domain allows for the first-principles description of spin dynamics. A numerical tool based on the Hirshfeld partitioning scheme is developed to analyze the time-dependent spin magnetization. In this work, we also introduce the coupling between electron spin and a homogenous magnetic field into the TDN2C framework to simulate the response of the electronic spin degrees of freedom to an external magnetic field. This is illustrated for several model systems, including the spin-frustrated Li{sub 3} molecule. Exact agreement is found between numerical and analytic results for Larmor precession of hydrogen and lithium atoms. The TDN2C method paves the way for the ab initio description of molecular spin transport and spintronics in the time domain.

  4. Coupling of temperature with pressure induced initial decomposition mechanisms of two molecular crystals: An ab initio molecular dynamics study

    QIONG WU; DONG XIANG; GUOLIN XIONG; WEIHUA ZHU; HEMING XIAO

    2016-05-01

    Ab initio molecular dynamics simulations were performed to study the initiation of decompositionand formation of first products of two molecular crystals pentaerythritol tetranitrate (PETN) and 5-nitro-2,4-dihydro-1,2,4-triazole-3-one (NTO) under thermal decomposition temperature (475 K for PETN and 531 Kfor NTO) coupled with different pressures (1-5 GPa). The pressure effects on the initial decomposition stepsand initially generated products on PETN and NTO were very different. PETN was triggered by C-H... O intermolecular hydrogen transfer. The initial decomposition mechanism was independent of the pressure. ForNTO, two different initial decomposition mechanisms were found. At 1, 2, and 3 GPa, it was triggered by NH....O intermolecular hydrogen transfer, while at 4 and 5 GPa, it was triggered by N-H.....N intermolecularhydrogen transfer. This indicates that the initial decomposition mechanism was dependent on the pressure.Our study may provide new insights into initial mechanisms and decomposition reactions of molecular crystalexplosives under thermal decomposition temperature coupled with different pressures with details at atomiclevel.

  5. Ab initio studies on [bmim][PF6]–CO2 mixture and CO2 clusters

    B L Bhargava; M Saharay; S Balasubramanian

    2008-06-01

    Ab initio molecular dynamics studies have been carried out on the room temperature ionic liquid, 1,n-butyl,3-methylimidazolium hexafluorophosphate ([bmim][PF6]) and supercritical carbon dioxide mixture at room temperature and experimental density. Partial radial distribution functions (RDF) for different sites have been computed to see the organization of CO2 molecules around the ionic liquid. Several partial RDFs around the carbon atom of CO2 molecule are compared to find out that the CO2 has specific interaction with a carbon atom present in the imidazolium ring. The CO2 is also found to be very well organized around the terminal carbon atom of the butyl chain. The partial RDFs for the oxygen atoms around oxygen and carbon atoms of the CO2 suggests that there is very good organization of CO2 molecules around themselves even in the [bmim][PF6] – CO2 mixture. The instantaneous quadrupole moment tensor has been calculated for the anion and the cation. The ensemble average of diagonal components of quadrupole moment tensor of the cation have finite values, whereas the off-diagonal components of the cation and both the diagonal and off-diagonal components of the anion have the value of zero with a large standard deviation. The CPMD studies performed on CO2 clusters reveals the greater tendency of the clusters with more CO2 units, to deviate from the linear geometry.

  6. Vibrational lifetimes of hydrogen on lead films: An ab initio molecular dynamics with electronic friction (AIMDEF) study

    Saalfrank, Peter [Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Strasse 24-25, D-14476 Potsdam (Germany); Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián (Spain); Juaristi, J. I. [Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián (Spain); Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián (Spain); Departamento de Física de Materiales, Facultad de Químicas UPV/EHU, Apartado 1072, 20018 Donostia-San Sebastián (Spain); Alducin, M.; Muiño, R. Díez [Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián (Spain); Centro de Física de Materiales CFM/MPC (CSIC-UPV/EHU), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián (Spain); Blanco-Rey, M. [Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián (Spain); Departamento de Física de Materiales, Facultad de Químicas UPV/EHU, Apartado 1072, 20018 Donostia-San Sebastián (Spain)

    2014-12-21

    Using density functional theory and Ab Initio Molecular Dynamics with Electronic Friction (AIMDEF), we study the adsorption and dissipative vibrational dynamics of hydrogen atoms chemisorbed on free-standing lead films of increasing thickness. Lead films are known for their oscillatory behaviour of certain properties with increasing thickness, e.g., energy and electron spillout change in discontinuous manner, due to quantum size effects [G. Materzanini, P. Saalfrank, and P. J. D. Lindan, Phys. Rev. B 63, 235405 (2001)]. Here, we demonstrate that oscillatory features arise also for hydrogen when chemisorbed on lead films. Besides stationary properties of the adsorbate, we concentrate on finite vibrational lifetimes of H-surface vibrations. As shown by AIMDEF, the damping via vibration-electron hole pair coupling dominates clearly over the vibration-phonon channel, in particular for high-frequency modes. Vibrational relaxation times are a characteristic function of layer thickness due to the oscillating behaviour of the embedding surface electronic density. Implications derived from AIMDEF for frictional many-atom dynamics, and physisorbed species will also be given.

  7. Ab initio study of anisotropic mechanical properties of LiCoO2 during lithium intercalation and deintercalation process

    The mechanical properties of LixCoO2 under various Li concentrations and associated anisotropy have been systematically studied using the first principles method. During the lithium intercalation process, the Young's modulus, bulk modulus, shear modulus, and ultimate strength increase with increasing lithium concentration. Strong anisotropy of mechanical properties between a-axis and c-axis in LixCoO2 is identified at low lithium concentrations, and the anisotropy decreases with increasing lithium concentration. The observed lithium concentration dependence and anisotropy are explained by analyzing the charge transfer using Bader charge analysis, bond order analysis, and bond strength by investigating partial density of states and charge density difference. With the decrease of Li concentration, the charge depletion in the bonding regions increases, indicating a weaker Co-O bond strength. Additionally, the Young's modulus, bulk modulus, shear modulus, and toughness are obtained by simulating ab initio tensile tests. From the simulated stress-strain curves, LixCoO2 shows the highest toughness, which is in contraction with Pugh criterion prediction based on elastic properties only

  8. 42214 layered Fe-based superconductors: An ab initio study of their structural, magnetic, and electronic properties

    Bucci, F.; Sanna, A.; Continenza, A.; Katrych, S.; Karpinski, J.; Gross, E. K. U.; Profeta, G.

    2016-01-01

    As a follow-up to the discovery of a new family of Fe-based superconductors, namely, the RE4Fe2As2Te1 -xO4 (42214) (RE = Pr, Sm, and Gd), we present a detailed ab initio study of these compounds highlighting the role of rare-earth (RE) atoms, external pressure, and Te content on their physical properties. Modifications of the structural, magnetic, and electronic properties of the pure (e.g., x =0.0 ) 42214 compounds and their possible correlations with the observed superconducting properties are calculated and discussed. The careful analysis of the results obtained shows that (i) changing the RE atoms allows one to tune the internal pressure acting on the As height with respect to the Fe planes; (ii) similarly to other Fe pnictides, the 42214 pure compounds show an antiferromagnetic-stripe magnetic ground state phase joined by an orthorhombic distortion (not experimentally found yet); (iii) smaller RE atoms increase the magnetic instability of the compounds possibly favoring the onset of the superconducting state; (iv) external pressure induces the vanishing of the magnetic order with a transition to the tetragonal phase and can be a possible experimental route towards higher superconducting critical temperature (Tc) ; and (v) Te vacancies act on the structural parameters, changing the As height and affecting the stability of the magnetic phase.

  9. Ab initio study of the one-dimensional H-bonded ferroelectric CsH2PO4

    Lasave, J.; Abufager, P.; Koval, S.

    2016-04-01

    We studied the microscopic mechanism of the paraelectric-ferroelectric (PE-FE) phase transition of CsH2PO4 (CDP) by means of first-principles electronic structure calculations. The calculated structural parameters in the PE and FE phases as well as the total spontaneous polarization Ps obtained with the Berry phase formalism for CDP are in good agreement with experiments. The main contribution to Ps originates from a large y x component of the calculated Born effective-charge tensor for the disordered protons. Moreover, this component is ≈2.7 times larger than the z x component of the proton effective-charge tensor relevant to the polarization in the H-bonded FE KH2PO4 (KDP). This is the main feature that compensates the different number of protons per formula unit involved in the phase transitions for CDP and KDP to give close values for their measured Ps. Correlations among protons and heavy atoms along chains in the b direction lead to larger instabilities for the global and local FE distortions in CDP and its deuterated counterpart DCDP. We conclude that the tunneling particle is a dressed proton (deuteron) in agreement with experiments and with recent ab initio results for KDP.

  10. Ab Initio Study of Water Adsorption and Reactivity on the (211) Surface of Anatase TiO2

    Xu, Jing; Xu, Li-Fang; Li, Zhen-Zhen; Wang, Jian-Tao; Lin, Zhe-Shuai; Liu, Kai; Cao, Yong-Ge; Selloni, Annabella

    2016-06-01

    The reactivity of the anatase TiO2 (211) surface is systematically studied by ab initio calculations of the surface energy and water-adsorption energy. We find that anatase (211) has a high surface energy of 0.97 J /m2 , close to that of the (001) surface, and the unsaturated fourfold-coordinated Ti4 atom is more reactive than the unsaturated fivefold-coordinated Ti5 atom. Accordingly, for water adsorption on the (211) surface, a dissociative form is favored on Ti4 sites, with a large adsorption energy Δ HH ,OH˜1.28 eV , while a nondissociative molecular form is favored on Ti5 sites, with a smaller adsorption energy Δ HH2O˜0.78 eV . Such distinct surface properties lead to a mixed dissociative and molecular adsorption configuration when the coverage is increased from 1 /3 to 1 monolayer. These results suggest that, similar to the (001) surface, the anatase (211) surface exhibits high reactivity and may be useful in catalytic and photocatalytic applications as well.

  11. A comparative ab initio study of superconductivity in the body centered tetragonal YC2 and LaC2

    Ab initio studies of the electronic band structure and phonon dispersion relations, using the planewave pseudopotential method and the density functional theory, have been made for the superconducting materials YC2 and LaC2. Differences in the phonon spectrum and density of states both in the acoustical and optical ranges between these materials are investigated and discussed. By integrating the Eliashberg spectral function α2F(ω), the average electron-phonon coupling parameter is found to be λ = 0.55 for YC2 and 0.54 for LaC2, indicating these to be weak-coupling BCS superconductors. It is established that about 60% of λ is contributed by acoustic phonons in both materials. Using a reasonable value of μ* = 0.13 for the effective Coulomb repulsion parameter, the superconducting critical temperature Tc is found to be 3.81 K for YC2 and 2.44 K for LaC2, in good agreement with values reported from experimental measurements

  12. Vibrational lifetimes of hydrogen on lead films: An ab initio molecular dynamics with electronic friction (AIMDEF) study

    Using density functional theory and Ab Initio Molecular Dynamics with Electronic Friction (AIMDEF), we study the adsorption and dissipative vibrational dynamics of hydrogen atoms chemisorbed on free-standing lead films of increasing thickness. Lead films are known for their oscillatory behaviour of certain properties with increasing thickness, e.g., energy and electron spillout change in discontinuous manner, due to quantum size effects [G. Materzanini, P. Saalfrank, and P. J. D. Lindan, Phys. Rev. B 63, 235405 (2001)]. Here, we demonstrate that oscillatory features arise also for hydrogen when chemisorbed on lead films. Besides stationary properties of the adsorbate, we concentrate on finite vibrational lifetimes of H-surface vibrations. As shown by AIMDEF, the damping via vibration-electron hole pair coupling dominates clearly over the vibration-phonon channel, in particular for high-frequency modes. Vibrational relaxation times are a characteristic function of layer thickness due to the oscillating behaviour of the embedding surface electronic density. Implications derived from AIMDEF for frictional many-atom dynamics, and physisorbed species will also be given

  13. Ab initio study of Kr in hcp Ti: Diffusion, formation and stability of small Kr-vacancy clusters

    Ab initio electronic structure calculations have been performed to study the formation and migration of Kr impurities, and the stability of small Kr-vacancy clusters for clusters with up to four vacancies and four Kr atoms, in hcp Ti. Both the substitutional and the interstitial configurations of Kr are found to be stable. The octahedral configuration is however found to be more stable than the tetrahedral. Interstitial Kr atoms are shown to have attractive interactions and a low migration barrier, suggesting that, at low temperature, Kr bubble formation is possible, even in the absence of vacancies. We also find vacancy clusters to be stable. The binding energies of an interstitial Kr atom and a vacancy to a Kr-vacancy cluster are obtained from the calculated formation energies of the clusters. The stability of small-vacancy clusters is found to be dependent on Kr-vacancy ratio. The trends of the calculated binding energies are discussed in terms of providing further insights on the behaviour of Kr in implanted Ti.

  14. Ab initio quantum mechanical study for the photolysis and unimolecular decomposition reactions in the atmosphere of CF3OH

    The electronic transitions from the ground state to low-lying excited states of CF3OH have been investigated using high level ab initio quantum mechanical techniques. Also the possible photo dissociation procedures of CF3OH have been considered. The highest level employed in this study is TZP CCSD(T) level of theory. The possible four low-lying excited states can result by the excitation of the lone pair electron (n) in oxygen to σ molecular orbital in C-O or O-H bond. The vertical transition (n→σ) energy is predicted to be 220.5 kcal/mol (130 nm) at TZ2P CISD level to theory. The bond dissociation energies of CF3OH to CF3O + H and CF3+OH have been predicted to be 119.5 kcal/mol 114.1 kcal/mol, respectively, at TZP CCSD level of theory. In addition, the transition state for the unimolecular decomposition of CF3OH into CF2O + HF has been examined. The activation energy and energy separation for this decomposition have been computed to be 43.6 kcal/mol and 5.0 kcal/mol including zero-point vibrational energy corrections at TZP CCSD(T) level of theory

  15. Ab initio and atomistic study of generalized stacking fault energies in Mg and Mg–Y alloys

    Magnesium–yttrium alloys show significantly improved room temperature ductility when compared with pure Mg. We study this interesting phenomenon theoretically at the atomic scale employing quantum-mechanical (so-called ab initio) and atomistic modeling methods. Specifically, we have calculated generalized stacking fault energies for five slip systems in both elemental magnesium (Mg) and Mg–Y alloys using (i) density functional theory and (ii) a set of embedded-atom-method (EAM) potentials. These calculations predict that the addition of yttrium results in a reduction in the unstable stacking fault energy of basal slip systems. Specifically in the case of an I2 stacking fault, the predicted reduction of the stacking fault energy due to Y atoms was verified by experimental measurements. We find a similar reduction for the stable stacking fault energy of the {11 2-bar 2} non-basal slip system. On the other hand, other energies along this particular γ-surface profile increase with the addition of Y. In parallel to our quantum-mechanical calculations, we have also developed a new EAM Mg–Y potential and thoroughly tested its performance. The comparison of quantum-mechanical and atomistic results indicates that the new potential is suitable for future large-scale atomistic simulations. (paper)

  16. Kopplung von Dichtefunktional- und ab-initio-Methoden

    Goll, Erich

    2008-01-01

    Im Rahmen der Doktorarbeit wurde untersucht, inwieweit die Kopplung von Dichtefunktionalmethoden und ab-initio-Korrelationsmethoden der Quantenchemie eine Verbesserung bezüglich beider Grenzmethoden erbringt. Die Kopplung erfolgt durch eine Aufspaltung des interelektronischen Hamiltonoperators (abstoßende Coulombwechselwirkung). Die kurzreichweitige Wechselwirkung wird mit Dichtefunktionaltheorie behandelt, die langreichweitige mit Hilfe von ab-initio-Methoden. Diese Aufteilung soll dazu dien...

  17. Exploring the impact of semicore level electronic relaxation on polaron dynamics: An adiabatic ab initio study of FePO4

    Wang, Zi; Bevan, Kirk H.

    2016-01-01

    In the present work, we study the effects of the electronic relaxation of semicore levels on polaron activation energies and dynamics. Within the framework of adiabatic ab initio theory, we utilize both static transition state theory and molecular dynamics methods for an in-depth study of polaronic hopping in delithiated LiFePO4 (FePO4). Our results show that electronic relaxation of semicore states is significant in FePO4, resulting in a lower activation barrier and kinetics that is one to two orders faster compared to the result of calculations that do not incorporate semicore states. In general, the results suggest that the relaxation of states far below the Fermi energy could dramatically impact the ab initio polaronic barrier estimates for many transition metal oxides and phosphates.

  18. Theoretical studies of H2--H2 collisions. IV. Ab initio calculations of anisotropic transport phenomena in para-hydrogen gas

    The temperature dependence of the effective Waldmann--Snider cross sections determining the Senftleben--Beenakker effects of viscosity and heat conductivity has been studied for pH2 gas between 10 and 200 K. From ab initio nonspherical potentials of H2--H2, scattering matrices have been determined in close-coupling calculations. From these, the elements of the scattering amplitude matrix have been obtained and used as input quantities for the evaluation of the various Waldmann--Snider collision integrals. The results of these first ab initio numerical calculations of anisotropic transport coefficients show excellent agreement of calculated and measured effective cross sections, especially for the most recent improved version of the interaction potential. In addition, it has been shown that the polarization production cross sections are quite sensitive to the potential anisotropy

  19. AB INITIO STUDY OF CHEMICAL ACTIVATION AND HYDROGENATION OF WHITE PHOSPHORUS IN REACTION WITH RHODIUM TRIHYDRIDE COMPLEX

    Iolanta I. Balan; Natalia N. Gorinchoy

    2011-01-01

    The four-stage mechanism of reaction of the rhodium trihydride complex [(triphos)RhH3] (triphos=1,1,1-tris(diphenylphosphanylmethyl)ethane) with the white phosphorus molecule resulting in the phosphane and the cyclo-P3 complex [(triphos)M(η3-P3] is analyzed on the basis of ab initio calculations of reactants, products, and intermediate complexes of reaction. It is shown that generation of the transient complex [(triphos)RhH(η1:η1-P4)] followed by intramolecular hydrogen atom migration from t...

  20. Ab initio calculations as a quantitative tool in the inelastic neutron scattering study of a single-molecule magnet analogue.

    Vonci, Michele; Giansiracusa, Marcus J; Gable, Robert W; Van den Heuvel, Willem; Latham, Kay; Moubaraki, Boujemaa; Murray, Keith S; Yu, Dehong; Mole, Richard A; Soncini, Alessandro; Boskovic, Colette

    2016-02-01

    Ab initio calculations carried out on the Tb analogue of the single-molecule magnet family Na9[Ln(W5O18)2] (Ln = Nd, Gd, Ho and Er) have allowed interpretation of the inelastic neutron scattering spectra. The combined experimental and theoretical approach sheds new light on the sensitivity of the electronic structure of the Tb(III) ground and excited states to small structural distortions from axial symmetry, thus revealing the subtle relationship between molecular geometry and magnetic properties of the two isostructural species that comprise the sample. PMID:26690503

  1. Stability of nucleic acid base pairs in organic solvents: Molecular dynamics, molecular dynamics/quenching, and correlated ab initio study

    Zendlová, Lucie; Hobza, Pavel; Kabeláč, Martin

    2007-01-01

    Roč. 111, č. 10 (2007), s. 2591-2609. ISSN 1520-6106 R&D Projects: GA MŠk LC512; GA ČR GA203/05/0009; GA ČR(CZ) GD203/05/H001; GA AV ČR KJB400550518 Institutional research plan: CEZ:AV0Z40550506 Keywords : ab initio calculation * organic solvent * nucleobases * H-bonding * stacking Subject RIV: CF - Physical ; Theoretical Chemistry Impact factor: 4.086, year: 2007

  2. Orbital free ab initio molecular dynamics simulation study of some static and dynamic properties of liquid noble metals

    G.M. Bhuiyan

    2012-10-01

    Full Text Available Several static and dynamic properties of liquid Cu, Ag and Au at thermodynamic states near their respective melting points, have been evaluated by means of the orbital free ab-initio molecular dynamics simulation method. The calculated static structure shows good agreement with the available X-ray and neutron diffraction data. As for the dynamic properties, the calculated dynamic structure factors point to the existence of collective density excitations along with a positive dispersion for l-Cu and l-Ag. Several transport coefficients have been obtained which show a reasonable agreement with the available experimental data.

  3. Solvation structure of glucosamine in aqueous solution as studied by Monte Carlo simulation using ab initio fitted potential

    Siraleartmukul, Krisana; Siriwong, Khatcharin; Remsungnen, Tawun; Muangsin, Nongnuj; Udomkichdecha, Werasak; Hannongbua, Supot

    2004-09-01

    The solvation structure of glucosamine in aqueous solution was investigated using Monte Carlo simulation at 298 K. The MCY rigid water model and ab initio glucosamine-water fitted potential were applied. The first hydration shell appears at 4.6 Å from the center of glucosamine with a coordination number of seven water molecules where one water lies in the ligand's plane while two and four of them are about 2-4 Å above and below the plane, respectively. Furthermore, the mobility distribution and orientation of the water molecules around the ligand have been intensively investigated and reported.

  4. Magnetism in zigzag and armchair CuO nanotubes: Ab-initio study

    Paudel, Samir; Dandeliya, Sushmita; Chaurasiya, Rajneesh; Srivastava, Anurag; Kaphle, Gopi Chandra

    2016-05-01

    The structural stability, electronic band structure and magnetic properties of zigzag (4≤n≤12) and armchair (3≤n≤8) copper oxide nanotubes have been analyzed by employing a standard Density Functional Theory based ab-intio approach using spin polarized generalized gradient approximation with revised Perdew Burke Ernzerhoff type parameterization. The binding energy of both zigzag as well as armchair CuO nanotubes increases with increasing diameter. Whereas, calculated total magnetic moment decreases with increasing diameter. Among the considered chiralities, (5, 5) armchair (diameter 9.08 Å) and (9, 0) zigzag (diameter 9.47 Å) CuO nanotubes show highest degree of spin polarization, however total magnetic moment is found to be highest for (4, 0) zigzag and (3, 3) armchair CuO nanotubes. The computed electronic properties of considered CuO Nanotube, confirms the metallic nature of these nanotubes.

  5. Ab initio simulation of helium inside carbon nanotubes

    In present work we consider the complex behaviour of quantum liquids like liquid He-4 inside carbon nanotubes. Interactions between helium atoms and carbon atoms of the short-length atomistic model and model with periodical boundary conditions of carbon nanotube were studied via ab initio quantum simulations. Effects of geometrical confinement of the tube on the He behaviour inside CNT (13,0) have been explored. Nanotubes with typical average diameter of 10 angstroms are under consideration.

  6. Ab Initio Modeling of Ecosystems with Artificial Life

    Adami, C.

    2002-01-01

    Artificial Life provides the opportunity to study the emergence and evolution of simple ecosystems in real time. We give an overview of the advantages and limitations of such an approach, as well as its relation to individual-based modeling techniques. The Digital Life system Avida is introduced and prospects for experiments with ab initio evolution (evolution "from scratch"), maintenance, as well as stability of ecosystems are discussed.

  7. Ab-initio calculations on melting of thorium

    Mukherjee, D.; Sahoo, B. D.; Joshi, K. D.; Kaushik, T. C.; Gupta, Satish C.

    2016-05-01

    Ab-initio molecular dynamics study has been performed on face centered cubic structured thorium to determine its melting temperature at room pressure. The ion-electron interaction potential energy calculated as a function of temperature for three volumes (a0)3 and (1.02a0)3 and (1.04a0)3 increases gradually with temperature and undergoes a sharp jump at ~2200 K, ~2100 K and ~1800 K, respectively. Here, a0 = 5.043 Å is the equilibrium lattice parameter at 0 K obtained from ab-initio calculations. These jumps in interaction energy are treated as due to the onset of melting and corresponding temperatures as melting point. The melting point of 2100 K is close to the experimental value of 2023K. Further, the same has been verified by plotting the atomic arrangement evolved at various temperatures and corresponding pair correlation functions.

  8. Augmented wave ab initio EFG calculations: some methodological warnings

    Errico, Leonardo A. [Departamento de Fisica-IFLP (CONICET), Facultad de Ciencias Exactas, Universidad Nacional de La Plata, CC67 (1900) La Plata (Argentina); Renteria, Mario [Departamento de Fisica-IFLP (CONICET), Facultad de Ciencias Exactas, Universidad Nacional de La Plata, CC67 (1900) La Plata (Argentina); Petrilli, Helena M. [Instituto de Fisica-DFMT, Universidade de Sao Paulo, C.P. 66318, 05315-970 Sao Paulo, SP (Brazil)]. E-mail: hmpetril@macbeth.if.usp.br

    2007-02-01

    We discuss some accuracy aspects inherent to ab initio electronic structure calculations in the understanding of nuclear quadrupole interactions. We use the projector augmented wave method to study the electric-field gradient (EFG) at both Sn and O sites in the prototype cases SnO and SnO{sub 2}. The term ab initio is used in the standard context of the also called first principles methods in the framework of the Density Functional Theory. As the main contributions of EFG calculations to problems in condensed matter physics are related to structural characterizations on the atomic scale, we discuss the 'state of the art' on theoretical EFG calculations and make a brief critical review on the subject, calling attention to some fundamental theoretical aspects.

  9. Augmented wave ab initio EFG calculations: some methodological warnings

    We discuss some accuracy aspects inherent to ab initio electronic structure calculations in the understanding of nuclear quadrupole interactions. We use the projector augmented wave method to study the electric-field gradient (EFG) at both Sn and O sites in the prototype cases SnO and SnO2. The term ab initio is used in the standard context of the also called first principles methods in the framework of the Density Functional Theory. As the main contributions of EFG calculations to problems in condensed matter physics are related to structural characterizations on the atomic scale, we discuss the 'state of the art' on theoretical EFG calculations and make a brief critical review on the subject, calling attention to some fundamental theoretical aspects

  10. Ab Initio Path to Heavy Nuclei

    Binder, Sven; Calci, Angelo; Roth, Robert

    2014-01-01

    We present the first ab initio calculations of nuclear ground states up into the domain of heavy nuclei, spanning the range from 16-O to 132-Sn based on two- plus three-nucleon interactions derived within chiral effective field theory. We employ the similarity renormalization group for preparing the Hamiltonian and use coupled-cluster theory to solve the many-body problem for nuclei with closed sub-shells. Through an analysis of theoretical uncertainties resulting from various truncations in this framework, we identify and eliminate the technical hurdles that previously inhibited the step beyond medium-mass nuclei, allowing for reliable validations of nuclear Hamiltonians in the heavy regime. Following this path we show that chiral Hamiltonians qualitatively reproduce the systematics of nuclear ground-state energies up to the neutron-rich Sn isotopes.