Ab initio potential energy surface for Ar 3 +

The potential energy surface (PES) of linear Ar ₃ ⁺ is calculated at the MP4/6-31G* level including all single, double, triple and quadruple excitations. The results show that the PES of the linear Ar ₃ ⁺ has a very flat valley along the asymmetric stretching vibration normal mode, ν₃. A higher level quadratic configuration interaction calculation including single, double and triple substitutions QCISD (T) along this flat valley suggests that an asymmetric geometry energy minimum reported earlier based on MP2 [1] is due to symmetry breaking in UHF. The global minimum of the PES is found to be for the symmetric geometry atR _ab =R _bc =2.66±0.01 Å, which is in good agreement with the MRD-CI calculation [2] and expectations from our earlier photodissociation experiments [3]. The calculational results are compared with other theoretical calculations, and are discussed in the context of the photodissociation and dynamics of dissociation experiments conducted on Ar ₃ ⁺ .     

Ab initio potential energy surface and bound states of the Xe-CO complex

The first two-dimensional potential energy surface for the Xe-CO van der Waals interaction is calculated by the single and double excitation coupled-cluster theory with noniterative treatment of triple excitations. Mixed basis sets, aug-cc-pVQZ for the C and O atoms, and aug-cc-pVQZ-PP for the Xe atom, with an additional (3s3p2d2f1g) set of midbond functions, are used. Our potential energy surface has a single, nearly T-shaped minimum of -131.87 cm(-1) at R(e)=7.80a(0) and theta(e)=102.5 degrees. Based on the potential, the bound state energies are calculated for seven isotopomers of the Xe-(12)C(16)O complex, seven isotopomers of the Xe-(13)C(16)O complex, and three isotopomers of the Xe-(13)C(18)O complex. Compared with available experimental data, the predicted transition frequencies and spectroscopic constants are in good agreement with the experimental results.     

Ab initio potential energy surface and bound states for the Kr-OCS complex

The first ab initio potential energy surface of the Kr-OCS complex is developed using the coupled-cluster singles and doubles with noniterative inclusion of connected triples [CCSD(T)]. The mixed basis sets, aug-cc-pVTZ for the O, C, and S atom, and aug-cc-pVQZ-PP for the Kr atom, with an additional (3s3p2d1f) set of midbond functions are used. A potential model is represented by an analytical function whose parameters are fitted numerically to the single point energies computed at 228 configurations. The potential has a T-shaped global minimum and a local linear minimum. The global minimum occurs at R = 7.146 a(0), θ = 105.0° with energy of -270.73 cm(-1). Bound state energies up to J = 9 are calculated for three isotopomers (82)Kr-OCS, (84)Kr-OCS, and (86)Kr-OCS. Analysis of the vibrational wavefunctions and energies suggests the complex can exist in two isomeric forms: T-shaped and quasi-linear. The calculated transition frequencies and spectroscopic constants of the three isotopomers are in good agreement with the experimental values.     

Ab initio potential energy surfaces and nonadiabatic interactions in the H+ +NO collision system

Ab initio calculations on the H(+)+NO system have been carried out in Jacobi coordinates at the multireference configuration interaction level employing Dunning's correlation-consistent polarized valence triple zeta basis set to analyze the role of low-lying electronic excited states in influencing the collision dynamics relevant to the experimental collision energy range of 9.5-30 eV. The lowest two adiabatic potential energy surfaces, asymptotically correlating to H(+)+NO(X (2)Pi) and H((2)S)+NO(+)(X (1)Sigma(+)), have been obtained. Using ab initio procedures, the (radial) nonadiabatic couplings and the mixing angle between the lowest two electronic states (1 (2)A' and 2 (2)A') have been obtained to yield the corresponding quasidiabatic potential energy matrix. The strengths of the computed vibrational coupling matrix elements reflect a similar trend, as has been observed experimentally in the magnitudes of the state-to-state transition probability for the inelastic vibrational excitations [J. Krutein and F. Linder, J. Chem. Phys. 71, 559 (1979); F. A. Gianturco et al., J. Phys. B 14, 667 (1981)].     

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.     

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 study of the potential energy surface for the OH+CO-->H+CO2 reaction

Potential energy surface for the reaction OH+CO-->H+CO2 has been calculated using the complete active space self-consistent-field and multireference configuration interaction methods with the correlation consistent triple-, quadruple-, and quintuple-zeta basis sets. A specific- reaction-parameters density functional theory has been suggested, in which the B3LYP functional is reoptimized to give the highly accurate potential energy surface with less computational efforts.     

Ab initio potential energy surface of CH and reaction dynamics of H + CH+

This work presents a new ground state potential energy surface (PES) for CH. The potential is tested using quasi classical trajectory (QCT) and quantum reactive scattering methods for the H + CH(+) reaction. Cross sections and rate coefficients for all reaction channels up to 300 K are calculated. The abstraction rate coefficients follow the expected slightly decreasing behaviour above 90 K, but have a positive gradient with lower temperatures. The inelastic collision and exchange reaction rate constants are increasing monotonically with temperature. The rate coefficients of the exchange reaction differ significantly between QCT and quantum reactive scattering, due to intrinsic shortcomings of the QCT final state distributions.     

Ab initio prediction of the potential energy surface and vibration-rotation energy levels of CaCl2

The equilibrium structure and potential energy surface of calcium dichloride (CaCl2) have been determined from accurate ab initio calculations using the coupled-cluster method, CCSD(T), in conjunction with basis sets of quadruple- and quintuple-zeta quality. The CaCl2 molecule was found to be linear at equilibrium. The vibration-rotation energy levels of various CaCl2 isotopomers were predicted by the variational method. The calculated spectroscopic constants could be used to guide future high-resolution spectroscopic experiments on calcium dichloride.     

Ab initio prediction of the potential energy surface and vibration-rotation energy levels of BeH2

The equilibrium structure and potential energy surface of beryllium dihydride BeH(2) in its ground electronic state have been determined from highly accurate ab initio calculations. The vibration-rotation energy levels of three isotopomers BeH(2), BeD(2), and BeHD were predicted using the variational method. The calculated spectroscopic constants are in remarkably good agreement with the existing experimental data (sub-cm(-1) accuracy) and should be useful in a further analysis of high-resolution vibration-rotation spectra of all three isotopomers.     

Ab initio potential energy surface and intermolecular vibrations of the naphthalene-argon van der Waals complex

The intermolecular potential energy surface (PES) of the naphthalene-argon (NpAr) complex is constructed using an ab initio method. The molecule-argon interaction energy is computed at the level of the second-order M?ller-Plesset (MP2) theory combined with the augmented correlation consistent polarized valence double-ζ basis set. The analytical PES fitted to a large set of single energy values is further improved with the help of correction functions determined by calculations of the interaction energy at the coupled cluster level including single and double excitations supplemented by triple excitations performed for a limited set of intermolecular configurations. The PES determined is very flat near its four equivalent global minima of -493 cm(-1) located from both sides of the Np plane at a distance of 3.435 A? and shifted from the center of Np by ±0.43 A? along its long symmetry axis. The large-amplitude motion of Ar in the complex is investigated, and dynamical consequence of a strong intermode coupling is discovered in the excited vibrational states. The theoretical results obtained allow for the reassignment of the spectral bands observed in the electronic transition S(1) ← S(0) of the NpAr complex.     

Ab initio studies of the ground-state potential energy surface of formamide

Ab initio calculations employing a standard double-zeta basis set augmented with various polarization functions have been used to investigate the lowest energy region of the ground-state potential energy surface of the formamide molecule. Hartree-Fock calculations with only d polarization functions on the nonhydrogen atoms located two stable minima, that with geometry distorted from planarity having slightly lower energy; only one stable minimum with planar structure is found when p polarization functions on the hydrogens are included. In contrast optimizations, which account approximately for the correlation energy using second-order Møller-Plesset perturbation theory consistently favor a single slightly nonplanar minimum energy geometry.     

Ab initio potential energy surface and rovibrational spectrum of Ar-HCCCN

We report an ab initio intermolecular potential energy surface of the Ar-HCCCN complex using a supermolecular method. The calculations were performed using the fourth-order M?ller-Plesset theory with the full counterpoise correction for the basis set superposition error and a large basis set including bond functions. The complex was found to have a planar T-shaped structure minimum and a linear minimum with the Ar atom facing the H atom. The T-shaped minimum is the global minimum with the well depth of 236.81 cm(-1). A potential barrier separating the two minima is located at R=5.57 A and theta=20.39 degrees with the height of 151.59 cm(-1). The two-dimensional discrete variable representation was employed to calculate the rovibrational energy levels for Ar-HCCCN. The rovibrational spectra including intensities for the ground state and the first excited intermolecular vibrational state are also presented. The results show that the spectra are mostly b-type (Delta K(a)=+/-1) transitions with weak a-type (Delta K(a)=0) transitions in structure, which are in good agreement with the recent experimental results [A. Huckauf, W. Jager, P. Botschwina, and R. Oswald, J. Chem. Phys. 119, 7749 (2003)].     

Ab initio potential energy surfaces of the ion-molecule reaction: C2H2 + O+

High level ab initio calculations using complete active space self-consistent field and multi reference single and double excitation configuration interaction methods with cc-pVDZ (correlation consistent polarized valence double zeta) and cc-pVTZ (triple zeta) basis sets have been performed to elucidate the reaction mechanism of the ion-molecule reaction, C2H2(1Sigmag+) + O+(4S), for which collision experiment has been performed by Chiu et al. [J. Chem. Phys. 109, 5300 (1998)]. The minor low-energy process leading to the weak spin-forbidden product C2H2+ (2Piu) + O(1D) has been studied previously and will not be discussed here. The major pathways to form charge-transfer (CT) products, C2H2+ (2Piu) + O(3P) (CT1) and C2H2+ (4A2) + O(3P) (CT2), and the covalently bound intermediates are investigated. The approach of the oxygen atom cation to acetylene goes over an energy barrier TS1 of 29 kcal/mol (relative to the reactant) and adiabatically leads the CT2 product or a weakly bound intermediate Int1 between CT2 products. This transition state TS1 is caused by the avoided crossing between the reactant and CT2 electronic states. As the C-O distance becomes shorter beyond the above intermediate, the C1 reaction pathway is energetically more favorable than the Cs pathway and goes over the second transition state TS2 of a relative energy of 39 kcal/mol. Although this TS connects diabatically to the covalent intermediate Int2, there are many states that interact adiabatically with this diabatic state and these lead to the other charge-transfer product CT1 via either of several nonadiabatic transitions. These findings are consistent with the experiment, in which charge transfer and chemical reaction products are detected above 35 and 39 kcal/mol collision energies, respectively.     

Ab initio potential energy surface for HeF2 in its ground electronic state

The ground state potential energy surface for He–F₂ has been generated using the coupled-cluster singles and doubles excitation approach with perturbative treatment of triple excitations [CCSD(T)] and multi-reference configuration interaction (MRCI) methodologies, with augmented correlation consistent quadruple zeta basis set and diffused functions. Both the CCSD(T) and MRCI surfaces are compared and the results analyzed. The CCSD(T) surface exhibits van der Waals minima at different distances for different orientations of He approaching F₂ and is adequate to describe accurately only in the region around the equilibrium bond distance of F₂. The MRCI surface, on the other hand, yields reliable results for a wider range of F–F bond distances leading to the correct asymptote. Davidson correction to the MRCI surface makes it purely repulsive over the regions investigated.     

Ab initio molecular orbital studies of the structure and potential energy surface of the LiAlF4 complex

Ab initio molecular orbital calculations have carried out on various structures of LiAlF₄ complex using minimal and extended basis sets. A C_2v structure with two fluorines in the bridge was found to be more stable than structures with one and three fluorines in the bridge. Migration of the Li⁺ in the complex is found to be relatively easy and the AlF^?₄ anion is found to be distorted from tetrahedral symmetry.     

Ab initio potential energy surface and spectrum of the B(3Pi) state of the HeI2 complex

The three-dimensional interaction potential for I2(B 3Pi0u+)+He is computed using accurate ab initio methods and a large basis set. Scalar relativistic effects are accounted for by large-core relativistic pseudopotentials for the iodine atoms. Using multireference configuration interaction calculations with subsequent treatment of spin-orbit coupling, it is shown for linear and perpendicular structures of the complex that the interaction potential for I2(B 3Pi0u+)+He is very well approximated by the average of the 3A' and 3A" interaction potentials obtained without spin-orbit coupling. The three-dimensional 3A' and 3A" interaction potentials are computed at the unrestricted open-shell coupled-cluster level of theory using large basis sets. Bound state calculations based on the averaged surface are carried out and binding energies, vibrationally averaged structures, and frequencies are determined. These results are found to be in excellent accord with recent experimental measurements from laser-induced fluorescence and action spectra of HeI2. Furthermore, in combination with a recent X-state potential, the spectral blueshift is obtained and compared with available experimental values.     

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