Ab initio calculations, potential representation and vibrational dynamics of He2Br2 van der Waals complex

An intermolecular potential energy surface for He(2)Br(2) complex in the ground state is calculated at the levels of fourth-order (MP4) Moller-Plesset and coupled-cluster [CCSD(T)] approximations, using large-core pseudopotential for Br atoms and the aug-cc-pV5Z basis set for He. The surface is characterized by three minima and the minimum energy pathways through them. The global minimum corresponds to a linear He-Br(2)-He configuration, while the two other ones to "police-nightstick" and tetrahedral structures. The corresponding well depths are -90.39/-89.18, -81.23/-80.78 and -74.40/-74.02 cm(-1), respectively, at MP4/CCSD(T) levels of theory. It is found that results obtained by summing three-body parametrized HeBr(2) interactions and the He-He interaction are in very good accord with the corresponding MP4/CSSD(T) configuration energies of the He(2)Br(2). Variational calculations using a sum of three-body interactions are presented to study the bound states of the vdW He(2)Br(2) complex. The binding energy D(0) and the corresponding vibrationally averaged structure are determined for different isomers of the cluster and their comparison with the available experimental data is discussed.     

Theoretical investigation of the He-I2(E3Πg) ion-pair state: ab initio intermolecular potential and vibrational levels

We present a theoretical study on the potential energy surface and vibrational bound states of the E electronic excited state of the HeI(2) van der Waals system. The interaction energies are computed using accurate ab initio methods and large basis sets. Relativistic small-core effective core potentials in conjunction with a quintuple-zeta quality basis set are employed for the heavy iodine atoms in multireference configuration interaction calculations for the (3)A' and (3)A" states. For the representation of the potential energy surface we used a general interpolation technique for constructing potential surfaces from ab initio data based on the reproducing kernel Hilbert space method. The surface presents global and local minima for T-shaped configurations with well-depths of 33.2 and 4.6 cm(-1), respectively. Vibrational energies and states are computed through variational quantum mechanical calculations. We found that the binding energy of the HeI(2)(E) T-shaped isomer is 16.85 cm(-1), in excellent agreement with recent experimental measurements. In lieu of more experimental data we also report our predictions on higher vibrational levels and we analyze the influence of the underlying surface on them. This is the first attempt to represent the potential surface of such a highly excited electronic state of a van der Waals complex, and it demonstrates the capability of the ab initio technology to provide accurate results for carrying out reliable studies to model experimental data.     

van der Waals corrections to density functional theory calculations: Methane,ethane, ethylene,benzene, formaldehyde,ammonia, water,PBE, and CPMD

Parameters are developed for a practical application of the empirical van der Waals (vdW) correction infrastructure available in the CPMD density functional theory (DFT) code. The binding energy, geometry, and potential energy surface (PES) are examined for methane, ethane, ethylene, formaldehyde, ammonia, three benzene dimer geometries, and three benzene–water geometries. The vdW corrected results compare favorably with MP2 and CCSD(T) calculations near the complete basis set limits, and with experimental results where they are available.     

Can a single-reference approach provide a balanced description of ground and excited states? A comparison of the completely renormalized equation-of-motion coupled-cluster method with multireference quasidegenerate perturbation theory near a conical intersection and along a photodissociation coordinate in ammonia

We calculated the two lowest electronically adiabatic potential energy surfaces of ammonia in the region of the conical intersection and at a sequence of geometries along which one of the N-H bonds is broken. We employed both a multireference (MR) method and a single-reference (SR) method. The MR calculations are based on multiconfiguration quasidegenerate perturbation theory (MC-QDPT) with a 6-311+G(3df,3pd) basis set. The SR calculations, carried out with the same basis, employ the completely renormalized equation-of-motion coupled-cluster method with singles and doubles, and a noniterative treatment of triples, denoted CR-EOMCCSD(T). At 91 geometries used for comparison, including geometries near a conical intersection, the surfaces agree to 7% on average.     

A hierarchical construction scheme for accurate potential energy surface generation: an application to the F+H2 reaction

We present a hierarchical construction scheme for accurate ab initio potential energy surface generation. The scheme is based on the observation that when molecular configuration changes, the variation in the potential energy difference between different ab initio methods is much smaller than the variation for potential energy itself. This means that it is easier to numerically represent energy difference to achieve a desired accuracy. Because the computational cost for ab initio calculations increases very rapidly with the accuracy, one can gain substantial saving in computational time by constructing a high accurate potential energy surface as a sum of a low accurate surface based on extensive ab initio data points and an energy difference surface for high and low accuracy ab initio methods based on much fewer data points. The new scheme was applied to construct an accurate ground potential energy surface for the FH(2) system using the coupled-cluster method and a very large basis set. The constructed potential energy surface is found to be more accurate on describing the resonance states in the FH(2) and FHD systems than the existing surfaces.     

A six-dimensional H(2)-H(2) potential energy surface for bound state spectroscopy

We present a six-dimensional potential energy surface for the (H(2))(2) dimer based on coupled-cluster electronic structure calculations employing large atom-centered Gaussian basis sets and a small set of midbond functions at the dimer's center of mass. The surface is intended to describe accurately the bound and quasibound states of the dimers (H(2))(2), (D(2))(2), and H(2)-D(2) that correlate with H(2) or D(2) monomers in the rovibrational levels (v,j)=(0,0), (0,2), (1,0), and (1,2). We employ a close-coupled approach to compute the energies of these bound and quasibound dimer states using our potential energy surface, and compare the computed energies for infrared and Raman transitions involving these states with experimentally measured transition energies. We use four of the experimentally measured dimer transition energies to make two empirical adjustments to the ab initio potential energy surface; the adjusted surface gives computed transition energies for 56 experimentally observed transitions that agree with experiment to within 0.036 cm(-1). For 26 of the 56 transitions, the agreement between the computed and measured transition energies is within the quoted experimental uncertainty. Finally, we use our potential energy surface to predict the energies of another 34 not-yet-observed infrared and Raman transitions for the three dimers.     

Computational study of the rovibrational spectrum of (OCS)2

In this paper, we report a new intermolecular potential energy surface and rovibrational transition frequencies and line strengths computed for the OCS dimer. The potential is made by fitting energies obtained from explicitly correlated coupled-cluster calculations and fit using an interpolating moving least squares method. The rovibrational Schroedinger equation is solved with a symmetry-adapted Lanczos algorithm and an uncoupled product basis set. All four intermolecular coordinates are included in the calculation. On the potential energy surface we find, previously unknown, cross-shaped isomers and also polar and non-polar isomers. The associated wavefunctions and energy levels are presented. To identify polar and cross states we use both calculations of line strengths and vibrational parent analysis. Calculated rotational constants differ from their experimental counterparts by less than 0.001 cm(-1).     

A theoretical study of He2ICl van der Waals cluster

The structure, energetics, and dynamics of He2ICl complex in its ground state are studied by means of ab initio electronic structure and quantum-mechanical calculations. Interaction energies for selected He2ICl configurations are calculated at the coupled-cluster [CCSD(T)] level of theory using a large-core pseudopotential for the I atom and the aug-cc-pVTZ and aug-cc-pV5Z basis sets for the Cl and He atoms, respectively. The surface is characterized around its lower five minima and the minimum energy pathways through them. The global minimum of the potential corresponds to a "police-nightstick (1)" configuration, the second one to a linear, the next one to tetrahedral configuration, and the following two to "bifork" and "police-nightstick (2)" structures, with well depths of -99.12, -97.42, -88.32, -85.84, and -78.54 cm(-1), respectively. An analytical form based on the sum of the three-body parametrized HeICl interactions plus the He-He interaction is found to represent very well the tetra-atomic CSSD(T) results. The present potential expression is employed to perform variational five-dimensional quantum-mechanical calculations to study the vibrational bound states of the van der Waals He2ICl complex. Results for total angular momentum J = 0 provide the binding energy D0 and the corresponding vibrationally averaged structure for different isomers of the cluster. Comparison of these results with recent experimental observations further justifies the potential used in this work.     

A comparative basis-set study of NeH+ using coupled-cluster techniques

Unrestricted Hartree-Fock, coupled-cluster calculations are reported for the ground state of NeH⁺ using atomic basis sets of increasing size and accuracy for both Ne and H. The goal is to determine the basis set and coupled-cluster level of calculation needed to obtain a NeH⁺ potential energy curve of known accuracy. Here, it is shown that calculations using a quintuple zeta basis at the coupled-cluster singles and doubles level with noniterative triples, CCSD(T) , predict a Ne—H bond dissociation energy that is within about 0.01 eV of the exact Born–Oppenheimer molecular electronic structure result. Spectroscopic constants determined using the Simons–Parr–Finlan procedure are found to be in very good agreement with the experimental results. Calculations at the augmented quadruple zeta level for the two lowest triplet excited states of the NeH⁺ species are presented. Both of these states separate into ground-state Ne⁺ and H(1s). The resulting potential curves predict stable minima at the SCF, CCSD, and CCSD(T) levels with dissociation energies of about 0.07 eV. Spectroscopic constants from the potential curves and dissociation constants are reported. © 1994 John Wiley & Sons, Inc.     

Potential energy surface and rovibrational states of the ground Ar-HI complex

A potential energy surface for the ground electronic state of the Ar-HI van der Waals complex is calculated at the coupled-cluster with single and double excitations and a noniterative perturbation treatment of triple excitations [CCSD(T)] level of theory. Calculations are performed using for the iodine atom a correlation consistent triple-zeta valence basis set in conjunction with large-core Stuttgart-Dresden-Bonn relativistic pseudopotential, whereas specific augmented correlation consistent basis sets are employed for the H and Ar atoms supplemented with an additional set of bond functions. In agreement with previous studies, the equilibrium structure is found to be linear Ar-I-H, with a well depth of 205.38 cm(-1). Another two secondary minima are also predicted at a linear and bent Ar-H-I configurations with well depths of 153.57 and 151.57 cm(-1), respectively. The parametrized CCSD(T) potential is used to calculate rovibrational bound states of Ar-HI/Ar-DI complexes, and the vibrationally averaged structures of the different isomers are determined. Spectroscopic constants are also computed from the CCSD(T) surface and their comparison with available experimental data demonstrates the quality of the present surface in the corresponding configuration regions.     

Theoretical characterization of intermolecular vibrational states through the multi-configuration time dependent Hartree approach: the He2,3ICl clusters

Benchmark, full-dimensional calculations on the ground and excited vibrational states for the tetra-, and penta-atomic weakly bound He(2,3)ICl complexes are reported. The representation of the potential energy surfaces includes three-body HeICl potentials parameterized to coupled-cluster singles, doubles, and perturbative triples ab initio data. These terms are important in accurately describing the interactions of such highly floppy systems. The corresponding 6D/9D computations are performed with the multi-configuration time dependent Hartree method, using natural potential fits, and a mode combination scheme to optimize the computational effort in the improved relaxation calculations. For these complexes several low-lying vibrational states are computed, and their binding energies and radial/angular probability density distributions are obtained. We found various isomers which are assigned to different structural models related with combinations of the triatomic isomers, like linear, T-shaped, and antilinear ones. Comparison of these results with recent experimental data is presented, and the quantitative deviations found with respect to the experiment are discussed.     

Prediction of physicochemical properties of organic molecules using van der Waals surface electrostatic potentials

The generalized interaction properties function (GIPF) methodology developed by Politzer and coworkers, which calculated molecular surface electrostatic potential (MSESP) on a density envelope surface, was modified by calculating the MSESP on a much simpler van der Waals (vdW) surface of a molecule. In this work, vdW molecular surfaces were obtained from the fully optimized structures confirmed by frequency calculations at B3LYP/6-31G(d) level of theory. Multiple linear regressions for normal boiling point, heats of vaporization, heats of sublimation, heats of fusion, liquid density, and solid density were performed using GIPF variables from vdW model surface. Results from our model are compared with those from Politzer and coworkers. The surface-dependent beta (and gamma) values are dependent on the surface models but the surface-independent alpha and regression coefficients (r) are constant when vdW surface and density surface with 0.001 a.u. contour value are compared. This interesting phenomenon is explained by linear dependencies of GIPF variables.     

van der Waals corrected density functional calculations of the adsorption of benzene on the Cu (111) surface

We investigate the performance of several van der Waals (vdW) functionals at calculating the interactions between benzene and the copper (111) surface, using the local orbital approach in the SIESTA code. We demonstrate the importance of using surface optimized basis sets to calculate properties of pure surfaces, including surface energies and the work function. We quantify the errors created using (3 × 3) supercells to study adsorbate interactions using much larger supercells, and show non‐negligible errors in the binding energies and separation distances. We examine the eight high‐symmetry orientations of benzene on the Cu (111) surface, reporting the binding energies, separation distance, and change in work function. The optimized vdW‐DF(optB88‐vdW) functional provides superior results to the vdW‐DF(revPBE) and vdW‐DF2(rPW86) functionals, and closely matches the experimental and experimentally deduced values. This work demonstrates that local orbital methods using appropriate basis sets combined with a vdW functional can model adsorption between metal surfaces and organic molecules.     

Ab initio calculation of the low-lying vibrational states of C2H2(A) in full dimensionality

We report full-dimensional calculations of vibrational energies of trans-C2H2(A) using the code MULTIMODE and with a full-dimensional potential energy surface obtained by fitting singles and doubles coupled-cluster equations-of-motion (EOM-CCSD) energies using a [3s 2p 1d] atomic natural orbital basis. The EOM-CCSD calculations were done with the code "ACES II". We compare the properties of the potential surface to previous calculations at the trans minimum and also compare the vibrational energies to experimental ones.     

An ab initio potential energy surface and predissociative resonances of HArF

A three-dimensional potential energy surface of the ground electronic state HArF is constructed from more than 2000 ab initio points at the multireference averaged quadratic coupled-cluster level employing an augmented large basis set. The calculations indicate that the linear HArF molecule is metastable with a barrier of 0.643 eV in the atomization (HArF --> H + Ar + F) channel and a barrier of 1.017 eV in the dissociation (HArF --> Ar + HF) channel. Variational calculations of low-lying predissociative resonances of both HArF and DArF are performed on the three-dimensional potential energy surface using a complex-symmetric Lanczos propagation method, which yields both positions and widths of the resonance states. The resonance lifetime generally decreases with energy, but strong mode selectivity exists. Reasonably good agreement with experiment confirms the accuracy of our potential. These calculations provide valuable information on the stability and dynamics of HArF/DArF in its ground electronic state.     

The calculated infrared spectrum of Cl- H2O using a full dimensional ab initio potential surface and dipole moment surface

We report a full dimensional ab initio-based global potential energy surface (PES) and dipole moment surface (DMS) for Cl(-)H(2)O. Both surfaces are symmetric with respect to interchange of the H atoms. The PES is a fit to thousands of electronic energies calculated using the coupled-cluster method (CCSD(T)) with a moderately large basis (aug-cc-pVTZ). The infrared spectrum and vibrational dynamics are reported and compared to experiment. These results are analyzed by examination of wave function and the dipole surface.     

Ab initio prediction of the potential energy surface and vibrational-rotational energy levels of dialuminum monoxide, Al2O

The equilibrium structure and potential energy surface of dialuminum monoxide, Al(2)O, have been determined from large-scale ab initio calculations using the coupled-cluster method, CCSD(T), in conjunction with basis sets of triple- through quintuple-zeta quality. The effects of core-electron correlation on the calculated molecular parameters were investigated. The vibrational-rotational energy levels of the Al(2) (16)O and Al(2) (18)O isotopic species were calculated by a variational approach. The predicted energy levels are in remarkably good agreement with the available experimental spectroscopic data (from laser-induced fluorescence), demonstrating that the Al(2)O molecule is linear at equilibrium in its ground electronic state. The reported theoretical data settle controversies between the experimental studies about the equilibrium structure and assignment of vibrational fundamentals of the Al(2)O molecule.