An Accurate Potential Energy Surface for H2O
NASA Technical Reports Server (NTRS)
Schwenke, David W.; Partridge, Harry; Langhoff, Stephen R. (Technical Monitor)
1997-01-01
We have carried out extensive high quality ab initio electronic structure calculations of the ground state potential energy surface (PES) and dipole moment function (DMF) for H2O. A small adjustment is made to the PES to improve the agreement of line positions from theory and experiment. The theoretical line positions are obtained from variational ro-vibrational calculations using the exact kinetic energy operator. For the lines being fitted, the root-mean-square error was reduced from 6.9 to 0.08 /cm. We were then able to match 30,092 of the 30,117 lines from the HITRAN 96 data base to theoretical lines, and 80% of the line positions differed less than 0.1 /cm. About 3% of the line positions in the experimental data base appear to be incorrect. Theory predicts the existence of many additional weak lines with intensities above the cutoff used in the data base. To obtain results of similar accuracy for HDO, a mass dependent correction to the PH is introduced and is parameterized by simultaneously fitting line positions for HDO and D2O. The mass dependent PH has good predictive value for T2O and HTO. Nonadiabatic effects are not explicitly included. Line strengths for vibrational bands summed over rotational levels usually agree well between theory and experiment, but individual line strengths can differ greatly. A high temperature line list containing about 380 million lines has been generated using the present PES and DMF
NASA Astrophysics Data System (ADS)
Rhee, Young Min
2000-10-01
A modified method to construct an accurate potential energy surface by interpolation is presented. The modification is based on the use of Cartesian coordinates in the weighting function. The translational and rotational invariance of the potential is incorporated by a proper definition of the distance between two Cartesian configurations. A numerical algorithm to find the distance is developed. It is shown that the present method is more exact in describing a planar system compared to the previous methods with weightings in internal coordinates. The applicability of the method to reactive systems is also demonstrated by performing classical trajectory simulations on the surface.
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.
Highly accurate potential energy surface for the He-H2 dimer.
Bakr, Brandon W; Smith, Daniel G A; Patkowski, Konrad
2013-10-14
A new highly accurate interaction potential is constructed for the He-H2 van der Waals complex. This potential is fitted to 1900 ab initio energies computed at the very large-basis coupled-cluster level and augmented by corrections for higher-order excitations (up to full configuration interaction level) and the diagonal Born-Oppenheimer correction. At the vibrationally averaged H-H bond length of 1.448736 bohrs, the well depth of our potential, 15.870 ± 0.065 K, is nearly 1 K larger than the most accurate previous studies have indicated. In addition to constructing our own three-dimensional potential in the van der Waals region, we present a reparameterization of the Boothroyd-Martin-Peterson potential surface [A. I. Boothroyd, P. G. Martin, and M. R. Peterson, J. Chem. Phys. 119, 3187 (2003)] that is suitable for all configurations of the triatomic system. Finally, we use the newly developed potentials to compute the properties of the lone bound states of (4)He-H2 and (3)He-H2 and the interaction second virial coefficient of the hydrogen-helium mixture. PMID:24116617
Tung, Wei-Cheng; Adamowicz, Ludwik
2014-03-28
Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH(+) ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations. PMID:24697449
NASA Astrophysics Data System (ADS)
Tung, Wei-Cheng; Adamowicz, Ludwik
2014-03-01
Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH+ ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations.
Fedorov, Dmitry A.; Varganov, Sergey A.; Derevianko, Andrei
2014-05-14
We calculate the potential energy curves, the permanent dipole moment curves, and the lifetimes of the ground and excited vibrational states of the heteronuclear alkali dimers XY (X, Y = Li, Na, K, Rb, Cs) in the X{sup 1}Σ{sup +} electronic state using the coupled cluster with singles doubles and triples method. All-electron quadruple-ζ basis sets with additional core functions are used for Li and Na, and small-core relativistic effective core potentials with quadruple-ζ quality basis sets are used for K, Rb, and Cs. The inclusion of the coupled cluster non-perturbative triple excitations is shown to be crucial for obtaining the accurate potential energy curves. A large one-electron basis set with additional core functions is needed for the accurate prediction of permanent dipole moments. The dissociation energies are overestimated by only 14 cm{sup −1} for LiNa and by no more than 114 cm{sup −1} for the other molecules. The discrepancies between the experimental and calculated harmonic vibrational frequencies are less than 1.7 cm{sup −1}, and the discrepancies for the anharmonic correction are less than 0.1 cm{sup −1}. We show that correlation between atomic electronegativity differences and permanent dipole moment of heteronuclear alkali dimers is not perfect. To obtain the vibrational energies and wave functions the vibrational Schrödinger equation is solved with the B-spline basis set method. The transition dipole moments between all vibrational states, the Einstein coefficients, and the lifetimes of the vibrational states are calculated. We analyze the decay rates of the vibrational states in terms of spontaneous emission, and stimulated emission and absorption induced by black body radiation. In all studied heteronuclear alkali dimers the ground vibrational states have much longer lifetimes than any excited states.
NASA Astrophysics Data System (ADS)
Fedorov, Dmitry A.; Derevianko, Andrei; Varganov, Sergey A.
2014-05-01
We calculate the potential energy curves, the permanent dipole moment curves, and the lifetimes of the ground and excited vibrational states of the heteronuclear alkali dimers XY (X, Y = Li, Na, K, Rb, Cs) in the X1Σ+ electronic state using the coupled cluster with singles doubles and triples method. All-electron quadruple-ζ basis sets with additional core functions are used for Li and Na, and small-core relativistic effective core potentials with quadruple-ζ quality basis sets are used for K, Rb, and Cs. The inclusion of the coupled cluster non-perturbative triple excitations is shown to be crucial for obtaining the accurate potential energy curves. A large one-electron basis set with additional core functions is needed for the accurate prediction of permanent dipole moments. The dissociation energies are overestimated by only 14 cm-1 for LiNa and by no more than 114 cm-1 for the other molecules. The discrepancies between the experimental and calculated harmonic vibrational frequencies are less than 1.7 cm-1, and the discrepancies for the anharmonic correction are less than 0.1 cm-1. We show that correlation between atomic electronegativity differences and permanent dipole moment of heteronuclear alkali dimers is not perfect. To obtain the vibrational energies and wave functions the vibrational Schrödinger equation is solved with the B-spline basis set method. The transition dipole moments between all vibrational states, the Einstein coefficients, and the lifetimes of the vibrational states are calculated. We analyze the decay rates of the vibrational states in terms of spontaneous emission, and stimulated emission and absorption induced by black body radiation. In all studied heteronuclear alkali dimers the ground vibrational states have much longer lifetimes than any excited states.
Temperature dependent effective potential method for accurate free energy calculations of solids
NASA Astrophysics Data System (ADS)
Hellman, Olle; Steneteg, Peter; Abrikosov, I. A.; Simak, S. I.
2013-03-01
We have developed a thorough and accurate method of determining anharmonic free energies, the temperature dependent effective potential technique (TDEP). It is based on ab initio molecular dynamics followed by a mapping onto a model Hamiltonian that describes the lattice dynamics. The formalism and the numerical aspects of the technique are described in detail. A number of practical examples are given, and results are presented, which confirm the usefulness of TDEP within ab initio and classical molecular dynamics frameworks. In particular, we examine from first principles the behavior of force constants upon the dynamical stabilization of the body centered phase of Zr, and show that they become more localized. We also calculate the phase diagram for 4He modeled with the Aziz potential and obtain results which are in favorable agreement both with respect to experiment and established techniques.
Communication: An accurate global potential energy surface for the ground electronic state of ozone
Dawes, Richard E-mail: hguo@unm.edu; Lolur, Phalgun; Li, Anyang; Jiang, Bin; Guo, Hua E-mail: hguo@unm.edu
2013-11-28
We report a new full-dimensional and global potential energy surface (PES) for the O + O{sub 2} → O{sub 3} ozone forming reaction based on explicitly correlated multireference configuration interaction (MRCI-F12) data. It extends our previous [R. Dawes, P. Lolur, J. Ma, and H. Guo, J. Chem. Phys. 135, 081102 (2011)] dynamically weighted multistate MRCI calculations of the asymptotic region which showed the widely found submerged reef along the minimum energy path to be the spurious result of an avoided crossing with an excited state. A spin-orbit correction was added and the PES tends asymptotically to the recently developed long-range electrostatic model of Lepers et al. [J. Chem. Phys. 137, 234305 (2012)]. This PES features: (1) excellent equilibrium structural parameters, (2) good agreement with experimental vibrational levels, (3) accurate dissociation energy, and (4) most-notably, a transition region without a spurious reef. The new PES is expected to allow insight into the still unresolved issues surrounding the kinetics, dynamics, and isotope signature of ozone.
Very accurate potential energy curve of the LiH molecule
NASA Astrophysics Data System (ADS)
Tung, Wei-Cheng; Pavanello, Michele; Adamowicz, Ludwik
2011-02-01
We present very accurate calculations of the ground-state potential energy curve (PEC) of the LiH molecule performed with all-electron explicitly correlated Gaussian functions with shifted centers. The PEC is generated with the variational method involving simultaneous optimization of all Gaussians with an approach employing the analytical first derivatives of the energy with respect to the Gaussian nonlinear parameters (i.e., the exponents and the coordinates of the shifts). The LiH internuclear distance is varied between 1.8 and 40 bohrs. The absolute accuracy of the generated PEC is estimated as not exceeding 0.3 cm-1. The adiabatic corrections for the four LiH isotopologues, i.e., 7LiH, 6LiH, 7LiD, and 6LiD, are also calculated and added to the LiH PEC. The aforementioned PECs are then used to calculate the vibrational energies for these systems. The maximum difference between the computed and the experimental vibrational transitions is smaller than 0.9 cm-1. The contribution of the adiabatic correction to the dissociation energy of 7LiH molecule is 10.7 cm-1. The magnitude of this correction shows its importance in calculating the LiH spectroscopic constants. As the estimated contribution of the nonadiabatic and relativistic effects to the ground state dissociation energy is around 0.3 cm-1, their inclusion in the LiH PEC calculation seems to be the next most important contribution to evaluate in order to improve the accuracy achieved in this work.
Very accurate potential energy curve of the LiH molecule.
Tung, Wei-Cheng; Pavanello, Michele; Adamowicz, Ludwik
2011-02-14
We present very accurate calculations of the ground-state potential energy curve (PEC) of the LiH molecule performed with all-electron explicitly correlated Gaussian functions with shifted centers. The PEC is generated with the variational method involving simultaneous optimization of all Gaussians with an approach employing the analytical first derivatives of the energy with respect to the Gaussian nonlinear parameters (i.e., the exponents and the coordinates of the shifts). The LiH internuclear distance is varied between 1.8 and 40 bohrs. The absolute accuracy of the generated PEC is estimated as not exceeding 0.3 cm(-1). The adiabatic corrections for the four LiH isotopologues, i.e., (7)LiH, (6)LiH, (7)LiD, and (6)LiD, are also calculated and added to the LiH PEC. The aforementioned PECs are then used to calculate the vibrational energies for these systems. The maximum difference between the computed and the experimental vibrational transitions is smaller than 0.9 cm(-1). The contribution of the adiabatic correction to the dissociation energy of (7)LiH molecule is 10.7 cm(-1). The magnitude of this correction shows its importance in calculating the LiH spectroscopic constants. As the estimated contribution of the nonadiabatic and relativistic effects to the ground state dissociation energy is around 0.3 cm(-1), their inclusion in the LiH PEC calculation seems to be the next most important contribution to evaluate in order to improve the accuracy achieved in this work. PMID:21322671
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.
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).
NASA Astrophysics Data System (ADS)
Lu-Lu, Zhang; Yu-Zhi, Song; Shou-Bao, Gao; Yuan, Zhang; Qing-Tian, Meng
2016-05-01
A globally accurate single-sheeted double many-body expansion potential energy surface is reported for the first excited state of HS2 by fitting the accurate ab initio energies, which are calculated at the multireference configuration interaction level with the aug-cc-pVQZ basis set. By using the double many-body expansion-scaled external correlation method, such calculated ab initio energies are then slightly corrected by scaling their dynamical correlation. A grid of 2767 ab initio energies is used in the least-square fitting procedure with the total root-mean square deviation being 1.406 kcal·mol‑1. The topographical features of the HS2(A2A‧) global potential energy surface are examined in detail. The attributes of the stationary points are presented and compared with the corresponding ab initio results as well as experimental and other theoretical data, showing good agreement. The resulting potential energy surface of HS2(A2A‧) can be used as a building block for constructing the global potential energy surfaces of larger S/H molecular systems and recommended for dynamic studies on the title molecular system. Project supported by the National Natural Science Foundation of China (Grant No. 11304185), the Taishan Scholar Project of Shandong Province, China, the Shandong Provincial Natural Science Foundation, China (Grant No. ZR2014AM022), the Shandong Province Higher Educational Science and Technology Program, China (Grant No. J15LJ03), the China Postdoctoral Science Foundation (Grant No. 2014M561957), and the Post-doctoral Innovation Project of Shandong Province, China (Grant No. 201402013).
Tung, Wei-Cheng; Adamowicz, Ludwik
2014-03-28
Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH{sup +} ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations.
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.}
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Dateo, Christopher E.; Schwenke, David W.; Chaban, Galina M.
2005-01-01
Accurate quartic force fields have been determined for the CCH- and NH2- molecular anions using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T). Very large one-particle basis sets have been used including diffuse functions and up through g-type functions. Correlation of the nitrogen and carbon core electrons has been included, as well as other "small" effects, such as the diagonal Born-Oppenheimer correction, and basis set extrapolation, and corrections for higher-order correlation effects and scalar relativistic effects. Fundamental vibrational frequencies have been computed using standard second-order perturbation theory as well as variational methods. Comparison with the available experimental data is presented and discussed. The implications of our research for the astronomical observation of molecular anions will be discussed.
Accurate double many-body expansion potential energy surface for the 2(1)A' state of N2O.
Li, Jing; Varandas, António J C
2014-08-28
An accurate double many-body expansion potential energy surface is reported for the 2(1)A' state of N2O. The new double many-body expansion (DMBE) form has been fitted to a wealth of ab initio points that have been calculated at the multi-reference configuration interaction level using the full-valence-complete-active-space wave function as reference and the cc-pVQZ basis set, and subsequently corrected semiempirically via double many-body expansion-scaled external correlation method to extrapolate the calculated energies to the limit of a complete basis set and, most importantly, the limit of an infinite configuration interaction expansion. The topographical features of the novel potential energy surface are then examined in detail and compared with corresponding attributes of other potential functions available in the literature. Exploratory trajectories have also been run on this DMBE form with the quasiclassical trajectory method, with the thermal rate constant so determined at room temperature significantly enhancing agreement with experimental data. PMID:25173014
Accurate double many-body expansion potential energy surface for the 21A' state of N_2O
NASA Astrophysics Data System (ADS)
Li, Jing; Varandas, António J. C.
2014-08-01
An accurate double many-body expansion potential energy surface is reported for the 21A' state of N_2O. The new double many-body expansion (DMBE) form has been fitted to a wealth of ab initio points that have been calculated at the multi-reference configuration interaction level using the full-valence-complete-active-space wave function as reference and the cc-pVQZ basis set, and subsequently corrected semiempirically via double many-body expansion-scaled external correlation method to extrapolate the calculated energies to the limit of a complete basis set and, most importantly, the limit of an infinite configuration interaction expansion. The topographical features of the novel potential energy surface are then examined in detail and compared with corresponding attributes of other potential functions available in the literature. Exploratory trajectories have also been run on this DMBE form with the quasiclassical trajectory method, with the thermal rate constant so determined at room temperature significantly enhancing agreement with experimental data.
Zou, Lindong; Li, Jun; Wang, Hui; Ma, Jianyi; Guo, Hua
2015-07-16
Full-dimensional quantum dynamics studies of the photodetachment of HCO2(-) and DCO2(-) are reported using a wave-packet method on an accurate global potential energy surface of the neutral HOCO/HCO2 system. The calculated photoelectron spectra reproduced both the positions and widths of the main HCO2 and DCO2 peaks observed in experiment. Specifically, both the (2)A1 and (2)B2 resonance peaks of the neutral radicals were identified in our simulations thanks to the adiabatic PES that captures both the (2)A1 and (2)B2 minima. The narrow widths and isotope effect of the lowest resonances are indicative of tunneling-facilitated predissociation. Furthermore, the dissociation product CO2 was found to be excited in both its symmetric stretching and bending modes, which are coupled via a strong Fermi resonance, but rotationally cold, in good agreement with the recent photoelectron-photodetachment coincidence experiments. PMID:25607218
NASA Astrophysics Data System (ADS)
Delahaye, Thibault; Nikitin, Andrei; Rey, Michaël; Szalay, Péter G.; Tyuterev, Vladimir G.
2014-09-01
In this paper we report a new ground state potential energy surface for ethylene (ethene) C2H4 obtained from extended ab initio calculations. The coupled-cluster approach with the perturbative inclusion of the connected triple excitations CCSD(T) and correlation consistent polarized valence basis set cc-pVQZ was employed for computations of electronic ground state energies. The fit of the surface included 82 542 nuclear configurations using sixth order expansion in curvilinear symmetry-adapted coordinates involving 2236 parameters. A good convergence for variationally computed vibrational levels of the C2H4 molecule was obtained with a RMS(Obs.-Calc.) deviation of 2.7 cm-1 for fundamental bands centers and 5.9 cm-1 for vibrational bands up to 7800 cm-1. Large scale vibrational and rotational calculations for 12C2H4, 13C2H4, and 12C2D4 isotopologues were performed using this new surface. Energy levels for J = 20 up to 6000 cm-1 are in a good agreement with observations. This represents a considerable improvement with respect to available global predictions of vibrational levels of 13C2H4 and 12C2D4 and rovibrational levels of 12C2H4.
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.
Dybeck, Eric C; Schieber, Natalie P; Shirts, Michael R
2016-08-01
We examine the free energies of three benzene polymorphs as a function of temperature in the point-charge OPLS-AA and GROMOS54A7 potentials as well as the polarizable AMOEBA09 potential. For this system, using a polarizable Hamiltonian instead of the cheaper point-charge potentials is shown to have a significantly smaller effect on the stability at 250 K than on the lattice energy at 0 K. The benzene I polymorph is found to be the most stable crystal structure in all three potentials examined and at all temperatures examined. For each potential, we report the free energies over a range of temperatures and discuss the added value of using full free energy methods over the minimized lattice energy to determine the relative crystal stability at finite temperatures. The free energies in the polarizable Hamiltonian are efficiently calculated using samples collected in a cheaper point-charge potential. The polarizable free energies are estimated from the point-charge trajectories using Boltzmann reweighting with MBAR. The high configuration-space overlap necessary for efficient Boltzmann reweighting is achieved by designing point-charge potentials with intramolecular parameters matching those in the expensive polarizable Hamiltonian. Finally, we compare the computational cost of this indirect reweighted free energy estimate to the cost of simulating directly in the expensive polarizable Hamiltonian. PMID:27341280
NASA Astrophysics Data System (ADS)
Shi, De-Heng; Liu, Yu-Fang; Sun, Jin-Feng; Zhu, Zun-Lue; Yang, Xiang-Dong
2006-12-01
The reasonable dissociation limit of the second excited singlet state B1Π of 7LiH molecule is obtained. The accurate dissociation energy and equilibrium geometry of the B1Π state are calculated using a symmetry-adapted-cluster configuration-interaction method in full active space. The whole potential energy curve for the B1Π state is obtained over the internuclear distance ranging from about 0.10 nm to 0.54 nm, and has a least-square fit to the analytic Murrell-Sorbie function form. The vertical excitation energy is calculated from the ground state to the B1Π state and compared with previous theoretical results. The equilibrium internuclear distance obtained by geometry optimization is found to be quite different from that obtained by single-point energy scanning under the same calculation condition. Based on the analytic potential energy function, the harmonic frequency value of the B1Π state is estimated. A comparison of the theoretical calculations of dissociation energies, equilibrium interatomic distances and the analytic potential energy function with those obtained by previous theoretical results clearly shows that the present work is more comprehensive and in better agreement with experiments than previous theories, thus it is an improvement on previous theories.
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)
Shi, De-Heng; Liu, Yu-Fang; Sun, Jin-Feng; Yang, Xiang-Dong; Zhu, Zun-Lue
2006-05-01
The reasonable dissociation limit of the A1Σ+ state 7LiH molecule is obtained. The accurate dissociation energy and the equilibrium geometry of this state are calculated using a symmetry-adapted-cluster configuration-interaction method in complete active space for the first time. The whole potential energy curve and the dipole moment function for the A1Σ+ state are calculated over a wide internuclear separation range from about 0.1 to 1.4 nm. The calculated equilibrium geometry and dissociation energy of this potential energy curve are of Re=0.2487 nm and De=1.064 eV, respectively. The unusual negative values of the anharmonicity constant and the vibration-rotational coupling constant are of ωeχe=-4.7158cm-1 and αe=-0.08649cm-1, respectively. The vertical excitation energy from the ground to the A1Σ+ state is calculated and the value is of 3.613 eV at 0.15875 nm (the equilibrium position of the ground state). The highly anomalous shape of this potential energy curve, which is exceptionally flat over a wide radial range around the equilibrium position, is discussed in detail. The harmonic frequency value of 502.47cm-1 about this state is approximately estimated. Careful comparison of the theoretical determinations with those obtained by previous theories about the A1Σ+ state dissociation energy clearly shows that the present calculations are much closer to the experiments than previous theories, thus represents an improvement.
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.}
A highly accurate interatomic potential for argon
NASA Astrophysics Data System (ADS)
Aziz, Ronald A.
1993-09-01
A modified potential based on the individually damped model of Douketis, Scoles, Marchetti, Zen, and Thakkar [J. Chem. Phys. 76, 3057 (1982)] is presented which fits, within experimental error, the accurate ultraviolet (UV) vibration-rotation spectrum of argon determined by UV laser absorption spectroscopy by Herman, LaRocque, and Stoicheff [J. Chem. Phys. 89, 4535 (1988)]. Other literature potentials fail to do so. The potential also is shown to predict a large number of other properties and is probably the most accurate characterization of the argon interaction constructed to date.
Mizukami, Wataru Tew, David P.; Habershon, Scott
2014-10-14
We present a new approach to semi-global potential energy surface fitting that uses the least absolute shrinkage and selection operator (LASSO) constrained least squares procedure to exploit an extremely flexible form for the potential function, while at the same time controlling the risk of overfitting and avoiding the introduction of unphysical features such as divergences or high-frequency oscillations. Drawing from a massively redundant set of overlapping distributed multi-dimensional Gaussian functions of inter-atomic separations we build a compact full-dimensional surface for malonaldehyde, fit to explicitly correlated coupled cluster CCSD(T)(F12*) energies with a root mean square deviations accuracy of 0.3%–0.5% up to 25 000 cm{sup −1} above equilibrium. Importance-sampled diffusion Monte Carlo calculations predict zero point energies for malonaldehyde and its deuterated isotopologue of 14 715.4(2) and 13 997.9(2) cm{sup −1} and hydrogen transfer tunnelling splittings of 21.0(4) and 3.2(4) cm{sup −1}, respectively, which are in excellent agreement with the experimental values of 21.583 and 2.915(4) cm{sup −1}.
NASA Astrophysics Data System (ADS)
Ishikawa, Atsushi; Nakai, Hiromi
2016-04-01
Gibbs free energy of hydration of a proton and standard hydrogen electrode potential were evaluated using high-level quantum chemical calculations. The solvent effect was included using the cluster-continuum model, which treated short-range effects by quantum chemical calculations of proton-water complexes, and the long-range effects by a conductor-like polarizable continuum model. The harmonic solvation model (HSM) was employed to estimate enthalpy and entropy contributions due to nuclear motions of the clusters by including the cavity-cluster interactions. Compared to the commonly used ideal gas model, HSM treatment significantly improved the contribution of entropy, showing a systematic convergence toward the experimental data.
Huang, Xinchuan; Schwenke, David W; Lee, Timothy J
2014-03-21
A purely ab initio potential energy surface (PES) was refined with selected (32)S(16)O2 HITRAN data. Compared to HITRAN, the root-mean-squares error (σ(RMS)) for all J = 0-80 rovibrational energy levels computed on the refined PES (denoted Ames-1) is 0.013 cm(-1). Combined with a CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), an infrared (IR) line list (denoted Ames-296K) has been computed at 296 K and covers up to 8000 cm(-1). Compared to the HITRAN and CDMS databases, the intensity agreement for most vibrational bands is better than 85%-90%. Our predictions for (34)S(16)O2 band origins, higher energy (32)S(16)O2 band origins and missing (32)S(16)O2 IR bands have been verified by most recent experiments and available HITRAN data. We conclude that the Ames-1 PES is able to predict (32/34)S(16)O2 band origins below 5500 cm(-1) with 0.01-0.03 cm(-1) uncertainties, and the Ames-296K line list provides continuous, reliable and accurate IR simulations. The K(a)-dependence of both line position and line intensity errors is discussed. The line list will greatly facilitate SO2 IR spectral experimental analysis, as well as elimination of SO2 lines in high-resolution astronomical observations. PMID:24655184
NASA Technical Reports Server (NTRS)
Huang, Xinchuan; Schwenke, David W.; Lee, Timothy J.
2014-01-01
A purely ab initio potential energy surface (PES) was refined with selected (32)S(16)O2 HITRAN data. Compared to HITRAN, the root-mean-squares error (RMS) error for all J=0-80 rovibrational energy levels computed on the refined PES (denoted Ames-1) is 0.013 cm(exp -1). Combined with a CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), an infrared (IR) line list (denoted Ames-296K) has been computed at 296K and covers up to 8,000 cm(exp -1). Compared to the HITRAN and CDMS databases, the intensity agreement for most vibrational bands is better than 85-90%. Our predictions for (34)S(16)O2 band origins, higher energy (32)S(16)O2 band origins and missing (32)S(16)O2 IR bands have been verified by most recent experiments and available HITRAN data. We conclude that the Ames-1 PES is able to predict (32/34)S(16)O2 band origins below 5500 cm(exp -1) with 0.01-0.03 cm(exp -1) uncertainties, and the Ames-296K line list provides continuous, reliable and accurate IR simulations. The Ka-dependence of both line position and line intensity errors is discussed. The line list will greatly facilitate SO2 IR spectral experimental analysis, as well as elimination of SO2 lines in high-resolution astronomical observations.
Liu, Tianhui; Zhang, Zhaojun; Fu, Bina; Yang, Xueming; Zhang, Dong H
2016-03-16
The mode-specific dynamics for the dissociative chemisorption of H2O on Cu(111) is first investigated by seven-dimensional quantum dynamics calculations, based on an accurately fitted potential energy surface (PES) recently developed by neural network fitting to DFT energy points. It is indicated that excitations in all three vibrational modes have a significant impact on reactivity, which are more efficacious than increasing the translational energy in promoting the reaction, with the largest enhancement for the excitation in the asymmetric stretching mode. There is large discrepancy between the six-dimensional reactivities with fixed azimuthal angles and seven-dimensional results, revealing that the 6D "flat surface" model cannot accurately characterize the reaction dynamics. The azimuthal angle-averaging approach is validated for vibrational excited states of the reactant, where the 7D mode-specific probability can be well reproduced by averaging the 6D azimuthal angle-fixed probabilities over 18 angles. PMID:26941197
Accurate potential energy functions, non-adiabatic and spin-orbit couplings in the ZnH+ system
NASA Astrophysics Data System (ADS)
Liang, Guiying; Liu, Xiaoting; Zhang, Xiaomei; Xu, Haifeng; Yan, Bing
2016-03-01
A high-level ab initio calculation on the ZnH+ cation has been carried out with the multi-reference configuration interaction method plus Davison correction (MRCI + Q). The scalar relativistic effect is included by using the Douglas-Kroll-Hess (DKH) method. The calculated potential energy curves (PECs) of the 7 Λ-S states are associated with the dissociation limits of Zn+(2Sg) + H(2Sg), Zn(1Sg) + H+(1Sg), and Zn+(2Pu) + H(2Sg), respectively (The Λ-S state is labeled as 2S + 1Λ, in which Λ is the quantum number for the projection along the internuclear axis of the total electronic orbital angular momentum and S is the total electron spin). The spectroscopic constants of the bound states are determined and in good agreement with the available theoretical and experimental results. The permanent dipole moments (PDMs) of Λ-S states and the spin-orbit (SO) matrix elements between Λ-S states are also computed. The results show that the abrupt changes of the PDMs and SO matrix elements come into being for the reason of the avoided crossing between the states with the same symmetry. In addition, the non-adiabatic couplings matrix elements between Λ-S states are also evaluated. Finally, the spin-orbit couplings (SOCs) for the low-lying states are considered with Breit-Pauli operator. The SOC effect makes the 7 Λ-S states of the ZnH+ cation split into 12 Ω states (Ω = Λ + Sz, in which Sz is projection of the total electron spin S along the internuclear Z-axis). For the (3)0+ state, the two energy minima exhibit in the potential, which could be attributed to the formation of the new avoided crossing point. The transition dipole moments (TDMs), Franck-Condon factors, and the radiative lifetimes of the selected transitions (2)0+-X0+, (3)0+-X0+, (2)1-X0+ and (3)1-X0+ have been reported.
Accurate potential energy functions, non-adiabatic and spin-orbit couplings in the ZnH(+) system.
Liang, Guiying; Liu, Xiaoting; Zhang, Xiaomei; Xu, Haifeng; Yan, Bing
2016-03-01
A high-level ab initio calculation on the ZnH(+) cation has been carried out with the multi-reference configuration interaction method plus Davison correction (MRCI+Q). The scalar relativistic effect is included by using the Douglas-Kroll-Hess (DKH) method. The calculated potential energy curves (PECs) of the 7 Λ-S states are associated with the dissociation limits of Zn(+)((2)Sg)+H((2)Sg), Zn((1)Sg)+H(+)((1)Sg), and Zn(+)((2)Pu)+H((2)Sg), respectively (The Λ-S state is labeled as (2S+1)Λ, in which Λ is the quantum number for the projection along the internuclear axis of the total electronic orbital angular momentum and S is the total electron spin). The spectroscopic constants of the bound states are determined and in good agreement with the available theoretical and experimental results. The permanent dipole moments (PDMs) of Λ-S states and the spin-orbit (SO) matrix elements between Λ-S states are also computed. The results show that the abrupt changes of the PDMs and SO matrix elements come into being for the reason of the avoided crossing between the states with the same symmetry. In addition, the non-adiabatic couplings matrix elements between Λ-S states are also evaluated. Finally, the spin-orbit couplings (SOCs) for the low-lying states are considered with Breit-Pauli operator. The SOC effect makes the 7 Λ-S states of the ZnH(+) cation split into 12 Ω states (Ω=Λ+Sz, in which Sz is projection of the total electron spin S along the internuclear Z-axis). For the (3)0(+) state, the two energy minima exhibit in the potential, which could be attributed to the formation of the new avoided crossing point. The transition dipole moments (TDMs), Franck-Condon factors, and the radiative lifetimes of the selected transitions (2)0(+)-X0(+), (3)0(+)-X0(+), (2)1-X0(+) and (3)1-X0(+) have been reported. PMID:26637984
NASA Astrophysics Data System (ADS)
Li, Yong-Qing; Song, Yu-Zhi; Joaquim de Campos Varandas, António
2015-01-01
An accurate single-sheeted double many-body expansion potential energy surface is reported for the title system. It is obtained by using the aug-cc-pVTZ and aug-cc-pVQZ basis sets with extrapolation of the electron correlation energy to the complete basis set limit, plus extrapolation to the complete basis set limit of the complete-active-space self-consistent field energy. The collinear and bending barrier heights of the new global potential energy surface is 2.301 and 1.768 kcal mol-1, in very good agreement with the values of 2.222 and 1.770 kcal mol-1 from the current best potential energy surface. In particular, the new potential energy surface describes well the important van der Waals interactions which is very useful for investigating the dynamics of the title system. Thus, the new potential energy surface can both be recommended for dynamics studies of the F + H2 reaction and as building block for constructing the potential energy surfaces of larger fluorine/hydrogen containing systems. Based on the new potential energy surface, a preliminary theoretical study of the reaction F(2P) + H2 (X1 Σg+) → FH(X1Σ+) + H(2S) has been carried out with the methods of quasi-classical trajectory and quantum mechanical. The results have shown that the new PES is suitable for any kind of dynamics studies. Supplementary material in the form of one pdf file available from the Journal web page at http://dx.doi.org/10.1140/epjd/e2014-50445-3
NASA Astrophysics Data System (ADS)
Nabok, Dmitrii; Gulans, Andris; Draxl, Claudia
2016-07-01
The G W approach of many-body perturbation theory has become a common tool for calculating the electronic structure of materials. However, with increasing number of published results, discrepancies between the values obtained by different methods and codes become more and more apparent. For a test set of small- and wide-gap semiconductors, we demonstrate how to reach the numerically best electronic structure within the framework of the full-potential linearized augmented plane-wave (FLAPW) method. We first evaluate the impact of local orbitals in the Kohn-Sham eigenvalue spectrum of the underlying starting point. The role of the basis-set quality is then further analyzed when calculating the G0W0 quasiparticle energies. Our results, computed with the exciting code, are compared to those obtained using the projector-augmented plane-wave formalism, finding overall good agreement between both methods. We also provide data produced with a typical FLAPW basis set as a benchmark for other G0W0 implementations.
NASA Astrophysics Data System (ADS)
Meng, Qingyong; Chen, Jun; Zhang, Dong H.
2015-09-01
The ring polymer molecular dynamics (RPMD) calculations are performed to calculate rate constants for the title reaction on the recently constructed potential energy surface based on permutation invariant polynomial (PIP) neural-network (NN) fitting [J. Li et al., J. Chem. Phys. 142, 204302 (2015)]. By inspecting convergence, 16 beads are used in computing free-energy barriers at 300 K ≤ T ≤ 1000 K, while different numbers of beads are used for transmission coefficients. The present RPMD rates are in excellent agreement with quantum rates computed on the same potential energy surface, as well as with the experimental measurements, demonstrating further that the RPMD is capable of producing accurate rates for polyatomic chemical reactions even at rather low temperatures.
Song, Y Z; Varandas, A J C
2009-04-01
A single-sheeted potential energy surface is reported for the electronic ground-state of H(2)S by fitting accurate multireference configuration interaction energies calculated using aug-cc-pVTZ and aug-cc-pVQZ basis sets with extrapolation of the electron correlation energy to the complete basis set limit, plus extrapolation to the complete basis set limit of the complete-active-space self-consistent field energy. A switching function formalism has been used to warrant the correct behavior at the H(2)(X (1)Sigma(g) (+))+S((1)D) and SH(X (2)Pi)+H((2)S) dissociation limits. The topographical features of the novel global potential energy surface are examined in detail, with the former being used for exploratory quasiclassical trajectory calculations of the thermal rate constant for the S((1)D)+H(2), S((1)D)+D(2), and S((1)D)+HD reactions at room temperature. A comparison with other available potential energy surfaces as well as kinetics data is also provided. PMID:19355742
Li, Yongqing; Yuan, Jiuchuang; Chen, Maodu; Ma, Fengcai; Sun, Mengtao
2013-07-15
An accurate single-sheeted double many-body expansion potential energy surface is reported for the title system. A switching function formalism has been used to warrant the correct behavior at the H2(X1Σg+)+N(2D) and NH (X3Σ-)+H(2S) dissociation channels involving nitrogen in the ground N(4S) and first excited N(2D) states. The topographical features of the novel 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. The novel surface can be using to treat well the Renner-Teller degeneracy of the 12A″ and 12A' states of NH 2. Such a work can both be recommended for dynamics studies of the N(2D)+H2 reaction and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen-containing systems. In turn, a test theoretical study of the reaction N(2D)+H2(X1Σg+)(ν=0,j=0)→NH (X3Σ-)+H(2S) has been carried out with the method of quantum wave packet on the new potential energy surface. Reaction probabilities, integral cross sections, and differential cross sections have been calculated. Threshold exists because of the energy barrier (68.5 meV) along the minimum energy path. On the curve of reaction probability for total angular momentum J = 0, there are two sharp peaks just above threshold. The value of integral cross section increases quickly from zero to maximum with the increase of collision energy, and then stays stable with small oscillations. The differential cross section result shows that the reaction is a typical forward and backward scatter in agreement with experimental measurement result. PMID:23666848
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)
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.
Bytautas, Laimutis; Matsunaga, Nikita; Ruedenberg, Klaus
2010-02-21
In the first paper of this series, a very accurate ab initio potential energy curve of the (3)Sigma(g)(-) ground state of O(2) has been determined in the approximation that all valence shell electron correlations were calculated at the complete basis set limit. In the present study, the corrections arising from core electron correlations and relativity effects, viz., spin-orbit coupling and scalar relativity, are determined and added to the potential energy curve. From the 24 points calculated on this curve, an analytical expression in terms of even-tempered Gaussian functions is determined and, from it, the vibrational and rotational energy levels are calculated by means of the discrete variable representation. We find 42 vibrational levels. Experimental data (from the Schumann-Runge band system) only yield the lowest 36 levels due to significant reduction in the transition intensities of higher levels. For the 35 term values G(v), the mean absolute deviation between theoretical and experimental data is 12.8 cm(-1). The dissociation energy with respect to the lowest vibrational energy is calculated within 25 cm(-1) of the experimental value of 41,268.2+/-3 cm(-1). The theoretical crossing between the (3)Sigma(g)(-) state and the (1)Sigma(g)(+) state is found to occur at 2.22 A and the spin-orbit coupling in this region is analyzed. PMID:20170227
Li, Jun; Guo, Hua
2015-12-14
A globally accurate full-dimensional potential energy surface (PES) for the OH + CH4 → H2O + CH3 reaction is developed using the permutation invariant polynomial-neural network approach based on ∼135,000 points at the level of correlated coupled cluster singles, doubles, and perturbative triples level with the augmented correlation consistent polarized valence triple-zeta basis set. The total root mean square fitting error is only 3.9 meV or 0.09 kcal/mol. This PES is shown to reproduce energies, geometries, and harmonic frequencies of stationary points along the reaction path. Kinetic and dynamical calculations on the PES indicated a good agreement with the available experimental data. PMID:26671351
NASA Astrophysics Data System (ADS)
Li, Jun; Guo, Hua
2015-12-01
A globally accurate full-dimensional potential energy surface (PES) for the OH + CH4 → H2O + CH3 reaction is developed using the permutation invariant polynomial-neural network approach based on ˜135 000 points at the level of correlated coupled cluster singles, doubles, and perturbative triples level with the augmented correlation consistent polarized valence triple-zeta basis set. The total root mean square fitting error is only 3.9 meV or 0.09 kcal/mol. This PES is shown to reproduce energies, geometries, and harmonic frequencies of stationary points along the reaction path. Kinetic and dynamical calculations on the PES indicated a good agreement with the available experimental data.
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
Jiang, Bin; Hu, Xixi; Lin, Sen; Xie, Daiqian; Guo, Hua
2015-09-28
Cobalt is a widely used catalyst for many heterogeneous reactions, including the Fischer-Tropsch (FT) process, which converts syngas (H2 and CO) to higher hydrocarbons. As a result, a better understanding of the key chemical steps on the Co surface, such as the dissociative chemisorption of H2 as an initial step of the FT process, is of fundamental importance. Here, we report an accurate full-dimensional global potential energy surface for the dissociative chemisorption of H2 on the rigid Co(0001) surface constructed from more than 3000 density functional theory points. The high-fidelity potential energy surface was obtained using the permutation invariant polynomial-neural network method, which preserves both the permutation symmetry of H2 and translational symmetry of the Co(0001) surface. The reaction path features a very low barrier on the top site. Full-dimensional quantum dynamical calculations provide insights into the dissociation dynamics and influence of the initial vibrational, rotational, and orientational degrees of freedom. PMID:26286861
Janke, Svenja M; Auerbach, Daniel J; Wodtke, Alec M; Kandratsenka, Alexander
2015-09-28
We have constructed a potential energy surface (PES) for H-atoms interacting with fcc Au(111) based on fitting the analytic form of the energy from Effective Medium Theory (EMT) to ab initio energy values calculated with density functional theory. The fit used input from configurations of the H-Au system with Au atoms at their lattice positions as well as configurations with the Au atoms displaced from their lattice positions. It reproduces the energy, in full dimension, not only for the configurations used as input but also for a large number of additional configurations derived from ab initio molecular dynamics (AIMD) trajectories at finite temperature. Adiabatic molecular dynamics simulations on this PES reproduce the energy loss behavior of AIMD. EMT also provides expressions for the embedding electron density, which enabled us to develop a self-consistent approach to simulate nonadiabatic electron-hole pair excitation and their effect on the motion of the incident H-atoms. For H atoms with an energy of 2.7 eV colliding with Au, electron-hole pair excitation is by far the most important energy loss pathway, giving an average energy loss ≈3 times that of the adiabatic case. This increased energy loss enhances the probability of the H-atom remaining on or in the Au slab by a factor of 2. The most likely outcome for H-atoms that are not scattered also depends prodigiously on the energy transfer mechanism; for the nonadiabatic case, more than 50% of the H-atoms which do not scatter are adsorbed on the surface, while for the adiabatic case more than 50% pass entirely through the 4 layer simulation slab. PMID:26429033
Janke, Svenja M.; Auerbach, Daniel J.; Kandratsenka, Alexander; Wodtke, Alec M.
2015-09-28
We have constructed a potential energy surface (PES) for H-atoms interacting with fcc Au(111) based on fitting the analytic form of the energy from Effective Medium Theory (EMT) to ab initio energy values calculated with density functional theory. The fit used input from configurations of the H–Au system with Au atoms at their lattice positions as well as configurations with the Au atoms displaced from their lattice positions. It reproduces the energy, in full dimension, not only for the configurations used as input but also for a large number of additional configurations derived from ab initio molecular dynamics (AIMD) trajectories at finite temperature. Adiabatic molecular dynamics simulations on this PES reproduce the energy loss behavior of AIMD. EMT also provides expressions for the embedding electron density, which enabled us to develop a self-consistent approach to simulate nonadiabatic electron-hole pair excitation and their effect on the motion of the incident H-atoms. For H atoms with an energy of 2.7 eV colliding with Au, electron-hole pair excitation is by far the most important energy loss pathway, giving an average energy loss ≈3 times that of the adiabatic case. This increased energy loss enhances the probability of the H-atom remaining on or in the Au slab by a factor of 2. The most likely outcome for H-atoms that are not scattered also depends prodigiously on the energy transfer mechanism; for the nonadiabatic case, more than 50% of the H-atoms which do not scatter are adsorbed on the surface, while for the adiabatic case more than 50% pass entirely through the 4 layer simulation slab.
Palma, Juliana; Manthe, Uwe
2015-12-17
Quasiclassical trajectories (QCT) have been employed to elucidate the effect of exciting the C-H bond in F + CHD3 collisions. The calculations were performed on a new potential energy surface that accurately describes the van der Waals complexes in the entrance channel of the reaction. It was found that exciting the C-H bond significantly enhances the yield of HF + CD3, whereas it has a minor effect on the production of DF + CHD2. Therefore, the net effect is that the total reactivity increases upon excitation. This result strongly contradicts recent experimental findings. Significant differences in regard to the yield of each product channel were also found between QCT results calculated with the new surface and those obtained with the surface previously developed by Czakó et al. This shows that relatively small variations in the topography of the entrance channel can result in huge discrepancies in the predicted DF/HF branching ratio. However, in regard to other attributes of the reaction, the agreement between QCT results computed with different surfaces, and between them and experimental results, is good. For the F + CHD3 → HF + CD3 reaction, at a collisional energy of 9.0 kcal/mol, experiments and QCT calculations agree, indicating that the extra energy deposited in the C-H bond is channelled into the HF product. In addition, the angular distribution of CD3 is backward oriented and is not sensitive to the excitation of the C-H bond. PMID:26270126
NASA Astrophysics Data System (ADS)
Huang, Xinchuan; Schwenke, David W.; Lee, Timothy J.
2009-06-01
A further refined, global potential energy surface (PES) is computed for the C_{3v} symmetry isotopologues of ammonia, including ^{14}NH_3, ^{15}NH_3, ^{14}ND_3 and ^{14}NT_3. The refinement procedure was similar to that used in our previously reported PES, but now extends to higher J energy levels and other isotopologues. Both the diagonal Born-Oppenheimer correction and the non-adiabatic correction were included. J=0-6 rovibrational energy levels and transition frequencies of ^{14}NH_3 computed on this PES are in excellent agreement with HITRAN data. Statistics on nearly 4100 transitions and more than 1000 energy levels demonstrate the accuracy achieved by the state-of-the-art "Best Theory + Experiment" strategy. Most transition frequencies are of ±0.01-0.02 cm^{-1} accuracy. Similar accuracy has been found on ^{15}NH_3 J=0-3 rovibrational energy levels. Several transitions and energy levels in HITRAN have been identified as unreliable or suspicious, and some have been re-assigned. For ^{14}ND_3 and ^{14}NT_3, J=0-3 calculations have been performed. Agreement for pure rotation-inversion transitions is nearly perfect, with more reliable energy levels presented. On the other hand, our J=0 results suggest a re-analysis on the ^{14}ND_3 ν_1 band origin is needed. Finally, we will discuss possible future refinements leading to an even better final PES for Ammonia. X. Huang, D.W. Schwenke, and T.J. Lee, J. Chem. Phys. 129, 214304 (2008).
Homayoon, Zahra
2014-09-28
A new, full (nine)-dimensional potential energy surface and dipole moment surface to describe the NO{sup +}(H{sub 2}O) cluster is reported. The PES is based on fitting of roughly 32 000 CCSD(T)-F12/aug-cc-pVTZ electronic energies. The surface is a linear least-squares fit using a permutationally invariant basis with Morse-type variables. The PES is used in a Diffusion Monte Carlo study of the zero-point energy and wavefunction of the NO{sup +}(H{sub 2}O) and NO{sup +}(D{sub 2}O) complexes. Using the calculated ZPE the dissociation energies of the clusters are reported. Vibrational configuration interaction calculations of NO{sup +}(H{sub 2}O) and NO{sup +}(D{sub 2}O) using the MULTIMODE program are performed. The fundamental, a number of overtone, and combination states of the clusters are reported. The IR spectrum of the NO{sup +}(H{sub 2}O) cluster is calculated using 4, 5, 7, and 8 modes VSCF/CI calculations. The anharmonic, coupled vibrational calculations, and IR spectrum show very good agreement with experiment. Mode coupling of the water “antisymmetric” stretching mode with the low-frequency intermolecular modes results in intensity borrowing.
NASA Astrophysics Data System (ADS)
Szabó, István; Telekes, Hajnalka; Czakó, Gábor
2015-06-01
We develop a full-dimensional global analytical potential energy surface (PES) for the F- + CH3F reaction by fitting about 50 000 energy points obtained by an explicitly correlated composite method based on the second-order Møller-Plesset perturbation-F12 and coupled-cluster singles, doubles, and perturbative triples-F12a methods and the cc-pVnZ-F12 [n = D, T] basis sets. The PES accurately describes the (a) back-side attack Walden inversion mechanism involving the pre- and post-reaction (b) ion-dipole and (c) hydrogen-bonded complexes, the configuration-retaining (d) front-side attack and (e) double-inversion substitution pathways, as well as (f) the proton-abstraction channel. The benchmark quality relative energies of all the important stationary points are computed using the focal-point analysis (FPA) approach considering electron correlation up to coupled-cluster singles, doubles, triples, and perturbative quadruples method, extrapolation to the complete basis set limit, core-valence correlation, and scalar relativistic effects. The FPA classical(adiabatic) barrier heights of (a), (d), and (e) are -0.45(-0.61), 46.07(45.16), and 29.18(26.07) kcal mol-1, respectively, the dissociation energies of (b) and (c) are 13.81(13.56) and 13.73(13.52) kcal mol-1, respectively, and the endothermicity of (f) is 42.54(38.11) kcal mol-1. Quasiclassical trajectory computations of cross sections, scattering (θ) and initial attack (α) angle distributions, as well as translational and internal energy distributions are performed for the F- + CH3F(v = 0) reaction using the new PES. Apart from low collision energies (Ecoll), the SN2 excitation function is nearly constant, the abstraction cross sections rapidly increase with Ecoll from a threshold of ˜40 kcal mol-1, and retention trajectories via double inversion are found above Ecoll = ˜ 30 kcal mol-1, and at Ecoll = ˜ 50 kcal mol-1, the front-side attack cross sections start to increase very rapidly. At low Ecoll, the
Szabó, István; Telekes, Hajnalka; Czakó, Gábor
2015-06-28
We develop a full-dimensional global analytical potential energy surface (PES) for the F(-) + CH3F reaction by fitting about 50 000 energy points obtained by an explicitly correlated composite method based on the second-order Møller-Plesset perturbation-F12 and coupled-cluster singles, doubles, and perturbative triples-F12a methods and the cc-pVnZ-F12 [n = D, T] basis sets. The PES accurately describes the (a) back-side attack Walden inversion mechanism involving the pre- and post-reaction (b) ion-dipole and (c) hydrogen-bonded complexes, the configuration-retaining (d) front-side attack and (e) double-inversion substitution pathways, as well as (f) the proton-abstraction channel. The benchmark quality relative energies of all the important stationary points are computed using the focal-point analysis (FPA) approach considering electron correlation up to coupled-cluster singles, doubles, triples, and perturbative quadruples method, extrapolation to the complete basis set limit, core-valence correlation, and scalar relativistic effects. The FPA classical(adiabatic) barrier heights of (a), (d), and (e) are -0.45(-0.61), 46.07(45.16), and 29.18(26.07) kcal mol(-1), respectively, the dissociation energies of (b) and (c) are 13.81(13.56) and 13.73(13.52) kcal mol(-1), respectively, and the endothermicity of (f) is 42.54(38.11) kcal mol(-1). Quasiclassical trajectory computations of cross sections, scattering (θ) and initial attack (α) angle distributions, as well as translational and internal energy distributions are performed for the F(-) + CH3F(v = 0) reaction using the new PES. Apart from low collision energies (Ecoll), the SN2 excitation function is nearly constant, the abstraction cross sections rapidly increase with Ecoll from a threshold of ∼40 kcal mol(-1), and retention trajectories via double inversion are found above Ecoll = ∼ 30 kcal mol(-1), and at Ecoll = ∼ 50 kcal mol(-1), the front-side attack cross sections start to increase very rapidly. At
NASA Astrophysics Data System (ADS)
Meng, Qingyong; Chen, Jun; Zhang, Dong H.
2016-04-01
To fast and accurately compute rate coefficients of the H/D + CH4 → H2/HD + CH3 reactions, we propose a segmented strategy for fitting suitable potential energy surface (PES), on which ring-polymer molecular dynamics (RPMD) simulations are performed. On the basis of recently developed permutation invariant polynomial neural-network approach [J. Li et al., J. Chem. Phys. 142, 204302 (2015)], PESs in local configuration spaces are constructed. In this strategy, global PES is divided into three parts, including asymptotic, intermediate, and interaction parts, along the reaction coordinate. Since less fitting parameters are involved in the local PESs, the computational efficiency for operating the PES routine is largely enhanced by a factor of ˜20, comparing with that for global PES. On interaction part, the RPMD computational time for the transmission coefficient can be further efficiently reduced by cutting off the redundant part of the child trajectories. For H + CH4, good agreements among the present RPMD rates and those from previous simulations as well as experimental results are found. For D + CH4, on the other hand, qualitative agreement between present RPMD and experimental results is predicted.
Huang, Xinchuan; Schwenke, David W; Tashkun, Sergey A; Lee, Timothy J
2012-03-28
An isotopic-independent, highly accurate potential energy surface (PES) has been determined for CO(2) by refining a purely ab initio PES with selected, purely experimentally determined rovibrational energy levels. The purely ab initio PES is denoted Ames-0, while the refined PES is denoted Ames-1. Detailed tests are performed to demonstrate the spectroscopic accuracy of the Ames-1 PES. It is shown that Ames-1 yields σ(rms) (root-mean-squares error) = 0.0156 cm(-1) for 6873 J = 0-117 (12)C(16)O(2) experimental energy levels, even though less than 500 (12)C(16)O(2) energy levels were included in the refinement procedure. It is also demonstrated that, without any additional refinement, Ames-1 yields very good agreement for isotopologues. Specifically, for the (12)C(16)O(2) and (13)C(16)O(2) isotopologues, spectroscopic constants G(v) computed from Ames-1 are within ±0.01 and 0.02 cm(-1) of reliable experimentally derived values, while for the (16)O(12)C(18)O, (16)O(12)C(17)O, (16)O(13)C(18)O, (16)O(13)C(17)O, (12)C(18)O(2), (17)O(12)C(18)O, (12)C(17)O(2), (13)C(18)O(2), (13)C(17)O(2), (17)O(13)C(18)O, and (14)C(16)O(2) isotopologues, the differences are between ±0.10 and 0.15 cm(-1). To our knowledge, this is the first time a polyatomic PES has been refined using such high J values, and this has led to new challenges in the refinement procedure. An initial high quality, purely ab initio dipole moment surface (DMS) is constructed and used to generate a 296 K line list. For most bands, experimental IR intensities are well reproduced for (12)C(16)O(2) using Ames-1 and the DMS. For more than 80% of the bands, the experimental intensities are reproduced with σ(rms)(ΔI) < 20% or σ(rms)(ΔI∕δ(obs)) < 5. A few exceptions are analyzed and discussed. Directions for future improvements are discussed, though it is concluded that the current Ames-1 and the DMS should be useful in analyzing and assigning high-resolution laboratory or astronomical spectra. PMID:22462861
An accurate model potential for alkali neon systems.
Zanuttini, D; Jacquet, E; Giglio, E; Douady, J; Gervais, B
2009-12-01
We present a detailed investigation of the ground and lowest excited states of M-Ne dimers, for M=Li, Na, and K. We show that the potential energy curves of these Van der Waals dimers can be obtained accurately by considering the alkali neon systems as one-electron systems. Following previous authors, the model describes the evolution of the alkali valence electron in the combined potentials of the alkali and neon cores by means of core polarization pseudopotentials. The key parameter for an accurate model is the M(+)-Ne potential energy curve, which was obtained by means of ab initio CCSD(T) calculation using a large basis set. For each MNe dimer, a systematic comparison with ab initio computation of the potential energy curve for the X, A, and B states shows the remarkable accuracy of the model. The vibrational analysis and the comparison with existing experimental data strengthens this conclusion and allows for a precise assignment of the vibrational levels. PMID:19968334
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.
An accurate analytic representation of the water pair potential.
Cencek, Wojciech; Szalewicz, Krzysztof; Leforestier, Claude; van Harrevelt, Rob; van der Avoird, Ad
2008-08-28
The ab initio water dimer interaction energies obtained from coupled cluster calculations and used in the CC-pol water pair potential (Bukowski et al., Science, 2007, 315, 1249) have been refitted to a site-site form containing eight symmetry-independent sites in each monomer and denoted as CC-pol-8s. Initially, the site-site functions were assumed in a B-spline form, which allowed a precise optimization of the positions of the sites. Next, these functions were assumed in the standard exponential plus inverse powers form. The root mean square error of the CC-pol-8s fit with respect to the 2510 ab initio points is 0.10 kcal mol(-1), compared to 0.42 kcal mol(-1) of the CC-pol fit (0.010 kcal mol(-1) compared to 0.089 kcal mol(-1) for points with negative interaction energies). The energies of the stationary points in the CC-pol-8s potential are considerably more accurate than in the case of CC-pol. The water dimer vibration-rotation-tunneling spectrum predicted by the CC-pol-8s potential agrees substantially and systematically better with experiment than the already very accurate spectrum predicted by CC-pol, while specific features that could not be accurately predicted previously now agree very well with experiment. This shows that the uncertainties of the fit were the largest source of error in the previous predictions and that the present potential sets a new standard of accuracy in investigations of the water dimer. PMID:18688514
Towards an accurate dissociative potential for water
NASA Astrophysics Data System (ADS)
Akin-Ojo, Omololu
2014-03-01
Most models of water describe the molecule as rigid, i.e., with fixed bond angles and bond lengths, or as flexible in which the bond angles and bond lengths vary but the chemical bonds cannot be broken. In this work we present our progress in the development of a water model which allows for the breaking and formation of chemical bonds. The force field was obtained by fitting ab initio (not DFT) energies, forces, and molecular properties. The ability of the model to predict properties of water at ambient and extreme conditions will be presented. We will also report on the modeling of small clusters of water using the dissociative force field.
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
Chen, Jun; Sun, Zhigang E-mail: zhangdh@dicp.ac.cn; Zhang, Dong H. E-mail: zhangdh@dicp.ac.cn
2015-01-14
A three dimensional potential energy surface for the F + H{sub 2} → HF + H reaction has been computed by the spin unrestricted coupled cluster method with singles, doubles, triples, and perturbative quadruples [UCCSDT(2){sub Q}] using the augmented correlation-consistent polarised valence quadruple zeta basis set for the fluorine atom and the correlation-consistent polarised valence quadruple zeta basis set for the hydrogen atom. All the calculations are based on the restricted open-shell Hartree-Fock orbitals, together with the frozen core approximations, and the UCCSD(T)/complete basis set (CBS) correction term was included. The global potential energy surface was calculated by fitting the sampled ab initio points without any scaling factor for the correlation energy part using a neutral network function method. Extensive dynamics calculations have been carried out on the potential energy surface. The reaction rate constants, integral cross sections, product rotational states distribution, and forward and backward scattering as a function of collision energy of the F + HD → HF + D, F + HD → DF + H, and F + H{sub 2} reaction, were calculated by the time-independent quantum dynamics scattering theory using the new surface. The satisfactory agreement with the reported experimental observations previously demonstrates the accuracy of the new potential energy surface.
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.
Araujo, Rafael B.; Almeida, J. de S; Ferreira da Silva, A.; Ahuja, Rajeev
2015-09-28
The main goals of this paper are to investigate the accuracy of the Tran-Blaha modified Becke Johnson (TB-mBJ) potential to predict the electronic structure of lithium iron phosphate and the related redox reaction energy with the lithium deintercalation process. The computed electronic structures show that the TB-mBJ method is able to partially localize Fe-3d electrons in LiFePO{sub 4} and FePO{sub 4} which usually is a problem for the generalized gradient approximation (GGA) due to the self interaction error. The energy band gap is also improved by the TB-mBJ calculations in comparison with the GGA results. It turned out, however, that the redox reaction energy evaluated by the TB-mBJ technique is not in good agreement with the measured one. It is speculated that this disagreement in the computed redox energy and the experimental value is due to the lack of a formal expression to evaluate the exchange and correlation energy. Therefore, the TB-mBJ is an efficient method to improve the prediction of the electronic structures coming form the standard GGA functional in LiFePO{sub 4} and FePO{sub 4}. However, it does not appear to have the same efficiency for evaluating the redox reaction energies for the investigated system.
NASA Astrophysics Data System (ADS)
B. Araujo, Rafael; S. de Almeida, J.; Ferreira da Silva, A.; Ahuja, Rajeev
2015-09-01
The main goals of this paper are to investigate the accuracy of the Tran-Blaha modified Becke Johnson (TB-mBJ) potential to predict the electronic structure of lithium iron phosphate and the related redox reaction energy with the lithium deintercalation process. The computed electronic structures show that the TB-mBJ method is able to partially localize Fe-3d electrons in LiFePO4 and FePO4 which usually is a problem for the generalized gradient approximation (GGA) due to the self interaction error. The energy band gap is also improved by the TB-mBJ calculations in comparison with the GGA results. It turned out, however, that the redox reaction energy evaluated by the TB-mBJ technique is not in good agreement with the measured one. It is speculated that this disagreement in the computed redox energy and the experimental value is due to the lack of a formal expression to evaluate the exchange and correlation energy. Therefore, the TB-mBJ is an efficient method to improve the prediction of the electronic structures coming form the standard GGA functional in LiFePO4 and FePO4. However, it does not appear to have the same efficiency for evaluating the redox reaction energies for the investigated system.
NASA Astrophysics Data System (ADS)
Seo, Dong-Hwa; Urban, Alexander; Ceder, Gerbrand
2015-09-01
Transition-metal (TM) oxides play an increasingly important role in technology today, including applications such as catalysis, solar energy harvesting, and energy storage. In many of these applications, the details of their electronic structure near the Fermi level are critically important for their properties. We propose a first-principles-based computational methodology for the accurate prediction of oxygen charge transfer in TM oxides and lithium TM (Li-TM) oxides. To obtain accurate electronic structures, the Heyd-Scuseria-Ernzerhof (HSE06) hybrid functional is adopted, and the amount of exact Hartree-Fock exchange (mixing parameter) is adjusted to reproduce reference band gaps. We show that the HSE06 functional with optimal mixing parameter yields not only improved electronic densities of states, but also better energetics (Li-intercalation voltages) for LiCo O2 and LiNi O2 as compared to the generalized gradient approximation (GGA), Hubbard U corrected GGA (GGA +U ), and standard HSE06. We find that the optimal mixing parameters for TM oxides are system specific and correlate with the covalency (ionicity) of the TM species. The strong covalent (ionic) nature of TM-O bonding leads to lower (higher) optimal mixing parameters. We find that optimized HSE06 functionals predict stronger hybridization of the Co 3 d and O 2 p orbitals as compared to GGA, resulting in a greater contribution from oxygen states to charge compensation upon delithiation in LiCo O2 . We also find that the band gaps of Li-TM oxides increase linearly with the mixing parameter, enabling the straightforward determination of optimal mixing parameters based on GGA (α =0.0 ) and HSE06 (α =0.25 ) calculations. Our results also show that G0W0@GGA +U band gaps of TM oxides (M O ,M =Mn ,Co ,Ni ) and LiCo O2 agree well with experimental references, suggesting that G0W0 calculations can be used as a reference for the calibration of the mixing parameter in cases when no experimental band gap has been
Accurate ionization potential of semiconductors from efficient density functional calculations
NASA Astrophysics Data System (ADS)
Ye, Lin-Hui
2016-07-01
Despite its huge successes in total-energy-related applications, the Kohn-Sham scheme of density functional theory cannot get reliable single-particle excitation energies for solids. In particular, it has not been able to calculate the ionization potential (IP), one of the most important material parameters, for semiconductors. We illustrate that an approximate exact-exchange optimized effective potential (EXX-OEP), the Becke-Johnson exchange, can be used to largely solve this long-standing problem. For a group of 17 semiconductors, we have obtained the IPs to an accuracy similar to that of the much more sophisticated G W approximation (GWA), with the computational cost of only local-density approximation/generalized gradient approximation. The EXX-OEP, therefore, is likely as useful for solids as for finite systems. For solid surfaces, the asymptotic behavior of the vx c has effects similar to those of finite systems which, when neglected, typically cause the semiconductor IPs to be underestimated. This may partially explain why standard GWA systematically underestimates the IPs and why using the same GWA procedures has not been able to get an accurate IP and band gap at the same time.
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
Boschen, Jeffery S.; Theis, Daniel; Ruedenberg, Klaus; Windus, Theresa L.
2013-12-07
The diatomic carbon molecule has a complex electronic structure with a large number of low-lying electronic excited states. In this work, the potential energy curves (PECs) of the four lowest lying singlet states (X-1 Sigma(+)(g), A(1)Pi(u), B-1 Delta(g), and B'(1)Sigma(+)(g)) were obtained by high-level ab initio calculations. Valence electron correlation was accounted for by the correlation energy extrapolation by intrinsic scaling (CEEIS) method. Additional corrections to the PECs included core-valence correlation and relativistic effects. Spin-orbit corrections were found to be insignificant. The impact of using dynamically weighted reference wave functions in conjunction with CEEIS was examined and found to give indistinguishable results from the even weighted method. The PECs showed multiple curve crossings due to the B-1 Delta(g) state as well as an avoided crossing between the two (1)Sigma(+)(g) states. Vibrational energy levels were computed for each of the four electronic states, as well as rotational constants and spectroscopic parameters. Comparison between the theoretical and experimental results showed excellent agreement overall. Equilibrium bond distances are reproduced to within 0.05 %. The dissociation energies of the states agree with experiment to within similar to 0.5 kcal/mol, achieving "chemical accuracy." Vibrational energy levels show average deviations of similar to 20 cm(-1) or less. The B-1 Delta(g) state shows the best agreement with a mean absolute deviation of 2.41 cm(-1). Calculated rotational constants exhibit very good agreement with experiment, as do the spectroscopic constants.
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.
Shi, Deheng; Li, Peiling; Sun, Jinfeng; Zhu, Zunlue
2014-01-01
The potential energy curves (PECs) of 28 Ω states generated from 9 Λ-S states (X(2)Π, 1(4)Π, 1(6)Π, 1(2)Σ(+), 1(4)Σ(+), 1(6)Σ(+), 1(4)Σ(-), 2(4)Π and 1(4)Δ) are studied for the first time using an ab initio quantum chemical method. All the 9 Λ-S states correlate to the first two dissociation limits, N((4)Su)+Se((3)Pg) and N((4)Su)+Se((3)Dg), of NSe radical. Of these Λ-S states, the 1(6)Σ(+), 1(4)Σ(+), 1(6)Π, 2(4)Π and 1(4)Δ are found to be rather weakly bound states. The 1(2)Σ(+) is found to be unstable and has double wells. And the 1(6)Σ(+), 1(4)Σ(+), 1(4)Π and 1(6)Π are found to be the inverted ones with the SO coupling included. The PEC calculations are made by the complete active space self-consistent field method, which is followed by the internally contracted multireference configuration interaction approach with the Davidson modification. The spin-orbit coupling is accounted for by the state interaction approach with the Breit-Pauli Hamiltonian. The convergence of the present calculations is discussed with respect to the basis set and the level of theory. Core-valence correlation corrections are included with a cc-pCVTZ basis set. Scalar relativistic corrections are calculated by the third-order Douglas-Kroll Hamiltonian approximation at the level of a cc-pV5Z basis set. All the PECs are extrapolated to the complete basis set limit. The variation with internuclear separation of spin-orbit coupling constants is discussed in brief for some Λ-S states with one shallow well on each PEC. The spectroscopic parameters of 9 Λ-S and 28 Ω states are determined by fitting the first ten vibrational levels whenever available, which are calculated by solving the rovibrational Schrödinger equation with Numerov's method. The splitting energy in the X(2)Π Λ-S state is determined to be about 864.92 cm(-1), which agrees favorably with the measurements of 891.80 cm(-1). Moreover, other spectroscopic parameters of Λ-S and Ω states involved here are
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
Meng, Qingyong Chen, Jun Zhang, Dong H.
2015-09-14
The ring polymer molecular dynamics (RPMD) calculations are performed to calculate rate constants for the title reaction on the recently constructed potential energy surface based on permutation invariant polynomial (PIP) neural-network (NN) fitting [J. Li et al., J. Chem. Phys. 142, 204302 (2015)]. By inspecting convergence, 16 beads are used in computing free-energy barriers at 300 K ≤ T ≤ 1000 K, while different numbers of beads are used for transmission coefficients. The present RPMD rates are in excellent agreement with quantum rates computed on the same potential energy surface, as well as with the experimental measurements, demonstrating further that the RPMD is capable of producing accurate rates for polyatomic chemical reactions even at rather low temperatures.
Accurate momentum transfer cross section for the attractive Yukawa potential
Khrapak, S. A.
2014-04-15
Accurate expression for the momentum transfer cross section for the attractive Yukawa potential is proposed. This simple analytic expression agrees with the numerical results better than to within ±2% in the regime relevant for ion-particle collisions in complex (dusty) plasmas.
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)
Shi, De-Heng; Liu, Qionglan; Yu, Wei; Sun, Jinfeng; Zhu, Zunlue
2014-05-01
The potential energy curves (PECs) of 23 Ω states generated from the 12 electronic states (X1 Σ +, 21 Σ +, 11 Σ -, 11 Π, 21 Π, 11 Δ, 13 Σ +, 23 Σ +, 13 Σ -, a3 Π, 23 Π and 13 Δ) are studied for the first time. All the states correlate to the first dissociation channel of the SiBr+ cation. Of these electronic states, the 23 Σ + is the repulsive one without the spin-orbit coupling, whereas it becomes the bound one with the spin-orbit coupling added. On the one hand, without the spin-orbit coupling, the 11 Π, 21 Π and 23 Π are the rather weakly bound states, and only the 11 Π state possesses the double well; on the other hand, with the spin-orbit coupling included, the a3 Π and 11 Π states possess the double well, and the 13 Σ + and 13 Σ - are the inverted states. The PECs are calculated by the CASSCF method, which is followed by the internally contracted MRCI approach with the Davidson modification. Scalar relativistic correction is calculated by the third-order Douglas-Kroll Hamiltonian approximation with a cc-pVTZ-DK basis set. Core-valence correlation correction is included with a cc-pCVTZ basis set. The spin-orbit coupling is accounted for by the state interaction method with the Breit-Pauli Hamiltonian using the all-electron aug-cc-pCVTZ basis set. All the PECs are extrapolated to the complete basis set limit. The variation with internuclear separation of the spin-orbit coupling constant is discussed in brief. The spectroscopic parameters are evaluated for the 11 bound electronic states and the 23 bound Ω states, and are compared with available measurements. Excellent agreement has been found between the present results and the experimental data. It demonstrates that the spectroscopic parameters reported here can be expected to be reliably predicted ones. The Franck-Condon factors and radiative lifetimes of the transitions from the a3 Π 0 + and a3 Π 1 states to the X1 Σ + 0+ state are calculated for several low vibrational levels, and
Shi, Deheng; Liu, Qionglan; Sun, Jinfeng; Zhu, Zunlue
2014-03-25
The potential energy curves (PECs) of 28 Ω states generated from the 12 states (X(4)Σ(-), 1(2)Π, 1(2)Σ(-), 1(2)Δ, 1(2)Σ(+), 2(2)Π, A(4)Π, B(4)Σ(-), 3(2)Π, 1(6)Σ(-), 2(2)Σ(-) and 1(6)Π) of the BN(+) cation are studied for the first time for internuclear separations from about 0.1 to 1.0 nm using an ab initio quantum chemical method. All the Λ-S states correlate to the first four dissociation channels. The 1(6)Σ(-), 3(2)Π and A(4)Π states are found to be the inverted ones. The 1(2)Σ(+), 2(2)Π, 3(2)Π and 2(2)Σ(-) states are found to possess the double well. The PECs are calculated by the complete active space self-consistent field method, which is followed by the internally contracted multireference configuration interaction approach with the Davidson correction. Core-valence correlation correction is included by a cc-pCV5Z basis set. Scalar relativistic correction is calculated by the third-order Douglas-Kroll Hamiltonian approximation at the level of a cc-pV5Z basis set. The convergent behavior of present calculations is discussed with respect to the basis set and level of theory. The spin-orbit coupling is accounted for by the state interaction approach with the Breit-Pauli Hamiltonian using the all-electron cc-pCV5Z basis set. All the PECs are extrapolated to the complete basis set limit. The spectroscopic parameters are obtained, and the vibrational properties of 1(2)Σ(+), 2(2)Π, 3(2)Π and 2(2)Σ(-) states are evaluated. Analyses demonstrate that the spectroscopic parameters reported here can be expected to be reliably predicted ones. The conclusion is gained that the effect of spin-orbit coupling on the spectroscopic parameters are not obvious almost for all the Λ-S states involved in the present paper. PMID:24334021
Huang, Xinchuan E-mail: Timothy.J.Lee@nasa.gov; Schwenke, David W.; Lee, Timothy J. E-mail: Timothy.J.Lee@nasa.gov
2014-03-21
A purely ab initio potential energy surface (PES) was refined with selected {sup 32}S{sup 16}O{sub 2} HITRAN data. Compared to HITRAN, the root-mean-squares error (σ{sub RMS}) for all J = 0–80 rovibrational energy levels computed on the refined PES (denoted Ames-1) is 0.013 cm{sup −1}. Combined with a CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), an infrared (IR) line list (denoted Ames-296K) has been computed at 296 K and covers up to 8000 cm{sup −1}. Compared to the HITRAN and CDMS databases, the intensity agreement for most vibrational bands is better than 85%–90%. Our predictions for {sup 34}S{sup 16}O{sub 2} band origins, higher energy {sup 32}S{sup 16}O{sub 2} band origins and missing {sup 32}S{sup 16}O{sub 2} IR bands have been verified by most recent experiments and available HITRAN data. We conclude that the Ames-1 PES is able to predict {sup 32/34}S{sup 16}O{sub 2} band origins below 5500 cm{sup −1} with 0.01–0.03 cm{sup −1} uncertainties, and the Ames-296K line list provides continuous, reliable and accurate IR simulations. The K{sub a}-dependence of both line position and line intensity errors is discussed. The line list will greatly facilitate SO{sub 2} IR spectral experimental analysis, as well as elimination of SO{sub 2} lines in high-resolution astronomical observations.
Liu, Hui; Shi, Deheng; Sun, Jinfeng; Zhu, Zunlue; Shulin, Zhang
2014-04-24
The potential energy curves (PECs) of 54 spin-orbit states generated from the 22 electronic states of O2 molecule are investigated for the first time for internuclear separations from about 0.1 to 1.0nm. Of the 22 electronic states, the X(3)Σg(-), A(')(3)Δu, A(3)Σu(+), B(3)Σu(-), C(3)Πg, a(1)Δg, b(1)Σg(+), c(1)Σu(-), d(1)Πg, f(1)Σu(+), 1(5)Πg, 1(3)Πu, 2(3)Σg(-), 1(5)Σu(-), 2(1)Σu(-) and 2(1)Δg are found to be bound, whereas the 1(5)Σg(+), 2(5)Σg(+), 1(1)Πu, 1(5)Δg, 1(5)Πu and 2(1)Πu are found to be repulsive ones. The B(3)Σu(-) and d(1)Πg states possess the double well. And the 1(3)Πu, C(3)Πg, A'(3)Δu, 1(5)Δg and 2(5)Σg(+) states are the inverted ones when the spin-orbit coupling is included. The PEC calculations are done by the complete active space self-consistent field (CASSCF) method, which is followed by the internally contracted multireference configuration interaction (icMRCI) approach with the Davidson correction. Core-valence correlation and scalar relativistic corrections are taken into account. The convergence of present calculations is evaluated with respect to the basis set and level of theory. The vibrational properties are discussed for the 1(5)Πg, 1(3)Πu, d(1)Πg and 1(5)Σu(-) states and for the second well of the B(3)Σu(-) state. The spin-orbit coupling effect is accounted for by the state interaction method with the Breit-Pauli Hamiltonian. The PECs of all the electronic states and spin-orbit states are extrapolated to the complete basis set limit. The spectroscopic parameters are obtained, and compared with available experimental and other theoretical results. Analyses demonstrate that the spectroscopic parameters reported here can be expected to be reliably predicted ones. The conclusion is obtained that the effect of spin-orbit coupling on the spectroscopic parameters are small almost for all the electronic states involved in this paper except for the 1(5)Σu(-), 1(5)Πg and 1(3)Πu. PMID:24486866
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({sup 2}D) + H{sub 2} 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({sup 2}D)+H{sub 2}(X{sup 1}Σ{sub g}{sup +})(ν=0,j=0)→NH(a{sup 1}Δ)+H({sup 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.
Accurate ab initio energy gradients in chemical compound space.
Anatole von Lilienfeld, O
2009-10-28
Analytical potential energy derivatives, based on the Hellmann-Feynman theorem, are presented for any pair of isoelectronic compounds. Since energies are not necessarily monotonic functions between compounds, these derivatives can fail to predict the right trends of the effect of alchemical mutation. However, quantitative estimates without additional self-consistency calculations can be made when the Hellmann-Feynman derivative is multiplied with a linearization coefficient that is obtained from a reference pair of compounds. These results suggest that accurate predictions can be made regarding any molecule's energetic properties as long as energies and gradients of three other molecules have been provided. The linearization coefficent can be interpreted as a quantitative measure of chemical similarity. Presented numerical evidence includes predictions of electronic eigenvalues of saturated and aromatic molecular hydrocarbons. PMID:19894922
Accurate complex scaling of three dimensional numerical potentials
Cerioni, Alessandro; Genovese, Luigi; Duchemin, Ivan; Deutsch, Thierry
2013-05-28
The complex scaling method, which consists in continuing spatial coordinates into the complex plane, is a well-established method that allows to compute resonant eigenfunctions of the time-independent Schroedinger operator. Whenever it is desirable to apply the complex scaling to investigate resonances in physical systems defined on numerical discrete grids, the most direct approach relies on the application of a similarity transformation to the original, unscaled Hamiltonian. We show that such an approach can be conveniently implemented in the Daubechies wavelet basis set, featuring a very promising level of generality, high accuracy, and no need for artificial convergence parameters. Complex scaling of three dimensional numerical potentials can be efficiently and accurately performed. By carrying out an illustrative resonant state computation in the case of a one-dimensional model potential, we then show that our wavelet-based approach may disclose new exciting opportunities in the field of computational non-Hermitian quantum mechanics.
NASA Astrophysics Data System (ADS)
Garrett, Bruce C.; Truhlar, Donald G.
1980-03-01
Rate constants are calculated for the reactions D+H2→DH+H and H+D2→HD+D and compared to measured values. An accurate potential energy surface, based on the ab initio calculations of Liu and Siegbahn, was used. Rates were calculated using both conventional transition state theory and canonical variational theory. In the former, the generalized transition state dividing surface is located at the saddle point; in the latter it is located to maximize the generalized free energy of activation. We show that, in the absence of tunneling corrections, locating the generalized-transition-state dividing surface variationally has an important quantitative effect on the predicted rate constants for these systems and that, when tunneling is included, most of the effect of using a better dividing surface can be included in conventional transition state theory for these systems by using a consistent transmission coefficient for quantal scattering by the vibrationally adiabatic potential energy curve. Tunneling effects are important for these reactions even for temperatures larger than 400 K. We show how to separate classical recrossing effects from quantal corrections on reaction-coordinate motion in both the transmission coefficients and the kinetic isotope effects. Our most complete calculations are in excellent agreement with most of the measured rate constants and kinetic isotope effects.
Accurate calculation of diffraction-limited encircled and ensquared energy.
Andersen, Torben B
2015-09-01
Mathematical properties of the encircled and ensquared energy functions for the diffraction-limited point-spread function (PSF) are presented. These include power series and a set of linear differential equations that facilitate the accurate calculation of these functions. Asymptotic expressions are derived that provide very accurate estimates for the relative amount of energy in the diffraction PSF that fall outside a square or rectangular large detector. Tables with accurate values of the encircled and ensquared energy functions are also presented. PMID:26368873
NASA Astrophysics Data System (ADS)
Czakó, Gábor; Bowman, Joel M.
2012-01-01
We report a high-quality, ab initio, full-dimensional global potential energy surface (PES) for the Cl(2P, 2P3/2) + CH4 reaction, which describes both the abstraction (HCl + CH3) and substitution (H + CH3Cl) channels. The analytical PES is a least-squares fit, using a basis of permutationally invariant polynomials, to roughly 16 000 ab initio energy points, obtained by an efficient composite method, including counterpoise and spin-orbit corrections for the entrance channel. This composite method is shown to provide accuracy almost equal to all-electron CCSD(T)/aug-cc-pCVQZ results, but at much lower computational cost. Details of the PES, as well as additional high-level benchmark characterization of structures and energetics are reported. The PES has classical barrier heights of 2650 and 15 060 cm-1 (relative to Cl(2P3/2) + CH4(eq)), respectively, for the abstraction and substitution reactions, in good agreement with the corresponding new computed benchmark values, 2670 and 14 720 cm-1. The PES also accurately describes the potential wells in the entrance and exit channels for the abstraction reaction. Quasiclassical trajectory calculations using the PES show that (a) the inclusion of the spin-orbit corrections in the PES decreases the cross sections by a factor of 1.5-2.5 at low collision energies (Ecoll); (b) at Ecoll ≈ 13 000 cm-1 the substitution channel opens and the H/HCl ratio increases rapidly with Ecoll; (c) the maximum impact parameter (bmax) for the abstraction reaction is ˜6 bohr; whereas bmax is only ˜2 bohr for the substitution; (d) the HCl and CH3 products are mainly in the vibrational ground state even at very high Ecoll; and (e) the HCl rotational distributions are cold, in excellent agreement with experiment at Ecoll = 1280 cm-1.
Accurate ab initio vibrational energies of methyl chloride
Owens, Alec; Yurchenko, Sergei N.; Yachmenev, Andrey; Tennyson, Jonathan; Thiel, Walter
2015-06-28
Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH{sub 3}{sup 35}Cl and CH{sub 3}{sup 37}Cl. The respective PESs, CBS-35{sup HL}, and CBS-37{sup HL}, are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY {sub 3}Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35{sup HL} and CBS-37{sup HL} PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm{sup −1}, respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH{sub 3}Cl without empirical refinement of the respective PESs.
Accurate ab initio vibrational energies of methyl chloride.
Owens, Alec; Yurchenko, Sergei N; Yachmenev, Andrey; Tennyson, Jonathan; Thiel, Walter
2015-06-28
Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH3 (35)Cl and CH3 (37)Cl. The respective PESs, CBS-35( HL), and CBS-37( HL), are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY 3Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35( HL) and CBS-37( HL) PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm(-1), respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH3Cl without empirical refinement of the respective PESs. PMID:26133427
Accurate ab initio vibrational energies of methyl chloride
NASA Astrophysics Data System (ADS)
Owens, Alec; Yurchenko, Sergei N.; Yachmenev, Andrey; Tennyson, Jonathan; Thiel, Walter
2015-06-01
Two new nine-dimensional potential energy surfaces (PESs) have been generated using high-level ab initio theory for the two main isotopologues of methyl chloride, CH335Cl and CH337Cl. The respective PESs, CBS-35 HL, and CBS-37 HL, are based on explicitly correlated coupled cluster calculations with extrapolation to the complete basis set (CBS) limit, and incorporate a range of higher-level (HL) additive energy corrections to account for core-valence electron correlation, higher-order coupled cluster terms, scalar relativistic effects, and diagonal Born-Oppenheimer corrections. Variational calculations of the vibrational energy levels were performed using the computer program TROVE, whose functionality has been extended to handle molecules of the form XY 3Z. Fully converged energies were obtained by means of a complete vibrational basis set extrapolation. The CBS-35 HL and CBS-37 HL PESs reproduce the fundamental term values with root-mean-square errors of 0.75 and 1.00 cm-1, respectively. An analysis of the combined effect of the HL corrections and CBS extrapolation on the vibrational wavenumbers indicates that both are needed to compute accurate theoretical results for methyl chloride. We believe that it would be extremely challenging to go beyond the accuracy currently achieved for CH3Cl without empirical refinement of the respective PESs.
Chemically accurate energy barriers of small gas molecules moving through hexagonal water rings.
Hjertenæs, Eirik; Trinh, Thuat T; Koch, Henrik
2016-07-21
We present chemically accurate potential energy curves of CH4, CO2 and H2 moving through hexagonal water rings, calculated by CCSD(T)/aug-cc-pVTZ with counterpoise correction. The barriers are extracted from a potential energy surface obtained by allowing the water ring to expand while the gas molecule diffuses through. State-of-the-art XC-functionals are evaluated against the CCSD(T) potential energy surface. PMID:27345929
Szabó, István; Telekes, Hajnalka; Czakó, Gábor
2015-06-28
We develop a full-dimensional global analytical potential energy surface (PES) for the F{sup −} + CH{sub 3}F reaction by fitting about 50 000 energy points obtained by an explicitly correlated composite method based on the second-order Møller–Plesset perturbation-F12 and coupled-cluster singles, doubles, and perturbative triples-F12a methods and the cc-pVnZ-F12 [n = D, T] basis sets. The PES accurately describes the (a) back-side attack Walden inversion mechanism involving the pre- and post-reaction (b) ion-dipole and (c) hydrogen-bonded complexes, the configuration-retaining (d) front-side attack and (e) double-inversion substitution pathways, as well as (f) the proton-abstraction channel. The benchmark quality relative energies of all the important stationary points are computed using the focal-point analysis (FPA) approach considering electron correlation up to coupled-cluster singles, doubles, triples, and perturbative quadruples method, extrapolation to the complete basis set limit, core-valence correlation, and scalar relativistic effects. The FPA classical(adiabatic) barrier heights of (a), (d), and (e) are −0.45(−0.61), 46.07(45.16), and 29.18(26.07) kcal mol{sup −1}, respectively, the dissociation energies of (b) and (c) are 13.81(13.56) and 13.73(13.52) kcal mol{sup −1}, respectively, and the endothermicity of (f) is 42.54(38.11) kcal mol{sup −1}. Quasiclassical trajectory computations of cross sections, scattering (θ) and initial attack (α) angle distributions, as well as translational and internal energy distributions are performed for the F{sup −} + CH{sub 3}F(v = 0) reaction using the new PES. Apart from low collision energies (E{sub coll}), the S{sub N}2 excitation function is nearly constant, the abstraction cross sections rapidly increase with E{sub coll} from a threshold of ∼40 kcal mol{sup −1}, and retention trajectories via double inversion are found above E{sub coll} = ∼ 30 kcal mol{sup −1}, and at E{sub coll} =
Thompson, A.P.; Swiler, L.P.; Trott, C.R.; Foiles, S.M.; Tucker, G.J.
2015-03-15
We present a new interatomic potential for solids and liquids called Spectral Neighbor Analysis Potential (SNAP). The SNAP potential has a very general form and uses machine-learning techniques to reproduce the energies, forces, and stress tensors of a large set of small configurations of atoms, which are obtained using high-accuracy quantum electronic structure (QM) calculations. The local environment of each atom is characterized by a set of bispectrum components of the local neighbor density projected onto a basis of hyperspherical harmonics in four dimensions. The bispectrum components are the same bond-orientational order parameters employed by the GAP potential [1]. The SNAP potential, unlike GAP, assumes a linear relationship between atom energy and bispectrum components. The linear SNAP coefficients are determined using weighted least-squares linear regression against the full QM training set. This allows the SNAP potential to be fit in a robust, automated manner to large QM data sets using many bispectrum components. The calculation of the bispectrum components and the SNAP potential are implemented in the LAMMPS parallel molecular dynamics code. We demonstrate that a previously unnoticed symmetry property can be exploited to reduce the computational cost of the force calculations by more than one order of magnitude. We present results for a SNAP potential for tantalum, showing that it accurately reproduces a range of commonly calculated properties of both the crystalline solid and the liquid phases. In addition, unlike simpler existing potentials, SNAP correctly predicts the energy barrier for screw dislocation migration in BCC tantalum.
Accurate nuclear radii and binding energies from a chiral interaction
Ekstrom, Jan A.; Jansen, G. R.; Wendt, Kyle A.; Hagen, Gaute; Papenbrock, Thomas F.; Carlsson, Boris; Forssen, Christian; Hjorth-Jensen, M.; Navratil, Petr; Nazarewicz, Witold
2015-05-01
With the goal of developing predictive ab initio capability for light and medium-mass nuclei, two-nucleon and three-nucleon forces from chiral effective field theory are optimized simultaneously to low-energy nucleon-nucleon scattering data, as well as binding energies and radii of few-nucleon systems and selected isotopes of carbon and oxygen. Coupled-cluster calculations based on this interaction, named NNLOsat, yield accurate binding energies and radii of nuclei up to 40Ca, and are consistent with the empirical saturation point of symmetric nuclear matter. In addition, the low-lying collective Jπ=3- states in 16O and 40Ca are described accurately, while spectra for selected p- and sd-shellmore » nuclei are in reasonable agreement with experiment.« less
Accurate nuclear radii and binding energies from a chiral interaction
Ekstrom, Jan A.; Jansen, G. R.; Wendt, Kyle A.; Hagen, Gaute; Papenbrock, Thomas F.; Carlsson, Boris; Forssen, Christian; Hjorth-Jensen, M.; Navratil, Petr; Nazarewicz, Witold
2015-05-01
With the goal of developing predictive ab initio capability for light and medium-mass nuclei, two-nucleon and three-nucleon forces from chiral effective field theory are optimized simultaneously to low-energy nucleon-nucleon scattering data, as well as binding energies and radii of few-nucleon systems and selected isotopes of carbon and oxygen. Coupled-cluster calculations based on this interaction, named NNLO_{sat}, yield accurate binding energies and radii of nuclei up to ^{40}Ca, and are consistent with the empirical saturation point of symmetric nuclear matter. In addition, the low-lying collective J^{π}=3^{-} states in ^{16}O and ^{40}Ca are described accurately, while spectra for selected p- and sd-shell nuclei are in reasonable agreement with experiment.
Highly accurate eigenvalues for the distorted Coulomb potential
NASA Astrophysics Data System (ADS)
Ixaru, L. Gr.; de Meyer, H.; vanden Berghe, G.
2000-03-01
We consider the eigenvalue problem for the radial Schrödinger equation with potentials of the form V(r)=S(r)/r+R(r) where S(r) and R(r) are well behaved functions which tend to some (not necessarily equal) constants when r-->0 and r-->∞. Formulas (14.4.5)-(14.4.8) of Abramowitz and Stegun [Handbook of Mathematical Functions, 8th ed. (Dover, New York, 1972)], corresponding to the pure Coulomb case, are here generalized for this distorted case. We also present a complete procedure for the numerical solution of the problem. Our procedure is robust, very economic and particularly suited for very large n. Numerical illustrations for n up to 2000 are given.
A fast, time-accurate unsteady full potential scheme
NASA Technical Reports Server (NTRS)
Shankar, V.; Ide, H.; Gorski, J.; Osher, S.
1985-01-01
The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi, obtained by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.
Development and application of accurate analytical models for single active electron potentials
NASA Astrophysics Data System (ADS)
Miller, Michelle; Jaron-Becker, Agnieszka; Becker, Andreas
2015-05-01
The single active electron (SAE) approximation is a theoretical model frequently employed to study scenarios in which inner-shell electrons may productively be treated as frozen spectators to a physical process of interest, and accurate analytical approximations for these potentials are sought as a useful simulation tool. Density function theory is often used to construct a SAE potential, requiring that a further approximation for the exchange correlation functional be enacted. In this study, we employ the Krieger, Li, and Iafrate (KLI) modification to the optimized-effective-potential (OEP) method to reduce the complexity of the problem to the straightforward solution of a system of linear equations through simple arguments regarding the behavior of the exchange-correlation potential in regions where a single orbital dominates. We employ this method for the solution of atomic and molecular potentials, and use the resultant curve to devise a systematic construction for highly accurate and useful analytical approximations for several systems. Supported by the U.S. Department of Energy (Grant No. DE-FG02-09ER16103), and the U.S. National Science Foundation (Graduate Research Fellowship, Grants No. PHY-1125844 and No. PHY-1068706).
Toward Hamiltonian Adaptive QM/MM: Accurate Solvent Structures Using Many-Body Potentials.
Boereboom, Jelle M; Potestio, Raffaello; Donadio, Davide; Bulo, Rosa E
2016-08-01
Adaptive quantum mechanical (QM)/molecular mechanical (MM) methods enable efficient molecular simulations of chemistry in solution. Reactive subregions are modeled with an accurate QM potential energy expression while the rest of the system is described in a more approximate manner (MM). As solvent molecules diffuse in and out of the reactive region, they are gradually included into (and excluded from) the QM expression. It would be desirable to model such a system with a single adaptive Hamiltonian, but thus far this has resulted in distorted structures at the boundary between the two regions. Solving this long outstanding problem will allow microcanonical adaptive QM/MM simulations that can be used to obtain vibrational spectra and dynamical properties. The difficulty lies in the complex QM potential energy expression, with a many-body expansion that contains higher order terms. Here, we outline a Hamiltonian adaptive multiscale scheme within the framework of many-body potentials. The adaptive expressions are entirely general, and complementary to all standard (nonadaptive) QM/MM embedding schemes available. We demonstrate the merit of our approach on a molecular system defined by two different MM potentials (MM/MM'). For the long-range interactions a numerical scheme is used (particle mesh Ewald), which yields energy expressions that are many-body in nature. Our Hamiltonian approach is the first to provide both energy conservation and the correct solvent structure everywhere in this system. PMID:27332140
Śmiga, Szymon; Della Sala, Fabio; Buksztel, Adam; Grabowski, Ireneusz; Fabiano, Eduardo
2016-08-15
One important property of Kohn-Sham (KS) density functional theory is the exact equality of the energy of the highest occupied KS orbital (HOMO) with the negative ionization potential of the system. This exact feature is out of reach for standard density-dependent semilocal functionals. Conversely, accurate results can be obtained using orbital-dependent functionals in the optimized effective potential (OEP) approach. In this article, we investigate the performance, in this context, of some advanced OEP methods, with special emphasis on the recently proposed scaled-opposite-spin OEP functional. Moreover, we analyze the impact of the so-called HOMO condition on the final quality of the HOMO energy. Results are compared to reference data obtained at the CCSD(T) level of theory. © 2016 Wiley Periodicals, Inc. PMID:27357413
Inverter Modeling For Accurate Energy Predictions Of Tracking HCPV Installations
NASA Astrophysics Data System (ADS)
Bowman, J.; Jensen, S.; McDonald, Mark
2010-10-01
High efficiency high concentration photovoltaic (HCPV) solar plants of megawatt scale are now operational, and opportunities for expanded adoption are plentiful. However, effective bidding for sites requires reliable prediction of energy production. HCPV module nameplate power is rated for specific test conditions; however, instantaneous HCPV power varies due to site specific irradiance and operating temperature, and is degraded by soiling, protective stowing, shading, and electrical connectivity. These factors interact with the selection of equipment typically supplied by third parties, e.g., wire gauge and inverters. We describe a time sequence model accurately accounting for these effects that predicts annual energy production, with specific reference to the impact of the inverter on energy output and interactions between system-level design decisions and the inverter. We will also show two examples, based on an actual field design, of inverter efficiency calculations and the interaction between string arrangements and inverter selection.
Accurate energy levels for singly ionized platinum (Pt II)
NASA Technical Reports Server (NTRS)
Reader, Joseph; Acquista, Nicolo; Sansonetti, Craig J.; Engleman, Rolf, Jr.
1988-01-01
New observations of the spectrum of Pt II have been made with hollow-cathode lamps. The region from 1032 to 4101 A was observed photographically with a 10.7-m normal-incidence spectrograph. The region from 2245 to 5223 A was observed with a Fourier-transform spectrometer. Wavelength measurements were made for 558 lines. The uncertainties vary from 0.0005 to 0.004 A. From these measurements and three parity-forbidden transitions in the infrared, accurate values were determined for 28 even and 72 odd energy levels of Pt II.
Accurate bond dissociation energies (D 0) for FHF- isotopologues
NASA Astrophysics Data System (ADS)
Stein, Christopher; Oswald, Rainer; Sebald, Peter; Botschwina, Peter; Stoll, Hermann; Peterson, Kirk A.
2013-09-01
Accurate bond dissociation energies (D 0) are determined for three isotopologues of the bifluoride ion (FHF-). While the zero-point vibrational contributions are taken from our previous work (P. Sebald, A. Bargholz, R. Oswald, C. Stein, P. Botschwina, J. Phys. Chem. A, DOI: 10.1021/jp3123677), the equilibrium dissociation energy (D e ) of the reaction ? was obtained by a composite method including frozen-core (fc) CCSD(T) calculations with basis sets up to cardinal number n = 7 followed by extrapolation to the complete basis set limit. Smaller terms beyond fc-CCSD(T) cancel each other almost completely. The D 0 values of FHF-, FDF-, and FTF- are predicted to be 15,176, 15,191, and 15,198 cm-1, respectively, with an uncertainty of ca. 15 cm-1.
Spectroscopically Accurate Calculations of the Rovibrational Energies of Diatomic Hydrogen
NASA Astrophysics Data System (ADS)
Perry, Jason
2005-05-01
The Born-Oppenheimer approximation has been used to calculate the rotational and vibrational states of diatomic hydrogen. Because it is an approximation, our group now wants to use a Born-Oppenheimer potential to calculate the electronic energy that has been corrected to match closely with spectroscopic results. We are using a code that has corrections for adiabatic, relativistic, radiative, and non-adiabatic effects. The rovibrational energies have now been calculated for both bound and quasi-bound states. We also want to compute quadrupole transition probabilities for diatomic hydrogen. These calculations aspire to investigate diatomic hydrogen in astrophysical environments.
Golibrzuch, Kai; Shirhatti, Pranav R.; Kandratsenka, Alexander; Wodtke, Alec M.; Bartels, Christof; Max Planck Institute for Biophysical Chemistry, Göttingen 37077 ; Rahinov, Igor; Auerbach, Daniel J.; Max Planck Institute for Biophysical Chemistry, Göttingen 37077; Department of Chemistry and Biochemistry, University of California Santa Barbara, Santa Barbara, California 93106
2014-01-28
We present a combined experimental and theoretical study of NO(v = 3 → 3, 2, 1) scattering from a Au(111) surface at incidence translational energies ranging from 0.1 to 1.2 eV. Experimentally, molecular beam–surface scattering is combined with vibrational overtone pumping and quantum-state selective detection of the recoiling molecules. Theoretically, we employ a recently developed first-principles approach, which employs an Independent Electron Surface Hopping (IESH) algorithm to model the nonadiabatic dynamics on a Newns-Anderson Hamiltonian derived from density functional theory. This approach has been successful when compared to previously reported NO/Au scattering data. The experiments presented here show that vibrational relaxation probabilities increase with incidence energy of translation. The theoretical simulations incorrectly predict high relaxation probabilities at low incidence translational energy. We show that this behavior originates from trajectories exhibiting multiple bounces at the surface, associated with deeper penetration and favored (N-down) molecular orientation, resulting in a higher average number of electronic hops and thus stronger vibrational relaxation. The experimentally observed narrow angular distributions suggest that mainly single-bounce collisions are important. Restricting the simulations by selecting only single-bounce trajectories improves agreement with experiment. The multiple bounce artifacts discovered in this work are also present in simulations employing electronic friction and even for electronically adiabatic simulations, meaning they are not a direct result of the IESH algorithm. This work demonstrates how even subtle errors in the adiabatic interaction potential, especially those that influence the interaction time of the molecule with the surface, can lead to an incorrect description of electronically nonadiabatic vibrational energy transfer in molecule-surface collisions.
NASA Astrophysics Data System (ADS)
Thompson, Aidan; Foiles, Stephen; Schultz, Peter; Swiler, Laura; Trott, Christian; Tucker, Garritt
2013-03-01
Molecular dynamics (MD) is a powerful condensed matter simulation tool for bridging between macroscopic continuum models and quantum models (QM) treating a few hundred atoms, but is limited by the accuracy of available interatomic potentials. Sound physical and chemical understanding of these interactions have resulted in a variety of concise potentials for certain systems, but it is difficult to extend them to new materials and properties. The growing availability of large QM data sets has made it possible to use more automated machine-learning approaches. Bartók et al. demonstrated that the bispectrum of the local neighbor density provides good regression surrogates for QM models. We adopt a similar bispectrum representation within a linear regression scheme. We have produced potentials for silicon and tantalum, and we are currently extending the method to III-V compounds. Results will be presented demonstrating the accuracy of these potentials relative to the training data, as well as their ability to accurately predict material properties not explicitly included in the training data. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Dept. of Energy Nat. Nuclear Security Admin. under Contract DE-AC04-94AL85000.
NASA Astrophysics Data System (ADS)
Colalongo, Luigi; Ghittorelli, Matteo; Torricelli, Fabrizio; Kovács-Vajna, Zsolt Miklos
2015-12-01
Surface-potential-based mathematical models are among the most accurate and physically based compact models of Thin-Film Transistors (TFTs) and, in turn, of Organic Thin-Film Transistors (OTFTs), available today. However, the need for iterative computations of the surface potential limits their computational efficiency and diffusion in CAD applications. The existing closed-form approximations of the surface potential are based on regional approximations and empirical smoothing functions that could result not enough accurate to model OTFTs and, in particular, transconductances and transcapacitances. In this paper we present an accurate and computationally efficient closed-form approximation of the surface potential, based on the Lagrange Reversion Theorem, that can be exploited in advanced surface-potential-based OTFTs and TFTs device models.
Navratil, P; Caurier, E
2003-10-14
The authors calculate properties of A = 6 system using the accurate charge-dependent nucleon-nucleon (NN) potential at fourth order of chiral perturbation theory. By application of the ab initio no-core shell model (NCSM) and a variational calculation in the harmonic oscillator basis with basis size up to 16 {h_bar}{Omega} they obtain the {sup 6}Li binding energy of 28.5(5) MeV and a converged excitation spectrum. Also, they calculate properties of {sup 10}B using the same NN potential in a basis space of up to 8 {h_bar}{Omega}. The results are consistent with results obtained by standard accurate NN potentials and demonstrate a deficiency of Hamiltonians consisting of only two-body terms. At this order of chiral perturbation theory three-body terms appear. It is expected that inclusion of such terms in the Hamiltonian will improve agreement with experiment.
The Calculation of Accurate Metal-Ligand Bond Energies
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W.; Partridge, Harry, III; Ricca, Alessandra; Arnold, James O. (Technical Monitor)
1997-01-01
The optimization of the geometry and calculation of zero-point energies are carried out at the B3LYP level of theory. The bond energies are determined at this level, as well as at the CCSD(T) level using very large basis sets. The successive OH bond energies to the first row transition metal cations are reported. For most systems there has been an experimental determination of the first OH. In general, the CCSD(T) values are in good agreement with experiment. The bonding changes from mostly covalent for the early metals to mostly electrostatic for the late transition metal systems.
High energy laser testbed for accurate beam pointing control
NASA Astrophysics Data System (ADS)
Kim, Dojong; Kim, Jae Jun; Frist, Duane; Nagashima, Masaki; Agrawal, Brij
2010-02-01
Precision laser beam pointing is a key technology in High Energy Laser systems. In this paper, a laboratory High Energy Laser testbed developed at the Naval Postgraduate School is introduced. System identification is performed and a mathematical model is constructed to estimate system performance. New beam pointing control algorithms are designed based on this mathematical model. It is shown in both computer simulation and experiment that the adaptive filter algorithm can improve the pointing performance of the system.
Alaska's renewable energy potential.
Not Available
2009-02-01
This paper delivers a brief survey of renewable energy technologies applicable to Alaska's climate, latitude, geography, and geology. We first identify Alaska's natural renewable energy resources and which renewable energy technologies would be most productive. e survey the current state of renewable energy technologies and research efforts within the U.S. and, where appropriate, internationally. We also present information on the current state of Alaska's renewable energy assets, incentives, and commercial enterprises. Finally, we escribe places where research efforts at Sandia National Laboratories could assist the state of Alaska with its renewable energy technology investment efforts.
Baer, M.R.; Hobbs, M.L.; McGee, B.C.
1998-11-03
Exponential-13,6 (EXP-13,6) potential pammeters for 750 gases composed of 48 elements were determined and assembled in a database, referred to as the JCZS database, for use with the Jacobs Cowperthwaite Zwisler equation of state (JCZ3-EOS)~l) The EXP- 13,6 force constants were obtained by using literature values of Lennard-Jones (LJ) potential functions, by using corresponding states (CS) theory, by matching pure liquid shock Hugoniot data, and by using molecular volume to determine the approach radii with the well depth estimated from high-pressure isen- tropes. The JCZS database was used to accurately predict detonation velocity, pressure, and temperature for 50 dif- 3 Accurate predictions were also ferent explosives with initial densities ranging from 0.25 glcm3 to 1.97 g/cm . obtained for pure liquid shock Hugoniots, static properties of nitrogen, and gas detonations at high initial pressures.
Accurate energies of the He atom with undergraduate quantum mechanics
NASA Astrophysics Data System (ADS)
Massé, Robert C.; Walker, Thad G.
2015-08-01
Estimating the energies and splitting of the 1s2s singlet and triplet states of helium is a classic exercise in quantum perturbation theory but yields only qualitatively correct results. Using a six-line computer program, the 1s2s energies calculated by matrix diagonalization using a seven-state basis improve the results to 0.4% error or better. This is an effective and practical illustration of the quantitative power of quantum mechanics, at a level accessible to undergraduate students.
Accurate Mass Determinations in Decay Chains with Missing Energy
Cheng, H.-C; Gunion, John F.; Han Zhenyu; Engelhardt, Dalit; McElrath, Bob
2008-06-27
Many beyond the standard model theories include a stable dark matter candidate that yields missing or invisible energy in collider detectors. If observed at the CERN Large Hadron Collider, we must determine if its mass and other properties (and those of its partners) predict the correct dark matter relic density. We give a new procedure for determining its mass with small error.
Accurate mass determinations in decay chains with missing energy.
Cheng, Hsin-Chia; Engelhardt, Dalit; Gunion, John F; Han, Zhenyu; McElrath, Bob
2008-06-27
Many beyond the standard model theories include a stable dark matter candidate that yields missing or invisible energy in collider detectors. If observed at the CERN Large Hadron Collider, we must determine if its mass and other properties (and those of its partners) predict the correct dark matter relic density. We give a new procedure for determining its mass with small error. PMID:18643654
Waldrop, Jonathan M; Song, Bo; Patkowski, Konrad; Wang, Xiaopo
2015-05-28
A new highly accurate potential energy curve for the krypton dimer was constructed using coupled-cluster calculations up to the singles, doubles, triples, and perturbative quadruples level, including corrections for core-core and core-valence correlation and for relativistic effects. The ab initio data points were fitted to an analytic potential which was used to compute the most important transport properties of the krypton gas. The viscosity, thermal conductivity, self-diffusion coefficient, and thermal diffusion factor were calculated by the kinetic theory at low density and temperatures from 116 to 5000 K. The comparisons with literature experimental data as well as with values from other pair potentials indicate that our new potential is superior to all previous ones. The transport property values computed in this work are recommended as standard values over the complete temperature range. PMID:26026447
The free energy cost of accurate biochemical oscillations
Cao, Yuansheng; Wang, Hongli; Ouyang, Qi; Tu, Yuhai
2015-01-01
Oscillation is an important cellular process that regulates timing of different vital life cycles. However, in the noisy cellular environment, oscillations can be highly inaccurate due to phase fluctuations. It remains poorly understood how biochemical circuits suppress phase fluctuations and what is the incurred thermodynamic cost. Here, we study three different types of biochemical oscillations representing three basic oscillation motifs shared by all known oscillatory systems. In all the systems studied, we find that the phase diffusion constant depends on the free energy dissipation per period following the same inverse relation parameterized by system specific constants. This relationship and its range of validity are shown analytically in a model of noisy oscillation. Microscopically, we find that the oscillation is driven by multiple irreversible cycles that hydrolyze the fuel molecules such as ATP; the number of phase coherent periods is proportional to the free energy consumed per period. Experimental evidence in support of this general relationship and testable predictions are also presented. PMID:26566392
Hansen, Katja; Biegler, Franziska; Ramakrishnan, Raghunathan; Pronobis, Wiktor; von Lilienfeld, O. Anatole; Müller, Klaus -Robert; Tkatchenko, Alexandre
2015-06-04
Simultaneously accurate and efficient prediction of molecular properties throughout chemical compound space is a critical ingredient toward rational compound design in chemical and pharmaceutical industries. Aiming toward this goal, we develop and apply a systematic hierarchy of efficient empirical methods to estimate atomization and total energies of molecules. These methods range from a simple sum over atoms, to addition of bond energies, to pairwise interatomic force fields, reaching to the more sophisticated machine learning approaches that are capable of describing collective interactions between many atoms or bonds. In the case of equilibrium molecular geometries, even simple pairwise force fields demonstratemore » prediction accuracy comparable to benchmark energies calculated using density functional theory with hybrid exchange-correlation functionals; however, accounting for the collective many-body interactions proves to be essential for approaching the “holy grail” of chemical accuracy of 1 kcal/mol for both equilibrium and out-of-equilibrium geometries. This remarkable accuracy is achieved by a vectorized representation of molecules (so-called Bag of Bonds model) that exhibits strong nonlocality in chemical space. The same representation allows us to predict accurate electronic properties of molecules, such as their polarizability and molecular frontier orbital energies.« less
Hansen, Katja; Biegler, Franziska; Ramakrishnan, Raghunathan; Pronobis, Wiktor; von Lilienfeld, O. Anatole; Müller, Klaus -Robert; Tkatchenko, Alexandre
2015-06-04
Simultaneously accurate and efficient prediction of molecular properties throughout chemical compound space is a critical ingredient toward rational compound design in chemical and pharmaceutical industries. Aiming toward this goal, we develop and apply a systematic hierarchy of efficient empirical methods to estimate atomization and total energies of molecules. These methods range from a simple sum over atoms, to addition of bond energies, to pairwise interatomic force fields, reaching to the more sophisticated machine learning approaches that are capable of describing collective interactions between many atoms or bonds. In the case of equilibrium molecular geometries, even simple pairwise force fields demonstrate prediction accuracy comparable to benchmark energies calculated using density functional theory with hybrid exchange-correlation functionals; however, accounting for the collective many-body interactions proves to be essential for approaching the “holy grail” of chemical accuracy of 1 kcal/mol for both equilibrium and out-of-equilibrium geometries. This remarkable accuracy is achieved by a vectorized representation of molecules (so-called Bag of Bonds model) that exhibits strong nonlocality in chemical space. The same representation allows us to predict accurate electronic properties of molecules, such as their polarizability and molecular frontier orbital energies.
Theoretical studies of potential energy surfaces
Harding, L.B.
1993-12-01
The goal of this program is to calculate accurate potential energy surfaces (PES) for both reactive and nonreactive systems. To do this the electronic Schrodinger equation must be solved. Our approach to this problem starts with multiconfiguration self-consistent field (MCSCF) reference wavefunctions. These reference wavefunctions are designed to be sufficiently flexible to accurately describe changes in electronic structure over a broad range of geometries. Electron correlation effects are included via multireference, singles and doubles configuration interaction (MRSDCI) calculations. With this approach, the authors are able to provide useful predictions of the energetics for a broad range of systems.
Computed potential energy surfaces for chemical reactions
NASA Technical Reports Server (NTRS)
Walch, Stephen P.
1990-01-01
The objective was to obtain accurate potential energy surfaces (PES's) for a number of reactions which are important in the H/N/O combustion process. The interest in this is centered around the design of the SCRAM jet engine for the National Aerospace Plane (NASP), which was envisioned as an air-breathing hydrogen-burning vehicle capable of reaching velocities as large as Mach 25. Preliminary studies indicated that the supersonic flow in the combustor region of the scram jet engine required accurate reaction rate data for reactions in the H/N/O system, some of which was not readily available from experiment. The most important class of combustion reactions from the standpoint of the NASP project are radical recombinaton reactions, since these reactions result in most of the heat release in the combustion process. Theoretical characterizations of the potential energy surfaces for these reactions are presented and discussed.
Taylor, Mark P; Ye, Yuting; Adhikari, Shishir R
2015-11-28
The conformation of a polymer chain in solution is coupled to the local structure of the surrounding solvent and can undergo large changes in response to variations in solvent density and temperature. The many-body effects of solvent on the structure of an n-mer polymer chain can be formally mapped to an exact n-body solvation potential. Here, we use a pair decomposition of this n-body potential to construct a set of two-body potentials for a Lennard-Jones (LJ) polymer chain in explicit LJ solvent. The solvation potentials are built from numerically exact results for 5-mer chains in solvent combined with an approximate asymptotic expression for the solvation potential between sites that are distant along the chain backbone. These potentials map the many-body chain-in-solvent problem to a few-body single-chain problem and can be used to study a chain of arbitrary length, thereby dramatically reducing the computational complexity of the polymer chain-in-solvent problem. We have constructed solvation potentials at a large number of state points across the LJ solvent phase diagram including the vapor, liquid, and super-critical regions. We use these solvation potentials in single-chain Monte Carlo (MC) simulations with n ≤ 800 to determine the size, intramolecular structure, and scaling behavior of chains in solvent. To assess our results, we have carried out full chain-in-solvent MC simulations (with n ≤ 100) and find that our solvation potential approach is quantitatively accurate for a wide range of solvent conditions for these chain lengths. PMID:26627969
Accurate potential drop sheet resistance measurements of laser-doped areas in semiconductors
Heinrich, Martin; Kluska, Sven; Binder, Sebastian; Hameiri, Ziv; Hoex, Bram; Aberle, Armin G.
2014-10-07
It is investigated how potential drop sheet resistance measurements of areas formed by laser-assisted doping in crystalline Si wafers are affected by typically occurring experimental factors like sample size, inhomogeneities, surface roughness, or coatings. Measurements are obtained with a collinear four point probe setup and a modified transfer length measurement setup to measure sheet resistances of laser-doped lines. Inhomogeneities in doping depth are observed from scanning electron microscope images and electron beam induced current measurements. It is observed that influences from sample size, inhomogeneities, surface roughness, and coatings can be neglected if certain preconditions are met. Guidelines are given on how to obtain accurate potential drop sheet resistance measurements on laser-doped regions.
Intermolecular potentials and the accurate prediction of the thermodynamic properties of water.
Shvab, I; Sadus, Richard J
2013-11-21
The ability of intermolecular potentials to correctly predict the thermodynamic properties of liquid water at a density of 0.998 g∕cm(3) for a wide range of temperatures (298-650 K) and pressures (0.1-700 MPa) is investigated. Molecular dynamics simulations are reported for the pressure, thermal pressure coefficient, thermal expansion coefficient, isothermal and adiabatic compressibilities, isobaric and isochoric heat capacities, and Joule-Thomson coefficient of liquid water using the non-polarizable SPC∕E and TIP4P∕2005 potentials. The results are compared with both experiment data and results obtained from the ab initio-based Matsuoka-Clementi-Yoshimine non-additive (MCYna) [J. Li, Z. Zhou, and R. J. Sadus, J. Chem. Phys. 127, 154509 (2007)] potential, which includes polarization contributions. The data clearly indicate that both the SPC∕E and TIP4P∕2005 potentials are only in qualitative agreement with experiment, whereas the polarizable MCYna potential predicts some properties within experimental uncertainty. This highlights the importance of polarizability for the accurate prediction of the thermodynamic properties of water, particularly at temperatures beyond 298 K. PMID:24320337
Intermolecular potentials and the accurate prediction of the thermodynamic properties of water
NASA Astrophysics Data System (ADS)
Shvab, I.; Sadus, Richard J.
2013-11-01
The ability of intermolecular potentials to correctly predict the thermodynamic properties of liquid water at a density of 0.998 g/cm3 for a wide range of temperatures (298-650 K) and pressures (0.1-700 MPa) is investigated. Molecular dynamics simulations are reported for the pressure, thermal pressure coefficient, thermal expansion coefficient, isothermal and adiabatic compressibilities, isobaric and isochoric heat capacities, and Joule-Thomson coefficient of liquid water using the non-polarizable SPC/E and TIP4P/2005 potentials. The results are compared with both experiment data and results obtained from the ab initio-based Matsuoka-Clementi-Yoshimine non-additive (MCYna) [J. Li, Z. Zhou, and R. J. Sadus, J. Chem. Phys. 127, 154509 (2007)] potential, which includes polarization contributions. The data clearly indicate that both the SPC/E and TIP4P/2005 potentials are only in qualitative agreement with experiment, whereas the polarizable MCYna potential predicts some properties within experimental uncertainty. This highlights the importance of polarizability for the accurate prediction of the thermodynamic properties of water, particularly at temperatures beyond 298 K.
Intermolecular potentials and the accurate prediction of the thermodynamic properties of water
Shvab, I.; Sadus, Richard J.
2013-11-21
The ability of intermolecular potentials to correctly predict the thermodynamic properties of liquid water at a density of 0.998 g/cm{sup 3} for a wide range of temperatures (298–650 K) and pressures (0.1–700 MPa) is investigated. Molecular dynamics simulations are reported for the pressure, thermal pressure coefficient, thermal expansion coefficient, isothermal and adiabatic compressibilities, isobaric and isochoric heat capacities, and Joule-Thomson coefficient of liquid water using the non-polarizable SPC/E and TIP4P/2005 potentials. The results are compared with both experiment data and results obtained from the ab initio-based Matsuoka-Clementi-Yoshimine non-additive (MCYna) [J. Li, Z. Zhou, and R. J. Sadus, J. Chem. Phys. 127, 154509 (2007)] potential, which includes polarization contributions. The data clearly indicate that both the SPC/E and TIP4P/2005 potentials are only in qualitative agreement with experiment, whereas the polarizable MCYna potential predicts some properties within experimental uncertainty. This highlights the importance of polarizability for the accurate prediction of the thermodynamic properties of water, particularly at temperatures beyond 298 K.
Dolgounitcheva, O; Díaz-Tinoco, Manuel; Zakrzewski, V G; Richard, Ryan M; Marom, Noa; Sherrill, C David; Ortiz, J V
2016-02-01
Comparison of ab initio electron-propagator predictions of vertical ionization potentials and electron affinities of organic, acceptor molecules with benchmark calculations based on the basis set-extrapolated, coupled cluster single, double, and perturbative triple substitution method has enabled identification of self-energy approximations with mean, unsigned errors between 0.1 and 0.2 eV. Among the self-energy approximations that neglect off-diagonal elements in the canonical, Hartree-Fock orbital basis, the P3 method for electron affinities, and the P3+ method for ionization potentials provide the best combination of accuracy and computational efficiency. For approximations that consider the full self-energy matrix, the NR2 methods offer the best performance. The P3+ and NR2 methods successfully identify the correct symmetry label of the lowest cationic state in two cases, naphthalenedione and benzoquinone, where some other methods fail. PMID:26730459
Haley, William E.; Ibrahim, El-Sayed H.; Qu, Mingliang; Cernigliaro, Joseph G.; Goldfarb, David S.; McCollough, Cynthia H.
2015-01-01
Dual-energy computed tomography (DECT) has recently been suggested as the imaging modality of choice for kidney stones due to its ability to provide information on stone composition. Standard postprocessing of the dual-energy images accurately identifies uric acid stones, but not other types. Cystine stones can be identified from DECT images when analyzed with advanced postprocessing. This case report describes clinical implications of accurate diagnosis of cystine stones using DECT. PMID:26688770
Haley, William E; Ibrahim, El-Sayed H; Qu, Mingliang; Cernigliaro, Joseph G; Goldfarb, David S; McCollough, Cynthia H
2015-01-01
Dual-energy computed tomography (DECT) has recently been suggested as the imaging modality of choice for kidney stones due to its ability to provide information on stone composition. Standard postprocessing of the dual-energy images accurately identifies uric acid stones, but not other types. Cystine stones can be identified from DECT images when analyzed with advanced postprocessing. This case report describes clinical implications of accurate diagnosis of cystine stones using DECT. PMID:26688770
Theoretical studies of potential energy surfaces.
Harding, L. B.
2006-01-01
The goal of this program is to calculate accurate potential energy surfaces for both reactive and nonreactive systems. To do this the electronic Schroedinger equation must be solved. Our approach starts with multiconfiguration self-consistent field (MCSCF) reference wave functions. These reference wavefunctions are designed to be sufficiently flexible to accurately describe changes in electronic structure over a broad range of geometries. Dynamical electron correlation effects are included via multireference, singles and doubles configuration interaction (MRCI) calculations. With this approach, we are able to provide chemically useful predictions of the energetics for many systems. A second aspect of this program is the development of techniques to fit multi-dimensional potential surfaces to convenient, global, analytic functions that can then be used in dynamics calculations.
Accurate variational calculations and analysis of the HOCl vibrational energy spectrum
Skokov, S.; Qi, J.; Bowman, J.M.; Yang, C.; Gray, S.K.; Peterson, K.A. |; Mandelshtam, V.A.
1998-12-01
Large scale variational calculations for the vibrational states of HOCl are performed using a recently developed, accurate {ital ab initio} potential energy surface. Three different approaches for obtaining vibrational states are employed and contrasted; a truncation/recoupling scheme with direct diagonalization, the Lanczos method, and Chebyshev iteration with filter diagonalization. The complete spectrum of bound states for nonrotating HOCl is computed and analyzed within a random matrix theory framework. This analysis indicates almost entirely regular dynamics with only a small degree of chaos. The nearly regular spectral structure allows us to make assignments for the most significant part of the spectrum, based on analysis of coordinate expectation values and eigenfunctions. Ground state dipole moments and dipole transition probabilities are also calculated using accurate {ital ab initio} data. Computed values are in good agreement with available experimental data. Some exact rovibrational calculations for J=1, including Coriolis coupling, are performed. The exact results are nearly identical with those obtained from the adiabatic rotation approximation and very close to those from the centrifugal sudden approximation, thus indicating a very small degree of asymmetry and Coriolis coupling for the HOCl molecule. {copyright} {ital 1998 American Institute of Physics.}
Accurate Energies and Orbital Description in Semi-Local Kohn-Sham DFT
NASA Astrophysics Data System (ADS)
Lindmaa, Alexander; Kuemmel, Stephan; Armiento, Rickard
2015-03-01
We present our progress on a scheme in semi-local Kohn-Sham density-functional theory (KS-DFT) for improving the orbital description while still retaining the level of accuracy of the usual semi-local exchange-correlation (xc) functionals. DFT is a widely used tool for first-principles calculations of properties of materials. A given task normally requires a balance of accuracy and computational cost, which is well achieved with semi-local DFT. However, commonly used semi-local xc functionals have important shortcomings which often can be attributed to features of the corresponding xc potential. One shortcoming is an overly delocalized representation of localized orbitals. Recently a semi-local GGA-type xc functional was constructed to address these issues, however, it has the trade-off of lower accuracy of the total energy. We discuss the source of this error in terms of a surplus energy contribution in the functional that needs to be accounted for, and offer a remedy for this issue which formally stays within KS-DFT, and, which does not harshly increase the computational effort. The end result is a scheme that combines accurate total energies (e.g., relaxed geometries) with an improved orbital description (e.g., improved band structure).
NASA Astrophysics Data System (ADS)
Smith, R.; Flynn, C.; Candlish, G. N.; Fellhauer, M.; Gibson, B. K.
2015-04-01
We present accurate models of the gravitational potential produced by a radially exponential disc mass distribution. The models are produced by combining three separate Miyamoto-Nagai discs. Such models have been used previously to model the disc of the Milky Way, but here we extend this framework to allow its application to discs of any mass, scalelength, and a wide range of thickness from infinitely thin to near spherical (ellipticities from 0 to 0.9). The models have the advantage of simplicity of implementation, and we expect faster run speeds over a double exponential disc treatment. The potentials are fully analytical, and differentiable at all points. The mass distribution of our models deviates from the radial mass distribution of a pure exponential disc by <0.4 per cent out to 4 disc scalelengths, and <1.9 per cent out to 10 disc scalelengths. We tabulate fitting parameters which facilitate construction of exponential discs for any scalelength, and a wide range of disc thickness (a user-friendly, web-based interface is also available). Our recipe is well suited for numerical modelling of the tidal effects of a giant disc galaxy on star clusters or dwarf galaxies. We consider three worked examples; the Milky Way thin and thick disc, and a discy dwarf galaxy.
Holton, James M; Classen, Scott; Frankel, Kenneth A; Tainer, John A
2014-01-01
In macromolecular crystallography, the agreement between observed and predicted structure factors (Rcryst and Rfree) is seldom better than 20%. This is much larger than the estimate of experimental error (Rmerge). The difference between Rcryst and Rmerge is the R-factor gap. There is no such gap in small-molecule crystallography, for which calculated structure factors are generally considered more accurate than the experimental measurements. Perhaps the true noise level of macromolecular data is higher than expected? Or is the gap caused by inaccurate phases that trap refined models in local minima? By generating simulated diffraction patterns using the program MLFSOM, and including every conceivable source of experimental error, we show that neither is the case. Processing our simulated data yielded values that were indistinguishable from those of real data for all crystallographic statistics except the final Rcryst and Rfree. These values decreased to 3.8% and 5.5% for simulated data, suggesting that the reason for high R-factors in macromolecular crystallography is neither experimental error nor phase bias, but rather an underlying inadequacy in the models used to explain our observations. The present inability to accurately represent the entire macromolecule with both its flexibility and its protein-solvent interface may be improved by synergies between small-angle X-ray scattering, computational chemistry and crystallography. The exciting implication of our finding is that macromolecular data contain substantial hidden and untapped potential to resolve ambiguities in the true nature of the nanoscale, a task that the second century of crystallography promises to fulfill. Database Coordinates and structure factors for the real data have been submitted to the Protein Data Bank under accession 4tws. PMID:25040949
NASA Astrophysics Data System (ADS)
Trautt, Zachary T.; Tavazza, Francesca; Becker, Chandler A.
2015-10-01
The Materials Genome Initiative seeks to significantly decrease the cost and time of development and integration of new materials. Within the domain of atomistic simulations, several roadblocks stand in the way of reaching this goal. While the NIST Interatomic Potentials Repository hosts numerous interatomic potentials (force fields), researchers cannot immediately determine the best choice(s) for their use case. Researchers developing new potentials, specifically those in restricted environments, lack a comprehensive portfolio of efficient tools capable of calculating and archiving the properties of their potentials. This paper elucidates one solution to these problems, which uses Python-based scripts that are suitable for rapid property evaluation and human knowledge transfer. Calculation results are visible on the repository website, which reduces the time required to select an interatomic potential for a specific use case. Furthermore, property evaluation scripts are being integrated with modern platforms to improve discoverability and access of materials property data. To demonstrate these scripts and features, we will discuss the automation of stacking fault energy calculations and their application to additional elements. While the calculation methodology was developed previously, we are using it here as a case study in simulation automation and property calculations. We demonstrate how the use of Python scripts allows for rapid calculation in a more easily managed way where the calculations can be modified, and the results presented in user-friendly and concise ways. Additionally, the methods can be incorporated into other efforts, such as openKIM.
Matsui, Toru; Kitagawa, Yasutaka; Okumura, Mitsutaka; Shigeta, Yasuteru
2015-01-15
We computationally evaluated the standard hydrogen electrode (SHE) potential in aqueous phase and the Gibbs energy of a proton from the experimental pKa values of alcohol molecules. From the "golden standard" CCSD(T)/aug-cc-pVTZ level calculation, we estimated the SHE potential as 4.48 V, which is very close to the IUPAC-recommended experimental value of 4.44 V. As applications to the Gaussian-3 (G3) methods, which also reproduce the "golden standard" level calculations, we computed various pKa values and redox potentials for a vitamin series. For vitamin C, we support the experimental result of +0.35 V and predict the pKa value of d-ascorbic acid to be 3.7-3.9. Using a model molecule for nicotinamide adenine dinucleotide (NAD), we reproduced the redox potential and determined the order of the proton/electron addition, based on both the proton affinity and redox potential. PMID:25514626
Knight, Joseph W; Wang, Xiaopeng; Gallandi, Lukas; Dolgounitcheva, Olga; Ren, Xinguo; Ortiz, J Vincent; Rinke, Patrick; Körzdörfer, Thomas; Marom, Noa
2016-02-01
The performance of different GW methods is assessed for a set of 24 organic acceptors. Errors are evaluated with respect to coupled cluster singles, doubles, and perturbative triples [CCSD(T)] reference data for the vertical ionization potentials (IPs) and electron affinities (EAs), extrapolated to the complete basis set limit. Additional comparisons are made to experimental data, where available. We consider fully self-consistent GW (scGW), partial self-consistency in the Green's function (scGW0), non-self-consistent G0W0 based on several mean-field starting points, and a "beyond GW" second-order screened exchange (SOSEX) correction to G0W0. We also describe the implementation of the self-consistent Coulomb hole with screened exchange method (COHSEX), which serves as one of the mean-field starting points. The best performers overall are G0W0+SOSEX and G0W0 based on an IP-tuned long-range corrected hybrid functional with the former being more accurate for EAs and the latter for IPs. Both provide a balanced treatment of localized vs delocalized states and valence spectra in good agreement with photoemission spectroscopy (PES) experiments. PMID:26731609
Renewable energy potential in Colombia
NASA Astrophysics Data System (ADS)
Correa Guzman, Jose Luis
2008-12-01
Renewable energy flows are very large in comparison with humankind's use of energy. In principle, all our energy needs, both now and into the future, can be met by energy from renewable sources. After many years trying to develop the alternative energy potential of Colombia, a major effort is principally being made since 2000 to explore and assess the renewable resources of the entire country. Until 2000, the availability of conventional energy sources in Colombia prevented renewable energy exploration from reaching a higher level. However, the extreme energy crisis of 1992 - 1993 alerted the authorities and the community to the necessity for exploring alternative energy sources. This energy study is a general approach to the current and future renewable energy scenario of Colombia. It was prepared in response to the increased interest around the world and in particular in Colombia to develop its non-fossil energy prospective. It, therefore, represents a working document giving an initial impression of the possible scale of the main renewables sources as a response to the concern about energy security and fossil fuel dependence problems. The assumptions made and calculations reported may therefore be subject to revision as more information becomes available. The aim of this dissertation is not only to improve the public understanding and discussion of renewable energy matters in Colombia but also to stimulate the development and application of renewable energy, wherever they have prospects of economic viability and environmental acceptability. To achieve such goal this paper reviews several renewable technologies, their availability, contribution and feasibility in Colombia.
Sun, Y. Y.; Kim, Y. H.; Lee, K.; Zhang, S. B.
2008-01-01
Density functional theory (DFT) in the commonly used local density or generalized gradient approximation fails to describe van der Waals (vdW) interactions that are vital to organic, biological, and other molecular systems. Here, we propose a simple, efficient, yet accurate local atomic potential (LAP) approach, named DFT+LAP, for including vdW interactions in the framework of DFT. The LAPs for H, C, N, and O are generated by fitting the DFT+LAP potential energy curves of small molecule dimers to those obtained from coupled cluster calculations with single, double, and perturbatively treated triple excitations, CCSD(T). Excellent transferability of the LAPs is demonstrated by remarkable agreement with the JSCH-2005 benchmark database [P. Jurecka et al. Phys. Chem. Chem. Phys. 8, 1985 (2006)], which provides the interaction energies of CCSD(T) quality for 165 vdW and hydrogen-bonded complexes. For over 100 vdW dominant complexes in this database, our DFT+LAP calculations give a mean absolute deviation from the benchmark results less than 0.5 kcal/mol. The DFT+LAP approach involves no extra computational cost other than standard DFT calculations and no modification of existing DFT codes, which enables straightforward quantum simulations, such as ab initio molecular dynamics, on biomolecular systems, as well as on other organic systems.
Accurate exchange-correlation energies for the warm dense electron gas
NASA Astrophysics Data System (ADS)
Malone, Fionn; Blunt, Nicholas; Shepherd, James; Lee, Derek; Spencer, James; Foulkes, Matthew
The accurate treatment of matter at high temperatures and densities is of increasing importance to many fields in physics and chemistry, with applications ranging from planetary physics to inertial confinement fusion and plasmonic catalysis. Faithfully including the effects of temperature in density functional theory simulations of warm dense matter requires accurate results for the uniform electron gas (UEG) across the whole temperature-density plane. While accurate ground state quantum Monte Carlo data have existed for over 30 years, there remains significant disagreement between results obtained using different path integral Monte Carlo methods at finite temperature. To resolve this disagreement, we use the systematically improvable density matrix quantum Monte Carlo method to calculate the exchange-correlation energy of the UEG. We also demonstrate how the evaluation of free energies emerges naturally from our method.
An accurate dynamical electron diffraction algorithm for reflection high-energy electron diffraction
NASA Astrophysics Data System (ADS)
Huang, J.; Cai, C. Y.; Lv, C. L.; Zhou, G. W.; Wang, Y. G.
2015-12-01
The conventional multislice method (CMS) method, one of the most popular dynamical electron diffraction calculation procedures in transmission electron microscopy, was introduced to calculate reflection high-energy electron diffraction (RHEED) as it is well adapted to deal with the deviations from the periodicity in the direction parallel to the surface. However, in the present work, we show that the CMS method is no longer sufficiently accurate for simulating RHEED with the accelerating voltage 3-100 kV because of the high-energy approximation. An accurate multislice (AMS) method can be an alternative for more accurate RHEED calculations with reasonable computing time. A detailed comparison of the numerical calculation of the AMS method and the CMS method is carried out with respect to different accelerating voltages, surface structure models, Debye-Waller factors and glancing angles.
California Industrial Energy Efficiency Potential
Coito, Fred; Worrell, Ernst; Price, Lynn; Masanet, Eric; RafaelFriedmann; Rufo, Mike
2005-06-01
This paper presents an overview of the modeling approach andhighlights key findings of a California industrial energy efficiencypotential study. In addition to providing estimates of technical andeconomic potential, the study examines achievable program potential undervarious program-funding scenarios. The focus is on electricity andnatural gas savings for manufacturing in the service territories ofCalifornia's investor-owned utilities (IOUs). The assessment is conductedby industry type and by end use. Both crosscutting technologies andindustry-specific process measures are examined. Measure penetration intothe marketplace is modeled as a function of customer awareness, measurecost effectiveness, and perceived market barriers. Data for the studycomes from a variety of sources, including: utility billing records, theEnergy Information Association (EIA) Manufacturing Energy ConsumptionSurvey (MECS), state-sponsored avoided cost studies, energy efficiencyprogram filings, and technology savings and cost data developed throughLawrence Berkeley National Laboratory (LBNL). The study identifies 1,706GWh and 47 Mth (million therms) per year of achievable potential over thenext twelve years under recent levels of program expenditures, accountingfor 5.2 percent of industrial electricity consumption and 1.3 percent ofindustrial natural gas consumption. These estimates grow to 2,748 GWh and192 Mth per year if all cost-effective and achievable opportunities arepursued. Key industrial electricity end uses, in terms of energy savingspotential, include compressed air and pumping systems that combine toaccount for about half of the total achievable potential estimates. Fornatural gas, savings are concentrated in the boiler and process heatingend uses, accounting for over 99 percent to total achievablepotential.
An accurate nucleon-nucleon potential with charge-independence breaking
Wiringa, R.B.; Stoks, V.G.J.; Schiavilla, R.
1995-08-01
We constructed a new NN potential, designated Argonne v{sub 18}, with explicit charge-independence breaking. It supersedes our older v{sub 14} model, which was our standard nonrelativistic NN potential for most of the last decade. The main part of the new potential is charge-independent, like the old v{sub 14} model, with 14 components, each consisting of a radial function v{sub p}(r{sub 12}) multiplied by an operator: 1, {sigma}{sub 1}{center_dot}{sigma}{sub 2}, S{sub 12}, L{center_dot}S, L{sup 2}, L{sup 2}{sigma}{sub 1}{center_dot}{sigma}{sub 2}, and (L{center_dot}S){sup 2}, and each of these times {tau}{sub l}{center_dot}{tau}{sub 2}. Three charge-dependent and one charge-asymmetric operators are added along with a complete electromagnetic interaction, resulting in a model that fits pp, np, and nn data simultaneously. The charge-dependent operators are obtained by multiplying the spin operators 1, {sigma}{sub 1}{center_dot}{sigma}{sub 2}, and S{sub 12} by the isotensor T{sub 12} = 3{tau}{sub 1z}{tau}{sub 2z} - {tau}{sub 1}{center_dot}{tau}{sub 2}, which differentiates between np and pp or nn T = 1 states. A major source of charge dependence in NN interactions is the mass difference of the charged and neutral pions, which is carefully treated in the new model. The charge-asymmetric operator is {tau}{sub 1z}+{tau}{sub 2z} which splits pp and nn states; it is constrained by the difference between nn and pp scattering lengths. The electromagnetic interaction includes Coulomb, Darwin-Foldy, vacuum polarization, and magnetic moment terms. The potential was fit directly to the Nijmegen pp and np scattering database as well as the nn scattering length and deuteron binding energy. With {approximately}40 adjustable parameters it gives an excellent {chi}{sup 2}/degree of freedom of 1.09 for 4301 pp and np data in the range 0-350 MeV. A consistent set of two-body charge and current operators has also been derived to evaluate the deuteron electromagnetic form factors.
Computed potential energy surfaces for chemical reactions
NASA Technical Reports Server (NTRS)
Walch, Stephen P.; Levin, Eugene
1993-01-01
A new global potential energy surface (PES) is being generated for O(P-3) + H2 yields OH + H. This surface is being fit using the rotated Morse oscillator method, which was used to fit the previous POL-CI surface. The new surface is expected to be more accurate and also includes a much more complete sampling of bent geometries. A new study has been undertaken of the reaction N + O2 yields NO + O. The new studies have focused on the region of the surface near a possible minimum corresponding to the peroxy form of NOO. A large portion of the PES for this second reaction has been mapped out. Since state to state cross sections for the reaction are important in the chemistry of high temperature air, these studies will probably be extended to permit generation of a new global potential for reaction.
The potential of renewable energy
Not Available
1990-03-01
On June 27 and 28, 1989, the US Department of Energy (DOE) national laboratories were convened to discuss plans for the development of a National Energy Strategy (NES) and, in particular, the analytic needs in support of NES that could be addressed by the laboratories. As a result of that meeting, interlaboratory teams were formed to produce analytic white papers on key topics, and a lead laboratory was designated for each core laboratory team. The broad-ranging renewables assignment is summarized by the following issue statement from the Office of Policy, Planning and Analysis: to what extent can renewable energy technologies contribute to diversifying sources of energy supply What are the major barriers to greater renewable energy use and what is the potential timing of widespread commercialization for various categories of applications This report presents the results of the intensive activity initiated by the June 1989 meeting to produce a white paper on renewable energy. Scores of scientists, analysts, and engineers in the five core laboratories gave generously of their time over the past eight months to produce this document. Their generous, constructive efforts are hereby gratefully acknowledged. 126 refs., 44 figs., 32 tabs.
A Cavity Corrected 3D-RISM Functional for Accurate Solvation Free Energies
2014-01-01
We show that an Ng bridge function modified version of the three-dimensional reference interaction site model (3D-RISM-NgB) solvation free energy method can accurately predict the hydration free energy (HFE) of a set of 504 organic molecules. To achieve this, a single unique constant parameter was adjusted to the computed HFE of single atom Lennard-Jones solutes. It is shown that 3D-RISM is relatively accurate at predicting the electrostatic component of the HFE without correction but requires a modification of the nonpolar contribution that originates in the formation of the cavity created by the solute in water. We use a free energy functional with the Ng scaling of the direct correlation function [Ng, K. C. J. Chem. Phys.1974, 61, 2680]. This produces a rapid, reliable small molecule HFE calculation for applications in drug design. PMID:24634616
A Cavity Corrected 3D-RISM Functional for Accurate Solvation Free Energies.
Truchon, Jean-François; Pettitt, B Montgomery; Labute, Paul
2014-03-11
We show that an Ng bridge function modified version of the three-dimensional reference interaction site model (3D-RISM-NgB) solvation free energy method can accurately predict the hydration free energy (HFE) of a set of 504 organic molecules. To achieve this, a single unique constant parameter was adjusted to the computed HFE of single atom Lennard-Jones solutes. It is shown that 3D-RISM is relatively accurate at predicting the electrostatic component of the HFE without correction but requires a modification of the nonpolar contribution that originates in the formation of the cavity created by the solute in water. We use a free energy functional with the Ng scaling of the direct correlation function [Ng, K. C. J. Chem. Phys. 1974, 61, 2680]. This produces a rapid, reliable small molecule HFE calculation for applications in drug design. PMID:24634616
NASA Astrophysics Data System (ADS)
Moore, Keith; McLaughlin, Brendan M.; Lane, Ian C.
2016-04-01
BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio 2Σ+ potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron effective core-potential and even-tempered augmented polarized core-valence basis sets (aug-pCVnZ-PP, n = 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy De for the X2Σ+ state (extrapolated to the CBS limit) is 16 895.12 cm-1 (2.094 eV), which agrees within 0.1% of a revised experimental value of <16 910.6 cm-1, while the calculated re is within 0.03 pm of the experimental result.
Molecular Multipole Potential Energy Functions for Water.
Tan, Ming-Liang; Tran, Kelly N; Pickard, Frank C; Simmonett, Andrew C; Brooks, Bernard R; Ichiye, Toshiko
2016-03-01
Water is the most common liquid on this planet, with many unique properties that make it essential for life as we know it. These properties must arise from features in the charge distribution of a water molecule, so it is essential to capture these features in potential energy functions for water to reproduce its liquid state properties in computer simulations. Recently, models that utilize a multipole expansion located on a single site in the water molecule, or "molecular multipole models", have been shown to rival and even surpass site models with up to five sites in reproducing both the electrostatic potential around a molecule and a variety of liquid state properties in simulations. However, despite decades of work using multipoles, confusion still remains about how to truncate the multipole expansions efficiently and accurately. This is particularly important when using molecular multipole expansions to describe water molecules in the liquid state, where the short-range interactions must be accurate, because the higher order multipoles of a water molecule are large. Here, truncation schemes designed for a recent efficient algorithm for multipoles in molecular dynamics simulations are assessed for how well they reproduce results for a simple three-site model of water when the multipole moments and Lennard-Jones parameters of that model are used. In addition, the multipole analysis indicates that site models that do not account for out-of-plane electron density overestimate the stability of a non-hydrogen-bonded conformation, leading to serious consequences for the simulated liquid. PMID:26562223
Oyeyemi, Victor B.; Krisiloff, David B.; Keith, John A.; Libisch, Florian; Pavone, Michele; Carter, Emily A.
2014-01-28
Oxygenated hydrocarbons play important roles in combustion science as renewable fuels and additives, but many details about their combustion chemistry remain poorly understood. Although many methods exist for computing accurate electronic energies of molecules at equilibrium geometries, a consistent description of entire combustion reaction potential energy surfaces (PESs) requires multireference correlated wavefunction theories. Here we use bond dissociation energies (BDEs) as a foundational metric to benchmark methods based on multireference configuration interaction (MRCI) for several classes of oxygenated compounds (alcohols, aldehydes, carboxylic acids, and methyl esters). We compare results from multireference singles and doubles configuration interaction to those utilizing a posteriori and a priori size-extensivity corrections, benchmarked against experiment and coupled cluster theory. We demonstrate that size-extensivity corrections are necessary for chemically accurate BDE predictions even in relatively small molecules and furnish examples of unphysical BDE predictions resulting from using too-small orbital active spaces. We also outline the specific challenges in using MRCI methods for carbonyl-containing compounds. The resulting complete basis set extrapolated, size-extensivity-corrected MRCI scheme produces BDEs generally accurate to within 1 kcal/mol, laying the foundation for this scheme's use on larger molecules and for more complex regions of combustion PESs.
Oyeyemi, Victor B; Krisiloff, David B; Keith, John A; Libisch, Florian; Pavone, Michele; Carter, Emily A
2014-01-28
Oxygenated hydrocarbons play important roles in combustion science as renewable fuels and additives, but many details about their combustion chemistry remain poorly understood. Although many methods exist for computing accurate electronic energies of molecules at equilibrium geometries, a consistent description of entire combustion reaction potential energy surfaces (PESs) requires multireference correlated wavefunction theories. Here we use bond dissociation energies (BDEs) as a foundational metric to benchmark methods based on multireference configuration interaction (MRCI) for several classes of oxygenated compounds (alcohols, aldehydes, carboxylic acids, and methyl esters). We compare results from multireference singles and doubles configuration interaction to those utilizing a posteriori and a priori size-extensivity corrections, benchmarked against experiment and coupled cluster theory. We demonstrate that size-extensivity corrections are necessary for chemically accurate BDE predictions even in relatively small molecules and furnish examples of unphysical BDE predictions resulting from using too-small orbital active spaces. We also outline the specific challenges in using MRCI methods for carbonyl-containing compounds. The resulting complete basis set extrapolated, size-extensivity-corrected MRCI scheme produces BDEs generally accurate to within 1 kcal/mol, laying the foundation for this scheme's use on larger molecules and for more complex regions of combustion PESs. PMID:25669533
NASA Astrophysics Data System (ADS)
Oyeyemi, Victor B.; Krisiloff, David B.; Keith, John A.; Libisch, Florian; Pavone, Michele; Carter, Emily A.
2014-01-01
Oxygenated hydrocarbons play important roles in combustion science as renewable fuels and additives, but many details about their combustion chemistry remain poorly understood. Although many methods exist for computing accurate electronic energies of molecules at equilibrium geometries, a consistent description of entire combustion reaction potential energy surfaces (PESs) requires multireference correlated wavefunction theories. Here we use bond dissociation energies (BDEs) as a foundational metric to benchmark methods based on multireference configuration interaction (MRCI) for several classes of oxygenated compounds (alcohols, aldehydes, carboxylic acids, and methyl esters). We compare results from multireference singles and doubles configuration interaction to those utilizing a posteriori and a priori size-extensivity corrections, benchmarked against experiment and coupled cluster theory. We demonstrate that size-extensivity corrections are necessary for chemically accurate BDE predictions even in relatively small molecules and furnish examples of unphysical BDE predictions resulting from using too-small orbital active spaces. We also outline the specific challenges in using MRCI methods for carbonyl-containing compounds. The resulting complete basis set extrapolated, size-extensivity-corrected MRCI scheme produces BDEs generally accurate to within 1 kcal/mol, laying the foundation for this scheme's use on larger molecules and for more complex regions of combustion PESs.
Toward Accurate Modeling of the Effect of Ion-Pair Formation on Solute Redox Potential.
Qu, Xiaohui; Persson, Kristin A
2016-09-13
A scheme to model the dependence of a solute redox potential on the supporting electrolyte is proposed, and the results are compared to experimental observations and other reported theoretical models. An improved agreement with experiment is exhibited if the effect of the supporting electrolyte on the redox potential is modeled through a concentration change induced via ion pair formation with the salt, rather than by only considering the direct impact on the redox potential of the solute itself. To exemplify the approach, the scheme is applied to the concentration-dependent redox potential of select molecules proposed for nonaqueous flow batteries. However, the methodology is general and enables rational computational electrolyte design through tuning of the operating window of electrochemical systems by shifting the redox potential of its solutes; including potentially both salts as well as redox active molecules. PMID:27500744
Convective Available Potential Energy of World Ocean
NASA Astrophysics Data System (ADS)
Su, Z.; Ingersoll, A. P.; Thompson, A. F.
2012-12-01
Here, for the first time, we propose the concept of Ocean Convective Available Potential Energy (OCAPE), which is the maximum kinetic energy (KE) per unit seawater mass achievable by ocean convection. OCAPE occurs through a different mechanism from atmospheric CAPE, and involves the interplay of temperature and salinity on the equation of state of seawater. The thermobaric effect, which arises because the thermal coefficient of expansion increases with depth, is an important ingredient of OCAPE. We develop an accurate algorithm to calculate the OCAPE for a given temperature and salinity profile. We then validate our calculation of OCAPE by comparing it with the conversion of OCAPE to KE in a 2-D numerical model. We propose that OCAPE is an important energy source of ocean deep convection and contributes to deep water formation. OCAPE, like Atmospheric CAPE, can help predict deep convection and may also provide a useful constraint for modelling deep convection in ocean GCMs. We plot the global distribution of OCAPE using data from the World Ocean Atlas 2009 (WOA09) and see many important features. These include large values of OCAPE in the Labrador, Greenland, Weddell and Mediterranean Seas, which are consistent with our present observations and understanding, but also identify some new features like the OCAPE pattern in the Antarctic Circumpolar Current (ACC). We propose that the diagnosis of OCAPE can improve our understanding of global patterns of ocean convection and deep water formation as well as ocean stratification, the meridional overturning circulation and mixed layer processes. The background of this work is briefly introduced as below. Open-ocean deep convection can significantly modify water properties both at the ocean surface and throughout the water column (Gordon 1982). Open-ocean convection is also an important mechanism for Ocean Deep Water formation and the transport of heat, freshwater and nutrient (Marshall and Schott 1999). Open
Puzzarini, Cristina; Ali, Ashraf; Biczysko, Malgorzata; Barone, Vincenzo
2014-09-10
An accurate spectroscopic characterization of protonated oxirane has been carried out by means of state-of-the-art computational methods and approaches. The calculated spectroscopic parameters from our recent computational investigation of oxirane together with the corresponding experimental data available were used to assess the accuracy of our predicted rotational and IR spectra of protonated oxirane. We found an accuracy of about 10 cm{sup –1} for vibrational transitions (fundamentals as well as overtones and combination bands) and, in relative terms, of 0.1% for rotational transitions. We are therefore confident that the spectroscopic data provided herein are a valuable support for the detection of protonated oxirane not only in Titan's atmosphere but also in the interstellar medium.
Puzzarini, Cristina; Ali, Ashraf; Biczysko, Malgorzata; Barone, Vincenzo
2015-01-01
An accurate spectroscopic characterization of protonated oxirane has been carried out by means of state-of-the-art computational methods and approaches. The calculated spectroscopic parameters from our recent computational investigation of oxirane together with the corresponding experimental data available were used to assess the accuracy of our predicted rotational and IR spectra of protonated oxirane. We found an accuracy of about 10 cm−1 for vibrational transitions (fundamentals as well as overtones and combination bands) and, in relative terms, of 0.1% for rotational transitions. We are therefore confident that the spectroscopic data provided herein are a valuable support for the detection of protonated oxirane not only in Titan’s atmosphere but also in the interstellar medium. PMID:26543241
Apparatus for use in rapid and accurate controlled-potential coulometric analysis
Frazzini, Thomas L.; Holland, Michael K.; Pietri, Charles E.; Weiss, Jon R.
1981-01-01
An apparatus for controlled-potential coulometric analysis of a solution includes a cell to contain the solution to be analyzed and a plurality of electrodes to contact the solution in the cell. Means are provided to stir the solution and to control the atmosphere above it. A potentiostat connected to the electrodes controls potential differences among the electrodes. An electronic circuit connected to the potentiostat provides analog-to-digital conversion and displays a precise count of charge transfer during a desired chemical process. This count provides a measure of the amount of an unknown substance in the solution.
Finite domain simulations with adaptive boundaries: accurate potentials and nonequilibrium movesets.
Wagoner, Jason A; Pande, Vijay S
2013-12-21
We extend the theory of hybrid explicit/implicit solvent models to include an explicit domain that grows and shrinks in response to a solute's evolving configuration. The goal of this model is to provide an appropriate but not excessive amount of solvent detail, and the inclusion of an adjustable boundary provides a significant computational advantage for solutes that explore a range of configurations. In addition to the theoretical development, a successful implementation of this method requires (1) an efficient moveset that propagates the boundary as a new coordinate of the system, and (2) an accurate continuum solvent model with parameters that are transferable to an explicit domain of any size. We address these challenges and develop boundary updates using Monte Carlo moves biased by nonequilibrium paths. We obtain the desired level of accuracy using a "decoupling interface" that we have previously shown to remove boundary artifacts common to hybrid solvent models. Using an uncharged, coarse-grained solvent model, we then study the efficiency of nonequilibrium paths that a simulation takes by quantifying the dissipation. In the spirit of optimization, we study this quantity over a range of simulation parameters. PMID:24359359
Schütt, Heiko H; Harmeling, Stefan; Macke, Jakob H; Wichmann, Felix A
2016-05-01
The psychometric function describes how an experimental variable, such as stimulus strength, influences the behaviour of an observer. Estimation of psychometric functions from experimental data plays a central role in fields such as psychophysics, experimental psychology and in the behavioural neurosciences. Experimental data may exhibit substantial overdispersion, which may result from non-stationarity in the behaviour of observers. Here we extend the standard binomial model which is typically used for psychometric function estimation to a beta-binomial model. We show that the use of the beta-binomial model makes it possible to determine accurate credible intervals even in data which exhibit substantial overdispersion. This goes beyond classical measures for overdispersion-goodness-of-fit-which can detect overdispersion but provide no method to do correct inference for overdispersed data. We use Bayesian inference methods for estimating the posterior distribution of the parameters of the psychometric function. Unlike previous Bayesian psychometric inference methods our software implementation-psignifit 4-performs numerical integration of the posterior within automatically determined bounds. This avoids the use of Markov chain Monte Carlo (MCMC) methods typically requiring expert knowledge. Extensive numerical tests show the validity of the approach and we discuss implications of overdispersion for experimental design. A comprehensive MATLAB toolbox implementing the method is freely available; a python implementation providing the basic capabilities is also available. PMID:27013261
Sansone, Giuseppe; Maschio, Lorenzo; Usvyat, Denis; Schütz, Martin; Karttunen, Antti
2016-01-01
The black phosphorus (black-P) crystal is formed of covalently bound layers of phosphorene stacked together by weak van der Waals interactions. An experimental measurement of the exfoliation energy of black-P is not available presently, making theoretical studies the most important source of information for the optimization of phosphorene production. Here, we provide an accurate estimate of the exfoliation energy of black-P on the basis of multilevel quantum chemical calculations, which include the periodic local Møller-Plesset perturbation theory of second order, augmented by higher-order corrections, which are evaluated with finite clusters mimicking the crystal. Very similar results are also obtained by density functional theory with the D3-version of Grimme's empirical dispersion correction. Our estimate of the exfoliation energy for black-P of -151 meV/atom is substantially larger than that of graphite, suggesting the need for different strategies to generate isolated layers for these two systems. PMID:26651397
Accurate nonrelativistic ground-state energies of 3d transition metal atoms
Scemama, A.; Applencourt, T.; Giner, E.; Caffarel, M.
2014-12-28
We present accurate nonrelativistic ground-state energies of the transition metal atoms of the 3d series calculated with Fixed-Node Diffusion Monte Carlo (FN-DMC). Selected multi-determinantal expansions obtained with the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) method and including the most prominent determinants of the full configuration interaction expansion are used as trial wavefunctions. Using a maximum of a few tens of thousands determinants, fixed-node errors on total DMC energies are found to be greatly reduced for some atoms with respect to those obtained with Hartree-Fock nodes. To the best of our knowledge, the FN-DMC/(CIPSI nodes) ground-state energies presented here are the lowest variational total energies reported so far. They differ from the recently recommended non-variational values of McCarthy and Thakkar [J. Chem. Phys. 136, 054107 (2012)] only by a few percents of the correlation energy. Thanks to the variational property of FN-DMC total energies, our results provide exact lower bounds for the absolute value of all-electron correlation energies, |E{sub c}|.
Do Bond Functions Help for the Calculation of Accurate Bond Energies?
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Arnold, James (Technical Monitor)
1998-01-01
The bond energies of 8 chemically bound diatomics are computed using several basis sets with and without bond functions (BF). The bond energies obtained using the aug-pVnZ+BF basis sets (with a correction for basis set superposition error, BSSE) tend to be slightly smaller that the results obtained using the aug-pV(n+I)Z basis sets, but slightly larger than the BSSE corrected aug-pV(n+I)Z results. The aug-cc-pVDZ+BF and aug-cc-pVTZ+BF basis sets yield reasonable estimates of bond energies, but, in most cases, these results cannot be considered highly accurate. Extrapolation of the results obtained with basis sets including bond functions appears to be inferior to the results obtained by extrapolation using atom-centered basis sets. Therefore bond functions do not appear to offer a path for obtaining highly accurate results for chemically bound systems at a lower computational cost than atom centered basis sets.
Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations
NASA Technical Reports Server (NTRS)
Anderson, W. Kyle; Batina, John T.
1988-01-01
Parameter studies are conducted using the Euler and potential flow equation models for steady and unsteady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux vector splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flow results are made, as well as with experimental data where available.
Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations
NASA Technical Reports Server (NTRS)
Anderson, W. Kyle; Batina, John T.
1988-01-01
Parameter studies are conducted using the Euler and potential flow equation models for unsteady and steady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux-vector-splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flows results are made as well as with experimental data where available.
Highly accurate analytical energy of a two-dimensional exciton in a constant magnetic field
NASA Astrophysics Data System (ADS)
Hoang, Ngoc-Tram D.; Nguyen, Duy-Anh P.; Hoang, Van-Hung; Le, Van-Hoang
2016-08-01
Explicit expressions are given for analytically describing the dependence of the energy of a two-dimensional exciton on magnetic field intensity. These expressions are highly accurate with the precision of up to three decimal places for the whole range of the magnetic field intensity. The results are shown for the ground state and some excited states; moreover, we have all formulae to obtain similar expressions of any excited state. Analysis of numerical results shows that the precision of three decimal places is maintained for the excited states with the principal quantum number of up to n=100.
Accurate thermochemistry for larger molecules : gaussian-2 theory with bond separation energies.
Raghavachari, K.; Stefanov, B. B.; Curtiss, L. A.; Lucent Tech.
1997-04-22
Gaussian-2 (G2) theory is combined with isodesmic bond separation reaction energies to yield accurate thermochemistry for larger molecules. For a test set of 40 molecules composed of H, C, O, and N, our method yields enthalpies of formation, {Delta}H{sub f}{sup 0}(298 K), with a mean absolute deviation from experiment of only 0.5 kcal/mol. This is an improvement of a factor of three over the deviation of 1.5 kcal/mol seen in standard G2 theory.
Amacher, David E
2016-07-01
Context The silencing or activation of cancer-associated genes by epigenetic mechanisms can ultimately lead to the clonal expansion of cancer cells. Objective The aim of this review is to summarize all relevant epigenetic biomarkers that have been proposed to date for the diagnosis of some prevalent human cancers. Methods A Medline search for the terms epigenetic biomarkers, human cancers, DNA methylation, histone modifications and microRNAs was performed. Results One hundred fifty-seven relevant publications were found and reviewed. Conclusion To date, a significant number of potential epigenetic cancer biomarkers of human cancer have been investigated, and some have advanced to clinical implementation. PMID:26983778
Exact kinetic energy enables accurate evaluation of weak interactions by the FDE-vdW method
Sinha, Debalina; Pavanello, Michele
2015-08-28
The correlation energy of interaction is an elusive and sought-after interaction between molecular systems. By partitioning the response function of the system into subsystem contributions, the Frozen Density Embedding (FDE)-vdW method provides a computationally amenable nonlocal correlation functional based on the adiabatic connection fluctuation dissipation theorem applied to subsystem density functional theory. In reproducing potential energy surfaces of weakly interacting dimers, we show that FDE-vdW, either employing semilocal or exact nonadditive kinetic energy functionals, is in quantitative agreement with high-accuracy coupled cluster calculations (overall mean unsigned error of 0.5 kcal/mol). When employing the exact kinetic energy (which we term the Kohn-Sham (KS)-vdW method), the binding energies are generally closer to the benchmark, and the energy surfaces are also smoother.
Santra, Robin; Christ, Kevin V.; Greene, Chris H.
2004-04-01
The first three electronically excited states in the alkaline-earth-metal atoms magnesium, calcium, and strontium comprise the (nsnp){sup 3}P{sub J}{sup o}(J=0,1,2) fine-structure manifold. All three states are metastable and are of interest for optical atomic clocks as well as for cold-collision physics. An efficient technique--based on a physically motivated potential that models the presence of the ionic core--is employed to solve the Schroedinger equation for the two-electron valence shell. In this way, radiative lifetimes, laser-induced clock shifts, and long-range interaction parameters are calculated for metastable Mg, Ca, and Sr.
Fang, Tao; Li, Wei; Gu, Fangwei; Li, Shuhua
2015-01-13
We extend the generalized energy-based fragmentation (GEBF) approach to molecular crystals under periodic boundary conditions (PBC), and we demonstrate the performance of the method for a variety of molecular crystals. With this approach, the lattice energy of a molecular crystal can be obtained from the energies of a series of embedded subsystems, which can be computed with existing advanced molecular quantum chemistry methods. The use of the field compensation method allows the method to take long-range electrostatic interaction of the infinite crystal environment into account and make the method almost translationally invariant. The computational cost of the present method scales linearly with the number of molecules in the unit cell. Illustrative applications demonstrate that the PBC-GEBF method with explicitly correlated quantum chemistry methods is capable of providing accurate descriptions on the lattice energies and structures for various types of molecular crystals. In addition, this approach can be employed to quantify the contributions of various intermolecular interactions to the theoretical lattice energy. Such qualitative understanding is very useful for rational design of molecular crystals. PMID:26574207
Zeller, Rudolf
2013-03-13
Although the full-potential Korringa-Kohn-Rostoker Green function method yields accurate results for many physical properties, the convergence of calculated total energies with respect to the angular momentum cutoff is usually considered to be less satisfactory. This is surprising because accurate single-particle energies are expected if they are calculated by Lloyd's formula and because accurate densities and hence accurate double-counting energies should result from the total energy variational principle. It is shown how the concept of projection potentials can be used as a tool to analyse the convergence behaviour. The key factor blocking fast convergence is identified and it is illustrated how total energies can be improved with only a modest increase of computing time. PMID:23396831
Geothermal Energy: Tapping the Potential
ERIC Educational Resources Information Center
Johnson, Bill
2008-01-01
Ground source geothermal energy enables one to tap into the earth's stored renewable energy for heating and cooling facilities. Proper application of ground-source geothermal technology can have a dramatic impact on the efficiency and financial performance of building energy utilization (30%+). At the same time, using this alternative energy…
Fast and accurate modeling of molecular atomization energies with machine learning.
Rupp, Matthias; Tkatchenko, Alexandre; Müller, Klaus-Robert; von Lilienfeld, O Anatole
2012-02-01
We introduce a machine learning model to predict atomization energies of a diverse set of organic molecules, based on nuclear charges and atomic positions only. The problem of solving the molecular Schrödinger equation is mapped onto a nonlinear statistical regression problem of reduced complexity. Regression models are trained on and compared to atomization energies computed with hybrid density-functional theory. Cross validation over more than seven thousand organic molecules yields a mean absolute error of ∼10 kcal/mol. Applicability is demonstrated for the prediction of molecular atomization potential energy curves. PMID:22400967
Accurate Complete Basis Set Extrapolation of Direct Random Phase Correlation Energies.
Mezei, Pál D; Csonka, Gábor I; Ruzsinszky, Adrienn
2015-08-11
The direct random phase approximation (dRPA) is a promising way to obtain improvements upon the standard semilocal density functional results in many aspects of computational chemistry. In this paper, we address the slow convergence of the calculated dRPA correlation energy with the increase of the quality and size of the popular Gaussian-type Dunning's correlation consistent aug-cc-pVXZ split valence atomic basis set family. The cardinal number X controls the size of the basis set, and we use X = 3-6 in this study. It is known that even the very expensive X = 6 basis sets lead to large errors for the dRPA correlation energy, and thus complete basis set extrapolation is necessary. We study the basis set convergence of the dRPA correlation energies on a set of 65 hydrocarbon isomers from CH4 to C6H6. We calculate the iterative density fitted dRPA correlation energies using an efficient algorithm based on the CC-like form of the equations using the self-consistent HF orbitals. We test the popular inverse cubic, the optimized exponential, and inverse power formulas for complete basis set extrapolation. We have found that the optimized inverse power based extrapolation delivers the best energies. Further analysis showed that the optimal exponent depends on the molecular structure, and the most efficient two-point energy extrapolations that use X = 3 and 4 can be improved considerably by considering the atomic composition and hybridization states of the atoms in the molecules. Our results also show that the optimized exponents that yield accurate X = 3 and 4 extrapolated dRPA energies for atoms or small molecules might be inaccurate for larger molecules. PMID:26574475
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.
Accurate calculation of binding energies for molecular clusters - Assessment of different models
NASA Astrophysics Data System (ADS)
Friedrich, Joachim; Fiedler, Benjamin
2016-06-01
In this work we test different strategies to compute high-level benchmark energies for medium-sized molecular clusters. We use the incremental scheme to obtain CCSD(T)/CBS energies for our test set and carefully validate the accuracy for binding energies by statistical measures. The local errors of the incremental scheme are <1 kJ/mol. Since they are smaller than the basis set errors, we obtain higher total accuracy due to the applicability of larger basis sets. The final CCSD(T)/CBS benchmark values are ΔE = - 278.01 kJ/mol for (H2O)10, ΔE = - 221.64 kJ/mol for (HF)10, ΔE = - 45.63 kJ/mol for (CH4)10, ΔE = - 19.52 kJ/mol for (H2)20 and ΔE = - 7.38 kJ/mol for (H2)10 . Furthermore we test state-of-the-art wave-function-based and DFT methods. Our benchmark data will be very useful for critical validations of new methods. We find focal-point-methods for estimating CCSD(T)/CBS energies to be highly accurate and efficient. For foQ-i3CCSD(T)-MP2/TZ we get a mean error of 0.34 kJ/mol and a standard deviation of 0.39 kJ/mol.
Lee, M.W.; Meuwly, M.
2013-01-01
The evaluation of hydration free energies is a sensitive test to assess force fields used in atomistic simulations. We showed recently that the vibrational relaxation times, 1D- and 2D-infrared spectroscopies for CN(-) in water can be quantitatively described from molecular dynamics (MD) simulations with multipolar force fields and slightly enlarged van der Waals radii for the C- and N-atoms. To validate such an approach, the present work investigates the solvation free energy of cyanide in water using MD simulations with accurate multipolar electrostatics. It is found that larger van der Waals radii are indeed necessary to obtain results close to the experimental values when a multipolar force field is used. For CN(-), the van der Waals ranges refined in our previous work yield hydration free energy between -72.0 and -77.2 kcal mol(-1), which is in excellent agreement with the experimental data. In addition to the cyanide ion, we also study the hydroxide ion to show that the method used here is readily applicable to similar systems. Hydration free energies are found to sensitively depend on the intermolecular interactions, while bonded interactions are less important, as expected. We also investigate in the present work the possibility of applying the multipolar force field in scoring trajectories generated using computationally inexpensive methods, which should be useful in broader parametrization studies with reduced computational resources, as scoring is much faster than the generation of the trajectories.
Lee, Myung Won; Meuwly, Markus
2013-12-14
The evaluation of hydration free energies is a sensitive test to assess force fields used in atomistic simulations. We showed recently that the vibrational relaxation times, 1D- and 2D-infrared spectroscopies for CN(-) in water can be quantitatively described from molecular dynamics (MD) simulations with multipolar force fields and slightly enlarged van der Waals radii for the C- and N-atoms. To validate such an approach, the present work investigates the solvation free energy of cyanide in water using MD simulations with accurate multipolar electrostatics. It is found that larger van der Waals radii are indeed necessary to obtain results close to the experimental values when a multipolar force field is used. For CN(-), the van der Waals ranges refined in our previous work yield hydration free energy between -72.0 and -77.2 kcal mol(-1), which is in excellent agreement with the experimental data. In addition to the cyanide ion, we also study the hydroxide ion to show that the method used here is readily applicable to similar systems. Hydration free energies are found to sensitively depend on the intermolecular interactions, while bonded interactions are less important, as expected. We also investigate in the present work the possibility of applying the multipolar force field in scoring trajectories generated using computationally inexpensive methods, which should be useful in broader parametrization studies with reduced computational resources, as scoring is much faster than the generation of the trajectories. PMID:24170171
Potential of renewable and alternative energy sources
NASA Astrophysics Data System (ADS)
Konovalov, V.; Pogharnitskaya, O.; Rostovshchikova, A.; Matveenko, I.
2015-11-01
The article deals with application potential of clean alternative renewable energy sources. By means of system analysis the forecast for consumption of electrical energy in Tomsk Oblast as well as main energy sources of existing energy system have been studied up to 2018. Engineering potential of renewable and alternative energy sources is evaluated. Besides, ranking in the order of their efficiency descending is performed. It is concluded that Tomsk Oblast has high potential of alternative and renewable energy sources, among which the most promising development perspective is implementation of gasification stations to save fuel consumed by diesel power stations as well as building wind-power plants.
Potential energy of a dinuclear system
Adamyan, G.G.; Antonenko, N.V.; Jolos, R.V.; Ivanova, S.P.; Mel`nikova, O.I.
1994-11-01
An effective method for calculating the potential energy of a dinuclear system is proposed. Analytic expressions for calculating the nuclear part of the nucleus-nucleus potential in the double-folding form are obtained. A relationship between this potential and the proximity potential is found. Effects of deformation and mutual orientation of nuclei on the interaction potential are studied. It follows from a comparison of calculated potential energies with nuclear binding energies that excited states of some nuclei can be treated as dinuclear or trinuclear systems. 27 refs., 7 figs.
Brandenburg, Jan Gerit; Caldeweyher, Eike; Grimme, Stefan
2016-06-21
We extend the recently introduced PBEh-3c global hybrid density functional [S. Grimme et al., J. Chem. Phys., 2015, 143, 054107] by a screened Fock exchange variant based on the Henderson-Janesko-Scuseria exchange hole model. While the excellent performance of the global hybrid is maintained for small covalently bound molecules, its performance for computed condensed phase mass densities is further improved. Most importantly, a speed up of 30 to 50% can be achieved and especially for small orbital energy gap cases, the method is numerically much more robust. The latter point is important for many applications, e.g., for metal-organic frameworks, organic semiconductors, or protein structures. This enables an accurate density functional based electronic structure calculation of a full DNA helix structure on a single core desktop computer which is presented as an example in addition to comprehensive benchmark results. PMID:27240749
Accurate calculation of the dissociation energy of the highly anharmonic system ClHCl(-).
Stein, Christopher; Oswald, Rainer; Botschwina, Peter; Peterson, Kirk A
2015-05-28
Accurate bond dissociation energies (D0) are reported for different isotopologues of the highly anharmonic system ClHCl(-). The mass-independent equilibrium dissociation energy De was obtained by a composite method with frozen-core (fc) CCSD(T) as the basic contribution. Basis sets as large as aug-cc-pV8(+d)Z were employed, and extrapolation to the complete basis set (CBS) limit was carried out. Explicitly correlated calculations with the CCSD(T)-F12b method were also performed to support the conventionally calculated values. Core-core and core-valence correlation, scalar relativity, and higher-order correlation were considered as well. Two mass-dependent contributions, namely, the diagonal Born-Oppenheimer correction and the difference in zero-point energies between the complex and the HCl fragment, were then added in order to arrive at precise D0 values. Results for (35)ClH(35)Cl(-) and (35)ClD(35)Cl(-) are 23.81 and 23.63 kcal/mol, respectively, with estimated uncertainties of 0.05 kcal/mol. In contrast to FHF(-) ( Stein , C. ; Oswald , R. ; Sebald , P. ; Botschwina , P. ; Stoll , H. , Peterson , K. A. Mol. Phys. 2013 , 111 , 2647 - 2652 ), the D0 values of the bichloride species are larger than their De counterparts, which is an unusual situation in hydrogen-bonded systems. PMID:25405989
NASA Astrophysics Data System (ADS)
Mead, A. J.; Heymans, C.; Lombriser, L.; Peacock, J. A.; Steele, O. I.; Winther, H. A.
2016-06-01
We present an accurate non-linear matter power spectrum prediction scheme for a variety of extensions to the standard cosmological paradigm, which uses the tuned halo model previously developed in Mead et al. We consider dark energy models that are both minimally and non-minimally coupled, massive neutrinos and modified gravitational forces with chameleon and Vainshtein screening mechanisms. In all cases, we compare halo-model power spectra to measurements from high-resolution simulations. We show that the tuned halo-model method can predict the non-linear matter power spectrum measured from simulations of parametrized w(a) dark energy models at the few per cent level for k < 10 h Mpc-1, and we present theoretically motivated extensions to cover non-minimally coupled scalar fields, massive neutrinos and Vainshtein screened modified gravity models that result in few per cent accurate power spectra for k < 10 h Mpc-1. For chameleon screened models, we achieve only 10 per cent accuracy for the same range of scales. Finally, we use our halo model to investigate degeneracies between different extensions to the standard cosmological model, finding that the impact of baryonic feedback on the non-linear matter power spectrum can be considered independently of modified gravity or massive neutrino extensions. In contrast, considering the impact of modified gravity and massive neutrinos independently results in biased estimates of power at the level of 5 per cent at scales k > 0.5 h Mpc-1. An updated version of our publicly available HMCODE can be found at https://github.com/alexander-mead/hmcode.
Determination of Multidimensional Intermolecular Potential Energy Surfaces
NASA Astrophysics Data System (ADS)
Cohen, Ronald Carl
High resolution spectroscopy of the low frequency van der Waals vibrations (also referred to as Vibration -Rotation-Tunneling (VRT) spectroscopy) in weakly bound complexes provides the means to probe intermolecular forces with unprecedented detail and precision. We present an overview of the experimental information on intermolecular forces and intermolecular dynamics which has been obtained by far infrared VRT spectroscopy of 18 complexes. We then turn to a detailed examination of the Ar-H_2O complex, a simple prototype for the study of intermolecular forces. The measurement and analysis of 9 VRT bands is described. These data are first used to obtain a qualitative description of the intermolecular potential energy surface (IPS). A new simple and efficient method for calculating the eigenvalues of the multidimensional intermolecular dynamics on the IPS has been developed. This algorithm (an adaptation of the Collocation Method) was then used in a direct fit to obtain an accurate and detailed description of the intermolecular forces acting within the Ar-H_2O complex.
Exam Question Exchange: Potential Energy Surfaces.
ERIC Educational Resources Information Center
Alexander, John J., Ed.
1988-01-01
Presents three examination questions, graded in difficulty, that explore the topic of potential energy surfaces using a diagrammatic approach. Provides and discusses acceptable solutions including diagrams. (CW)
Esque, Jeremy; Cecchini, Marco
2015-04-23
The calculation of the free energy of conformation is key to understanding the function of biomolecules and has attracted significant interest in recent years. Here, we present an improvement of the confinement method that was designed for use in the context of explicit solvent MD simulations. The development involves an additional step in which the solvation free energy of the harmonically restrained conformers is accurately determined by multistage free energy perturbation simulations. As a test-case application, the newly introduced confinement/solvation free energy (CSF) approach was used to compute differences in free energy between conformers of the alanine dipeptide in explicit water. The results are in excellent agreement with reference calculations based on both converged molecular dynamics and umbrella sampling. To illustrate the general applicability of the method, conformational equilibria of met-enkephalin (5 aa) and deca-alanine (10 aa) in solution were also analyzed. In both cases, smoothly converged free-energy results were obtained in agreement with equilibrium sampling or literature calculations. These results demonstrate that the CSF method may provide conformational free-energy differences of biomolecules with small statistical errors (below 0.5 kcal/mol) and at a moderate computational cost even with a full representation of the solvent. PMID:25807150
NASA Astrophysics Data System (ADS)
Komsa, Hannu-Pekka; Berseneva, Natalia; Krasheninnikov, Arkady V.; Nieminen, Risto M.
2014-07-01
Impurities and defects frequently govern materials properties, with the most prominent example being the doping of bulk semiconductors where a minute amount of foreign atoms can be responsible for the operation of the electronic devices. Several computational schemes based on a supercell approach have been developed to get insights into types and equilibrium concentrations of point defects, which successfully work in bulk materials. Here, we show that many of these schemes cannot directly be applied to two-dimensional (2D) systems, as formation energies of charged point defects are dominated by large spurious electrostatic interactions between defects in inhomogeneous environments. We suggest two approaches that solve this problem and give accurate formation energies of charged defects in 2D systems in the dilute limit. Our methods, which are applicable to all kinds of charged defects in any 2D system, are benchmarked for impurities in technologically important h-BN and MoS2 2D materials, and they are found to perform equally well for substitutional and adatom impurities.
No Galaxy Left Behind: Accurate Measurements with the Faintest Objects in the Dark Energy Survey
Suchyta, E.
2015-07-29
Accurate statistical measurement with large imaging surveys has traditionally required throwing away a sizable fraction of the data. This is because most measurements have have relied on selecting nearly complete samples, where variations in the composition of the galaxy population with seeing, depth, or other survey characteristics are small. We introduce a new measurement method that aims to minimize this wastage, allowing precision measurement for any class of stars or galaxies detectable in an imaging survey. We have implemented our proposal in Balrog, a software package which embeds fake objects in real imaging in order to accurately characterize measurement biases. We also demonstrate this technique with an angular clustering measurement using Dark Energy Survey (DES) data. We first show that recovery of our injected galaxies depends on a wide variety of survey characteristics in the same way as the real data. We then construct a flux-limited sample of the faintest galaxies in DES, chosen specifically for their sensitivity to depth and seeing variations. Using the synthetic galaxies as randoms in the standard LandySzalay correlation function estimator suppresses the effects of variable survey selection by at least two orders of magnitude. Now our measured angular clustering is found to be in excellent agreement with that of a matched sample drawn from much deeper, higherresolution space-based COSMOS imaging; over angular scales of 0.004° < θ < 0.2 ° , we find a best-fit scaling amplitude between the DES and COSMOS measurements of 1.00 ± 0.09. We expect this methodology to be broadly useful for extending the statistical reach of measurements in a wide variety of coming imaging surveys.
No galaxy left behind. Accurate measurements with the faintest objects in the Dark Energy Survey
Suchyta, E.; Huff, E. M.; Aleksić, J.; Melchior, P.; Jouvel, S.; MacCrann, N.; Ross, A. J.; Crocce, M.; Gaztanaga, E.; Honscheid, K.; et al
2016-01-27
Accurate statistical measurement with large imaging surveys has traditionally required throwing away a sizable fraction of the data. This is because most measurements have relied on selecting nearly complete samples, where variations in the composition of the galaxy population with seeing, depth, or other survey characteristics are small. Here, we introduce a new measurement method that aims to minimize this wastage, allowing precision measurement for any class of detectable stars or galaxies. Moreover, our proposal was implemented in BALROG, software which embeds fake objects in real imaging to accurately characterize measurement biases. We demonstrate this technique with an angular clusteringmore » measurement using Dark Energy Survey (DES) data. We first show that recovery of our injected galaxies depends on a variety of survey characteristics in the same way as the real data. We then construct a flux-limited sample of the faintest galaxies in DES, chosen specifically for their sensitivity to depth and seeing variations. Using the synthetic galaxies as randoms in the Landy–Szalay estimator suppresses the effects of variable survey selection by at least two orders of magnitude. With this correction, our measured angular clustering is found to be in excellent agreement with that of a matched sample from much deeper, higher resolution space-based Cosmological Evolution Survey (COSMOS) imaging; over angular scales of 0.°004 < θ < 0.°2, we find a best-fitting scaling amplitude between the DES and COSMOS measurements of 1.00 ± 0.09. We expect this methodology to be broadly useful for extending measurements’ statistical reach in a variety of upcoming imaging surveys.« less
No Galaxy Left Behind: Accurate Measurements with the Faintest Objects in the Dark Energy Survey
Suchyta, E.
2016-01-27
Accurate statistical measurement with large imaging surveys has traditionally required throwing away a sizable fraction of the data. This is because most measurements have have relied on selecting nearly complete samples, where variations in the composition of the galaxy population with seeing, depth, or other survey characteristics are small. We introduce a new measurement method that aims to minimize this wastage, allowing precision measurement for any class of stars or galaxies detectable in an imaging survey. We have implemented our proposal in Balrog, a software package which embeds fake objects in real imaging in order to accurately characterize measurement biases.more » We also demonstrate this technique with an angular clustering measurement using Dark Energy Survey (DES) data. We first show that recovery of our injected galaxies depends on a wide variety of survey characteristics in the same way as the real data. We then construct a flux-limited sample of the faintest galaxies in DES, chosen specifically for their sensitivity to depth and seeing variations. Using the synthetic galaxies as randoms in the standard LandySzalay correlation function estimator suppresses the effects of variable survey selection by at least two orders of magnitude. Now our measured angular clustering is found to be in excellent agreement with that of a matched sample drawn from much deeper, higherresolution space-based COSMOS imaging; over angular scales of 0.004° < θ < 0.2 ° , we find a best-fit scaling amplitude between the DES and COSMOS measurements of 1.00 ± 0.09. We expect this methodology to be broadly useful for extending the statistical reach of measurements in a wide variety of coming imaging surveys.« less
No galaxy left behind: accurate measurements with the faintest objects in the Dark Energy Survey
NASA Astrophysics Data System (ADS)
Suchyta, E.; Huff, E. M.; Aleksić, J.; Melchior, P.; Jouvel, S.; MacCrann, N.; Ross, A. J.; Crocce, M.; Gaztanaga, E.; Honscheid, K.; Leistedt, B.; Peiris, H. V.; Rykoff, E. S.; Sheldon, E.; Abbott, T.; Abdalla, F. B.; Allam, S.; Banerji, M.; Benoit-Lévy, A.; Bertin, E.; Brooks, D.; Burke, D. L.; Rosell, A. Carnero; Kind, M. Carrasco; Carretero, J.; Cunha, C. E.; D'Andrea, C. B.; da Costa, L. N.; DePoy, D. L.; Desai, S.; Diehl, H. T.; Dietrich, J. P.; Doel, P.; Eifler, T. F.; Estrada, J.; Evrard, A. E.; Flaugher, B.; Fosalba, P.; Frieman, J.; Gerdes, D. W.; Gruen, D.; Gruendl, R. A.; James, D. J.; Jarvis, M.; Kuehn, K.; Kuropatkin, N.; Lahav, O.; Lima, M.; Maia, M. A. G.; March, M.; Marshall, J. L.; Miller, C. J.; Miquel, R.; Neilsen, E.; Nichol, R. C.; Nord, B.; Ogando, R.; Percival, W. J.; Reil, K.; Roodman, A.; Sako, M.; Sanchez, E.; Scarpine, V.; Sevilla-Noarbe, I.; Smith, R. C.; Soares-Santos, M.; Sobreira, F.; Swanson, M. E. C.; Tarle, G.; Thaler, J.; Thomas, D.; Vikram, V.; Walker, A. R.; Wechsler, R. H.; Zhang, Y.; DES Collaboration
2016-03-01
Accurate statistical measurement with large imaging surveys has traditionally required throwing away a sizable fraction of the data. This is because most measurements have relied on selecting nearly complete samples, where variations in the composition of the galaxy population with seeing, depth, or other survey characteristics are small. We introduce a new measurement method that aims to minimize this wastage, allowing precision measurement for any class of detectable stars or galaxies. We have implemented our proposal in BALROG, software which embeds fake objects in real imaging to accurately characterize measurement biases. We demonstrate this technique with an angular clustering measurement using Dark Energy Survey (DES) data. We first show that recovery of our injected galaxies depends on a variety of survey characteristics in the same way as the real data. We then construct a flux-limited sample of the faintest galaxies in DES, chosen specifically for their sensitivity to depth and seeing variations. Using the synthetic galaxies as randoms in the Landy-Szalay estimator suppresses the effects of variable survey selection by at least two orders of magnitude. With this correction, our measured angular clustering is found to be in excellent agreement with that of a matched sample from much deeper, higher resolution space-based Cosmological Evolution Survey (COSMOS) imaging; over angular scales of 0.004° < θ < 0.2°, we find a best-fitting scaling amplitude between the DES and COSMOS measurements of 1.00 ± 0.09. We expect this methodology to be broadly useful for extending measurements' statistical reach in a variety of upcoming imaging surveys.
Doblhoff-Dier, Katharina; Meyer, Jörg; Hoggan, Philip E; Kroes, Geert-Jan; Wagner, Lucas K
2016-06-14
Transition metals and transition metal compounds are important to catalysis, photochemistry, and many superconducting systems. We study the performance of diffusion Monte Carlo (DMC) applied to transition metal containing dimers (TMCDs) using single-determinant Slater-Jastrow trial wavefunctions and investigate the possible influence of the locality and pseudopotential errors. We find that the locality approximation can introduce nonsystematic errors of up to several tens of kilocalories per mole in the absolute energy of Cu and CuH if Ar or Mg core pseudopotentials (PPs) are used for the 3d transition metal atoms. Even for energy differences such as binding energies, errors due to the locality approximation can be problematic if chemical accuracy is sought. The use of the Ne core PPs developed by Burkatzki et al. (J. Chem. Phys. 2008, 129, 164115), the use of linear energy minimization rather than unreweighted variance minimization for the optimization of the Jastrow function, and the use of large Jastrow parametrizations reduce the locality errors. In the second section of this article, we study the general performance of DMC for 3d TMCDs using a database of binding energies of 20 TMCDs, for which comparatively accurate experimental data is available. Comparing our DMC results to these data for our results that compare best with experiment, we find a mean unsigned error (MUE) of 4.5 kcal/mol. This compares well with the achievable accuracy in CCSDT(2)Q (MUE = 4.6 kcal/mol) and the best all-electron DFT results (MUE = 4.5 kcal/mol) for the same set of systems (Truhlar et al. J. Chem. Theory Comput. 2015, 11, 2036-2052). The mean errors in DMC depend less on the exchange-correlation functionals used to generate the trial wavefunction than the corresponding mean errors in the underlying DFT calculations. Furthermore, the QMC results obtained for each molecule individually vary less with the functionals used. These observations are relevant for systems such as
Improved DFT Potential Energy Surfaces via Improved Densities.
Kim, Min-Cheol; Park, Hansol; Son, Suyeon; Sim, Eunji; Burke, Kieron
2015-10-01
Density-corrected DFT is a method that cures several failures of self-consistent semilocal DFT calculations by using a more accurate density instead. A novel procedure employs the Hartree-Fock density to bonds that are more severely stretched than ever before. This substantially increases the range of accurate potential energy surfaces obtainable from semilocal DFT for many heteronuclear molecules. We show that this works for both neutral and charged molecules. We explain why and explore more difficult cases, for example, CH(+), where density-corrected DFT results are even better than sophisticated methods like CCSD. We give a simple criterion for when DC-DFT should be more accurate than self-consistent DFT that can be applied for most cases. PMID:26722874
Energy potential of modern landfills
Bogner, J.E.
1990-01-01
Methane produced by refuse decomposition in a sanitary landfill can be recovered for commercial use. Landfill methane is currently under-utilized, with commercial recovery at only a small percentage of US landfills. New federal regulations mandating control of landfill gas migration and atmospheric emissions are providing impetus to methane recovery schemes as a means of recovering costs for increased environmental control. The benefits of landfill methane recovery include utilization of an inexpensive renewable energy resource, removal of explosive gas mixtures from the subsurface, and mitigation of observed historic increases in atmospheric methane. Increased commercial interest in landfill methane recovery is dependent on the final form of Clean Air Act amendments pertaining to gaseous emissions from landfills; market shifts in natural gas prices; financial incentives for development of renewable energy resources; and support for applied research and development to develop techniques for increased control of the gas generation process in situ. This paper will discuss the controls on methane generation in landfills. In addition, it will address how landfill regulations affect landfill design and site management practices which, in turn, influence decomposition rates. Finally, future trends in landfilling, and their relationship to gas production, will be examined. 19 refs., 2 figs., 3 tabs.
Geothermal Energy Potential in Western United States
ERIC Educational Resources Information Center
Pryde, Philip R.
1977-01-01
Reviews types of geothermal energy sources in the western states, including hot brine systems and dry steam systems. Conversion to electrical energy is a major potential use of geothermal energy, although it creates environmental disruptions such as noise, corrosion, and scaling of equipment. (AV)
LiF TLD-100 as a Dosimeter in High Energy Proton Beam Therapy-Can It Yield Accurate Results?
Zullo, John R. Kudchadker, Rajat J.; Zhu, X. Ronald; Sahoo, Narayan; Gillin, Michael T.
2010-04-01
In the region of high-dose gradients at the end of the proton range, the stopping power ratio of the protons undergoes significant changes, allowing for a broad spectrum of proton energies to be deposited within a relatively small volume. Because of the potential linear energy transfer dependence of LiF TLD-100 (thermolumescent dosimeter), dose measurements made in the distal fall-off region of a proton beam may be less accurate than those made in regions of low-dose gradients. The purpose of this study is to determine the accuracy and precision of dose measured using TLD-100 for a pristine Bragg peak, particularly in the distal fall-off region. All measurements were made along the central axis of an unmodulated 200-MeV proton beam from a Probeat passive beam-scattering proton accelerator (Hitachi, Ltd., Tokyo, Japan) at varying depths along the Bragg peak. Measurements were made using TLD-100 powder flat packs, placed in a virtual water slab phantom. The measurements were repeated using a parallel plate ionization chamber. The dose measurements using TLD-100 in a proton beam were accurate to within {+-}5.0% of the expected dose, previously seen in our past photon and electron measurements. The ionization chamber and the TLD relative dose measurements agreed well with each other. Absolute dose measurements using TLD agreed with ionization chamber measurements to within {+-} 3.0 cGy, for an exposure of 100 cGy. In our study, the differences in the dose measured by the ionization chamber and those measured by TLD-100 were minimal, indicating that the accuracy and precision of measurements made in the distal fall-off region of a pristine Bragg peak is within the expected range. Thus, the rapid change in stopping power ratios at the end of the range should not affect such measurements, and TLD-100 may be used with confidence as an in vivo dosimeter for proton beam therapy.
Wave Energy Potential in the Latvian EEZ
NASA Astrophysics Data System (ADS)
Beriņš, J.; Beriņš, J.; Kalnačs, J.; Kalnačs, A.
2016-06-01
The present article deals with one of the alternative forms of energy - sea wave energy potential in the Latvian Exclusice Economic Zone (EEZ). Results have been achieved using a new method - VEVPP. Calculations have been performed using the data on wave parameters over the past five years (2010-2014). We have also considered wave energy potential in the Gulf of Riga. The conclusions have been drawn on the recommended methodology for the sea wave potential and power calculations for wave-power plant pre-design stage.
A global potential energy surface for ArH2
NASA Technical Reports Server (NTRS)
Schwenke, David W.; Walch, Stephen P.; Taylor, Peter R.
1992-01-01
We describe a simple analytic representation of the ArH2 potential energy surface which well reproduces the results of extensive ab initio electronic structure calculations. The analytic representation smoothly interpolates between the dissociated H2 and strong bonding limits. In the fitting process, emphasis is made on accurately reproducing regions of the potential expected to be important for high temperature (ca. 3000 K) collision processes. Overall, the anisotropy and H2 bond length dependence of the analytic representation well reproduce the input data.
Theoretical studies of potential energy surfaces and computational methods.
Shepard, R.
2006-01-01
This project involves the development, implementation, and application of theoretical methods for the calculation and characterization of potential energy surfaces (PES) involving molecular species that occur in hydrocarbon combustion. These potential energy surfaces require an accurate and balanced treatment of reactants, intermediates, and products. Most of our work focuses on general multiconfiguration self-consistent-field (MCSCF) and multireference single- and double-excitation configuration interaction (MRSDCI) methods. In contrast to the more common single-reference electronic structure methods, this approach is capable of describing accurately molecular systems that are highly distorted away from their equilibrium geometries, including reactant, fragment, and transition-state geometries, and of describing regions of the potential surface that are associated with electronic wave functions of widely varying nature. The MCSCF reference wave functions are designed to be sufficiently flexible to describe qualitatively the changes in the electronic structure over the broad range of molecular geometries of interest. The necessary mixing of ionic, covalent, and Rydberg contributions, along with the appropriate treatment of the different electron-spin components (e.g. closed shell, high-spin open-shell, low-spin open shell, radical, diradical, etc.) of the wave functions are treated correctly at this level. Further treatment of electron correlation effects is included using large scale multireference CI wave functions, particularly including the single and double excitations relative to the MCSCF reference space. This leads to the most flexible and accurate large-scale MRSDCI wave functions that have been used to date in global PES studies.
Theoretical studies of potential energy surfaces and computational methods
Shepard, R.
1993-12-01
This project involves the development, implementation, and application of theoretical methods for the calculation and characterization of potential energy surfaces involving molecular species that occur in hydrocarbon combustion. These potential energy surfaces require an accurate and balanced treatment of reactants, intermediates, and products. This difficult challenge is met with general multiconfiguration self-consistent-field (MCSCF) and multireference single- and double-excitation configuration interaction (MRSDCI) methods. In contrast to the more common single-reference electronic structure methods, this approach is capable of describing accurately molecular systems that are highly distorted away from their equilibrium geometries, including reactant, fragment, and transition-state geometries, and of describing regions of the potential surface that are associated with electronic wave functions of widely varying nature. The MCSCF reference wave functions are designed to be sufficiently flexible to describe qualitatively the changes in the electronic structure over the broad range of geometries of interest. The necessary mixing of ionic, covalent, and Rydberg contributions, along with the appropriate treatment of the different electron-spin components (e.g. closed shell, high-spin open-shell, low-spin open shell, radical, diradical, etc.) of the wave functions, are treated correctly at this level. Further treatment of electron correlation effects is included using large scale multireference CI wave functions, particularly including the single and double excitations relative to the MCSCF reference space. This leads to the most flexible and accurate large-scale MRSDCI wave functions that have been used to date in global PES studies.
Calculation of molecular free energies in classical potentials
NASA Astrophysics Data System (ADS)
Farhi, Asaf; Singh, Bipin
2016-02-01
Free energies of molecules can be calculated by quantum chemistry computations or by normal mode classical calculations. However, the first can be computationally impractical for large molecules and the second is based on the assumption of harmonic dynamics. We present a novel, accurate and complete calculation of molecular free energies in standard classical potentials. In this method we transform the molecule by relaxing potential terms which depend on the coordinates of a group of atoms in that molecule and calculate the free energy difference associated with the transformation. Then, since the transformed molecule can be treated as non-interacting systems, the free energy associated with these atoms is analytically or numerically calculated. This two-step calculation can be applied to calculate free energies of molecules or free energy difference between (possibly large) molecules in a general environment. We demonstrate the method in free energy calculations for methanethiol and butane molecules in vacuum and solvent. We suggest the potential application of free energy calculation of chemical reactions in classical molecular simulations.
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.
Global potential energy hypersurface for dynamical studies of energy transfer in HF--HF collisions
Redmon, M.J.; Binkley, J.S.
1987-07-15
The interaction energy of two HF molecules at 1332 individual points has been calculated with Moeller--Plesset (many--body) perturbation theory at the MP4-SDTQ level using a 6-311G** basis set. 293 of the points correspond to stretching of one HF molecule from its equilibrium geometry. No attempt was made to use a sufficiently fine grid to accurately describe the well region corresponding to hydrogen bonding. However, the location and minimum energy are consistent with experiment and other accurate theoretical results. An extensive global fit (rms error of 1 kcal/mol) is reported of 1319 points (below 10 eV of potential energy) using a modified London potential with corrections obtained using polynomials through four-body interactions. A model electrostatic potential represents the long-range interaction. In addition, the use of an expansion in products of three Legendre functions is discussed. It is shown that the latter approach, although accurately fitting the ab initio data, has difficulties interpolating in regions of the surface exhibiting diverse magnitudes of potential energy, and therefore must be used with caution. This surface should be useful for studies of T--V--R processes in this system.
Wind energy in China: Estimating the potential
NASA Astrophysics Data System (ADS)
Yuan, Jiahai
2016-07-01
Persistent and significant curtailment has cast concern over the prospects of wind power in China. A comprehensive assessment of the production of energy from wind has identified grid-integrated wind generation potential at 11.9–14% of China's projected energy demand by 2030.
Potential of energy production from conserved forages
Technology Transfer Automated Retrieval System (TEKTRAN)
Forages have a potential role in meeting the demand for energy. Perennial forages are attractive for various reasons. One, both the monetary and energy cost of planting is spread over many years. Two, we already have the equipment for harvesting, storing and transporting this source of biomass. Thre...
Theoretical studies of potential energy surfaces
Harding, L.B.
1995-07-01
MRCI (configuration interaction) calculations were used to examine possible pathways for the O{sub 2} + CCH reaction. The H{sub 2} + CN potential surface was examined. An initial survey was made of the HCl + CN potential energy surface at a low level of theory.
Economic Energy Savings Potential in Federal Buildings
Brown, Daryl R.; Dirks, James A.; Hunt, Diane M.
2000-09-04
The primary objective of this study was to estimate the current life-cycle cost-effective (i.e., economic) energy savings potential in Federal buildings and the corresponding capital investment required to achieve these savings, with Federal financing. Estimates were developed for major categories of energy efficiency measures such as building envelope, heating system, cooling system, and lighting. The analysis was based on conditions (building stock and characteristics, retrofit technologies, interest rates, energy prices, etc.) existing in the late 1990s. The potential impact of changes to any of these factors in the future was not considered.
Energy expenditure during level human walking: seeking a simple and accurate predictive solution.
Ludlow, Lindsay W; Weyand, Peter G
2016-03-01
Accurate prediction of the metabolic energy that walking requires can inform numerous health, bodily status, and fitness outcomes. We adopted a two-step approach to identifying a concise, generalized equation for predicting level human walking metabolism. Using literature-aggregated values we compared 1) the predictive accuracy of three literature equations: American College of Sports Medicine (ACSM), Pandolf et al., and Height-Weight-Speed (HWS); and 2) the goodness-of-fit possible from one- vs. two-component descriptions of walking metabolism. Literature metabolic rate values (n = 127; speed range = 0.4 to 1.9 m/s) were aggregated from 25 subject populations (n = 5-42) whose means spanned a 1.8-fold range of heights and a 4.2-fold range of weights. Population-specific resting metabolic rates (V̇o2 rest) were determined using standardized equations. Our first finding was that the ACSM and Pandolf et al. equations underpredicted nearly all 127 literature-aggregated values. Consequently, their standard errors of estimate (SEE) were nearly four times greater than those of the HWS equation (4.51 and 4.39 vs. 1.13 ml O2·kg(-1)·min(-1), respectively). For our second comparison, empirical best-fit relationships for walking metabolism were derived from the data set in one- and two-component forms for three V̇o2-speed model types: linear (∝V(1.0)), exponential (∝V(2.0)), and exponential/height (∝V(2.0)/Ht). We found that the proportion of variance (R(2)) accounted for, when averaged across the three model types, was substantially lower for one- vs. two-component versions (0.63 ± 0.1 vs. 0.90 ± 0.03) and the predictive errors were nearly twice as great (SEE = 2.22 vs. 1.21 ml O2·kg(-1)·min(-1)). Our final analysis identified the following concise, generalized equation for predicting level human walking metabolism: V̇o2 total = V̇o2 rest + 3.85 + 5.97·V(2)/Ht (where V is measured in m/s, Ht in meters, and V̇o2 in ml O2·kg(-1)·min(-1)). PMID:26679617
Beste, Ariana; Harrison, Robert J; Yanai, Takeshi
2006-01-01
Chemists are mainly interested in energy differences. In contrast, most quantum chemical methods yield the total energy which is a large number compared to the difference and has therefore to be computed to a higher relative precision than would be necessary for the difference alone. Hence, it is desirable to compute energy differences directly, thereby avoiding the precision problem. Whenever it is possible to find a parameter which transforms smoothly from an initial to a final state, the energy difference can be obtained by integrating the energy derivative with respect to that parameter (c.f., thermodynamic integration or adiabatic connection methods). If the dependence on the parameter is predominantly linear, accurate results can be obtained by single-point integration. In density functional theory (DFT) and Hartree-Fock, we applied the formalism to ionization potentials, excitation energies, and chemical bond breaking. Example calculations for ionization potentials and excitation energies showed that accurate results could be obtained with a linear estimate. For breaking bonds, we introduce a non-geometrical parameter which gradually turns the interaction between two fragments of a molecule on. The interaction changes the potentials used to determine the orbitals as well as constraining the orbitals to be orthogonal.
NASA Astrophysics Data System (ADS)
Beste, A.; Harrison, R. J.; Yanai, T.
2006-08-01
Chemists are mainly interested in energy differences. In contrast, most quantum chemical methods yield the total energy which is a large number compared to the difference and has therefore to be computed to a higher relative precision than would be necessary for the difference alone. Hence, it is desirable to compute energy differences directly, thereby avoiding the precision problem. Whenever it is possible to find a parameter which transforms smoothly from an initial to a final state, the energy difference can be obtained by integrating the energy derivative with respect to that parameter (cf. thermodynamic integration or adiabatic connection methods). If the dependence on the parameter is predominantly linear, accurate results can be obtained by single-point integration. In density functional theory and Hartree-Fock, we applied the formalism to ionization potentials, excitation energies, and chemical bond breaking. Example calculations for ionization potentials and excitation energies showed that accurate results could be obtained with a linear estimate. For breaking bonds, we introduce a nongeometrical parameter which gradually turns the interaction between two fragments of a molecule on. The interaction changes the potentials used to determine the orbitals as well as the constraint on the orbitals to be orthogonal.
Technology Transfer Automated Retrieval System (TEKTRAN)
The three evapotranspiration (ET) measurement/retrieval techniques used in this study, lysimeter, scintillometer and remote sensing vary in their level of complexity, accuracy, resolution and applicability. The lysimeter with its point measurement is the most accurate and direct method to measure ET...
Rucker, Robert; Oelschlaeger, Peter; Warshel, Arieh
2010-01-01
DNA polymerase β (pol β) is a small eukaryotic enzyme with the ability to repair short single-stranded DNA gaps that has found use as a model system for larger replicative DNA polymerases. For all DNA polymerases, the factors determining their catalytic power and fidelity are the interactions between the bases of the base pair, amino acids near the active site, and the two magnesium ions. In this report, we study effects of all three aspects on human pol β transition state (TS) binding free energies by reproducing a consistent set of experimentally determined data for different structures. Our calculations comprise the combination of four different base pairs (incoming pyrimidine nucleotides incorporated opposite both matched and mismatched purines) with four different pol β structures (wild type and three separate mutations of ionized residues to alanine). We decompose the incoming deoxynucleoside 5′-triphosphate-TS, and run separate calculations for the neutral base part and the highly charged triphosphate part, using different dielectric constants in order to account for the specific electrostatic environments. This new approach improves our ability to predict the effect of matched and mismatched base pairing and of mutations in DNA polymerases on fidelity and may be a useful tool in studying the potential of DNA polymerase mutations in the development of cancer. It also supports our point of view with regards to the origin of the structural control of fidelity, allowing for a quantified description of the fidelity of DNA polymerases. PMID:19842163
NREL's Field Data Repository Supports Accurate Home Energy Analysis (Fact Sheet)
None, None
2012-02-01
This fact sheet discusses NREL's work to develop a repository of research-level residential building characteristics and historical energy use data to support ongoing efforts to improve the accuracy of residential energy analysis tools and the efficiency of energy assessment processes. The objective of this project is to create a robust empirical data source to support the research goals of the Department of Energy's Building America program, which is to improve the efficiency of existing U.S. homes by 30% to 50%. Researchers can use this data source to test the accuracy of building energy simulation software and energy audit procedures, ultimately leading to more credible and less expensive energy analysis.
Laurens, Lieve M L; Quinn, Matthew; Van Wychen, Stefanie; Templeton, David W; Wolfrum, Edward J
2012-04-01
In the context of algal biofuels, lipids, or better aliphatic chains of the fatty acids, are perhaps the most important constituents of algal biomass. Accurate quantification of lipids and their respective fuel yield is crucial for comparison of algal strains and growth conditions and for process monitoring. As an alternative to traditional solvent-based lipid extraction procedures, we have developed a robust whole-biomass in situ transesterification procedure for quantification of algal lipids (as fatty acid methyl esters, FAMEs) that (a) can be carried out on a small scale (using 4-7 mg of biomass), (b) is applicable to a range of different species, (c) consists of a single-step reaction, (d) is robust over a range of different temperature and time combinations, and (e) tolerant to at least 50% water in the biomass. Unlike gravimetric lipid quantification, which can over- or underestimate the lipid content, whole biomass transesterification reflects the true potential fuel yield of algal biomass. We report here on the comparison of the yield of FAMEs by using different catalysts and catalyst combinations, with the acid catalyst HCl providing a consistently high level of conversion of fatty acids with a precision of 1.9% relative standard deviation. We investigate the influence of reaction time, temperature, and biomass water content on the measured FAME content and profile for 4 different samples of algae (replete and deplete Chlorella vulgaris, replete Phaeodactylum tricornutum, and replete Nannochloropsis sp.). We conclude by demonstrating a full mass balance closure of all fatty acids around a traditional lipid extraction process. PMID:22349344
Laurens, L. M. L.; Quinn, M.; Van Wychen, S.; Templeton, D. W.; Wolfrum, E. J.
2012-04-01
In the context of algal biofuels, lipids, or better aliphatic chains of the fatty acids, are perhaps the most important constituents of algal biomass. Accurate quantification of lipids and their respective fuel yield is crucial for comparison of algal strains and growth conditions and for process monitoring. As an alternative to traditional solvent-based lipid extraction procedures, we have developed a robust whole-biomass in situ transesterification procedure for quantification of algal lipids (as fatty acid methyl esters, FAMEs) that (a) can be carried out on a small scale (using 4-7 mg of biomass), (b) is applicable to a range of different species, (c) consists of a single-step reaction, (d) is robust over a range of different temperature and time combinations, and (e) tolerant to at least 50% water in the biomass. Unlike gravimetric lipid quantification, which can over- or underestimate the lipid content, whole biomass transesterification reflects the true potential fuel yield of algal biomass. We report here on the comparison of the yield of FAMEs by using different catalysts and catalyst combinations, with the acid catalyst HCl providing a consistently high level of conversion of fatty acids with a precision of 1.9% relative standard deviation. We investigate the influence of reaction time, temperature, and biomass water content on the measured FAME content and profile for 4 different samples of algae (replete and deplete Chlorella vulgaris, replete Phaeodactylum tricornutum, and replete Nannochloropsis sp.). We conclude by demonstrating a full mass balance closure of all fatty acids around a traditional lipid extraction process.
Computed rotational rainbows from realistic potential energy surfaces
Gianturco, F.A.; Palma, A.
1985-08-01
The quantal IOS approximation in here employed to study interference structures in the rotationally inelastic, state-to-state differential cross sections for polar diatomic targets (LiH, FH, and CO) interacting with He atoms. Quite realistic expressions are used to describe the relevant potential energy surfaces (PES) which were taken from previous works that tested them against accurate experimental findings for total and partial differential cross sections. Specific features like short-range anisotropy and well depth, long-range attractive regions and overall range of action for each potential employed are analyzed and discussed in relation to their influence on rotational rainbows appearance and on the possible observation of cross section extrema in rotational energy distributions.
Terahertz absorption spectra and potential energy distribution of liquid crystals.
Chen, Zezhang; Jiang, Yurong; Jiang, Lulu; Ma, Heng
2016-01-15
In this work, the terahertz (THz) absorption spectra of a set of nematic liquid crystals were studied using the density functional theories (DFT). An accurate assignment of the vibrational modes corresponding to absorption frequencies were performed using potential energy distribution (PED) in a frequency range of 0-3 THz. The impacts of different core structures on THz absorption spectra were discussed. The results indicate that scope of application must be considered in the LC-based THz device designing. This proposed work may give a useful suggestion on the design of novel liquid crystal material in THz wave. PMID:26476072
Potential energy function for the hydroperoxyl radical
Lemon, W.J.; Hase, W.L.
1987-03-12
A switching function formalism is used to derive an analytic potential energy surface for the O + OH in equilibrium HO/sub 2/ in equilibrium H + O/sub 2/ reactive system. Both experimental and ab initio data are used to derive parameters for the potential energy surface. Trajectory calculations for highly excited HO/sub 2/ are performed on this surface. From these trajectories quasi-periodic eigentrajectories are found for vibrational levels near the HO/sub 2/ dissociation threshold with small amounts of quanta in the OH stretch mode and large amounts of quanta in the OO stretch mode.
Geothermal energy and its potential. [Utah
Berge, C.W.
1980-06-01
A brief review of geothermal energy and its potential as a future energy source is presented. The type of geothermal systems and their geologic occurrence is discussed, and the Phillips Petroleum Company's exploration and drilling programs in the Roosevelt Hot Springs area in parts of Iron, Beaver, and Millard Counties, Utah are detailed. A section on the rock behavior and mechanical properties of rocks in the Roosevelt Hot Springs area is included. (JMT)
Intermolecular potential energy surface and thermophysical properties of ethylene oxide
Crusius, Johann-Philipp Hassel, Egon; Hellmann, Robert; Bich, Eckard
2014-10-28
A six-dimensional potential energy hypersurface (PES) for two interacting rigid ethylene oxide (C{sub 2}H{sub 4}O) molecules was determined from high-level quantum-chemical ab initio calculations. The counterpoise-corrected supermolecular approach at the MP2 and CCSD(T) levels of theory was utilized to determine interaction energies for 10178 configurations of two molecules. An analytical site-site potential function with 19 sites per ethylene oxide molecule was fitted to the interaction energies and fine tuned to agree with data for the second acoustic virial coefficient from accurate speed of sound measurements. The PES was validated by computing the second virial coefficient, shear viscosity, and thermal conductivity. The values of these properties are substantiated by the best experimental data as they tend to fall within the uncertainty intervals and also obey the experimental temperature functions, except for viscosity, where experimental data are insufficient. Due to the lack of reliable data, especially for the transport properties, our calculated values are currently the most accurate estimates for these properties of ethylene oxide.
Intermolecular potential energy surface and thermophysical properties of ethylene oxide.
Crusius, Johann-Philipp; Hellmann, Robert; Hassel, Egon; Bich, Eckard
2014-10-28
A six-dimensional potential energy hypersurface (PES) for two interacting rigid ethylene oxide (C2H4O) molecules was determined from high-level quantum-chemical ab initio calculations. The counterpoise-corrected supermolecular approach at the MP2 and CCSD(T) levels of theory was utilized to determine interaction energies for 10178 configurations of two molecules. An analytical site-site potential function with 19 sites per ethylene oxide molecule was fitted to the interaction energies and fine tuned to agree with data for the second acoustic virial coefficient from accurate speed of sound measurements. The PES was validated by computing the second virial coefficient, shear viscosity, and thermal conductivity. The values of these properties are substantiated by the best experimental data as they tend to fall within the uncertainty intervals and also obey the experimental temperature functions, except for viscosity, where experimental data are insufficient. Due to the lack of reliable data, especially for the transport properties, our calculated values are currently the most accurate estimates for these properties of ethylene oxide. PMID:25362314
Biomass resource potential using energy crops
Wright, L.L.; Cushman, J.H.; Martin, S.A.
1993-09-01
Biomass energy crops can provide a significant and environmentally beneficial source of renewable energy feedstocks for the future. They can revitalize the agricultural sector of the US economy by providing profitable uses for marginal cropland. Energy crops include fast-growing trees, perennial grasses, and annual grasses, all capable of collecting solar energy and storing it as cellulosic compounds for several months to several years. Once solar energy is thus captured, it can be converted by means of currently available technologies to a wide variety of energy products such as electricity, heat, liquid transportation fuels, and gases. Experimental results from field trials have generated optimism that selected and improved energy crops, established on cropland with moderate limitations for crop production, have the potential for producing high yields. Both trees and grasses, under very good growing conditions, have produced average annual yields of 20 to 40 dry Mg ha{sup {minus}1} year{sup {minus}1}. Sorghum has shown especially high yields in the Midwest. Hybrids between sugar cane and its wild relatives, called energy cane, have yielded as much as 50 dry Mg ha{sup {minus}1} year{sup {minus}1} in Florida. These experimental results demonstrate that some species have the genetic potential for very rapid growth rates. New wood energy crop systems developed by the Department of Energy`s Biofuels Feedstock Development Program offer, at a minimum, a 100% increase in biomass production rates over the 2 to 4 Mg ha{sup {minus}1} year{sup {minus}1} of dry leafless woody biomass produced by most natural forest systems. Experimental data indicate that short rotation wood crops established on cropland with moderate limitations are capable of producing biomass yields of 8--20 dry Mg ha{sup {minus}1} year{sup {minus}1} with a present average about 11 dry Mg ha{sup {minus}1} year{sup {minus}1} on typical cropland sites.
Potential reduction of DSN uplink energy cost
NASA Technical Reports Server (NTRS)
Dolinsky, S.; Degroot, N. F.
1982-01-01
DSN Earth stations typically transmit more power than that required to meet minimum specifications for uplink performance. Energy and cost savings that could result from matching the uplink power to the amount required for specified performance are studied. The Galileo mission was selected as a case study. Although substantial reduction in transmitted energy is possible, potential savings in source energy (oil or electricity) savings are much less. This is because of the rising inefficiency in power conversion and radio frequency power generation that accompanies reduced power output.
Potential energy studies on silane dimers
NASA Astrophysics Data System (ADS)
Mahlanen, Riina; Pakkanen, Tapani A.
2011-04-01
Intermolecular interactions and parameters for use in MD studies of large molecule systems have earlier been determined for hydrocarbons, carbon tetrahalides and sulfur. The paper reports a model representing nonbonding interactions between silane molecules, which were examined in the same way as hydrocarbons in an earlier (neopentane, isopropane, propane, and ethane) study. Intermolecular potentials were determined for 11 combinations of silane compound pairs (silane SiH 4, disilane Si 2H 6, trisilane Si 3H 8, isotetrasilane Si 4H 10 and neopentasilane Si 5H 12) with MP2/aug(df)-6-311G ∗ab initio calculations. The most stable dimer configurations were identified. With use of the modified Morse potential model to represent the interactions, 276 new potential energy surfaces were generated for silane dimers. Separate and generic pair potentials were calculated for the silanes. The pair potentials can be used in MD studies of silanes.
NASA Astrophysics Data System (ADS)
Chiarelli, Antonio M.; Maclin, Edward L.; Low, Kathy A.; Mathewson, Kyle E.; Fabiani, Monica; Gratton, Gabriele
2016-03-01
Diffuse optical tomography (DOT) provides data about brain function using surface recordings. Despite recent advancements, an unbiased method for estimating the depth of absorption changes and for providing an accurate three-dimensional (3-D) reconstruction remains elusive. DOT involves solving an ill-posed inverse problem, requiring additional criteria for finding unique solutions. The most commonly used criterion is energy minimization (energy constraint). However, as measurements are taken from only one side of the medium (the scalp) and sensitivity is greater at shallow depths, the energy constraint leads to solutions that tend to be small and superficial. To correct for this bias, we combine the energy constraint with another criterion, minimization of spatial derivatives (Laplacian constraint, also used in low resolution electromagnetic tomography, LORETA). Used in isolation, the Laplacian constraint leads to solutions that tend to be large and deep. Using simulated, phantom, and actual brain activation data, we show that combining these two criteria results in accurate (error <2 mm) absorption depth estimates, while maintaining a two-point spatial resolution of <24 mm up to a depth of 30 mm. This indicates that accurate 3-D reconstruction of brain activity up to 30 mm from the scalp can be obtained with DOT.
Sahu, Nityananda; Singh, Gurmeet; Nandi, Apurba; Gadre, Shridhar R
2016-07-21
Owing to the steep scaling behavior, highly accurate CCSD(T) calculations, the contemporary gold standard of quantum chemistry, are prohibitively difficult for moderate- and large-sized water clusters even with the high-end hardware. The molecular tailoring approach (MTA), a fragmentation-based technique is found to be useful for enabling such high-level ab initio calculations. The present work reports the CCSD(T) level binding energies of many low-lying isomers of large (H2O)n (n = 16, 17, and 25) clusters employing aug-cc-pVDZ and aug-cc-pVTZ basis sets within the MTA framework. Accurate estimation of the CCSD(T) level binding energies [within 0.3 kcal/mol of the respective full calculation (FC) results] is achieved after effecting the grafting procedure, a protocol for minimizing the errors in the MTA-derived energies arising due to the approximate nature of MTA. The CCSD(T) level grafting procedure presented here hinges upon the well-known fact that the MP2 method, which scales as O(N(5)), can be a suitable starting point for approximating to the highly accurate CCSD(T) [that scale as O(N(7))] energies. On account of the requirement of only an MP2-level FC on the entire cluster, the current methodology ultimately leads to a cost-effective solution for the CCSD(T) level accurate binding energies of large-sized water clusters even at the complete basis set limit utilizing off-the-shelf hardware. PMID:27351269
Pseudospectral Gaussian quantum dynamics: Efficient sampling of potential energy surfaces
NASA Astrophysics Data System (ADS)
Heaps, Charles W.; Mazziotti, David A.
2016-04-01
Trajectory-based Gaussian basis sets have been tremendously successful in describing high-dimensional quantum molecular dynamics. In this paper, we introduce a pseudospectral Gaussian-based method that achieves accurate quantum dynamics using efficient, real-space sampling of the time-dependent basis set. As in other Gaussian basis methods, we begin with a basis set expansion using time-dependent Gaussian basis functions guided by classical mechanics. Unlike other Gaussian methods but characteristic of the pseudospectral and collocation methods, the basis set is tested with N Dirac delta functions, where N is the number of basis functions, rather than using the basis function as test functions. As a result, the integration for matrix elements is reduced to function evaluation. Pseudospectral Gaussian dynamics only requires O ( N ) potential energy calculations, in contrast to O ( N 2 ) evaluations in a variational calculation. The classical trajectories allow small basis sets to sample high-dimensional potentials. Applications are made to diatomic oscillations in a Morse potential and a generalized version of the Henon-Heiles potential in two, four, and six dimensions. Comparisons are drawn to full analytical evaluation of potential energy integrals (variational) and the bra-ket averaged Taylor (BAT) expansion, an O ( N ) approximation used in Gaussian-based dynamics. In all cases, the pseudospectral Gaussian method is competitive with full variational calculations that require a global, analytical, and integrable potential energy surface. Additionally, the BAT breaks down when quantum mechanical coherence is particularly strong (i.e., barrier reflection in the Morse oscillator). The ability to obtain variational accuracy using only the potential energy at discrete points makes the pseudospectral Gaussian method a promising avenue for on-the-fly dynamics, where electronic structure calculations become computationally significant.
Pseudospectral Gaussian quantum dynamics: Efficient sampling of potential energy surfaces.
Heaps, Charles W; Mazziotti, David A
2016-04-28
Trajectory-based Gaussian basis sets have been tremendously successful in describing high-dimensional quantum molecular dynamics. In this paper, we introduce a pseudospectral Gaussian-based method that achieves accurate quantum dynamics using efficient, real-space sampling of the time-dependent basis set. As in other Gaussian basis methods, we begin with a basis set expansion using time-dependent Gaussian basis functions guided by classical mechanics. Unlike other Gaussian methods but characteristic of the pseudospectral and collocation methods, the basis set is tested with N Dirac delta functions, where N is the number of basis functions, rather than using the basis function as test functions. As a result, the integration for matrix elements is reduced to function evaluation. Pseudospectral Gaussian dynamics only requires O(N) potential energy calculations, in contrast to O(N(2)) evaluations in a variational calculation. The classical trajectories allow small basis sets to sample high-dimensional potentials. Applications are made to diatomic oscillations in a Morse potential and a generalized version of the Henon-Heiles potential in two, four, and six dimensions. Comparisons are drawn to full analytical evaluation of potential energy integrals (variational) and the bra-ket averaged Taylor (BAT) expansion, an O(N) approximation used in Gaussian-based dynamics. In all cases, the pseudospectral Gaussian method is competitive with full variational calculations that require a global, analytical, and integrable potential energy surface. Additionally, the BAT breaks down when quantum mechanical coherence is particularly strong (i.e., barrier reflection in the Morse oscillator). The ability to obtain variational accuracy using only the potential energy at discrete points makes the pseudospectral Gaussian method a promising avenue for on-the-fly dynamics, where electronic structure calculations become computationally significant. PMID:27131532
Potential energy savings from aquifer thermal energy storage
Anderson, M.R.; Weijo, R.O.
1988-07-01
Pacific Northwest Laboratory researchers developed an aggregate-level model to estimate the short- and long-term potential energy savings from using aquifer thermal storage (ATES) in the United States. The objectives of this effort were to (1) develop a basis from which to recommend whether heat or chill ATES should receive future research focus and (2) determine which market sector (residential, commercial, or industrial) offers the largest potential energy savings from ATES. Information was collected on the proportion of US land area suitable for ATES applications. The economic feasibility of ATES applications was then evaluated. The potential energy savings from ATES applications was calculated. Characteristic energy use in the residential, commercial, and industrial sectors was examined, as was the relationship between waste heat production and consumption by industrial end-users. These analyses provided the basis for two main conclusions: heat ATES applications offer higher potential for energy savings than do chill ATES applications; and the industrial sector can achieve the highest potential energy savings for the large consumption markets. Based on these findings, it is recommended that future ATES research and development efforts be directed toward heat ATES applications in the industrial sector. 11 refs., 6 figs., 9 tabs.
Global Potential of Energy Efficiency Standards and Labeling Programs
McNeil, Michael A; McNeil, Michael A.; Letschert, Virginie; de la Rue du Can, Stephane
2008-06-15
This report estimates the global potential reductions in greenhouse gas emissions by 2030 for energy efficiency improvements associated with equipment (appliances, lighting, and HVAC) in buildings by means of energy efficiency standards and labels (EES&L). A consensus has emerged among the world's scientists and many corporate and political leaders regarding the need to address the threat of climate change through emissions mitigation and adaptation. A further consensus has emerged that a central component of these strategies must be focused around energy, which is the primary generator of greenhouse gas emissions. Two important questions result from this consensus: 'what kinds of policies encourage the appropriate transformation to energy efficiency' and 'how much impact can these policies have'? This report aims to contribute to the dialogue surrounding these issues by considering the potential impacts of a single policy type, applied on a global scale. The policy addressed in this report is Energy Efficient Standards and Labeling (EES&L) for energy-consuming equipment, which has now been implemented in over 60 countries. Mandatory energy performance standards are important because they contribute positively to a nation's economy and provide relative certainty about the outcome (both timing and magnitudes). Labels also contribute positively to a nation's economy and importantly increase the awareness of the energy-consuming public. Other policies not analyzed here (utility incentives, tax credits) are complimentary to standards and labels and also contribute in significant ways to reducing greenhouse gas emissions. We believe the analysis reported here to be the first systematic attempt to evaluate the potential of savings from EES&L for all countries and for such a large set of products. The goal of the analysis is to provide an assessment that is sufficiently well-quantified and accurate to allow comparison and integration with other strategies under
NASA Astrophysics Data System (ADS)
Artrith, Nongnuch; Kolpak, Alexie
2014-03-01
The shape, size, and composition of catalyst nanoparticles can have a significant influence on their catalytic activity. Understanding such structure-reactivity relationships is crucial for the optimization of industrial catalysts and the design of novel catalysts with enhanced properties. In this work, we investigate the equilibrium shape and surface structure/composition of Au/Cu nanoparticles in solution, which have recently been shown to be stable and efficient catalysts for CO2 reduction. Using a combination of density functional theory calculations and large-scale Monte-Carlo and molecular dynamics simulations with reactive atomistic potentials, we determine how the nanoparticle shape, surface structure, and surface stoichiometry (i.e., fraction of Au at the surface relative to overall composition), evolve as a function of varying catalytic conditions. We discuss the effects of these changes on the surface electronic structure and binding energies of CO2, H2, and CH3OH. Our results emphasize the important relationships between catalytic environment (e.g., solvent effects), catalyst structure, and catalytic activity. We thank the Schlumberger Foundation Faculty for the Future for financial support. Computing time at XSEDE and NERSC clusters are gratefully acknowledged.
Nuclear binding energy and symmetry energy of nuclear matter with modern nucleon-nucleon potentials
Hassaneen, Kh.S.A.; Abo-Elsebaa, H.M.; Sultan, E.A.; Mansour, H.M.M.
2011-03-15
Research Highlights: > The nuclear matter is studied within the Brueckner-Hartree-Fock (BHF) approach employing the most recent accurate nucleon-nucleon potentials. > The results come out by approximating the single particle self-consistent potential with a parabolic form. > We discuss the current status of the Coester line, i.e., density and energy of the various saturation points being strongly linearly correlated. > The nuclear symmetry energy is calculated as the difference between the binding energy of pure neutron matter and that of symmetric nuclear matter. - Abstract: The binding energy of nuclear matter at zero temperature in the Brueckner-Hartree-Fock approximation with modern nucleon-nucleon potentials is studied. Both the standard and continuous choices of single particle energies are used. These modern nucleon-nucleon potentials fit the deuteron properties and are phase shifts equivalent. Comparison with other calculations is made. In addition we present results for the symmetry energy obtained with different potentials, which is of great importance in astrophysical calculation.
An Energy-Efficient Strategy for Accurate Distance Estimation in Wireless Sensor Networks
Tarrío, Paula; Bernardos, Ana M.; Casar, José R.
2012-01-01
In line with recent research efforts made to conceive energy saving protocols and algorithms and power sensitive network architectures, in this paper we propose a transmission strategy to minimize the energy consumption in a sensor network when using a localization technique based on the measurement of the strength (RSS) or the time of arrival (TOA) of the received signal. In particular, we find the transmission power and the packet transmission rate that jointly minimize the total consumed energy, while ensuring at the same time a desired accuracy in the RSS or TOA measurements. We also propose some corrections to these theoretical results to take into account the effects of shadowing and packet loss in the propagation channel. The proposed strategy is shown to be effective in realistic scenarios providing energy savings with respect to other transmission strategies, and also guaranteeing a given accuracy in the distance estimations, which will serve to guarantee a desired accuracy in the localization result. PMID:23202218
Asmadi, Aldi; Neumann, Marcus A; Kendrick, John; Girard, Pascale; Perrin, Marc-Antoine; Leusen, Frank J J
2009-12-24
In the 2007 blind test of crystal structure prediction hosted by the Cambridge Crystallographic Data Centre (CCDC), a hybrid DFT/MM method correctly ranked each of the four experimental structures as having the lowest lattice energy of all the crystal structures predicted for each molecule. The work presented here further validates this hybrid method by optimizing the crystal structures (experimental and submitted) of the first three CCDC blind tests held in 1999, 2001, and 2004. Except for the crystal structures of compound IX, all structures were reminimized and ranked according to their lattice energies. The hybrid method computes the lattice energy of a crystal structure as the sum of the DFT total energy and a van der Waals (dispersion) energy correction. Considering all four blind tests, the crystal structure with the lowest lattice energy corresponds to the experimentally observed structure for 12 out of 14 molecules. Moreover, good geometrical agreement is observed between the structures determined by the hybrid method and those measured experimentally. In comparison with the correct submissions made by the blind test participants, all hybrid optimized crystal structures (apart from compound II) have the smallest calculated root mean squared deviations from the experimentally observed structures. It is predicted that a new polymorph of compound V exists under pressure. PMID:19950907
Potential energy hypersurface and molecular flexibility
NASA Astrophysics Data System (ADS)
Koča, Jaroslav
1993-02-01
The molecular flexibility phenomenon is discussed from the conformational potential energy(hyper) surface (PES) point of view. Flexibility is considered as a product of three terms: thermodynamic, kinetic and geometrical. Several expressions characterizing absolute and relative molecular flexibility are introduced, depending on a subspace studied of the entire conformational space, energy level E of PES as well as absolute temperature. Results obtained by programs DAISY, CICADA and PANIC in conjunction with molecular mechanics program MMX for flexibility analysis of isopentane, 2,2-dimethylpentane and isohexane molecules are introduced.
Oyeyemi, Victor B.; Pavone, Michele; Carter, Emily A.
2011-11-03
Quantum chemistry has become one of the most reliable tools for characterizing the thermochemical underpinnings of reactions, such as bond dissociation energies (BDEs). The accurate prediction of these particular properties (BDEs) are challenging for ab initio methods based on perturbative corrections or coupled cluster expansions of the single-determinant Hartree-Fock wave function: the processes of bond breaking and forming are inherently multi-configurational and require an accurate description of non-dynamical electron correlation. To this end, we present a systematic ab initio approach for computing BDEs that is based on three components: (1) multi-reference single and double excitation configuration interaction (MRSDCI) for the electronic energies; (2) a two-parameter scheme for extrapolating MRSDCI energies to the complete basis set limit; and (3) DFT-B3LYP calculations of minimumenergy structures and vibrational frequencies to account for zero point energy and thermal corrections. We validated our methodology against a set of reliable experimental BDE values of C*C and C*H bonds of hydrocarbons. The goal of chemical accuracy is achieved, on average, without applying any empirical corrections to the MRSDCI electronic energies. We then use this composite scheme to make predictions of BDEs in a large number of hydrocarbon molecules for which there are no experimental data, so as to provide needed thermochemical estimates for fuel molecules.
Computed potential energy surfaces for chemical reactions
NASA Technical Reports Server (NTRS)
Walch, Stephen P.
1994-01-01
Quantum mechanical methods have been used to compute potential energy surfaces for chemical reactions. The reactions studied were among those believed to be important to the NASP and HSR programs and included the recombination of two H atoms with several different third bodies; the reactions in the thermal Zeldovich mechanism; the reactions of H atom with O2, N2, and NO; reactions involved in the thermal De-NO(x) process; and the reaction of CH(squared Pi) with N2 (leading to 'prompt NO'). These potential energy surfaces have been used to compute reaction rate constants and rates of unimolecular decomposition. An additional application was the calculation of transport properties of gases using a semiclassical approximation (and in the case of interactions involving hydrogen inclusion of quantum mechanical effects).
NASA Astrophysics Data System (ADS)
Balabin, Roman M.; Lomakina, Ekaterina I.
2009-08-01
Artificial neural network (ANN) approach has been applied to estimate the density functional theory (DFT) energy with large basis set using lower-level energy values and molecular descriptors. A total of 208 different molecules were used for the ANN training, cross validation, and testing by applying BLYP, B3LYP, and BMK density functionals. Hartree-Fock results were reported for comparison. Furthermore, constitutional molecular descriptor (CD) and quantum-chemical molecular descriptor (QD) were used for building the calibration model. The neural network structure optimization, leading to four to five hidden neurons, was also carried out. The usage of several low-level energy values was found to greatly reduce the prediction error. An expected error, mean absolute deviation, for ANN approximation to DFT energies was 0.6±0.2 kcal mol-1. In addition, the comparison of the different density functionals with the basis sets and the comparison of multiple linear regression results were also provided. The CDs were found to overcome limitation of the QD. Furthermore, the effective ANN model for DFT/6-311G(3df,3pd) and DFT/6-311G(2df,2pd) energy estimation was developed, and the benchmark results were provided.
Global Expression for Representing Diatomic Potential-Energy Curves
NASA Technical Reports Server (NTRS)
Ferrante, John; Schlosser, Herbert; Smith, John R.
1991-01-01
A three-parameter expression that gives an accurate fit to diatomic potential curves over the entire range of separation for charge transfers between 0 and 1. It is based on a generalization of the universal binding-energy relation of Smith et al. (1989) with a modification that describes the crossover from a partially ionic state to the neutral state at large separations. The expression is tested by comparison with first-principles calculations of the potential curves ranging from covalently bonded to ionically bonded. The expression is also used to calculate spectroscopic constants form a curve fit to the first-principles curves. A comparison is made with experimental values of the spectroscopic constants.
Potential Energy Curves of Hydrogen Fluoride
NASA Technical Reports Server (NTRS)
Fallon, Robert J.; Vanderslice, Joseph T.; Mason, Edward A.
1960-01-01
Potential energy curves for the X(sup 1)sigma+ and V(sup 1)sigma+ states of HF and DF have been calculated by the Rydberg-Klein-Rees method. The results calculated from the different sets of data for HF and DF are found to be in very good agreement. The theoretical results of Karo are compared to the experimental results obtained here.
Energy Savings Potential of Radiative Cooling Technologies
Fernandez, Nicholas; Wang, Weimin; Alvine, Kyle J.; Katipamula, Srinivas
2015-11-30
Pacific Northwest National Laboratory (PNNL), with funding from the U.S. Department of Energy’s (DOE’s) Building Technologies Program (BTP), conducted a study to estimate, through simulation, the potential cooling energy savings that could be achieved through novel approaches to capturing free radiative cooling in buildings, particularly photonic ‘selective emittance’ materials. This report documents the results of that study.
The Wind Energy Potential of Iceland
NASA Astrophysics Data System (ADS)
Nawri, Nikolai; Nína Petersen, Guðrún; Bjornsson, Halldór; Hahmann, Andrea N.; Jónasson, Kristján; Bay Hasager, Charlotte; Clausen, Niels-Erik
2014-05-01
While Iceland has an abundant wind energy resource, its use for electrical power production has so far been limited. Electricity in Iceland is generated primarily from hydro- and geothermal sources, and adding wind energy has so far not been considered practical or even necessary. However, wind energy is becoming a more viable option, as opportunities for new hydro- or geothermal power installations become limited. In order to obtain an estimate of the wind energy potential of Iceland, a wind atlas has been developed as part of the joint Nordic project 'Improved Forecast of Wind, Waves and Icing' (IceWind). Downscaling simulations performed with the Weather Research and Forecasting (WRF) model were used to determine the large-scale wind energy potential of Iceland. Local wind speed distributions are represented by Weibull statistics. The shape parameter across Iceland varies between 1.2 and 3.6, with the lowest values indicative of near-exponential distributions at sheltered locations, and the highest values indicative of normal distributions at exposed locations in winter. Compared with summer, average power density in winter is increased throughout Iceland by a factor of 2.0 - 5.5. In any season, there are also considerable spatial differences in average wind power density. Relative to the average value within 10 km of the coast, power density across Iceland varies between 50 - 250%, excluding glaciers, or between 300 - 1500 W m-2 at 50 m above ground level in winter. At intermediate elevations of 500 - 1000 m above mean sea level, power density is independent of the distance to the coast. In addition to seasonal and spatial variability, differences in average wind speed and power density also exist for different wind directions. Along the coast in winter, power density of onshore winds is higher by 100 - 700 W m-2 than that of offshore winds. The regions with the highest average wind speeds are impractical for wind farms, due to the distances from road
Potential energy surfaces of Polonium isotopes
NASA Astrophysics Data System (ADS)
Nerlo-Pomorska, B.; Pomorski, K.; Schmitt, C.; Bartel, J.
2015-11-01
The evolution of the potential energy landscape is analysed in detail for ten even-even polonium isotopes in the mass range 188\\lt A\\lt 220 as obtained within the macroscopic-microscopic approach, relying on the Lublin-Strasbourg drop model and the Yukawa-folded single-particle energies for calculating the microscopic shell and pairing corrections. A variant of the modified Funny-Hills nuclear shape parametrization is used to efficiently map possible fission paths. The approach explains the main features of the fragment partition as measured in low-energy fission along the polonium chain. The latter lies in a transitional region of the nuclear chart, and will be essential to consistently understand the evolution of fission properties from neutron-deficient mercury to heavy actinides. The ability of our method to predict fission observables over such an extended region looks promising.
NASA Astrophysics Data System (ADS)
Sellers, Michael; Lisal, Martin; Brennan, John
2015-06-01
Investigating the ability of a molecular model to accurately represent a real material is crucial to model development and use. When the model simulates materials in extreme conditions, one such property worth evaluating is the phase transition point. However, phase transitions are often overlooked or approximated because of difficulty or inaccuracy when simulating them. Techniques such as super-heating or super-squeezing a material to induce a phase change suffer from inherent timescale limitations leading to ``over-driving,'' and dual-phase simulations require many long-time runs to seek out what frequently results in an inexact location of phase-coexistence. We present a compilation of methods for the determination of solid-solid and solid-liquid phase transition points through the accurate calculation of the chemical potential. The methods are applied to the Smith-Bharadwaj atomistic potential's representation of cyclotrimethylene trinitramine (RDX) to accurately determine its melting point (Tm) and the alpha to gamma solid phase transition pressure. We also determine Tm for a coarse-grain model of RDX, and compare its value to experiment and atomistic counterpart. All methods are employed via the LAMMPS simulator, resulting in 60-70 simulations that total 30-50 ns. Approved for public release. Distribution is unlimited.
Communication: Accurate hydration free energies at a wide range of temperatures from 3D-RISM.
Misin, Maksim; Fedorov, Maxim V; Palmer, David S
2015-03-01
We present a new model for computing hydration free energies by 3D reference interaction site model (3D-RISM) that uses an appropriate initial state of the system (as suggested by Sergiievskyi et al.). The new adjustment to 3D-RISM theory significantly improves hydration free energy predictions for various classes of organic molecules at both ambient and non-ambient temperatures. An extensive benchmarking against experimental data shows that the accuracy of the model is comparable to (much more computationally expensive) molecular dynamics simulations. The calculations can be readily performed with a standard 3D-RISM algorithm. In our work, we used an open source package AmberTools; a script to automate the whole procedure is available on the web (https://github.com/MTS-Strathclyde/ISc). PMID:25747054
NASA Astrophysics Data System (ADS)
Yoshidome, Takashi; Ekimoto, Toru; Matubayasi, Nobuyuki; Harano, Yuichi; Kinoshita, Masahiro; Ikeguchi, Mitsunori
2015-05-01
The hydration free energy (HFE) is a crucially important physical quantity to discuss various chemical processes in aqueous solutions. Although an explicit-solvent computation with molecular dynamics (MD) simulations is a preferable treatment of the HFE, huge computational load has been inevitable for large, complex solutes like proteins. In the present paper, we propose an efficient computation method for the HFE. In our method, the HFE is computed as a sum of
NASA Technical Reports Server (NTRS)
Stallcop, James R.; Partridge, Harry; Levin, Eugene; Arnold, Jim (Technical Monitor)
2001-01-01
The potential energy surfaces for H2-N and N2-N interactions are calculated by accurate ab initio methods and applied to determine transport data. The results confirm that an effective potential energy for accurately determining transport properties can be calculated using a single orientation. A simple method is developed to determine the dispersion coefficients of effective potential energies Effective potential energies required for O2-O collisions are determ=ined. The H2-N, N2-N, O2-H, and O2-O collision integrals are calculated and tabulated for a large range of temperatures. The theoretical values of the N2-N and O2-O diffusion coefficients compare well with measured data available at room temperature.
NASA Astrophysics Data System (ADS)
Sun, Yuansheng; Periasamy, Ammasi
2010-03-01
Förster resonance energy transfer (FRET) microscopy is commonly used to monitor protein interactions with filter-based imaging systems, which require spectral bleedthrough (or cross talk) correction to accurately measure energy transfer efficiency (E). The double-label (donor+acceptor) specimen is excited with the donor wavelength, the acceptor emission provided the uncorrected FRET signal and the donor emission (the donor channel) represents the quenched donor (qD), the basis for the E calculation. Our results indicate this is not the most accurate determination of the quenched donor signal as it fails to consider the donor spectral bleedthrough (DSBT) signals in the qD for the E calculation, which our new model addresses, leading to a more accurate E result. This refinement improves E comparisons made with lifetime and spectral FRET imaging microscopy as shown here using several genetic (FRET standard) constructs, where cerulean and venus fluorescent proteins are tethered by different amino acid linkers.
Accurate hydrogen bond energies within the density functional tight binding method.
Domínguez, A; Niehaus, T A; Frauenheim, T
2015-04-01
The density-functional-based tight-binding (DFTB) approach has been recently extended by incorporating one-center exchange-like terms in the expansion of the multicenter integrals. This goes beyond the Mulliken approximation and leads to a scheme which treats in a self-consistent way the fluctuations of the whole dual density matrix and not only its diagonal elements (Mulliken charges). To date, only the performance of this new formalism to reproduce excited-state properties has been assessed (Domínguez et al. J. Chem. Theory Comput., 2013, 9, 4901-4914). Here we study the effect of our corrections on the computation of hydrogen bond energies for water clusters and water-containing systems. The limitations of traditional DFTB to reproduce hydrogen bonds has been acknowledged often. We compare our results for a set of 22 small water clusters and water-containing systems as well as for five water hexadecamers to those obtained with the DFTB3 method. Additionally, we combine our extension with a third-order energy expansion in the charge fluctuations. Our results show that the new formalisms significantly improve upon original DFTB. PMID:25763597
Accurate human tissue characterization for energy-efficient wireless on-body communications.
Vallejo, Mónica; Recas, Joaquín; del Valle, Pablo García; Ayala, José L
2013-01-01
The demand for Wireless Body Sensor Networks (WBSNs) is rapidly increasing due to the revolution in wearable systems demonstrated by the penetration of on-the-body sensors in hospitals, sports medicine and general health-care practices. In WBSN, the body acts as a communication channel for the propagation of electromagnetic (EM) waves, where losses are mainly due to absorption of power in the tissue. This paper shows the effects of the dielectric properties of biological tissues in the signal strength and, for the first time, relates these effects with the human body composition. After a careful analysis of results, this work proposes a reactive algorithm for power transmission to alleviate the effect of body movement and body type. This policy achieves up to 40.8% energy savings in a realistic scenario with no performance overhead. PMID:23752565