Accurate intermolecular ground state potential of the Ar-N2 van der Waals complex

After carrying out a systematic basis set convergence study, we evaluate several ground state potential energy surfaces of the Ar-N(2) van der Waals complex at the coupled cluster singles and doubles model including connected triples corrections. We use the aug-cc-pVXZ (X=5,Q,D) and the daug-cc-pVQZ basis sets augmented with a set of 3s3p2d1f1g (denoted 33211) and 3s3p2d2f1g (denoted 33221) midbond functions, respectively. aug-cc-pVTZ-33211 results were available in the literature. The aug-cc-pV5Z-33211 (daug-cc-pVQZ-33221) surface is characterized by a T-shaped minimum at R(e)=3.709 (3.701) A and of 99.01 (102.50) cm(-1), and a linear saddle point at 4.260 (4.257) A and D(e)=75.28 (79.73) cm(-1). These results are compared with the values provided by the semiempirical potentials available, and those of previous theoretical studies. The basis set convergence of the intermolecular potentials is also analyzed. From the potentials the rovibronic spectroscopic properties are determined. We study the basis set convergence of the rotational frequencies. The binding parameters that characterized the aug-cc-pVTZ-33211 surface are reasonable, but the surface is not good enough to evaluate the microwave spectra. The aug-cc-pVQZ-33211 basis set results considerably improve the triple zeta and are close to the aug-cc-pV5Z-33211. Considering the small differences between the quadruple and the quintuple zeta surfaces, the latter results can be expected to be close to convergence. At this level the differences with respect to the accurate experimental frequencies are in the order of 0.7%. In the case of the daug-cc-pVXZ-33211,33221 (X=5,Q,T,D) series, the convergence of the interaction energies with respect to basis set improvement is not so smooth. The errors in the frequencies obtained with the daug-cc-pVQZ-33221 basis set with respect to experiment are in the order of 0.4%.     

Basis sets for the evaluation of van der Waals complex interaction energies: Ne–N2 intermolecular potential and microwave spectrum

In order to obtain efficient basis sets for the evaluation of van der Waals complex intermolecular potentials, we carry out systematic basis set studies. For this, interaction energies at representative geometries on the potential energy surfaces are evaluated using the CCSD(T) correlation method and large polarized LPol‐n and augmented polarization‐consistent aug‐pc‐2 basis sets extended with different sets of midbond functions. On the basis of the root mean square errors calculated with respect to the values for the most accurate potentials available, basis sets are selected for fitting the corresponding interaction energies and getting analytical potentials. In this work, we study the Ne–N₂ van der Waals complex and after the above procedure, the aug‐pc‐2–3321 and the LPol‐ds‐33221 basis set results are fitted. The obtained potentials are characterized by T‐shaped global minima at distances between the Ne atom and the N₂ center of mass of 3.39 Å, with interaction energies of ?49.36 cm^?1 for the aug‐pc‐2–3321 surface and ?50.28 cm^?1 for the LPol‐ds‐33221 surface. Both sets of results are in excellent agreement with the reference surface. To check the potentials further microwave transition frequencies are calculated that agree well with the experimental and the aV5Z‐33221 values. The success of this study suggests that it is feasible to carry out similar accurate calculations of interaction energies and ro‐vibrational spectra at reduced cost for larger complexes than has been possible hitherto. © 2013 Wiley Periodicals, Inc.     

Accurate intermolecular ground-state potential-energy surfaces of the HCCH-He, Ne, and Ar van der Waals complexes

Accurate ground-state intermolecular potential-energy surfaces are obtained for the HCCH-He, Ne, and Ar van der Waals complexes. The interaction energies are calculated at the coupled cluster singles and doubles including connected triple excitations level and fitted to analytic functions. For the three complexes we start with systematic basis set studies carried out at several intermolecular geometries, and using augmented correlation consistent polarized valence basis sets x-aug-cc-pVXZ (x=-,d; X=D,T,Q,5), also extended with a set of 3s3p2d1f1g midbond functions. The aug-cc-pVQZ-33211 surfaces of HCCH-He, Ne, and Ar complexes are characterized by absolute minima of -24.22, -50.20, and -122.17 cm(-1) at distances R between the rare-gas atom and the HCCH centers of mass of 4.35, 3.95, and 3.99 A, respectively; and at angles between the vector R and the HCCH main symmetry axis of 0 degrees , 43.3 degrees , and 60.6 degrees . The results are compared and considerably improve those previously available.     

Ab initio study of the structure, bonding, vibrational spectra, and energetics of XBS+ (where X=H, F, and Cl)

High level ab initio electronic structure calculations at the CCSD(T) level with augmented correlation-consistent basis set extrapolated to complete basis set limit have been performed on XBS and XBS+ for X=H, F, and Cl. The geometries have been optimized up through the aug-cc-pV5Z level and the vibrational frequencies have been calculated with the aug-cc-pVQZ basis sets. Analysis of the bonding in XBS and XBS+ using natural bond orbital analysis shows that the BS bond in XBS is a triple bond, while in XBS+ it is a double bond. The energetic properties of XBS cation and its first excited state are reported. The calculated adiabatic ionization potential is 11.11+/-0.01 eV as compared to the experimental value of 11.11+/-0.03 eV for HBS. The adiabatic ionization potentials for FBS and CIBS are 10.89+/-0.01 and 10.57+/-0.01 eV, respectively.     

Interaction potentials for Br(-)-Rg (Rg=He-Rn): spectroscopy and transport coefficients

High-level ab initio CCSD(T) calculations are performed in order to obtain accurate interaction potentials for the Br(-) anion interacting with each rare gas (Rg) atom. For the Rg atoms from He to Ar, two approaches are taken. The first one implements a relativistic core potential and an aug-cc-pVQZ basis set for bromine, an aug-cc-pV5Z basis set for Rg, and a set of bond functions placed at the midpoint of the Rg-Br distance. The second one uses the all-electron approximation with aug-cc-pV5Z bases further augmented by an extra diffuse function in each shell. Comparison reveals close similarity between both sets of results, so for Rg atoms from Kr to Rn only the second approach is exploited. Calculated potentials are assessed against the previous empirical, semiempirical, and ab initio potentials, and against available beam scattering data, zero electron kinetic energy spectroscopic data, and various sets of the measured ion mobilities and diffusion coefficients. This multiproperty analysis leads to the conclusion that the present potentials are consistently good for the whole series of Br(-)-Rg pairs over the whole range of internuclear distances covered.     

Interaction of CO with Kr: Potential energy surface and bound states

The first ab initio potential energy surface of the Kr-CO complex is developed using single and double excitation coupled-cluster theory with noniterative treatment of triple excitations. Mixed basis sets, aug-cc-pVQZ for the C and O atoms and aug-cc-pVQZ-PP for the Kr atom, with an additional (3s3p2d2f1g) set of midbond functions are used. The computed interaction energies in 336 configurations are analytically fitted to a two-dimensional potential model by a least squares fit. The potential has a minimum of -119.68 cm(-1) with Re=7.35a 0 at an approximate T-shaped geometry (theta e=98.5 degrees ). Bound state energies are calculated up to J=12, thus enabling a comprehensive comparison between theory and available experimental data as well as providing detailed guidance for future spectroscopic investigations of higher-lying states. The predicted transition frequencies and spectroscopic constants are in good agreement with the experimental results.     

Vibrational energies of PH3 calculated variationally at the complete basis set limit

The potential energy surface for the electronic ground state of PH(3) was calculated at the CCSD(T) level using aug-cc-pV(Q+d)Z and aug-cc-pVQZ basis sets for P and H, respectively, with scalar relativistic corrections included. A parametrized function was fitted through these ab initio points, and one parameter of this function was empirically adjusted. This analytical PES was employed in variational calculations of vibrational energies with the newly developed program TROVE. The convergence of the calculated vibrational energies with increasing vibrational basis set size was improved by means of an extrapolation scheme analogous to the complete basis set limit schemes used in ab initio electronic structure calculations. The resulting theoretical energy values are in excellent agreement with the available experimentally derived values.     

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

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

Comment on: “Estimating the Hartree–Fock limit from finite basis set calculations” [Jensen F (2005) Theor Chem Acc 113:267]

We demonstrate that a minor modification of the extrapolation proposed by Jensen [(2005): Theor Chem Acc 113: 267] yields very reliable estimates of the Hartree–Fock limit in conjunction with correlation consistent basis sets. Specifically, a two-point extrapolation of the form yields HF limits E _HF,∞ with an RMS error of 0.1 millihartree using aug-cc-pVQZ and aug-cc-pV5Z basis sets, and of 0.01 millihartree using aug-cc-pV5Z and aug-cc-pV6Z basis sets.     

Structural and spectroscopic study of the van der Waals complex of CO with HCO+ and the isoelectronic complex of CS with HCS+

This work reports the results of high level ab initio calculations of the OC-HCO(+) complex and the SC-HCS(+) complex and their hydrogen migration transition states. Geometry optimizations are performed at the CCSD(T)/aug-cc-pV5Z level of theory. Subsequent frequency calculations are carried out at the CCSD(T)/aug-cc-pVQZ level of theory. Additional geometry optimizations and harmonic frequency calculations for all the species involved in this study have been done with the explicitly correlated CCSD(T)-F12 method with the aug-cc-pVTZ and VTZ-F12 basis set. The geometries, rotational constants, harmonic vibrational frequencies, and energetics of the species involved in the complex are reported. These methods result in accurate computational predictions that have mean deviations for bond lengths, rotational constants, and vibrational frequencies of 0.001 A?, 163 MHz, and 46 cm(-1), respectively. These results provide essential spectroscopic properties for the complexes that can facilitate both laboratory and interstellar observations, and they also provide a comparison between oxygen and sulfur complex observability based on thermodynamic stability.     

Ab initio ground- and excited-state intermolecular potential energy surfaces for the NO-Ne and NO-Ar van der Waals complexes

The ground- [NO(X(2)Π)] and excited-state [NO(A(2)Σ(+))] intermolecular potential energy surfaces (IPESs) of the NO-Ne and NO-Ar van der Waals complexes are evaluated using the RCCSD(T) spin-restricted coupled cluster method and d-aug-cc-pVQZ basis set extended with a set of 3s3p2d1f1g midbond functions. These bases are selected from the results of a systematic basis-set convergence study carried out for the NO(A(2)Σ(+))-Ar state. We fit the interaction energies to analytic functions and compare the results to those previously available. The NO-Ar (NO-Ne) IPESs are characterized by absolute minima of -120 and -75 cm(-1) (-58 and -5 cm(-1)) at the ground and first excited state, respectively, located close to the T-shaped geometries for the ground states and at linear dispositions in the case of the excited states. The potentials are further used in the evaluation of the rovibrational spectra of the complexes, and the results are compared to those available in the literature.     

An ab initio benchmark study of hydrogen bonded formamide dimers

The five singly and doubly hydrogen bonded dimers of formamide are calculated at the correlated level by using resolution of identity M?ller-Plesset second-order perturbation theory (RIMP2) and the coupled cluster with singles, doubles, and perturbative triples [CCSD(T)] method. All structures are optimized with the Dunning aug-cc-pVTZ and aug-cc-pVQZ basis sets. The binding energies are extrapolated to the complete basis set (CBS) limit by using the aug-cc-pVXZ (X = D, T, Q) basis set series. The effect of extending the basis set to aug-cc-pV5Z on the geometries and binding energies is studied for the centrosymmetric doubly N-H...O bonded dimer FA1 and the doubly C-H...O bonded dimer FA5. The MP2 CBS limits range from -5.19 kcal/mol for FA5 to -14.80 kcal/mol for the FA1 dimer. The DeltaCCSD(T) corrections to the MP2 CBS limit binding energies calculated with the 6-31+G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ basis sets are mutually consistent to within < or =0.03 kcal/mol. The DeltaCCSD(T) correction increases the binding energy of the C-H...O bonded FA5 dimer by 0.4 kcal/mol or approximately 9% over the distance range +/-0.5 Angstrom relative to the potential minimum. This implies that the ubiquitous long-range C-H...O interactions in proteins are stronger than hitherto calculated.     

Microsolvation of protonated methane: structures and energetics of CH5(+)(H2)n

Effects of microsolvating CH(5)(+) with up to four H(2) molecules have been investigated in terms of structures and energies. For the smaller complexes, benchmark calculations have been carried out using MP2 and CCSD(T) with basis sets up to aug-cc-pV5Z quality and energies have been extrapolated to the infinite basis set limit. It is found that MP2 calculations using the aug-cc-pVQZ basis set or better yield robust reference data for both structures and energies. More than 30 stationary points including minima and first-order as well as second-order stationary points have been characterized by this method and are discussed in terms of solvation motifs. Finally, the performance of several density functionals has been assessed for this very demanding case. Popular GGA functionals such as BLYP and PBE fail, whereas the TPSS meta-GGA functional captures many structural and energetic aspects of microsolvation satisfactorily.     

Systematically convergent basis sets for explicitly correlated wavefunctions: the atoms H, He, B-Ne, and Al-Ar

Correlation consistent basis sets have been optimized for use with explicitly correlated F12 methods. The new sets, denoted cc-pVnZ-F12 (n=D,T,Q), are similar in size and construction to the standard aug-cc-pVnZ and aug-cc-pV(n+d)Z basis sets, but the new sets are shown in the present work to yield much improved convergence toward the complete basis set limit in MP2-F12/3C calculations on several small molecules involving elements of both the first and second row. For molecules containing only first row atoms, the smallest cc-pVDZ-F12 basis set consistently recovers nearly 99% of the MP2 valence correlation energy when combined with the MP2-F12/3C method. The convergence with basis set for molecules containing second row atoms is slower, but the new DZ basis set still recovers 97%-99% of the frozen core MP2 correlation energy. The accuracy of the new basis sets for relative energetics is demonstrated in benchmark calculations on a set of 15 chemical reactions.     

Density functional theory and the correlation consistent basis sets: the tight d effect on HSO and HOS

The HSO and HOS isomers have been revisited using the DFT functionals, B3LYP, B3PW91, and PBE, in combination with tight d-augmented correlation consistent basis sets, cc-pV(x+d)Z and aug-cc-pV(x+d)Z for second-row atoms. Structures, vibrationally averaged structures, relative energies, harmonic and anharmonic frequencies, enthalpies of formation of HSO and HOS, and the barrier for the HSO/HOS isomerization have been determined. These results were compared with results from previous DFT and ab initio studies in which the standard correlation consistent basis sets were used. The relative energies of the two isomers converge more rapidly and smoothly with respect to increasing basis set size for the tight d-augmented sets than for the standard basis sets. Our best calculations, B3PW91/aug-cc-pV(5+d)Z, for the relative energy of the isomers are in excellent agreement with previous CCSD(T) results given by Wilson and Dunning.     

A simple and efficient CCSD(T)-F12 approximation

A new explicitly correlated CCSD(T)-F12 approximation is presented and tested for 23 molecules and 15 chemical reactions. The F12 correction strongly improves the basis set convergence of correlation and reaction energies. Errors of the Hartree-Fock contributions are effectively removed by including MP2 single excitations into the auxiliary basis set. Using aug-cc-pVTZ basis sets the CCSD(T)-F12 calculations are more accurate and two orders of magnitude faster than standard CCSD(T)/aug-cc-pV5Z calculations.     

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.     

Anchoring the potential energy surface of the cyclic water trimer

Six cyclic stationary points on the water trimer potential energy surface have been fully optimized at the MP2 level with the aug-cc-pVQZ basis set. In agreement with previous work, harmonic vibrational frequencies indicate that two structures are minima, three are transition states connecting minima on the surface while the remaining stationary point is a higher-order saddle point. The 1- and n-particle limits of the electronic energies of each of these six structures were estimated by systematically varying both the basis sets and theoretical methods. The former limit was approached with the cc-pVXZ and aug-cc-pVXZ families of basis sets (X=2-7) while MP2, CCSD(T), and BD(TQ) calculations helped examine the latter. Core correlation effects have also been assessed at the MP2 level with the cc-pCVXZ series of basis sets (X=2-5). These data have been combined to provide highly accurate relative energies and dissociation energies for these stationary points.     

Calibration study of the CCSD(T)-F12a/b methods for C2 and small hydrocarbons

Explicitly correlated CCSD(T)-F12a/b methods combined with basis sets specifically designed for this technique have been tested for their ability to reproduce standard CCSD(T) benchmark data covering 16 small molecules composed of hydrogen and carbon. The standard method calibration set was obtained with very large one-particle basis sets, including some aug-cc-pV7Z and aug-cc-pV8Z results. Whenever possible, the molecular properties (atomization energies, structures, and harmonic frequencies) were extrapolated to the complete basis set limit in order to facilitate a direct comparison of the standard and explicitly correlated approaches without ambiguities arising from the use of different basis sets. With basis sets of triple-ζ quality or better, the F12a variant was found to overshoot the presumed basis set limit, while the F12b method converged rapidly and uniformly. Extrapolation of F12b energies to the basis set limit was found to be very effective at reproducing the best standard method atomization energies. Even extrapolations based on the small cc-pVDZ-F12/cc-pVTZ-F12 combination proved capable of a mean absolute deviation of 0.20 kcal/mol. The accuracy and simultaneous cost savings of the F12b approach are such that it should enable high quality property calculations to be performed on chemical systems that are too large for standard CCSD(T).