A new ab initio intermolecular potential energy surface and predicted rotational spectra of the Ar-H2S complex

We report a reliable three-dimensional ab initio intermolecular potential energy surface for the Ar-H(2)S complex with H(2)S monomer fixed at its experimental average structure. The potential energies were evaluated using the supermolecular approach at the coupled-cluster level with a large basis set including bond functions. The full counterpoise procedure was used to correct the basis set superposition error. The potential has a planar T-shaped global minimum with a well depth of 177.48 cm(-1) at the intermolecular distance of 3.72 ?. An additional planar local minimum is also found and is separated from the global minimum with an energy barrier with a height of 47.46 cm(-1). The combined radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm were employed to calculate the rovibrational energy levels for three isotopic species of Ar-H(2)S complexes (Ar-H(2)(32)S, Ar-H(2)(33)S, and Ar-H(2)(34)S). The rotational transition frequencies and structural parameters for the three isotopomers were also determined for the ground and the first excited states, which are all in good agreement with the available experimental values.     

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 intermolecular potential energy surface, bound states, and microwave spectra for the van der Waals complex Ne-HCCCN

We report two ab initio intermolecular potential energy surfaces for Ne-HCCCN using a supermolecular method. The calculations were performed at the fourth-order M?ller-Plesset (MP4) and the coupled cluster singles-and-doubles with noniterative inclusion of connected triples [CCSD(T)] levels 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 Ne atom facing the H atom. The two-dimensional discrete variable representation method was employed to calculate the rovibrational bound states. In addition, the microwave spectra including intensities for the ground vibrational state were predicted. The results show that the spectrum is dominated by b-type (DeltaK(a) = +/-1) transitions with very weak a-type (DeltaK(a) = 0) transitions. The calculated results at the CCSD(T) potential are in good agreement with those at MP4 potential.     

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.     

Spectroscopy and potential energy surface of the H2-CO2 van der Waals complex: experimental and theoretical studies

A 4-D ab initio potential energy surface is calculated for the intermolecular interaction of hydrogen and carbon dioxide, using the CCSD(T) method with a large basis set. The surface has a global minimum with a well depth of 212 cm(-1) and an intermolecular distance of 2.98 A for a planar configuration with both the O-C-O and H-H axes perpendicular to the intermolecular axis. Bound state calculations are performed for the H(2)-CO(2) van der Waals complex with H(2) in both the para and ortho spin states, and the binding energy of paraH(2)-CO(2)(50.4 cm(-1)) is found to be significantly less than that of orthoH(2)-CO(2)(71.7 cm(-1)). The surface supports 7 bound intermolecular vibrational states for paraH(2)-CO(2) and 19 for orthoH(2)-CO(2), and the lower rotational levels with J< or = 4 follow an asymmetric rotor pattern. The calculated infrared spectrum of paraH(2)-CO(2) agrees well with experiment. For orthoH(2)-CO(2), the ground state rotational levels allowed by symmetry are found to have (K(a), K(c))=(even, odd) or (odd, even). This somewhat unexpected fact enables the previously observed experimental spectrum to be assigned for the first time, in good agreement with theory, and indicates that the orientation of hydrogen is perpendicular to the intermolecular axis in the ground state of the orthoH(2)-CO(2) complex.     

Intermolecular potential and second virial coefficient of the water-nitrogen complex

The authors construct a rigid-body (five-dimensional) potential energy surface for the water-nitrogen complex using the systematic intermolecular potential extrapolation routine. The intermolecular potential is then extrapolated to the limit of a complete basis set. An analytic fit of this surface is obtained, and, using this, the global minimum energy is found. The minimum is located in an arrangement in which N2 is near the H atom of H2O, almost collinear with the OH bond. The best estimate of the binding energy is 441 cm-1 (1 cm-1 approximately 1.986 43x10(-23) J). The extrapolated potential is then used to calculate the second cross virial coefficient over a wide temperature range (100-3000 K). These calculated second virial coefficients are generally consistent with experimental data, but for the most part the former have smaller uncertainties.     

Ab initio study of the intermolecular potential energy surface in the ion-induced-dipole hydrogen-bonded O2(-)(X2Πg)-H2(X1Σg(+)) complex

This work presents the first investigation on the intermolecular potential energy surface of the ground electronic state of the O2(-)(2Πg)-H2(1Σg(+)) complex. High level correlated ab initio calculations were carried out using the Hartree-Fock spin-unrestricted coupled cluster singles and doubles including perturbative triples correction [RHF-UCCSD(T)]/aug-cc-pVXZ levels of calculations, where XZ = DZ, TZ, QZ, and 5Z. Results of full geometry optimization and the intermolecular potential energy surface (IPES) calculations show four equivalent minimum energy structures of L-shaped geometry with Cs symmetry at equilibrium along the 2A″ surface of the complex. For these equilibrium minimum energy structures, the most accurate value for the dissociation energy (De) was calculated as 1407.7 cm(-1), which was obtained by extrapolating the counterpoise (CP) corrected De values to the complete basis set (CBS) limit. This global minimum energy structure is stabilized by an ion-induced-dipole hydrogen bond. Detailed investigations of the IPES show that the collinear structure is unstable, while the C2v geometries present saddle points along the 2A″ surface. The barrier height between the two equivalent structures that differs in whether the hydrogen-bonded hydrogen atom is above or below the axis that connects centers of masses of the H2 and O2(-) moieties within the complex was calculated as 70 cm(-1). This suggests that the complex exhibits large amplitude motion. The barrier height to rotation of the H2 moiety by 180° within the complex is 1020 cm(-1). Anharmonic oscillator calculations predicted a strong H-H stretch fundamental transition at 3807 cm(-1). Results of the current work are expected to stimulate further theoretical and experimental investigations on the nature of intermolecular interactions in complexes that contain the superoxide radical and various closed-shell molecules that model atmospheric and biological molecules. These studies are fundamental to understanding the role of the O2(-) anion in chemistry in the atmosphere and in biological systems.     

Potential energy surface, van der Waals motions, and vibronic transitions in phenol-argon complex

The structure and intermolecular vibrational energy levels of the phenol-Ar complex are calculated from its potential energy surface. This surface is constructed from a large set of the interaction energy values computed using second-order Moller-Plesset perturbation theory with the augmented correlation consistent polarized valence double-zeta basis set. The global minimum in the potential energy surface corresponds to a cluster structure with Ar located over the geometric center of the phenol ring at a distance of 3.510 A and shifted by 0.1355 A towards oxygen. The calculated dissociation energy of 371 cm(-1) is in accordance with the experiment. Additional local minima higher in energy are with Ar placed in the phenol plane. However, they are too shallow to form the bound states corresponding to planar isomers. The deformation of the potential energy surface shape, created by the interaction of Ar with the phenolic oxygen, is responsible for a pronounced intermode mixing. As a result, a set of hybrid stretching-bending states appears which cannot be described in terms of the standard models. The intermode coupling is reflected in the vibronic structure of the S1-S0 electronic transition. The intensities of the vibronic bands are calculated from the electronic transition dipole moment surfaces determined using the ab initio single-excitation configuration interaction method. They allow us to correct and complete the assignment of the spectra observed in phenol-Ar, as well as in the analogous complexes of phenol with Kr and Xe.     

A new ab initio potential energy surface for the NeCO complex with the vibrational coordinate dependence

A new high quality three-dimensional potential energy surface for the Ne-CO van der Waals complex is developed using the CCSD(T) method and avqz∕avqz+33221 basis set. The ab initio calculation is performed in a total of 1365 configurations with supermolecule method. There is a single global minimum located in a nearly T-shaped geometry. The global minimum energy is -49.4090 cm(-1) at R(e)=6.40a(0) and θ(e)=82.5(°) for V(00). Using the three-dimensional potential energy surface, we have calculated bound rovibrational energy levels up to J = 10 including the Coriolis coupling terms. Compared with the experimental transition frequencies, the theoretical results are in good agreement with the experimental results.     

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.     

Ab initio investigation of the NH(X)-N2 van der Waals complex

The NH-N(2) van der Waals complex has been examined at the CCSD(T) level of theory using aug-cc-pVDZ and aug-cc-pVTZ basis sets. The full basis set superposition error correction was applied. Two minimum energy structures were located for the electronic ground state. The global minimum corresponds to a linear geometry of the complex (NH-N-N), with D(e)=236 cm(-1) and R(c.m.)=4.22 A. The secondary minimum corresponds to a T-shaped geometry of C(2v) symmetry, where the nitrogen atom of the H-N moiety points toward the center of mass of the N(2) unit, aligned with the a-inertial axis of the complex. The binding energy and R(c.m.) value for the secondary minimum were 144 cm(-1) and 3.63 A, respectively. This potential energy surface is consistent with the properties of matrix-isolated NH-N(2), and it is predicted that linear NH-N(2) will be a stable complex in the gas phase at low temperatures.     

Xe-N2O复合物的分子间势能面和振转光谱的理论研究

采用CCSD(T)方法研究了范德华分子体系Xe—N2O复合物的势能面和振转光谱性质,研究表明,该势能面有两个极小点,分别对应T构型和线性Xe—ONN构型,采用离散变量表象和Lanczos算法计算了体系的振转能级,计算结果表明,CCSD(T)势能面支持97个振动束缚态,并对能级进行了指认,计算得到的Xe—N2O转动跃迁频率与实验值吻合得很好。     

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.     

Rg-HX分子间势的精确从头计算研究

在用非迭代的三重激发项来校正CCSD的CCSD(T)理论水平下,采用aug cc pVQZ基函数对He HF的分子间势进行了系统的研究.结果表明：He HF以线型结构存在.在极限基的情况下,复合物两种线型极小点结构He H F和He F H势阱深分别为46.614 cm－1和25.026 cm－1,对应He原子到HF分子质心的距离Rm分别为0.3149 nm和0.3012 nm.讨论了不同的基函数和理论方法在研究此类弱束缚态复合物的分子间势时的可靠性及其对结果的影响,并研究了HF分子中H－F键长的改变对势能的影响,同时也给出了势函数的解析形式.     

The fluorobenzene-argon S(1) excited-state intermolecular potential energy surface

We evaluate the first excited-state (S1) intermolecular potential energy surface for the fluorobenzene-Ar van der Waals complex using the coupled cluster method and the augmented correlation-consistent polarized valence double-zeta basis set extended with a set of 3s3p2d1f1g midbond functions. To calculate the S(1) interaction energies, we use ground-state interaction energies evaluated with the same basis set and the coupled cluster singles and doubles (CCSD) including connected triple excitations [CCSD(T)] model and interaction and excitation energies evaluated at the CCSD level. The surface minima are characterized by the Ar atom located above and below the fluorobenzene ring at a distance of 3.5060 A with respect to the fluorobenzene center of mass and at an angle of 5.89 degrees with respect to the axis perpendicular to the fluorobenzene plane. The corresponding interaction energy is -425.226 cm(-1). The surface is used in the evaluation of the intermolecular level structure of the complex, and the results are compared to the experimental data available and to those found in previous theoretical papers on ground-state potentials for similar complexes.     

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 and analytic intermolecular potentials for Ar-CH(3)OH

Ab initio calculations at the CCSD(T)/aug-cc-pVTZ level of theory were used to characterize the Ar-CH(3)OH intermolecular potential energy surface (PES). Potential energy curves were calculated for four different Ar + CH(3)OH orientations and used to derive an analytic function for the intermolecular PES. A sum of Ar-C, Ar-O, Ar-H(C), and Ar-H(O) two-body potentials gives an excellent fit to these potential energy curves up to 100 kcal mol(-1), and adding an additional r(-n) term to the Buckingham two-body potential results in only a minor improvement in the fit. Three Ar-CH(3)OH van der Waals minima were found from the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ calculations. The structure of the global minimum is in overall good agreement with experiment (X.-C. Tan, L. Sun and R. L. Kuczkowski, J. Mol. Spectrosc., 1995, 171, 248). It is T-shaped with the hydroxyl H-atom syn with respect to Ar. Extrapolated to the complete basis set (CBS) limit, the global minimum has a well depth of 0.72 kcal mol(-1) with basis set superposition error (BSSE) correction. The aug-cc-pVTZ basis set gives a well depth only 0.10 kcal mol(-1) smaller than this value. The well depths of the other two minima are within 0.16 kcal mol(-1) of the global minimum. The analytic Ar-CH(3)OH intermolecular potential also identifies these three minima as the only van der Waals minima and the structures predicted by the analytic potential are similar to the ab initio structures. The analytic potential identifies the same global minimum and the predicted well depths for the minima are within 0.05 kcal mol(-1) of the ab initio values. Combining this Ar-CH(3)OH intermolecular potential with a potential for a OH-terminated alkylthiolate self-assembled monolayer surface (i.e., HO-SAM) provides a potential to model Ar + HO-SAM collisions.     

Five-dimensional ab initio potential energy surface and predicted infrared spectra of H2-CO2 van der Waals complexes

The authors present a new five-dimensional potential energy surface for H2-CO2 including the Q3 normal mode for the nu3 antisymmetric stretching vibration of the CO2 molecule. The potential energies were calculated using the supermolecular approach with the full counterpoise correction at the CCSD(T) level with an aug-cc-pVTZ basis set supplemented with bond functions. The global minimum is at two equivalent T-shaped coplanar configurations with a well depth of 219.68 cm-1. The rovibrational energy levels for four species of H2-CO2 (paraH2-, orthoH2-, paraD2-, and orthoD2-CO2) were calculated employing the discrete variable representation (DVR) for radial variables and finite basis representation (FBR) for angular variables and the Lanczos algorithm. Our calculations showed that the off-diagonal intra- and intermolecular vibrational coupling could be neglected, and separation of the intramolecular vibration by averaging the total Hamiltonian with the wave function of a specific vibrational state of CO2 should be a good approximation with high accuracy. The calculated band origin shift in the infrared spectra in the nu3 region of CO2 is -0.113 cm-1 for paraH2-CO2 and -0.099 cm-1 for orthoH2-CO2, which agrees well with the observed values of -0.198 and -0.096 cm-1. The calculated rovibrational spectra for H2-CO2 are consistent with the available experimental spectra. For D2-CO2, it is predicted that only a-type transitions occur for paraD2-CO2, while both a-type and b-type transitions are significant for orthoD2-CO2.     

Three-dimensional ab initio potential-energy surface and rovibrational spectra of the H2-Kr complex

We report a three-dimensional ab initio potential-energy surface for the H2-Kr complex calculated using a supermolecular method. The electronic calculations were performed at the coupled-cluster singles and doubles level with noniterative inclusion of connected triples levels with a large basis set including midbond functions and the full counterpoise correction for the basis-set superposition error. The intermolecular potential energy between the H2 molecule and the Kr atom were evaluated at five potential-optimized discrete variable representation (DVR) grid points generated from the potential-energy curve of H2. The potential for other bond lengths of H2 could be deduced using polynomial interpolations. The complex is found to have a linear preferred structure with a rather flat energy barrier. The three-dimensional DVR method and the Lanczos propagation algorithm were employed to calculate the rovibrational states without separating the inter- and intramolecular nuclear motions. In addition, the rovibrational spectra from the H2 fundamental vibrational band were calculated. The calculated shift for the band origin is -1.50 cm-1, which is in good agreement with the experimental value of -1.706 cm-1, and the calculated transition frequencies in Q1(0) and S1(0) bands are within 3% of the observed values.     

Intermolecular potential energy surface and rovibrational spectra of the He-N2O complex from ab initio calculations

We report an ab initio intermolecular potential energy surface calculation on the He-N(2)O complex with N(2)O at its ground state using a supermolecular approach. The calculation was performed at the coupled-cluster [CCSD(T)] level, with the full counterpoise correction for the basis set superposition error and a large basis set including midpoint bond functions. The CCSD(T) potential is found to have two minima corresponding to the T-shaped and linear He-ONN structures. The T-shaped minimum is the global minimum. The two-dimensional discrete variable representation method was employed to calculate the rovibrational energy levels for (4)He-N(2)O and (3)He-N(2)O with N(2)O at its ground and nu(3) excited states. The results indicate that the CCSD(T) potential supports five and four vibrational bound states for the (4)He-N(2)O and (3)He-N(2)O, respectively. Moreover, the calculations on the line intensities of the rotational transitions in the nu(3) region of N(2)O for the ground vibrational state shows that the (3)He-N(2)O spectrum is dominated by a-type transitions (DeltaK(a)=0), while the (4)He-N(2)O spectrum is contributed by both the a-type and b-type (DeltaK(a)=+/-1) transitions. The calculated transition frequencies and the intensities are in good agreement with the observed results.