An ab initio potential energy surface and vibrational states of MgH2(1(1)A')

A three-dimensional global potential energy surface for the ground electronic state of MgH(2) is constructed from more than 3000 ab initio points calculated using the internally contracted multireference configuration interaction method with the Davidson correction at the complete basis set limit. Low-lying vibrational energy levels of MgH(2) and MgD(2) are calculated using the Lanczos algorithm, and found to be in good agreement with known experimental band origins. The majority of the vibrational energy levels up to 8000 cm(-1) are assigned with normal mode quantum numbers. However, our results indicate a gradual transition from a normal mode regime for the stretching vibrations at low energies to a local mode regime near 7400 cm(-1), as evidenced by a decreasing energy gap between the (n(1),0,0) and (n(1)-1,0,1) vibrational states and bifurcation of the corresponding wave functions.     

Accurate ab initio based multisheeted double many-body expansion potential energy surface for the three lowest electronic singlet states of H3+

The authors present diabatic and adiabatic potential energy surfaces for the three lowest electronic singlet states of H3+. The modeling of the surfaces is based on the multi-sheeted double many-body expansion method which consists of dressing the various matrix elements of the diatomics-in-molecules potential matrix with three-body terms. The avoided crossing between the two lowest states and the conical intersection between the second and the third state are accurately represented by construction.     

Ab initio calculation and quasi‐classical dynamics study of the two lowest potential energy surfaces of the O(1D)+HBr system*

The lowest singlet 1¹A′ and 1¹A″ potential energy surfaces (PES) of the O(¹D)+HBr system have been ab initio computed. The complete active space self‐consistent field (CASSCF) method was used in most of the calculations, considering all the valence orbitals as active. The calculations were complemented with both analytical gradient calculations to characterize the stationary points and multireference configuration interaction (MRCI) calculations at selected nuclear geometries to improve the determination of the barrier heights and of the energetics. Electronic energy values for both PESs were then independently fitted by polynomial expansions in bond order coordinates. On the fitted surfaces quasi‐classical trajectories were separately run. Single‐surface calculations behave qualitatively different for the ground and the excited PES at low collision energies. A satisfactory agreement with existing experimental data was obtained by using the ground PES while calculations performed on the excited 1¹A″ PES worsened the agreement. However, when collision energy is increased, detailed experimental distributions are less well reproduced by calculations on the ground PES. This may imply the participation via nonadiabatic transitions of the 2¹A′ PES at higher energies while the adiabatic ground singlet PES well describes reactive scattering at low collision energy. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Ab initio adiabatic and quasidiabatic potential energy surfaces of lowest four electronic states of the H+ + O2 system

Ab initio global adiabatic and quasidiabatic potential energy surfaces of lowest four electronic (1-4 (3)A(")) states of the H(+)+O(2) system have been computed in the Jacobi coordinates (R,r,γ) using Dunning's cc-pVTZ basis set at the internally contracted multireference (single and double) configuration interaction level of accuracy, which are relevant to the dynamics studies of inelastic vibrational and charge transfer processes observed in the scattering experiments. The computed equilibrium geometry parameters of the bound [HO(2)](+) ion in the ground electronic state and other parameters for the transition state for the isomerization process, HOO(+)?OOH(+) are in good quantitative agreement with those available from the high level ab initio calculations, thus lending credence to the accuracy of the potential energy surfaces. The nonadiabatic couplings between the electronic states have been analyzed in both the adiabatic and quasidiabatic frameworks by computing the nonadiabatic coupling matrix elements and the coupling potentials, respectively. It is inferred that the dynamics of energy transfer processes in the scattering experiments carried out in the range of 9.5-23 eV would involve all the four electronic states.     

High-dimensional ab initio potential energy surfaces for reaction dynamics calculations

There has been great progress in the development of potential energy surfaces (PESs) for reaction dynamics that are fits to ab initio energies. The fitting techniques described here explicitly represent the invariance of the PES with respect to all permutations of like atoms. A review of a subset of dynamics calculations using such PESs (currently 16 such PESs exist) is then given. Bimolecular reactions of current interest to the community, namely, H + CH(4) and F + CH(4), are focused on. Unimolecular reactions are then reviewed, with a focus on the photodissociation dynamics of H(2)CO and CH(3)CHO, where so-called "roaming" pathways have been discovered. The challenges for electronically non-adiabatic reactions, and associated PESs, are presented with a focus on the OH* + H(2) reaction. Finally, some thoughts on future directions and challenges are given.     

Rate constant for OH(2 Pi)+O(3P)-->H(2S)+O2(3 Sigma g-) reaction on an improved ab initio potential energy surface and implications for the interstellar oxygen problem

The authors report a global potential energy surface for the ground electronic state of HO(2)(X (2)A(")), which improves upon the XXZLG potential [Xu and et al., J. Chem. Phys. 122, 244305 (2005)] with additional high-level ab initio points for the long-range interaction potential in the O+OH channel. Exact J=0 quantum mechanical reaction probabilities were calculated on the new potential and the rate constant for the title reaction was obtained using a J-shifting method. The calculated rate constant is in good agreement with available experimental values and our results predict a significantly lower rate at temperature range below 30 K, offering a possible explanation for the "interstellar oxygen problem."     

CVRQD ab initio ground-state adiabatic potential energy surfaces for the water molecule

The high accuracy ab initio adiabatic potential energy surfaces (PESs) of the ground electronic state of the water molecule, determined originally by Polyansky et al. [Science 299, 539 (2003)] and called CVRQD, are extended and carefully characterized and analyzed. The CVRQD potential energy surfaces are obtained from extrapolation to the complete basis set of nearly full configuration interaction valence-only electronic structure computations, augmented by core, relativistic, quantum electrodynamics, and diagonal Born-Oppenheimer corrections. We also report ab initio calculations of several quantities characterizing the CVRQD PESs, including equilibrium and vibrationally averaged (0 K) structures, harmonic and anharmonic force fields, harmonic vibrational frequencies, vibrational fundamentals, and zero-point energies. They can be considered as the best ab initio estimates of these quantities available today. Results of first-principles computations on the rovibrational energy levels of several isotopologues of the water molecule are also presented, based on the CVRQD PESs and the use of variational nuclear motion calculations employing an exact kinetic energy operator given in orthogonal internal coordinates. The variational nuclear motion calculations also include a simplified treatment of nonadiabatic effects. This sophisticated procedure to compute rovibrational energy levels reproduces all the known rovibrational levels of the water isotopologues considered, H(2) (16)O, H(2) (17)O, H(2) (18)O, and D(2) (16)O, to better than 1 cm(-1) on average. Finally, prospects for further improvement of the ground-state adiabatic ab initio PESs of water are discussed.     

A quasiclassical trajectory study for the N(4S)+O2(X 3Sigmag-)-->NO(X 2Pi)+O(3P) reaction on the new 2A' and 4A' potential-energy surfaces

A quasiclassical trajectory study with the sixth-order explicit symplectic algorithm of the N(4S)+O2(X 3Sigmag-)-->NO(X 2Pi)+O(3P) atmospheric reaction has been performed by employing the new 2A' and 4A' potential-energy surfaces reported by Sayos et al. [J. Chem. Phys. 117, 670 (2002)]. For the translational temperature considered up to 10,000 K, the larger relative translational energy and the higher rovibrational levels of O2 molecule with respect to the previous works have been taken into account, and a clearer database about the character of the total reaction cross section has been presented in this work. The dependence of microscopic rate constants on the vibrational level of O2 molecule at T=3000, 5000, and 10,000 K has been exhibited, and we can see that the values of log10 k(T,v,J) vary almost linearly with the vibrational level of O2 molecule. The thermal rate constants at the translational temperature between 300 and 10,000 K have been evaluated and compared with the experimental and previous theoretical results. It is found that the thermal rate constants determined in this work have a better agreement with the experimental data and can provide a more valid theoretical reference at the translational temperature considered for the title reaction.     

Ab initio potential energy surfaces, bound states, and electronic spectrum of the Ar-SH complex

New ab initio potential energy surfaces for the (2)Pi ground electronic state of the Ar-SH complex are presented, calculated at the RCCSD(T)/aug-cc-pV5Z level. Weakly bound rotation-vibration levels are calculated using coupled-channel methods that properly account for the coupling between the two electronic states. The resulting wave functions are analyzed and a new adiabatic approximation including spin-orbit coupling is proposed. The ground-state wave functions are combined with those obtained for the excited (2)Sigma(+) state [D. M. Hirst, R. J. Doyle, and S. R. Mackenzie, Phys. Chem. Chem. Phys. 6, 5463 (2004)] to produce transition dipole moments. Modeling the transition intensities as a combination of these dipole moments and calculated lifetime values [A. B. McCoy, J. Chem. Phys. 109, 170 (1998)] leads to a good representation of the experimental fluorescence excitation spectrum [M.-C. Yang, A. P. Salzberg, B.-C. Chang, C. C. Carter, and T. A. Miller, J. Chem. Phys. 98, 4301 (1993)].     

O(3P)+HBr(DBr)反应的含时量子散射计算

基于LEPS势能面, 用三维含时量子波包法对O(³P)+HBr(DBr)反应进行了准确的动力学计算. 计算的结果表明, 振动激发对这个反应是有效的, 而转动激发在某一能量范围内具有方位效应. 计算得到了该反应的速率常数和反应截面, 速率常数kO+HBr的计算值同实验值符合得很好. 通过对相应结果的对比, 可以发现这个反应具有比较明显的同位素效应.     

Quantum dynamical study of the O(1D)+HCl reaction employing three electronic state potential energy surfaces

Quantum dynamical calculations are reported for the title reaction, for both product arrangement channels and using potential energy surfaces corresponding to the three electronic states, 1 1A', 2 1A', and 1 1A", which correlate with both reactants and products. The calculations have been performed for J=0 using the time-dependent real wavepacket approach by Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. Reaction probabilities for both product arrangement channels on all three potential energy surfaces are presented for total energies between 0.1 and 1.1 eV. Product vibrational state distributions at two total energies, 0.522 and 0.722 eV, are also presented for both channels and all three electronic states. Product rotational quantum state distributions are presented for both product arrangement channels and all three electronic states for the first six product vibrational states.     

Near-infrared spectra and rovibrational dynamics on a four-dimensional ab initio potential energy surface of (HBr)2

Supersonic jet investigations of the (HBr)(2) dimer have been carried out using a tunable diode laser spectrometer to provide accurate data for comparison with results from a four-dimensional (4-D) ab initio potential energy surface (PES). The near-infrared nu(1) (+/-), nu(2) (+/-), and (nu(1)+nu(4))(-) bands of (H (79)Br)(2), (H (79)Br-H (81)Br), and (H (81)Br)(2) isotopomers have been recorded in the range 2500-2600 cm(-1) using a CW slit jet expansion with an upgraded near-infrared diode laser spectrometer. The 4-D PES has been calculated for (HBr)(2) using second-order M?ller-Plesset perturbation theory with an augmented and polarized 6-311G basis set. The potential is characterized by a global minimum occurring at the H bond structure with the distance between the center of masses (CM) of the monomer being R(CM)=4.10 A with angles theta(A)=10 degrees, theta(B)=100 degrees and a well depth of 692.2 cm(-1), theta(A) is the angle the HBr bond of monomer A makes with the vector from the CM of A to the CM of B, and theta(B) is the corresponding angle monomer B makes with the same CM-CM vector. The barrier for the H interchange occurs at the closed C(2h) structure for which R(CM)=4.07 A, theta(A)=45 degrees, theta(B)=135 degrees, and the barrier height is 73.9 cm(-1). The PES was fitted using a linear-least squares method and the rovibrational energy levels of the complex were calculated by a split pseudospectral method. The spectroscopic data provide accurate molecular parameters for the dimer that are then compared with the results predicted on the basis of the 4-D ab initio PES.     

HN2(2A') electronic manifold. I. A global ab initio study of first two states

A detailed ab initio multireference configuration interaction calculation with a standard aug-cc-pVTZ basis set is reported for the 1(2)A' and 2(2)A' states of the title system. The aim is to establish the dissociation scheme of all channels, while revealing the 2(2)A'/3(2)A' seam of conical intersections consistent with the crossings in the diatomic fragments. An ab initio mapping of linear NNH and T-shaped and linear NHN loci of conical intersections is also reported, jointly with a discussion of the topological features associated to a newly reported 2(2)A'/3(2)A' crossing seam.     

O(~3P)+HBr(DBr)反应的含时量子散射计算(英文)

基于LEPS势能面, 用三维含时量子波包法对O(3~P)+HBr(DBr)反应进行了准确的动力学计算. 计算的结果表明, 振动激发对这个反应是有效的, 而转动激发在某一能量范围内具有方位效应. 计算得到了该反应的速率常数和反应截面, 速率常数kO+HBr的计算值同实验值符合得很好. 通过对相应结果的对比, 可以发现这个反应具有比较明显的同位素效应.     

An ab initio study of the low-lying electronic states of S3

Accurate calculations of the low-lying singlet and triplet electronic states of thiozone, S(3), have been carried out using large multireference configuration interaction wave functions. Cuts of the full potential energy surfaces along the stretching and bending coordinates have been presented, together with the vertical excitation spectra. The strong experimentally observed absorption around 395 nm is assigned to the 1 (1)B(2) state, which correlates to ground state products. Absorption at wavelengths shorter than 260 nm is predicted to lead to singlet excited state products, S(2) (a (1)Delta(g))+S((1)D). The spectroscopic properties of the X (3)Sigma(g) (-), a (1)Delta(g), and b (1)Sigma(g) (+) electronic states of the S(2) radical have also been accurately characterized in this work. The investigations of the low-lying electronic states were accompanied by accurate ground state coupled cluster calculations of the thermochemistry of both S(2) and S(3) using large correlation consistent basis sets with corrections for core-valence correlation, scalar relativity, and atomic spin-orbit effects. Resulting values for D(0)(S(2)+S) and SigmaD(0) for S(3) are predicted to be 61.3 and 162.7 kcal/mol, respectively, with conservative uncertainties of +/-1 kcal/mol. Analogous calculations predict the C(2v)-D(3h) (open-cyclic) isomerization energy of S(3) to be 4.4+/-0.5 kcal/mol.     

Accurate ab initio determination of the adiabatic potential energy function and the Born-Oppenheimer breakdown corrections for the electronic ground state of LiH isotopologues

High level ab initio potential energy functions have been constructed for LiH in order to predict vibrational levels up to dissociation. After careful tests of the parameters of the calculation, the final adiabatic potential energy function has been composed from: (a) an ab initio nonrelativistic potential obtained at the multireference configuration interaction with singles and doubles level including a size-extensivity correction and quintuple-sextuple ζ extrapolations of the basis, (b) a mass-velocity-Darwin relativistic correction, and (c) a diagonal Born-Oppenheimer (BO) correction. Finally, nonadiabatic effects have also been considered by including a nonadiabatic correction to the kinetic energy operator of the nuclei. This correction is calculated from nonadiabatic matrix elements between the ground and excited electronic states. The calculated vibrational levels have been compared with those obtained from the experimental data [J. A. Coxon and C. S. Dickinson, J. Chem. Phys. 134, 9378 (2004)]. It was found that the calculated BO potential results in vibrational levels which have root mean square (rms) deviations of about 6-7 cm(-1) for LiH and ～3 cm(-1) for LiD. With all the above mentioned corrections accounted for, the rms deviation falls down to ～1 cm(-1). These results represent a drastic improvement over previous theoretical predictions of vibrational levels for all isotopologues of LiH.     

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.     

Ab initio potential energy surfaces, total absorption cross sections, and product quantum state distributions for the low-lying electronic states of N(2)O

Adiabatic potential energy surfaces for the six lowest singlet electronic states of N(2)O (X (1)A('), 2 (1)A('), 3 (1)A('), 1 (1)A("), 2 (1)A(") and 3 (1)A(")) have been computed using an ab initio multireference configuration interaction (MRCI) method and a large orbital basis set (aug-cc-pVQZ). The potential energy surfaces display several symmetry related and some nonsymmetry related conical intersections. Total photodissociation cross sections and product rotational state distributions have been calculated for the first ultraviolet absorption band of the system using the adiabatic ab initio potential energy and transition dipole moment surfaces corresponding to the lowest three excited electronic states. In the Franck-Condon region the potential energy curves corresponding to these three states lie very close in energy and they all contribute to the absorption cross section in the first ultraviolet band. The total angular momentum is treated correctly in both the initial and final states. The total photodissociation spectra and product rotational distributions are determined for N(2)O initially in its ground vibrational state (0,0,0) and in the vibrationally excited (0,1,0) (bending) state. The resulting total absorption spectra are in good quantitative agreement with the experimental results over the region of the first ultraviolet absorption band, from 150 to 220 nm. All of the lowest three electronically excited states [(1)Sigma(-)(1 (1)A(")), (1)Delta(2 (1)A(')), and (1)Delta(2 (1)A("))] have zero transition dipole moments from the ground state [(1)Sigma(+)(1 (1)A('))] in its equilibrium linear configuration. The absorption becomes possible only through the bending motion of the molecule. The (1)Delta(2 (1)A('))<--X (1)Sigma(+)((1)A(')) absorption dominates the absorption cross section with absorption to the other two electronic states contributing to the shape and diffuse structure of the band. It is suggested that absorption to the bound (1)Delta(2 (1)A(")) state makes an important contribution to the experimentally observed diffuse structure in the first ultraviolet absorption band. The predicted product rotational quantum state distribution at 203 nm agrees well with experimental observations.