Theoretical studies of the N2O van der Waals dimer: ab initio potential energy surface, intermolecular vibrations and rotational transition frequencies

Theoretical studies of the potential energy surface and bound states were performed for the N(2)O dimer. A four-dimensional intermolecular potential energy surface (PES) was constructed at the CCSD(T) level with aug-cc-pVTZ basis set supplemented with bond functions. Three co-planar local minima were found on this surface. They correspond to a nonpolar isomer with slipped-antiparallel planar structure and two equivalent polar isomers with slipped-parallel planar structures. The nonpolar isomer is energetically more stable than the polar ones by 162 cm(-1). To assign the fundamental vibrational frequencies for both isomers, more than 150 vibrational bound states were calculated based on this PES. The orientation of the nodal surface of the wave functions plays an important role in the assignment of disrotation and conrotation vibrational modes. The calculated vibrational frequencies are in good agreement with the available experimental data. We have also found a quantum tunneling effect between the two equivalent polar structures in the higher vibrational excited states. Rotational transition frequencies of the polar structure were also calculated. The accuracy of the PES is validated by the good agreement between theoretical and experimental results for the transition frequencies and spectroscopic parameters.     

Accurate intermolecular ground state potential of the Ne-HCl van der Waals complex

From an accurate ground state intermolecular potential energy surface we evaluate the rovibrational spectrum of the Ne-HCl van der Waals complex. The intermolecular potential is obtained by fitting a considerable number of interaction energies obtained at the coupled cluster singles and doubles including connected triple excitations level and with the augmented correlation consistent polarized valence quintuple zeta basis set extended with a set of 3s3p2d1f1g midbond functions. This basis set is selected after a systematic basis set study carried out at geometries close to those of the three main surface stationary points. The surface is characterized by two linear minima, i.e. Ne-ClH and Ne-HCl, with distances from the Ne atom to the HCl center of mass of 3.398 and 3.833 angstroms, respectively; and binding energies of -65.10 and -66.85 cm(-1), respectively. These results agree well with the experimental data available in contrast to previous theoretical results. The rovibrational spectra calculated for the different isotopic species are also compared to the experiments.     

Ab initio ground state phenylacetylene-argon intermolecular potential energy surface and rovibrational spectrum

We evaluate the phenylacetylene-argon intermolecular potential energy surface by fitting a representative number of ab initio interaction energies to an analytic function. These energies are calculated at a grid of intermolecular geometries, using the CCSD(T) method and the aug-cc-pVDZ basis set extended with a series of 3s3p2d1f1g midbond functions. The potential is characterized by two equivalent global minima where the Ar atom is located above and below the phenylacetylene plane at a distance of 3.5781 A? from the molecular center of mass and at an angle of 9.08° with respect to the axis perpendicular to the phenylacetylene plane and containing the center of mass. The calculated interaction energy is -418.9 cm(-1). To check further the potential, we obtain the rovibrational spectrum of the complex and the results are compared to the available experimental data.     

The Al+ -H(2) cation complex: rotationally resolved infrared spectrum, potential energy surface, and rovibrational calculations

The infrared spectrum of the Al(+)-H(2) complex is recorded in the H-H stretch region (4075-4110 cm(-1)) by monitoring Al(+) photofragments. The H-H stretch band is centered at 4095.2 cm(-1), a shift of -66.0 cm(-1) from the Q(1)(0) transition of the free H(2) molecule. Altogether, 47 rovibrational transitions belonging to the parallel K(a)=0-0 and 1-1 subbands were identified and fitted using a Watson A-reduced Hamiltonian, yielding effective spectroscopic constants. The results suggest that Al(+)-H(2) has a T-shaped equilibrium configuration with the Al(+) ion attached to a slightly perturbed H(2) molecule, but that large-amplitude intermolecular vibrational motions significantly influence the rotational constants derived from an asymmetric rotor analysis. The vibrationally averaged intermolecular separation in the ground vibrational state is estimated as 3.03 A, decreasing by 0.03 A when the H(2) subunit is vibrationally excited. A three-dimensional potential energy surface for Al(+)-H(2) is calculated ab initio using the coupled cluster CCSD(T) method and employed for variational calculations of the rovibrational energy levels and wave functions. Effective dissociation energies for Al(+)-H(2)(para) and Al(+)-H(2)(ortho) are predicted, respectively, to be 469.4 and 506.4 cm(-1), in good agreement with previous measurements. The calculations reproduce the experimental H-H stretch frequency to within 3.75 cm(-1), and the calculated B and C rotational constants to within approximately 2%. Agreement between experiment and theory supports both the accuracy of the ab initio potential energy surface and the interpretation of the measured spectrum.     

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)].     

Potential energy surface and rovibrational calculations for the Mg+-H2 and Mg+-D2 complexes

A three-dimensional potential energy surface is developed to describe the structure and dynamical behavior of the Mg(+)-H(2) and Mg(+)-D(2) complexes. Ab initio points calculated using the RCCSD(T) method and aug-cc-pVQZ basis set (augmented by bond functions) are fitted using a reproducing kernel Hilbert space method [Ho and Rabitz, J. Chem. Phys. 104, 2584 (1996)] to generate an analytical representation of the potential energy surface. The calculations confirm that Mg(+)-H(2) and Mg(+)-D(2) essentially consist of a Mg(+) atomic cation attached, respectively, to a moderately perturbed H(2) or D(2) molecule in a T-shaped configuration with an intermolecular separation of 2.62 A? and a well depth of D(e) = 842 cm(-1). The barrier for internal rotation through the linear configuration is 689 cm(-1). Interaction with the Mg(+) ion is predicted to increase the H(2) molecule's bond-length by 0.008 A?. Variational rovibrational energy level calculations using the new potential energy surface predict a dissociation energy of 614 cm(-1) for Mg(+)-H(2) and 716 cm(-1) for Mg(+)-D(2). The H-H and D-D stretch band centers are predicted to occur at 4059.4 and 2929.2 cm(-1), respectively, overestimating measured values by 3.9 and 2.6 cm(-1). For Mg(+)-H(2) and Mg(+)-D(2), the experimental B and C rotational constants exceed the calculated values by ～1.3%, suggesting that the calculated potential energy surface slightly overestimates the intermolecular separation. An ab initio dipole moment function is used to simulate the infrared spectra of both complexes.     

Higher energy states in the CO dimer: millimeter-wave spectra and rovibrational calculations

New extensive millimeter-wave measurements of the 12C16O dimer have been made, and more than 300 new spectral transitions have been observed in the frequency range 81-135 GHz. A joint analysis of these and previous millimeter-wave data yielded the precise location of 33 new energy levels of A+ symmetry and 20 levels of A- symmetry. These energy levels are located at 8-18 cm(-1) above the zero-point level. Some of them belong to already known stacks, and others make up 9 new stacks of the dimer. Newly determined stacks have K=0, 1, and, for the first time, 2, where K is the projection of the total angular momentum on the intermolecular axis. The energy levels from accompanying rovibrational calculations with the use of a recently developed hybrid CCSD(T)/DFT-SAPT potential are in very good agreement with experiment. Analysis of the calculated wave functions revealed that two new stacks of A+ symmetry with K=2 correspond to overall rotation of the dimer while the other newly observed stacks belong to the geared bend overtone modes. The ground vibrational states of the two "isomers" found are more or less localized at the two minima in the potential surface, whereas all the geared bend excited states show a considerable amount of delocalization.     

Spectroscopy of Ar-SH and Ar-SD. II. Determination of the three-dimensional intermolecular potential-energy surface

All the pure rotational transitions reported in the previous studies [J. Chem. Phys. 113, 10121 (2000); J. Mol. Spectrosc. 222, 22 (2003)] and newly observed rotation-vibration transitions, P = 1/2 <-- 3/2, for Ar-SH and Ar-SD [J. Chem. Phys. (2005), the preceding paper] have been simultaneously analyzed to determine a new intermolecular potential-energy surface of Ar-SH in the ground state. A Schrodinger equation considering the three-dimensional freedom of motion for an atom-diatom complex in the Jacobi coordinate, R, theta, and r, was numerically solved to obtain energies of the rovibrational levels using the discrete variable representation method. A three-dimensional potential-energy surface is determined by a least-squares fitting with initial values of the parameters for the potential obtained by ab initio calculations at the RCCSD(T)/aug-cc-pVTZ level of theory. The potential well reproduces all the observed data in the microwave and millimeter wave regions with parity doublings and hyperfine splittings. Several low-lying rovibrational energies are calculated using the new potential-energy surface. The dependence of the interaction energy between Ar and SH(2pi(i)) on the bond length of the SH monomer is discussed.     

Nitrous oxide dimer: a new potential energy surface and rovibrational spectrum of the nonpolar isomer

The spectrum of nitrous oxide dimer was investigated by constructing new potential energy surfaces using coupled-cluster theory and solving the rovibrational Schro?dinger equation with a Lanczos algorithm. Two four-dimensional (rigid monomer) global ab initio potential energy surfaces (PESs) were made using an interpolating moving least-squares (IMLS) fitting procedure specialized to describe the interaction of two linear fragments. The first exploratory fit was made from 1646 CCSD(T)/3ZaP energies. Isomeric minima and connecting transition structures were located on the fitted surface, and the energies of those geometries were benchmarked using complete basis set (CBS) extrapolations, counterpoise (CP) corrections, and explicitly correlated (F12b) methods. At the geometries tested, the explicitly correlated F12b method produced energies in close agreement with the estimated CBS limit. A second fit to 1757 data at the CCSD(T)-F12b/VTZ-F12 level was constructed with an estimated fitting error of less than 1.5?cm(-1). The second surface has a global nonpolar O-in minimum, two T-shaped N-in minima, and two polar minima. Barriers between these minima are small and some wave functions have amplitudes in several wells. Low-lying rovibrational wave functions and energy levels up to about 150?cm(-1) were computed on the F12b PES using a discrete variable representation/finite basis representation method. Calculated rotational constants and intermolecular frequencies are in very close agreement with experiment.     

Exploring the new three-dimensional ab initio interaction energy surface of the Ar-HF complex: rovibrational calculations for Ar-HF and Ar-DF with vibrationally excited diatoms

Several features and the performance of the recently published [P. Jankowski and M. Ziolkowski, Mol. Phys. 104, 2293 (2006)] three-dimensional intermolecular potential energy surface for the Ar-HF complex have been investigated. This full-dimensional surface has been obtained using the method of the local expansion of the exact interaction energy surface [P. Jankowski, J. Chem. Phys. 121, 1655 (2004)] in the Taylor series with respect to intramolecular coordinates. The interaction energies have been calculated with the coupled-cluster supermolecular method with single, double, and noniterative triple excitations. The convergence of the interaction energy with respect to the size of the basis set is discussed. The two-dimensional surfaces resulting from averaging of the full-dimensional surface over the intramolecular vibration of HF have been obtained and directly compared to the empirical H6(4,3,2) set of surfaces proposed by Hutson [J. Chem. Phys. 96, 6752 (1992)]. A very good agreement has been observed. The averaged potentials have been used to calculate the rovibrational energy levels of the Ar-HF and Ar-DF complexes and compared to the experimental data. The accuracy of rovibrational calculations achieved with the new surface is much better than with any of the ab initio surfaces available so far. Predictions of the rovibrational energy levels and spectroscopic constants have also been done for Ar-HF with HF in the v=4,5 vibrational states, and for Ar-DF and DF in the v=3,4 states. The full-dimensional surface studied in this paper is the first ab initio surface which is fully compatible with the empirical H6(4,3,2) surface proposed by Hutson.     

Isotope effects in the CO dimer: millimeter wave spectrum and rovibrational calculations of (12C18O)2

The millimeter wave spectrum of the isotopically substituted CO dimer, (12C18O)2, was studied with the Orotron jet spectrometer, confirming and extending a previous infrared study [A. R. W. McKellar, J. Mol. Spectrosc. 226, 190 (2004)]. A very dilute gas mixture of CO in Ne was used, which resulted in small consumption of 12C18O sample gas and produced cold and simple spectra. Using the technique of combination differences together with the data from the infrared work, six transitions in the 84-127 GHz region have been assigned. They belong to two branches, which connect four low levels of A+ symmetry to three previously unknown levels of A- symmetry. The discovery of the lowest state of A- symmetry, which corresponds to the projection K=0 of the total angular momentum J onto the intermolecular axis, identifies the geared bending mode of the 12C18O dimer at 3.607 cm(-1). Accompanying rovibrational calculations using a recently developed hybrid potential from ab initio coupled cluster [CCSD(T)] and symmetry-adapted perturbation theory calculations [G. W. M. Vissers et al., J. Chem. Phys. 122, 054306 (2005)] gave very good agreement with experiment. The isotopic dependence of the A+/A- energy splitting, the intermolecular separation R, and the energy difference of two ground state isomers, which change significantly when 18O or 13C are substituted into the normal (12C16O)2 isotopolog [L. A. Surin et al., J. Mol. Spectrosc. 223, 132 (2004)], was explained by these calculations. It turns out that the change in anisotropy of the intermolecular potential with respect to the shifted monomer centers of mass is particularly significant.     

Rotationally resolved infrared spectrum of the Li+-D2 cation complex

The infrared spectrum of mass selected Li(+)-D(2) cations is recorded in the D-D stretch region (2860-2950 cm(-1)) in a tandem mass spectrometer by monitoring Li(+) photofragments. The D-D stretch vibration of Li(+)-D(2) is shifted by -79 cm(-1) from that of the free D(2) molecule indicating that the vibrational excitation of the D(2) subunit strengthens the effective Li(+)cdots, three dots, centeredD(2) intermolecular interaction. Around 100 rovibrational transitions, belonging to parallel K(a)=0-0, 1-1, and 2-2 subbands, are fitted to a Watson A-reduced Hamiltonian to yield effective molecular parameters. The infrared spectrum shows that the complex consists of a Li(+) ion attached to a slightly perturbed D(2) molecule with a T-shaped equilibrium configuration and a 2.035 A vibrationally averaged intermolecular separation. Comparisons are made between the spectroscopic data and data obtained from rovibrational calculations using a recent three dimensional Li(+)-D(2) potential energy surface [R. Martinazzo, G. Tantardini, E. Bodo, and F. Gianturco, J. Chem. Phys. 119, 11241 (2003)].     

A new ab initio interaction energy surface and high-resolution spectra of the H2-CO van der Waals complex

A new four-dimensional intermolecular potential-energy surface for the H(2)-CO complex is presented. The ab initio points have been computed on a five-dimensional grid including the dependence on the H-H separation (the C-O separation was fixed). The surface has then been obtained by averaging over the intramolecular vibration of H(2). The coupled-cluster supermolecular method with single, double, and noniterative triple excitations has been used to calculate the interaction energy. The correlation part of the interaction energy has been obtained from extrapolations based on calculations in a series of basis sets. An analytical fit of the ab initio potential-energy surface has the global minimum of -93.049 cm(-1) at the intermolecular separation of 7.92 bohr for the linear geometry with the C atom pointing toward the H(2) molecule. For the other linear geometry, with the O atom pointing toward H(2), the local minimum of -72.741 cm(-1) has been found for the intermolecular separation of 7.17 bohr. The potential has been used to calculate the rovibrational energy levels of the para-H(2)-CO complex. The results agree very well with those observed by McKellar [A. R. W. McKellar J. Chem. Phys. 108, 1811 (1998)]: the discrepancies are smaller than 0.1 cm(-1). The calculated dissociation energy is equal to 19.527 cm(-1) and significantly smaller than the value of 22 cm(-1) estimated from the experiment. Predictions of rovibrational energy levels for ortho-H(2)-CO have also been done and can serve as a guidance to assign recorded experimental spectra. The interaction second virial coefficient has been calculated and compared with the experimental data.     

Quasiclassical trajectory calculations of the OH+NO2 association reaction on a global potential energy surface

We report a full-dimensional potential energy surface (PES) for the OH+NO(2) reaction based on fitting more than 55,000 energies obtained with density functional theory-B3LYP6-311G(d,p) calculations. The PES is invariant with respect to permutation of like nuclei and describes all isomers of HOONO, HONO(2), and the fragments OH+NO(2) and HO(2)+NO. Detailed comparison of the structures, energies, and harmonic frequencies of various stationary points on the PES are made with previous and present high-level ab initio calculations. Two hydrogen-bond complexes are found on the PES and confirmed by new ab initio CASPT2 calculations. Quasiclassical trajectory calculations of the cross sections for ground rovibrational OH+NO(2) association reactions to form HOONO and HONO(2) are done using this PES. The cross section to form HOONO is larger than the one to form HONO(2) at low collision energies but the reverse is found at higher energies. The enhancement of the HOONO complex at low collision energies is shown to be due, in large part, to the transient formation of a H-bond complex, which decays preferentially to HOONO. The association cross sections are used to obtain rate constants for formation of HOONO and HONO(2) for the ground rovibrational states in the high-pressure limit.     

Intermolecular interaction in an open-shell pi-bound cationic complex: IR spectrum and coupled cluster calculations for C2H2+-Ar

The intermolecular potential energy surface (PES) of Ar interacting with the acetylene cation in its (2)Pi(u) ground electronic state is characterized by infrared photodissociation (IRPD) spectroscopy and quantum chemical calculations. In agreement with the theoretical predictions, the rovibrational analysis of the IRPD spectrum of C(2)H(2) (+)-Ar recorded in the vicinity of the antisymmetric CH stretching fundamental (nu(3)) is consistent with a vibrationally averaged T-shaped structure and a ground-state center-of-mass separation of R(c.m.) = 2.86 +/- 0.09 A. The nu(3) band experiences a blueshift of 16.7 cm(-1) upon complexation, indicating that vibrational excitation slightly reduces the interaction strength. The two-dimensional intermolecular PES of C(2)H(2) (+)-Ar, obtained from coupled cluster calculations with a large basis set, features strong angular-radial coupling and supports in addition to a global pi-bound minimum also two shallow side wells with linear H-bound geometries. Bound state rovibrational energy level calculations are carried out for rotational angular momentum J = 0-10 (both parities) employing a discrete variable representation-distributed Gaussian basis method. Effective spectroscopic constants are determined for the vibrational ground state by fitting the calculated rotational energies to the standard Watson A-type Hamiltonian for a slightly asymmetric prolate top.     

Polar isomer of formic acid dimers formed in helium nanodroplets

The infrared spectrum of formic acid dimers in helium nanodroplets has been observed corresponding to excitation of the "free" OH and CH stretches. The experimental results are consistent with a polar acyclic structure for the dimer. The formation of this structure in helium, as opposed to the much more stable cyclic isomer with two O-H...O hydrogen bonds, is attributed to the unique growth conditions that exist in helium droplets, at a temperature of 0.37 K. Theoretical calculations are also reported to aid in the interpretation of the experimental results. At long range the intermolecular interaction between the two monomers is dominated by the dipole-dipole interaction, which favors the formation of a polar dimer. By following the minimum-energy path, the calculations predict the formation of an acyclic dimer having one O-H...O and one C-H...O contact. This structure corresponds to a local minimum on the potential energy surface and differs significantly from the structure observed in the gas phase.     

Six-dimensional potential energy surface of the dissociative chemisorption of HCl on Au(111) using neural networks

We constructed a six-dimensional potential energy surface(PES)for the dissociative chemisorption of HCl on Au(111)using the neural networks method based on roughly 70000 energies obtained from extensive density functional theory(DFT)calculations.The resulting PES is accurate and smooth,based on the small fitting errors and good agreement between the fitted PES and the direct DFT calculations.Time-dependent wave packet calculations show that the potential energy surface is very well converged with respect to the number of DFT data points,as well as to the fitting process.The dissociation probabilities of HCl initially in the ground rovibrational state from six-dimensional quantum dynamical calculations are quite diferent from the four-dimensional fixed-site calculations,indicating it is essential to perform full-dimensional quantum dynamical studies for the title molecule-surface interaction system.     

Multimode calculations of rovibrational energies and dipole transition intensities for polyatomic molecules with torsional motion: application to H2O2

We report rigorous calculations of rovibrational energies and dipole transition intensities for hydrogen peroxide using a new version of MULTIMODE as applied to molecules with torsional (reaction path) motion. The key features which permit such calculations for moderately sized polyatomic molecules of this general type are briefly described. A previous, accurate potential energy surface and a new high-level ab initio dipole moment surface are employed in these calculations. Detailed comparisons are made with high-resolution experimental spectral intensities from the HITRAN database.     

Ne-HCl势能面和振转光谱的理论研究

利用量子化学计算方法CCSD(T)和大基组aug-cc-pVTZ加键函数3_s3_p2_d对Ne-HCl体系的分子间势能面进行了理论研究.结果表明,势能面上有两个势阱,分别对应于线性Ne-ClH和Ne-HCl构型.通过精确求解核运动方程发现,该从头算势能面分别支持5个(对Ne-HCl)和7个(Ne-DCl)振动束缚态.计算得到的振转跃迁频率与实值值吻合.     

Ab initio calculations of the rotationally resolved infrared spectrum of KNa 2 +

Summary Ab initio variational calculations were performed on the rotationally resolved infrared spectrum of KNa ₂ ⁺ . A discrete potential energy surface was generated using the configuration interaction ansatz coupled with the frozen core approximation, from which an analytical representation was obtained using a power series expansion employing a Dunham expansion variable. This force field was embedded in an Eckart-Watson rovibrational Hamiltonian, from which eigenfunctions and eigenenergies were calculated. An SCF dipole moment surface was generated and used to calculate absolute line intensities and square dipole matrix elements between the vibrational ground state and the lowest-lying excited states for some of the most intense transitions within the P, Q and R branches.