Basis set convergence of explicitly correlated double-hybrid density functional theory calculations

The basis set convergence of explicitly correlated double-hybrid density functional theory (DFT) is investigated using the B2GP-PLYP functional. As reference values, we use basis set limit B2GP-PLYP-F12 reaction energies extrapolated from the aug(')-cc-pV(Q+d)Z and aug(')-cc-pV(5+d)Z basis sets. Explicitly correlated double-hybrid DFT calculations converge significantly faster to the basis set limit than conventional calculations done with basis sets saturated up to the same angular momentum (typically, one "gains" one angular momentum in the explicitly correlated calculations). In explicitly correlated F12 calculations the VnZ-F12 basis sets converge faster than the orbital A(')VnZ basis sets. Furthermore, basis set convergence of the MP2-F12 component is apparently faster than that of the underlying Kohn-Sham calculation. Therefore, the most cost-effective approach consists of combining the MP2-F12 correlation energy from a comparatively small basis set such as VDZ-F12 with a DFT energy from a larger basis set such as aug(')-cc-pV(T+d)Z.     

Optimized accurate auxiliary basis sets for RI-MP2 and RI-CC2 calculations for the atoms Rb to Rn

The introduction of the resolution-of-the-identity (RI) approximation for electron repulsion integrals in quantum chemical calculations requires in addition to the orbital basis so-called auxiliary or fitting basis sets. We report here such auxiliary basis sets optimized for second-order Møller–Plesset perturbation theory for the recently published (Weigend and Ahlrichs Phys Chem Chem Phys, 2005, 7, 3297–3305) segmented contracted Gaussian basis sets of split, triple-ζ and quadruple-ζ valence quality for the atoms Rb–Rn (except lanthanides). These basis sets are designed for use in connection with small-core effective core potentials including scalar relativistic corrections. Hereby accurate resolution-of-the-identity calculations with second-order Møller–Plesset perturbation theory (MP2) and related methods can now be performed for molecules containing elements from H to Rn. The error of the RI approximation has been evaluated for a test set of 385 small and medium sized molecules, which represent the common oxidation states of each element, and is compared with the one-electron basis set error, estimated based on highly accurate explicitly correlated MP2–R12 calculations. With the reported auxiliary basis sets the RI error for MP2 correlation energies is typically two orders of magnitude smaller than the one-electron basis set error, independent on the position of the atoms in the periodic table.     

Toward a less costly but accurate calculation of the CCSD(T)/CBS noncovalent interaction energy

The popular method of calculating the noncovalent interaction energies at the coupled-cluster single-, double-, and perturbative triple-excitations [CCSD(T)] theory level in the complete basis set (CBS) limit was to add a CCSD(T) correction term to the CBS second-order Møller-Plesset perturbation theory (MP2). The CCSD(T) correction term is the difference between the CCSD(T) and MP2 interaction energies evaluated in a medium basis set. However, the CCSD(T) calculations with the medium basis sets are still very expensive for systems with more than 30 atoms. Comparatively, the domain-based local pair natural orbital coupled-cluster method [DLPNO-CCSD(T)] can be applied to large systems with over 1,000 atoms. Considering both the computational accuracy and efficiency, in this work, we propose a new scheme to calculate the CCSD(T)/CBS interaction energies. In this scheme, the MP2/CBS term keeps intact and the CCSD(T) correction term is replaced by a DLPNO-CCSD(T) correction term which is the difference between the DLPNO-CCSD(T) and DLPNO-MP2 interaction energies evaluated in a medium basis set. The interaction energies of the noncovalent systems in the S22, HSG, HBC6, NBC10, and S66 databases were recalculated employing this new scheme. The consistent and tight settings of the truncation parameters for DLPNO-CCSD(T) and DLPNO-MP2 in this noncanonical CCSD(T)/CBS calculations lead to the maximum absolute deviation and root-mean-square deviation from the canonical CCSD(T)/CBS interaction energies of less than or equal to 0.28 kcal/mol and 0.09 kcal/mol, respectively. The high accuracy and low cost of this new computational scheme make it an excellent candidate for the study of large noncovalent systems.     

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.     

Auxiliary basis sets for density-fitted correlated wavefunction calculations: weighted core-valence and ECP basis sets for post-d elements

We report optimised auxiliary basis sets for the resolution-of-the-identity (or density-fitting) approximation of two-electron integrals in second-order M?ller-Plesset perturbation theory (MP2) and similar electronic structure calculations with correlation-consistent basis sets for the post-d elements Ga-Kr, In-Xe, and Tl-Rn. The auxiliary basis sets are optimised such that the density-fitting error is negligible compared to the one-electron basis set error. To check to which extent this criterion is fulfilled we estimated for a test set of 80 molecules the basis set limit of the correlation energy at the MP2 level and evaluated the remaining density-fitting and the one-electron basis set errors. The resulting auxiliary basis sets are only 2-6 times larger than the corresponding one-electron basis sets and lead in MP2 calculations to speed-ups of the integral evaluation by one to three orders of magnitude. The density-fitting errors in the correlation energy are at least hundred times smaller than the one-electron basis set error, i.e. in the order of only 1-100 μH per atom.     

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.     

Comparison of ab initio and DFT electronic structure methods for peptides containing an aromatic ring: effect of dispersion and BSSE

We establish that routine B3LYP and MP2 methods give qualitatively wrong conformations for flexible organic systems containing pi systems and that recently developed methods can overcome the known inadequacies of these methods. This is illustrated for a molecule (a conformer of the Tyr-Gly dipeptide) for which B3LYP/6-31+G(d) and MP2/6-31+G(d) geometry optimizations yield strikingly different structures [Mol. Phys. 2006, 104, 559-570]: MP2 predicts a folded "closed-book" conformer with the glycine residue located above the tyrosine ring, whereas B3LYP predicts a more open conformation. By employing different levels of theory, including the local electron correlation methods LMP2 (local MP2) and LCCSD(T0) (local coupled cluster with single, double, and noniterative local triple excitations) and large basis sets (aug-cc-pVnZ, n=D, T, Q), it is shown that the folded MP2 minimum is an artifact caused by large intramolecular BSSE (basis set superposition error) effects in the MP2/6-31+G(d) calculations. The B3LYP functional gives the correct minimum, but the potential energy apparently rises too steeply when the glycine and tyrosine residues approach each other, presumably due to missing dispersion effects in the B3LYP calculations. The PWB6K and M05-2X functionals, designed to give good results for weak interactions, remedy this to some extent. The reduced BSSE in the LMP2 calculations leads to faster convergence with increasing basis set quality, and accurate results can be obtained with smaller basis sets as compared to canonical MP2. We propose LMP2 as a suitable method to study interactions with pi-electron clouds.     

Analysis of the errors in explicitly correlated electronic structure theory

The explicitly correlated second order M?ller-Plesset (MP2-R12) methods perform well in reproducing the last detail of the correlation cusp, allowing higher accuracy than can be accessed through conventional means. Nevertheless in basis sets that are practical for calculations on larger systems (i.e., around triple- or perhaps quadruple-zeta) MP2-R12 fails to bridge the divide between conventional MP2 and the MP2 basis set limit. In this contribution we analyse the sources of error in MP2-R12 calculations in such basis sets. We conclude that the main source of error is the choice of the correlation factor r12. Sources of error that must be avoided for accurate quantum chemistry include the neglect of exchange commutators and the extended Brillouin condition. The generalized Brillouin condition is found not to lead to significant errors.     

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.     

Quantum dynamics of vibrationally activated OH-CO reactant complexes

The MP2 complete basis set (CBS) limit for the binding energy of the two low-lying water octamer isomers of D2d and S4 symmetry is estimated at -72.7+/-0.4 kcal/mol using the family of augmented correlation-consistent orbital basis sets of double through quintuple zeta quality. The largest MP2 calculation with the augmented quintuple zeta (aug-cc-pV5Z) basis set produced binding energies of -73.70 (D2d) and -73.67 kcal/mol (S4). The effects of higher correlation, computed at the CCSD(T) level of theory, are estimated at <0.1 kcal/mol. The newly established MP2/CBS limit for the water octamer is reproduced quite accurately by the newly developed all atom polarizable, flexible interaction potential (TTM2-F). The TTM2-F binding energies of -73.21 (D2d) and -73.24 kcal/mol (S4) for the two isomers are just 0.5 kcal/mol (or 0.7%) larger than the MP2/CBS limit.     

Extrapolation of electron correlation energies to finite and complete basis set targets

The electron correlation energy of two-electron atoms is known to converge asymptotically as approximately (L+1)(-3) to the complete basis set limit, where L is the maximum angular momentum quantum number included in the basis set. Numerical evidence has established a similar asymptotic convergence approximately X(-3) with the cardinal number X of correlation-consistent basis sets cc-pVXZ for coupled cluster singles and doubles (CCSD) and second order perturbation theory (MP2) calculations of molecules. The main focus of this article is to probe for deviations from asymptotic convergence behavior for practical values of X by defining a trial function X(-beta) that for an effective exponent beta=beta(eff)(X,X+1,X+N) provides the correct energy E(X+N), when extrapolating from results for two smaller basis sets, E(X) and E(X+1). This analysis is first applied to "model" expansions available from analytical theory, and then to a large body of finite basis set results (X=D,T,Q,5,6) for 105 molecules containing H, C, N, O, and F, complemented by a smaller set of 14 molecules for which accurate complete basis set limits are available from MP2-R12 and CCSD-R12 calculations. beta(eff) is generally found to vary monotonically with the target of extrapolation, X+N, making results for large but finite basis sets a useful addition to the limited number of cases where complete basis set limits are available. Significant differences in effective convergence behavior are observed between MP2 and CCSD (valence) correlation energies, between hydrogen-rich and hydrogen-free molecules, and, for He, between partial-wave expansions and correlation-consistent basis sets. Deviations from asymptotic convergence behavior tend to get smaller as X increases, but not always monotonically, and are still quite noticeable even for X=5. Finally, correlation contributions to atomization energies (rather than total energies) exhibit a much larger variation of effective convergence behavior, and extrapolations from small basis sets are found to be particularly erratic for molecules containing several electronegative atoms. Observed effects are discussed in the light of results known from analytical theory. A carefully calibrated protocol for extrapolations to the complete basis set limit is presented, based on a single "optimal" exponent beta(opt)(X,X+1,infinity) for the entire set of molecules, and compared to similar approaches reported in the literature.     

A local second-order Møller-Plesset method with localized orbitals: a parallelized efficient electron correlation method

Using orthogonal localized occupied orbitals we have developed and implemented a parallelized local second-order M?ller-Plesset (MP2) method based on the idea developed by Head-Gordon and co-workers. A subset of nonorthogonal correlation functions (the orbital domain) was assigned to each of the localized occupied orbitals using a distance criterion and excitations from localized occupied orbitals that were arranged into subsets. The correlation energy was estimated using a partial diagonalization and an iterative efficient method for solving large-scale linear equations. Some illustrative calculations are provided for molecules with up to 1484 Cartesian basis sets. The orbital domain sizes were found to be independent of the molecular size, and the present local MP2 method covered about 98%-99% of the correlation energy of the conventional canonical MP2 method.     

Perturbation expansion theory corrected from basis set superposition error. I. Locally projected excited orbitals and single excitations

The locally projected self-consistent field molecular orbital method for molecular interaction (LP SCF MI) is reformulated for multifragment systems. For the perturbation expansion, two types of the local excited orbitals are defined; one is fully local in the basis set on a fragment, and the other has to be partially delocalized to the basis sets on the other fragments. The perturbation expansion calculations only within single excitations (LP SE MP2) are tested for water dimer, hydrogen fluoride dimer, and colinear symmetric ArM+ Ar (M = Na and K). The calculated binding energies of LP SE MP2 are all close to the corresponding counterpoise corrected SCF binding energy. By adding the single excitations, the deficiency in LP SCF MI is thus removed. The results suggest that the exclusion of the charge-transfer effects in LP SCF MI might indeed be the cause of the underestimation for the binding energy.     

Explicitly correlated RMP2 for high-spin open-shell reference states

We present an explicitly correlated version of the high-spin open-shell RMP2 method. The theory is derived in a unitarily invariant form, which is suitable for the insertion of local approximations. It is demonstrated that the rapid basis set convergence of closed-shell MP2-F12 is also achieved in RMP2-F12, and similar Ansatze and approximations can be employed. All integrals are computed using efficient density fitting approximations, and many-electron integrals are avoided using resolution of the identity approximations. The performance of the method is demonstrated by benchmark calculations on a large set of ionization potentials, electron affinities and atomization energies. Using triple-zeta basis sets RMP2-F12 yields results that are closer to the basis set limit than standard RMP2 with augmented quintuple-zeta basis sets for all properties. Different variants of perturbative corrections for the open-shell Hartree-Fock treatment are described and tested.     

Treating dispersion effects in extended systems by hybrid MP2:DFT calculations--protonation of isobutene in zeolite ferrierite

We propose use of a hybrid method to study problems that involve both bond rearrangements and van-der-Waals interactions. The method combines second-order M?ller-Plesset perturbation theory (MP2) calculations for the reaction site with density functional theory (DFT) calculations for a large system under periodic boundary conditions. Hybrid MP2:DFT structure optimisation for a cluster embedded in the periodic model is the first of three steps in a multi-level approach. The second step is extrapolation of the MP2 energy to the complete basis set limit. The third step is extrapolating the high-level (MP2) correction to the limiting case of the full periodic structure. This is done by calculating the MP2 correction for a series of cluster models of increasing size, fitting an analytic expression to these energy corrections, and applying the fitted expression to the full periodic structure. We assume that, up to a constant, the high-level correction is described by a damped dispersion expression. Combining the results of all three steps yields an estimate of the MP2 reaction energy for the full periodic system at the complete basis set level. The method is designed for a reaction between a small or medium sized substrate molecule and a very large chemical system. For adsorption of isobutene in zeolite H-ferrierite, the energies obtained for the formation of different structures, the pi-complex, the isobutoxide, the tert-butoxide, and the tert-butyl carbenium ion, are -78, -73, -48, and -21 kJ mol(-1), respectively. This corresponds to corrections of the pure DFT (PBE functional) results by -62, -70, -67, and -29 kJ mol(-1), respectively. Hence, the MP2 corrections are substantial and, perhaps more importantly, not the same for the different hydrocarbon species in the zeolite. Coupled-cluster (CCSD(T)) calculations change the MP2 energies by -4 kJ mol(-1) (tert-butyl cation) or less (below +/-1 kJ mol(-1) for the other species).     

Polarization-consistent versus correlation-consistent basis sets in predicting molecular and spectroscopic properties

Compared to the correlation-consistent basis sets, it is not known if polarization-consistent pc-n basis sets, which were initially developed for HF and DFT calculations, can provide a monotonic and faster convergence toward the basis-set limit for results at correlated levels as well as better accuracy for a similar number of basis functions. It is also not known whether the pc-n basis sets can compute second derivatives of energy, such as nuclear magnetic shielding tensors, efficiently. To address these questions, the pc-n (n = 1-4), cc-pVxZ, and/or aug-cc-pVxZ (x = D, T, Q, 5, and 6) basis sets were used to compute the molecular and/or spectroscopic parameters of H2, H2O, and NH3 at the RHF, B3-LYP, MP2, and/or CCSD(T) levels of theory. The results show that compared to the cc-pVxZ and/or aug-cc-pVxZ basis sets the pc-n basis sets yield faster convergence toward the basis-set limit but equivalent molecular and/or spectroscopic parameters in the basis-set limit at the RHF, DFT, MP2, and CCSD(T) levels. Because the pc-n basis sets show faster convergence, fewer basis-set functions are needed to reach the accuracy obtained with the aug-cc-pVxZ basis sets, enabling faster calculations and less computer storage space. The results also show that the pc-n basis sets, in conjunction with the "locally dense" basis-set approach, could be applied to predict accurate parameters; thus, they could be used to estimate accurate molecular or spectroscopic properties (e.g., NMR parameters) for larger systems such as the active site of enzymes.     

On geometries of stacked and H-bonded nucleic acid base pairs determined at various DFT, MP2, and CCSD(T) levels up to the CCSD(T)/complete basis set limit level

The geometries and interaction energies of stacked and hydrogen-bonded uracil dimers and a stacked adeninecdots, three dots, centeredthymine pair were studied by means of high-level quantum chemical calculations. Specifically, standard as well as counterpoise-corrected optimizations were performed at second-order Moller-Plesset (MP2) and coupled cluster level of theory with single, double, and perturbative triple excitations [CCSD(T)] levels with various basis sets up to the complete basis set limit. The results can be summarized as follows: (i) standard geometry optimization with small basis set (e.g., 6-31G(*)) provides fairly reasonable intermolecular separation; (ii) geometry optimization with extended basis sets at the MP2 level underestimates the intermolecular distances compared to the reference CCSD(T) results, whereas the MP2/cc-pVTZ counterpoise-corrected optimization agrees well with the reference geometries and, therefore, is recommended as a next step for improving MP2/cc-pVTZ geometries; (iii) the stabilization energy of stacked nucleic acids base pairs depends considerably on the method used for geometry optimization, so the use of reliable geometries, such as counterpoise-corrected MP2/cc-pVTZ ones, is recommended; (iv) the density functional theory methods fail completely in locating the energy minima for stacked structures and when the geometries from MP2 calculations are used, the resulting stabilization energies are strongly underestimated; (v) the self-consistent charges-density functional tight binding method, with inclusion of the empirical dispersion energy, accurately reproduces interaction energies and geometries of dispersion-bonded (stacked) complexes; this method can thus be recommended for prescanning the potential energy surfaces of van der Waals complexes.     

Can extrapolation to the basis set limit be an alternative to the counterpoise correction? A study on the helium dimer

Configuration interaction and coupled cluster calculations are reported for He₂ using various orbital basis sets of the d-aug-AVXZ type, with the results being extrapolated to the one electron basis set limit both with counterpoise and without counterpoise correction. A generalized uniform singlet- and triplet-pair extrapolation scheme has been utilized for such a purpose. Using appropriate corrections to mimic full configuration interaction, the energies were predicted in excellent agreement with the best available estimates. The results also suggest that extrapolation to the complete basis set limit may be a general alternative to the counterpoise correction that yields a more accurate potential energy while being more economical.     

Is spin-component scaled second-order Møller-Plesset perturbation theory an appropriate method for the study of noncovalent interactions in molecules?

Testing of the spin-component scaled second-order M?ller-Plesset (SCS-MP2) method for the computation of noncovalent interaction energies is done with a database of 165 biologically relevant complexes. The effects of the spin-scaling procedure (i.e., MP2 vs SCS-MP2), the basis set size, and the corrections for basis set superposition error (BSSE) are systematically examined. When using two-point basis set extrapolations for the correlation energy, augmentation of the atomic orbital basis with computationally costly diffuse functions is found to be obsolete. In general, SCS-MP2 also improves results for noncovalent interactions statistically on MP2, and significant outliers are removed. Moreover, it is shown that effects of BSSE and one-particle basis set incompleteness almost cancel each other in the case of triple-zeta sets (SCS-MP2/TZVPP or SCS-MP2/cc-pVTZ without counterpoise correction), which opens a practical route to efficient computations for large systems. We recommend SCS-MP2 as the preferred quantum chemical wave function based method for the noncovalent interactions in large biologically relevant systems when reasonable coupled-cluster with single and double and perturbative triple excitations (CCSD(T)) calculations cannot be performed anymore. A comparison to MP2 and CCSD(T) interaction energies for n-alkane dimers, however, indicates (and this also holds to a lesser extent for hydrogen-bonded systems) limitations of SCS-MP2 when treating chemically "saturated" interactions. The different behavior of second-order perturbation theory for saturated and for stacked pi-systems is discussed.