X

Search Results

Searching....

6.7 Attenuated MP2

6.7.1 Overview

(July 14, 2022)

MP2(attenuator, basis) approximates MP2 by splitting the Coulomb operator in two pieces and preserving only short-range two-electron interactions, akin to the CASE approximation, 302 Dombroski J. P., Taylor S. W., Gill P. M. W.
J. Phys. Chem.
(1996), 100, pp. 6272.
Link
, 23 Adamson R. D., Dombroski J. P., Gill P. M. W.
Chem. Phys. Lett.
(1996), 254, pp. 329.
Link
but without modification of the underlying SCF calculation. While MP2 is a comparatively efficient method for estimating the correlation energy, it converges slowly with basis set size — and, even in the complete basis limit, contains fundamentally inaccurate physics for long-range interactions. Basis set superposition error and the MP2-level treatment of long-range interactions both typically artificially increase correlation energies for non-covalent interactions. Attenuated MP2 improves upon MP2 for inter- and intramolecular interactions, with significantly better performance for relative and binding energies of non-covalent complexes, frequently outperforming complete basis set estimates of MP2. 402 Goldey M., Head-Gordon M.
J. Phys. Chem. Lett.
(2012), 3, pp. 3592.
Link
, 401 Goldey M., Dutoi A., Head-Gordon M.
Phys. Chem. Chem. Phys.
(2013), 15, pp. 15869.
Link

Attenuated MP2, denoted MP2(attenuator, basis) is implemented in Q-Chem based on the complementary terf function, below:

s(r)=terfc(r,r0)=12{erfc[ω(r-r0)]+erfc[ω(r+r0)]} (6.26)

By choosing the terfc short-range operator, we optimally preserve the short-range behavior of the Coulomb operator while smoothly and rapidly switching off around the distance r0. Since this directly addresses basis set superposition error, parameterization must be done for specific basis sets. This has been performed for the basis sets, aug-cc-pVDZ 402 Goldey M., Head-Gordon M.
J. Phys. Chem. Lett.
(2012), 3, pp. 3592.
Link
and aug-cc-pVTZ. 401 Goldey M., Dutoi A., Head-Gordon M.
Phys. Chem. Chem. Phys.
(2013), 15, pp. 15869.
Link
Other basis sets are not recommended for general use until further testing has been done.

Energies and gradients are functional with and without the resolution of the identity approximation using correlation keywords ATTMP2 and ATTRIMP2.

Example 6.20  Example of RI-MP2(terfc, aug-cc-pVDZ) energy evaluation

$molecule
   0 1
   O    -1.551007   -0.114520    0.000000
   H    -1.934259    0.762503    0.000000
   H    -0.599677    0.040712    0.000000
$end

$rem
   METHOD          attrimp2
   BASIS           aug-cc-pvdz
   AUX_BASIS       rimp2-aug-cc-pvdz
$end

Example 6.21  Example of MP2(terfc, aug-cc-pVTZ) geometry optimization

$molecule
   0 1
   H    0.0   0.0   0.0
   H    0.0   0.0   0.9
$end

$rem
   JOBTYPE         opt
   METHOD          attmp2
   BASIS           aug-cc-pvtz
$end