p-Difluorobenzene-argon ground state intermolecular potential energy surface

The ground state intermolecular potential energy surface for the p-difluorobenzene-Ar van der Waals complex is evaluated using the coupled cluster singles and doubles including connected triple excitations [CCSD(T)] model and the augmented correlation consistent polarized valence double-zeta basis set extended with a set of 3s3p2d1f1g midbond functions. The surface minima are characterized by the Ar atom located above and below the difluorobenzene center of mass at a distance of 3.5290 A. The corresponding binding energy is -398.856 cm(-1). The surface is used in the evaluation of the intermolecular level structure of the complex. The results clearly improve previously available data and show the importance of using a good correlation method and basis set when dealing with van der Waals complexes.     

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     

The chlorobenzene-argon ground state intermolecular potential energy surface

Using the coupled cluster singles and doubles including connected triple excitations model with the augmented correlation consistent polarized valence double zeta basis set extended with a set of 3s3p2d1f1g midbond functions, we evaluate the ground state intermolecular potential energy surface of the chlorobenzene-argon van der Waals complex. The minima of 420 cm(-1) are characterized by Ar atom position vectors of the length 3.583 A, forming an angle of 9.87 degrees with respect to the axis perpendicular to the chlorobenzene plane. These results are compared to those obtained for similar complexes and to the experimental data available. From the potential the three-dimensional vibrational eigenfunctions and eigenvalues are calculated and the results allow to correct and complete the experimental assignment.     

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.     

The p-difluorobenzene-argon S1 excited state intermolecular potential energy surface

The first excited state (S1) intermolecular potential energy surface for the p-difluorobenzene-Ar van der Waals complex is evaluated using the coupled-cluster method and the augmented correlation consistent polarized valence double-zeta basis set extended with a set of 3s3p2d1f1g midbond functions. In order to calculate the S1 interaction energies we use the ground state surface 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 p-difluorobenzene center of mass at a distance of 3.4736 A. The corresponding interaction energy is -435.233 cm-1. The surface is used in the evaluation of the intermolecular level structure of the complex.     

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 intermolecular potential energy surface and thermophysical properties of hydrogen sulfide

A six-dimensional potential energy hypersurface (PES) for two interacting rigid hydrogen sulfide molecules was determined from high-level quantum-mechanical ab initio computations. A total of 4016 points for 405 different angular orientations of two molecules were calculated utilizing the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory and extrapolating the calculated interaction energies to the complete basis set limit. An analytical site-site potential function with eleven sites per hydrogen sulfide molecule was fitted to the interaction energies. The PES has been validated by computing the second pressure virial coefficient, shear viscosity, thermal conductivity and comparing with the available experimental data. The calculated values of volume viscosity were not used to validate the potential as the low accuracy of the available data precluded such an approach. The second pressure virial coefficient was evaluated by means of the Takahashi and Imada approach, while the transport properties, in the dilute limit, were evaluated by utilizing the classical trajectory method. In general, the agreement with the primary experimental data is within the experimental error for temperatures higher than 300 K. For lower temperatures the lack of reliable data indicates that the values of the second pressure virial coefficient and of the transport properties calculated in this work are currently the most accurate estimates for the thermophysical properties of hydrogen sulfide.     

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 ground-state potential energy function of PO+

The ground-state potential energy function of PO⁺ has been calculated from the set of molecular constants B _e, ω_e, a _i (i = 1, … , 5), R _e, D _e and C₄ in the form of generalized potential energy function previously suggested by us for solving the inverse spectroscopic problem.     

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.     

Ab initio intermolecular potential energy surface and second pressure virial coefficients of methane

A six-dimensional potential energy hypersurface (PES) for two interacting rigid methane molecules was determined from high-level quantum-mechanical ab initio computations. A total of 272 points for 17 different angular orientations on the PES were calculated utilizing the counterpoise-corrected supermolecular approach at the CCSD(T) level of theory with basis sets of aug-cc-pVTZ and aug-cc-pVQZ qualities. The calculated interaction energies were extrapolated to the complete basis set limit. An analytical site-site potential function with nine sites per methane molecule was fitted to the interaction energies. In addition, a semiempirical correction to the analytical potential function was introduced to take into account the effects of zero-point vibrations. This correction includes adjustments of the dispersion coefficients and of a single-parameter within the fit to the measured values of the second virial coefficient B(T) at room temperature. Quantitative agreement was then obtained with the measured B values over the whole temperature range of the measurements. The calculated B values should definitely be more reliable at very low temperatures (T<150 K) than values extrapolated using the currently recommended equation of state.     

An ab initio intermolecular potential energy surface for the F(2) dimer

Two analytical representations for the potential energy surface of the F(2) dimer were constructed on the basis of ab initio calculations up to the fourth-order of M?ller-Plesset (MP) perturbation theory. The best estimate of the complete basis set limit of interaction energy was derived for analysis of basis set incompleteness errors. At the MP4/aug-cc-pVTZ level of theory, the most stable structure of the dimer was obtained at R = 6.82 au, theta(a) = 12.9 degrees , theta(b) = 76.0 degrees , and phi = 180 degrees , with a well depth of 716 microE(h). Two other minima were found for canted and X-shaped configurations with potential energies around -596 and -629 microE(h), respectively. Hexadecapole moments of monomers play an important role in the anisotropy of interaction energy that is highly R-dependent at intermediate intermolecular distances. The quality of potentials was tested by computing values of the second virial coefficient. The fitted MP4 potential has a more reasonable agreement with experimental values.     

A new form for intermolecular potential energy functions

A new intermolecular potential energy function is presented for Monte Carlo and molecular dynamics simulations of liquids, solutions and other molecular assemblies. The potential energy function is expressed in terms of intermolecular overlap integrals over localized molecular orbitais of isolated molecules and Coulomb potentials between fractional point charges placed on the nuclei. The potential function is easy to generate and is applicable to a wide range of molecules. As examples the potential functions are generated for the water, ammonia and hydrogen fluoride dimers.     

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.     

Determination of the intermolecular potential energy surface for (HCl)2 from vibration--rotation--tunneling spectra

An accurate and detailed semiempirical intermolecular potential energy surface for (HCl)2 has been determined by a direct nonlinear least-squares fit to 33 microwave, far-infrared and near-infrared spectroscopic quantities using the analytical potential model of Bunker et al. [J. Mol. Spectrosc. 146, 200 (l99l)] and a rigorous four-dimensional dynamical method (described in the accompanying paper). The global minimum (De= -692 cm-1) is located near the hydrogen-bonded L-shaped geometry (R=3.746 angstroms, theta1=9 degrees, theta2=89.8 degrees, and phi=180 degrees). The marked influence of anisotropic repulsive forces is evidenced in the radial dependence of the donor-acceptor interchange tunneling pathway. The minimum energy pathway in this low barrier (48 cm-1) process involves a contraction of 0.1 angstroms in the center of mass distance (R) at the C2h symmetry barrier position. The new surface is much more accurate than either the ab initio formulation of Bunker et al. or a previous semiempirical surface [J. Chem. Phys. 78, 6841 (1983)].     

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcalmol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively.     

Structure of the adiabatic ground-state potential surface for the ozone molecule

Institute for Physical and Organic Chemistry, Rostov University. Translated from Zhurnal Strukturnoi Khimii, Vol. 32, No. 6, pp. 15–20, November–December, 1991.     

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.     

The ab initio ground-state potential energy function of beryllium monohydride, BeH

The accurate ground-state potential energy function of beryllium monohydride, BeH, has been determined from large-scale ab initio calculations using the multi-reference averaged coupled-pair functional (MR-ACPF) method in conjunction with the correlation-consistent core-valence basis sets up to septuple-zeta quality. The effects of electron correlation beyond the MR-ACPF level of approximation were taken into account. The scalar relativistic and adiabatic (the diagonal correction) effects, as well as some of the nonadiabatic effects, were also discussed. The vibration-rotation energy levels of three isotopologues, BeH, BeD, and BeT, were predicted to sub-cm(-1) accuracy.