A fully ab initio potential energy surface for C1H_2 reactive system

An ab initio analytical potential energy surface called BW3 for the CIH2 reactive system is presented. The fit of this surface is based on about 1 200 ab initio energy points, computed with multi-reference configuration interaction(MRCI) and scaling external correlation (SEC) method and a very large basis set. The precision in the fit is very high. The BW3 surface could reproduce correctly the dissociation energy of H2 and HCI, and the endothermicity of the Cl H2 abstraction reaction. For the Cl H2 abstraction reaction, the saddle point of BW3 lies in collinear geometries, and the barrier height is 32.84 kJ/mol; for the H CIH exchange reaction, the barrier of BW3 is also linear, with a height of 77.40 kJ/mol.     

A fully ab initio potential energy surface for ClH2 reactive system

Anab initio analytical potential energy surface called BW3 for the ClH₂ reactive system is presented. The fit of this surface is based on about 1 200ab initio energy points, computed with multi-reference configuration interaction(MRCI) and scaling external correlation (SEC) method and a very large basis set. The precision in the fit is very high. The BW3 surface could reproduce correctly the dissociation energy of H₂ and HCl, and the endothermicity of the Cl + H₂ abstraction reaction. For the Cl + H₂ abstraction reaction, the saddle point of BW3 lies in collinear geometries, and the barrier height is 32.84 kJ/mol; for the H + ClH exchange reaction, the barrier of BW3 is also linear, with a height of 77.40 kJ/mol.     

General methodology in two dimensions for classical simulation of reactive and nonreactive events on ab initio potential energy surfaces

A completely general two-dimensional (2D) methodology for the classical simulation of reactive and nonreactive events on ab initio potential energy surfaces is introduced and tested. The methodology requires the minimum amount of information given a priori—geometries and energies at these geometries. From a list of ab initio geometries and energies, simulations may be executed and a distribution of outcomes obtained. The method introduced attempts a local approach at simulating the dynamics of the system, rather than a global analytic fit to the potential energy surface. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1431–1444, 1998     

A fullyab initio potential energy surface for ClH2 reactive system

Anab initio analytical potential energy surface called BW3 for the ClH₂ reactive system is presented. The fit of this surface is based on about 1 200ab initio energy points, computed with multi-reference configuration interaction(MRCI) and scaling external correlation (SEC) method and a very large basis set. The precision in the fit is very high. The BW3 surface could reproduce correctly the dissociation energy of H₂ and HCl, and the endothermicity of the Cl + H₂ abstraction reaction. For the Cl + H₂ abstraction reaction, the saddle point of BW3 lies in collinear geometries, and the barrier height is 32.84 kJ/mol; for the H + ClH exchange reaction, the barrier of BW3 is also linear, with a height of 77.40 kJ/mol.     

An ab initio potential energy surface and vibrational energy levels of ZnH2

A three‐dimensional potential energy surface of the electronic ground state of ZnH₂ (${X}^1\sum _g^ +$ ) molecule is constructed from more than 7500 ab initio points calculated at the internally contracted multireference configuration interaction with the Davidson correction (icMRCI+Q) level employing large basis sets. The calculated relative energies of various dissociation reactions are in good agreement with the previous theoretical/experimental values. Low‐lying vibrational energy levels of ZnH₂, ZnD₂, and HZnD are calculated on the three‐dimensional potential energy surface using the Lanczos algorithm, and found to be in good agreement with the available experimental band origins and the previous theoretical values. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Efficient evaluation of the accuracy of molecular quantum dynamics on an approximate analytical or interpolated ab initio potential energy surface

Ab initio electronic structure methods have reached a satisfactory accuracy for the calculation of static properties, but remain too expensive for quantum dynamical calculations. Recently, an efficient semiclassical method was proposed to evaluate the accuracy of quantum dynamics on an approximate potential without having to perform the expensive quantum dynamics on the accurate potential. Here, this method is applied for the first time to evaluate the accuracy of quantum dynamics on an approximate analytical or interpolated potential in comparison to the quantum dynamics on an accurate potential obtained by an ab initio electronic structure method. Specifically, the vibrational dynamics of H₂ on a Morse potential is compared with that on the full CI potential, and the photodissociation dynamics of CO₂ on a LEPS potential with that on the excited ¹Π surface computed at the EOM‐CCSD/aug‐cc‐pVDZ level of theory. Finally, the effect of discretization of a potential energy surface on the quantum dynamics is evaluated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2426–2435, 2010     

An accurate three-dimensional potential energy surface for the He-Na_2 complex

An accurate three-dimensional potential energy surface(PES) for the He-Na2 van der Waals comple was calculated at the coupled cluster singles-and-doubles with noniterative inclusion of connecte triple(CCSD(T)) level of theory.A mixed basis set,aug-cc-pVQZ for the He atom and cc-pCVQZ for th sodium atom,and an additional(3s3p2d1f) set of midbond functions were used.The computed inte action energies in 819 configurations were fitted to a 96-parameter analytic potential model by leas squares fitting.The PES has two shallow wells corresponding to the T-shaped structure and the linea configuration,which are located at 12.5a0 and 14 a0 with depths of 1.769 and 1.684 cm-1,respectivel The whole potential energy surface exhibits weak anisotropy.Based on the fitted PES,state-to-stat differential cross sections were calculated.     

Complete basis set ab initio and hybrid density functional theory exploration of the potential energy surface in the reaction between an amino radical and nitrogen oxide

The potential energy surface for the reaction involving NH₂ plus NO was explored with a quadratic complete basis set ab initio approach and three hybrid density functional theory methods, the target being to accurately estimate activation barriers and the relative stability of the nitrogen–oxygen isomers. The computational results were compared with previously performed ab initio calculations and new, more accredited values for the NH₂NO rearrangement to HNNOH and for the HNNOH decomposition reaction were suggested. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 409–414, 1998     

An ab initio MO study on the direct and migratory reaction mechanism for the formation of HF in H+ClF system

ＰｏｌａｎｙｉａｎｄｃｏｗｏｒｋｅｒｓｈａｖｅｓｔｕｄｉｅｄａｓｅｒｉｅｓｏｆｒｅａｃｔｉｏｎｓｂｅｔｗｅｅｎＨａｔｏｍｓａｎｄｉｎｔｅｒｈａｌｏｇｅｎｓＸＹ(Ｙｉｓｔｈｅｍｏｒｅｅｌｅｃｔｒｏｎｅｇａｔｉｖｅｈａｌｏｇｅｎａｔｏｍ)ｕｓｉｎｇｔｈｅｉｎｆｒａｒｅｄｃｈｅｍｉｌｕｍｉｎｅｓｃｅｎｃｅｍｅｔｈｏｄ[1,2]ａｎｄｈａｖｅｏｂｓｅｒｖｅｄｂｉｍｏｄａｌｅｎｅｒｇｙｄｉｓｔｒｉｂｕｔｉｏｎｓｆｏｒｔｈｅＨＹｐｒｏｄｕｃｔ;ｔｈａｔｉｓｔｏｓａｙ,ｔｈｅｔｏｔａｌａｖａｉｌａｂｌｅｅｎｅｒｇｙｆ…     

Accurate ab initio potential energy surface and vibration‐rotation energy levels of hydrogen peroxide

The accurate ground‐state potential energy surface of hydrogen peroxide, H₂O₂, has been determined from ab initio calculations using the coupled‐cluster approach in conjunction with the correlation‐consistent basis sets up to septuple‐zeta quality. Results obtained with the conventional and explicitly correlated coupled‐cluster methods were compared. The core–electron correlation, scalar relativistic, and higher‐order valence–electron correlation effects were taken into account. The adiabatic effects were also discussed. The vibration–rotation energy levels of the H₂O₂, D₂O₂, and HOOD isotopologues were predicted, and the experimental vibrational fundamental wavenumbers were reproduced to 1 cm^?1 (“spectroscopic”) accuracy. © 2012 Wiley Periodicals, Inc.     

F2 dimer: Improved intermolecular potential energy surface using ab initio calculations

A computational study on the intermolecular potential energy of 44 different orientations of F₂ dimers is presented. Basis set superposition error (BSSE) corrected potential energy surface is calculated using the supermolecular approach at CCSD(T) and QCISD(T) levels of theory. The interaction energies obtained using the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets are extrapolated to the complete basis set limit using the latest extrapolation scheme. The basis set effect is checked and it is found that the extrapolated intermolecular energies provide the best compromise between the accuracy and computational cost. Among 1320 energy points of F₂–F₂ system covering more relative orientations, the most stable structure of the dimers was obtained with a well depth of ?146.62 cm^?1 that related to cross configuration, and the most unstable structure is related to linear orientation with a well depth of ?52.63 cm^?1. The calculated second virial coefficients are in good agreement with experimental data. The latest extrapolation scheme of the complete basis set limit at the CCSD(T) level of theory is used to determine the intermolecular potential energy surface of the F₂ dimer. Comparing the results obtained by the latest scheme with those by older schemes show that the new approach provides the best compromise between accuracy and computational cost.     

An ab initio Study of Intermolecular Potential for Ne-HBr Complex

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｒａｒｅ ｇａｓｈｙｄｒｏｇｅｎｈａｌｉｄｅ (Ｍ—ＨＸ)ｃｏｍｐｌｅｘｅｓｈａｖｅｌｏｎｇｂｅｅｎｓｔｕｄｉｅｄｆｏｒｉｎｓｉｇｈｔｓｉｎｔｏｔｈｅｎａｔｕｒｅｏｆｉｎ ｔｅｒｍｏｌｅｃｕｌａｒｆｏｒｃｅｓａｎｄｄｅｔａｉｌｓｏｆｉｎｔｅｒ ａｎｄｉｎｔｒａｍｏｌｅｃｕｌａｒｄｙｎａｍｉｃｓ .1,2 ＳｙｓｔｅｍａｔｉｃｉｎｖｅｓｔｉｇａｔｉｏｎｓｏｆｔｈｅｃｏｍｐｌｅｘｅｓｏｆＸｅ ,Ｋｒ,Ａｒ,ａｎｄｌａｔｅｌｙＮｅｗｉｔｈＨＦ ,ＨＣｌ,ＨＢｒａｎｄｔｈｅｉｒｄｅｕｔｅｒａ…     

DNA base trimers: empirical and quantum chemical ab initio calculations versus experiment in vacuo

A complete scan of the potential‐energy surfaces for selected DNA base trimers has been performed by a molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The resulting most stable/populated structures were then reoptimized at a correlated ab initio level by employing resolution of the identity, Møller–Plesset second‐order perturbation theory (RI‐MP2). A systematic study of these trimers at such a complete level of electronic structure theory is presented for the first time. We show that prior experimental and theoretical interpretations were incorrect in assuming that the most stable structures of the methylated trimers corresponded to planar systems characterized by cyclic intermolecular hydrogen bonding. We found that stacked structures of two bases with the third base in a T‐shape arrangement are the global minima in all of the methylated systems: they are more stable than the cyclic planar structures by about 10 kcal mol^?1. The different behaviors of nonmethylated and methylated trimers is also discussed. The high‐level geometries and interaction energies computed for the trimers serve also as a reference for the testing of recently developed density functional theory (DFT) functionals with respect to their ability to correctly describe the balance between the electrostatic and dispersion contributions that bind these trimers together. The recently reported M052X functional with a polarized triple‐zeta basis set predicts 11 uracil trimer interaction energies with a root‐mean‐square error of 2.3 kcal mol^?1 relative to highly correlated ab initio theoretical calculations.     

Ab initio quantum chemical studies of reaction mechanism for CN with CH2CO

Six product channels have been found in the association reaction of CN + CH₂CO, and a variety of possible complexes and saddle points along the minimum energy reaction paths have been characterized at the UMP2(full)/6‐31G(d) level. The dominant reaction channels are the production of CH₂CN + CO and CH₂NC + CO. The isomerization and dissociation reactions of the major products of CH₂CN and CH₂NC have been investigated using the G2MP2 level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Ab initio intermolecular potential energy surface of He-LiH

The intermolecular potential energy surface of He-LiH complex was studied using the full-electronic complete forth-order Miller-Plesset perturbation (MPPT) method.In ab initio calculations,the bond length of LiH was fixed at 0 159 5 nm.The potential has two local minima of Vm=-179.93 cm for the linear He LiH geormetrv at Rm=0.227 nm and Vm=-10.44 cm-1 for the linear He-HL1 geometry at Rm=0.516 nm The potemal exhibits strong anisotropy The analytic potential function with 31 parameters was determined by fitting to the calculated ab,mtio potentials The influence of variation of LiH bond length on the potential energy surface was also studied     

Ab initio MO study of potential energy surface of NH2 with CN reaction

An extensive quantum chemical study of the potential energy surfaces (PES) for the association reaction of NH₂ with CN and the subsequent isomerization and dissociation reactions has been carried out using density functional theory (DFT)/B3LYP/6‐311++G(3df,2p) level of theory on both singlet and triplet states. The reaction mechanism on the triplet surface is more complicated than that on the singlet surface. A total of 19 isomers and 46 transition states have been identified and characterized on the triplet PES. Among them, IM2 (IM2a), IM3 (IM3a, IM3b), and IM10 are the lowest‐lying isomers with thermodynamic stability. Twenty available dissociation channels, depending on the different initial isomers, have been identified. On the singlet surface, only 12 isomers and 16 transition states have been found, and among them IM1(S) and IM2(S) are the lowest‐lying isomers. The higher isomerization and dissociation barriers on the singlet surface indicate that the addition and the subsequent reactions of NH₂+CN are most likely to occur on the triplet PES because of the lower barriers. A prediction can be made for the possible mechanism explaining the production of H+HNCN. Besides HNCN, other major products are NH+HCN and NH+HNC, which are produced by direct dissociation reactions from triplet IM2 and IM3, respectively. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Singlet-Triplet Excited-State Inversion in Heptazine and Related Molecules: Assessment of TD-DFT and ab initio Methods

We have investigated the origin of the S₁-T₁ energy levels inversion for heptazine, and other N-doped π-conjugated hydrocarbons, leading thus to an unusually negative singlet-triplet energy gap ( ). Since this inversion might rely on substantial doubly-excited configurations to the S₁ and/or T₁ wavefunctions, we have systematically applied multi-configurational SA-CASSCF and SC-NEVPT2 methods, SCS-corrected CC2 and ADC(2) approaches, and linear-response TD-DFT, to analyze if the latter method could also face this challenging issue. We have also extended the study to B-doped π-conjugated systems, to see the effect of chemical composition on the results. For all the systems studied, an intricate interplay between the singlet-triplet exchange interaction, the influence of doubly-excited configurations, and the impact of dynamic correlation effects, serves to explain the values found for most of the compounds, which is not predicted by TD-DFT.