Q-Chem 4.4 User’s Manual

5.8 Noniterative Corrections to Coupled Cluster Energies

5.8.1 (T) Triples Corrections

To approach chemical accuracy in reaction energies and related properties, it is necessary to account for electron correlation effects that involve three electrons simultaneously, as represented by triple substitutions relative to the mean field single determinant reference, which arise in MP4. The best standard methods for including triple substitutions are the CCSD(T) [300] and QCISD(T) methods [291] The accuracy of these methods is well-documented for many cases [301], and in general is a very significant improvement relative to the starting point (either CCSD or QCISD). The cost of these corrections scales with the 7th power of molecule size (or the 4th power of the number of basis functions, for a fixed molecule size), although no additional disk resources are required relative to the starting coupled-cluster calculation. Q-Chem supports the evaluation of CCSD(T) and QCISD(T) energies, as well as the corresponding OD(T) correction to the optimized doubles method discussed in the previous subsection. Gradients and properties are not yet available for any of these (T) corrections.

5.8.2 (2) Triples and Quadruples Corrections

While the (T) corrections discussed above have been extraordinarily successful, there is nonetheless still room for further improvements in accuracy, for at least some important classes of problems. They contain judiciously chosen terms from 4th- and 5th-order Møller-Plesset perturbation theory, as well as higher order terms that result from the fact that the converged cluster amplitudes are employed to evaluate the 4th- and 5th-order order terms. The (T) correction therefore depends upon the bare reference orbitals and orbital energies, and in this way its effectiveness still depends on the quality of the reference determinant. Since we are correcting a coupled-cluster solution rather than a single determinant, this is an aspect of the (T) corrections that can be improved. Deficiencies of the (T) corrections show up computationally in cases where there are near-degeneracies between orbitals, such as stretched bonds, some transition states, open shell radicals, and diradicals.

Prof. Steve Gwaltney, while working at Berkeley with Martin Head-Gordon, has suggested a new class of non iterative correction that offers the prospect of improved accuracy in problem cases of the types identified above [302]. Q-Chem contains Gwaltney’s implementation of this new method, for energies only. The new correction is a true second order correction to a coupled-cluster starting point, and is therefore denoted as (2). It is available for two of the cluster methods discussed above, as OD(2) and CCSD(2) [302, 303]. Only energies are available at present.

The basis of the (2) method is to partition not the regular Hamiltonian into perturbed and unperturbed parts, but rather to partition a similarity-transformed Hamiltonian, defined as $\tilde{H}=e^{-\hat{T}}\hat{H}e^{\hat{T}}$. In the truncated space (call it the $p$-space) within which the cluster problem is solved (e.g., singles and doubles for CCSD), the coupled-cluster wave function is a true eigenvalue of $\tilde{H}$. Therefore we take the zero order Hamiltonian, $\tilde{H}^{(0)}$, to be the full $\tilde{H}$ in the p-space, while in the space of excluded substitutions (the q-space) we take only the one-body part of $ \tilde{H}$ (which can be made diagonal). The fluctuation potential describing electron correlations in the $q$-space is $\tilde{H}-\tilde{H}^{(0)}$, and the (2) correction then follows from second order perturbation theory.

The new partitioning of terms between the perturbed and unperturbed Hamiltonians inherent in the (2) correction leads to a correction that shows both similarities and differences relative to the existing (T) corrections. There are two types of higher correlations that enter at second order: not only triple substitutions, but also quadruple substitutions. The quadruples are treated with a factorization ansatz, that is exact in 5th order Møller-Plesset theory [304], to reduce their computational cost from $N^{9}$ to $N^{6}$. For large basis sets this can still be larger than the cost of the triples terms, which scale as the 7th power of molecule size, with a factor twice as large as the usual (T) corrections.

These corrections are feasible for molecules containing between four and ten first row atoms, depending on computer resources, and the size of the basis set chosen. There is early evidence that the (2) corrections are superior to the (T) corrections for highly correlated systems [302]. This shows up in improved potential curves, particularly at long range and may also extend to improved energetic and structural properties at equilibrium in problematical cases. It will be some time before sufficient testing on the new (2) corrections has been done to permit a general assessment of the performance of these methods. However, they are clearly very promising, and for this reason they are available in Q-Chem.

5.8.3 (dT) and (fT) corrections

Alternative inclusion of noniterative $N^7$ triples corrections is described in Section 6.7.20. These methods called (dT) and (fT) are of similar accuracy to other triples corrections. CCSD(dT) and CCSD(fT) are equivalent to the CR-CCSD(T)$_ L$ and CR-CCSD(T)$_2$ methods of Piecuch and co-workers [305]

5.8.4 Job Control Options

The evaluation of a non iterative (T) or (2) correction after a coupled-cluster singles and doubles level calculation (either CCSD, QCISD or OD) is controlled by the correlation keyword, and the specification of any frozen orbitals via N_FROZEN_CORE (and possibly N_FROZEN_VIRTUAL).

There is only one additional job control option. For the (2) correction, it is possible to apply the frozen core approximation in the reference coupled cluster calculation, and then correlate all orbitals in the (2) correction. This is controlled by CC_INCL_CORE_CORR, described below.

The default is to include core and core-valence correlation automatically in the CCSD(2) or OD(2) correction, if the reference CCSD or OD calculation was performed with frozen core orbitals. The reason for this choice is that core correlation is economical to include via this method (the main cost increase is only linear in the number of core orbitals), and such effects are important to account for in accurate calculations. This option should be made false if a job with explicitly frozen core orbitals is desired. One good reason for freezing core orbitals in the correction is if the basis set is physically inappropriate for describing core correlation (e.g., standard Pople basis sets, and Dunning cc-pV$x$Z basis sets are designed to describe valence-only correlation effects). Another good reason is if a direct comparison is desired against another method such as CCSD(T) which is always used in the same orbital window as the CCSD reference.

CC_INCL_CORE_CORR

Whether to include the correlation contribution from frozen core orbitals in non iterative (2) corrections, such as OD(2) and CCSD(2).


TYPE:

LOGICAL


DEFAULT:

TRUE


OPTIONS:

TRUE

FALSE


RECOMMENDATION:

Use the default unless no core-valence or core correlation is desired (e.g., for comparison with other methods or because the basis used cannot describe core correlation).


5.8.5 Example

Example 5.86  Two jobs that compare the correlation energy calculated via the standard CCSD(T) method with the new CCSD(2) approximation, both using the frozen core approximation. This requires that CC_INCL_CORE_CORR must be specified as FALSE in the CCSD(2) input.

$molecule
   0  2
   O
   H  O  0.97907
$end

$rem
   METHOD          ccsd(t)
   BASIS           cc-pvtz
   N_FROZEN_CORE   fc
$end

@@@

$molecule
   read
$end

$rem
   METHOD              ccsd(2)
   BASIS               cc-pvtz
   N_FROZEN_CORE       fc
   CC_INCL_CORE_CORR   false
$end

Example 5.87  Water: Ground state CCSD(dT) calculation using RI

$molecule
0 1
O
H1 O OH
H2 O OH H1 HOH

OH  = 0.957
HOH = 104.5
$end

$rem
JOBTYPE       SP
BASIS         cc-pvtz
AUX_BASIS     rimp2-cc-pvtz
METHOD        CCSD(dT)
$end