Q-Chem 5.1 User’s Manual

6.4 Exact MP2 Methods

6.4.1 Algorithm

Second-order Møller-Plesset theory[Møller and Plesset(1934)] (MP2) is probably the simplest useful wave function-based electron correlation method. Revived in the mid-1970s, it remains highly popular today, because it offers systematic improvement in optimized geometries and other molecular properties relative to Hartree-Fock (HF) theory.[Hehre et al.(1986)Hehre, Radom, v. R. Schleyer, and Pople] Indeed, in a recent comparative study of small closed-shell molecules,[Helgaker et al.(1997)Helgaker, Gauss, Jorgensen, and Olsen] MP2 outperformed much more expensive singles and doubles coupled-cluster theory for such properties! Relative to state-of-the-art Kohn-Sham density functional theory (DFT) methods, which are the most economical methods to account for electron correlation effects, MP2 has the advantage of properly incorporating long-range dispersion forces. The principal weaknesses of MP2 theory are for open shell systems, and other cases where the HF determinant is a poor starting point.

Q-Chem contains an efficient conventional semi-direct method to evaluate the MP2 energy and gradient.[Head-Gordon(1999)] These methods require $OVN$ memory ($O$, $V$, $N$ are the numbers of occupied, virtual and total orbitals, respectively), and disk space which is bounded from above by $OVN^{2}/2$. The latter can be reduced to $IVN^2/2$ by treating the occupied orbitals in batches of size $I$, and re-evaluating the two-electron integrals $O/I$ times. This approach is tractable on modern workstations for energy and gradient calculations of at least 500 basis functions or so, or molecules of between 15 and 30 first row atoms, depending on the basis set size. The computational cost increases between the 3rd and 5th power of the size of the molecule, depending on which part of the calculation is time-dominant.

The algorithm and implementation in Q-Chem is improved over earlier methods,[Head-Gordon et al.(1988)Head-Gordon, Pople, and Frisch, Frisch et al.(1990)Frisch, Head-Gordon, and Pople] particularly in the following areas:

The implementation offers the user the following alternatives:

The semi-direct algorithm is the only choice for gradient calculations. It is also normally the most efficient choice for energy calculations. There are two classes of exceptions:

There are three important options that should be wisely chosen by the user in order to exploit the full efficiency of Q-Chem’s direct and semi-direct MP2 methods (as discussed above, the LOCAL_OCCUPIED method has different requirements).

6.4.2 Algorithm Control and Customization

The direct and semi-direct integral transformation algorithms used by Q-Chem (e.g., MP2, CIS(D)) are limited by available disk space, $D$, and memory, $C$, the number of basis functions, $N$, the number of virtual orbitals, $V$ and the number of occupied orbitals, $O$, as discussed above. The generic description of the key $rem variables are:

MEM_STATIC

Sets the memory for Fortran AO integral calculation and transformation modules.


TYPE:

INTEGER


DEFAULT:

64

corresponding to 64 Mb.


OPTIONS:

$n$

User-defined number of megabytes.


RECOMMENDATION:

For direct and semi-direct MP2 calculations, this must exceed OVN + requirements for AO integral evaluation (32–160 Mb), as discussed above.


MEM_TOTAL

Sets the total memory available to Q-Chem, in megabytes.


TYPE:

INTEGER


DEFAULT:

2000

Corresponding to 2000 Mb.


OPTIONS:

$n$

User-defined number of megabytes.


RECOMMENDATION:

Use the default, or set equal to the physical memory of your machine. Note that if the memory allocation total more than 1 Gb for a CCMAN job, the memory is allocated as follows

12%

MEM_STATIC

50%

CC_MEMORY

35%

Other memory requirements:


AO2MO_DISK

Sets the amount of disk space (in megabytes) available for MP2 calculations.


TYPE:

INTEGER


DEFAULT:

2000

Corresponding to 2000 Mb.


OPTIONS:

$n$

User-defined number of megabytes.


RECOMMENDATION:

Should be set as large as possible, discussed in Section 6.4.1.


CD_ALGORITHM

Determines the algorithm for MP2 integral transformations.


TYPE:

STRING


DEFAULT:

Program determined.


OPTIONS:

DIRECT

Uses fully direct algorithm (energies only).

SEMI_DIRECT

Uses disk-based semi-direct algorithm.

LOCAL_OCCUPIED

Alternative energy algorithm (see 6.4.1).


RECOMMENDATION:

Semi-direct is usually most efficient, and will normally be chosen by default.


6.4.3 Example

Example 6.69  Example of an MP2/6-31G* calculation employing the frozen core approximation. Note that the EXCHANGE $rem variable will default to HF

$molecule
   0 1
   O
   H1  O  oh
   H2  O  oh  H1  hoh

   oh  = 1.01
   hoh = 105
$end

$rem
   METHOD          mp2
   BASIS           6-31g*
   N_FROZEN_CORE   fc
$end