## Abstract

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.

Original language | English (US) |
---|---|

Article number | 221106 |

Journal | Journal of Chemical Physics |

Volume | 127 |

Issue number | 22 |

DOIs | |

State | Published - 2007 |

## All Science Journal Classification (ASJC) codes

- Physics and Astronomy(all)
- Physical and Theoretical Chemistry

## Access to Document

## Other files and links

## Fingerprint

Dive into the research topics of 'A simple and efficient CCSD(T)-F12 approximation'. Together they form a unique fingerprint.## Cite this

- APA
- Author
- BIBTEX
- Harvard
- Standard
- RIS
- Vancouver

*Journal of Chemical Physics*,

*127*(22), [221106]. https://doi.org/10.1063/1.2817618

**A simple and efficient CCSD(T)-F12 approximation**. In: Journal of Chemical Physics. 2007 ; Vol. 127, No. 22.

}

*Journal of Chemical Physics*, vol. 127, no. 22, 221106. https://doi.org/10.1063/1.2817618

**A simple and efficient CCSD(T)-F12 approximation.**/ Adler, Thomas B.; Knizia, Gerald; Werner, Hans Joachim.

In: Journal of Chemical Physics, Vol. 127, No. 22, 221106, 2007.

Research output: Contribution to journal › Article › peer-review

TY - JOUR

T1 - A simple and efficient CCSD(T)-F12 approximation

AU - Adler, Thomas B.

AU - Knizia, Gerald

AU - Werner, Hans Joachim

N1 - Funding Information: Adler Thomas B. Knizia Gerald Werner Hans-Joachim a) Institut für Theoretische Chemie, Universität Stuttgart , Pfaffenwaldring 55, D-70569 Stuttgart, Germany a) Author to whom correspondence should be addressed. Electronic mail: werner@theochem.uni-stuttgart.de 14 12 2007 127 22 221106 19 10 2007 01 11 2007 14 12 2007 2007-12-14T15:41:38 2007 American Institute of Physics 0021-9606/2007/127(22)/221106/4/ $23.00 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. DFG The painfully slow convergence of the electron correlation energy with respect to the basis set size can be dramatically improved by including terms into the wave function that depend explicitly on the interelectronic distances. 1,2 In recent years, much progress has been made in improving the accuracy and efficiency of such explicitly correlated methods 3–18 and in making them applicable to rather large molecules. 15–17 Most of this work was restricted to second-order Møller-Plesset perturbation theory, but recently also practical explicitly correlated coupled-cluster (CC) methods have been developed. Based on earlier work of Noga et al. , 19,20 Fliegl et al. 21,22 proposed a CCSD(T)(F12) model, which differs from the exact theory by the neglect of terms that are quadratic in the explicitly correlated amplitudes. This year, Tew et al. 23 demonstrated that with augmented triple-zeta basis sets this method yields correlation energies and reaction energies that are as accurate as standard calculations with augmented quintuple-zeta basis sets. They stated that such calculations are about three times as expensive as conventional CCSD(T) calculations. In this communication we present a method that is even more efficient and accurate. The general ansatz for closed-shell CCSD-F12 wave functions with single and double excitations, Ψ CCSD = e T ̂ 1 + T ̂ 2 Ψ HF , (1) employs excitation operators T ̂ 1 = t a i E ̂ a i , (2) T ̂ 2 = T a b i j E ̂ a i E ̂ b j + T α β i j E ̂ α i E ̂ β j . (3) Here and in the following, summation over repeated dummy indices is implied. The indices i , j , … refer to occupied orbitals a , b to virtual (external) orbitals, and α , β to a complete orbital basis set. t a i and T a b i j are, respectively, the conventional singles and doubles amplitudes. In F12 theory, the additional amplitudes T α β i j are approximated as T α β i j = ⟨ α β ∣ Q ̂ 12 F ̂ 12 ∣ k l ⟩ T k l i j , (4) where the projector Q ̂ 12 = 1 − ∣ r s ⟩ ⟨ r s ∣ − ∣ m x ⟩ ⟨ m x ∣ − ∣ x m ⟩ ⟨ x m ∣ (5) ensures strong orthogonality of the explicitly correlated terms to the Hartree-Fock reference function as well as to the conventional double excitations. 3,6,17,24–26 The indices r , s denote the full orbital (MO) basis, and x the complementary auxiliary (CA) orbital basis. 6,17 Note that T α β i j = 0 if both α and β belong to the orbital basis, or if at least one of them corresponds to an occupied orbital. As first proposed by Ten-no 11 and subsequently used by several authors, 7,10,12–14,16,17 the correlation factor F ̂ 12 is taken to be a simple Slater function, F ( r 12 ) = exp ( − β r 12 ) . (6) In practice, we approximate the Slater function by a linear combination of Gaussians. 7,17 The above ansatz is orbital invariant. It suffers, however, from geminal basis set superposition errors, 17,27 and for larger molecules also numerical instabilities occur. These problems can be avoided by using the “diagonal” ansatz, 17,25 in which only the terms for k l = i j and k l = j i are kept. This ansatz is not orbital invariant and only size consistent when localized orbitals are used. With localized orbitals it yields more accurate results than the more general orbital invariant ansatz. 17 Even simpler is the diagonal fixed amplitudeansatz by Ten-no, 11 in which T i i i i = t s , T i j i j = 1 2 [ t s + t t ] , and T j i i j = 1 2 [ t s − t t ] ( i ≠ j ) . All other amplitudes T k l i j are zero. The fixed parameters t s = − 1 ∕ ( 2 β ) and t t = − 1 ∕ ( 4 β ) follow from the cusp conditions 28,29 for singlet and triplet pairs, respectively. This ansatz, which is orbital invariant, size consistent, and free of geminal superposition errors, has been used in the current work. For simplicity, we will in the following still write T k l i j , but it will always be assumed that k l is restricted as defined above. The CCSD-F12 residual equations for the doubles amplitudes T a b i j can be straightforwardly derived by using the ansatz in Eq. (3) . However, this requires resolution of the identity (RI) approximations in many places and does not lead to a numerically stable method. In addition, if the amplitudes T k l i j are fully optimized, an extended set of CC equations has to be solved. In order to overcome these problems, Fliegl et al. 21 have introduced various approximations; in particular, they kept only linear terms in T k l i j . In the current work, we propose even simpler approximations, which turn out to work equally well. First, if the diagonal fixed amplitude ansatz is used, no equations need to be solved for the amplitudes T k l i j . It remains, however, the question of how to compute the energy. For fully optimized amplitudes the CCSD energy is simply E corr = K a b i j [ 2 T ¯ a b i j − T ¯ b a i j ] + V k l i j [ 2 T k l i j − T l k i j ] , (7) with T ¯ a b i j = T a b i j + t a i t b j , K a b i j = ⟨ i j ∣ r 12 − 1 ∣ a b ⟩ , and V k l i j = ⟨ i j ∣ r 12 − 1 Q ̂ 12 F ̂ 12 ∣ k l ⟩ . This expression is linear in the amplitudes, and therefore errors of the amplitudes affect the energy in first order. This is critical if the fixed amplitude ansatz is applied in which the amplitudes T k l i j are not optimized. It turns out, however, that the error is very small if the energy is optimized with respect to the exponent β of the correlation factor. In other words, when β is optimized, the optimized amplitudes T i j i j and T j i i j are very close to the fixed values from the cusp conditions, and Eq. (7) should give accurate results. This is demonstrated in Fig. 1 for the H 2 O molecule. Unfortunately, the optimum value of β varies for different molecules and basis sets. For the aug-cc-pVTZ (AVTZ) basis the optimal values lie roughly between 1.1 a 0 − 1 for pure hydrocarbons and 1.5 a 0 − 1 for molecules containing mainly oxygen atoms. In order to reduce the dependence on β , we have chosen to replace the second term in Eq. (7) by the MP2-F12 Hylleraas functional, i.e., V k l i j is replaced by V k l i j + R k l i j , where R k l i j is the residual for the explicit part of the MP2-F12 equations. In the case of MP2-F12, this energy expression is accurate to second order in errors of the amplitudes T k l i j and, therefore, much less sensitive to the choice of β and T k l i j . This can also be seen in Fig. 1 . In our CCSD(T)-F12a approximation the residuals R k l i j are computed using the complementary auxiliary basis set approach 6 and approximation 3C. This approximation is based on the MP2-R12/C method by Kedžuch et al. 30 and our implementation is described in detail in Ref. 17 . The only difference to Ref. 17 is that the CCSD amplitudes T a b i j replace the MP2 amplitudes in the term C a b k l T a b i j . The second approximation concerns the contributions of the explicitly correlated terms to the CC equations for the singles and doubles amplitudes. The simplest approximation would be to neglect these entirely and to compute the correlation energy as sum of the CCSD(T) correlation energy and the MP2-F12 correction (i.e., the difference of the pure MP2-F12 and MP2 energies). It has been noted that, usually, the basis set dependence of MP2 is stronger than that of CCSD(T), and that, therefore, such a correction has limited accuracy. Nevertheless, already this approximation works quite well (see supplementary material 35 ). A better approximation can be made by including the most important explicitly correlated contributions into the CCSD residuals. The conventional CCSD residuals R a b i j can be written in matrix form as R CC i j = R MP 2 i j + K ( D i j ) + α i j , k l T ¯ k l + G i j + G j i † , (8) where R MP 2 i j is the MP2 residual, K ( D i j ) are external exchange operators which account for the contributions of the three- and four-external integrals (see Ref. 31 for definitions and details 39 ). The explicit form of the quantities α i j , k l and G i j involve products of integral and amplitude matrices. In our method we assume that the explicitly correlated terms contribute only to the terms that are linear in the amplitude matrices. Furthermore, we make a “hybrid approximation” 5,17 in which we neglect all contributions of the CA basis set { x } . This leads to the correction Δ R F 12 a i j = [ W k l − K ( F k l ) + C k l ] T k l i j , (9) where W a b k l = ⟨ k l ∣ F ̂ 12 r 12 − 1 ∣ a b ⟩ and [ K ( F k l ) ] a b = ⟨ a b ∣ r 12 − 1 ∣ r s ⟩ ⟨ r s ∣ F ̂ 12 ∣ k l ⟩ . (10) The matrices C k l arise from R MP 2 i j ; thus, they are the same as in MP2-F12 theory. 17 The remaining two terms arise from the explicitly correlated contributions to the external exchange operators K ( D i j ) , where we have approximated the projector in Eq. (5) as Q ̂ 12 = 1 − ∣ r s ⟩ ⟨ r s ∣ . As in theCCSD(T)(F12) method of Tew et al. , 23 no direct contributions of the explicitly correlated terms to the triples are included. A similar term as in Eq. (9) occurs in the CCSD-F12 residual R k l i j . If this is included and the energy contribution of R k l i j is still computed as [ 2 T k l i j − T l k i j ] R k l i j , one obtains an additional energy correction, Δ E = ( 2 T k l i j − T l k i j ) tr [ ( W k l − K ( F k l ) ) D j i ] . (11) We denote the approximation that includes this correction as CCSD(T)-F12b. This contribution roughly doubles the effect of the coupling of the conventional and explicit parts. It is not a priori clear whether this will improve the accuracy, since many further terms in the R k l i j residual are neglected, and the resulting energy functional is an approximation by itself. 40 In fact, as will be shown below, it leads to a worsening of the convergence of the correlation energies with basis set size. The contributions in Eqs. (9) and (10) can be computed with negligible additional effort together with the operators K ( D i j ) as K ( D i j − F ̃ i j ) , where F ̃ i j = F k l T k l i j . The internal-external parts of the modified external exchange operators also contribute to the singles residuals. They are evaluated in the AO basis, 31 and therefore no additional transformed integrals are required. Thus, as compared to a conventional CCSD(T) calculation, the only additional effort is to perform a MP2-F12 calculation beforehand. All required integrals for the MP2-F12 calculation are evaluated using efficient density-fitting (DF) approximations, as described in previous works. 4,7,16,17 In the Hartree-Fock and CCSD(T) calculations no density-fitting approximations were used. In order to test the method we have studied the same 15 reactions as Tew et al. 23 These involve 23 molecules with first- and second-row atoms. Only the valence electrons were correlated and MP2/aug-cc-pVTZ geometries were used. Reference values were computed using the standard CCSD(T) and augmented correlation consistent basis sets (aug-cc-pVnZ). 32 The results up to quintuple zeta are identical to the values reported by Tew et al. In addition, for the smaller molecules we were also able to use aug-cc-pV6Z basis sets. For the second-row atoms, the aug - cc - pV ( n + d ) Z basis sets 33 with one additional tight d shell have been used. For simplicity, in the following we will denote these basis sets AVnZ. The basis set limits were estimated by two-point fits of the CCSD(T) correlation energies to the extrapolation formula 34 E n = E CBS + A n − 3 . Results obtained with n = 4 , 5 and 5, 6 will be denoted CBS[45] and CBS[56], respectively. The convergence of the valence-shell correlation energy for H 2 O is shown in Table I for the AVnZ basis sets and β = 1.3 a 0 − 1 . The convergence of the CCSD(T)-F12a method is found to be faster than that of the CCSD(T)(F12) approximation of Tew et al. 23 and more similar to the basis set convergence of MP2-F12. The CCSD(T)-F12a AVQZ and AV5Z values are close to the CBS[45] estimate ( 308.97 m E h ) but slightly overestimate the CBS[56] value ( 308.34 m E h ) . On the other hand, the CCSD(T)-F12b approximation appears to underestimate the correlation energies. Using the AVTZ basis the CCSD(T)-F12a and CCSD(T)-F12b approximations yield 99.5% and 98.1% of the CBS[56] basis set limit. On average, 99.2% and 97.6%, respectively, of the CBS[45] values were obtained for the 23 molecules . The mean and maximum absolute errors of the Hartree-Fock (HF) and correlation energy contributions to the reaction energies are presented in Table II (for individual values see supplementary material 35 ). As HF reference values we took the AV5Z results. It is found that for the AVTZ basis set the mean and maximum HF errors are significantly larger than those of the CCSD(T)-F12a correlation contributions. However, the HF basis set error can be almost entirely removed by adding an MP2 singles energy correction Δ E s = 2 t α i f i α , obtained by allowing single excitations into the CA orbital space. This means solving f α i = t α k f k i − f α β t β i , where α , β run over the virtual and CA orbitals (this is entirely distinct from the CCSD calculation). If this correction is applied, the accuracy of the HF contribution is again comparable to that of the correlation energy contributions. For the correlation contributions, the CBS[56] values have been taken as reference whenever available, otherwise the CBS[45] ones (using the CBS[56] values slightly reduces the errors). The parameter β has been chosen to be 1.0 a 0 − 1 . In the current tests this value yielded best results, but the dependence of the errors on the choice of β was found to be rather weak on the average. 35 The results in Table II demonstrate that our CCSD(T)-F12a approximation yields even better accuracy than the CCSD(T)(F12) method of Tew et al. , 23 and the reaction energies obtained with the AVTZ basis are better than quintuple-zeta quality. The mean absolute error of the total CCSD(T)-F12a reaction energies (including HF and Δ E s ) amounts to 0.4 kJ ∕ mol , the maximum error to 1.0 kJ ∕ mol . For the CCSD(T)-F12b approximation the corresponding values are 0.6 and 1.4 kJ ∕ mol , respectively. Finally, we mention that the MP2-F12 calculation, which precedes the CCSD(T) one, takes only a relatively small amount of time. For the largest case in the present study, namely, the C H 3 C H O molecule, the total CPU time (using MOLPRO , 36 including integral evaluation and HF) was 3856 s (AMD Opteron 2.6 GHz ; no symmetry has been used). Of this time, the DF-MP2-F12 took 334 sec , i.e., less than 10%. Thus, for almost no additional cost, results of at least CCSD(T)/AV5Z quality are obtained. The corresponding standard CCSD(T)/AV5Z calculations would require about 100 times more time and disk space. This work was funded in the priority program 1145 by the Deutsche Forschungsgemeinschaft. One of the authors (H.J.W.) is also grateful for generous support by the Fonds der Chemischen Industrie. We thank Fred Manby for providing the integral routines used in the DF-MP2-F12 program. Table I. Convergence of the valence-shell CCSD(T) correlation energy of H 2 O as a function of the basis set (in − m E h ) Basis Orbital (F12) a F12a b F12b b AVDZ 232.71 296.09 299.58 293.15 AVTZ 282.04 304.56 306.82 302.42 AVQZ 297.93 307.12 308.78 306.20 AV5Z 303.32 307.81 308.93 307.40 AV6Z 305.43 a CCSD(T)(F12) values from Ref. 23 , β = 1.3 a 0 − 1 . b β = 1.3 a 0 − 1 . The V5Z/JKFIT (Ref. 37 ) and AV5Z/MP2FIT (Ref. 38 ) bases have been used for RI and DF, respectively. Table II. Mean absolute (MAD), root mean square (RMS), and maximum absolute (MAX) errors of Hartree-Fock and CCSD(T) correlation energy contributions to the reaction energies in kJ/mol. The AVTZ basis set was used for all F12 calculations. Error HF CCSD(T) AVTZ AVTZ a AV5Z (F12) b F12a c F12b c MAD 1.77 0.28 0.58 0.68 0.38 0.54 RMS 2.75 0.36 0.86 0.90 0.55 0.61 MAX 6.42 0.86 2.15 1.62 1.69 1.01 a Including singles corrections, see text. b CCSD(T)(F12) values from Ref. 23 , β = 1.3 a 0 − 1 . c β = 1.0 a 0 − 1 . The VQZ/JKFIT (Ref. 37 ) and AVQZ/MP2FIT (Ref. 38 ) bases have been used for RI and DF, respectively. FIG. 1. Dependence of the MP2-F12 energy correction of H 2 O on the choice of β (upper panel) and the fixed amplitudes t S = − t / β and t T = − t ∕ ( 2 β ) for β = 1.5 a 0 − 1 (lower panel); basis AVTZ. The full lines represent the energy computed by the Hylleraas functional, the dashed lines (open symbols) the first-order approximation (see text).

PY - 2007

Y1 - 2007

N2 - 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.

AB - 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.

UR - http://www.scopus.com/inward/record.url?scp=37149036190&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=37149036190&partnerID=8YFLogxK

U2 - 10.1063/1.2817618

DO - 10.1063/1.2817618

M3 - Article

AN - SCOPUS:37149036190

SN - 0021-9606

VL - 127

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

IS - 22

M1 - 221106

ER -