Towards SiC Surface Functionalization: An Ab Initio Study
Cicero, G; Catellani, A
2005-01-28
We present a microscopic model of the interaction and adsorption mechanism of simple organic molecules on SiC surfaces as obtained from ab initio molecular dynamics simulations. Our results open the way to functionalization of silicon carbide, a leading candidate material for bio-compatible devices.
Ab initio study of hydrogen on beryllium surfaces
NASA Astrophysics Data System (ADS)
Bachurin, D. V.; Vladimirov, P. V.
2015-11-01
Static ab initio calculations were performed for five principal hexagonal close-packed beryllium surfaces: basal, prismatic (type I and II) and pyramidal (type I and II). The basal plane was found to be the most energetically favorable, while the energies of the prismatic (type I) and pyramidal (type I) planes were slightly higher followed by the type II planes. Beryllium is known to show extreme interlayer distance relaxation near the surface. Up to five outermost atomic layers were involved in surface relaxation. The presence of hydrogen on the beryllium surfaces led to a noticeable reduction of the surface energy.
Ab initio Potential Energy Surface for H-H2
NASA Technical Reports Server (NTRS)
Partridge, Harry; Bauschlicher, Charles W., Jr.; Stallcop, James R.; Levin, Eugene
1993-01-01
Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- (mu)E(sub h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(sub 0) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.
Ab initio potential energy surface and rovibrational states of HBO
NASA Astrophysics Data System (ADS)
Ha, Tae-Kyu; Makarewicz, Jan
1999-01-01
The potential energy surface describing the large-amplitude motion of H around the BO core in the HBO molecule has been determined from ab initio calculations. This surface has been sampled by a set of 170 grid points from a two-dimensional space defined by the stretching and the bending coordinates of the H nucleus. At each grid point, the BO bond length has been optimized using the second-order Møller-Plesset perturbation theory with the basis set aug-cc-pVTZ. The surface has a local minimum for the linear as well as the bent configuration of HBO. A low energy barrier to the linear configuration BOH causes a large-amplitude motion and a strong rovibrational interaction in the molecule. Its rovibrational dynamics is different from the dynamics in bent or quasilinear triatomic molecules.
Ab Initio Potential Energy Surface for H-H2
NASA Technical Reports Server (NTRS)
Patridge, Harry; Bauschlicher, Charles W., Jr.; Stallcop, James R.; Levin, Eugene
1993-01-01
Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- 3 micro E(h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces (25-70 kcal/mol above the H-H2 asymptote) at small interatomic separations; the Boothroyd, Keogh, Martin, and Peterson (BKMP) potential energy surface is found to agree with results of the present calculations within the expected uncertainty (+/- 1 kcal/mol) of the fit. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(0)) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.
An Ab Initio Based Potential Energy Surface for Water
NASA Technical Reports Server (NTRS)
Partridge, Harry; Schwenke, David W.; Langhoff, Stephen R. (Technical Monitor)
1996-01-01
We report a new determination of the water potential energy surface. A high quality ab initio potential energy surface (PES) and dipole moment function of water have been computed. This PES is empirically adjusted to improve the agreement between the computed line positions and those from the HITRAN 92 data base. The adjustment is small, nonetheless including an estimate of core (oxygen 1s) electron correlation greatly improves the agreement with experiment. Of the 27,245 assigned transitions in the HITRAN 92 data base for H2(O-16), the overall root mean square (rms) deviation between the computed and observed line positions is 0.125/cm. However the deviations do not correspond to a normal distribution: 69% of the lines have errors less than 0.05/cm. Overall, the agreement between the line intensities computed in the present work and those contained in the data base is quite good, however there are a significant number of line strengths which differ greatly.
Exploring the free energy surface using ab initio molecular dynamics.
Samanta, Amit; Morales, Miguel A; Schwegler, Eric
2016-04-28
Efficient exploration of configuration space and identification of metastable structures in condensed phase systems are challenging from both computational and algorithmic perspectives. In this regard, schemes that utilize a set of pre-defined order parameters to sample the relevant parts of the configuration space [L. Maragliano and E. Vanden-Eijnden, Chem. Phys. Lett. 426, 168 (2006); J. B. Abrams and M. E. Tuckerman, J. Phys. Chem. B 112, 15742 (2008)] have proved useful. Here, we demonstrate how these order-parameter aided temperature accelerated sampling schemes can be used within the Born-Oppenheimer and the Car-Parrinello frameworks of ab initio molecular dynamics to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways. We have used these methods to identify the metastable structures and reaction pathways in SiO2 and Ti. In addition, we have used the string method [W. E, W. Ren, and E. Vanden-Eijnden, Phys. Rev. B 66, 052301 (2002); L. Maragliano et al., J. Chem. Phys. 125, 024106 (2006)] within the density functional theory to study the melting pathways in the high pressure cotunnite phase of SiO2 and the hexagonal closed packed to face centered cubic phase transition in Ti. PMID:27131525
Exploring the free energy surface using ab initio molecular dynamics
NASA Astrophysics Data System (ADS)
Samanta, Amit; Morales, Miguel A.; Schwegler, Eric
2016-04-01
Efficient exploration of configuration space and identification of metastable structures in condensed phase systems are challenging from both computational and algorithmic perspectives. In this regard, schemes that utilize a set of pre-defined order parameters to sample the relevant parts of the configuration space [L. Maragliano and E. Vanden-Eijnden, Chem. Phys. Lett. 426, 168 (2006); J. B. Abrams and M. E. Tuckerman, J. Phys. Chem. B 112, 15742 (2008)] have proved useful. Here, we demonstrate how these order-parameter aided temperature accelerated sampling schemes can be used within the Born-Oppenheimer and the Car-Parrinello frameworks of ab initio molecular dynamics to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways. We have used these methods to identify the metastable structures and reaction pathways in SiO2 and Ti. In addition, we have used the string method [W. E, W. Ren, and E. Vanden-Eijnden, Phys. Rev. B 66, 052301 (2002); L. Maragliano et al., J. Chem. Phys. 125, 024106 (2006)] within the density functional theory to study the melting pathways in the high pressure cotunnite phase of SiO2 and the hexagonal closed packed to face centered cubic phase transition in Ti.
CVRQD ab initio ground-state adiabatic potential energy surfaces for the water molecule.
Barletta, Paolo; Shirin, Sergei V; Zobov, Nikolai F; Polyansky, Oleg L; Tennyson, Jonathan; Valeev, Edward F; Császár, Attila G
2006-11-28
The high accuracy ab initio adiabatic potential energy surfaces (PESs) of the ground electronic state of the water molecule, determined originally by Polyansky et al. [Science 299, 539 (2003)] and called CVRQD, are extended and carefully characterized and analyzed. The CVRQD potential energy surfaces are obtained from extrapolation to the complete basis set of nearly full configuration interaction valence-only electronic structure computations, augmented by core, relativistic, quantum electrodynamics, and diagonal Born-Oppenheimer corrections. We also report ab initio calculations of several quantities characterizing the CVRQD PESs, including equilibrium and vibrationally averaged (0 K) structures, harmonic and anharmonic force fields, harmonic vibrational frequencies, vibrational fundamentals, and zero-point energies. They can be considered as the best ab initio estimates of these quantities available today. Results of first-principles computations on the rovibrational energy levels of several isotopologues of the water molecule are also presented, based on the CVRQD PESs and the use of variational nuclear motion calculations employing an exact kinetic energy operator given in orthogonal internal coordinates. The variational nuclear motion calculations also include a simplified treatment of nonadiabatic effects. This sophisticated procedure to compute rovibrational energy levels reproduces all the known rovibrational levels of the water isotopologues considered, H(2) (16)O, H(2) (17)O, H(2) (18)O, and D(2) (16)O, to better than 1 cm(-1) on average. Finally, prospects for further improvement of the ground-state adiabatic ab initio PESs of water are discussed. PMID:17144700
NASA Astrophysics Data System (ADS)
Geng, Hua Y.
2015-02-01
A multilevel approach to sample the potential energy surface in a path integral formalism is proposed. The purpose is to reduce the required number of ab initio evaluations of energy and forces in ab initio path integral molecular dynamics (AI-PIMD) simulation, without compromising the overall accuracy. To validate the method, the internal energy and free energy of an Einstein crystal are calculated and compared with the analytical solutions. As a preliminary application, we assess the performance of the method in a realistic model-the FCC phase of dense atomic hydrogen, in which the calculated result shows that the acceleration rate is about 3 to 4-fold for a two-level implementation, and can be increased up to 10 times if extrapolation is used. With only 16 beads used for the ab initio potential sampling, this method gives a well converged internal energy. The residual error in pressure is just about 3 GPa, whereas it is about 20 GPa for a plain AI-PIMD calculation with the same number of beads. The vibrational free energy of the FCC phase of dense hydrogen at 300 K is also calculated with an AI-PIMD thermodynamic integration method, which gives a result of about 0.51 eV/proton at a density of rs = 0.912.
Geng, Hua Y.
2015-02-15
A multilevel approach to sample the potential energy surface in a path integral formalism is proposed. The purpose is to reduce the required number of ab initio evaluations of energy and forces in ab initio path integral molecular dynamics (AI-PIMD) simulation, without compromising the overall accuracy. To validate the method, the internal energy and free energy of an Einstein crystal are calculated and compared with the analytical solutions. As a preliminary application, we assess the performance of the method in a realistic model—the FCC phase of dense atomic hydrogen, in which the calculated result shows that the acceleration rate is about 3 to 4-fold for a two-level implementation, and can be increased up to 10 times if extrapolation is used. With only 16 beads used for the ab initio potential sampling, this method gives a well converged internal energy. The residual error in pressure is just about 3 GPa, whereas it is about 20 GPa for a plain AI-PIMD calculation with the same number of beads. The vibrational free energy of the FCC phase of dense hydrogen at 300 K is also calculated with an AI-PIMD thermodynamic integration method, which gives a result of about 0.51 eV/proton at a density of r{sub s}=0.912.
Enhancing mechanical toughness of aluminum surfaces by nano-boron implantation: An ab initio study
NASA Astrophysics Data System (ADS)
Zhu, Zhen; Kwon, Dae-Gyeon; Kwon, Young-Kyun; Tománek, David
2015-01-01
Searching for ways to enhance surface hardness of aluminum, we study the equilibrium structure, stability, elastic properties and formation dynamics of a boron-enriched surface using ab initio density functional calculations. We used molecular dynamics simulations to model the implantation of energetic boron nanoparticles in Al and identify structural arrangements that optimize the formation of strong covalent Bsbnd Al bonds. Nano-indentation simulations based on constrained optimization suggest that presence of boron nanostructures in the subsurface region enhances significantly the mechanical hardness of aluminum surfaces.
Ab initio study of tungsten defects near the surface
NASA Astrophysics Data System (ADS)
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.
Surface Segregation Energies of BCC Binaries from Ab Initio and Quantum Approximate Calculations
NASA Technical Reports Server (NTRS)
Good, Brian S.
2003-01-01
We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy method. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameterization. Quantum approximate segregation energies are computed with and without atomistic relaxation. The ab initio calculations are performed without relaxation for the most part, but predicted relaxations from quantum approximate calculations are used in selected cases to compute approximate relaxed ab initio segregation energies. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with other quantum approximate and ab initio theoretical work, and available experimental results.
Exploring the free energy surface using ab initio molecular dynamics
Samanta, Amit; Morales, Miguel A.; Schwegler, Eric
2016-04-22
Efficient exploration of the configuration space and identification of metastable structures are challenging from both computational as well as algorithmic perspectives. Here, we extend the recently proposed orderparameter aided temperature accelerated sampling schemes to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways within the framework of density functional theory based molecular dynamics. The sampling method is applied to explore the relevant parts of the configuration space in prototypical materials SiO2 and Ti to identify the different metastable structures corresponding to different phases in these materials. In addition, we use the string method inmore » collective variables to study the melting pathways in the high pressure cotunnite phase of SiO2 and the hcp to fcc phase transition in Ti.« less
Gor, Gennady Y; Bernstein, Noam
2016-05-31
Adsorption-induced deformation is expansion or contraction of a solid due to adsorption on its surface. This phenomenon is important for a wide range of applications, from chemomechanical sensors to methane recovery from geological formations. The strain of the solid is driven by the change of the surface stress due to adsorption. Using ab initio molecular dynamics, we calculate the surface stresses for the dry α-quartz surfaces, and investigate how these stresses change when the surfaces are exposed to water. We find that the nonhydroxylated surface shows small and approximately isotropic changes in stress, while the hydroxylated surface, which interacts more strongly with the polar water molecules, shows larger and qualitatively anisotropic (opposite sign in xx and yy) surface stress changes. All of these changes are several times larger than the surface tension of water itself. The anisotropy and possibility of positive surface stress change can explain experimentally observed surface area contraction due to adsorption. PMID:27159032
Phosphine adsorption and dissociation on the Si(001) surface: An ab initio survey of structures
NASA Astrophysics Data System (ADS)
Warschkow, O.; Wilson, H. F.; Marks, N. A.; Schofield, S. R.; Curson, N. J.; Smith, P. V.; Radny, M. W.; McKenzie, D. R.; Simmons, M. Y.
2005-09-01
We report a comprehensive ab initio survey of possible dissociation intermediates of phosphine (PH3) on the Si(001) surface. We assign three scanning tunneling microscopy (STM) features, commonly observed in room-temperature dosing experiments, to PH2+H , PH+2H , and P+3H species, respectively, on the basis of calculated energetics and STM simulation. These assignments and a time series of STM images which shows these three STM features converting into another, allow us to outline a mechanism for the complete dissociation of phosphine on the Si(001) surface. This mechanism closes an important gap in the understanding of the doping process of semiconductor devices.
Ho, T.; Rabitz, H.
1996-02-01
A general interpolation method for constructing smooth molecular potential energy surfaces (PES{close_quote}s) from {ital ab} {ital initio} data are proposed within the framework of the reproducing kernel Hilbert space and the inverse problem theory. The general expression for an {ital a} {ital posteriori} error bound of the constructed PES is derived. It is shown that the method yields globally smooth potential energy surfaces that are continuous and possess derivatives up to second order or higher. Moreover, the method is amenable to correct symmetry properties and asymptotic behavior of the molecular system. Finally, the method is generic and can be easily extended from low dimensional problems involving two and three atoms to high dimensional problems involving four or more atoms. Basic properties of the method are illustrated by the construction of a one-dimensional potential energy curve of the He{endash}He van der Waals dimer using the exact quantum Monte Carlo calculations of Anderson {ital et} {ital al}. [J. Chem. Phys. {bold 99}, 345 (1993)], a two-dimensional potential energy surface of the HeCO van der Waals molecule using recent {ital ab} {ital initio} calculations by Tao {ital et} {ital al}. [J. Chem. Phys. {bold 101}, 8680 (1994)], and a three-dimensional potential energy surface of the H{sup +}{sub 3} molecular ion using highly accurate {ital ab} {ital initio} calculations of R{umlt o}hse {ital et} {ital al}. [J. Chem. Phys. {bold 101}, 2231 (1994)]. In the first two cases the constructed potentials clearly exhibit the correct asymptotic forms, while in the last case the constructed potential energy surface is in excellent agreement with that constructed by R{umlt o}hse {ital et} {ital al}. using a low order polynomial fitting procedure. {copyright} {ital 1996 American Institute of Physics.}
An Ab Initio Approach Towards Engineering Fischer-Tropsch Surface Chemistry
Matthew Neurock; David A. Walthall
2006-05-07
One of the greatest societal challenges over the next decade is the production of cheap, renewable energy for the 10 billion people that inhabit the earth. This will require the development of various different energy sources potentially including fuels derived from methane, coal, and biomass and alternatives sources such as solar, wind and nuclear energy. One approach will be to synthesize gasoline and other fuels from simpler hydrocarbons such as CO derived from methane or other U.S. based sources such as coal. Syngas (CO and H{sub 2}) can be readily converted into higher molecular weight hydrocarbons through Fischer-Tropsch synthesis. Fischer-Tropsch synthesis involves the initiation or activation of CO and H{sub 2} bonds, the subsequent propagation steps including hydrogenation and carbon-carbon coupling, followed by chain termination reactions. Commercially viable catalysts include supported Co and Co-alloys. Over the first two years of this project we have used ab initio methods to determine the adsorption energies for all reactants, intermediates, and products along with the overall reaction energies and their corresponding activation barriers over the Co(0001) surface. Over the third year of the project we developed and advanced an ab initio-based kinetic Monte Carlo simulation code to simulate Fischer Tropsch synthesis. This report details our work over the last year which has focused on the derivation of kinetic parameters for the elementary steps involved in FT synthesis from ab initio density functional theoretical calculations and the application of the kinetic Monte Carlo algorithm to simulate the initial rates of reaction for FT over the ideal Co(0001) surface. The results from our simulations over Co(0001) indicate the importance of stepped surfaces for the activation of adsorbed CO. In addition, they demonstrate that the dominant CH{sub x}* surface intermediate under steady state conditions is CH*. This strongly suggests that hydrocarbon coupling
Full-dimensional (15-dimensional) ab initio analytical potential energy surface for the H7+ cluster
NASA Astrophysics Data System (ADS)
Barragán, Patricia; Prosmiti, Rita; Wang, Yimin; Bowman, Joel M.
2012-06-01
Full-dimensional ab initio potential energy surface is constructed for the H_7^+ cluster. The surface is a fit to roughly 160 000 interaction energies obtained with second-order MöllerPlesset perturbation theory and the cc-pVQZ basis set, using the invariant polynomial method [B. J. Braams and J. M. Bowman, Int. Rev. Phys. Chem. 28, 577 (2009), 10.1080/01442350903234923]. We employ permutationally invariant basis functions in Morse-type variables for all the internuclear distances to incorporate permutational symmetry with respect to interchange of H atoms into the representation of the surface. We describe how different configurations are selected in order to create the database of the interaction energies for the linear least squares fitting procedure. The root-mean-square error of the fit is 170 cm-1 for the entire data set. The surface dissociates correctly to the H_5^+ + H2 fragments. A detailed analysis of its topology, as well as comparison with additional ab initio calculations, including harmonic frequencies, verify the quality and accuracy of the parameterized potential. This is the first attempt to present an analytical representation of the 15-dimensional surface of the H_7^+ cluster for carrying out dynamics studies.
The ab initio potential energy surface and spectroscopic constants of HOCl
NASA Astrophysics Data System (ADS)
Koput, Jacek; Peterson, Kirk A.
1998-02-01
The potential energy surface of hypochlorous acid, HOCl, has been determined from large-scale ab initio calculations using the coupled-cluster method CCSD(T), with basis sets of quadruple- and quintuple-zeta quality. The effect of core-electron correlation on the calculated structural parameters has been investigated. The vibrational-rotational energy levels of the three isotopic species of HOCl have then been calculated using the variational method and have been further characterized by the spectroscopic constants determined using the perturbational approach. The spectroscopic constants determined, are found to be in excellent agreement with experimental data.
A new ab initio potential energy surface for the Ne-H 2 interaction
NASA Astrophysics Data System (ADS)
Lique, François
2009-03-01
A new accurate three-dimensional potential energy surface for the Ne-H 2 system, which explicitly takes into account the r-dependence of the H 2 vibration, was determined from ab initio calculations. It was obtained with the single and double excitation coupled-cluster method with noniterative perturbational treatment of triple excitation [CCSD(T)]. Calculations was been performed using the augmented correlation-consistent polarized quintuple zeta basis set (aug-cc-pV5Z) for the three atoms. We checked the accuracy of the present ab initio calculations. We have determined, using the new Ne-H 2 potential energy surface, differential cross-sections for the rotational excitation of the H 2 and D 2 molecules in collision with Ne and we have compared them with experimental results of Faubel et al. [M. Faubel, F.A. Gianturco, F. Ragnetti, L.Y. Rusin, F. Sondermann, U. Tappe, J.P. Toennies, J. Chem. Phys. 101 (1994) 8800]. The overall agreement confirms that the new potential energy surface can be used for the simulation of molecular collisions and/or molecular spectroscopy of the van der Waals complex Ne-H 2.
NASA Astrophysics Data System (ADS)
Majumder, Moumita; Dawes, Richard; Wang, Xiao-Gang; Carrington, Tucker; Li, Jun; Guo, Hua; Manzhos, Sergei
2014-06-01
New potential energy surfaces for methane were constructed, represented as analytic fits to about 100,000 individual high-level ab initio data. Explicitly-correlated multireference data (MRCI-F12(AE)/CVQZ-F12) were computed using Molpro [1] and fit using multiple strategies. Fits with small to negligible errors were obtained using adaptations of the permutation-invariant-polynomials (PIP) approach [2,3] based on neural-networks (PIP-NN) [4,5] and the interpolative moving least squares (IMLS) fitting method [6] (PIP-IMLS). The PESs were used in full-dimensional vibrational calculations with an exact kinetic energy operator by representing the Hamiltonian in a basis of products of contracted bend and stretch functions and using a symmetry adapted Lanczos method to obtain eigenvalues and eigenvectors. Very close agreement with experiment was produced from the purely ab initio PESs. References 1- H.-J. Werner, P. J. Knowles, G. Knizia, 2012.1 ed. 2012, MOLPRO, a package of ab initio programs. see http://www.molpro.net. 2- Z. Xie and J. M. Bowman, J. Chem. Theory Comput 6, 26, 2010. 3- B. J. Braams and J. M. Bowman, Int. Rev. Phys. Chem. 28, 577, 2009. 4- J. Li, B. Jiang and Hua Guo, J. Chem. Phys. 139, 204103 (2013). 5- S Manzhos, X Wang, R Dawes and T Carrington, JPC A 110, 5295 (2006). 6- R. Dawes, X-G Wang, A.W. Jasper and T. Carrington Jr., J. Chem. Phys. 133, 134304 (2010).
Ab Initio Potential Energy Surfaces and the Calculation of Accurate Vibrational Frequencies
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Dateo, Christopher E.; Martin, Jan M. L.; Taylor, Peter R.; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within plus or minus 8 cm(exp -1) on average, and molecular bond distances are accurate to within plus or minus 0.001-0.003 Angstroms, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as vibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy will be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.
Ab initio calculations of the reaction pathways for methane decomposition over the Cu (111) surface
NASA Astrophysics Data System (ADS)
Gajewski, Grzegorz; Pao, Chun-Wei
2011-08-01
Growth of large-area, few-layer graphene has been reported recently through the catalytic decomposition of methane (CH4) over a Cu surface at high temperature. In this study, we used ab initio calculations to investigate the minimum energy pathways of successive dehydrogenation reactions of CH4 over the Cu (111) surface. The geometries and energies of all the reaction intermediates and transition states were identified using the climbing image nudged elastic band method. The activation barriers for CH4 decomposition over this Cu surface are much lower than those in the gas phase; furthermore, analysis of electron density differences revealed significant degrees of charge transfer between the adsorbates and the Cu atoms along the reaction path; these features reveal the role of Cu as the catalytic material for graphene growth. All the dehydrogenation reactions are endothermic, except for carbon dimer (C2) formation, which is, therefore, the most critical step for subsequent graphene growth, in particular, on Cu (111) surface.
Global ab initio ground-state potential energy surface of N4
NASA Astrophysics Data System (ADS)
Paukku, Yuliya; Yang, Ke R.; Varga, Zoltan; Truhlar, Donald G.
2013-07-01
We present a global ground-state potential energy surface for N4 suitable for treating high-energy vibrational-rotational energy transfer and collision-induced dissociation in N2-N2 collisions. To obtain the surface, complete active space second-order perturbation theory calculations were performed for the ground singlet state with an active space of 12 electrons in 12 orbitals and the maug-cc-pVTZ triple zeta basis set. About 17 000 ab initio data points have been calculated for the N4 system, distributed along nine series of N2 + N2 geometries and three series of N3 + N geometries. The six-dimensional ground-state potential energy surface is fitted using least-squares fits to the many-body component of the electronic energies based on permutationally invariant polynomials in bond order variables.
High-level ab initio potential energy surfaces and vibrational energies of H2CS
NASA Astrophysics Data System (ADS)
Yachmenev, Andrey; Yurchenko, Sergei N.; Ribeyre, Tristan; Thiel, Walter
2011-08-01
Six-dimensional (6D) potential energy surfaces (PESs) of H2CS have been generated ab initio using the recently proposed explicitly correlated (F12) singles and doubles coupled cluster method including a perturbational estimate of connected triple excitations, CCSD(T)-F12b [T. B. Adler, G. Knizia, and H.-J. Werner, J. Chem. Phys. 127, 221106 (2007)] in conjunction with F12-optimized correlation consistent basis sets. Core-electron correlation, high-order correlation, scalar relativistic, and diagonal Born-Oppenheimer terms were included as additive high-level (HL) corrections. The resulting 6D PESs were represented by analytical functions which were used in variational calculations of the vibrational term values below 5000 cm-1. The best PESs obtained with and without the HL corrections, VQZ-F12* HL and VQZ-F12*, reproduce the fundamental vibrational wavenumbers with mean absolute deviations of 1.13 and 1.22 cm-1, respectively. A detailed analysis of the effects of the HL corrections shows how the VQZ-F12 results benefit from error cancellation. The present purely ab initio PESs will be useful as starting points for empirical refinements towards an accurate "spectroscopic" PES of H2CS.
Han, Huixian; Li, Anyang; Guo, Hua
2014-12-28
A new full-dimensional global potential energy surface (PES) for the acetylene-vinylidene isomerization on the ground (S{sub 0}) electronic state has been constructed by fitting ∼37 000 high-level ab initio points using the permutation invariant polynomial-neural network method with a root mean square error of 9.54 cm{sup −1}. The geometries and harmonic vibrational frequencies of acetylene, vinylidene, and all other stationary points (two distinct transition states and one secondary minimum in between) have been determined on this PES. Furthermore, acetylene vibrational energy levels have been calculated using the Lanczos algorithm with an exact (J = 0) Hamiltonian. The vibrational energies up to 12 700 cm{sup −1} above the zero-point energy are in excellent agreement with the experimentally derived effective Hamiltonians, suggesting that the PES is approaching spectroscopic accuracy. In addition, analyses of the wavefunctions confirm the experimentally observed emergence of the local bending and counter-rotational modes in the highly excited bending vibrational states. The reproduction of the experimentally derived effective Hamiltonians for highly excited bending states signals the coming of age for the ab initio based PES, which can now be trusted for studying the isomerization reaction.
An Ab Initio Approach Towards Engineering Fischer-Tropsch Surface Chemistry
Matthew Neurock
2005-06-13
As petroleum prices continue to rise and the United States seeks to reduce its dependency on foreign oil, there is a renewed interest in the research and development of more efficient and alternative energy sources, such as fuel cells. One approach is to utilize processes that can produce long-chain hydrocarbons from other sources. One such reaction is Fischer-Tropsch synthesis. Fischer-Tropsch synthesis is a process by which syngas (CO and H{sub 2}) is converted to higher molecular weight hydrocarbons. The reaction involves a complex set of bond-breaking and bond-making reactions, such as CO and H{sub 2} activation, hydrocarbon hydrogenation reactions, and hydrocarbon coupling reactions. This report details our initial construction of an ab initio based kinetic Monte Carlo code that can be used to begin to simulate Fischer-Tropsch synthesis over model Co(0001) surfaces. The code is based on a stochastic kinetic formalism that allows us to explicitly track the transformation of all reactants, intermediates and products. The intrinsic kinetics for the simulations were derived from the ab initio results that we reported in previous year summaries.
A coupled channel study of HN2 unimolecular decay based on a global ab initio potential surface
NASA Technical Reports Server (NTRS)
Koizumi, Hiroyasu; Schatz, George C.; Walch, Stephen P.
1991-01-01
The unimolecular decay lifetimes of several vibrational states of HN2 are determined on the basis of an accurate coupled channel dynamics study using a global analytical potential surface. The surface reproduces the ab initio points with an rms error of 0.08 kcal/mol for energies below 20 kcal/mol. Modifications to the potential that describe the effect of improving the basis set in the ab initio calculations are provided. Converged coupled channel calculations are performed for the ground rotational state of HN2 to determine the lifetimes of the lowest ten vibrational states. Only the ground vibrational state (000) and first excited bend (001) are found to have lifetimes longer than 1 ps. The lifetimes of these states are estimated at 3 x 10 to the -9th and 2 x 10 to the -10th s, respectively. Variation of these results with quality of the ab initio calculations is not more than a factor of 5.
Yu, Dequan; Chen, Jun; Cong, Shulin; Sun, Zhigang
2015-12-17
The FH2– anion has a stable structure that resembles a configuration in the vicinity of the transition state for neutral reaction F + H2 → HF + H. Electron photodetachment spectra of the FH2– anion reveal the neutral reaction dynamics in the critical transition-state region. Accurate quantum dynamics simulations of the photodetachment spectra using highly accurate new ab initio potential energy surfaces for both anionic and neutral FH2 are performed and compared with all available experimental results. The results provide reliable interpretations for the experimental observations of FH2– photoelectron detachment and reveal a detailed picture of the molecular dynamics around the transition state of the F + H2 reaction. The latest high-resolution photoelectron detachment spectra [Kim et al. Science, 2015, 349, 510-513] confirm the high accuracy of our new potential energy surface for describing the resonance-enhanced reactivity of the neutral F + H2 reaction. PMID:26550683
NASA Astrophysics Data System (ADS)
Tokizaki, Chihiro; Yoshida, Takahiko; Takayanagi, Toshiyuki
2016-05-01
The cyclooctatetraene (COT) anion has a stable D4h structure that is similar to the transition state configurations of the neutral C-C bond-alternation (D4h ↔ D8h ↔ D4h) and ring-inversion (D2d ↔ D4h ↔ D2d) unimolecular reactions. The previously measured photodetachment spectrum of COT- revealed the reaction dynamics in the vicinity of the two transition states on the neutral potential energy surface. In this work, the photodetachment spectrum is calculated quantum mechanically on ab initio-level potential energy surfaces within a three degree-of-freedom reduced-dimensionality model. Very good agreement has been obtained between theory and experiment, providing reliable interpretations for the experimental spectrum. A detailed picture of the reactive molecular dynamics of the COT unimolecular reaction in the transition state region is also discussed.
Rotational Energy Transfer of N2 Determined Using a New Ab Initio Potential Energy Surface
NASA Technical Reports Server (NTRS)
Huo, Winifred M.; Stallcop, James R.; Partridge, Harry; Langhoff, Stephen R. (Technical Monitor)
1997-01-01
A new N2-N2 rigid-rotor surface has been determined using extensive Ab Initio quantum chemistry calculations together with recent experimental data for the second virial coefficient. Rotational energy transfer is studied using the new potential energy surface (PES) employing the close coupling method below 200 cm(exp -1) and coupled state approximation above that. Comparing with a previous calculation based on the PES of van der Avoird et al.,3 it is found that the new PES generally gives larger cross sections for large (delta)J transitions, but for small (delta)J transitions the cross sections are either comparable or smaller. Correlation between the differences in the cross sections and the two PES will be attempted. The computed cross sections will also be compared with available experimental data.
STM and ab initio study of holmium nanowires on a Ge(111) surface
NASA Astrophysics Data System (ADS)
Eames, C.; Bonet, C.; Probert, M. I. J.; Tear, S. P.; Perkins, E. W.
2006-11-01
A nanorod structure has been observed on the Ho/Ge(111) surface using scanning tunneling microscopy (STM). The rods do not require patterning of the surface or defects such as step edges in order to grow as is the case for nanorods on Si(111). At low holmium coverage the nanorods exist as isolated nanostructures while at high coverage they form a periodic 5×1 structure. We propose a structural model for the 5×1 unit cell and show using an ab initio calculation that the STM profile of our model structure compares favorably to that obtained experimentally for both filled and empty states sampling. The calculated local density of states shows that the nanorod is metallic in character.
Surface tension of ab initio liquid water at the water-air interface
NASA Astrophysics Data System (ADS)
Nagata, Yuki; Ohto, Tatsuhiko; Bonn, Mischa; Kühne, Thomas D.
2016-05-01
We report calculations on the surface tension of the water-air interface using ab initio molecular dynamics (AIMD) simulations. We investigate the influence of the cell size on surface tension of water from force field molecular dynamics simulations. We find that the calculated surface tension increases with increasing simulation cell size, thereby illustrating that a correction for finite size effects is essential for small systems that are customary in AIMD simulations. Moreover, AIMD simulations reveal that the use of a double-ζ basis set overestimates the experimentally measured surface tension due to the Pulay stress while more accurate triple and quadruple-ζ basis sets give converged results. We further demonstrate that van der Waals corrections critically affect the surface tension. AIMD simulations without the van der Waals correction substantially underestimate the surface tension while the van der Waals correction with the Grimme's D2 technique results in a value for the surface tension that is too high. The Grimme's D3 van der Waals correction provides a surface tension close to the experimental value. Whereas the specific choices for the van der Waals correction and basis sets critically affect the calculated surface tension, the surface tension is remarkably insensitive to the details of the exchange and correlation functionals, which highlights the impact of long-range interactions on the surface tension. Our simulated values provide important benchmarks, both for improving van der Waals corrections and AIMD simulations of aqueous interfaces.
NASA Astrophysics Data System (ADS)
Khaliullin, Rustam Z.; Eshet, Hagai; Kühne, Thomas D.; Behler, Jörg; Parrinello, Michele
2010-03-01
An interatomic potential for the diamond and graphite phases of carbon has been created using a neural-network (NN) representation of the ab initio potential energy surface. The NN potential combines the accuracy of a first-principles description of both phases with the efficiency of empirical force fields and allows one to perform a molecular-dynamics study, of ab initio quality, of the thermodynamics of graphite-diamond coexistence. Good agreement between the experimental and calculated coexistence curves is achieved if nuclear quantum effects are included in the simulation.
NASA Astrophysics Data System (ADS)
George, D. X. F.; Kumar, Sanjay
2010-08-01
Ab initio global adiabatic as well as quasidiabatic potential energy surfaces for the ground and the first excited electronic states of the H + + CO system have been computed as a function of the Jacobi coordinates ( R, r, γ) using Dunning's cc-pVTZ basis set at the internally contracted multi-reference (single and double) configuration interaction level of accuracy. In addition, nonadiabatic coupling matrix elements arising from radial motion, mixing angle and coupling potential have been computed using the ab initio procedure [Simah et al. (1999) [66
NASA Technical Reports Server (NTRS)
Jaffe, Richard L.; Pattengill, Merle D.; Schwenke, David W.
1989-01-01
Strategies for constructing global potential energy surfaces from a limited number of accurate ab initio electronic energy calculations are discussed. Generally, these data are concentrated in small regions of configuration space (e.g., in the vicinity of saddle points and energy minima) and difficulties arise in generating a potential function that is globally well-behaved. Efficient computer codes for carrying out classical trajectory calculations on vector and parallel processors are also described. Illustrations are given from recent work on the following chemical systems: Ca + HF yields CaF + H, H + H + H2 yields H2 + H2, N + O2 yields NO + O and O + N2 yields NO + N. The dynamics and kinetics of metathesis, dissociation, recombination, energy transfer and complex formation processes will be discussed.
Ab initio intermolecular potential energy surface and thermophysical properties of nitrous oxide
Crusius, Johann-Philipp Hassel, Egon; Hellmann, Robert Bich, Eckard
2015-06-28
We present an analytical intermolecular potential energy surface (PES) for two rigid nitrous oxide (N{sub 2}O) molecules derived from high-level quantum-chemical ab initio calculations. Interaction energies for 2018 N{sub 2}O–N{sub 2}O configurations were computed utilizing the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory using basis sets up to aug-cc-pVQZ supplemented with bond functions. A site-site potential function with seven sites per N{sub 2}O molecule was fitted to the pair interaction energies. We validated our PES by computing the second virial coefficient as well as shear viscosity and thermal conductivity in the dilute-gas limit. The values of these properties are substantiated by the best experimental data.
Ab initio intermolecular potential energy surface and thermophysical properties of nitrous oxide.
Crusius, Johann-Philipp; Hellmann, Robert; Hassel, Egon; Bich, Eckard
2015-06-28
We present an analytical intermolecular potential energy surface (PES) for two rigid nitrous oxide (N2O) molecules derived from high-level quantum-chemical ab initio calculations. Interaction energies for 2018 N2O-N2O configurations were computed utilizing the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory using basis sets up to aug-cc-pVQZ supplemented with bond functions. A site-site potential function with seven sites per N2O molecule was fitted to the pair interaction energies. We validated our PES by computing the second virial coefficient as well as shear viscosity and thermal conductivity in the dilute-gas limit. The values of these properties are substantiated by the best experimental data. PMID:26133428
Surface electron density models for accurate ab initio molecular dynamics with electronic friction
NASA Astrophysics Data System (ADS)
Novko, D.; Blanco-Rey, M.; Alducin, M.; Juaristi, J. I.
2016-06-01
Ab initio molecular dynamics with electronic friction (AIMDEF) is a valuable methodology to study the interaction of atomic particles with metal surfaces. This method, in which the effect of low-energy electron-hole (e-h) pair excitations is treated within the local density friction approximation (LDFA) [Juaristi et al., Phys. Rev. Lett. 100, 116102 (2008), 10.1103/PhysRevLett.100.116102], can provide an accurate description of both e-h pair and phonon excitations. In practice, its applicability becomes a complicated task in those situations of substantial surface atoms displacements because the LDFA requires the knowledge at each integration step of the bare surface electron density. In this work, we propose three different methods of calculating on-the-fly the electron density of the distorted surface and we discuss their suitability under typical surface distortions. The investigated methods are used in AIMDEF simulations for three illustrative adsorption cases, namely, dissociated H2 on Pd(100), N on Ag(111), and N2 on Fe(110). Our AIMDEF calculations performed with the three approaches highlight the importance of going beyond the frozen surface density to accurately describe the energy released into e-h pair excitations in case of large surface atom displacements.
Surface magnetism of the carbon foam: An ab-initio theoretical study
Yu, Guodong; Jiang, Liwei; Zheng, Yisong
2014-08-11
By means of ab-initio calculations, we study the surface magnetism of carbon foam, a carbon allotrope constructed by both the sp{sup 2} and sp{sup 3} hybridized carbon atoms. It is found that the sp{sup 3}-terminated surface and sp{sup 2}-terminated surface with hydrogen saturation (sp{sup 2}H) have the antiferromagnetic and ferromagnetic ground states, respectively. And the corresponding phase transition temperatures of these two surfaces are estimated to exceed 100 K. In addition, the magnetic moment of the sp{sup 2}H surface can be effectively modulated by strain. More interestingly, at the level of generalized gradient approximation, as the lattice are stretched, a bulk bandgap will open. And in such a bandgap, the sp{sup 2}H surface states can form spin polarized Dirac cones at K(K′) point. Thus, in the vicinity of the Fermi level, only the electron with a specific spin behaves like a Dirac Fermion.
Ab initio calculations for dissociative hydrogen adsorption on lithium oxide surfaces
Sutjianto, A. |; Tam, S.W.; Curtiss, L.A.; Johnson, C.E.; Pandey, R.
1994-12-01
Lithium ceramics are one class of materials being considered as tritium breeders for fusion technology,and hydrogen is known to enhance the release of tritium from lithium ceramic materials. Dissociative hydrogen chemisorption on the Li{sub 2}O surfaces of the (100), (110), and (111) planes has been investigated with ab initio Hartree-Fock calculations. Calculations for unrelaxed crystal Li{sub 2}O structures indicated that except for the (100) surface, the (110) and (111) surfaces are stable. Results on the heterolytic sites of n-layer (110) (where n {ge} 2) slabs and three-layer (111) slabs suggest that dissociative hydrogen chemisorption is endothermic. For a one-layer (110) slab at 100% surface coverage, the dissociative hydrogen chemisorption is exothermic, forming OH{sup {minus}} and Li{sup +}H{sup {minus}}Li{sup +}. The results also indicate that the low coordination environment in surface step structures, such as kinks and ledges, may plan an important role in the hydrogen chemisorption process. On the homolytic sites of the (110) and (111) surfaces, there is no hydrogen chemisorption.
NASA Astrophysics Data System (ADS)
Liu, Shi-Yu; Liu, Shiyang; Li, De-Jun; Wang, Sanwu; Guo, Jing; Shen, Yaogen
2015-02-01
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.
Ab initio thermodynamic evaluation of Pd atom interaction with CeO(2) surfaces.
Mayernick, Adam D; Janik, Michael J
2009-08-28
Palladium supported on ceria is an effective catalytic material for three-way automotive catalysis, catalytic combustion, and solid-oxide fuel cell (SOFC) anodes. The morphology, oxidation state, and particle size of Pd on ceria affect catalytic activity and are a function of experimental conditions. This work utilizes ab initio thermodynamics using density functional theory (DFT) (DFT+U) methods to evaluate the stability of Pd atoms, PdO(x) species, and small Pd particles in varying configurations on CeO(2) (111), (110), and (100) single crystal surfaces. Over specific oxygen partial pressure and temperature ranges, palladium incorporation to form a mixed surface oxide is thermodynamically favorable versus other single Pd atom states on each ceria surface. For example, Pd atoms may incorporate into Ce fluorite lattice positions in a Pd(4+) oxidation state on the CeO(2) (111) surface. The ceria support shifts the transition between formal Pd oxidation states (Pd(0), Pd(2+), Pd(4+)) relative to bulk palladium and stabilizes certain oxidized palladium species on each surface. We show that temperature, oxygen pressure, and cell potential in a SOFC can influence the stable states of palladium supported on ceria surfaces, providing insight into structural stability during catalytic operation. PMID:19725615
Liu, Shi-Yu; Liu, Shiyang; Li, De-Jun; Wang, Sanwu; Guo, Jing; Shen, Yaogen
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.
An Ab Initio Approach Towards Engineering Fischer-Tropsch Surface Chemistry
Matthew Neurock
2002-09-11
As the US seeks to develop an energy strategy that reduces the reliance on foreign oil, there is a renewed interest in research and development of the Fischer Tropsch synthesis of converting syngas into long chain hydrocarbon products. This report investigates some of the basic elementary steps for Fischer-Tropsch synthesis over ideal Co and Ru metal surfaces by using ab initio density functional theoretical calculations. This includes activation of CO of CO, the hydrogenation of CH{sub x} intermediates, and the adsorption and dissociation of water. The activation of CO is studied in detail showing a strong dependence on the surface coverage, defect sites and Co-Ru alloy formation. The barriers for CO activation over the ideal (0001) surfaces are quite high making CO activation at the terrace sites unlikely under operating conditions. The calculations for the overall reaction energies at the step edges indicate that these sites are much more reactive. The hydrogenation of the CHx intermediates occurs in a sequential fashion. CH1 was found to be the most stable intermediate over various surfaces. The barriers to form both CH* as well as CH{sub 4} are both found to be highly activated and potentially difficult steps. Water which is a reaction product was found to be weakly adsorbed on Co. Analysis of the microscopic reverse reaction of water activation indicates that this process has a very low activation barrier. Consequently, any water which forms desorbs or is activated to form surface hydroxyl intermediates.
Ab initio calculations of the reaction pathways for methane decomposition over the Cu (111) surface.
Gajewski, Grzegorz; Pao, Chun-Wei
2011-08-14
Growth of large-area, few-layer graphene has been reported recently through the catalytic decomposition of methane (CH(4)) over a Cu surface at high temperature. In this study, we used ab initio calculations to investigate the minimum energy pathways of successive dehydrogenation reactions of CH(4) over the Cu (111) surface. The geometries and energies of all the reaction intermediates and transition states were identified using the climbing image nudged elastic band method. The activation barriers for CH(4) decomposition over this Cu surface are much lower than those in the gas phase; furthermore, analysis of electron density differences revealed significant degrees of charge transfer between the adsorbates and the Cu atoms along the reaction path; these features reveal the role of Cu as the catalytic material for graphene growth. All the dehydrogenation reactions are endothermic, except for carbon dimer (C(2)) formation, which is, therefore, the most critical step for subsequent graphene growth, in particular, on Cu (111) surface. PMID:21842949
Probing the Si(001) surface with a Si tip: An ab initio study
NASA Astrophysics Data System (ADS)
Kantorovich, Lev; Hobbs, Chris
2006-06-01
Topographic noncontact atomic force microscopy (NC-AFM) images of the p(2×1) and c(4×2) reconstructions of the Si(001) surface are simulated for the cases of weak and strong tip-surface interactions and various temperatures using ab initio density functional theory. In the simulations the surface is imaged by a sharp silicon tip with a single dangling bond at its apex. At a very close approach to the surface, the tip flips a surface dimer when positioned close to its lower atom. The energy barriers for an individual flipped surface dimer to regain its initial configuration are calculated to be ˜0.1eV , implying that the surface should be able to “heal” itself at all but extremely low temperatures during one oscillation cycle of the cantilever. Thus, at small enough temperatures, T⩽70K , and large frequency shifts, the imaging process is dominated by tip induced dimer flip events resulting in a permanent deformation of the surface and an apparent p(2×1) symmetric phase to be observed. No dissipation is expected as the tip oscillations are conservative at these conditions. At intermediate temperatures, 70K⩽T⩽200K , the flipped dimers are able to return to the ground state during each tip oscillation, resulting in continuous healing of the surface and thus large dissipation is expected. At T⩾200K dimers flip back and forth easily resulting in an apparent symmetric p(2×1) phase and noticeable dissipation. At small frequency shifts the dimers do not flip, still the upper dimer atoms are imaged as bright so that surface reconstruction can easily be determined. The possibility of manipulating the orientation of dimers at low temperatures and large frequency shifts by means of preprogrammed scan directions, is also discussed.
MoS2 on an amorphous HfO2 surface: An ab initio investigation
NASA Astrophysics Data System (ADS)
Scopel, W. L.; Miwa, R. H.; Schmidt, T. M.; Venezuela, P.
2015-05-01
The energetic stability, electronic and structural properties of MoS2 adsorbed on an amorphous a-HfO2 surface (MoS2/HfO2) are examined through ab initio theoretical investigations. Our total energy results indicate that the formation of MoS2/HfO2 is an exothermic process with an adsorption energy of 34 meV/Å2, which means that it is more stable than similar systems like graphene/HfO2 and MoS2/SiO2. There are no chemical bonds at the MoS2-HfO2 interface. Upon formation of MoS2/HfO2, the electronic charge distribution is mostly localized at the interface region with no net charge transfer between the adsorbed MoS2 sheet and -HfO2 surface. However, the MoS2 sheet becomes n-type doped when there are oxygen vacancies in the HfO2 surface. Further investigation of the electronic distribution reveals that there are no electron- and hole-rich regions (electron-hole puddles) on the MoS2 sheet, which makes this system promising for use in high-speed nanoelectronic devices.
Optical properties of surfaces with supercell ab initio calculations: Local-field effects
NASA Astrophysics Data System (ADS)
Tancogne-Dejean, Nicolas; Giorgetti, Christine; Véniard, Valérie
2015-12-01
Surface optical and electronic properties are crucial for material science and have implications in fields as various as nanotechnology, nonlinear optics, and spectroscopies. In particular, the huge variation of electronic density perpendicular to the surface is expected to play a key role in absorption due to local-field effects. Numerous state-of-the-art theoretical and numerical ab initio formalisms developed for studying these properties are based on supercell approaches, in reciprocal space, due to their efficiency. In this paper, we show that the standard scheme fails for the out-of-plane optical response of the surface. This response is interpreted using the "effective-medium theory" with vacuum and also in terms of interaction between replicas, as the supercell approach implies a periodicity which is absent in the real system. We propose an alternative formulation, also based on the supercell, for computing the macroscopic dielectric function. Application to the clean Si(001) 2 ×1 surface allows us to present the effect of the local fields for both peak positions and line shape on the bulk and surface contributions. It shows how local fields built up for the in-plane and out-of-plane dielectric responses of the surface. In addition to their conceptual impact, our results explain why the standard approach gives reliable predictions for the in-plane components, leading to correct reflectance anisotropy spectra. Our scheme can be further generalized to other low-dimensional geometries, such as clusters or nanowires, and open the way to nonlinear optics for surfaces.
Potential energy surface and second virial coefficient of methane-water from ab initio calculations
NASA Astrophysics Data System (ADS)
Akin-Ojo, Omololu; Szalewicz, Krzysztof
2005-10-01
Six-dimensional intermolecular potential energy surfaces (PESs) for the interaction of CH4 with H2O are presented, obtained from ab initio calculations using symmetry-adapted perturbation theory (SAPT) at two different levels of intramonomer correlation and the supermolecular approach at three different levels of electron correlation. Both CH4 and H2O are assumed to be rigid molecules with interatomic distances and angles fixed at the average values in the ground-state vibration. A physically motivated analytical expression for each PES has been developed as a sum of site-site functions. The PES of the CH4-H2O dimer has only two symmetry-distinct minima. From the SAPT calculations, the global minimum has an energy of -1.03kcal /mol at a geometry where H2O is the proton donor, HO -H⋯CH4, with the O-H-C angle of 165°, while the secondary minimum, with an energy of -0.72kcal/mol, has CH4 in the role of the proton donor (H3C -H⋯OH2). We estimated the complete basis set limit of the SAPT interaction energy at the global minimum to be -1.06kcal/mol. The classical cross second virial coefficient B12(T) has been calculated for the temperature range 298-653K. Our best results agree well with some experiments, allowing an evaluation of the quality of experimental results.
Direct molecular simulation of nitrogen dissociation based on an ab initio potential energy surface
Valentini, Paolo Schwartzentruber, Thomas E. Bender, Jason D. Nompelis, Ioannis Candler, Graham V.
2015-08-15
The direct molecular simulation (DMS) approach is used to predict the internal energy relaxation and dissociation dynamics of high-temperature nitrogen. An ab initio potential energy surface (PES) is used to calculate the dynamics of two interacting nitrogen molecules by providing forces between the four atoms. In the near-equilibrium limit, it is shown that DMS reproduces the results obtained from well-established quasiclassical trajectory (QCT) analysis, verifying the validity of the approach. DMS is used to predict the vibrational relaxation time constant for N{sub 2}–N{sub 2} collisions and its temperature dependence, which are in close agreement with existing experiments and theory. Using both QCT and DMS with the same PES, we find that dissociation significantly depletes the upper vibrational energy levels. As a result, across a wide temperature range, the dissociation rate is found to be approximately 4–5 times lower compared to the rates computed using QCT with Boltzmann energy distributions. DMS calculations predict a quasi-steady-state distribution of rotational and vibrational energies in which the rate of depletion of high-energy states due to dissociation is balanced by their rate of repopulation due to collisional processes. The DMS approach simulates the evolution of internal energy distributions and their coupling to dissociation without the need to precompute rates or cross sections for all possible energy transitions. These benchmark results could be used to develop new computational fluid dynamics models for high-enthalpy flow applications.
Liu, Hanchao; Wang, Yimin; Bowman, Joel M.
2015-05-21
The calculation and characterization of the IR spectrum of liquid water have remained a challenge for theory. In this paper, we address this challenge using a combination of ab initio approaches, namely, a quantum treatment of IR spectrum using the ab initio WHBB water potential energy surface and a refined ab initio dipole moment surface. The quantum treatment is based on the embedded local monomer method, in which the three intramolecular modes of each embedded H{sub 2}O monomer are fully coupled and also coupled singly to each of six intermolecular modes. The new dipole moment surface consists of a previous spectroscopically accurate 1-body dipole moment surface and a newly fitted ab initio intrinsic 2-body dipole moment. A detailed analysis of the new dipole moment surface in terms of the coordinate dependence of the effective atomic charges is done along with tests of it for the water dimer and prism hexamer double-harmonic spectra against direct ab initio calculations. The liquid configurations are taken from previous molecular dynamics calculations of Skinner and co-workers, using the TIP4P plus E3B rigid monomer water potential. The IR spectrum of water at 300 K in the range of 0–4000 cm{sup −1} is calculated and compared with experiment, using the ab initio WHBB potential and new ab initio dipole moment, the q-TIP4P/F potential, which has a fixed-charged description of the dipole moment, and the TTM3-F potential and dipole moment surfaces. The newly calculated ab initio spectrum is in very good agreement with experiment throughout the above spectral range, both in band positions and intensities. This contrasts to results with the other potentials and dipole moments, especially the fixed-charge q-TIP4P/F model, which gives unrealistic intensities. The calculated ab initio spectrum is analyzed by examining the contribution of various transitions to each band.
NASA Astrophysics Data System (ADS)
Liu, Hanchao; Wang, Yimin; Bowman, Joel M.
2015-05-01
The calculation and characterization of the IR spectrum of liquid water have remained a challenge for theory. In this paper, we address this challenge using a combination of ab initio approaches, namely, a quantum treatment of IR spectrum using the ab initio WHBB water potential energy surface and a refined ab initio dipole moment surface. The quantum treatment is based on the embedded local monomer method, in which the three intramolecular modes of each embedded H2O monomer are fully coupled and also coupled singly to each of six intermolecular modes. The new dipole moment surface consists of a previous spectroscopically accurate 1-body dipole moment surface and a newly fitted ab initio intrinsic 2-body dipole moment. A detailed analysis of the new dipole moment surface in terms of the coordinate dependence of the effective atomic charges is done along with tests of it for the water dimer and prism hexamer double-harmonic spectra against direct ab initio calculations. The liquid configurations are taken from previous molecular dynamics calculations of Skinner and co-workers, using the TIP4P plus E3B rigid monomer water potential. The IR spectrum of water at 300 K in the range of 0-4000 cm-1 is calculated and compared with experiment, using the ab initio WHBB potential and new ab initio dipole moment, the q-TIP4P/F potential, which has a fixed-charged description of the dipole moment, and the TTM3-F potential and dipole moment surfaces. The newly calculated ab initio spectrum is in very good agreement with experiment throughout the above spectral range, both in band positions and intensities. This contrasts to results with the other potentials and dipole moments, especially the fixed-charge q-TIP4P/F model, which gives unrealistic intensities. The calculated ab initio spectrum is analyzed by examining the contribution of various transitions to each band.
NASA Astrophysics Data System (ADS)
Seller, Harrell; Shustorovich, Evgeny
1996-02-01
Binding energies for sulfur oxides, SO x, x = 1-3, have been determined for several coordination modes on silver, gold and palladium surfaces employing ab initio quantum chemical methods and the bond order conservation Morse potential (BOC-MP) method. SO 2 coordination was studied in the most detail. In general the agreement between the BOC-MP and ab initio binding energies is good for the (111) surfaces of silver and palladium with both methods predicting that, in the zero coverage limit, di-coordination via S,O and O,O will be more favorable energetically than mono-coordination via S. In the case of chemisorption on the Pd (110) surface the two methods agree well for the cases in which there are formulas for the BOC-MP binding energies. In going from the (111) surfaces to the (110) surfaces of silver and palladium the ab initio calculations predict that the preferred chemisorption site shifts from the bridge site to the hollow site. On the silver surfaces the net charge transferred to the adsorbate as judged from the Mulliken populations correlates roughly with the binding energy. No significant charge transfer was found on the palladium surfaces. Our SO 2 chemisorption calculations indicate that the work functions of the metal surfaces examined should increase upon mono-S adsorption, increase to a lesser extent upon di S,O adsorption and may even decrease upon di O,O adsorption. Ab initio calculations provide evidence of the existence of SO 2 surface dimers. The binding energy predicted by the BOC-MP model for SO 3 in the bridging site agrees well with the ab initio result for SO 3 di-coordinated in the long bridge of the Ag(110) surface. The methods yield similar predictions for the case of SO on silver. Our modeling provides a coherent picture consistent with many aspects of the experimental literature. We present some model predictions, particularly the di O,O coordination mode for SO 2, that require verification experimentally.
An Ab Initio Approach Towards Engineering Fischer-Tropsch Surface Chemistry
Matthew Neurock; Siddharth Chopra
2003-09-11
As the US seeks to develop an energy strategy that reduces the reliance on foreign oil, there is a renewed interest in the research and development of the Fischer Tropsch synthesis for converting syngas into long chain hydrocarbon products. This report investigates some of the basic elementary steps for Fischer-Tropsch synthesis over ideal Pt, Ru and carbon-covered Pt and Ru metal surfaces by using ab initio density functional theoretical calculations. We examine in detail the adsorption sites as well as the binding energies for C, CH, CH{sub 2}, CH3 and CH4 on Pt(111), Ru(0001), 2x2-C-Pt(111) and 2x2-C-Ru(0001). The results indicate that the binding energies increase with decreasing the hydrogen in the fragment molecule, i.e. CH{sub 4} < CH{sub 3} < CH{sub 2} < CH < C. More specifically the work analyzes the elementary steps involved in the activation of methane. This is simply the reverse set of steps necessary for the hydrogenation of C to CH{sub 4}. The results indicate that these hydrocarbon intermediates bind more strongly to Ru than Pt. The introduction of co-adsorbed carbon atoms onto both Ru(0001) as well as Pt(111) significantly increased the overall energies as well as the activation barriers for C-H bond activation. The results suggest that Ru may be so active that it initially can initially activate CH4 into CH or C but ultimately it dies because the CH and C intermediates poison the surface and thus kill its activity. Methane can dissociate on Pt but subsequent hydrocarbon coupling reactions act to remove the surface carbon.
NASA Astrophysics Data System (ADS)
Vlahos, Vasilios; Lee, Yueh-Lin; Booske, John H.; Morgan, Dane; Turek, Ladislav; Kirshner, Mark; Kowalczyk, Richard; Wilsen, Craig
2009-05-01
Scandate cathodes (BaxScyOz on W) are important thermionic electron emission materials whose emission mechanism remains unclear. Ab initio modeling is used to investigate the surface properties of both scandate and traditional B-type (Ba-O on W) cathodes. We demonstrate that the Ba-O dipole surface structure believed to be present in active B-type cathodes is not thermodynamically stable, suggesting that a nonequilibrium steady state dominates the active cathode's surface structure. We identify a stable, low work function BaxScyOz surface structure, which may be responsible for some scandate cathode properties and demonstrate that multicomponent surface coatings can lower cathode work functions.
An ab initio potential energy surface and dynamics of the Ar+H2+ → ArH + + H reaction
NASA Astrophysics Data System (ADS)
Liu, Xinguo; Liu, Huirong; Zhang, Qinggang
2011-04-01
An ab initio potential energy surface (PES) for the ground state (1 2A') of the chemical reaction Ar+H2+ → ArH + + H has been constructed from a set of accurate ab initio data, which we have computed using the coupled-cluster theory including all single and double excitations plus perturbative corrections for the triples UCCSD(T) with a large orbital basis set of aug-cc-pV5Z. The new PES has a root-mean-square (rms) error of 0.5341 kcal/mol. The total integral reaction cross-sections have been calculated at three collision energies by means of the quasi-classical trajectory (QCT) calculation based on the new PES and compared with previous TSH results.
Quantal Study of the Exchange Reaction for N + N2 using an ab initio Potential Energy Surface
NASA Technical Reports Server (NTRS)
Wang, Dunyou; Stallcop, James R.; Huo, Winifred M.; Dateo, Christopher E.; Schwenke, David W.; Partridge, Harry; Kwak, Dochan (Technical Monitor)
2002-01-01
The N + N2 exchange rate is calculated using a time-dependent quantum dynamics method on a newly determined ab initio potential energy surface (PES) for the ground A" state. This ab initio PES shows a double barrier feature in the interaction region with the barrier height at 47.2 kcal/mol, and a shallow well between these two barriers, with the minimum at 43.7 kcal/mol. A quantum dynamics wave packet calculation has been carried out using the fitted PES to compute the cumulative reaction probability for the exchange reaction of N + N2(J=O). The J - K shift method is then employed to obtain the rate constant for this reaction. The calculated rate constant is compared with experimental data and a recent quasi-classical calculation using a LEPS PES. Significant differences are found between the present and quasiclassical results. The present rate calculation is the first accurate 3D quantal dynamics study for N + N2 reaction system and the ab initio PES reported here is the first such surface for N3.
A computational ab initio study of surface diffusion of sulfur on the CdTe (111) surface
NASA Astrophysics Data System (ADS)
Naderi, Ebadollah; Ghaisas, S. V.
2016-08-01
In order to discern the formation of epitaxial growth of CdS shell over CdTe nanocrystals, kinetics related to the initial stages of the growth of CdS on CdTe is investigated using ab-initio methods. We report diffusion of sulfur adatom on the CdTe (111) A-type (Cd-terminated) and B-type (Te-terminated) surfaces within the density functional theory (DFT). The barriers are computed by applying the climbing Nudge Elastic Band (c-NEB) method. From the results surface hopping emerges as the major mode of diffusion. In addition, there is a distinct contribution from kick-out type diffusion in which a CdTe surface atom is kicked out from its position and is replaced by the diffusing sulfur atom. Also, surface vacancy substitution contributes to the concomitant dynamics. There are sites on the B- type surface that are competitively close in terms of the binding energy to the lowest energy site of epitaxy on the surface. The kick-out process is more likely for B-type surface where a Te atom of the surface is displaced by a sulfur adatom. Further, on the B-type surface, subsurface migration of sulfur is indicated. Furthermore, the binding energies of S on CdTe reveal that on the A-type surface, epitaxial sites provide relatively higher binding energies and barriers than on B-type.
Skokov, S.; Peterson, K.A.; Bowman, J.M.
1998-08-01
Accurate {ital ab initio} multireference configuration interaction (CI) calculations with large correlation-consistent basis sets are performed for HOCl. After extrapolation to the complete basis set limit, the {ital ab initio} data are precisely fit to give a semiglobal three-dimensional potential energy surface to describe HOCl{r_arrow}Cl+OH from high overtone excitation of the OH-stretch. The average absolute deviation between the {ital ab initio} and fitted energies is 4.2thinspcm{sup {minus}1} for energies up to 60 kcal/mol relative to the HOCl minimum. Vibrational energies of HOCl including the six overtones of the OH-stretch are computed using a vibrational-Cl method on the fitted potential and also on a slightly adjusted potential. Near-spectroscopic accuracy is obtained using the adjusted potential; the average absolute deviation between theory and experiment for 19 experimentally reported states is 4.8thinspcm{sup {minus}1}. Very good agreement with experiment is also obtained for numerous rotational energies for the ground vibrational state, the ClO-stretch fundamental, and the fifth overtone of the OH-stretch. {copyright} {ital 1998 American Institute of Physics.}
An Ab Initio Approach Towards Engineering Fischer-Tropsch Surface Chemistry
Matthew Neurock
2006-09-11
One of the greatest societal challenges over the next decade is the production of cheap, renewable energy for the 10 billion people that inhabit the earth. This will require the development of various energy sources which will likely include fuels derived from methane, coal, and biomass and alternatives sources such as solar, wind and nuclear energy. One approach will be to synthesize gasoline and other fuels from simpler hydrocarbons such as CO derived from methane or other U.S. based sources such as coal. Syngas (CO and H{sub 2}) can be readily converted into higher molecular weight hydrocarbons through Fischer-Tropsch synthesis. Fischer-Tropsch (FT) synthesis involves the adsorption and the activation of CO and H{sub 2}, the subsequent propagation steps including hydrogenation and carbon-carbon coupling, followed by chain termination reactions. The current commercial catalysts are supported Co and Co-alloys particles. This project set out with the following objectives in mind: (1) understand the reaction mechanisms that control FT kinetics, (2) predict how the intrinsic metal-adsorbate bond affects the sequence of elementary steps in FT, (3) establish the effects of the reaction environment on catalytic activity and selectivity, (4) construct a first-principles based algorithm that can incorporate the detailed atomic surface structure and simulate the kinetics for the myriad of elementary pathways that make up FT chemistry, and (5) suggest a set of optimal features such as alloy composition and spatial configuration, oxide support, distribution of defect sites. As part of this effort we devoted a significant portion of time to develop an ab initio based kinetic Monte Carlo simulation which can be used to follow FT surface chemistry over different transition metal and alloy surfaces defined by the user. Over the life of this program, we have used theory and have developed and applied stochastic Monte Carlo simulations in order to establish the fundamental
Zhang, Chunfang; Fu, Mingkai; Shen, Zhitao; Ma, Haitao; Bian, Wensheng
2014-06-21
A new global ab initio potential energy surface (called ZMB-a) for the 1(1)A' state of the C((1)D)+H2 reactive system has been constructed. This is based upon ab initio calculations using the internally contracted multireference configuration interaction approach with the aug-cc-pVQZ basis set, performed at about 6300 symmetry unique geometries. Accurate analytical fits are generated using many-body expansions with the permutationally invariant polynomials, except that the fit of the deep well region is taken from our previous fit. The ZMB-a surface is unique in the accurate description of the regions around conical intersections (CIs) and of van der Waals (vdW) interactions. The CIs between the 1(1)A' and 2(1)A' states cause two kinds of barriers on the ZMB-a surface: one is in the linear H-CH dissociation direction with a barrier height of 9.07 kcal/mol, which is much higher than those on the surfaces reported before; the other is in the C((1)D) collinearly attacking H2 direction with a barrier height of 12.39 kcal/mol. The ZMB-a surface basically reproduces our ab initio calculations in the vdW interaction regions, and supports a linear C-HH vdW complex in the entrance channel, and two vdW complexes in the exit channel, at linear CH-H and HC-H geometries, respectively. PMID:24952535
Ab initio atomic recombination reaction energetics on model heat shield surfaces
NASA Technical Reports Server (NTRS)
Senese, Fredrick; Ake, Robert
1992-01-01
Ab initio quantum mechanical calculations on small hydration complexes involving the nitrate anion are reported. The self-consistent field method with accurate basis sets has been applied to compute completely optimized equilibrium geometries, vibrational frequencies, thermochemical parameters, and stable site labilities of complexes involving 1, 2, and 3 waters. The most stable geometries in the first hydration shell involve in-plane waters bridging pairs of nitrate oxygens with two equal and bent hydrogen bonds. A second extremely labile local minimum involves out-of-plane waters with a single hydrogen bond and lies about 2 kcal/mol higher. The potential in the region of the second minimum is extremely flat and qualitatively sensitive to changes in the basis set; it does not correspond to a true equilibrium structure.
Wong, Kim F.; Sonnenberg, Jason L.; Paesani, Francesco; Yamamoto, Takeshi; Vaníček, Jiří; Zhang, Wei; Schlegel, H. Bernhard; Case, David A.; Cheatham, Thomas E.; Miller, William H.; Voth, Gregory A.
2010-01-01
The rates of intramolecular proton transfer are calculated on a full-dimensional reactive electronic potential energy surface that incorporates high level ab initio calculations along the reaction path and by using classical Transition State theory, Path-Integral Quantum Transition State Theory, and the Quantum Instanton approach. The specific example problem studied is malonaldehyde. Estimates of the kinetic isotope effect using the latter two methods are found to be in reasonable agreement with each other. Improvements and extensions of this practical, yet chemically accurate framework for the calculations of quantized, reactive dynamics are also discussed. PMID:21116485
NASA Astrophysics Data System (ADS)
Ri, Mun-Hyok; Yu, Chol-Jun; Jang, Yong-Man; Kim, Song-Un
2016-03-01
We report a computational study of the adsorption of zoledronic acid molecule on hydroxyapatite (001) surface within ab initio density functional theory. The systematic study has been performed, from hydroxyapatite bulk and surface, and zoledronic acid molecule to the adsorption of the molecule on the surface. The optimized bond lengths and bond angles were obtained and analyzed, giving an evidence of structural similarity between subjects under study. The formation energies of hydroxyapatite (001) surfaces with two kinds of terminations were computed as about 1.2 and 1.5 J/m^2 with detailed atomistic structural information. We determined the adsorption energies of zoledronic acid molecule on the surfaces, which are -260 kJ/mol at 0.25 ML and -400 kJ/mol at 0.5 ML. An atomistic insight of strong binding affinity of zoledronic acid to the hydroxyapatite surface was given and discussed.
AB INITIO AND CALPHAD THERMODYNAMICS OF MATERIALS
Turchi, P A
2004-04-14
Ab initio electronic structure methods can supplement CALPHAD in two major ways for subsequent applications to stability in complex alloys. The first one is rather immediate and concerns the direct input of ab initio energetics in CALPHAD databases. The other way, more involved, is the assessment of ab initio thermodynamics {acute a} la CALPHAD. It will be shown how these results can be used within CALPHAD to predict the equilibrium properties of multi-component alloys.
NASA Astrophysics Data System (ADS)
Song, Yu-Zhi; Li, Yong-Qing; Gao, Shou-Bao; Meng, Qing-Tian
2014-01-01
A globally accurate potential energy surface is reported for the electronic ground-state HLi2 by fitting ab initio energies to double many-body expansion formalism. The total 3726 ab initio energies used to map the HLi2 potential energy surface are calculated using the multi-reference configuration interaction method, with their dynamical correlation being semiempirically corrected by the double many-body expansion-scaled external correlation method. The current potential energy surface generates an excellent fit of the ab initio energies, showing a small root-mean squared derivation of 0.636 kcal mol-1. The topographical features of the HLi2 potential energy surface are examined in detail, which concludes that the H + Li2(X 1Σg) → Li + LiH(X 1Σ) reaction is essentially barrierless and the exothermicity is calculated to be 33.668 kcal mol-1, thus corroborates the available experimental and theoretical results.
NASA Astrophysics Data System (ADS)
Xavier, F. George D.; Kumar, Sanjay
2010-10-01
Ab initio global adiabatic and quasidiabatic potential energy surfaces of lowest four electronic (1-4 A3″) states of the H++O2 system have been computed in the Jacobi coordinates (R,r,γ) using Dunning's cc-pVTZ basis set at the internally contracted multireference (single and double) configuration interaction level of accuracy, which are relevant to the dynamics studies of inelastic vibrational and charge transfer processes observed in the scattering experiments. The computed equilibrium geometry parameters of the bound [HO2]+ ion in the ground electronic state and other parameters for the transition state for the isomerization process, HOO+⇌OOH+ are in good quantitative agreement with those available from the high level ab initio calculations, thus lending credence to the accuracy of the potential energy surfaces. The nonadiabatic couplings between the electronic states have been analyzed in both the adiabatic and quasidiabatic frameworks by computing the nonadiabatic coupling matrix elements and the coupling potentials, respectively. It is inferred that the dynamics of energy transfer processes in the scattering experiments carried out in the range of 9.5-23 eV would involve all the four electronic states.
Yousaf, Masood; Shin, Dongbin; Ruoff, Rodney; Park, Noejung
2015-12-17
We used ab initio molecular dynamics (AIMD) to investigate the effect of a monochromatic oscillating electric field in resonance with a particular molecular vibration on surfaces. As a case study, AIMD simulations were carried out for hydroxyl functional groups on graphene. When the frequency of the applied field matches with the C-OH vibration frequency, the amplitude is monotonically amplified, leading to a complete desorption from the surface, overcoming the substantial barrier. This suggests the possibility of activating a particular bond without damaging the remaining surface. We extended this work to the case of the amination of sp(2)-bonded carbon surfaces and discussed the general perspective that, in general, an unfavorable chemical process can be activated by applying an external electric field with an appropriate resonance frequency. PMID:26634785
Feng, Ya-Juan; Huang, Teng; Wang, Chao; Liu, Yi-Rong; Jiang, Shuai; Miao, Shou-Kui; Chen, Jiao; Huang, Wei
2016-07-14
Molecular level insight into the interaction between volatile organic compounds (VOCs) and aerosols is crucial for improvement of atmospheric chemistry models. In this paper, the interaction between adsorbed toluene, one of the most significant VOCs in the urban atmosphere, and the aqueous surface of aerosols was studied by means of combined molecular dynamics simulations and ab initio quantum chemistry calculations. It is revealed that toluene can be stably adsorbed on the surface of aqueous droplets via hydroxyl-π hydrogen bonding between the H atoms of the water molecules and the C atoms in the aromatic ring. Further, significant modifications on the electrostatic potential map and frontier molecular orbital are induced by the solvation effect of surface water molecules, which would affect the reactivity and pathway of the atmospheric photooxidation of toluene. This study demonstrates that the surface interactions should be taken into consideration in the atmospheric chemical models on oxidation of aromatics. PMID:27280740
NASA Astrophysics Data System (ADS)
Molina Sanchez, Alejandro; Wirtz, Ludger
2014-03-01
The interaction of graphene with substrates can alter its electronic and vibrational properties and is relevant for the practical use of graphene. In this work, we describe the graphene-substrate interaction through the theoretical study of the vibrational properties. We focus on three paradigmatic cases where the interaction strength changes gradually: graphene@BN, graphene@Ir(111), and graphene@SiC. We use ab initio methods to obtain the phonon band structure, the density of states, and the strength of the electron-phonon coupling. Graphene on boron nitride exhibits a weak interaction but a non-negligible shift of the 2D Raman band. We explain this observation by connecting the increase of the dielectric screening with the softening of the electron-phonon interaction. Graphene on iridium, also displays weak interaction but the substrate is a metal. In this case the electron-electron interaction in graphene is screened by a metal electron gas. In the last case, we study the buffer layer of graphene on silicon carbide. The strong hybridization of graphene with silicon carbide changes substantially the electronic structure of graphene. All the calculations are compared to experimental data of Raman spectroscopy and angle-resolved inelastic electron scattering. University of Luxembourg.
Hellmann, Robert; Bich, Eckard; Vogel, Eckhard
2008-06-01
A six-dimensional potential energy hypersurface (PES) for two interacting rigid methane molecules was determined from high-level quantum-mechanical ab initio computations. A total of 272 points for 17 different angular orientations on the PES were calculated utilizing the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory with basis sets of aug-cc-pVTZ and aug-cc-pVQZ qualities. The calculated interaction energies were extrapolated to the complete basis set limit. An analytical site-site potential function with nine sites per methane molecule was fitted to the interaction energies. In addition, a semiempirical correction to the analytical potential function was introduced to take into account the effects of zero-point vibrations. This correction includes adjustments of the dispersion coefficients and of a single-parameter within the fit to the measured values of the second virial coefficient B(T) at room temperature. Quantitative agreement was then obtained with the measured B values over the whole temperature range of the measurements. The calculated B values should definitely be more reliable at very low temperatures (T<150 K) than values extrapolated using the currently recommended equation of state. PMID:18537418
Ab initio intermolecular potential energy surface and thermophysical properties of hydrogen sulfide.
Hellmann, Robert; Bich, Eckard; Vogel, Eckhard; Vesovic, Velisa
2011-08-14
A six-dimensional potential energy hypersurface (PES) for two interacting rigid hydrogen sulfide molecules was determined from high-level quantum-mechanical ab initio computations. A total of 4016 points for 405 different angular orientations of two molecules were calculated utilizing the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory and extrapolating the calculated interaction energies to the complete basis set limit. An analytical site-site potential function with eleven sites per hydrogen sulfide molecule was fitted to the interaction energies. The PES has been validated by computing the second pressure virial coefficient, shear viscosity, thermal conductivity and comparing with the available experimental data. The calculated values of volume viscosity were not used to validate the potential as the low accuracy of the available data precluded such an approach. The second pressure virial coefficient was evaluated by means of the Takahashi and Imada approach, while the transport properties, in the dilute limit, were evaluated by utilizing the classical trajectory method. In general, the agreement with the primary experimental data is within the experimental error for temperatures higher than 300 K. For lower temperatures the lack of reliable data indicates that the values of the second pressure virial coefficient and of the transport properties calculated in this work are currently the most accurate estimates for the thermophysical properties of hydrogen sulfide. PMID:21720616
Communication: A benchmark-quality, full-dimensional ab initio potential energy surface for Ar-HOCO
Conte, Riccardo E-mail: jmbowma@emory.edu; Bowman, Joel M. E-mail: jmbowma@emory.edu; Houston, Paul L.
2014-04-21
A full-dimensional, global ab initio potential energy surface (PES) for the Ar-HOCO system is presented. The PES consists of a previous intramolecular ab initio PES for HOCO [J. Li, C. Xie, J. Ma, Y. Wang, R. Dawes, D. Xie, J. M. Bowman, and H. Guo, J. Phys. Chem. A 116, 5057 (2012)], plus a new permutationally invariant interaction potential based on fitting 12 432 UCCSD(T)-F12a/aVDZ counterpoise-corrected energies. The latter has a total rms fitting error of about 25 cm{sup −1} for fitted interaction energies up to roughly 12 000 cm{sup −1}. Two additional fits are presented. One is a novel very compact permutational invariant representation, which contains terms only involving the Ar-atom distances. The rms fitting error for this fit is 193 cm{sup −1}. The other fit is the widely used pairwise one. The pairwise fit to the entire data set has an rms fitting error of 427 cm{sup −1}. All of these potentials are used in preliminary classical trajectory calculations of energy transfer with a focus on comparisons with the results using the benchmark potential.
NASA Astrophysics Data System (ADS)
Delahaye, Thibault; Rey, Michael; Tyuterev, Vladimir; Nikitin, Andrei V.; Szalay, Peter
2015-06-01
Hydrocarbons such as ethylene (C_2H_4) and methane (CH_4) are of considerable interest for the modeling of planetary atmospheres and other astrophysical applications. Knowledge of rovibrational transitions of hydrocarbons is of primary importance in many fields but remains a formidable challenge for the theory and spectral analysis. Essentially two theoretical approaches for the computation and prediction of spectra exist. The first one is based on empirically-fitted effective spectroscopic models. Several databases aim at collecting the corresponding data but the information about C_2H_4 spectrum present in these databases remains limited, only some spectral ranges around 1000, 3000 and 6000 cm-1 being available. Another way for computing energies, line positions and intensities is based on global variational calculations using ab initio surfaces. Although they do not yet reach the spectroscopic accuracy, they could provide reliable predictions which could be quantitatively accurate with respect to the precision of available observations and as complete as possible. All this thus requires extensive first-principles quantum mechanical calculations essentially based on two necessary ingredients: (i) accurate intramolecular potential energy surface and dipole moment surface components and (ii) efficient computational methods to achieve a good numerical convergence. We report predictions of vibrational and rovibrational energy levels of C_2H_4 using our new ground state potential energy surface obtained from extended ab initio calculations. Additionally we will introduce line positions and line intensities predictions based on a new dipole moment surface for ethylene. These results will be compared with previous works on ethylene and its isotopologues.
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Taylor, Peter R.
1988-01-01
Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F+H2 yields HF+H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces.
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Taylor, Peter R.
1989-01-01
Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F + H2 yields HF + H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces.
NASA Astrophysics Data System (ADS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Taylor, Peter R.
Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F + H2 yields HF + H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces.
NASA Astrophysics Data System (ADS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Taylor, Peter R.
Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F+H2 yields HF+H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces.
Wang, Zhongquan; Feng, Eryin; Yu, Haijun; Zhang, Chunzao; Du, Jianming
2011-01-01
A new high quality three-dimensional potential energy surface for the Ne–CO van der Waals complex is developed using the CCSD(T) method and avqz/avqz+33221 basis set. The ab initio calculation is performed in a total of 1365 configurations with supermolecule method. There is a single global minimum located in a nearly T-shaped geometry. The global minimum energy is −49.4090 cm−1 at Re =6.40a0 and θe =82.5° for V00 . Using the three-dimensional potential energy surface, we have calculated bound rovibrational energy levels up to J = 10 including the Coriolis coupling terms. Compared with the experimental transition frequencies, the theoretical results are in good agreement with the experimental results. PMID:21241112
Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface
NASA Astrophysics Data System (ADS)
McCoy, Anne B.; Huang, Xinchuan; Carter, Stuart; Landeweer, Marc Y.; Bowman, Joel M.
2005-02-01
We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh et al. [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.
Grenier, Romain; To, Quy-Dong; de Lara-Castells, María Pilar; Léonard, Céline
2015-07-01
Global potentials for the interaction between the Ar atom and gold surfaces are investigated and Ar-Au pair potentials suitable for molecular dynamics simulations are derived. Using a periodic plane-wave representation of the electronic wave function, the nonlocal van-der-Waals vdW-DF2 and vdW-OptB86 approaches have been proved to describe better the interaction. These global interaction potentials have been decomposed to produce pair potentials. Then, the pair potentials have been compared with those derived by combining the dispersionless density functional dlDF for the repulsive part with an effective pairwise dispersion interaction. These repulsive potentials have been obtained from the decomposition of the repulsive interaction between the Ar atom and the Au2 and Au4 clusters and the dispersion coefficients have been evaluated by means of ab initio calculations on the Ar+Au2 complex using symmetry adapted perturbation theory. The pair potentials agree very well with those evaluated through periodic vdW-DF2 calculations. For benchmarking purposes, CCSD(T) calculations have also been performed for the ArAu and Ar+Au2 systems using large basis sets and extrapolations to the complete basis set limit. This work highlights that ab initio calculations using very small surface clusters can be used either as an independent cross-check to compare the performance of state-of-the-art vdW-corrected periodic DFT approaches or, directly, to calculate the pair potentials necessary in further molecular dynamics calculations. PMID:26046588
NASA Astrophysics Data System (ADS)
Zhang, Peng; Irle, Stephan; Morokuma, Keiji; Tschumper, Gregory S.
2003-10-01
The mechanism of photodissociation of the vinyl radical, C2H3, starting from the first doublet excited (D1,Ã) state was studied with high-level ab initio methods as well as with ab initio direct dynamics. Geometry optimizations of stationary points and surface crossing seams were performed with complete active space self-consistent field (CASSCF) method, and the energies were re-evaluated with single-point multireference single and double excitation configuration interaction (MRCISD) calculations. Both internal conversion and intersystem crossing channels, which could bring the excited vinyl radical down to the ground state potential energy surface leading to dissociation on the ground state, have been identified within planar Cs, twisted Cs and C2v symmetry. Direct dynamics calculation indicates that the most feasible reaction channel is the direct internal conversion from D1 to the ground state (D0) within planar Cs symmetry, through a minimum of seam of crossing (conical intersection) at an energy of about 80 kcal/mol (with respect to the ground-state equilibrium geometry). The other internal conversions from D1 to D0 through conical intersections within twisted Cs symmetry require energies of about 80 and 75 kcal/mol at the two minima of seam of crossing, respectively, and they are not favored dynamically without initial out-of-plane vibrational excitation. The intersystem crossing channels between D1 and the lowest quartet state (Q1) and D0 and Q1 within twisted Cs and C2v symmetry are not efficient due to the high energy of the minima of seam of crossing as well as the small spin-orbit coupling.
Ab Initio Study of Water Adsorption and Reactivity on the (211) Surface of Anatase TiO2
NASA Astrophysics Data System (ADS)
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.
Dawes, R.; Thompson, D. L.; Guo, Y.; Wagner, A. F.; Minkoff, M.; Chemistry; Univ. of Missouri-Columbia; Oklahoma State Univ.
2007-05-11
A highly accurate and efficient method for molecular global potential energy surface (PES) construction and fitting is demonstrated. An interpolating-moving-least-squares (IMLS)-based method is developed using low-density ab initio Hessian values to compute high-density PES parameters suitable for accurate and efficient PES representation. The method is automated and flexible so that a PES can be optimally generated for classical trajectories, spectroscopy, or other applications. Two important bottlenecks for fitting PESs are addressed. First, high accuracy is obtained using a minimal density of ab initio points, thus overcoming the bottleneck of ab initio point generation faced in applications of modified-Shepard-based methods. Second, high efficiency is also possible (suitable when a huge number of potential energy and gradient evaluations are required during a trajectory calculation). This overcomes the bottleneck in high-order IMLS-based methods, i.e., the high cost/accuracy ratio for potential energy evaluations. The result is a set of hybrid IMLS methods in which high-order IMLS is used with low-density ab initio Hessian data to compute a dense grid of points at which the energy, Hessian, or even high-order IMLS fitting parameters are stored. A series of hybrid methods is then possible as these data can be used for neural network fitting, modified-Shepard interpolation, or approximate IMLS. Results that are indicative of the accuracy, efficiency, and scalability are presented for one-dimensional model potentials as well as for three-dimensional (HCN) and six-dimensional (HOOH) molecular PESs
Collective rotation from ab initio theory
NASA Astrophysics Data System (ADS)
Caprio, M. A.; Maris, P.; Vary, J. P.; Smith, R.
2015-08-01
Through ab initio approaches in nuclear theory, we may now seek to quantitatively understand the wealth of nuclear collective phenomena starting from the underlying internucleon interactions. No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, as evidenced by rotational patterns for excitation energies, electromagnetic moments and electromagnetic transitions. In this review, NCCI calculations of 7-9Be are used to illustrate and explore ab initio rotational structure, and the resulting predictions for rotational band properties are compared with experiment. We highlight the robustness of ab initio rotational predictions across different choices for the internucleon interaction.
NASA Astrophysics Data System (ADS)
Chu, Xingli; Zhang, Yanxing; Li, Shasha; Yang, Zongxian
2014-04-01
The first-principles method based on density functional theory (DFT) is used to investigate the reaction mechanism for the adsorption of H2S on the oxygen-enriched yttria-stabilized zirconia (YSZ + O) (111) surface. It is found that the H2S dissociation processes have low energy barriers (< 0.5 eV) and high exothermicities (2.5 eV), and the dissociative S atoms may result in the poisoning of the YSZ + O surface by forming the SO and the hyposulfite (SO22 -) species with very strong bonds to the surface. In addition, using the ab initio atomistic thermodynamics method, the surface regeneration or de-sulfurization process of a sulfur-poisoned (i.e. sulfur-covered) YSZ + O(111) surface is studied. According to the phase diagram, the adsorbed atomic sulfur can be oxidized to SO2 and removed from the YSZ + O surface by introducing oxidizing reagents, e.g. O2 and H2O.
Improved ab initio calculation of surface second-harmonic generation from Si(111)(1 ×1 ):H
NASA Astrophysics Data System (ADS)
Anderson, Sean M.; Tancogne-Dejean, Nicolas; Mendoza, Bernardo S.; Véniard, Valérie
2016-06-01
We carry out an improved ab initio calculation of surface second-harmonic generation (SSHG) from the Si(111)(1 ×1 ):H surface. This calculation includes three new features in one formulation: (i) the scissors correction, (ii) the contribution of the nonlocal part of the pseudopotentials, and (iii) the inclusion of a cut function to extract the surface response, all within the independent particle approximation. We apply these improvements on the Si(111)(1 ×1 ):H surface and compare with various experimental spectra from several different sources. We also revisit the three-layer model for the SSHG yield and demonstrate that it provides more accurate results over several, more common, two-layer models. We demonstrate the importance of using properly relaxed coordinates for the theoretical calculations. We conclude that this approach to the calculation of the second-harmonic spectra is versatile and accurate within this level of approximation. This well-characterized surface offers an excellent platform for comparison with theory and allows us to offer this study as an efficient benchmark for this type of calculation.
Ab initio nuclear structure theory
NASA Astrophysics Data System (ADS)
Negoita, Gianina Alina
Ab initio no core methods have become major tools for understanding the properties of light nuclei based on realistic nucleon-nucleon (NN) and three-nucleon (NNN) interactions. A brief description is provided for the inter-nucleon interactions that fit two-body scattering and bound state data, as well as NNN interactions. Major new progress, including the goal of applying these interactions to solve for properties of nuclei, is limited by convergence issues. That is, with the goal of obtaining high precision solutions of the nuclear many-body Hamiltonian with no core methods (all nucleons treated on the same footing), one needs to proceed to very large basis spaces to achieve a convergence pattern suitable for extrapolation to the exact result. This thesis investigates (1) the similarity renormalization group (SRG) approach to soften the interaction, while preserving its phase shift properties, and (2) adoption of a realistic basis space using Woods-Saxon (WS) single-particle wavefunctions. Both have their advantages and limitations, discussed here. For (1), SRG was demonstrated by applying it to a realistic NN interaction, JISP16, in a harmonic oscillator (HO) representation. The degree of interaction softening achieved through a regulator parameter is examined. For (2), new results are obtained with the realistic JISP16 NN interaction in ab initio calculations of light nuclei 4He, 6He and 12C, using a WS basis optimized to minimize the ground-state energy within the truncated no core shell model. These are numerically-intensive many-body calculations. Finally, to gain insight into the potential for no core investigations of heavier nuclei, an initial investigation was obtained for the odd mass A = 47 - 49 region nuclei straddling 48Ca. The motivation for selecting these nuclei stems from the aim of preparing for nuclear double beta-decay studies of 48Ca. In these heavier systems, phenomenological additions to the realistic NN interaction determined by previous
Mills, Jeffrey D; Ben-Nun, Michal; Rollin, Kyle; Bromley, Michael W J; Li, Jiabo; Hinde, Robert J; Winstead, Carl L; Sheehy, Jeffrey A; Boatz, Jerry A; Langhoff, Peter W
2016-08-25
Continuing attention has addressed incorportation of the electronically dynamical attributes of biomolecules in the largely static first-generation molecular-mechanical force fields commonly employed in molecular-dynamics simulations. We describe here a universal quantum-mechanical approach to calculations of the electronic energy surfaces of both small molecules and large aggregates on a common basis which can include such electronic attributes, and which also seems well-suited to adaptation in ab initio molecular-dynamics applications. In contrast to the more familiar orbital-product-based methodologies employed in traditional small-molecule computational quantum chemistry, the present approach is based on an "ex-post-facto" method in which Hamiltonian matrices are evaluated prior to wave function antisymmetrization, implemented here in the support of a Hilbert space of orthonormal products of many-electron atomic spectral eigenstates familiar from the van der Waals theory of long-range interactions. The general theory in its various forms incorporates the early semiempirical atoms- and diatomics-in-molecules approaches of Moffitt, Ellison, Tully, Kuntz, and others in a comprehensive mathematical setting, and generalizes the developments of Eisenschitz, London, Claverie, and others addressing electron permutation symmetry adaptation issues, completing these early attempts to treat van der Waals and chemical forces on a common basis. Exact expressions are obtained for molecular Hamiltonian matrices and for associated energy eigenvalues as sums of separate atomic and interaction-energy terms, similar in this respect to the forms of classical force fields. The latter representation is seen to also provide a long-missing general definition of the energies of individual atoms and of their interactions within molecules and matter free from subjective additional constraints. A computer code suite is described for calculations of the many-electron atomic eigenspectra and
Ab-initio study of hydrogen doping and oxygen vacancy at anatase TiO{sub 2} surface
Sotoudeh, M.; Mohammadizadeh, M. R.; Hashemifar, S. J.; Abbasnejad, M.
2014-02-15
Density functional–pseudopotential calculations were performed to study the effects of hydrogen doping and oxygen vacancy, both individually and together, on the electronic structure and stability of (001) surface of TiO{sub 2} in the anatase phase. Based on our calculations, O/Ti termination is the most stable one, and it appears that p-states of deep and surface O atoms and d-orbitals of surface Ti atoms have roles in the valence band and, the conduction band comes from the d-orbitals of deep Ti atoms. Although, no considerable change was seen during H doping, a mid-gap state appeared below the conduction band in the O vacancy configuration. In the framework of ab-initio atomistic thermodynamics, we argue that the anatase TiO{sub 2} prefers a defected O layer termination in the [001] direction. The obtained electronic structures indicate that H doping in the bulk creates the empty mid-gap state below the conduction band and hence decreases the band gap of the system. This phenomenon may explain the enhanced photocatalytic activity of the anatase TiO{sub 2} (001) surface after hydrogenation.
NASA Astrophysics Data System (ADS)
Zimmermann, Bernd; Mavropoulos, Phivos; Long, Nguyen H.; Gerhorst, Christian-Roman; Blügel, Stefan; Mokrousov, Yuriy
2016-04-01
The Fermi surfaces and Elliott-Yafet spin-mixing parameter (EYP) of several elemental metals are studied by ab initio calculations. We focus first on the anisotropy of the EYP as a function of the direction of the spin-quantization axis [B. Zimmermann et al., Phys. Rev. Lett. 109, 236603 (2012), 10.1103/PhysRevLett.109.236603]. We analyze in detail the origin of the gigantic anisotropy in 5 d hcp metals as compared to 5 d cubic metals by band structure calculations and discuss the stability of our results against an applied magnetic field. We further present calculations of light (4 d and 3 d ) hcp crystals, where we find a huge increase of the EYP anisotropy, reaching colossal values as large as 6000 % in hcp Ti. We attribute these findings to the reduced strength of spin-orbit coupling, which promotes the anisotropic spin-flip hot loops at the Fermi surface. In order to conduct these investigations, we developed an adapted tetrahedron-based method for the precise calculation of Fermi surfaces of complicated shape and accurate Fermi-surface integrals within the full-potential relativistic Korringa-Kohn-Rostoker Green function method.
Aruguete, Deborah A.; Marcus, Matthew A.; Li, Liang-shi; Williamson, Andrew; Fakra, Sirine; Gygi, Francois; Galli, Giulia; Alivisatos, A. Paul
2006-01-27
We report orientation-specific, surface-sensitive structural characterization of colloidal CdSe nanorods with extended X-ray absorption fine structure spectroscopy and ab-initio density functional theory calculations. Our measurements of crystallographically-aligned CdSe nanorods show that they have reconstructed Cd-rich surfaces. They exhibit orientation-dependent changes in interatomic distances which are qualitatively reproduced by our calculations. These calculations reveal that the measured interatomic distance anisotropy originates from the nanorod surface.
Ab initio phonon limited transport
NASA Astrophysics Data System (ADS)
Verstraete, Matthieu
We revisit the thermoelectric (TE) transport properties of two champion materials, PbTe and SnSe, using fully first principles methods. In both cases the performance of the material is due to subtle combinations of structural effects, scattering, and phase space reduction. In PbTe anharmonic effects are completely opposite to the predicted quasiharmonic evolution of phonon frequencies and to frequently (and incorrectly) cited extrapolations of experiments. This stabilizes the material at high T, but also tends to enhance its thermal conductivity, in a non linear manner, above 600 Kelvin. This explains why PbTe is in practice limited to room temperature applications. SnSe has recently been shown to be the most efficient TE material in bulk form. This is mainly due to a strongly enhanced carrier concentration and electrical conductivity, after going through a phase transition from 600 to 800 K. We calculate the transport coefficients as well as the defect concentrations ab initio, showing excellent agreement with experiment, and elucidating the origin of the double phase transition as well as the new charge carriers. AH Romero, EKU Gross, MJ Verstraete, and O Hellman PRB 91, 214310 (2015) O. Hellman, IA Abrikosov, and SI Simak, PRB 84 180301 (2011)
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-01
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined. PMID:25106565
MoS{sub 2} on an amorphous HfO{sub 2} surface: An ab initio investigation
Scopel, W. L.; Miwa, R. H. Schmidt, T. M.; Venezuela, P.
2015-05-21
The energetic stability, electronic and structural properties of MoS{sub 2} adsorbed on an amorphous a-HfO{sub 2} surface (MoS{sub 2}/HfO{sub 2}) are examined through ab initio theoretical investigations. Our total energy results indicate that the formation of MoS{sub 2}/HfO{sub 2} is an exothermic process with an adsorption energy of 34 meV/Å{sup 2}, which means that it is more stable than similar systems like graphene/HfO{sub 2} and MoS{sub 2}/SiO{sub 2}. There are no chemical bonds at the MoS{sub 2}-HfO{sub 2} interface. Upon formation of MoS{sub 2}/HfO{sub 2}, the electronic charge distribution is mostly localized at the interface region with no net charge transfer between the adsorbed MoS{sub 2} sheet and –HfO{sub 2} surface. However, the MoS{sub 2} sheet becomes n-type doped when there are oxygen vacancies in the HfO{sub 2} surface. Further investigation of the electronic distribution reveals that there are no electron- and hole-rich regions (electron-hole puddles) on the MoS{sub 2} sheet, which makes this system promising for use in high-speed nanoelectronic devices.
NASA Astrophysics Data System (ADS)
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-01
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined.
Jackson, Bret; Nattino, Francesco; Kroes, Geert-Jan
2014-08-07
The dissociative chemisorption of methane on metal surfaces is of great practical and fundamental importance. Not only is it the rate-limiting step in the steam reforming of natural gas, the reaction exhibits interesting mode-selective behavior and a strong dependence on the temperature of the metal. We present a quantum model for this reaction on Ni(100) and Ni(111) surfaces based on the reaction path Hamiltonian. The dissociative sticking probabilities computed using this model agree well with available experimental data with regard to variation with incident energy, substrate temperature, and the vibrational state of the incident molecule. We significantly expand the vibrational basis set relative to earlier studies, which allows reaction probabilities to be calculated for doubly excited initial vibrational states, though it does not lead to appreciable changes in the reaction probabilities for singly excited initial states. Sudden models used to treat the center of mass motion parallel to the surface are compared with results from ab initio molecular dynamics and found to be reasonable. Similar comparisons for molecular rotation suggest that our rotationally adiabatic model is incorrect, and that sudden behavior is closer to reality. Such a model is proposed and tested. A model for predicting mode-selective behavior is tested, with mixed results, though we find it is consistent with experimental studies of normal vs. total (kinetic) energy scaling. Models for energy transfer into lattice vibrations are also examined.
Ab initio Potential-Energy Surfaces and Electron-Spin-Exchange Cross Sections for H-O2 Interactions
NASA Technical Reports Server (NTRS)
Stallcop, James R.; Partridge, Harry; Levin, Eugene
1996-01-01
Accurate quartet- and doublet-state potential-energy surfaces for the interaction of a hydrogen atom and an oxygen molecule in their ground states have been determined from an ab initio calculation using large-basis sets and the internally contracted multireference configuration interaction method. These potential surfaces have been used to calculate the H-O2 electron-spin-exchange cross section; the square root of the cross section (in a(sub 0)), not taking into account inelastic effects, can be obtained approximately from the expressions 2.390E(sup -1/6) and 5.266-0.708 log10(E) at low and high collision energies E (in E(sub h)), respectively. These functional forms, as well as the oscillatory structure of the cross section found at high energies, are expected from the nature of the interaction energy. The mean cross section (the cross section averaged over a Maxwellian velocity distribution) agrees reasonably well with the results of measurements.
Zabidi, Noriza Ahmad; Kassim, Hasan Abu; Shrivastava, Keshav N.
2008-05-20
Polonium is the only element with a simple cubic (sc) crystal structure. Atoms in solid polonium sit at the corners of a simple cubic unit cell and no where else. Polonium has a valence electron configuration 6s{sup 2}6p{sup 4} (Z = 84). The low temperature {alpha}-phase transforms into the rhombohedral (trigonal) {beta} structure at {approx}348 K. The sc {alpha}-Po unit cell constant is a = 3.345 A. The beta form of polonium ({beta}-Po) has the lattice parameters, a{sub R} = 3.359 A and a rhombohedral angle 98 deg. 13'. We have performed an ab initio electronic structure calculation by using the density functional theory. We have performed the calculation with and without spin-orbit (SO) coupling by using both the LDA and the GGA for the exchange-correlations. The k-points in a simple cubic BZ are determined by R (0.5, 0.5, 0.5), {gamma} (0, 0, 0), X (0.5, 0, 0), M (0.5, 0.5, 0) and {gamma} (0, 0, 0). Other directions of k-points are {gamma} (0, 0, 0), X (0.5, 0, 0), R (0.5, 0.5, 0.5) and {gamma} (0, 0, 0). The SO splittings of p states at the {gamma} point in the GGA+SO scheme for {alpha}-Po are 0.04 eV and 0.02 eV while for the {beta}-Po these are 0.03 eV and 0.97 eV. We have also calculated the vibrational spectra for the unit cells in both the structures. We find that exchanging of a Po atom by Pb atom produces several more bands and destabilizes the {beta} phase.
Alexandrov, Vitali Y.; Rosso, Kevin M.
2015-01-01
Goethite (α-FeOOH) surfaces represent one of the most ubiquitous redox-active interfaces in the environment, playing an important role in biogeochemical metal cycling and contaminant residence in the subsurface. Fe(II)-catalyzed recrystallization of goethite is a fundamental process in this context, but the proposed Fe(II)aq-Fe(III)goethite electron and iron atom exchange mechanism of recrystallization remains poorly understood at the atomic level. We examine the adsorption of aqueous Fe(II) and subsequent interfacial electron transfer (ET) between adsorbed Fe(II) and structural Fe(III) at the (110) and (021) goethite surfaces using density functional theory calculations including Hubbard U corrections (DFT+U) aided by ab initio molecular dynamics simulations. We investigate various surface sites for the adsorption of Fe2+(H2O)6 in different coordination environments. Calculated energies for adsorbed complexes at both surfaces favor monodentate complexes with reduced 4- and 5-fold coordination over higher-dentate structures and 6- fold coordination. The hydrolysis of H2O ligands is observed for some pre-ET adsorbed Fe(II) configurations. ET from the adsorbed Fe(II) into the goethite lattice is calculated to be energetically uphill always, but simultaneous proton transfer from H2O ligands of the adsorbed complexes to the surface oxygen species stabilizes post-ET states. We find that surface defects such as oxygen vacancies near the adsorption site also can stabilize post-ET states, enabling the Fe(II)aq-Fe(III)goethite interfacial electron transfer reaction implied from experiments to proceed.
Marquardt, Roberto; Sagui, Kenneth; Zheng, Jingjing; Thiel, Walter; Luckhaus, David; Yurchenko, Sergey; Mariotti, Fabio; Quack, Martin
2013-08-15
The analytical, full-dimensional, and global representation of the potential energy surface of NH(3) in the lowest adiabatic electronic state developed previously (Marquardt, R.; et al. J. Phys. Chem. B 2005, 109, 8439–8451) is improved by adjustment of parameters to an enlarged set of electronic energies from ab initio calculations using the coupled cluster method with single and double substitutions and a perturbative treatment of connected triple excitations (CCSD(T)) and the method of multireference configuration interaction (MRCI). CCSD(T) data were obtained from an extrapolation of aug-cc-pVXZ results to the basis set limit (CBS), as described in a previous work (Yurchenko, S.N.; et al. J. Chem. Phys 2005, 123, 134308); they cover the region around the NH3 equilibrium structures up to 20,000 hc cm(–1). MRCI energies were computed using the aug-cc-pVQZ basis to describe both low lying singlet dissociation channels. Adjustment was performed simultaneously to energies obtained from the different ab initio methods using a merging strategy that includes 10,000 geometries at the CCSD(T) level and 500 geometries at the MRCI level. Characteristic features of this improved representation are NH3 equilibrium geometry r(eq)(NH(3)) ≈ 101.28 pm, α(eq)(NH(3)) ≈ 107.03°, the inversion barrier at r(inv)(NH(3)) ≈ 99.88 pm and 1774 hc cm(–1) above the NH(3) minimum, and dissociation channel energies 41,051 hc cm(–1) (for NH(3) → ((2)B(2))NH(2) + ((2)S(1/2))H) and 38,450 hc cm(–1) (for NH(3) → ((3)Σ(–))NH +((1)Σ(g)(+))H(2)); the average agreement between calculated and experimental vibrational line positions is 11 cm(–1) for (14)N(1)H(3) in the spectral region up to 5000 cm(–1). A survey of our current knowledge on the vibrational spectroscopy of ammonia and its isotopomers is also given. PMID:23688044
NASA Astrophysics Data System (ADS)
Churakov, S. V.; Parrinello, M.
2003-04-01
The carbonation of brucite (Mg(OH)2) has been considered as a potential technology for cleaning industrial carbon dioxide waste. The kinetics of the reaction Mg(OH)2 + CO2 -> MgCO3 + H2O have been studied experimentally at 573°C by Bearat at al. [1]. Their experiments suggest that the carbonation of magnesium hydroxide proceeds by the reaction Mg(OH)2 -> MgO + H2O followed by the adsorption of CO2 molecules on the dehydrated brucite surface. Due to the large difference in volumes between Mg(OH)2 and MgO, dehydration causes the formation of dislocations and cracks, allowing water molecules to leave the brucite surface and facilitating the advance of the carbonation front in the bulk solid. The detailed mechanism of this process is however unknown. We used the Car-Parrinello ab initio molecular dynamics method to study the structure and dynamics of the (0001), (1-100) and (11-20) surfaces of brucit and calculated the enthalpy and activation barrier of H2O nucleation and dehydration on different surfaces. The results obtained are in agreement with previous studies of brucite dehydration by Masini and Bernasconi [2]. The reactive Car-Parrinello molecular dynamics method [3] has been applied to investigate the detailed mechanism of the dehydration-carbonation reaction at the (1-100) interface of brucite with the gas phase. Based on the results of our MD simulations and the calculated enthalpy of CO2 adsorption on the dehydrated brucite surfaces we propose a mechanism for the dehydration/carbonation reaction. [1] Bearat H, McKelvy MJ, Chizmeshya AVG, Sharma R, Carpenter RW (2002) J. Amer. Ceram. Soc. 85(4):742 [2] Masini P and Bernasconi M (2001) J. Phys. Cond. Mat. 13: 1-12 [3] Iannuzzi M, Laio A and Parrinello M (2003) Phys. Rev. Lett. (submitted)
NASA Astrophysics Data System (ADS)
Pai, Sung Jin; Bae, Young Chan
2014-08-01
A six-dimensional intermolecular potential energy surface for a rigid methane (CH4) and carbon dioxide (CO2) dimer was developed from the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory. A total of 466 grid points distributed to 46 orientations were calculated from the complete basis set limit extrapolation based on up to aug-cc-pVQZ basis set. A modified site-site pair potential function was proposed for rapid representation of the high level ab initio calculations. A nonadditive three-body interaction was represented by the Axilrod-Teller-Muto expression for mixtures with the polarizability and the London dispersion constant of each molecule. Second to fourth virial coefficients of CH4 and CO2 mixtures were calculated using both the Mayer sampling Monte Carlo method and the present potential functions. The virial equation of state derived from these coefficients was used to predict the pVT values and showed good agreement with experimental data below 200 bar at 300 K. The vapor-liquid coexistence curves of pure CH4, CO2 and their mixtures were presented with the aid of Gibbs ensemble Monte Carlo simulations. The predicted tie lines agreed with the experimental data within the uncertainties up to near the critical point.
Ab initio potential energy surface and vibration-rotation energy levels of silicon dicarbide, SiC2.
Koput, Jacek
2016-10-01
The accurate ground-state potential energy surface of silicon dicarbide, SiC2 , has been determined from ab initio calculations using the coupled-cluster approach. Results obtained with the conventional and explicitly correlated coupled-cluster methods were compared. The core-electron correlation, higher-order valence-electron correlation, and scalar relativistic effects were taken into account. The potential energy barrier to the linear SiCC configuration was predicted to be 1782 cm(-1) . The vibration-rotation energy levels of the SiC2 , (29) SiC2 , (30) SiC2 , and SiC(13) C isotopologues were calculated using a variational method. The experimental vibration-rotation energy levels of the main isotopologue were reproduced to high accuracy. In particular, the experimental energy levels of the highly anharmonic vibrational ν3 mode of SiC2 were reproduced to within 6.7 cm(-1) , up to as high as the v3 = 16 state. PMID:27481562
Marinakis, Sarantos; Dean, Indigo Lily; Kłos, Jacek; Lique, François
2015-09-01
We present a new set of potential energy surfaces (PESs) for the CH(X(2)Π)-He van der Waals system. Ab initio calculations of the CH-He PES were carried out using the open-shell single- and double-excitation coupled cluster approach with non-iterative perturbational treatment of triple excitations [RCCSD(T)]. The augmented correlation-consistent polarized valence quadruple-zeta (aug-cc-pVQZ) basis set was employed augmented by mid-bond functions. Integral cross sections for the rotational excitation in CH-He collisions were calculated using the new PES and compared with available experimental results. The newly constructed PES reproduces the available experimental results for CH(X(2)Π, v = 0)-He collisions better than any previously available PES. Differential cross sections (DCS) are presented for the first time for this system and discussed within the context of rotational rainbows. Finally, our work provides the first rate thermal coefficients for this system that are crucially needed for astrochemical modelling and future anticipated experiments in CH(X(2)Π)-He collisions. PMID:26220835
Four-faceted nanowires generated from densely-packed TiO2 rutile surfaces: Ab initio calculations
NASA Astrophysics Data System (ADS)
Evarestov, R. A.; Zhukovskii, Yu. F.
2013-02-01
Two-dimensional (2D) slabs and monoperiodic (1D) nanowires orthogonal to the slab surface of rutile-based TiO2 structure terminated by densely-packed surfaces and facets, respectively, have been simulated in the current study. The procedure of structural generation of nanowires (NWs) from titania slabs (2D → 1D) is described. We have simulated: (i) (110), (100), (101) and (001) slabs of different thicknesses as well as (ii) [001]- and [110]-oriented nanowires of different diameters terminated by either four types of related {110} facets or alternating {11¯0} and {001} facets, respectively. Nanowires have been described using both the Ti atom-centered rotation axes as well as the hollow site-centered axes passing through the interstitial sites between the Ti and O atoms closest to the axes. For simulations on TiO2 slabs and NWs, we have performed large-scale ab initio Density Functional Theory (DFT) and hybrid DFT-Hartree Fock (DFT-HF) calculations with the total geometry optimization within the Generalized Gradient Approximation (GGA) in the form of the Perdew-Becke-Ernzenhof exchange-correlation functionals (PBE and PBE0, respectively), using the formalism of linear combination of localized atomic functions (LCAO) implemented in CRYSTAL09 code. Both structural and electronic properties of enumerated rutile-based titania slabs and nanowires have been calculated. According to the results of our surface energy calculations, the most stable rutile-based titania slab is terminated by (110) surfaces whereas the energetically favorable [001]-oriented NWs are also terminated by {110} facets only, thus confirming results of previous studies.
Dawes, R.; Wagner, A. F.; Thompson, D. L.; Chemical Sciences and Engineering Division; Univ. of Missouri at Columbia
2009-04-23
We report here calculated J = 0 vibrational frequencies for {sup 1}CH{sub 2} and HCN with root-mean-square error relative to available measurements of 2.0 cm{sup -1} and 3.2 cm{sup -1}, respectively. These results are obtained with DVR calculations with a dense grid on ab initio potential energy surfaces (PESs). The ab initio electronic structure calculations employed are Davidson-corrected MRCI calculations with double-, triple-, and quadruple-{zeta} basis sets extrapolated to the complete basis set (CBS) limit. In the {sup 1}CH{sub 2} case, Full CI tests of the Davidson correction at small basis set levels lead to a scaling of the correction with the bend angle that can be profitably applied at the CBS limit. Core-valence corrections are added derived from CCSD(T) calculations with and without frozen cores. Relativistic and non-Born-Oppenheimer corrections are available for HCN and were applied. CBS limit CCSD(T) and CASPT2 calculations with the same basis sets were also tried for HCN. The CCSD(T) results are noticeably less accurate than the MRCI results while the CASPT2 results are much poorer. The PESs were generated automatically using the local interpolative moving least-squares method (L-IMLS). A general triatomic code is described where the L-IMLS method is interfaced with several common electronic structure packages. All PESs were computed with this code running in parallel on eight processors. The L-IMLS method provides global and local fitting error measures important in automatically growing the PES from initial ab initio seed points. The reliability of this approach was tested for {sup 1}CH{sub 2} by comparing DVR-calculated vibrational levels on an L-IMLS ab initio surface with levels generated by an explicit ab initio calculation at each DVR grid point. For all levels ({approx}200) below 20000 cm{sup -1}, the mean unsigned difference between the levels of these two calculations was 0.1 cm{sup -1}, consistent with the L-IMLS estimated mean unsigned
NASA Astrophysics Data System (ADS)
Kalugina, Yulia N.; Lokshtanov, Sergei E.; Cherepanov, Victor N.; Vigasin, Andrey A.
2016-02-01
We present new three-dimensional potential energy surface (PES) and dipole moment surfaces (DMSs) for the CH4-Ar van der Waals system. Ab initio calculations of the PES and DMS were carried out using the closed-shell single- and double-excitation coupled cluster approach with non-iterative perturbative treatment of triple excitations. The augmented correlation-consistent aug-cc-pVXZ (X = D,T,Q) basis sets were employed, and the energies obtained were then extrapolated to the complete basis set limit. The dipole moment surface was obtained using aug-cc-pVTZ basis set augmented with mid-bond functions for better description of exchange interactions. The second mixed virial coefficient was calculated and compared to available experimental data. The equilibrium constant for true dimer formation was calculated using classical partition function based on the knowledge of ab initio PES. Temperature variations of the zeroth spectral moment of the rototranslational collision-induced band as well as its true dimer constituent were traced with the use of the Boltzmann-weighted squared induced dipole properly integrated over respective phase space domains. Height profiles for N2-N2, N2-H2, CH4-N2, (CH4)2, and CH4-Ar true bound dimers in Titan's atmosphere were calculated with the use of reliable ab initio PESs.
Peterson, Kirk A
2000-09-15
A global, analytical potential energy surface for the ground electronic state of HOBr has been determined using highly correlated multireference configuration interaction wave functions and explicit basis set extrapolations of large correlation consistent basis sets. The ab initio data have been fit to an analytical functional form that accurately includes both the HOBr and HBrO minima, as well as all dissociation asymptotes. Small adjustments to this surface are made based on the limited experimental data available and by indirectly taking into account the effects of spin-orbit coupling on the OH+Br dissociation channel. Vibrational energy levels are calculated variationally for both HOBr and HBrO up to the OH+Br dissociation limit using a truncation/recoupling method. The HOBr isomer is calculated to contain 708 bound vibrational energy levels, while the HBrO minimum lies above the OH+Br dissociation limit but is calculated to have 74 ''quasibound,'' localized eigenstates. Infrared intensities for all of these vibrational transitions are also calculated using MRCI dipole moment functions. The assignment of the HOBr states is complicated by strong stretch-bend resonances even at relatively low energies. In contrast to the HOCl case, these state mixings made it particularly difficult to assign the relatively intense OH overtone bands above v{sub 1}=2. The vibrational density of states of HOBr at the OH+Br dissociation limit is determined to be 0.16 states/cm-1. Comparisons to recent work on HOCl using similar methods are made throughout. (c) 2000 American Institute of Physics.
An Accurate Global Ab Initio Potential Energy Surface for the X(1)A' Electronic State of HOBr
Peterson, Kirk A.
1999-12-01
A global, analytical potential energy surface for the ground electronic state of HOBr has been determined using highly correlated multireference configuration interaction wave functions and explicit basis set extrapolations of large correlation consistent basis sets. The ab initio data have been fit to an analytical functional form that accurately includes both the HOBr and HBrO minima, as well as all dissociation asymptotes. Small adjustments to this surface are made based on the limited experimental data available and by indirectly taking into account the effects of spin-orbit coupling on the OH+Br dissociation channel. Vibrational energy levels are calculated variationally for both HOBr and HBrO up to the OH+Br dissociation limit using a truncation-recoupling method. The HOBr isomer is calculated to contain 708 bound vibrational energy levels, while the HBrO minimum lies above the OH+Br dissociation limit but is calculated to have 74 quasibound, localized eigenstates. Infrared intensities for all of these vibrational transitions are also calculated using MRCI dipole moment functions. The assignment of the HOBr states is complicated by strong stretch-bend resonances even at relatively low energies. In contrast to the HOCl case, these state mixings made it particularly difficult to assign the relatively intense OH overtone bands above v1=2. The vibrational density of states of HOBr at the OH+Br dissociation limit is determined to be 0.16 states/cm-1. Comparisons to recent work on HOCl using similar methods are made throughout.
Li, Jun E-mail: zhangdh@dicp.ac.cn; Chen, Jun; Zhao, Zhiqiang; Zhang, Dong H. E-mail: zhangdh@dicp.ac.cn; Xie, Daiqian; Guo, Hua
2015-05-28
We report a permutationally invariant global potential energy surface (PES) for the H + CH{sub 4} system based on ∼63 000 data points calculated at a high ab initio level (UCCSD(T)-F12a/AVTZ) using the recently proposed permutation invariant polynomial-neural network method. The small fitting error (5.1 meV) indicates a faithful representation of the ab initio points over a large configuration space. The rate coefficients calculated on the PES using tunneling corrected transition-state theory and quasi-classical trajectory are found to agree well with the available experimental and previous quantum dynamical results. The calculated total reaction probabilities (J{sub tot} = 0) including the abstraction and exchange channels using the new potential by a reduced dimensional quantum dynamic method are essentially the same as those on the Xu-Chen-Zhang PES [Chin. J. Chem. Phys. 27, 373 (2014)].
Ab Initio Simulations of Silica Alpha Quartz Surfaces and their Interaction with Water
NASA Astrophysics Data System (ADS)
Chen, Yun-Wen; Cheng, Hai-Ping
2008-03-01
Two types of dry silica alpha quartz (0001) surfaces have been investigated through density functional calculations. One is the dense surface proposed by Rignanese et al.[1] that has 3-fold flower like six-member rings on the top and the three-member rings underneath. The other one is a newfound surface that has zigzag-shape six-member rings on the top and the three-member rings underneath. It is found that the new one is energetically more stable than the dense surface. The interactions of the surfaces with water molecules have also been concerned. Both of the surfaces have the similar hydrophilic properties. We report our results from MD simulations at a temperature of 300K. [1]. G.-M. Rignanese, Alessandro De Vita, J.-C. Charlier, Roberto Car, Phys. Rev. B 61, 13250 (2000)
SurfKin: an ab initio kinetic code for modeling surface reactions.
Le, Thong Nguyen-Minh; Liu, Bin; Huynh, Lam K
2014-10-01
In this article, we describe a C/C++ program called SurfKin (Surface Kinetics) to construct microkinetic mechanisms for modeling gas-surface reactions. Thermodynamic properties of reaction species are estimated based on density functional theory calculations and statistical mechanics. Rate constants for elementary steps (including adsorption, desorption, and chemical reactions on surfaces) are calculated using the classical collision theory and transition state theory. Methane decomposition and water-gas shift reaction on Ni(111) surface were chosen as test cases to validate the code implementations. The good agreement with literature data suggests this is a powerful tool to facilitate the analysis of complex reactions on surfaces, and thus it helps to effectively construct detailed microkinetic mechanisms for such surface reactions. SurfKin also opens a possibility for designing nanoscale model catalysts. PMID:25111729
Sun, J.; Pohl, K.; Mikkelsen, A.; Fuglsang Jensen, M.; Hofmann, Ph.; Koroteev, Y. M.; Bihlmayer, G.; Chulkov, E. V.
2006-12-15
The surface structure of Bi(110) has been investigated by low-energy electron diffraction intensity analysis and by first-principles calculations. Diffraction patterns at a sample temperature of 110 K and normal incidence reveal a bulk truncated (1x1) surface without indication of any structural reconstruction despite the presence of dangling bonds on the surface layer. Good agreement is obtained between the calculated and measured diffraction intensities for this surface containing only one mirror-plane symmetry element and a buckled bilayer structure. No significant interlayer spacing relaxations are found. The Debye temperature for the surface layer is found to be lower than in the bulk, which is indicative of larger atomic vibrational amplitudes at the surface. Meanwhile, the second layer shows a Debye temperature close to the bulk value. The experimental results for the relaxations agree well with those of our first-principles calculation.
Chemisorption of a Molecular Oxygen On the UN(0 0 1) Surface: Ab Initio Calculations
Zhukovskii, Yuri F.; Bocharov, Dmitry; Kotomin, Eugene Alexej
2009-09-15
The results of DFT GGA calculations on oxygen molecules adsorbed upon the (0 0 1) surface of uranium mononitride (UN) are presented and discussed. We demonstrate that O2 molecules oriented parallel to the substrate can dissociate either (i) spontaneously when the molecular center lies above the surface hollow site or atop N ion, (ii) with the activation barrier when a molecule sits atop the surface U ion. This explains fast UN oxidation in air.
Ab Initio Studies of Anatase TiO2 (101) Surface-supported Au8 Clusters.
Mikolajczyk, A; Pinto, H P; Gajewicz, A; Puzyn, T; Leszczynski, J
2015-01-01
Supported transition metals on TiO2 surfaces have shown exceptional catalytic properties in many important process such as CO oxidation, selective propane oxidation, hydrogenation, water adsorption and other catalytic and photocatalytic oxidation reaction at low-temperature. Among the three polymorphs of TiO2, the anatase crystal is the more photoactive. The anatase (101) surface attracts more attention since it has lower surface energy relative to (001) and (100) surfaces and it is observed to adsorb small molecules on its surface. Using density-functional theory (DFT) with on-site Coulomb interactions corrections, we have computed the structural and electronic properties of selected Au8 clusters interacting with clean and reduced anatase TiO2(101) surfaces. The computed adsorption energies are suggesting that the considered Au8 clusters are only physisorbed onto pristine TiO2(101) surface. Oxygen vacancies are found to enhance the absorption of Au8 on the Ti2(101) surface. Accurate simulations required spin polarized DFT since the ground state of Au8 interacting with defective TiO2(101) shows magnetic solutions. The results show that Au8 clusters are chemically bonded to the surface around the locality of the oxygen vacancy. The surface oxygen vacancy is found to be energetically more favourable than sub-surface oxygen vacancy configuration. These vacancy sites may act as nucleation sites for small Au clusters or Au atoms. Finally, the computed electronic structure of all the Au8/TiO2(101) configurations considered in this work are analysed in the light of available experimental data. PMID:25961526
NASA Astrophysics Data System (ADS)
Wang, Yimin; Bowman, Joel M.; Kamarchik, Eugene
2016-03-01
We report full-dimensional, ab initio-based potentials and dipole moment surfaces for NaCl, NaF, Na+H2O, F-H2O, and Cl-H2O. The NaCl and NaF potentials are diabatic ones that dissociate to ions. These are obtained using spline fits to CCSD(T)/aug-cc-pV5Z energies. In addition, non-linear least square fits using the Born-Mayer-Huggins potential are presented, providing accurate parameters based strictly on the current ab initio energies. The long-range behavior of the NaCl and NaF potentials is shown to go, as expected, accurately to the point-charge Coulomb interaction. The three ion-H2O potentials are permutationally invariant fits to roughly 20 000 coupled cluster CCSD(T) energies (awCVTZ basis for Na+ and aVTZ basis for Cl- and F-), over a large range of distances and H2O intramolecular configurations. These potentials are switched accurately in the long range to the analytical ion-dipole interactions, to improve computational efficiency. Dipole moment surfaces are fits to MP2 data; for the ion-ion cases, these are well described in the intermediate- and long-range by the simple point-charge expression. The performance of these new fits is examined by direct comparison to additional ab initio energies and dipole moments along various cuts. Equilibrium structures, harmonic frequencies, and electronic dissociation energies are also reported and compared to direct ab initio results. These indicate the high fidelity of the new PESs.
Wang, Yimin; Bowman, Joel M; Kamarchik, Eugene
2016-03-21
We report full-dimensional, ab initio-based potentials and dipole moment surfaces for NaCl, NaF, Na(+)H2O, F(-)H2O, and Cl(-)H2O. The NaCl and NaF potentials are diabatic ones that dissociate to ions. These are obtained using spline fits to CCSD(T)/aug-cc-pV5Z energies. In addition, non-linear least square fits using the Born-Mayer-Huggins potential are presented, providing accurate parameters based strictly on the current ab initio energies. The long-range behavior of the NaCl and NaF potentials is shown to go, as expected, accurately to the point-charge Coulomb interaction. The three ion-H2O potentials are permutationally invariant fits to roughly 20,000 coupled cluster CCSD(T) energies (awCVTZ basis for Na(+) and aVTZ basis for Cl(-) and F(-)), over a large range of distances and H2O intramolecular configurations. These potentials are switched accurately in the long range to the analytical ion-dipole interactions, to improve computational efficiency. Dipole moment surfaces are fits to MP2 data; for the ion-ion cases, these are well described in the intermediate- and long-range by the simple point-charge expression. The performance of these new fits is examined by direct comparison to additional ab initio energies and dipole moments along various cuts. Equilibrium structures, harmonic frequencies, and electronic dissociation energies are also reported and compared to direct ab initio results. These indicate the high fidelity of the new PESs. PMID:27004880
Rong, Xi; Kolpak, Alexie M
2015-05-01
The design of efficient, stable, and inexpensive catalysts for oxygen evolution and reduction is crucial for the development of electrochemical energy conversion devices such as fuel cells and metal-air batteries. Currently, such design is limited by challenges in atomic-scale experimental characterization and computational modeling of solid-liquid interfaces. Here, we begin to address these issues by developing a general-, first-principles-, and electrochemical-principles-based framework for prediction of catalyst surface structure, stoichiometry, and stability as a function of pH, electrode potential, and aqueous cation concentration. We demonstrate the approach by determining the surface phase diagram of LaMnO3, which has been studied for oxygen evolution and reduction and computing the reaction overpotentials on the relevant surface phases. Our results illustrate the critical role of solvated cation species in governing the catalyst surface structure and stoichiometry, and thereby catalytic activity, in aqueous solution. PMID:26263350
Ab initio study of molecular and atomic oxygen on GeTe(111) surfaces
Deringer, Volker L.; Dronskowski, Richard
2014-11-07
Oxidation of the phase-change material germanium telluride (GeTe) is an atomic-scale process of fundamental importance, as it is detrimental to the stability of GeTe-based data-storage devices. Here, we present comprehensive density-functional theory simulations of molecular and atomic oxygen in contact with GeTe(111) surfaces. Molecular O{sub 2} is predicted to readily adsorb on the Ge-terminated (111) surface; the pristine Te-terminated counterpart, by contrast, appears quite inert. The coverage-dependent adsorption of O atoms is then investigated, and based on these data, a surface phase diagram for GeTe(111)/O is constructed. These results afford a detailed, atom-resolved picture of the initial surface oxidation of GeTe, and they harmonize well with a previous X-ray photoelectron spectroscopy study on this very topic.
Comparison of the bulk and surface properties of ceria and zirconia by ab initio investigations
Gennard, S.; Cora, F.; Catlow, C.R.A.
1999-11-18
Quantum mechanical (QM) calculations were made at a periodic Hartree-Fock (HF) level, on the bulk and surface properties of cubic CeO{sub 2} and ZrO{sub 2}. The authors have investigated the M-O bonding features, and established the high degree of ionicity of both materials, which is greater for CeO{sub 2} than ZrO{sub 2}. The calculated values for the C{sub 11}, C{sub 12}, and C{sub 44} elastic constants, are in close agreement with experiment; an extended oxygen basis set, containing d-orbital polarization functions, is essential to model accurately the symmetry lowering during the C{sub 44} distortion. In the surface studies, they have calculated the surface energies of the {l{underscore}brace}011{r{underscore}brace} and {l{underscore}brace}111{r{underscore}brace} faces of both ceria and zirconia. QM results are compared with interatomic potential-based (IP) methods to assess the accuracy of the latter. IP methods were found to provide a correct estimate of the surface relaxations and the correct order of stability of the two faces examined, with the energy difference between the {l{underscore}brace}011{r{underscore}brace} and the {l{underscore}brace}111{r{underscore}brace} surfaces being approximately 1 J/m{sup 2}, as in the QM study. However, IP methods do not always discriminate adequately between the properties of the two materials. Finally, geometric and electronic relaxations in the {l{underscore}brace}111{r{underscore}brace} surface are confined to the outermost oxygen ions, while in the {l{underscore}brace}011{r{underscore}brace} slabs they are more important and extend to the subsurface layers in a columnar way. The unsaturation of the surface ions in the {l{underscore}brace}011{r{underscore}brace} face may have important implications for the catalytic activity of the materials.
Nie, JL; Xiao, Haiyan J.; Gao, Fei; Zu, Xiaotao T.
2009-05-12
First-principles calculations based on density functional theory with the generalized gradient approximation have been performed to study the aluminum (Al) adsorption on the (001) surface of α-uranium (α-U). The geometric, electronic and magnetic properties have been investigated at coverages of 0.25 and 0.5 monolayer. The results show that the quasi-trigonal sites are preferred at both coverages. The bonding of Al with U is found to be metallic, which mainly arises from the mixing of Al 3sp and U 5f states. A ferromagnetic phase is determined for the bare α-U(001) surface, while the adsorption of Al on the surface significantly perturbs the spin arrangement pattern and reduces the local magnetic moment, leading to a ferrimagnetic phase on the α-U(001) surface at the coverage of 0.5 monolayer. However, the Al overlayer is paramagnetic. Generally, the spin polarization has negligible effects on the geometric and electronic properties of Al atoms on the α-U(001) surface.
Ab initio study of atomic disorder on As-rich GaAs(111)A surface
NASA Astrophysics Data System (ADS)
Romanyuk, O.; Mutombo, P.; Grosse, F.
2015-11-01
Mechanisms for the appearance of disorder on the As-rich GaAs(111)A surface were investigated employing density functional theory (DFT). Focus was given to the As trimer interactions by considering different surface symmetries and rest site occupations. The (2 × 2) and the c(4 × 2) structure models with As trimer and an As rest site were found the most energetically stable under the As-rich experimental conditions at T = 0 K. Low interactions between neighboring As trimers causes disorder in thermodynamic equilibrium at finite temperatures. A careful analysis of the configurational entropy contributions including the different statistics was carried out. The experimentally observed As-rich (2 × 2) structure was confirmed to be kinetically stabilized. The stabilization mechanism is discussed with respect to the As trimer migration on the surface, which is limited by a large diffusion barrier through the As rest sites.
Ab-initio calculation of C and CO adsorption on the Co (110) surface
NASA Astrophysics Data System (ADS)
Chin, Shin-Liang; Ionescu, Adrian; Reeve, Robert M.; Cheng, Jun; Barnes, Crispin H. W.
2013-02-01
The adsorption energies, structural, electrical and magnetic properties of adsorption of C and CO on the fcc Co (110) surface have been investigated using density functional theory. The preferential adsorption site for the fcc Co (110) has been calculated. For the case of C adsorption, the preferential adsorption site is the long-bridge for both the 0.5 and 1.0 monolayers (ML) coverages. Whilst for the CO case, the preferential adsorption sites are at the atop and short-bridge site for 0.5 and 1.0 ML coverages respectively. Structurally, the first two layers of the bare Co (110) surface expand whereas the second and third layers contract. Upon adsorption of either C or CO, however, the degree of expansion and compression reduces. Magnetically, the adsorbates were found to couple ferrimagnetically to the surface and suppress the magnetic moment of the Co layers beneath them. The C adsorbate has a much stronger suppression effect as compared to the CO adsorbate. At 0.5 ML coverage, the C adatom suppresses up to 47% of the magnetic moment in the surface layer compared to a clean Co (110), whereas the CO adsorbate only suppresses up to 16%. For the 1.0 ML coverage case, both the C and CO adsorbates suppress almost equivalently well at 68% and 63% respectively. We also report on a correlation between the amount of charge transfer and the degree of suppression of the surface magnetic moment. Finally, we observe that the electronic charge is shared in the [001] direction for the C adsorbate and in the [11¯0] direction for the CO adsorbate. This anisotropy of surface bonding, in conjunction with the ligand field theory, explains the mechanism behind the spin reorientation transition that occurs uniquely on the CO/Co(110) system.
Kokalj, Anton
2015-01-01
The bonding of benzotriazole-an outstanding corrosion inhibitor for copper-on reduced and oxidized copper surfaces is discussed on the basis of density functional theory (DFT) calculations. Calculations reveal that benzotriazole is able to bond with oxide-free and oxidized copper surfaces and on both of them it bonds significantly stronger to coordinatively unsaturated Cu sites. This suggests that benzotriazole is able to passivate the reactive under-coordinated surface sites that are plausible microscopic sites for corrosion attack. Benzotriazole can adsorb in a variety of different forms, yet it forms a strong molecule-surface bond only in deprotonated form. The bonding is even stronger when the deprotonated form is incorporated into organometallic adcomplexes. This is consistent with existing experimental evidence that benzotriazole inhibits corrosion by forming protective organometallic complexes. It is further shown that adsorption of benzotriazole considerably reduces the metal work function, which is a consequence of a large permanent molecular dipole and a properly oriented adsorption structure. It is argued that such a pronounced effect on the work function might be relevant for corrosion inhibition, because it should diminish the anodic corrosion reaction, which is consistent with existing experimental evidence that benzotriazole, although a mixed type inhibitor, predominantly affects the anodic reaction. PMID:25955130
Ab initio study of silane and disilane adsorption on Si(100)-(2x1) surface
NASA Astrophysics Data System (ADS)
Huang, Min; Randall, John; Chabal, Yves J.; Wallace, Robert M.; Cho, Kyeongjae
2010-03-01
Silane (SiH4) and disilane (Si2H6) are common precursors for the growth of Si and SiO2 thin films for microelectronic and photovoltaic devices. The adsorption of silane and disilane on Si(100)-(2x1) surface, which are important steps in the growth of Si films in atomic layer epitaxy (ALE), were investigated using density functional theory calculations. The silane molecule dissociates on the Si surface at the intra-dimer site with barrier energy of 0.22 eV. We investigate both Si-Si bond cleavage and Si-H bond cleavage mechanisms for adsorption of Si2H6 on Si (100) surface. A Si-H bond cleavage mechanism was found to be more favored than Si-Si bond cleavage mechanism due to the lower barrier energy of 0.04 eV. The lower barrier energy for Si2H6 dissociation than that of SiH4 agrees well with the experimental results showing that Si2H6 has higher sticking coefficient than SiH4 on Si surface at 300K. The vibration frequencies of Si2H5, SiH3, SiH2, SiH resulting from dissociation of silane and disilane were calculated and compared with experimental results available. The simulation results will facilitate the controlled ALE for atomically precise manufacturing applications.
Rotational Energy Transfer of N2 Gas Determined Using a New Ab Initio Potential Energy Surface
NASA Technical Reports Server (NTRS)
Huo, Winifred M.; Stallcop, James R.; Partridge, Harry; Langhoff, Stephen R. (Technical Monitor)
1997-01-01
Rotational energy transfer between two N2 molecules is a fundamental process of some importance. Exchange is expected to play a role, but its importance is somewhat uncertain. Rotational energy transfer cross sections of N2 also have applications in many other fields including modeling of aerodynamic flows, laser operations, and linewidth analysis in nonintrusive laser diagnostics. A number of N2-N2 rigid rotor potential energy surface (PES) has been reported in the literature.
Ab initio cluster model study of the chemisorptions of CO on low-index platinum surfaces
Curulla, D.; Clotet, A.; Ricart, J.M.; Illas, F.
1999-06-24
A systematic theoretical study of the adsorption of CO on the Pt{l_brace}100{r_brace}, Pt{l_brace}110{r_brace}, and Pt{l_brace}111{r_brace} surfaces is presented. The calculated equilibrium geometries and vibrational frequencies have been found to be rather independent of the cluster model chosen to represent the surface. However, calculated interaction energies are found to be very sensitive to the surface cluster model. The analysis of the chemisorption bond has been carried out by means of the constrained space orbital variation, CSOV, and of projection operator techniques. These analysis reveal that the bonding interactions are dominated by the {pi}-back-donation although {sigma}-donation plays a significant role. It is also clearly shown that all bonding mechanisms, other than Pauli repulsion, but specially {pi}-back-donation, contribute to the observed red shift. However, the {pi}-back-donation contribution to the red shift is very similar for CO on different sites. Hence, {pi}-back-donation cannot be the mechanism responsible for the observed difference for the CO vibrational frequency on on-top and bridge sites. The CSOV decomposition reveals that the leading term contributing to this difference in vibrational frequency of chemisorbed CO is the initial Pauli repulsion or wall effect; this is a new, important and unexpected conclusion.
Du Weina; Luo Cheng; Li Zesheng
2008-11-07
The dissociations of CH{sub 2}CHCN into CH{sub 2}CH+CN and CH{sub 2}C+HCN in the S{sub 0}, T{sub 1}, and {sup 1}{pi}{sub 2}{pi}{sub C{identical_to}}{sub N}* (definitions of {pi} orbitals can refer to computational details) states, have been explored at the complete active space self-consistent field level of theory employing the Dunning correlation consistent triple-zeta basis set. The lowest energy points of the surface crossing seams have been searched. Two conical intersections, from {sup 1}{pi}{sub C{identical_to}}{sub N}{pi}{sub 1}* to {sup 1}{pi}{sub 2}{pi}{sub 1}* (CI{sub 1}) and from {sup 1}{pi}{sub 2}{pi}{sub 1}* to S{sub 0} (CI{sub 2}), and one intersystem crossing point (T{sub 1}/S{sub 0}) have been located. The energies of all critical points have been recomputed with the multiconfigurational second-order perturbation method. At each conical intersection, derivative coupling and unscaled gradient difference vectors have been analyzed to determine the relaxation channels that the molecule may evolve in after nonradiative decay. Once the molecule is photoexcited to the {sup 1}{pi}{sub 2}{pi}{sub 1}* or {sup 1}{pi}{sub C{identical_to}}{sub N}{pi}{sub 1}* state, it would relax along the similar pathway: funneling through CI{sub 1} and then CI{sub 2}, and finally populate the ground state. Our results show that upon 193 nm photoexcitation, the most probable reaction channel is the ground-state HCN elimination following radiationless decays from excited states through surface crossings, which consists with experimental results J. Chem. Phys. 108, 5784 (1998). The investigated dissociation channels on the {sup 1}{pi}{sub 2}{pi}{sub C{identical_to}}{sub N}* surface, which are inaccessible upon 193 nm photoexcitation, may provide information for reactions induced by higher energy excitations.
NASA Astrophysics Data System (ADS)
Marqués, R.; Baena, J. D.; Beruete, M.; Falcone, F.; Lopetegi, T.; Sorolla, M.; Martín, F.; Garcia, J.
2005-02-01
Frequency selective surfaces (FSSs) made up of periodic arrays of split ring resonators (SRRs) are analysed. This analysis includes complementary screens, or complementary SRR-FSSs (CSRR-FSSs). It is shown that these FSSs show a dual behaviour, with a stop/pass band behaviour at the frequency of resonance of the SRRs/CSRRs. Cross-polarization effects in the SRR and the CSRR are considered, and it is shown that they permit resonance to occur for normally incident plane wave excitation. This latter property of SRRs and CSRRs also implies that the FSSs considered may act as polarizers and polarization converters as well. An analytical theory, valid for perfectly conducting and infinitely thin screens, is proposed for the SRR-FSSs and CSRR-FSSs. These approximations are valid in the microwave and millimetre-wave range, and up to the terahertz range.
Ab initio constructed diabatic surfaces of NO2 and the photodetachment spectra of its anion.
Mukherjee, Saikat; Mukherjee, Bijit; Sardar, Subhankar; Adhikari, Satrajit
2015-12-28
A thorough investigation has been performed for electronic structure, topological effect, and nuclear dynamics of NO2 molecule, where the adiabatic potential energy surfaces (PESs), conical intersections between the ground (X(2)A1) and the first excited state (A(2)B2), and the corresponding non-adiabatic coupling terms between those states are recalculated [Chem. Phys. 416, 11 (2013)] to achieve enough accuracy in dynamics. We employ beyond Born-Oppenheimer theory for these two state sub-Hilbert space to carry out adiabatic to diabatic transformation (ADT) to obtain the ADT angles and thereby, to construct single-valued, smooth, and continuous diabatic PESs. The analytic expressions for the adiabatic PESs and ADT angles are provided to represent a two-state three-mode diabatic Hamiltonian of NO2 for performing nuclear dynamics to calculate the photo-electron spectra of its anion. It appears that not only Jahn-Teller type coupling but also Renner-Teller interaction contributes significantly on the overall spectrum. The coupling between the electronic states (X(2)A1 and A(2)B2) of NO2 is essentially through the asymmetric stretching mode, where the functional form of such interaction is distinctly symmetric and non-linear. PMID:26723671
Ab initio constructed diabatic surfaces of NO2 and the photodetachment spectra of its anion
NASA Astrophysics Data System (ADS)
Mukherjee, Saikat; Mukherjee, Bijit; Sardar, Subhankar; Adhikari, Satrajit
2015-12-01
A thorough investigation has been performed for electronic structure, topological effect, and nuclear dynamics of NO2 molecule, where the adiabatic potential energy surfaces (PESs), conical intersections between the ground (X2A1) and the first excited state (A2B2), and the corresponding non-adiabatic coupling terms between those states are recalculated [Chem. Phys. 416, 11 (2013)] to achieve enough accuracy in dynamics. We employ beyond Born-Oppenheimer theory for these two state sub-Hilbert space to carry out adiabatic to diabatic transformation (ADT) to obtain the ADT angles and thereby, to construct single-valued, smooth, and continuous diabatic PESs. The analytic expressions for the adiabatic PESs and ADT angles are provided to represent a two-state three-mode diabatic Hamiltonian of NO2 for performing nuclear dynamics to calculate the photo-electron spectra of its anion. It appears that not only Jahn-Teller type coupling but also Renner-Teller interaction contributes significantly on the overall spectrum. The coupling between the electronic states (X2A1 and A2B2) of NO2 is essentially through the asymmetric stretching mode, where the functional form of such interaction is distinctly symmetric and non-linear.
Ab initio pseudopotential band calculation of organic conductors
Ishibashi, Shoji; Kohyama, Masanori
1999-12-01
The authors have calculated the band structures of organic conductors TTF-TCNQ and {beta}-(BEDT-TTF){sub 2}I{sub 3} using the ab initio plane-wave pseudopotential method within the local-density approximation (LDA). The Fermi-surface shape and the origin of bands near the Fermi level are investigated for each compound.
NASA Astrophysics Data System (ADS)
Yung, Y. Y.; Choi, B. H.; Tang, K. T.
1980-01-01
Three dimensional quantum mechanical calculations are carried out for the reactive scattering of D+H2→DH+H on the ab initio potential energy surface calculated by Liu and Siegbahn and fitted by Truhlar and Horowitz. The differential and total cross sections as well as the S matrix elements are obtained from the adiabatic distorted wave method. Threshold energy, cross sections and product distributions over final states are all in good agreement with experimental measurements. Results are also compared with the corresponding ones obtained on the Porter-Karplus and the Yates-Lester semi-empirical surfaces.
NASA Astrophysics Data System (ADS)
Ryabishchenkova, A. G.; Otrokov, M. M.; Kuznetsov, V. M.; Chulkov, E. V.
2015-09-01
Ab initio study of the adsorption, diffusion, and intercalation of alkali metal adatoms on the (0001) step surface of the topological insulator Bi2Se3 has been performed for the case of low coverage. The calculations of the activation energies of diffusion of adatoms on the surface and in van der Waals gaps near steps, as well as the estimate of diffusion lengths, have shown that efficient intercalation through steps is possible only for Li and Na. Data obtained for K, Rb, and Cs atoms indicate that their thermal desorption at high temperatures can occur before intercalation. The results have been discussed in the context of existing experimental data.
NASA Astrophysics Data System (ADS)
Bo, Tao; Lan, Jian-Hui; Zhao, Yao-Lin; Zhang, Yu-Juan; He, Chao-Hui; Chai, Zhi-Fang; Shi, Wei-Qun
2016-02-01
We have investigated the low-index (111), (110), and (100) surfaces of NpO2 as well as the adsorption and dissociation behaviors of water on these surfaces by using DFT + U method. The calculated surface energies for the (111), (110), and (100) surfaces are 0.81, 1.14, and 1.67 J m- 2, respectively. Based on our calculations, the presence of surface oxygen vacancy enhances the dissociation adsorption of water molecule. There is no significant coverage dependence for molecular adsorption of water on the NpO2 (111), (110), and (100) surfaces. The average adsorption energy of water in dissociative states increases slightly from 1/4 to 2/4 ML and rises significantly with the coverage increase to 1 ML on these surfaces. In addition, we have analyzed the water adsorption on the NpO2 surfaces for different temperatures and H2O partial pressures and have plotted the pressure-temperature phase diagrams by using the "ab initio atomistic thermodynamics" approach.
Kalugina, Yulia N; Lokshtanov, Sergei E; Cherepanov, Victor N; Vigasin, Andrey A
2016-02-01
We present new three-dimensional potential energy surface (PES) and dipole moment surfaces (DMSs) for the CH4-Ar van der Waals system. Ab initio calculations of the PES and DMS were carried out using the closed-shell single- and double-excitation coupled cluster approach with non-iterative perturbative treatment of triple excitations. The augmented correlation-consistent aug-cc-pVXZ (X = D,T,Q) basis sets were employed, and the energies obtained were then extrapolated to the complete basis set limit. The dipole moment surface was obtained using aug-cc-pVTZ basis set augmented with mid-bond functions for better description of exchange interactions. The second mixed virial coefficient was calculated and compared to available experimental data. The equilibrium constant for true dimer formation was calculated using classical partition function based on the knowledge of ab initio PES. Temperature variations of the zeroth spectral moment of the rototranslational collision-induced band as well as its true dimer constituent were traced with the use of the Boltzmann-weighted squared induced dipole properly integrated over respective phase space domains. Height profiles for N2-N2, N2-H2, CH4-N2, (CH4)2, and CH4-Ar true bound dimers in Titan's atmosphere were calculated with the use of reliable ab initio PESs. PMID:26851918
Phonocatalysis. An ab initio simulation experiment
NASA Astrophysics Data System (ADS)
Kim, Kwangnam; Kaviany, Massoud
2016-06-01
Using simulations, we postulate and show that heterocatalysis on large-bandgap semiconductors can be controlled by substrate phonons, i.e., phonocatalysis. With ab initio calculations, including molecular dynamic simulations, the chemisorbed dissociation of XeF6 on h-BN surface leads to formation of XeF4 and two surface F/h-BN bonds. The reaction pathway and energies are evaluated, and the sorption and reaction emitted/absorbed phonons are identified through spectral analysis of the surface atomic motion. Due to large bandgap, the atomic vibration (phonon) energy transfer channels dominate and among them is the match between the F/h-BN covalent bond stretching and the optical phonons. We show that the chemisorbed dissociation (the pathway activation ascent) requires absorption of large-energy optical phonons. Then using progressively heavier isotopes of B and N atoms, we show that limiting these high-energy optical phonons inhibits the chemisorbed dissociation, i.e., controllable phonocatalysis.
Li, Y. Q.; Zhang, P. Y.; Han, K. L.
2015-03-28
A global many-body expansion potential energy surface is reported for the electronic ground state of CH{sub 2}{sup +} by fitting high level ab initio energies calculated at the multireference configuration interaction level with the aug-cc-pV6Z basis set. The topographical features of the new global potential energy surface are examined in detail and found to be in good agreement with those calculated directly from the raw ab initio energies, as well as previous calculations available in the literature. In turn, in order to validate the potential energy surface, a test theoretical study of the reaction CH{sup +}(X{sup 1}Σ{sup +})+H({sup 2}S)→C{sup +}({sup 2}P)+H{sub 2}(X{sup 1}Σ{sub g}{sup +}) has been carried out with the method of time dependent wavepacket on the title potential energy surface. The total integral cross sections and the rate coefficients have been calculated; the results determined that the new potential energy surface can both be recommended for dynamics studies of any type and as building blocks for constructing the potential energy surfaces of larger C{sup +}/H containing systems.
Ab initio infrared and Raman spectra
NASA Astrophysics Data System (ADS)
Fredkin, Donald R.; Komornicki, Andrew; White, Steven R.; Wilson, Kent R.
1983-06-01
We discuss several ways in which molecular absorption and scattering spectra can be computed ab initio, from the fundamental constants of nature. These methods can be divided into two general categories. In the first, or sequential, type of approach, one first solves the electronic part of the Schrödinger equation in the Born-Oppenheimer approximation, mapping out the potential energy, dipole moment vector (for infrared absorption) and polarizability tensor (for Raman scattering) as functions of nuclear coordinates. Having completed the electronic part of the calculation, one then solves the nuclear part of the problem either classically or quantum mechanically. As an example of the sequential ab initio approach, the infrared and Raman rotational and vibrational-rotational spectral band contours for the water molecule are computed in the simplest rigid rotor, normal mode approximation. Quantum techniques are used to calculate the necessary potential energy, dipole moment, and polarizability information at the equilibrium geometry. A new quick, accurate, and easy to program classical technique involving no reference to Euler angles or special functions is developed to compute the infrared and Raman band contours for any rigid rotor, including asymmetric tops. A second, or simultaneous, type of ab initio approach is suggested for large systems, particularly those for which normal mode analysis is inappropriate, such as liquids, clusters, or floppy molecules. Then the curse of dimensionality prevents mapping out in advance the complete potential, dipole moment, and polarizability functions over the whole space of nuclear positions of all atoms, and a solution in which the electronic and nuclear parts of the Born-Oppenheimer approximation are simultaneously solved is needed. A quantum force classical trajectory (QFCT) molecular dynamic method, based on linear response theory, is described, in which the forces, dipole moment, and polarizability are computed quantum
M Kondo; T Mates; D Fischer; F Wudl; E Kramer
2011-12-31
Interfaces between phenylacetylene (PA) monolayers and two silicon surfaces, Si(111) and Si(100), are probed by X-ray photoelectron spectroscopy (XPS), ultraviolet photoelectron spectroscopy (UPS), and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy, and the results are analyzed using ab initio molecular orbital calculations. The monolayer systems are prepared via the surface hydrosilylation reaction between PA and hydrogen-terminated silicon surfaces. The following spectral features are obtained for both of the PA-Si(111) and PA-Si(100) systems: a broad {pi}-{pi}* shakeup peak at 292 eV (XPS), a broad first ionization peak at 3.8 eV (UPS), and a low-energy C 1s {yields} {pi}* resonance peak at 284.3 eV (NEXAFS). These findings are ascribed to a styrene-like {pi}-conjugated molecular structure at the PA-Si interface by comparing the experimental data with theoretical analysis results. A conclusion is drawn that the vinyl group can keep its {pi}-conjugation character on the hydrogen-terminated Si(100) [H:Si(100)] surface composed of the dihydride (SiH{sub 2}) groups as well as on hydrogen-terminated Si(111) having the monohydride (SiH) group. The formation mechanism of the PA-Si(100) interface is investigated within cluster ab initio calculations, and the possible structure of the H:Si(100) surface is discussed based on available data.
Song, Lingchun; Gao, Jiali
2009-01-01
A theoretical model is presented for deriving effective diabatic states based on ab initio self-consistent field valence bond (VBSCF) theory by reducing the multi-configurational VB Hamiltonian into an effective two-state model. We describe two computational approaches for the optimization of the effective diabatic configurations, resulting in two ways of interpreting such effective diabatic states. In the variational diabatic configuration (VDC) method, the energies of the individual diabatic states are variationally minimized. In the consistent diabatic configuration (CDC) method, both the configuration coefficients and orbital coefficients are simultaneously optimized to minimize the adiabatic ground-state energy in VBSCF calculations. In addition, we describe a mixed molecular orbital and valence bond (MOVB) approach to construct the CDC diabatic and adiabatic states for a chemical reaction, whereas the VDC-MOVB method has been described previously. Employing the symmetric SN2 reaction between NH3 and CH3NH3+ as a test system, we found that the results from ab initio VBSCF and from MOVB calculations are in good agreement, suggesting that the computationally efficient MOVB method is a reasonable model for VB simulations of condensed phase reactions. The results indicate that CDC and VDC diabatic states converge, respectively, to covalent and ionic states as the molecular geometries are distorted from the minimum of the respective diabatic state along the reaction coordinate. Furthermore, the resonance energy that stabilizes the energy of crossing between the two diabatic states, resulting in the transition state of the adiabatic ground-state reaction, has a strong dependence on the overlap integral between the two diabatic states and is a function of both the exchange integral and the total diabatic ground-state energy. PMID:18828577
Hua -Gen Yu; Han, Huixian; Guo, Hua
2016-03-29
Vibrational energy levels of the ammonium cation (NH4+) and its deuterated isotopomers are calculated using a numerically exact kinetic energy operator on a recently developed nine-dimensional permutation invariant semiglobal potential energy surface fitted to a large number of high-level ab initio points. Like CH4, the vibrational levels of NH4+ and ND4+ exhibit a polyad structure, characterized by a collective quantum number P = 2(v1 + v3) + v2 + v4. As a result, the low-lying vibrational levels of all isotopomers are assigned and the agreement with available experimental data is better than 1 cm–1.
Collective rotation from ab initio theory
NASA Astrophysics Data System (ADS)
Caprio, Mark A.; Maris, Pieter; Vary, James P.
2015-10-01
The challenge of ab initio nuclear theory is to quantitatively predict the complex and highly-correlated behavior of the nuclear many-body system, starting from the underlying internucleon interactions. We may now seek to understand the wealth of nuclear collective phenomena through ab initio approaches. No-core configuration interaction (NCCI) calculations for p-shell nuclei give rise to rotational bands, as evidenced by rotational patterns for excitation energies, electromagnetic moments, and electromagnetic transitions. In this talk, the intrinsic structure of these bands is discussed, and the predicted rotational bands are compared to experiment. Supported by the US DOE under Award Nos. DE-FG02-95ER-40934, DESC0008485 (SciDAC/NUCLEI), and DE-FG02-87ER40371 and the US NSF under Award No. 0904782. Computational resources provided by NERSC (US DOE Contract No. DE-AC02-05CH11231).
Ab initio infrared and Raman spectra
NASA Technical Reports Server (NTRS)
Fredkin, D. R.; White, S. R.; Wilson, K. R.; Komornicki, A.
1983-01-01
It is pointed out that with increased computer power and improved computational techniques, such as the gradients developed in recent years, it is becoming practical to compute spectra ab initio, from the fundamental constants of nature, for systems of increasing complexity. The present investigation has the objective to explore several possible ab initio approaches to spectra, giving particular attention to infrared and nonresonance Raman. Two approaches are discussed. The sequential approach, in which first the electronic part and then later the nuclear part of the Born-Oppenheimer approximation is solved, is appropriate for small systems. The simultaneous approach, in which the electronic and nuclear parts are solved at the same time, is more appropriate for many-atom systems. A review of the newer quantum gradient techniques is provided, and the infrared and Raman spectral band contours for the water molecule are computed.
NASA Astrophysics Data System (ADS)
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.
Ab initio Bogoliubov coupled cluster theory
NASA Astrophysics Data System (ADS)
Signoracci, Angelo; Hagen, Gaute; Duguet, Thomas
2014-09-01
Coupled cluster (CC) theory has become a standard method in nuclear theory for realistic ab initio calculations of medium mass nuclei, but remains limited by its requirement of a Slater determinant reference state which reasonably approximates the nuclear system of interest. Extensions of the method, such as equation-of-motion CC, permit the calculation of nuclei with one or two nucleons added or removed from a doubly magic core, yet still only a few dozen nuclei are accessible with modern computational restrictions. In order to extend the applicability of ab initio methods to open-shell systems, the superfluid nature of nuclei must be taken into account. By utilizing Bogoliubov algebra and employing spontaneous symmetry breaking with respect to particle number conservation, superfluid systems can be treated by a single reference state. An ab initio theory to include correlations on top of a Bogoliubov reference state has been developed in the guise of standard CC theory. The formalism and first results of this Bogoliubov coupled cluster theory will be presented to demonstrate the applicability of the method.
Makarewicz, Jan; Shirkov, Leonid
2016-05-28
The pyridine-Ar (PAr) van der Waals (vdW) complex is studied using a high level ab initio method. Its structure, binding energy, and intermolecular vibrational states are determined from the analytical potential energy surface constructed from interaction energy (IE) values computed at the coupled cluster level of theory with single, double, and perturbatively included triple excitations with the augmented correlation consistent polarized valence double-ζ (aug-cc-pVDZ) basis set complemented by midbond functions. The structure of the complex at its global minimum with Ar at a distance of 3.509 Å from the pyridine plane and shifted by 0.218 Å from the center of mass towards nitrogen agrees well with the corresponding equilibrium structure derived previously from the rotational spectrum of PAr. The PAr binding energy De of 392 cm(-1) is close to that of 387 cm(-1) calculated earlier at the same ab initio level for the prototypical benzene-Ar (BAr) complex. However, under an extension of the basis set, De for PAr becomes slightly lower than De for BAr. The ab initio vdW vibrational energy levels allow us to estimate the reliability of the methods for the determination of the vdW fundamentals from the rotational spectra. To disclose the character of the intermolecular interaction in PAr, the symmetry-adapted perturbation theory (SAPT) is employed for the analysis of different physical contributions to IE. It is found that SAPT components of IE can be approximately expressed in the binding region by only two of them: the exchange repulsion and dispersion energy. The total induction effect is negligible. The interrelations between various SAPT components found for PAr are fulfilled for a few other complexes involving aromatic molecules and Ar or Ne, which indicates that they are valid for all rare gas (Rg) atoms and aromatics. PMID:27250287
Li, Xiaoping; Mandal, Anirban; Miliordos, Evangelos; Hunt, Katharine L C
2012-01-28
We report new ab initio results for the interaction-induced dipole moments Δμ of hydrogen molecules colliding with helium atoms. These results are needed in order to calculate collision-induced absorption spectra at high temperatures; applications include modeling the radiative profiles of very cool white dwarf stars, with temperatures from 3500 K to 9000 K. We have evaluated the dipoles based on finite-field calculations, with coupled cluster methods in MOLPRO 2006 and aug-cc-pV5Z (spdfg) basis sets for both the H and He centers. We have obtained values of Δμ for eight H(2) bond lengths ranging from 0.942 a.u. to 2.801 a.u., for 15 intermolecular separations R ranging from 2.0 a.u. to 10.0 a.u., and for 19 different relative orientations. In general, our values agree well with earlier ab initio results, for the geometrical configurations that are treated in common, but we have determined more points on the collision-induced dipole surface by an order of magnitude. These results make it possible to calculate transition probabilities for molecules in excited vibrational states, overtones, and rotational transitions with ΔJ > 4. We have cast our results in the symmetry-adapted form needed for absorption line shape calculations, by expressing Δμ as a series in the spherical harmonics of the orientation angles of the intermolecular vector and of a unit vector along the H(2) bond axis. The expansion coefficients depend on the H(2) bond length and the intermolecular distance R. For large separations R, we show that the ab initio values of the leading coefficients converge to the predictions from perturbation theory, including both classical multipole polarization and dispersion effects. PMID:22299884
NASA Astrophysics Data System (ADS)
Makarewicz, Jan; Shirkov, Leonid
2016-05-01
The pyridine-Ar (PAr) van der Waals (vdW) complex is studied using a high level ab initio method. Its structure, binding energy, and intermolecular vibrational states are determined from the analytical potential energy surface constructed from interaction energy (IE) values computed at the coupled cluster level of theory with single, double, and perturbatively included triple excitations with the augmented correlation consistent polarized valence double-ζ (aug-cc-pVDZ) basis set complemented by midbond functions. The structure of the complex at its global minimum with Ar at a distance of 3.509 Å from the pyridine plane and shifted by 0.218 Å from the center of mass towards nitrogen agrees well with the corresponding equilibrium structure derived previously from the rotational spectrum of PAr. The PAr binding energy De of 392 cm-1 is close to that of 387 cm-1 calculated earlier at the same ab initio level for the prototypical benzene-Ar (BAr) complex. However, under an extension of the basis set, De for PAr becomes slightly lower than De for BAr. The ab initio vdW vibrational energy levels allow us to estimate the reliability of the methods for the determination of the vdW fundamentals from the rotational spectra. To disclose the character of the intermolecular interaction in PAr, the symmetry-adapted perturbation theory (SAPT) is employed for the analysis of different physical contributions to IE. It is found that SAPT components of IE can be approximately expressed in the binding region by only two of them: the exchange repulsion and dispersion energy. The total induction effect is negligible. The interrelations between various SAPT components found for PAr are fulfilled for a few other complexes involving aromatic molecules and Ar or Ne, which indicates that they are valid for all rare gas (Rg) atoms and aromatics.
Paul, Amit Kumar; Ray, Somrita; Mukhopadhyay, Debasis; Adhikari, Satrajit
2011-07-21
We perform ab initio calculation using quantum chemistry package (MOLPRO) on the excited states of Na(3) cluster and present the adiabatic PESs for the electronic states 2(2)E' and 1(2)A(1)', and the non-adiabatic coupling (NAC) terms among those states. Since the ab initio calculated NAC elements for the states 2(2)E' and 1(2)A(1)' demonstrate the numerical validity of so called "Curl Condition," such states closely form a sub-Hilbert space. For this subspace, we employ the NAC terms to solve the "adiabatic-diabatic transformation (ADT)" equations to obtain the functional form of the transformation angles and pave the way to construct the continuous and single valued diabatic potential energy surface matrix by exploiting the existing first principle based theoretical means on beyond Born-Oppenheimer treatment. Nuclear dynamics has been carried out on those diabatic surfaces to reproduce the experimental spectrum for system B of Na(3) cluster and thereby, to explore the numerical validity of the theoretical development on beyond Born-Oppenheimer approach for adiabatic to diabatic transformation. PMID:21786987
Skutterudites under pressure: An ab initio study
Ram, Swetarekha; Kanchana, V.; Valsakumar, M. C.
2014-03-07
Ab initio results on the band structure, density of states, and Fermi surface (FS) properties of LaRu{sub 4}X{sub 12} (X = P, As, Sb) are presented at ambient pressure as well as under compression. The analysis of density of states reveals the major contribution at the Fermi level to be mainly from the Ru-d and X-p states. We have a complicated Fermi surface with both electron and hole characters for all the three compounds which is derived mainly from the Ru-d and X-p states. There is also a simpler FS with hole character derived from the P-p{sub z} orbital for LaRu{sub 4}P{sub 12} and Ru-d{sub z{sup 2}} orbital in the case of As and Sb containing compounds. More interestingly, Fermi surface nesting feature is observed only in the case of the LaRu{sub 4}P{sub 12}. Under compression, we observe the topology of the complicated FS sheet of LaRu{sub 4}As{sub 12} to change around V/V{sub 0} = 0.85, leading to a behaviour similar to that of a multiband superconductor, and in addition, we have two more hole pockets centered around Γ at V/V{sub 0} = 0.8 for the same compound. Apart from this, we find the hole pocket to vanish at V/V{sub 0} = 0.8 in the case of LaRu{sub 4}Sb{sub 12} and the opening of the complicated FS sheet gets reduced. The de Haas van Alphen calculation shows the number of extremal orbits in the complicated sheet to change in As and Sb containing compounds under compression, where we also observe the FS topology to change.
NASA Astrophysics Data System (ADS)
Gonzalez, Dayana; Mebel, Alexander
2016-03-01
It has been recently shown that Titan provides a unique perspective in our solar system: its atmosphere is comparable to a model of prebiotic Earth's. Provided the organic cationic and anionic molecular species identified by the Cassini spacecraft, this research characterizes reaction pathways for the reactions of methyl derivatives of the cyclopropenyl cation, the methyl cation with methyl- and dimethyl-acetylene, and reactions of resonance structures of protonated acrylonitrile with CH2NH. Isomerization and dissociation reactions involving methyl-cyclopropenyl cations, the perinaphthenyl cation and anion, and cations of pyrimidine and purine precursors of nucleobases will be examined to locate reaction pathways, intermediates, transition states, and products of the reactions. Gaussian '09 software is used for ab initio calculations to map out the PES. Geometry optimizations and vibrational frequency computations are preformed via the double-hybrid density functional B2PLYP-D3. Single-point energies are refined by use of the explicitly-correlated coupled-cluster CCSD(T)-F12 method. Rate constants are calculated using microcanonical RRKM theory, and pressure effects evaluated used the Master Equation approach; these allow for prediction of absolute rate constants and product branching ratios at different pressures and temperatures.
NASA Astrophysics Data System (ADS)
Kozin, Igor N.; Jensen, Per; Li, Yan; Buenker, Robert J.; Hirsch, Gerhard; Klee, Stefan
1997-01-01
The present work reports an ab initioMRD-CI calculation of the dipole moment surfaces for the electronic ground state of the H 2Te molecule. Using the ab initioresults, we calculate the vibrational transition moments, and we simulate the far-infrared spectrum of H 2Te by means of the MORBID program system. We obtain the equilibrium value of the dipole moment from the ab initiocalculation as 0.377 Debye based on our initial theoretical treatment which was employed over a wide range of molecular geometries. However, the use of an improved AO basis at the equilibrium geometry of H 2Te lowers this result to 0.298 Debye. The comparison of our simulated far-infrared spectrum with the experimental spectrum suggests that this value is too large, and that the correct value is certainly larger than 0.19 Debye and very probably smaller than 0.26 Debye. From the ab initiodata, we predict many vibrational transition moments for H 2Te, D 2Te, and HDTe. We hope that these results will be of assistance in the interpretation of the rotation-vibration spectrum of these molecules.
Ab initio no core full configuration approach for light nuclei
NASA Astrophysics Data System (ADS)
Kim, Youngman; Shin, Ik Jae; Maris, Pieter; Vary, James P.; Forssén, Christian; Rotureau, Jimmy
2014-07-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.
Ab initio no core full configuration approach for light nuclei
NASA Astrophysics Data System (ADS)
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.
Ab Initio Electronic Relaxation Times and Transport in Noble Metals
NASA Astrophysics Data System (ADS)
Mustafa, Jamal I.; Bernardi, Marco; Neaton, Jeffrey B.; Louie, Steven G.
Relaxation times employed to study electron transport in metals are typically assumed to be constants and obtained empirically using the Drude model. Here, we employ ab initio calculations to compute the electron-phonon relaxation times of Cu, Ag, and Au, and find that they vary significantly on the Fermi surface, spanning ~15 -45 fs. We compute room temperature resistivities in excellent agreement with experiment by combining GW bandstructures, Wannier-interpolated band velocities, and ab initio relaxation times. Our calculations are compared to other approximations used for the relaxation times. Additionally, an importance sampling scheme is introduced to speed up the convergence of resistivity and transport calculations by sampling directly points on the Fermi surface. This work was supported by NSF Grant No. DMR15-1508412 and U.S. DOE under Contract No. DE-AC02-05CH11231. Computational resources have been provided by DOE at LBNL's NERSC facility.
GAUSSIAN 76: An ab initio Molecular Orbital Program
DOE R&D Accomplishments Database
Binkley, J. S.; Whiteside, R.; Hariharan, P. C.; Seeger, R.; Hehre, W. J.; Lathan, W. A.; Newton, M. D.; Ditchfield, R.; Pople, J. A.
1978-01-01
Gaussian 76 is a general-purpose computer program for ab initio Hartree-Fock molecular orbital calculations. It can handle basis sets involving s, p and d-type Gaussian functions. Certain standard sets (STO-3G, 4-31G, 6-31G*, etc.) are stored internally for easy use. Closed shell (RHF) or unrestricted open shell (UHF) wave functions can be obtained. Facilities are provided for geometry optimization to potential minima and for limited potential surface scans.
Ab initio melting curve of osmium
NASA Astrophysics Data System (ADS)
Burakovsky, L.; Burakovsky, N.; Preston, D. L.
2015-11-01
The melting curve of osmium up to a pressure P of 500 GPa is obtained from an extensive suite of ab initio quantum molecular dynamics (QMD) simulations using the Z method. The ab initio P =0 melting point of Os is 3370 ±75 K; this range encompasses all of the available data in the literature and corroborates the conclusion of J. W. Arblaster [Platinum Metals Rev. 49, 166 (2005)], 10.1595/147106705X70264 that the melting temperature of pure Os is 3400 ±50 K and that the 3300 K typically quoted in the literature is the melting point of impure Os. The T =0 equation of state (EOS) of Os and the P dependence of the optimized c /a ratio for the hexagonal unit cell, both to pressures ˜900 GPa, are obtained in the ab initio approach as validation of its use. Although excellent agreement with the available experimental data (P ≲80 GPa) is found, it is the third-order Birch-Murnaghan EOS with B0'=5 rather than the more widely accepted B0'=4 that describes the QMD data to higher pressures, in agreement with the more recent experimental EOS by Godwal et al. The theoretical melting curve of Os obtained earlier by Joshi et al. is shown to be inconsistent with our QMD results, and the possible reason for this discrepancy is suggested. Regularities in the melting curves of Os and five other third-row transition metals (Ta, W, Re, Pt, Au) could be used to estimate the currently unknown melting curves of Hf and Ir.
Ab initio quantum chemistry: Methodology and applications
Friesner, Richard A.
2005-01-01
This Perspective provides an overview of state-of-the-art ab initio quantum chemical methodology and applications. The methods that are discussed include coupled cluster theory, localized second-order Moller–Plesset perturbation theory, multireference perturbation approaches, and density functional theory. The accuracy of each approach for key chemical properties is summarized, and the computational performance is analyzed, emphasizing significant advances in algorithms and implementation over the past decade. Incorporation of a condensed-phase environment by means of mixed quantum mechanical/molecular mechanics or self-consistent reaction field techniques, is presented. A wide range of illustrative applications, focusing on materials science and biology, are discussed briefly. PMID:15870212
NASA Astrophysics Data System (ADS)
Otrokov, Mikhail; Ryabishchenkova, Anastasia; Gosalvez, Miguel Angel; Kuznetsov, Vladimir; Chulkov, Evgueni
We present the results of an ab initio study of the adsorption, diffusion, and intercalation of alkali metal adatoms on the (0001) stepped surface of the topological insulator Bi2Se3 for the case of low coverage. The calculations of the activation energies of the adatoms diffusion on the surface and in the van der Waals gaps near the steps, as well as the estimation of diffusion lengths, show that efficient intercalation through the steps is possible only for Li and Na. Data obtained for K, Rb, and Cs atoms indicate that their thermal desorption at high temperatures can start before intercalation. These results are discussed in the context of the experimental data available.
NASA Astrophysics Data System (ADS)
Saheer, V. C.; Kumar, Sanjay
2016-01-01
The global ground and first three excited electronic state adiabatic as well as the corresponding quasidiabatic potential energy surfaces is reported as a function of nuclear geometries in the Jacobi coordinates ( R → , r → , γ ) using Dunning's cc-pVTZ basis set at the internally contracted multi-reference (single and double) configuration interaction level of accuracy. Nonadiabatic couplings, arising out of relative motion of proton and the vibrational motion of CO, are also reported in terms of coupling potentials. The quasidiabatic potential energy surfaces and the coupling potentials have been obtained using the ab initio procedure [Simah et al., J. Chem. Phys. 111, 4523 (1999)] for the purpose of dynamics studies.
Peterson, K.A.; Skokov, S.; Bowman, J.M.
1999-10-01
A new, global analytical potential energy surface is constructed for the X&hthinsp;{sup 1}A{sup {prime}} electronic ground state of HOCl that accurately includes the HClO isomer. The potential is obtained by using accurate {ital ab initio} data from a previously published surface [Skokov {ital et al.}, J. Chem. Phys. {bold 109}, 2662 (1998)], as well as a significant number of new data for the HClO region of the surface at the same multireference configuration interaction, complete basis set limit level of theory. Vibrational energy levels and intensities are computed for both HOCl and HClO up to the OH+Cl dissociation limit and above the isomerization barrier. After making only minor adjustments to the {ital ab initio} surface, the errors with respect to experiment for HOCl are generally within a few cm{sup {minus}1} for 22 vibrational levels with the largest error being 26 cm{sup {minus}1}. A total of 813 bound vibrational states are calculated for HOCl. The HClO potential well supports 57 localized states, of which only the first 3 are bound. The strongest dipole transitions for HClO were computed for the fundamentals{emdash}33, 2.9, and 25 km/mol for {nu}{sub 1}, {nu}{sub 2}, and {nu}{sub 3}, respectively. From exact J=1 ro-vibrational calculations, state dependent rotational constants have been calculated for HClO. Lastly, resonance calculations with the new potential demonstrate that the presence of the HClO minimum has a negligible effect on the resonance states of HOCl near the dissociation threshold due to the relatively high and wide isomerization barrier. {copyright} {ital 1999 American Institute of Physics.}
Ab initio non-relativistic spin dynamics
Ding, Feizhi; Goings, Joshua J.; Li, Xiaosong; Frisch, Michael J.
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.
Ab initio non-relativistic spin dynamics
NASA Astrophysics Data System (ADS)
Ding, Feizhi; Goings, Joshua J.; Frisch, Michael J.; Li, Xiaosong
2014-12-01
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 Li3 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.
Ab Initio Neutron Drops with Chiral Hamiltonians
NASA Astrophysics Data System (ADS)
Potter, Hugh; Maris, Pieter; Vary, James
2015-04-01
Ab initio calculations for neutron drops are of interest for insights into neutron-rich nuclei and neutron star matter, and for examining the neutron-only sector of nucleon-nucleon and 3-nucleon interactions. I present ab initio results calculated using the no-core shell model with 2- and 3-body chiral Hamiltonians for neutron drops up to 20 neutrons confined in a 10 MeV harmonic trap. I discuss ground state energies, internal energies, radii, and evidence for pairing. In addition, excitation energies can be used to investigate the spin-orbit splittings in the p-shell and sd -shell. Prior Green's Function Monte Carlo calculations using the Argonne v8' potential with added 3-nucleon forces serve as a comparison. Supported by DOE Grants DESC0008485 (SciDAC/NUCLEI), DE-FG02-87ER40371, and NSF Grant 0904782; computational resources provided by the Oak Ridge Leadership Computing Facility (DOE Office of Science Contract DE-AC05-00OR22725) under an INCITE award.
NASA Astrophysics Data System (ADS)
Orellana, Walter
2014-07-01
The stability, electronic, and optical properties of (6,5) single-walled carbon nanotubes (CNTs) functionalized with free-base tetraphenylporphyrin (TPP) molecules through π-stacking interactions are studied by ab-initio calculations. The stability and optical response of the CNT-TPP compounds for increasing CNT-surface coverage are investigated. Our results show that four TPP molecules forming a ring around the CNT is the most stable configuration, showing strong binding energies of about 2.5 eV/TPP. However, this binding energy can increase even more after additional molecules assemble side by side along the CNT, favoring the formation of a full single layer of TPP, as experimentally suggested. The strong π-π attractive forces induce molecular distortions that move the TPP higher-occupied molecular orbital levels inside the CNT bandgap, changing the optical response of the TPP molecules stacked on the CNT.
NASA Astrophysics Data System (ADS)
Petrongolo, Carlo; Ranghino, Graziella; Scordamaglia, Raimondo
1980-01-01
The potential energy surface of the β-lactam + OH - reaction, related to the mode of action of β-lactam antibiotics, was investigated using the ab initio Hartree—Fock method with the STO-3G basis set. Three possible reaction paths for the B A C2 breaking of the amidic CN bond were obtained and discussed. The minimum-energy reaction path is characterized by the following processes: (1) the formation of a tetrahedral intermediate, ≈ 121 kcal mol -1 more stable than the reagents; (2) a barrier, ≈ 15 kcal mol -1 above the intermediate, which is mainly due to the partial breaking of the amidic bond; (3) the complete breaking of the amidic bond concerted with a proton transfer till the formation of the final product, ≈ 34 kcal mol -1 more stable than the intermediate. The evolution of some molecular orbitals and of the electron population along the reaction path was also discussed.
NASA Astrophysics Data System (ADS)
Novko, D.; Blanco-Rey, M.; Juaristi, J. I.; Alducin, M.
2015-11-01
The relaxation dynamics of hot H, N, and N2 on Pd(100), Ag(111), and Fe(110), respectively, is studied by means of ab initio molecular dynamics with electronic friction. This method is adapted here to account for the electron density changes caused by lattice vibrations, thus treating on an equal footing both electron-hole (e -h ) pair and phonon excitations. We find that even if the latter increasingly dominate the heavier is the hot species, the contribution of e -h pairs is by no means negligible in these cases because it gains relevance at the last stage of the relaxation process. The quantitative details of energy dissipation depend on the interplay of the potential energy surface, electronic structure, and kinetic factors.
Orellana, Walter
2014-07-14
The stability, electronic, and optical properties of (6,5) single-walled carbon nanotubes (CNTs) functionalized with free-base tetraphenylporphyrin (TPP) molecules through π-stacking interactions are studied by ab-initio calculations. The stability and optical response of the CNT-TPP compounds for increasing CNT-surface coverage are investigated. Our results show that four TPP molecules forming a ring around the CNT is the most stable configuration, showing strong binding energies of about 2.5 eV/TPP. However, this binding energy can increase even more after additional molecules assemble side by side along the CNT, favoring the formation of a full single layer of TPP, as experimentally suggested. The strong π-π attractive forces induce molecular distortions that move the TPP higher-occupied molecular orbital levels inside the CNT bandgap, changing the optical response of the TPP molecules stacked on the CNT.
Barron, Hector; Mottet, Christine; Weissker, Hans-Christian
2014-01-01
It is known that the surface-plasmon resonance (SPR) in small spherical Au nanoparticles of about 2 nm is strongly damped. We demonstrate that small Au nanorods of high aspect ratio develop a strong longitudinal SPR, of intensity comparable to that in Ag rods, as soon as the resonance energy drops below the onset of the interband transitions due to the geometry. We present ab initio calculations of time-dependent density-functional theory of rods with lengths of up to 7 nm. Changing length and width, not only the energy but also the character of the resonance in Au rods can be tuned. Moreover, the aspect ratio alone is not sufficient to predict even the character of the spectrum; the absolute size matters. PMID:24343038
Wu, Ying; Zhang, Chunfang; Cao, Jianwei; Bian, Wensheng
2014-06-19
Quasiclassical trajectory (QCT) calculations have been performed on a new global ab initio potential energy surface (PES) for the singlet ground state (1(1)A') of the CH2 reactive system. Our new PES can give a very good description of the well and asymptote regions, and particularly regions around conical intersections (CIs) and of van der Waals (vdW) interactions. The integral cross sections, differential cross sections, and product rovibrational state distributions for the C((1)D) + H2 → CH + H reaction have been investigated in a wide range of collision energies. The present integral cross sections are much larger than the previous QCT results at low collision energies, which can be attributed to the differences of the PESs in the regions around the CIs and vdW complexes. The thermal rate coefficients in the temperature range 200-1500 K have also been calculated and very good agreement with experiment is obtained. PMID:24878310
Timoshevskii, Vladimir; Feng, Zimin; Bevan, Kirk H; Zaghib, Karim
2015-08-26
The atomic and electronic structures of the LiFePO4 (LFP) surface, both bare and reconstructed upon possible oxygenation, are theoretically studied by ab initio methods. On the basis of total energy calculations, the atomic structure of the oxygenated surface is proposed, and the effect of surface reconstruction on the electronic properties of the surface is clarified. While bare LFP(010) surface is insulating, adsorption of oxygen leads to the emergence of semimetallic behavior by inducing the conducting states in the band gap of the system. The physical origin of these conducting states is investigated. We further demonstrate that deposition of Li2S layers on top of oxygenated LFP(010) surface leads to the formation of additional conducting hole states in the first layer of Li2S surface because of the charge transfer from sulfur p-states to the gap states of LFP surface. This demonstrates that oxygenated LFP surface not only provides conducting layers itself, but also induces conducting channels in the top layer of Li2S. These results help to achieve further understanding of potential role of LFP particles in improving the performance of Li-S batteries through emergent interface conductivity. PMID:26237114
Towards Accurate Ab Initio Predictions of the Spectrum of Methane
NASA Technical Reports Server (NTRS)
Schwenke, David W.; Kwak, Dochan (Technical Monitor)
2001-01-01
We have carried out extensive ab initio calculations of the electronic structure of methane, and these results are used to compute vibrational energy levels. We include basis set extrapolations, core-valence correlation, relativistic effects, and Born- Oppenheimer breakdown terms in our calculations. Our ab initio predictions of the lowest lying levels are superb.
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics
Makhov, Dmitry V.; Shalashilin, Dmitrii V.; Glover, William J.; Martinez, Todd J.
2014-08-07
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as “cloning,” in analogy to the “spawning” procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, “trains,” as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics
NASA Astrophysics Data System (ADS)
Makhov, Dmitry V.; Glover, William J.; Martinez, Todd J.; Shalashilin, Dmitrii V.
2014-08-01
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions.
Ab initio multiple cloning algorithm for quantum nonadiabatic molecular dynamics.
Makhov, Dmitry V; Glover, William J; Martinez, Todd J; Shalashilin, Dmitrii V
2014-08-01
We present a new algorithm for ab initio quantum nonadiabatic molecular dynamics that combines the best features of ab initio Multiple Spawning (AIMS) and Multiconfigurational Ehrenfest (MCE) methods. In this new method, ab initio multiple cloning (AIMC), the individual trajectory basis functions (TBFs) follow Ehrenfest equations of motion (as in MCE). However, the basis set is expanded (as in AIMS) when these TBFs become sufficiently mixed, preventing prolonged evolution on an averaged potential energy surface. We refer to the expansion of the basis set as "cloning," in analogy to the "spawning" procedure in AIMS. This synthesis of AIMS and MCE allows us to leverage the benefits of mean-field evolution during periods of strong nonadiabatic coupling while simultaneously avoiding mean-field artifacts in Ehrenfest dynamics. We explore the use of time-displaced basis sets, "trains," as a means of expanding the basis set for little cost. We also introduce a new bra-ket averaged Taylor expansion (BAT) to approximate the necessary potential energy and nonadiabatic coupling matrix elements. The BAT approximation avoids the necessity of computing electronic structure information at intermediate points between TBFs, as is usually done in saddle-point approximations used in AIMS. The efficiency of AIMC is demonstrated on the nonradiative decay of the first excited state of ethylene. The AIMC method has been implemented within the AIMS-MOLPRO package, which was extended to include Ehrenfest basis functions. PMID:25106573
Guiding ab initio calculations by alchemical derivatives.
to Baben, M; Achenbach, J O; von Lilienfeld, O A
2016-03-14
We assess the concept of alchemical transformations for predicting how a further and not-tested change in composition would change materials properties. This might help to guide ab initio calculations through multidimensional property-composition spaces. Equilibrium volumes, bulk moduli, and relative lattice stability of fcc and bcc 4d transition metals Zr, Nb, Mo, Tc, Ru, Rh, Pd, and Ag are calculated using density functional theory. Alchemical derivatives predict qualitative trends in lattice stability while equilibrium volumes and bulk moduli are predicted with less than 9% and 28% deviation, respectively. Predicted changes in equilibrium volume and bulk moduli for binary and ternary mixtures of Rh-Pd-Ag are in qualitative agreement even for predicted bulk modulus changes as large as +100% or -50%. Based on these results, it is suggested that alchemical transformations could be meaningful for enhanced sampling in the context of virtual high-throughput materials screening projects. PMID:26979677
Ab initio determination of light hadron masses.
Dürr, S; Fodor, Z; Frison, J; Hoelbling, C; Hoffmann, R; Katz, S D; Krieg, S; Kurth, T; Lellouch, L; Lippert, T; Szabo, K K; Vulvert, G
2008-11-21
More than 99% of the mass of the visible universe is made up of protons and neutrons. Both particles are much heavier than their quark and gluon constituents, and the Standard Model of particle physics should explain this difference. We present a full ab initio calculation of the masses of protons, neutrons, and other light hadrons, using lattice quantum chromodynamics. Pion masses down to 190 mega-electron volts are used to extrapolate to the physical point, with lattice sizes of approximately four times the inverse pion mass. Three lattice spacings are used for a continuum extrapolation. Our results completely agree with experimental observations and represent a quantitative confirmation of this aspect of the Standard Model with fully controlled uncertainties. PMID:19023076
Discovering chemistry with an ab initio nanoreactor
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-11-02
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 accelerated, first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor we show new pathways for glycine synthesis frommore » primitive compounds proposed to exist on the early Earth, providing new insight into the classic Urey-Miller experiment. Ultimately, these results highlight the emergence of theoretical and computational chemistry as a tool for discovery in addition to its traditional role of interpreting experimental findings.« less
Ab initio molar volumes and Gaussian radii.
Parsons, Drew F; Ninham, Barry W
2009-02-12
Ab initio molar volumes are calculated and used to derive radii for ions and neutral molecules using a spatially diffuse model of the electron distribution with Gaussian spread. The Gaussian radii obtained can be used for computation of nonelectrostatic ion-ion dispersion forces that underlie Hofmeister specific ion effects. Equivalent hard-sphere radii are also derived, and these are in reasonable agreement with crystalline ionic radii. The Born electrostatic self-energy is derived for a Gaussian model of the electronic charge distribution. It is shown that the ionic volumes used in electrostatic calculations of strongly hydrated cosmotropic ions ought best to include the first hydration shell. Ionic volumes for weakly hydrated chaotropic metal cations should exclude electron overlap (in electrostatic calculations). Spherical radii are calculated as well as nonisotropic ellipsoidal radii for nonspherical ions, via their nonisotropic static polarizability tensors. PMID:19140766
Discovering chemistry with an ab initio nanoreactor
NASA Astrophysics Data System (ADS)
Martinez, Todd
Traditional approaches for modeling chemical reaction networks such as those involved in combustion have focused on identifying individual reactions and using theoretical approaches to explore the underlying mechanisms. Recent advances involving graphical processing units (GPUs), commodity products developed for the videogaming industry, have made it possible to consider a distinct approach wherein one attempts to discover chemical reactions and mechanisms. We provide a brief summary of these developments and then discuss the concept behind the ``ab initio nanoreactor'' which explores the space of possible chemical reactions and molecular species for a given stoichiometry. The nanoreactor concept is exemplified with an example to the Urey-Miller reaction network which has been previously advanced as a potential model for prebiotic chemistry. We briefly discuss some of the future directions envisioned for the development of this nanoreactor concept.
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 accelerated, first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, providing new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery in addition to its traditional role of interpreting experimental findings. PMID:25411881
Guiding ab initio calculations by alchemical derivatives
NASA Astrophysics Data System (ADS)
to Baben, M.; Achenbach, J. O.; von Lilienfeld, O. A.
2016-03-01
We assess the concept of alchemical transformations for predicting how a further and not-tested change in composition would change materials properties. This might help to guide ab initio calculations through multidimensional property-composition spaces. Equilibrium volumes, bulk moduli, and relative lattice stability of fcc and bcc 4d transition metals Zr, Nb, Mo, Tc, Ru, Rh, Pd, and Ag are calculated using density functional theory. Alchemical derivatives predict qualitative trends in lattice stability while equilibrium volumes and bulk moduli are predicted with less than 9% and 28% deviation, respectively. Predicted changes in equilibrium volume and bulk moduli for binary and ternary mixtures of Rh-Pd-Ag are in qualitative agreement even for predicted bulk modulus changes as large as +100% or -50%. Based on these results, it is suggested that alchemical transformations could be meaningful for enhanced sampling in the context of virtual high-throughput materials screening projects.
Discovering chemistry with an ab initio nanoreactor.
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S; Martínez, Todd J
2014-12-01
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provide detailed physical insight. Although 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 accelerated first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor, we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, which provide new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery, in addition to its traditional role of interpreting experimental findings. PMID:25411881
Discovering chemistry with an ab initio nanoreactor
NASA Astrophysics Data System (ADS)
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-12-01
Chemical understanding is driven by the experimental discovery of new compounds and reactivity, and is supported by theory and computation that provide detailed physical insight. Although 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 accelerated first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor, we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, which provide new insight into the classic Urey-Miller experiment. These results highlight the emergence of theoretical and computational chemistry as a tool for discovery, in addition to its traditional role of interpreting experimental findings.
Ab Initio Modeling of Molecular Radiation
NASA Technical Reports Server (NTRS)
Jaffe, Richard; Schwenke, David
2014-01-01
Radiative emission from excited states of atoms and molecules can comprise a significant fraction of the total heat flux experienced by spacecraft during atmospheric entry at hypersonic speeds. For spacecraft with ablating heat shields, some of this radiative flux can be absorbed by molecular constituents in the boundary layer that are formed by the ablation process. Ab initio quantum mechanical calculations are carried out to predict the strengths of these emission and absorption processes. This talk will describe the methods used in these calculations using, as examples, the 4th positive emission bands of CO and the 1g+ 1u+ absorption in C3. The results of these calculations are being used as input to NASA radiation modeling codes like NeqAir, HARA and HyperRad.
Discovering chemistry with an ab initio nanoreactor
Wang, Lee-Ping; Titov, Alexey; McGibbon, Robert; Liu, Fang; Pande, Vijay S.; Martínez, Todd J.
2014-11-02
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 accelerated, first-principles molecular dynamics simulation of chemical reactions that discovers new molecules and mechanisms without preordained reaction coordinates or elementary steps. Using the nanoreactor we show new pathways for glycine synthesis from primitive compounds proposed to exist on the early Earth, providing new insight into the classic Urey-Miller experiment. Ultimately, these results highlight the emergence of theoretical and computational chemistry as a tool for discovery in addition to its traditional role of interpreting experimental findings.
Ab Initio Calculation of the Hoyle State
Epelbaum, Evgeny; Krebs, Hermann; Lee, Dean; Meissner, Ulf-G.
2011-05-13
The Hoyle state plays a crucial role in the helium burning of stars heavier than our Sun and in the production of carbon and other elements necessary for life. This excited state of the carbon-12 nucleus was postulated by Hoyle as a necessary ingredient for the fusion of three alpha particles to produce carbon at stellar temperatures. Although the Hoyle state was seen experimentally more than a half century ago nuclear theorists have not yet uncovered the nature of this state from first principles. In this Letter we report the first ab initio calculation of the low-lying states of carbon-12 using supercomputer lattice simulations and a theoretical framework known as effective field theory. In addition to the ground state and excited spin-2 state, we find a resonance at -85(3) MeV with all of the properties of the Hoyle state and in agreement with the experimentally observed energy.
NASA Astrophysics Data System (ADS)
Nattino, Francesco; Galparsoro, Oihana; Costanzo, Francesca; Díez Muiño, Ricardo; Alducin, Maite; Kroes, Geert-Jan
2016-06-01
Accurately modeling surface temperature and surface motion effects is necessary to study molecule-surface reactions in which the energy dissipation to surface phonons can largely affect the observables of interest. We present here a critical comparison of two methods that allow to model such effects, namely, the ab initio molecular dynamics (AIMD) method and the generalized Langevin oscillator (GLO) model, using the dissociation of N2 on W(110) as a benchmark. AIMD is highly accurate as the surface atoms are explicitly part of the dynamics, but this advantage comes with a large computational cost. The GLO model is much more computationally convenient, but accounts for lattice motion effects in a very approximate way. Results show that, despite its simplicity, the GLO model is able to capture the physics of the system to a large extent, returning dissociation probabilities which are in better agreement with AIMD than static-surface results. Furthermore, the GLO model and the AIMD method predict very similar energy transfer to the lattice degrees of freedom in the non-reactive events, and similar dissociation dynamics.
Ab initio alpha-alpha scattering
NASA Astrophysics Data System (ADS)
Elhatisari, Serdar; Lee, Dean; Rupak, Gautam; Epelbaum, Evgeny; Krebs, Hermann; Lähde, Timo A.; Luu, Thomas; Meißner, Ulf-G.
2015-12-01
Processes such as the scattering of alpha particles (4He), the triple-alpha reaction, and alpha capture play a major role in stellar nucleosynthesis. In particular, alpha capture on carbon determines the ratio of carbon to oxygen during helium burning, and affects subsequent carbon, neon, oxygen, and silicon burning stages. It also substantially affects models of thermonuclear type Ia supernovae, owing to carbon detonation in accreting carbon-oxygen white-dwarf stars. In these reactions, the accurate calculation of the elastic scattering of alpha particles and alpha-like nuclei—nuclei with even and equal numbers of protons and neutrons—is important for understanding background and resonant scattering contributions. First-principles calculations of processes involving alpha particles and alpha-like nuclei have so far been impractical, owing to the exponential growth of the number of computational operations with the number of particles. Here we describe an ab initio calculation of alpha-alpha scattering that uses lattice Monte Carlo simulations. We use lattice effective field theory to describe the low-energy interactions of protons and neutrons, and apply a technique called the ‘adiabatic projection method’ to reduce the eight-body system to a two-cluster system. We take advantage of the computational efficiency and the more favourable scaling with system size of auxiliary-field Monte Carlo simulations to compute an ab initio effective Hamiltonian for the two clusters. We find promising agreement between lattice results and experimental phase shifts for s-wave and d-wave scattering. The approximately quadratic scaling of computational operations with particle number suggests that it should be possible to compute alpha scattering and capture on carbon and oxygen in the near future. The methods described here can be applied to ultracold atomic few-body systems as well as to hadronic systems using lattice quantum chromodynamics to describe the interactions of
Ab initio alpha-alpha scattering.
Elhatisari, Serdar; Lee, Dean; Rupak, Gautam; Epelbaum, Evgeny; Krebs, Hermann; Lähde, Timo A; Luu, Thomas; Meißner, Ulf-G
2015-12-01
Processes such as the scattering of alpha particles ((4)He), the triple-alpha reaction, and alpha capture play a major role in stellar nucleosynthesis. In particular, alpha capture on carbon determines the ratio of carbon to oxygen during helium burning, and affects subsequent carbon, neon, oxygen, and silicon burning stages. It also substantially affects models of thermonuclear type Ia supernovae, owing to carbon detonation in accreting carbon-oxygen white-dwarf stars. In these reactions, the accurate calculation of the elastic scattering of alpha particles and alpha-like nuclei--nuclei with even and equal numbers of protons and neutrons--is important for understanding background and resonant scattering contributions. First-principles calculations of processes involving alpha particles and alpha-like nuclei have so far been impractical, owing to the exponential growth of the number of computational operations with the number of particles. Here we describe an ab initio calculation of alpha-alpha scattering that uses lattice Monte Carlo simulations. We use lattice effective field theory to describe the low-energy interactions of protons and neutrons, and apply a technique called the 'adiabatic projection method' to reduce the eight-body system to a two-cluster system. We take advantage of the computational efficiency and the more favourable scaling with system size of auxiliary-field Monte Carlo simulations to compute an ab initio effective Hamiltonian for the two clusters. We find promising agreement between lattice results and experimental phase shifts for s-wave and d-wave scattering. The approximately quadratic scaling of computational operations with particle number suggests that it should be possible to compute alpha scattering and capture on carbon and oxygen in the near future. The methods described here can be applied to ultracold atomic few-body systems as well as to hadronic systems using lattice quantum chromodynamics to describe the interactions of
Comas-Vives, Aleix
2016-03-14
In this contribution, realistic amorphous SiO2 models of 2.1 × 2.1 nm with silanol densities ranging 1.1-7.2 OH per nm(2) are obtained by means of ab initio calculations via the dehydroxylation of a fully hydroxylated silica surface. The dehydroxyation process is considered to take place via direct condensation of adjacent silanol groups and silica migration steps. The latter reconstructions are needed in order to obtain highly dehydroxylated silica surfaces with favorable energetics and without the formation of defects. The obtained surface phase diagram of different silica models as a function of temperature and PH2O is able to correctly describe the silanol density under different conditions, and the IR spectroscopic signatures of the silanols are in qualitative agreement with the experiment. The amorphous silica models presented here have a high degree of heterogeneity as found from the big variability obtained in the energetics of the dehydroxylation steps. It was also found that the resulting average Si-O distance of the newly formed siloxane bridges serves as a descriptor of the strain introduced in the silica surface. All these factors can be crucial in order to simulate the activity of catalysts grafted onto silica with different silanol densities, especially the one containing ca. 1 OH per nm(2), which can serve as a model for the SiO2 surface pretreated under high vacuum and at 700 °C. PMID:26898649
Zimmermann, Tomáš; Vaníček, Jiří
2014-10-07
We derive a somewhat crude, yet very efficient semiclassical approximation for computing nonadiabatic spectra. The resulting method, which is a generalization of the multiple-surface dephasing representation, includes quantum effects through interference of mixed quantum-classical trajectories and through quantum treatment of the collective electronic degree of freedom. The method requires very little computational effort beyond the fewest-switches surface hopping or Ehrenfest locally mean-field dynamics and is very easy to implement. The proposed approximation is tested by computing the absorption and time-resolved stimulated emission spectra of pyrazine using the four-dimensional three-surface model which allows for comparison with the numerically exact quantum spectra. As expected, the multiple-surface dephasing representation is not suitable for high-resolution linear spectra, yet it seems to capture all the important features of pump-probe spectra. Finally, the method is combined with on-the-fly ab initio evaluation of the electronic structure (i.e., energies, forces, electric-dipole, and nonadiabatic couplings) in order to compute fully dimensional nonadiabatic spectra of pyrazine without approximations inherent to analytical, including vibronic-coupling models. The Appendix provides derivations of perturbative expressions for linear and pump-probe spectra of arbitrary mixed states and for arbitrary laser pulse shapes.
Ab initio derivation of model energy density functionals
NASA Astrophysics Data System (ADS)
Dobaczewski, Jacek
2016-08-01
I propose a simple and manageable method that allows for deriving coupling constants of model energy density functionals (EDFs) directly from ab initio calculations performed for finite fermion systems. A proof-of-principle application allows for linking properties of finite nuclei, determined by using the nuclear nonlocal Gogny functional, to the coupling constants of the quasilocal Skyrme functional. The method does not rely on properties of infinite fermion systems but on the ab initio calculations in finite systems. It also allows for quantifying merits of different model EDFs in describing the ab initio results.
The implementation of ab initio quantum chemistry calculations on transporters.
Cooper, M D; Hillier, I H
1991-06-01
The RHF and geometry optimization sections of the ab initio quantum chemistry code, GAMESS, have been optimized for a network of parallel microprocessors, Inmos T800-20 transputers, using both indirect and direct SCF techniques. The results indicate great scope for implementation of such codes on small parallel computer systems, very high efficiencies having been achieved, particularly in the cases of direct SCF and geometry optimization with large basis sets. The work, although performed upon one particular parallel system, the Meiko Computing Surface, is applicable to a wide range of parallel systems with both shared and distributed memory. PMID:1919615
Li, Y Q; Ma, F C; Sun, M T
2013-10-21
A full three-dimensional global potential energy surface is reported first time for the title system, which is important for the photodissociation processes. It is obtained using double many-body expansion theory and an extensive set of accurate ab initio energies extrapolated to the complete basis set limit. Such a work can be recommended for dynamics studies of the N((2)D) + H2 reaction, a reliable theoretical treatment of the photodissociation dynamics and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen containing systems. In turn, a preliminary theoretical study of the reaction N((2)D)+H2(X(1)Σg (+))(ν=0,j=0)→NH(a(1)Δ)+H((2)S) has been carried out with the method of quasi-classical trajectory on the new potential energy surface. Integral cross sections and thermal rate constants have been calculated, providing perhaps the most reliable estimate of the integral cross sections and the rate constants known thus far for such a reaction. PMID:24160511
NASA Astrophysics Data System (ADS)
Almeida, A. L.; Martins, João B. L.; Taft, C. A.; Longo, E.; Lester, W. A.
1998-09-01
Ab initio and semiempirical calculations of large cluster models have been performed in order to study water adsorption and dissociation on pure, defective (vacancies) and doped (Li, Na, K, Ca, Fe) MgO (001) surfaces. The geometries of the adsorbed and dissociated molecules have been optimized preparatory to analysis of binding energies, stretching frequencies, charge transfers, preferential sites of interaction, and bond distances. We have used Mulliken, natural bond order, and electrostatic-derived atomic and overlap populations to analyze charge distributions in the clusters. We have also investigated transition structures, activation energies, energy gaps, HOMO, density of states, SCF orbital energies as well as the acid-base properties of our cluster model. Numerical results are compared, where possible, with experiment, interpreted in the framework of various analytical models, and correlated with site coordination numbers, corner and edge site preferential locations, and direction of charge transfer. A thorough charge analysis indicates substantial charge redistribution in the magnesium oxide crystal as a result of water adsorption and dissociation in pure, defective, and doped MgO crystals. The introduction of heavier impurities and vacancies could produce substantial changes in the physical and chemical properties of the catalyst and increase the binding and dissociation energies. Some of the largest changes originate from the introduction of vacancies. Two and three-dimensional potential energy surfaces are used to investigate activation energies of hydroxylation on the MgO surface. Stretching frequencies are correlated with magnesium and oxygen coordination numbers.
Brüssel, Marc; di Dio, Philipp J.; Muñiz, Kilian; Kirchner, Barbara
2011-01-01
We carried out ab initio molecular dynamic simulations in order to determine the free energy surfaces of two selected reactions including solvents, namely a rearrangement of a ruthenium oxoester in water and a carbon dioxide addition to a palladium complex in carbon dioxide. For the latter reaction we also investigated the gas phase reaction in order to take solvent effects into account. We used two techniques to reconstruct the free energy surfaces: thermodynamic integration and metadynamics. Furthermore, we gave a reasonable error estimation of the computed free energy surface. We calculated a reaction barrier of ΔF = 59.5 ± 8.5 kJ mol−1 for the rearrangement of a ruthenium oxoester in water from thermodynamic integration. For the carbon dioxide addition to the palladium complex in carbon dioxide we found a ΔF = 44.9 ± 3.3 kJ mol−1 from metadynamics simulations with one collective variable. The investigation of the same reactions in the gas phase resulted in ΔF = 24.9 ± 6.7 kJ mol−1 from thermodynamic integration, in ΔF = 26.7 ± 2.3 kJ mol−1 from metadynamics simulations with one collective variable, and in ΔF = 27.1 ± 5.9 kJ mol−1 from metadynamics simulations with two collective variables. PMID:21541065
NASA Astrophysics Data System (ADS)
Li, Y. Q.; Ma, F. C.; Sun, M. T.
2013-10-01
A full three-dimensional global potential energy surface is reported first time for the title system, which is important for the photodissociation processes. It is obtained using double many-body expansion theory and an extensive set of accurate ab initio energies extrapolated to the complete basis set limit. Such a work can be recommended for dynamics studies of the N(2D) + H2 reaction, a reliable theoretical treatment of the photodissociation dynamics and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen containing systems. In turn, a preliminary theoretical study of the reaction N(^2D)+H_2(X^1Σ _g^+)(ν =0,j=0)rArr NH(a^1Δ )+H(^2S) has been carried out with the method of quasi-classical trajectory on the new potential energy surface. Integral cross sections and thermal rate constants have been calculated, providing perhaps the most reliable estimate of the integral cross sections and the rate constants known thus far for such a reaction.
Kempisty, Paweł
2014-11-15
Adsorption of ammonia at NH{sub 3}/NH{sub 2}/H-covered GaN(0001) surface was analyzed using results of ab initio calculations. The whole configuration space of partially NH{sub 3}/NH{sub 2}/H-covered GaN(0001) surface was divided into zones of differently pinned Fermi level: at the Ga broken bond state for dominantly bare surface (region I), at the valence band maximum (VBM) for NH{sub 2} and H-covered surface (region II), and at the conduction band minimum (CBM) for NH{sub 3}-covered surface (region III). The electron counting rule (ECR) extension was formulated for the case of adsorbed molecules. The extensive ab intio calculations show the validity of the ECR in case of all mixed H-NH{sub 2}-NH{sub 3} coverages for the determination of the borders between the three regions. The adsorption was analyzed using the recently identified dependence of the adsorption energy on the charge transfer at the surface. For region I ammonia adsorbs dissociatively, disintegrating into a H adatom and a HN{sub 2} radical for a large fraction of vacant sites, while for region II adsorption of ammonia is molecular. The dissociative adsorption energy strongly depends on the Fermi level at the surface (pinned) and in the bulk (unpinned) while the molecular adsorption energy is determined by bonding to surface only, in accordance to the recently published theory. Adsorption of Ammonia in region III (Fermi level pinned at CBM) leads to an unstable configuration both molecular and dissociative, which is explained by the fact that broken Ga-bonds are doubly occupied by electrons. The adsorbing ammonia brings 8 electrons to the surface, necessitating the transfer of these two electrons from the Ga broken bond state to the Fermi level. This is an energetically costly process. Adsorption of ammonia at H-covered site leads to the creation of a NH{sub 2} radical at the surface and escape of H{sub 2} molecule. The process energy is close to 0.12 eV, thus not large, but the direct inverse
Ab Initio: And a New Era of Airline Pilot Training.
ERIC Educational Resources Information Center
Gesell, Laurence E.
1995-01-01
Expansion of air transportation and decreasing numbers seeking pilot training point to a shortage of qualified pilots. Ab initio training, in which candidates with no flight time are trained to air transport proficiency, could resolve the problem. (SK)
An ab initio-based approach to the stability of GaN(0 0 0 1) surfaces under Ga-rich conditions
NASA Astrophysics Data System (ADS)
Ito, Tomonori; Akiyama, Toru; Nakamura, Kohji
2009-05-01
Structural stability of GaN(0 0 0 1) under Ga-rich conditions is systematically investigated by using our ab initio-based approach. The surface phase diagram for GaN(0 0 0 1) including (2×2) and pseudo-(1×1) is obtained as functions of temperature and Ga beam equivalent pressure by comparing chemical potentials of Ga atom in the gas phase with that on the surface. The calculated results reveal that the pseudo-(1×1) appearing below 684-973 K changes its structure to the (2×2) with Ga adatom at higher temperatures beyond 767-1078 K via the newly found (1×1) with two adlayers of Ga. These results are consistent with the stable temperature range of both the pseudo-(1×1) and (2×2) with Ga adatom obtained experimentally. Furthermore, it should be noted that the structure with another coverage of Ga adatoms between the (1×1) and (2×2)-Ga does not appear as a stable structure of GaN(0 0 0 1). Furthermore, ghost island formation observed by scanning tunneling microscopy is discussed on the basis of the phase diagram.
Denis-Alpizar, Otoniel; Kalugina, Yulia; Stoecklin, Thierry; Vera, Mario Hernández; Lique, François
2013-12-14
We present a new four-dimensional potential energy surface for the collisional excitation of HCN by H{sub 2}. Ab initio calculations of the HCN–H{sub 2} van der Waals complex, considering both molecules as rigid rotors, were carried out at the explicitly correlated coupled cluster with single, double, and perturbative triple excitations [CCSD(T)-F12a] level of theory using an augmented correlation-consistent triple zeta (aVTZ) basis set. The equilibrium structure is linear HCN–H{sub 2} with the nitrogen pointing towards H{sub 2} at an intermolecular separation of 7.20 a{sub 0}. The corresponding well depth is −195.20 cm{sup −1}. A secondary minimum of −183.59 cm{sup −1} was found for a T-shape configuration with the H of HCN pointing to the center of mass of H{sub 2}. We also determine the rovibrational energy levels of the HCN–para-H{sub 2} and HCN–ortho-H{sub 2} complexes. The calculated dissociation energies for the para and ortho complexes are 37.79 cm{sup −1} and 60.26 cm{sup −1}, respectively. The calculated ro-vibrational transitions in the HCN–H{sub 2} complex are found to agree by more than 0.5% with the available experimental data, confirming the accuracy of the potential energy surface.
Ab initio two-component Ehrenfest dynamics
NASA Astrophysics Data System (ADS)
Ding, Feizhi; Goings, Joshua J.; Liu, Hongbin; Lingerfelt, David B.; Li, Xiaosong
2015-09-01
We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H2 and O2. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.
Ab Initio Studies of Calcium Carbonate Hydration.
Lopez-Berganza, Josue A; Diao, Yijue; Pamidighantam, Sudhakar; Espinosa-Marzal, Rosa M
2015-11-25
Ab initio simulations of large hydrated calcium carbonate clusters are challenging due to the existence of multiple local energy minima. Extensive conformational searches around hydrated calcium carbonate clusters (CaCO3·nH2O for n = 1-18) were performed to find low-energy hydration structures using an efficient combination of Monte Carlo searches, density-functional tight binding (DFTB+) method, and density-functional theory (DFT) at the B3LYP level, or Møller-Plesset perturbation theory at the MP2 level. This multilevel optimization yields several low-energy structures for hydrated calcium carbonate. Structural and energetics analysis of the hydration of these clusters revealed a first hydration shell composed of 12 water molecules. Bond-length and charge densities were also determined for different cluster sizes. The solvation of calcium carbonate in bulk water was investigated by placing the explicitly solvated CaCO3·nH2O clusters in a polarizable continuum model (PCM). The findings of this study provide new insights into the energetics and structure of hydrated calcium carbonate and contribute to the understanding of mechanisms where calcium carbonate formation or dissolution is of relevance. PMID:26505205
Ab initio two-component Ehrenfest dynamics
Ding, Feizhi; Goings, Joshua J.; Liu, Hongbin; Lingerfelt, David B.; Li, Xiaosong
2015-09-21
We present an ab initio two-component Ehrenfest-based mixed quantum/classical molecular dynamics method to describe the effect of nuclear motion on the electron spin dynamics (and vice versa) in molecular systems. The two-component time-dependent non-collinear density functional theory is used for the propagation of spin-polarized electrons while the nuclei are treated classically. We use a three-time-step algorithm for the numerical integration of the coupled equations of motion, namely, the velocity Verlet for nuclear motion, the nuclear-position-dependent midpoint Fock update, and the modified midpoint and unitary transformation method for electronic propagation. As a test case, the method is applied to the dissociation of H{sub 2} and O{sub 2}. In contrast to conventional Ehrenfest dynamics, this two-component approach provides a first principles description of the dynamics of non-collinear (e.g., spin-frustrated) magnetic materials, as well as the proper description of spin-state crossover, spin-rotation, and spin-flip dynamics by relaxing the constraint on spin configuration. This method also holds potential for applications to spin transport in molecular or even nanoscale magnetic devices.
Nickolaisen, S.L.; Roehl, C.M.; Blakeley, L.K.; Friedl, R.R.; Francisco, J.S.; Liu, R.; Sander, S.P.
2000-01-20
The kinetics of the HO{sub 2} + ClO reaction was studied using the flash photolysis/ultraviolet absorption technique over the temperature range 203--364 K and pressure range 50--700 Torr of N{sub 2}. In contrast to previous work, the temperature dependence displayed linear Arrhenius behavior over the entire temperature range. Ab initio calculations of intermediates and transition states have been carried out on the singlet and triplet potential energy surfaces. These calculations show that the reaction proceeds mainly through the ClO-HO{sub 2} complex on the triplet surface; however, collisionally stabilized HOOOCl formed on the singlet surface will possess an appreciable lifetime due to large barriers toward decomposition to HCl and HOCl. Termolecular rate calculations using ab initio parameters lead to a strong collision rate constant for HOOOCl formation. This intermediate may be important under both laboratory and atmospheric conditions.
Le, Hung M; Dinh, Thach S; Le, Hieu V
2011-10-13
The singlet-triplet transformation and molecular dissociation of ozone (O(3)) gas is investigated by performing quasi-classical molecular dynamics (MD) simulations on an ab initio potential energy surface (PES) with visible and near-infrared excitations. MP4(SDQ) level of theory with the 6-311g(2d,2p) basis set is executed for three different electronic spin states (singlet, triplet, and quintet). In order to simplify the potential energy function, an approximation is adopted by ignoring the spin-orbit coupling and allowing the molecule to switch favorably and instantaneously to the spin state that is more energetically stable (lowest in energy among the three spin states). This assumption has previously been utilized to study the SiO(2) system as reported by Agrawal et al. (J. Chem. Phys. 2006, 124 (13), 134306). The use of such assumption in this study probably makes the upper limits of computed rate coefficients the true rate coefficients. The global PES for ozone is constructed by fitting 5906 ab initio data points using a 60-neuron two-layer feed-forward neural network. The mean-absolute error and root-mean-squared error of this fit are 0.0446 eV (1.03 kcal/mol) and 0.0756 eV (1.74 kcal/mol), respectively, which reveal very good fitting accuracy. The parameter coefficients of the global PES are reported in this paper. In order to identify the spin state with high confidence, we propose the use of a pattern-recognition neural network, which is trained to predict the spin state of a given configuration (with a prediction accuracy being 95.6% on a set of testing data points). To enhance the prediction effectiveness, a buffer series of five points are validated to confirm the spin state during the MD process to gain better confidence. Quasi-classical MD simulations from 1.2 to 2.4 eV of total internal energy (including zero-point energy) result in rate coefficients of singlet-triplet transformation in the range of 0.027 ps(-1) to 1.21 ps(-1). Also, we find very
NASA Astrophysics Data System (ADS)
Czakó, Gábor
2013-04-01
Chemically accurate full-dimensional non-spin-orbit and spin-orbit (SO) ground-state potential energy surfaces (PESs) are obtained for the Br + CH4 → HBr + CH3 reaction by fitting 21 574 composite ab initio energy points. The composite method considers electron correlation methods up to CCSD(T), basis sets up to aug-cc-pwCVTZ-PP, correlation of the core electrons, scalar relativistic effects via an effective core potential (ECP), and SO corrections, thereby achieving an accuracy better than 0.5 kcal/mol. Benchmark structures and relative energies are computed for the stationary points using the ab initio focal-point analysis (FPA) scheme based on both ECP and Douglas-Kroll approaches providing all-electron relativistic CCSDT(Q)/complete-basis-set quality energies. The PESs accurately describe the saddle point of the abstraction reaction and the van der Waals complexes in the entrance and product channels. The SO-corrected PES provides a classical barrier height of 7285(7232 ± 50) cm-1, De values of 867(799 ± 10) and 399(344 ± 10) cm-1 for the complexes CH3-HBr and CH3-BrH, respectively, and reaction endothermicity of 7867(7857 ± 50) cm-1, in excellent agreement with the new, FPA-based benchmark data shown in parentheses. The difference between the Br + CH4 asymptotes of the non-SO and SO PESs is 1240 cm-1, in good agreement with the experiment (1228 cm-1). Quasiclassical trajectory calculations based on more than 13 million trajectories for the late-barrier Br + CH4(vk = 0, 1) [k = 1, 2, 3, 4] reactions show that the vibrational energy, especially the excitation of the stretching modes, activates the reaction much more efficiently than translational energy, in agreement with the extended Polanyi rules. Angular distributions show dominant backward scattering for the ground-state reaction and forward scattering for the stretching-excited reactions. The reactivity on the non-SO PES is about 3-5 times larger than that on the SO PES in a wide collision energy
Varandas, A J C
2013-04-01
The Combined-Hyperbolic-Inverse-Power-Representation method, which treats evenly both short- and long-range interactions, is used to fit an extensive set of ab initio points for HO2 previously utilized [Xu et al., J. Chem. Phys. 122, 244305 (2005)] to develop a spline interpolant. The novel form is shown to perform accurately when compared with others, while quasiclassical trajectory calculations of the O + OH reaction clearly pinpoint the role of long-range forces at low temperatures. PMID:23574218
Chen, Liuyang; Shao, Kejie; Chen, Jun; Yang, Minghui; Zhang, Dong H
2016-05-21
This work performs a time-dependent wavepacket study of the H2 + C2H → H + C2H2 reaction on a new ab initio potential energy surface (PES). The PES is constructed using neural network method based on 68 478 geometries with energies calculated at UCCSD(T)-F12a/aug-cc-pVTZ level and covers H2 + C2H↔H + C2H2, H + C2H2 → HCCH2, and HCCH2 radial isomerization reaction regions. The reaction dynamics of H2 + C2H → H + C2H2 are investigated using full-dimensional quantum dynamics method. The initial-state selected reaction probabilities are calculated for reactants in eight vibrational states. The calculated results showed that the H2 vibrational excitation predominantly enhances the reactivity while the excitation of bending mode of C2H slightly inhibits the reaction. The excitations of two stretching modes of C2H molecule have negligible effect on the reactivity. The integral cross section is calculated with J-shift approximation and the mode selectivity in this reaction is discussed. The rate constants over 200-2000 K are calculated and agree well with the experimental measured values. PMID:27208951
Song, Hongwei; Lu, Yunpeng; Li, Jun; Yang, Minghui; Guo, Hua
2016-04-28
An initial state selected time-dependent wave packet method is applied to study the dynamics of the OH + CHD3 reaction with a six-dimensional model on a newly developed full-dimensional ab initio potential energy surface (PES). This quantum dynamical (QD) study is complemented by full-dimensional quasi-classical trajectory (QCT) calculations on the same PES. The QD results indicate that both translational energy and the excitation of the CH stretching mode significantly promote the reaction while the excitation of the umbrella mode has a negligible effect on the reactivity. For this early barrier reaction, interestingly, the CH stretching mode is more effective than translational energy in promoting the reaction except at very low collision energies. These QD observations are supported by QCT results. The higher efficacy of the CH stretching model in promoting this early barrier reaction is inconsistent with the prediction of the naively extended Polanyi's rules, but can be rationalized by the recently proposed sudden vector projection model. PMID:27131546
NASA Astrophysics Data System (ADS)
Song, Hongwei; Lu, Yunpeng; Li, Jun; Yang, Minghui; Guo, Hua
2016-04-01
An initial state selected time-dependent wave packet method is applied to study the dynamics of the OH + CHD3 reaction with a six-dimensional model on a newly developed full-dimensional ab initio potential energy surface (PES). This quantum dynamical (QD) study is complemented by full-dimensional quasi-classical trajectory (QCT) calculations on the same PES. The QD results indicate that both translational energy and the excitation of the CH stretching mode significantly promote the reaction while the excitation of the umbrella mode has a negligible effect on the reactivity. For this early barrier reaction, interestingly, the CH stretching mode is more effective than translational energy in promoting the reaction except at very low collision energies. These QD observations are supported by QCT results. The higher efficacy of the CH stretching model in promoting this early barrier reaction is inconsistent with the prediction of the naively extended Polanyi's rules, but can be rationalized by the recently proposed sudden vector projection model.
Li, Jun; Jiang, Bin; Song, Hongwei; Ma, Jianyi; Zhao, Bin; Dawes, Richard; Guo, Hua
2015-05-21
We survey the recent advances in theoretical understanding of quantum state resolved dynamics, using the title reactions as examples. It is shown that the progress was made possible by major developments in two areas. First, an accurate analytical representation of many high-level ab initio points over a large configuration space can now be made with high fidelity and the necessary permutation symmetry. The resulting full-dimensional global potential energy surfaces enable dynamical calculations using either quasi-classical trajectory or more importantly quantum mechanical methods. The second advance is the development of accurate and efficient quantum dynamical methods, which are necessary for providing a reliable treatment of quantum effects in reaction dynamics such as tunneling, resonances, and zero-point energy. The powerful combination of the two advances has allowed us to achieve a quantitatively accurate characterization of the reaction dynamics, which unveiled rich dynamical features such as steric steering, strong mode specificity, and bond selectivity. The dependence of reactivity on reactant modes can be rationalized by the recently proposed sudden vector projection model, which attributes the mode specificity and bond selectivity to the coupling of reactant modes with the reaction coordinate at the relevant transition state. The deeper insights provided by these theoretical studies have advanced our understanding of reaction dynamics to a new level. PMID:25886142
NASA Astrophysics Data System (ADS)
Chen, Liuyang; Shao, Kejie; Chen, Jun; Yang, Minghui; Zhang, Dong H.
2016-05-01
This work performs a time-dependent wavepacket study of the H2 + C2H → H + C2H2 reaction on a new ab initio potential energy surface (PES). The PES is constructed using neural network method based on 68 478 geometries with energies calculated at UCCSD(T)-F12a/aug-cc-pVTZ level and covers H2 + C2H↔H + C2H2, H + C2H2 → HCCH2, and HCCH2 radial isomerization reaction regions. The reaction dynamics of H2 + C2H → H + C2H2 are investigated using full-dimensional quantum dynamics method. The initial-state selected reaction probabilities are calculated for reactants in eight vibrational states. The calculated results showed that the H2 vibrational excitation predominantly enhances the reactivity while the excitation of bending mode of C2H slightly inhibits the reaction. The excitations of two stretching modes of C2H molecule have negligible effect on the reactivity. The integral cross section is calculated with J-shift approximation and the mode selectivity in this reaction is discussed. The rate constants over 200-2000 K are calculated and agree well with the experimental measured values.
Role of chlorine in the dynamics of GaN(0 0 0 1) surface during HVPE GaN growth—Ab initio study
NASA Astrophysics Data System (ADS)
Krukowski, Stanisław; Kempisty, Paweł; Strąk, Paweł
2008-04-01
Atomistic processes, involving chlorine, occurring at GaN(0 0 0 1) surface during hydride vapor-phase epitaxy (HVPE) growth, were modeled, using ab-initio quantum mechanical density functional theory (QM DFT) SIESTA code. The process of adsorption of hydrogen chloride (HCl) at the surface was simulated. Both Ga-rich and N-rich surface conditions were considered. For both cases it was assumed that hydrogen is abundant, which leads to H and NH 2 coverage for Ga-rich and N-rich state of the surface, respectively. By modeling a HCl molecule, approaching the surface, it was shown that two states exist: physisorbed at far distance, and chemisorbed at closer location. For both the cases, i.e. N-rich and Ga-rich, the physical adsorption is barrierless. The transition to chemically adsorbed state encounters an energy barrier that severely depends on the state of the GaN(0 0 0 1) surface. In the case of the Ga-rich surface, the barrier is relatively small, of the order of 0.2 eV. In the case of the N-rich surface, a high-energy barrier exists, close to 2 eV. At the final i.e. chemisorbed state, the molecule disintegrates, leaving chlorine atom strongly attached to the topmost Ga atoms; the hydrogen atom is detached and paired with the other hydrogen atom. The adsorption of HCl molecule approaching Ga-H radical, located on top of NH 2 coverage, occurs via physic to chemi-sorbed state. The transition to chemically adsorbed state requires surmounting the energy barrier of about 0.6 eV: free H 2 molecule is left and the Ga-Cl admolecule remains at the surface. As before, H coverage is ineffective in protection of gallium against HCl attack. The desorption of chlorine atom from the GaCl admolecule, attached at the GaN(0 0 0 1) surface, was also analyzed. Single Cl atom desorption encounters a high-energy barrier, of the order of 4.5 eV. The desorption of HCl, involving the additional H atom, needs to overcome a significantly lower-energy barrier, of the order of 0.6 eV, creating
Ab initio water pair potential with flexible monomers.
Jankowski, Piotr; Murdachaew, Garold; Bukowski, Robert; Akin-Ojo, Omololu; Leforestier, Claude; Szalewicz, Krzysztof
2015-03-26
A potential energy surface for the water dimer with explicit dependence on monomer coordinates is presented. The surface was fitted to a set of previously published interaction energies computed on a grid of over a quarter million points in the 12-dimensional configurational space using symmetry-adapted perturbation theory and coupled-cluster methods. The present fit removes small errors in published fits, and its accuracy is critically evaluated. The minimum and saddle-point structures of the potential surface were found to be very close to predictions from direct ab initio optimizations. The computed second virial coefficients agreed well with experimental values. At low temperatures, the effects of monomer flexibility in the virial coefficients were found to be much smaller than the quantum effects. PMID:25687650
NASA Astrophysics Data System (ADS)
Mancini, John S.; Bowman, Joel M.
2013-03-01
We report a global, full-dimensional, ab initio potential energy surface describing the HCl-H2O dimer. The potential is constructed from a permutationally invariant fit, using Morse-like variables, to over 44 000 CCSD(T)-F12b/aug-cc-pVTZ energies. The surface describes the complex and dissociated monomers with a total RMS fitting error of 24 cm-1. The normal modes of the minima, low-energy saddle point and separated monomers, the double minimum isomerization pathway and electronic dissociation energy are accurately described by the surface. Rigorous quantum mechanical diffusion Monte Carlo (DMC) calculations are performed to determine the zero-point energy and wavefunction of the complex and the separated fragments. The calculated zero-point energies together with a De value calculated from CCSD(T) with a complete basis set extrapolation gives a D0 value of 1348 ± 3 cm-1, in good agreement with the recent experimentally reported value of 1334 ± 10 cm-1 [B. E. Casterline, A. K. Mollner, L. C. Ch'ng, and H. Reisler, J. Phys. Chem. A 114, 9774 (2010), 10.1021/jp102532m]. Examination of the DMC wavefunction allows for confident characterization of the zero-point geometry to be dominant at the C2v double-well saddle point and not the Cs global minimum. Additional support for the delocalized zero-point geometry is given by numerical solutions to the 1D Schrödinger equation along the imaginary-frequency out-of-plane bending mode, where the zero-point energy is calculated to be 52 cm-1 above the isomerization barrier. The D0 of the fully deuterated isotopologue is calculated to be 1476 ± 3 cm-1, which we hope will stand as a benchmark for future experimental work.
Long, Run; Prezhdo, Oleg V
2011-11-30
Following recent experiments [Science 2010, 328, 1543; PNAS 2011, 108, 965], we report an ab initio nonadiabatic molecular dynamics (NAMD) simulation of the ultrafast photoinduced electron transfer (ET) from a PbSe quantum dot (QD) into the rutile TiO(2) (110) surface. The system forms the basis for QD-sensitized semiconductor solar cells and demonstrates that ultrafast interfacial ET is instrumental for achieving high efficiencies in solar-to-electrical energy conversion. The simulation supports the observation that the ET successfully competes with energy losses due to electron-phonon relaxation. The ET proceeds by the adiabatic mechanism because of strong donor-acceptor coupling. High frequency polar vibrations of both QD and TiO(2) promote the ET, since these modes can rapidly influence the donor-acceptor state energies and coupling. Low frequency vibrations generate a distribution of initial conditions for ET, which shows a broad variety of scenarios at the single-molecule level. Compared to the molecule-TiO(2) interfaces, the QD-TiO(2) system exhibits pronounced differences that arise due to the larger size and higher rigidity of QDs relative to molecules. Both donor and acceptor states are more delocalized in the QD system, and the ET is promoted by optical phonons, which have relatively low frequencies in the QD materials composed of heavy elements. In contrast, in molecular systems, optical phonons are not thermally accessible under ambient conditions. Meanwhile, TiO(2) acceptor states resemble surface impurities due to the local influence of molecular chromophores. At the same time, the photoinduced ET at both QD-TiO(2) and molecule-TiO(2) interfaces is ultrafast and occurs by the adiabatic mechanism, as a result of strong donor-acceptor coupling. The reported state-of-the-art simulation generates a detailed time-domain atomistic description of the interfacial ET process that is fundamental to a wide variety of applications. PMID:22007727
Adsorption of ethanol on V2O5 (010) surface for gas-sensing applications: Ab initio investigation
NASA Astrophysics Data System (ADS)
Qin, Yuxiang; Cui, Mengyang; Ye, Zhenhua
2016-08-01
The adsorption of ethanol on V2O5 (010) surface was investigated by means of density functional theory (DFT) with a combined generalized gradient approximation (GGA) plus Hubbard U approach to exploit the potential sensing applications. The adsorption configurations were first constructed by considering different orientations of ethanol molecule to V and O sites on the "Hill"- and "Valley"-like regions of corrugated (010) surface. It is found that ethanol molecule can adsorb on whole surface in multiple stable configurations. Nevertheless the molecular adsorption on the "Hill"-like surface is calculated to occur preferentially, and the single coordinated oxygen on "Hill"-like surface (O1(H)) acting as the most energetically favorable adsorption site shows the strongest adsorption ability to ethanol molecule. Surface adsorption of ethanol tunes the electronic structure of V2O5 and cause an n-doping effect. As a consequence, the Fermi levels shift toward the conductive bond increasing the charge carrier concentration of electrons in adsorbed V2O5. The sensitive electronic structure and the multiple stable configurations to ethanol adsorption highlight the high adsorption activity and then the potential of V2O5 (010) surface applied to high sensitive sensor for ethanol vapor detection. Further Mulliken population and Natural bond orbital (NBO) calculations quantify the electron transfer from the adsorbed ethanol to the surface, and correlates the adsorption ability of surface sites with the charge donation and dispersion.
NASA Astrophysics Data System (ADS)
Dholabhai, Pratik P.; Atta-Fynn, Raymond; Ray, Asok K.
2008-12-01
In our continuing attempts to understand theoretically various surface properties such as corrosion and potential catalytic activity of actinide surfaces in the presence of environmental gases, we report here the first ab initio study of molecular adsorption on the double hexagonal close-packed (dhcp) americium (Am) (0 0 0 1) surface. Specifically, molecular oxygen adsorption on the (0 0 0 1) surface of dhcp Am has been studied in detail within the framework of density functional theory using a full-potential all-electron linearized augmented plane wave plus local orbitals (FP-LAPW+lo) method. Dissociative adsorption is found to be energetically more favorable compared to molecular adsorption. Chemisorption energies were optimized with respect to the distance of adsorbates from the surface for three approach positions at three adsorption sites, namely t1 (one-fold top), b2 (two-fold bridge), and h3 (three-fold hollow) sites. Chemisorption energies were computed at the scalar-relativistic-no-spin-orbit-coupling (SR-NSOC) and at the fully relativistic-with-spin-orbit-coupling (FR-SOC) levels of theory. The most stable configuration corresponds to a horizontal approach molecular dissociation with the oxygen atoms occupying neighboring h3 sites, with chemisorption energies at the NSOC and SOC theoretical levels being 9.395 and 9.886 eV, respectively. The corresponding distances of the oxygen molecule from the surface and oxygen-oxygen distance were found to be 0.953 and 3.731 Å, respectively. Overall our calculations indicate that chemisorption energies in cases with SOC are slightly more stable than those with NSOC in the 0.089-0.493 eV range. The work functions and net magnetic moments, respectively, increased and decreased in all cases compared to corresponding quantities of the bare dhcp-Am (0 0 0 1) surface. Adsorbate-substrate interactions have been analyzed in detail using partial charges inside muffin-tin spheres, difference charge density distributions, and
Ab Initio Studies of Stratospheric Ozone Depletion Chemistry
NASA Technical Reports Server (NTRS)
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.
THERMODYNAMICS OF MATERIALS: FROM AB INITIO TO PHENOMENOLOGY
Turchi, P A
2004-09-24
Quantum mechanical-based (or ab initio) methods are used to predict the stability properties of materials although their application is limited to relatively simple systems in terms of structures and number of alloy components. However thermodynamics of complex multi-component alloys requires a more versatile approach afforded within the CALPHAD formalism. Despite its success, the lack of experimental data very often prevents the design of robust thermodynamic databases. After a brief survey of ab initio methodologies and CALPHAD, it will be shown how ab initio electronic structure methods can supplement in two ways CALPHAD for subsequent applications. The first one is rather immediate and concerns the direct input of ab initio energetics in CALPHAD databases. The other way, more involved, is the assessment of ab initio thermodynamics '{acute a} la CALPHAD'. It will be shown how these results can be used within CALPHAD to predict the equilibrium properties of multi-component alloys. Finally, comments will be made on challenges and future prospects.
NASA Astrophysics Data System (ADS)
Chakrabarty, Aurab; Bouhali, Othmane; Mousseau, Normand; Becquart, Charlotte S.; El-Mellouhi, Fedwa
2016-08-01
Adsorption and dissociation of hydrocarbons on metallic surfaces represent crucial steps on the path to carburization, eventually leading to dusting corrosion. While adsorption of CO molecules on Fe surface is a barrier-less exothermic process, this is not the case for the dissociation of CO into C and O adatoms and the diffusion of C beneath the surface that are found to be associated with large energy barriers. In practice, these barriers can be affected by numerous factors that combine to favour the CO-Fe reaction such as the abundance of CO and other hydrocarbons as well as the presence of structural defects. From a numerical point of view, studying these factors is challenging and a step-by-step approach is necessary to assess, in particular, the influence of the finite box size on the reaction parameters for adsorption and dissociation of CO on metal surfaces. Here, we use density functional theory (DFT) total energy calculations with the climbing-image nudged elastic band method to estimate the adsorption energies and dissociation barriers for different CO coverages with surface supercells of different sizes. We further compute the effect of periodic boundary condition for DFT calculations and find that the contribution from van der Waals interaction in the computation of adsorption parameters is important as they contribute to correcting the finite-size error in small systems. The dissociation process involves carbon insertion into the Fe surface causing a lattice deformation that requires a larger surface system for unrestricted relaxation. We show that, in the larger surface systems associated with dilute CO-coverages, C-insertion is energetically more favourable, leading to a significant decrease in the dissociation barrier. This observation suggests that a large surface system with dilute coverage is necessary for all similar metal-hydrocarbon reactions in order to study their fundamental electronic mechanisms, as an isolated phenomenon, free from
Barrett, B R; Navratil, P; Vary, J P
2011-04-11
A long-standing goal of nuclear theory is to determine the properties of atomic nuclei based on the fundamental interactions among the protons and neutrons (i.e., nucleons). By adopting nucleon-nucleon (NN), three-nucleon (NNN) and higher-nucleon interactions determined from either meson-exchange theory or QCD, with couplings fixed by few-body systems, we preserve the predictive power of nuclear theory. This foundation enables tests of nature's fundamental symmetries and offers new vistas for the full range of complex nuclear phenomena. Basic questions that drive our quest for a microscopic predictive theory of nuclear phenomena include: (1) What controls nuclear saturation; (2) How the nuclear shell model emerges from the underlying theory; (3) What are the properties of nuclei with extreme neutron/proton ratios; (4) Can we predict useful cross sections that cannot be measured; (5) Can nuclei provide precision tests of the fundamental laws of nature; and (6) Under what conditions do we need QCD to describe nuclear structure, among others. Along with other ab initio nuclear theory groups, we have pursued these questions with meson-theoretical NN interactions, such as CD-Bonn and Argonne V18, that were tuned to provide high-quality descriptions of the NN scattering phase shifts and deuteron properties. We then add meson-theoretic NNN interactions such as the Tucson-Melbourne or Urbana IX interactions. More recently, we have adopted realistic NN and NNN interactions with ties to QCD. Chiral perturbation theory within effective field theory ({chi}EFT) provides us with a promising bridge between QCD and hadronic systems. In this approach one works consistently with systems of increasing nucleon number and makes use of the explicit and spontaneous breaking of chiral symmetry to expand the strong interaction in terms of a dimensionless constant, the ratio of a generic small momentum divided by the chiral symmetry breaking scale taken to be about 1 GeV/c. The resulting NN
Atomic and electronic structure of hydrogen on ZnO (1bar 100) surface: ab initio hybrid calculations
NASA Astrophysics Data System (ADS)
Usseinov, A. B.; Kotomin, E. A.; Zhukovskii, Yu F.; Purans, J.; Sorokin, A. V.; Akilbekov, A. T.
2013-12-01
Hydrogen atoms unavoidably incorporated into ZnO during growth of bulk samples and thin films considerably affect their electrical conductivity. The results of first principles hybrid LCAO calculations are discussed for hydrogen atoms in the bulk and on the non-polar ZnO (1bar 100) surface. The incorporation energy, the atomic relaxation, the electronic density redistribution and the electronic structure modifications are compared for the surface adsorption and bulk interstitial H positions. It is shown that hydrogen has a strong binding with the surface O ions (2.7 eV) whereas its incorporation into bulk is energetically unfavorable. Surface hydrogen atoms are very shallow donors, thus, contributing to the electronic conductivity.
NASA Astrophysics Data System (ADS)
Banerjee, Paramita; Das, G. P.
2016-05-01
In our effort to do first principles design of suitable materials for hydrogen storage, we have explored the interaction characteristics of a hydrogen molecule with pure as well as a 3d-transition metal (TM) atom doped Mg(0001) surface using density functional theory (DFT) based approach. Doping of a 3d-TM atom by creating a vacancy on the top most layer of Mg(0001) surface, enhances the molecular hydrogen adsorption efficiency of this surface by ~ 6 times. The TM atom gains some charge from the defected site of the Mg(0001) surface, becomes anionic and adsorbs the hydrogen molecule via Anti Kubas-type interaction. The interaction energy of this H2 molecule, including van der Waals dispersion correction, turns out to be ~ 0.4 eV, which falls in the right energy window between physisorption and chemisorption. On full coverage of this 3d-TM atom doped Mg(0001) surface with hydrogen molecules, the gravimetric density of hydrogen has been estimated to be ~ 5.6 wt %, thereby satisfying the criteria set by the department of energy (DOE) for efficient hydrogen storage.
Ab initio study of the p-hole magnetism at polar surfaces of ZnO: the role of correlations.
Fischer, Guntram; Sanchez, Nadiezhda; Adeagbo, Waheed; Szotek, Zdzislawa; Temmerman, Walter M; Ernst, Arthur; Hoffmann, Martin; Hergert, Wolfram; Muñoz, M Carmen
2016-01-13
A standard local density approximation and its self-interaction corrected version are applied to study spontaneous magnetization, promoted by localized p electron holes, of polar oxygen-terminated ZnO surfaces. The electronic properties and magnetic exchange interactions of three different facets are calculated. It is demonstrated that partially filled oxygen p orbitals of the polar surfaces exhibit magnetic moment formation and long range magnetic order leading to the occurrence of a ferromagnetic ground state. Monte Carlo simulations predict Curie temperatures above room temperature. In contrast to isolated defects in bulk materials, applying correlation corrections to the localized p-like surface states does not lead to a collapse of magnetic interaction: as the weakening of the magnetic interaction, caused by the reduced electronic overlap, is compensated by a strengthening due to an increase of the magnetic moments, the ferromagnetism can principally persist above room temperature, provided a large hole concentration exists. PMID:26657257
Taylor, Christopher D; Lillard, R Scott
2009-01-01
Density functional theory was applied to the initial steps of uranium hydriding: surface phenomena, absorption, bulk transport and trapping. H adsorbs exothermically to the (0 0 1) surface, yet H absorption into the bulk is endothermic, with off-center octahedral absorption having the lowest absorption energy of 0.39 eV, relative to molecular H{sub 2}. H absorption in interstitial sites causes a local softening of the bulk modulus. Diffusion of H in unstrained {alpha}-U has a barrier of 0.6 eV. The energy of H absorption adjacent to the chemical impurities C, S, Si was lowered by an amount proportional to the size of the impurity atom, and the resulting lattice strain Si > S > C. Thus, impurities may promote hydriding by providing surfaces or prestrained zones for H uptake.
NASA Astrophysics Data System (ADS)
Varandas, A. J. C.
2013-04-01
The Combined-Hyperbolic-Inverse-Power-Representation method, which treats evenly both short- and long-range interactions, is used to fit an extensive set of ab initio points for HO2 previously utilized [Xu et al., J. Chem. Phys. 122, 244305 (2005), 10.1063/1.1944290] to develop a spline interpolant. The novel form is shown to perform accurately when compared with others, while quasiclassical trajectory calculations of the O + OH reaction clearly pinpoint the role of long-range forces at low temperatures.
NASA Astrophysics Data System (ADS)
Boilleau, Corentin; Suaud, Nicolas; Guihéry, Nathalie
2012-12-01
In spin-crossover (SCO) compounds exhibiting a light induced excited spin state trapping (LIESST) effect, the thermodynamic T1/2 and kinetic T(LIESST) temperature values depend on the features of the potential energy surfaces (PES) of the two lowest singlet and quintet states but also on vibrational contributions, collective effects, such as electrostatics, for instance, spin-orbit couplings to a lesser extent, etc. In this work, the question of the link between the shape of the PES of SCO compounds exhibiting a LIESST effect and their first coordination sphere structure is addressed from wave function theory based ab initio calculations. Fe(II) complexes based on model ligands suited to reproduce the main characteristics of the PES of such compounds are distorted to emphasize selectively the role played by the metal-ligand distances and the ligand-metal-ligand angles. The studied angular deformations are those usually observed in many Fe(L)2(NCS)2 complexes. It is shown that the larger the deformation between the low spin and high spin equilibrium geometries, the higher the energy barrier from the high spin state and the weaker the energy difference between the bottom of the wells. These results corroborate observations made by experimentalists on a large number of complexes. While the PES features only constitutes one of the contributions to these temperatures, it is worth noticing that, relating T1/2 to the energy difference between the bottoms of the singlet and quintet wells and the T(LIESST) to the energy barrier from the quintet bottom well, the same slope of the empirical law T(LIESST) = -0.3T1/2+T0 is observed.
Teng, Wen-Shuang; Moskaleva, Lyudmila V; Chen, Hui-Lung; Lin, M C
2013-07-18
The reaction of NCN with H atoms has been investigated by ab initio MO and RRKM theory calculations. The mechanisms for formation of major products on the doublet and quartet potential energy surfaces have been predicted at the CCSD(T) level of theory with the complete basis set limit. In addition, the heat of formation for NCN predicted at this rigorous level and those from five isogyric reactions are in close agreement with the best value based on the isodesmic process, (3)CCO + N2 = (3)NCN + CO, 109.4 kcal/mol, which lies within the two existing experimental values. The rate constants for the three possible reaction channels, H + NCN → CH + N2 (k(P1)), HCN + (4)N (k(QP1)), and HNC + (4)N (k(QP2)), have been calculated in the temperature range 298-3000 K. The results show that k(P1) is significantly higher than k(QP1) and k(QP2) and that the total rate constant agrees well with available experimental values in the whole temperature range studied. The kinetics of the reverse CH + N2 reaction has also been revisited at the CCSD(T)/CBS level; the predicted total rate constants at 760 Torr Ar pressure can be represented by kr = 4.01 × 10(-15) T(0.90) exp(-17.42 kcal mol(-1)/RT) cm(3) molecule(-1) s(-1) at T = 800-4000 K. The result agrees closely with the most recent experimental data and the best theoretical result of Harding et al. (J. Phys. Chem. A 2008, 112, 522) as well as that of Moskaleva and Lin (Proc. Combust. Inst. 2000, 28, 2393) evaluated with a steady-state approximation after a coding error correction made in this study. PMID:23755876
NASA Astrophysics Data System (ADS)
Buchachenko, A. A.; Grinev, T. A.; Kłos, J.; Bieske, E. J.; Szczȩśniak, M. M.; Chałasiński, G.
2003-12-01
Three-dimensional potential energy and dipole moment surfaces of the Cl--H2 system are calculated ab initio by means of a coupled cluster method with single and double excitations and noniterative correction to triple excitations with augmented correlation consistent quadruple-zeta basis set supplemented with bond functions, and represented in analytical forms. Variational calculations of the energy levels up to the total angular momentum J=25 provide accurate estimations of the measured rotational spectroscopic constants of the ground van der Waals levels n=0 of the Cl-⋯H2/D2 complexes although they underestimate the red shifts of the mid-infrared spectra with v=0→v=1 vibrational excitation of the monomer. They also attest to the accuracy of effective radial interaction potentials extracted previously from experimental data using the rotational RKR procedure. Vibrational predissociation of the Cl-⋯H2/D2(v=1) complexes is shown to follow near-resonant vibrational-to-rotational energy transfer mechanism so that more than 97% of the product monomers are formed in the highest accessible rotational level. This mechanism explains the strong variation of the predissociation rate with isotopic content and nuclear spin form of the complex. Strong deviation of the observed relative abundances of ortho and para forms of the complexes from those of the monomers is qualitatively explained by the secondary ligand exchange reactions in the ionic beam, within the simple thermal equilibrium model. Positions and intensities of the hot v=0, n=1→v=1, n=1 and combination v=0, n=0→v=1, n=1 bands are predicted, and implications to the photoelectron spectroscopy of the complex are briefly discussed.
NASA Astrophysics Data System (ADS)
Bartolomei, Massimiliano; Carmona-Novillo, Estela; Hernández, Marta I.; Campos-Martínez, José; Hernandez-Lamoneda, Ramón
2008-06-01
A new potential energy surface (PES) for the quintet state of rigid O2(3Σg-)+O2(3Σg-) has been obtained using restricted coupled-cluster theory with singles, doubles, and perturbative triple excitations [RCCSD(T)]. A large number of relative orientations of the monomers (65) and intermolecular distances (17) have been considered. A spherical harmonic expansion of the interaction potential has been built from the ab initio data. It involves 29 terms, as a consequence of the large anisotropy of the interaction. The spherically averaged term agrees quite well with the one obtained from analysis of total integral cross sections. The absolute minimum of the PES corresponds to the crossed (D2d) structure (X shape) with an intermolecular distance of 6.224 bohrs and a well depth of 16.27 meV. Interestingly, the PES presents another (local) minimum close in energy (15.66 meV) at 6.50 bohrs and within a planar skewed geometry (S shape). We find that the origin of this second structure is due to the orientational dependence of the spin-exchange interactions which break the spin degeneracy and leads to three distinct intermolecular PESs with singlet, triplet, and quintet multiplicities. The lowest vibrational bound states of the O2-O2 dimer have been obtained and it is found that they reflect the above mentioned topological features of the PES: The first allowed bound state for the 16O isotope has an X structure but the next state is just 0.12 meV higher in energy and exhibits an S shape.
NASA Astrophysics Data System (ADS)
Veiga, R. G. A.; Miwa, R. H.; McLean, A. B.
2016-03-01
We report first-principles calculations of the energetic stability and electronic properties of metal-phthalocyanine (MPc) molecules (M = Cr, Mn, Fe, Co, Ni, Cu, and Zn) adsorbed on the δ -doped Si(111)-B (√{3 }×√{3 }) reconstructed surface. (i) It can be seen that CrPc, MnPc, FePc, and CoPc are chemically anchored to the topmost Si atom. (ii) Contrastingly, the binding of the NiPc, CuPc, and ZnPc molecules to the Si (111 ) -B (√{3 }×√{3 }) surface is exclusively ruled by van der Waals interactions, the main implication being that these molecules may diffuse and rearrange to form clusters and/or self-organized structures on this surface. The electronic structure calculations reveal that in point (i), owing to the formation of the metal-Si covalent bond, the net magnetic moment of the molecule is quenched by 1 μB , remaining unchanged in point (ii). In particular, the magnetic moment of CuPc (1 μB ) is preserved after adsorption. Finally, we verify that the formation of ZnPc, CuPc, and NiPc molecular (self-assembled) arrangements on the Si(111)-B (√{3 }×√{3 } ) surface is energetically favorable, in good agreement with recent experimental findings.
Photo-driven oxidation of water on α-Fe2O3 surfaces: an ab initio study.
Nguyen, Manh-Thuong; Seriani, Nicola; Piccinin, Simone; Gebauer, Ralph
2014-02-14
Adopting the theoretical scheme developed by the Nørskov group [see, for example, Nørskov et al., J. Phys. Chem. B 108, 17886 (2004)], we conducted a density functional theory study of photo-driven oxidation processes of water on various terminations of the clean hematite (α-Fe2O3) (0001) surface, explicitly taking into account the strong correlation among the 3d states of iron through the Hubbard U parameter. Six best-known terminations, namely, Fe−Fe−O3− (we call S1), O−Fe−Fe−(S2), O2−Fe−Fe−(S3), O3−Fe−Fe− (S4), Fe−O3−Fe− (S5), and O−Fe−O3−(S6), are first exposed to water, the stability of resulting surfaces is investigated under photoelectrochemical conditions by considering different chemical reactions (and their reaction free energies) that lead to surfaces covered by O atoms or/and OH groups. Assuming that the water splitting reaction is driven by the redox potential for photogenerated holes with respect to the normal hydrogen electrode, UVB, at voltage larger than UVB, most 3-oxygen terminated substrates are stable. These results thus suggest that the surface, hydroxylated in the dark, should release protons under illumination. Considering the surface free energy of all the possible terminations shows that O3–S5 and O3–S1 are the most thermodynamically stable. While water oxidation process on the former requires an overpotential of 1.22 V, only 0.84 V is needed on the latter. PMID:24527933
Sun, Geng; Jiang, Hong
2015-12-21
A comprehensive understanding of surface thermodynamics and kinetics based on first-principles approaches is crucial for rational design of novel heterogeneous catalysts, and requires combining accurate electronic structure theory and statistical mechanics modeling. In this work, ab initio molecular dynamics (AIMD) combined with the integrated tempering sampling (ITS) method has been explored to study thermodynamic and kinetic properties of elementary processes on surfaces, using a simple reaction CH2⇌CH+H on the Ni(111) surface as an example. By a careful comparison between the results from ITS-AIMD simulation and those evaluated in terms of the harmonic oscillator (HO) approximation, it is found that the reaction free energy and entropy from the HO approximation are qualitatively consistent with the results from ITS-AIMD simulation, but there are also quantitatively significant discrepancies. In particular, the HO model misses the entropy effects related to the existence of multiple adsorption configurations arising from the frustrated translation and rotation motion of adsorbed species, which are different in the reactant and product states. The rate constants are evaluated from two ITS-enhanced approaches, one using the transition state theory (TST) formulated in terms of the potential of mean force (PMF) and the other one combining ITS with the transition path sampling (TPS) technique, and are further compared to those based on harmonic TST. It is found that the rate constants from the PMF-based TST are significantly smaller than those from the harmonic TST, and that the results from PMF-TST and ITS-TPS are in a surprisingly good agreement. These findings indicate that the basic assumptions of transition state theory are valid in such elementary surface reactions, but the consideration of statistical averaging of all important adsorption configurations and reaction pathways, which are missing in the harmonic TST, are critical for accurate description of
ForceFit: a code to fit classical force fields to ab-initio potential energy surfaces
Henson, Neil Jon; Waldher, Benjamin; Kuta, Jadwiga; Clark, Aurora; Clark, Aurora E
2009-01-01
The ForceFit program package has been developed for fitting classical force field parameters based upon a force matching algorithm to quantum mechanical gradients of configurations that span the potential energy surface of the system. The program, which runs under Unix and is written in C++, is an easy to use, nonproprietary platform that enables gradient fitting of a wide variety of functional force field forms to quantum mechanical information obtained from an array of common electronic structure codes. All aspects of the fitting process are run from a graphical user interface, from the parsing of quantum mechanical data, assembling of a potential energy surface database, setting the force field and variables to be optimized, choosing a molecular mechanics code for comparison to the reference data, and finally, the initiation of a least squares minimization algorithm. Furthermore, the code is based on a modular templated code design that enables the facile addition of new functionality to the program.
Crowley, Jason M; Tahir-Kheli, Jamil; Goddard, William A
2015-10-01
It has been established experimentally that Bi2Te3 and Bi2Se3 are topological insulators, with zero band gap surface states exhibiting linear dispersion at the Fermi energy. Standard density functional theory (DFT) methods such as PBE lead to large errors in the band gaps for such strongly correlated systems, while more accurate GW methods are too expensive computationally to apply to the thin films studied experimentally. We show here that the hybrid B3PW91 density functional yields GW-quality results for these systems at a computational cost comparable to PBE. The efficiency of our approach stems from the use of Gaussian basis functions instead of plane waves or augmented plane waves. This remarkable success without empirical corrections of any kind opens the door to computational studies of real chemistry involving the topological surface state, and our approach is expected to be applicable to other semiconductors with strong spin-orbit coupling. PMID:26722872
Onida, Barbara; Borello, Luisa; Busco, Claudia; Ugliengo, Piero; Goto, Yasutomo; Inagaki, Shinji; Garrone, Edoardo
2005-06-23
Joint IR and computational results allow a detailed characterization of the surface properties of a mesoporous benzene-silica hybrid material with crystal-like wall structure. After outgassing at 450 degrees C, hydroxyl species mainly consist of noninteracting silanols, with both O-H and Si-O stretching modes at lower frequencies than those of SiOH in silica. Interaction with several probe molecules, followed both by experiment and calculus, shows that the aryl group in the coordination sphere of Si imparts a lesser acidity with respect to the isolated silanol in silica. In contrast, adsorption isotherms indicate that the interaction with acetone is stronger with benzene-silica than with silica: this is interpreted in terms of secondary interactions taking place between the slightly acidic CH in acetone and the electronic cloud in benzene-like rings. This suggests that both the inorganic component and the organic one play a role in dictating the surface behavior. PMID:16852474
Strak, Pawel; Sakowski, Konrad; Kempisty, Pawel
2015-09-07
Properties of bare and nitrogen-covered Al-terminated AlN(0001) surface were determined using density functional theory (DFT) calculations. At a low nitrogen coverage, the Fermi level is pinned by Al broken bond states located below conduction band minimum. Adsorption of nitrogen is dissociative with an energy gain of 6.05 eV/molecule at a H3 site creating an overlap with states of three neighboring Al surface atoms. During this adsorption, electrons are transferred from Al broken bond to topmost N adatom states. Accompanying charge transfer depends on the Fermi level. In accordance with electron counting rule (ECR), the DFT results confirm the Fermi level is not pinned at the critical value of nitrogen coverage θ{sub N}(1) = 1/4 monolayer (ML), but it is shifted from an Al-broken bond state to Np{sub z} state. The equilibrium thermodynamic potential of nitrogen in vapor depends drastically on the Fermi level pinning being shifted by about 4 eV for an ECR state at 1/4 ML coverage. For coverage above 1/4 ML, adsorption is molecular with an energy gain of 1.5 eV at a skewed on-top position above an Al surface atom. Electronic states of the admolecule are occupied as in the free molecule, no electron transfer occurs and adsorption of a N{sub 2} molecule does not depend on the Fermi level. The equilibrium pressure of molecular nitrogen above an AlN(0001) surface depends critically on the Fermi level position, being very low and very high for low and high coverage, respectively. From this fact, one can conclude that at typical growth conditions, the Fermi level is not pinned, and the adsorption and incorporation of impurities depend on the position of Fermi level in the bulk.
NASA Astrophysics Data System (ADS)
Buchachenko, A. A.; Kroupnov, A. A.; Kovalev, V. L.
2015-08-01
Elementary stage rate coefficients of the full system of kinetic equations describing heterogeneous catalytic recombination of the dissociated air on the surfaces of thermal protective ceramic coatings of β-cristobalite and α-Al2O3 are determined using the quantum-mechanical calculations within the framework of cluster models and literature data. Both the impact and associative recombination processes of adsorbed oxygen and nitrogen atoms are taken into account.
Kowalewski, Markus Mukamel, Shaul
2015-07-28
Femtosecond Stimulated Raman Spectroscopy (FSRS) signals that monitor the excited state conical intersections dynamics of acrolein are simulated. An effective time dependent Hamiltonian for two C—H vibrational marker bands is constructed on the fly using a local mode expansion combined with a semi-classical surface hopping simulation protocol. The signals are obtained by a direct forward and backward propagation of the vibrational wave function on a numerical grid. Earlier work is extended to fully incorporate the anharmonicities and intermode couplings.
NASA Astrophysics Data System (ADS)
Kowalewski, Markus; Mukamel, Shaul
2015-07-01
Femtosecond Stimulated Raman Spectroscopy (FSRS) signals that monitor the excited state conical intersections dynamics of acrolein are simulated. An effective time dependent Hamiltonian for two C—H vibrational marker bands is constructed on the fly using a local mode expansion combined with a semi-classical surface hopping simulation protocol. The signals are obtained by a direct forward and backward propagation of the vibrational wave function on a numerical grid. Earlier work is extended to fully incorporate the anharmonicities and intermode couplings.
Topological insulator Bi2Se3(111) surface doped with transition metals: An ab initio investigation
NASA Astrophysics Data System (ADS)
Abdalla, L. B.; Seixas, L.; Schmidt, T. M.; Miwa, R. H.; Fazzio, A.
2013-07-01
Based on first-principles calculations, we provide a map of the energetic stability and the electronic behavior of the topologically protected surface states of the topological insulator (TI), Bi2Se3, doped with transition metals (TMs), TM/Bi2Se3(111). We find that Fe, Mn, and Cr impurities are energetically more stable at the Bi substitutional sites, whereas we may find energetically stable substitutional as well as interstitial configurations for Co and Ni impurity atoms. Through scanning tunneling microscopy simulations, we verify that each TM atomic species and its position in the Bi2Se3(111) surface can be identified. The substitutional Fe and Cr impurities exhibit an energetic preference for the out-of-plane net magnetic moment, giving rise to a small energy gap at the Dirac point (DP), whereas the in-plane magnetic moment of substitutional Mn/Bi2Se3(111) promotes a shift of the DP from the center of the surface Brillouin zone, opening up a small energy gap. For the substitutional impurities, the shapes of the metallic bands are somewhat preserved compared with the energy bands of the pristine Bi2Se3(111) surface. Interstitial Co atoms also present an in-plane net magnetic moment, where we find the formation of metallic bands, suppressing the presence of the DP. For the Ni/Bi2Se3(111) impurity there is not a net magnetic moment, and therefore, the DP is preserved. Further formation energy results indicate other plausible (meta)stable impurity configurations, giving rise to quite different scenarios for the topological properties of TM/Bi2Se3(111), even for the same TM impurity.
NASA Astrophysics Data System (ADS)
Strak, Pawel; Sakowski, Konrad; Kempisty, Pawel; Krukowski, Stanislaw
2015-09-01
Properties of bare and nitrogen-covered Al-terminated AlN(0001) surface were determined using density functional theory (DFT) calculations. At a low nitrogen coverage, the Fermi level is pinned by Al broken bond states located below conduction band minimum. Adsorption of nitrogen is dissociative with an energy gain of 6.05 eV/molecule at a H3 site creating an overlap with states of three neighboring Al surface atoms. During this adsorption, electrons are transferred from Al broken bond to topmost N adatom states. Accompanying charge transfer depends on the Fermi level. In accordance with electron counting rule (ECR), the DFT results confirm the Fermi level is not pinned at the critical value of nitrogen coverage θN(1) = 1/4 monolayer (ML), but it is shifted from an Al-broken bond state to Npz state. The equilibrium thermodynamic potential of nitrogen in vapor depends drastically on the Fermi level pinning being shifted by about 4 eV for an ECR state at 1/4 ML coverage. For coverage above 1/4 ML, adsorption is molecular with an energy gain of 1.5 eV at a skewed on-top position above an Al surface atom. Electronic states of the admolecule are occupied as in the free molecule, no electron transfer occurs and adsorption of a N2 molecule does not depend on the Fermi level. The equilibrium pressure of molecular nitrogen above an AlN(0001) surface depends critically on the Fermi level position, being very low and very high for low and high coverage, respectively. From this fact, one can conclude that at typical growth conditions, the Fermi level is not pinned, and the adsorption and incorporation of impurities depend on the position of Fermi level in the bulk.