Coupled Cluster Methods¶

Code author: T. Daniel Crawford

Section author: T. Daniel Crawford

Module: Keywords, PSI Variables, CCENERGY

Module: Keywords, PSI Variables, CCEOM

Module: Keywords, PSI Variables, CCRESPONSE

Module: Keywords, PSI Variables, CCTRIPLES

Module: Keywords, CCDENSITY

Module: Keywords, CCHBAR

Module: Keywords, CCLAMBDA

Module: Keywords, CCSORT

The coupled cluster approach is one of the most accurate and reliable quantum chemical techniques for including the effects of electron correlation. Instead of the linear expansion of the wavefunction used by configuation interaction, coupled cluster uses an exponential expansion,

(1) $| \Psi \rangle &= e^{\hat{T}} | \Phi_0 \rangle \\ &= \left( 1 + {\hat{T}} + \frac{1}{2} {\hat{T}}^2 + \frac{1}{3!}{\hat{T}}^3 + \cdots \right) | \Phi_0 \rangle,$

where the cluster operator ${\hat{T}}$ is written as a sum of operators that generate singly-excited, doubly-excited, etc., determinants:

${\hat{T}} = {\hat{T}_1} + {\hat{T}_2} + {\hat{T}_3} + \cdots + {\hat{T}_N},$

with

${\hat T}_1 | \Phi_0 \rangle &= \sum_{i}^{\rm occ} \sum_a^{\rm vir} t_i^a | \Phi_i^a \rangle \\ {\hat T}_2 | \Phi_0 \rangle &= \sum_{i<j}^{\rm occ} \sum_{a<b}^{\rm vir} t_{ij}^{ab} | \Phi_{ij}^{ab} \rangle,$

etc. The popular coupled cluster singles and doubles (CCSD) model [Purvis:1982] truncates the expansion at ${\hat{T}} = {\hat{T}_1} + {\hat{T}_2}$ . This model has the same number of parameters as configuration interaction singles and doubles (CISD) but improves upon it by approximately accounting for higher-order terms using products of lower-order terms (e.g., the term ${\hat{T}_2}^2$ approximately accounts for quadruple excitations). The inclusion of such products makes coupled-cluster methods size extensive, meaning that the quality of the computation should not degrade for larger molecules. The computational cost for CCSD scales as ${\cal{O}}(o^2 v^4)$ , where $o$ is the number of occupied orbitals and $v$ is the number of virtual orbitals.

Improving upon CCSD, the CCSD(T) method [Raghavachari:1989] includes a perturbative estimate of the energy contributed by the ${\hat{T}_3}$ operator. The computational cost of this additional term scales as ${\cal{O}}(o^3 v^4)$ , making it rather expensive for molecules with more than a dozen heavy atoms or so. However, when this method is affordable, it provides very high quality results in most cases.

PSI4 is capable of computing energies and analytic gradients for a number of coupled cluster models. It can also compute linear response properties (such as static or frequency-dependent polarizability, or optical rotation angles) for some models. Excited states can also be computed by the CC2 and CC3 models, or by EOM-CCSD. Table CC Methods summarizes these capabilities. This section describes how to carry out coupled cluster calculations within PSI4. For higher-order coupled-cluster methods like CCSDT and CCSDTQ, PSI4 can interface to Kállay’s MRCC code (see MRCC).

Current coupled cluster capabilities of PSI4
Reference	Method	Energy	Gradient	Exc. Energies	LR Props
RHF	CC2	Y	N	Y	Y
UHF	CC2	Y	N	Y	N
ROHF	CC2	Y	N	Y	N
RHF	CCSD	Y	Y	Y	Y
RHF	CCSD(T)	Y	N	—	—
ROHF	CCSD	Y	Y	Y	N
ROHF	CCSD(T)	Y	N	—	—
UHF	CCSD	Y	Y	Y	N
UHF	CCSD(T)	Y	Y	—	—
RHF	CC3	Y	N	Y	N
UHF	CC3	Y	N	Y	N
ROHF	CC3	Y	N	Y	N
Brueckner	CCD	Y	N	N	N
Brueckner	CCD(T)	Y	N	—	—

The following wavefunctions are currently recognized by PSI4 as arguments to functions like energy(): 'ccsd', 'ccsd(t)', 'cc2', 'cc3', 'bccd' (CCD with Brueckner orbitals), 'bccd(t)' (CCD(T) with Brueckner orbitals), 'eom-ccsd', 'eom-cc2' (CC2 for excited states), 'eom-cc3' (CC3 for excited states). Response properties can be obtained by calling the function property() (instead of, for example, energy(), e.g., property('ccsd'). There are many sample coupled cluster inputs provided in psi4/samples.

Basic Keywords¶

A complete list of keywords related to coupled-cluster computations is provided in the appendices, with the majority of the relevant keywords appearing in Appendix CCENERGY. For a standard ground-state CCSD or CCSD(T) computation, the following keywords are common:

REFERENCE¶

Reference wavefunction type

Type: string

Possible Values: RHF, ROHF, UHF

Default: RHF

R_CONVERGENCE¶

Convergence criterion for wavefunction (change) in CC amplitude equations.

Type: conv double

Default: 1e-7

MAXITER¶

Maximum number of iterations to solve the CC equations

Type: integer

Default: 50

BRUECKNER_ORBS_R_CONVERGENCE¶

Convergence criterion for Breuckner orbitals. The convergence is determined based on the largest $T_1$ amplitude.

Type: conv double

Default: 1e-5

RESTART¶

Do restart the coupled-cluster iterations from old $t_1$ and $t_2$ amplitudes? For geometry optimizations, Brueckner calculations, etc. the iterative solution of the CC amplitude equations may benefit considerably by reusing old vectors as initial guesses. Assuming that the MO phases remain the same between updates, the CC codes will, by default, re-use old vectors, unless the user sets RESTART = false.

Type: boolean

Default: true

CACHELEVEL¶

Cacheing level for libdpd governing the storage of amplitudes, integrals, and intermediates in the CC procedure. A value of 0 retains no quantities in cache, while a level of 6 attempts to store all quantities in cache. For particularly large calculations, a value of 0 may help with certain types of memory problems. The default is 2, which means that all four-index quantites with up to two virtual-orbital indices (e.g., $\langle ij | ab \rangle>$ integrals) may be held in the cache.

Type: integer

Default: 2

CACHETYPE¶

Selects the priority type for maintaining the automatic memory cache used by the libdpd codes. A value of LOW selects a “low priority” scheme in which the deletion of items from the cache is based on pre-programmed priorities. A value of LRU selects a “least recently used” scheme in which the oldest item in the cache will be the first one deleted.

Type: string

Possible Values: LOW, LRU

Default: LOW

NUM_AMPS_PRINT¶

Number of important $t_1$ and $t_2$ amplitudes to print

Type: integer

Default: 10

MP2_AMPS_PRINT¶

Do print the MP2 amplitudes which are the starting guesses for RHF and UHF reference functions?

Type: boolean

Default: false

Larger Calculations¶

Here are a few recommendations for carrying out large-basis-set coupled cluster calculations with PSI4:

In most cases it is reasonable to set the memory keyword to 90% of the available physical memory, at most. There is a small amount of overhead associated with the coupled cluster modules that is not accounted for by the internal CC memory handling routines. Thus, the user should not specify the entire physical memory of the system, or swapping is likely. However, for especially large calculations, it is better to set the memory keyword to a value less than 16 GB.
Set the CACHELEVEL keyword to 0. This will turn off cacheing, which, for very large calculations, can lead to heap fragmentation and memory faults, even when sufficient physical memory exists.
Set the PRINT keyword to 2. This will help narrow where memory bottlenecks or other errors exist in the event of a crash.

Excited State Coupled Cluster Calculations¶

A complete list of keywords related to coupled cluster linear response is provided in Appendix CCEOM. The most important keywords associated with EOM-CC calculations are:

ROOTS_PER_IRREP¶

Number of excited states per irreducible representation for EOM-CC and CC-LR calculations. Irreps denote the final state symmetry, not the symmetry of the transition.

Type: array

Default: No Default

E_CONVERGENCE¶

Convergence criterion for excitation energy (change) in the Davidson algorithm for CC-EOM.

Type: conv double

Default: 1e-8

SINGLES_PRINT¶

Do print information on the iterative solution to the single-excitation EOM-CC problem used as a guess to full EOM-CC?

Type: boolean

Default: false

SCHMIDT_ADD_RESIDUAL_TOLERANCE¶

Minimum absolute value above which a guess vector to a root is added to the Davidson algorithm in the EOM-CC iterative procedure.

Type: conv double

Default: 1e-3

EOM_GUESS¶

Specifies a set of single-excitation guess vectors for the EOM-CC procedure. If EOM_GUESS = SINGLES, the guess will be taken from the singles-singles block of the similarity-transformed Hamiltonian, Hbar. If EOM_GUESS = DISK, guess vectors from a previous computation will be read from disk. If EOM_GUESS = INPUT, guess vectors will be specified in user input. The latter method is not currently available.

Type: string

Possible Values: SINGLES, DISK, INPUT

Default: SINGLES

Linear Response (CCLR) Calculations¶

Linear response computations are invoked like property('ccsd') or property('cc2'), along with a list of requested properties. A complete list of keywords related to coupled cluster linear response is provided in Appendix CCRESPONSE.

The most important keywords associated with CC-LR calculations are as follows.

PROPERTY¶

The response property desired. Acceptable values are POLARIZABILITY (default) for dipole-polarizabilities, ROTATION for specific rotations, ROA for Raman Optical Activity, and ALL for all of the above.

Type: string

Possible Values: POLARIZABILITY, ROTATION, ROA, ALL

Default: POLARIZABILITY

OMEGA¶

Array that specifies the desired frequencies of the incident radiation field in CCLR calculations. If only one element is given, the units will be assumed to be atomic units. If more than one element is given, then the units must be specified as the final element of the array. Acceptable units are HZ, NM, EV, and AU.

Type: array

Default: No Default

GAUGE¶

Specifies the choice of representation of the electric dipole operator. Acceptable values are LENGTH for the usual length-gauge representation, VELOCITY for the modified velocity-gauge representation in which the static-limit optical rotation tensor is subtracted from the frequency- dependent tensor, or BOTH. Note that, for optical rotation calculations, only the choices of VELOCITY or BOTH will yield origin-independent results.

Type: string

Possible Values: LENGTH, VELOCITY, BOTH

Default: LENGTH

Coupled Cluster Methods¶

Basic Keywords¶

REFERENCE¶

R_CONVERGENCE¶

MAXITER¶

BRUECKNER_ORBS_R_CONVERGENCE¶

RESTART¶

CACHELEVEL¶

CACHETYPE¶

NUM_AMPS_PRINT¶

MP2_AMPS_PRINT¶

Larger Calculations¶

Excited State Coupled Cluster Calculations¶

ROOTS_PER_IRREP¶

E_CONVERGENCE¶

SINGLES_PRINT¶

SCHMIDT_ADD_RESIDUAL_TOLERANCE¶

EOM_GUESS¶

Linear Response (CCLR) Calculations¶

PROPERTY¶

OMEGA¶

GAUGE¶

Table Of Contents

Previous topic

Next topic

This Page

Navigation

Coupled Cluster Methods¶

Basic Keywords¶

Larger Calculations¶

Excited State Coupled Cluster Calculations¶

Linear Response (CCLR) Calculations¶

Previous topic

Next topic

This Page

Quick search

Navigation