Q-Chem 5.1 User’s Manual

6.5 Local MP2 Methods

6.5.1 Local Triatomics in Molecules (TRIM) Model

The development of what may be called “fast methods” for evaluating electron correlation is a problem of both fundamental and practical importance, because of the unphysical increases in computational complexity with molecular size which afflict “exact” implementations of electron correlation methods. Ideally, the development of fast methods for treating electron correlation should not impact either model errors or numerical errors associated with the original electron correlation models. Unfortunately this is not possible at present, as may be appreciated from the following rough argument. Spatial locality is what permits re-formulations of electronic structure methods that yield the same answer as traditional methods, but faster. The one-particle density matrix decays exponentially with a rate that relates to the HOMO-LUMO gap in periodic systems. When length scales longer than this characteristic decay length are examined, sparsity will emerge in both the one-particle density matrix and also pair correlation amplitudes expressed in terms of localized functions. Very roughly, such a length scale is about 5 to 10 atoms in a line, for good insulators such as alkanes. Hence sparsity emerges beyond this number of atoms in 1-D, beyond this number of atoms squared in 2-D, and this number of atoms cubed in 3-D. Thus for three-dimensional systems, locality only begins to emerge for systems of between hundreds and thousands of atoms.

If we wish to accelerate calculations on systems below this size regime, we must therefore introduce additional errors into the calculation, either as numerical noise through looser tolerances, or by modifying the theoretical model, or perhaps both. Q-Chem’s approach to local electron correlation is based on modifying the theoretical models describing correlation with an additional well-defined local approximation. We do not attempt to accelerate the calculations by introducing more numerical error because of the difficulties of controlling the error as a function of molecule size, and the difficulty of achieving reproducible significant results. From this perspective, local correlation becomes an integral part of specifying the electron correlation treatment. This means that the considerations necessary for a correlation treatment to qualify as a well-defined theoretical model chemistry apply equally to local correlation modeling. The local approximations should be

To ensure that these model chemistry criteria are met, Q-Chem’s local MP2 methods[Lee et al.(2000)Lee, Maslen, and Head-Gordon, Head-Gordon et al.(1999a)Head-Gordon, Lee, and Maslen] express the double substitutions (i.e., the pair correlations) in a redundant basis of atom-labeled functions. The advantage of doing this is that local models satisfying model chemistry criteria can be defined by performing an atomic truncation of the double substitutions. A general substitution in this representation will then involve the replacement of occupied functions associated with two given atoms by empty (or virtual) functions on two other atoms, coupling together four different atoms. We can force one occupied to virtual substitution (of the two that comprise a double substitution) to occur only between functions on the same atom, so that only three different atoms are involved in the double substitution. This defines the triatomics in molecules (TRIM) local model for double substitutions. The TRIM model offers the potential for reducing the computational requirements of exact MP2 theory by a factor proportional to the number of atoms. We could also force each occupied to virtual substitution to be on a given atom, thereby defining a more drastic diatomics in molecules (DIM) local correlation model.

The simplest atom-centered basis that is capable of spanning the occupied space is a minimal basis of core and valence atomic orbitals on each atom. Such a basis is necessarily redundant because it also contains sufficient flexibility to describe the empty valence anti-bonding orbitals necessary to correctly account for non-dynamical electron correlation effects such as bond-breaking. This redundancy is actually important for the success of the atomic truncations because occupied functions on adjacent atoms to some extent describe the same part of the occupied space. The minimal functions we use to span the occupied space are obtained at the end of a large basis set calculation, and are called extracted polarized atomic orbitals (EPAOs).[Lee and Head-Gordon(2000b)] We discuss them briefly below. It is even possible to explicitly perform an SCF calculation in terms of a molecule-optimized minimal basis of polarized atomic orbitals (PAOs) (see Chapter 4). To span the virtual space, we use the full set of atomic orbitals, appropriately projected into the virtual space.

We summarize the situation. The number of functions spanning the occupied subspace will be the minimal basis set dimension, $M$, which is greater than the number of occupied orbitals, $O$, by a factor of up to about two. The virtual space is spanned by the set of projected atomic orbitals whose number is the atomic orbital basis set size $N$, which is fractionally greater than the number of virtuals $VNO$. The number of double substitutions in such a redundant representation will be typically three to five times larger than the usual total. This will be more than compensated by reducing the number of retained substitutions by a factor of the number of atoms, $A$, in the local triatomics in molecules model, or a factor of $A^2$ in the diatomics in molecules model.

The local MP2 energy in the TRIM and DIM models are given by the following expressions, which can be compared against the full MP2 expression given earlier in Eq. (6.13). First, for the DIM model:

  \begin{equation} \label{eq516} E_{\ensuremath{\mathrm{DIM\, MP2}}} = -\frac{1}{2}\sum _{\bar{P},\bar{Q}} \frac{ (\bar{P} | \bar{Q}) (\bar{P} || \bar{Q}) }{ \Delta _{\bar{P}} + \Delta _{\bar{Q}} } \end{equation}   (6.17)

The sums run over the linear number of atomic single excitations after they have been canonicalized. Each term in the denominator is thus an energy difference between occupied and virtual levels in this local basis. Similarly, the TRIM model corresponds to the following local MP2 energy:

  \begin{equation} \label{eq517} E_{\ensuremath{\mathrm{TRIM\, MP2}}} = -\sum _{\bar{P},jb} \frac{ (\bar{P} | jb) (\bar{P} || jb) }{ \Delta _{\bar{P}} + \varepsilon _ b - \varepsilon _ j } - E_{\ensuremath{\mathrm{DIM\, MP2}}} \end{equation}   (6.18)

where the sum is now mixed between atomic substitutions $\bar{P}$, and non-local occupied $j$ to virtual $b$ substitutions. See Refs. Lee:2000b,Head-Gordon:1999c for a full derivation and discussion.

The accuracy of the local TRIM and DIM models has been tested in a series of calculations.[Lee et al.(2000)Lee, Maslen, and Head-Gordon, Head-Gordon et al.(1999a)Head-Gordon, Lee, and Maslen] In particular, the TRIM model has been shown to be quite faithful to full MP2 theory via the following tests:

6.5.2 EPAO Evaluation Options

When a local MP2 job (requested by the LOCAL_MP2 option for CORRELATION) is performed, the first new step after the SCF calculation is converged is to extract a minimal basis of polarized atomic orbitals (EPAOs) that spans the occupied space. There are three valid choices for this basis, controlled by the PAO_METHOD and EPAO_ITERATE keywords described below.

PAO_METHOD

Controls the type of PAO calculations requested.


TYPE:

STRING


DEFAULT:

EPAO

For local MP2, EPAOs are chosen by default.


OPTIONS:

EPAO

Find EPAOs by minimizing delocalization function.

PAO

Do SCF in a molecule-optimized minimal basis.


RECOMMENDATION:

None


EPAO_ITERATE

Controls iterations for EPAO calculations (see PAO_METHOD).


TYPE:

INTEGER


DEFAULT:

0

Use non-iterated EPAOs based on atomic blocks of SPS.


OPTIONS:

$n$

Optimize the EPAOs for up to $n$ iterations.


RECOMMENDATION:

Use the default. For molecules that are not too large, one can test the sensitivity of the results to the type of minimal functions by the use of optimized EPAOs in which case a value of $n=500$ is reasonable.


EPAO_WEIGHTS

Controls algorithm and weights for EPAO calculations (see PAO_METHOD).


TYPE:

INTEGER


DEFAULT:

115

Standard weights, use 1$^{\ensuremath{\mathrm{st}}}$ and 2$^{\ensuremath{\mathrm{nd}}}$ order optimization


OPTIONS:

15 

Standard weights, with 1$^{\ensuremath{\mathrm{st}}}$ order optimization only.


RECOMMENDATION:

Use the default, unless convergence failure is encountered.


6.5.3 Algorithm Control and Customization

A local MP2 calculation (requested by the LOCAL_MP2 option for CORRELATION) consists of the following steps:

Details of the efficient implementation of the local MP2 method described above are reported in the recent thesis of Dr. Michael Lee.[Lee(2000)] Here we simply summarize the capabilities of the program. The computational advantage associated with these local MP2 methods varies depending upon the size of molecule and the basis set. As a rough general estimate, TRIM MP2 calculations are feasible on molecule sizes about twice as large as those for which conventional MP2 calculations are feasible on a given computer, and this is their primary advantage. Our implementation is well suited for large basis set calculations. The AO basis two-electron integrals are evaluated four times. DIM MP2 calculations are performed as a by-product of TRIM MP2 but no separately optimized DIM algorithm has been implemented.

The resource requirements for local MP2 calculations are as follows:

The evaluation of the local MP2 energy does not require any further customization. An adequate amount of MEM_STATIC (80 to 160 Mb) should be specified to permit efficient AO basis two-electron integral evaluation, but all large scratch arrays are allocated from MEM_TOTAL.

6.5.4 Examples

Example 6.70  A relative energy evaluation using the local TRIM model for MP2 with the 6-311G** basis set. The energy difference is the internal rotation barrier in propenal, with the first geometry being planar trans, and the second the transition structure.

$molecule
   0 1
   C
   C  1  1.32095
   C  2  1.47845  1  121.19  
   O  3  1.18974  2  123.83  1   180.00
   H  1  1.07686  2  121.50  3     0.00
   H  1  1.07450  2  122.09  3   180.00
   H  2  1.07549  1  122.34  3   180.00
   H  3  1.09486  2  115.27  4   180.00
$end

$rem
   METHOD        local_mp2
   BASIS         6-311g**
$end

@@@

$molecule
   0 1
   C
   C  1  1.31656 
   C  2  1.49838  1  123.44 
   O  3  1.18747  2  123.81  1   92.28
   H  1  1.07631  2  122.03  3   -0.31
   H  1  1.07484  2  121.43  3  180.28
   H  2  1.07813  1  120.96  3  180.34
   H  3  1.09387  2  115.87  4  179.07
$end

$rem
   CORRELATION   local_mp2
   BASIS         6-311g**
$end