Basis set convergence of molecular correlation energy differences within the random phase approximation

The basis set convergence of energy differences obtained from the random phase approximation (RPA) to the correlation energy is investigated for a wide range of molecular interactions. For dispersion bound systems the basis set incompleteness error is most pronounced, as shown for the S22 benchmark [P. Jurecka et al., Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The use of very large basis sets (> quintuple-zeta) or extrapolation to the complete basis set (CBS) limit is necessary to obtain a reliable estimate of the binding energy for these systems. Counterpoise corrected results converge to the same CBS limit, but counterpoise correction without extrapolation is insufficient. Core-valence correlations do not play a significant role. For medium- and short-range correlation, quadruple-zeta results are essentially converged, as demonstrated for relative alkane conformer energies, reaction energies dominated by intramolecular dispersion, isomerization energies, and reaction energies of small organic molecules. Except for weakly bound systems, diffuse augmentation almost universally slows down basis set convergence. For most RPA applications, quadruple-zeta valence basis sets offer a good balance between accuracy and efficiency.     

Improved efficiency of focal point conformational analysis with truncated correlation consistent basis sets

It has been suggested that the computational cost of correlated ab initio calculations could be reduced efficiently by using truncated basis sets on hydrogen atoms (Mintz et al., J Chem Phys 2004, 121, 5629). We now explore this proposal in the context of conformational analysis of small molecules, such as hydrogen peroxide, dimethyl ether, ethyl methyl ether, formic acid, methyl formate, and several small alcohols. It is found that truncated correlation consistent basis sets that lack certain higher angular momentum functions on hydrogen atoms offer accuracy similar to traditional Dunning's basis sets for conformational analysis. Combination of such basis sets with the basis set extrapolation technique to estimate Hartree-Fock and M?ller-Plesset second order energies provides composite extrapolation model chemistries that are significantly more accurate and faster than analogous single point calculations with traditional correlation consistent basis sets. Root mean square errors of best composite extrapolation model chemistries on the used set of molecules are within 0.03 kcal/mol of traditional focal point conformational energies. The applicability of composite extrapolation methods is illustrated by performing conformational analysis of tert-butanol and cyclohexanol. For comparison, conformational energies calculated with popular molecular mechanics force fields are also given.     

Extrapolating to the one-electron basis-set limit in electronic structure calculations

A simple, yet reliable, scheme based on treating uniformly singlet-pair and triplet-pair interactions is suggested to extrapolate atomic and molecular electron correlation energies calculated at two basis-set levels of ab initio theory to the infinite one-electron basis-set limit. The novel dual-level method is first tested on extrapolating the full correlation in single-reference coupled-cluster singles and doubles energies for the closed-shell systems CH2((1)A1), H2O, HF, N2, CO, Ne, and F2 with correlation-consistent basis sets of the type cc-pVXZ (X=D,T,Q,5,6) reported by Klopper [Mol. Phys. 6, 481 (2001)] against his own benchmark calculations with large uncontracted basis sets obtained from explicit correlated singles and doubles coupled-cluster theory. Comparisons are also reported for the same data set but using both single-reference Moller-Plesset and coupled-cluster doubles methods. The results show a similar, often better, accordance with the target results than Klopper's extrapolations where singlet-pair and triplet-pair energies are extrapolated separately using the popular X(-3) and X(-5) dual-level laws, respectively. Applications to the extrapolation of the dynamical correlation in multireference configuration interaction calculations carried out anew for He, H2, HeH+, He2 ++, H3+(1 (1)A'), H3+(1 (3)A'), BH, CH, NH, OH, FH, B2, C2, N2, O2, F2, BO, CO, NO, BN, CN, SH, H2O, and NH3 with standard augmented correlation-consistent basis sets of the type aug-cc-pVXZ (X=D,T,Q,5,6) are also reported. Despite lacking accurate theoretical or experimental data for comparison in the case of most diatomic systems, the new method also shows in this case a good performance when judged from the results obtained with the traditional schemes which extrapolate using the two largest affordable basis sets. For the Hartree-Fock and complete-active space self-consistent field energies, a simple pragmatic extrapolation rule is examined whose results are shown to compare well with the ones obtained from the best reported schemes.     

Accuracy of basis-set extrapolation schemes for DFT-RPA correlation energies in molecular calculations

We construct a reference benchmark set for atomic and molecular random phase approximation (RPA) correlation energies in a density functional theory framework at the complete basis-set limit. This set is used to evaluate the accuracy of some popular extrapolation schemes for RPA all-electron molecular calculations. The results indicate that for absolute energies, accurate results, clearly outperforming raw data, are achievable with two-point extrapolation schemes based on quintuple- and sextuple-zeta basis sets. Moreover, we show that results in good agreement with the benchmark can also be obtained by using a semiempirical extrapolation procedure based on quadruple- and quintuple-zeta basis sets. Finally, we analyze the performance of different extrapolation schemes for atomization energies.     

Accurate calculations of intermolecular interaction energies using explicitly correlated wave functions

Explicitly correlated second-order M?ller-Plesset (MP2-F12) calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented and compared with standard MP2 results. The MP2 complete basis set limits are estimated using basis set extrapolation and augmented quadruple-zeta and quintuple-zeta basis sets. Already with augmented double-zeta basis sets the MP2-F12 interaction energies are found to be closer to the complete basis set limits than standard MP2 calculations with augmented quintuple-zeta basis sets. Various possible approximations in the MP2-F12 method are systematically tested. Best results are obtained with localized orbitals and the diagonal MP2-F12/C(D) ansatz. Hybrid approximations, in which some contributions of the auxiliary basis set are neglected and which considerably reduce the computational cost, have a negligible effect on the interaction energies. Also the orbital-invariant fixed-amplitude approximation of Ten-no leads to only slightly less accurate results. Preliminary results for the neon and benzene dimers, obtained with the recently proposed CCSD(T)-F12a approximation, indicate that the CCSD(T) basis set limits can also be very closely approached using augmented triple-zeta basis sets.     

Intermolecular potentials of the methane dimer calculated with Møller-Plesset perturbation theory and density functional theory

We have calculated the intermolecular interaction potentials of the methane dimer at the minimum-energy D(3d) conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with the Perdew-Wang (PW91) functional as the exchange or the correlation part. The HF calculations yield unbound potentials largely due to the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater-type orbitals fitted with Gaussian functions (STO-nG) (n=3-6) [Quantum Theory of Molecular and Solids: The Self-Consistent Field for Molecular and Solids (McGraw-Hill, New York, 1974), Vol. 4], Pople's medium size basis sets of Krishnan et al. [J. Chem. Phys. 72, 650 (1980)] [up to 6-311++G(3df,3pd)] to Dunning's correlation consistent basis sets [J. Chem. Phys. 90, 1007 (1989)] (cc-pVXZ and aug-cc-pVXZ) (X=D, T, and Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(2d,2p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy (approximately 0.01 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the destined C6 value from molecular polarizability calculations. The slow convergence could indicate the inefficacy of using the MP2 calculations with Gaussian-type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the destined potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energy calculated using the PW91PW91 functional and the equilibrium bond length calculated using the PW91VP86 functional are close to the MP2 results at the basis set limit.     

Quantum Monte Carlo calculations of the dissociation energy of the water dimer

We report diffusion quantum Monte Carlo (DMC) calculations of the equilibrium dissociation energy D(e) of the water dimer. The dissociation energy measured experimentally, D(0), can be estimated from D(e) by adding a correction for vibrational effects. Using the measured dissociation energy and the modern value of the vibrational energy Mas et al., [J. Chem. Phys. 113, 6687 (2000)] leads to D(e)=5.00+/-0.7 kcal mol(-1), although the result Curtiss et al., [J. Chem. Phys. 71, 2703 (1979)] D(e)=5.44+/-0.7 kcal mol(-1), which uses an earlier estimate of the vibrational energy, has been widely quoted. High-level coupled cluster calculations Klopper et al., [Phys. Chem. Chem. Phys. 2, 2227 (2000)] have yielded D(e)=5.02+/-0.05 kcal mol(-1). In an attempt to shed new light on this old problem, we have performed all-electron DMC calculations on the water monomer and dimer using Slater-Jastrow wave functions with both Hartree-Fock approximation (HF) and B3LYP density functional theory single-particle orbitals. We obtain equilibrium dissociation energies for the dimer of 5.02+/-0.18 kcal mol(-1) (HF orbitals) and 5.21+/-0.18 kcal mol(-1) (B3LYP orbitals), in good agreement with the coupled cluster results.     

State-to-state reactive differential cross sections for the H+H2-->H2+H reaction on five different potential energy surfaces employing a new quantum wavepacket computer code: DIFFREALWAVE

State-to-state differential cross sections have been calculated for the hydrogen exchange reaction, H+H2-->H2+H, using five different high quality potential energy surfaces with the objective of examining the sensitivity of these detailed cross sections to the underlying potential energy surfaces. The calculations were performed using a new parallel computer code, DIFFREALWAVE. The code is based on the real wavepacket approach of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)]. The calculations are parallelized over the helicity quantum number Omega' (i.e., the quantum number for the body-fixed z component of the total angular momentum) and wavepackets for each J,Omega' set are assigned to different processors, similar in spirit to the Coriolis-coupled processors approach of Goldfield and Gray [Comput. Phys. Commun. 84, 1 (1996)]. Calculations for J=0-24 have been performed to obtain converged state-to-state differential cross sections in the energy range from 0.4 to 1.2 eV. The calculations employ five different potential energy surfaces, the BKMP2 surface and a hierarchical family of four new ab initio surfaces [S. L. Mielke, et al., J. Chem. Phys. 116, 4142 (2002)]. This family of four surfaces has been calculated using three different hierarchical sets of basis functions and also an extrapolation to the complete basis set limit, the so called CCI surface. The CCI surface is the most accurate surface for the H3 system reported to date. Our calculations of differential cross sections are the first to be reported for the A2, A3, A4, and CCI surfaces. They show that there are some small differences in the cross sections obtained from the five different surfaces, particularly at higher energies. The calculations also show that the BKMP2 performs well and gives cross sections in very good agreement with the results from the CCI surface, displaying only small divergences at higher energies.     

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.     

On the validity of complete basis set extrapolation formula optimized for the equilibrium distance applied to the potential energy surface for the correlation energy of the helium dimer

Full configuration interaction calculations are performed for He₂ using various orbital basis sets of the aug‐cc‐pVXZ type, with the correlation energies being extrapolated to the complete basis set (CBS) limit. A two‐point CBS extrapolation formula has been utilized for such a purpose. It is shown that the extrapolation formula with the offset parameter k(R) optimized for the equilibrium distance is not uniformly applicable to He He distances in the very short region of the potential energy curve. The offset parameter k(R) in the repulsive region of the potential energy curve can be largely different with the one in the long‐range distances especially in the cases of basis‐sets with large cardinality number. It is also noticed that the accuracy of this extrapolation scheme may not be improved with the increasing of the cardinality number.     

Compatibility of correlation‐consistent basis sets with a hybrid Hartree–Fock/density functional method

The present study examines the feasibility of combining the correlation‐consistent basis sets developed by Dunning and coworkers with the hybrid Hartree–Fock/density functional method B3LYP. Furthermore, extrapolation to the complete basis set (CBS) limit minimizes errors due to the presence of an incomplete basis set and can act as a rigorous test of the limitations of the B3LYP method. Equilibrium geometries, energies, and harmonic vibrational frequencies were determined for a series of well‐studied, yet computationally challenging, small inorganics and their respective ions. The results were then extrapolated to the CBS limit, where applicable, and compared to experiment. It was found that a union between the hybrid Hartree–Fock/density functional B3LYP method and Dunning's augmented correlation‐consistent basis sets gave results that were comparable to molecular orbital methods that explicitly account for electron correlation. Furthermore, the minimum basis set necessary to attain reasonable results for the systems studied was aug‐cc‐pVTZ. Upgrading to the aug‐cc‐pVQZ level and subsequent extrapolation to the CBS limit further improved the overall agreement with the experiment. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 207–216, 1999     

The convergence of complete active space self-consistent-field energies to the complete basis set limit

Examination of the convergence of full valence complete active space self-consistent-field energies with expansion of the one-electron basis set reveals a pattern very similar to the convergence of single determinant Hartree-Fock energies. Calculations on 26 molecular examples with the sequence of ntuple-zeta augmented polarized (nZaP) basis sets (n=2, 3, 4, 5, and 6) are used to evaluate complete basis set extrapolation schemes. The most effective extrapolation reduces the rms one-electron basis set truncation errors from 3.03, 0.58, and 0.12 mhartree to 0.23, 0.05, and 0.014 mhartree for the 3ZaP, 4ZaP, and 5ZaP basis sets, respectively.     

Towards a complete basis set limit of Hartree–Fock method: correlation-consistent versus polarized-consistent basis sets

In this paper the convergence pattern of correlation-consistent (cc-pVxZ) and polarized-consistent (PC-n) hierarchies relative to the complete basis set limit have been considered in a small set of diatomic molecules. Using the sequence of these basis sets it was demonstrated that potential energy surfaces derived from basis-set-dependent solution of the Hartree–Fock equations achieves the exact numerical derived potential energy surfaces (PESs) in an ordered manner. So it was possible to compute the spectroscopic parameters in the complete basis set limit with considerable accuracy using the most extended members of both hierarchies. On the other hand, for the first time the detailed convergence patterns of total energies in three separate inter-nuclear distances have been considered in these molecules and it was demonstrated that the total energies arrive at microhartree accuracy at a considerable rate. Possible performance of extrapolation schemes is discussed and it was demonstrated that reliable extrapolation procedures indeed exist. A successful test of the proposed extrapolation method, using the three most extended members of polarized-consistent basis sets, has been accomplished on selected polyatomic molecules.     

Coupled-cluster dynamic polarizabilities including triple excitations

Dynamic polarizabilities for open- and closed-shell molecules were obtained by using coupled-cluster (CC) linear response theory with full treatment of singles, doubles, and triples (CCSDT-LR) with large basis sets utilizing the NWChem software suite. By using four approximate CC methods in conjunction with augmented cc-pVNZ basis sets, we are able to evaluate the convergence in both many-electron and one-electron spaces. For systems with primarily dynamic correlation, the results for CC3 and CCSDT are almost indistinguishable. For systems with significant static correlation, the CC3 tends to overestimate the triples contribution, while the PS(T) approximation [J. Chem. Phys. 127, 164105 (2007)] produces mixed results that are heavily dependent on the accuracies provided by noniterative approaches used to correct the equation-of-motion CCSD excitation energies. Our results for open-shell systems show that the choice of reference (restricted open-shell Hartree-Fock versus unrestricted Hartree-Fock) can have a significant impact on the accuracy of polarizabilities. A simple extrapolation based on pentuple-zeta CCSD calculations and triple-zeta CCSDT calculations reproduces experimental results with good precision in most cases.     

Estimated MP2 and CCSD(T) interaction energies of n-alkane dimers at the basis set limit: comparison of the methods of Helgaker et al. and Feller

The MP2 (the second-order M?ller-Plesset calculation) and CCSD(T) (coupled cluster calculation with single and double substitutions with noniterative triple excitations) interaction energies of all-trans n-alkane dimers were calculated using Dunning's [J. Chem. Phys. 90, 1007 (1989)] correlation consistent basis sets. The estimated MP2 interaction energies of methane, ethane, and propane dimers at the basis set limit [EMP2(limit)] by the method of Helgaker et al. [J. Chem. Phys. 106, 9639 (1997)] from the MP2/aug-cc-pVXZ (X=D and T) level interaction energies are very close to those estimated from the MP2/aug-cc-pVXZ (X=T and Q) level interaction energies. The estimated EMP2(limit) values of n-butane to n-heptane dimers from the MP2/cc-pVXZ (X=D and T) level interaction energies are very close to those from the MP2/aug-cc-pVXZ (X=D and T) ones. The EMP2(limit) values estimated by Feller's [J. Chem. Phys. 96, 6104 (1992)] method from the MP2/cc-pVXZ (X=D, T, and Q) level interaction energies are close to those estimated by the method of Helgaker et al. from the MP2/cc-pVXZ (X=T and Q) ones. The estimated EMP2(limit) values by the method of Helgaker et al. using the aug-cc-pVXZ (X=D and T) are close to these values. The estimated EMP2(limit) of the methane, ethane, propane, n-butane, n-pentane, n-hexane, n-heptane, n-octane, n-nonane, and n-decane dimers by the method of Helgaker et al. are -0.48, -1.35, -2.08, -2.97, -3.92, -4.91, -5.96, -6.68, -7.75, and -8.75 kcal/mol, respectively. Effects of electron correlation beyond MP2 are not large. The estimated CCSD(T) interaction energies of the methane, ethane, propane, and n-butane dimers at the basis set limit by the method of Helgaker et al. (-0.41, -1.22, -1.87, and -2.74 kcal/mol, respectively) from the CCSD(T)/cc-pVXZ (X=D and T) level interaction energies are close to the EMP2(limit) obtained using the same basis sets. The estimated EMP2(limit) values of the ten dimers were fitted to the form m0+m1X (X is 1 for methane, 2 for ethane, etc.). The obtained m0 and m1 (0.595 and -0.926 kcal/mol) show that the interactions between long n-alkane chains are significant. Analysis of basis set effects shows that cc-pVXZ (X=T, Q, or 5), aug-cc-pVXZ (X=D, T, Q, or 5) basis set, or 6-311G** basis set augmented with diffuse polarization function is necessary for quantitative evaluation of the interaction energies between n-alkane chains.     

Improved correlation energy extrapolation schemes based on local pair natural orbital methods

It is well-known that the basis set limit is difficult to reach in correlated post Hartree-Fock ab initio calculations. One possible route forward is to employ basis set extrapolation schemes. In order to avoid prohibitively expensive calculations, the highest level calculation (typically based on the "gold standard" coupled cluster theory with single, double, and perturbative triple excitations, CCSD(T)) is only performed with the smallest basis set, and the remaining basis set incompleteness is estimated at a lower level of theory, typically second-order M?ller-Plesset perturbation theory (MP2). In this work, we provide a comprehensive investigation of alternative schemes where the MP2 extrapolation is replaced by the coupled-electron pair approximation, version 1 (CEPA/1) or the local pair natural orbital version of this method (LPNO-CEPA/1). It is shown that the MP2 method achieves apparent accuracy only due to error cancellation. Systematically more accurate results at small additional computational cost are obtained if the MP2 step is replaced by LPNO-CEPA/1. The errors of LPNO-CEPA/1 relative to canonical CEPA/1 are negligible. Owing to the highly systematic nature of the deviations between canonical and LPNO methods, basis set extrapolation reduces the LPNO errors in the total energies by 1 order of magnitude (~0.2 kcal/mol) and errors in energy differences to essentially zero. Using the CCSD(T)/LPNO-CEPA/1-based extrapolation scheme, new reference values are proposed for the recently published S66 set of interaction energies. The deviations between the new values and the original interactions energies are mostly very small but reach values up to 0.3 kcal/mol.     

Variational formulation of perturbative explicitly-correlated coupled-cluster methods

We present a variational formulation of the recently-proposed CCSD(2)(R12) method [Valeev, Phys. Chem. Chem. Phys., 2008, 10, 106]. The centerpiece of this approach is the CCSD(2)(R12) Lagrangian obtained via L?wdin partitioning of the coupled-cluster singles and doubles (CCSD) Hamiltonian. Extremization of the Lagrangian yields the second-order basis set incompleteness correction for the CCSD energy. We also developed a simpler Hylleraas-type functional that only depends on one set of geminal amplitudes by applying screening approximations. This functional is used to develop a diagonal orbital-invariant version of the method in which the geminal amplitudes are fixed at the values determined by the first-order cusp conditions. Extension of the variational method to include perturbatively the effect of connected triples produces the method that approximates the complete basis-set limit of the standard CCSD plus perturbative triples [CCSD(T)] method. For a set of 20 small closed-shell molecules, the method recovered at least 94.5/97.3% of the CBS CCSD(T) correlation energy with the aug-cc-pVDZ/aug-cc-pVTZ orbital basis set. For 12 isogyric reactions involving these molecules, combining the aug-cc-pVTZ correlation energies with the aug-cc-pVQZ Hartree-Fock energies produces the electronic reaction energies with a mean absolute deviation of 1.4 kJ mol(-1) from the experimental values. The method has the same number of optimized parameters as the corresponding CCSD(T) model, does not require any modification of the coupled-cluster computer program, and only needs a small triple-zeta basis to match the precision of the considerably more expensive standard quintuple-zeta CCSD(T) computation.     

Estimating the Hartree—Fock limit from finite basis set calculations

It is demonstrated that the polarization-consistent basis sets, which are optimized for density functional methods, are also suitable for Hartree–Fock calculations, and can be used for estimating the Hartree–Fock basis set limit to within a few micro-hartree accuracy. Various two- and three-point extrapolation schemes are tested and exponential functions are found to be superior compared to functions depending on the inverse power of the highest angular momentum function in the basis set. Total energies can be improved by roughly an order of magnitude, but atomization energies are only marginally improved by extrapolation.     

Adsorption‐induced changes of intramolecular optical transitions: PTCDA/NaCl and PTCDA/KCl

Structural and optical properties of isolated perylene‐3,4,9,10‐tetracarboxylic acid dianhydride molecules adsorbed on (100) oriented NaCl and KCl surfaces were studied theoretically to analyze the recently observed red‐shift of the optical excitation spectrum after adsorption (Müller et al., Phys. Rev. B, 2011, 83, 241203; Paulheim et al. Phys. Chem. Chem. Phys., 2013, 15, 4906). The ground‐state structures were obtained by periodic dispersion‐corrected density functional theory (DFT) calculations. For the excited‐state calculations, nonperiodic time‐dependent DFT methods were applied for a cluster model embedded in point charges. The range‐separated hybrid functional CAM‐B3LYP was used. Correlation‐consistent basis sets were used and the calculated excitation energies were extrapolated to the complete basis set limit. The shift of the first optical excitation energy was analyzed in terms of electronic and geometric contributions. It was found that both the distortion of the molecule due to the interaction with the surface and the electrostatic potential of the surface play an important role. © 2015 Wiley Periodicals, Inc.