6.6 Auxiliary Basis (Resolution of the Identity) MP2 Methods

6.6.6 Brueckner CC2

(May 16, 2021)

Brueckner orbitals (BOs) are shown to largely ameliorate the artificial symmetry breaking that occurs at the Hartree-Fock (HF) level. 632 Lee J., Head-Gordon M.
Phys. Chem. Chem. Phys.
(2019), 21, pp. 47638.
Link
This can lead to improved results even in higher order wavefunction theories that typically use HF orbitals as a starting point. Unfortunately, the cheapest traditional Brueckner theory, Brueckner doubles (BD), is still quite computationally demanding. Orbital optimized MP2 (OOMP2) was proposed as a low scaling approximation to BD 690 Lochan R. C., Head-Gordon M.
J. Chem. Phys.
(2007), 126, pp. 164101.
Link
due to the similarity of orbital optimized coupled cluster doubles (OOCCD) and BD. Another possible source of approximate BD orbitals is Bruckner CC2 (BCC2). 31 Akinaga Y., Kawashima Y., Ten-no S.
Chem. Phys. Lett.
(2011), 506, pp. 276.
Link

Rather than optimizing the energy as in OOMP2, BCC2 optimizes the orbitals in order to reduce the CC2 t1 amplitudes to 0. In the absence of t1 amplitudes, the CC2 doubles amplitudes are exactly the same as MP2, and the singles amplitude equations are as follows:

Ωia=fai+jbfjbtijab+12jbcja||cbtijbc+12jkbjk||bitjkab (6.24)

BCC2 has a computational cost scaling as 𝒪(o2v3) - a significant improvement over OOCCD/BD, and can greatly improve the quality of orbitals over HF in open-shell molecules due to the reduction of artificial spin contamination.

Example 6.15  Example of BCC2 methodology

$molecule
   1 2
   F
   H 1 1.001
$end

$rem
   UNRESTRICTED      TRUE
   JOBTYPE           SP
   EXCHANGE          HF
   GEN_SCFMAN_FINAL  TRUE
   DO_BCC2           3                run BCC2
   SCF_ALGORITHM     DIIS
   SCF_GUESS         GWH
   BASIS             sto-3g
   AUX_BASIS         rimp2-vdz
   SCF_CONVERGENCE   8
   THRESH            14
   SYMMETRY          FALSE
   PURECART          1111
$end

View output

As CC2 is a perturbative approach similar to MP2, breakdowns occur in the presence of orbital degeneracy. While BCC2 does not explicitly optimize the energy and thereby drive orbitals toward degeneracy, results are still severely hindered when orbital energy differences become small. The κ regularizer originally proposed by Joonho Lee and Martin Head-Gordon 631 Lee J., Head-Gordon M.
J. Chem. Theory Comput.
(2018), 14, pp. 5203.
Link
can be used in to further improve the results of BCC2 in this case. The regularizer is applied by modifying the t amplitudes as follows:

tijab=-ab||ijΔijab(1-exp(-κΔijab))2 (6.25)

As in κ-OOMP2, the parameter κ sets the regularization strength, with κ=0 yielding HF, and κ= yielding unregularized BCC2. A value of κ=1.2 Eh-1 is suggested for most applications. κ-BCC2 runs through libgscf and libgmbpt. Currently, DIIS should be used to converge BCC2 orbitals, as the singles residual is not an orbital gradient and cannot be used with gradient-based algorithms. The BCC2 code can currently handle restricted (R) and unrestricted (U) orbital types.

Summary of rem variables relevant to run κ-BCC2:

CORRELATION None (default)
JOBTYPE sp (default) single point energy evaluation
BASIS user’s choice (standard or user-defined: GENERAL or MIXED)
GEN_SCFMAN_FINAL TRUE
SCF_ALGORITHM DIIS (gradient based algorithms currently unsupported)
AUX_BASIS corresponding auxiliary basis (standard or user-defined:
AUX_GENERAL or AUX_MIXED)
DO_BCC2 3 (run BCC2)
REGULARIZED_BCC2 0 (no regularizer; default)
1 (δ-regularizer)
2 (κ-regularizer; recommended)
3 (σ-regularizer)
REG_PARAMETER regularization parameter multiplied by 1e3; no default
1200 (Recommended value for κ-BCC2)
N_FROZEN_CORE 0 (frozen core currently unsupported)
N_FROZEN_VIRTUAL 0 (frozen core currently unsupported)
DO_S2 0 (default)
1 (Compute S2 at the MP2 level)

Example 6.16  Example of κ-BCC2 with the R orbital type applied to a water dimer

$molecule
   0  1
   O    -2.766559046    0.187082886    0.566917837
   H    -3.696304300    1.179189102   -0.642506882
   H    -3.395837846   -1.509891173    0.389283582
   O     2.587035064    0.275900014   -0.746441819
   H     3.579141280    0.918406897    0.633058252
   H     0.852266482    0.311804811   -0.156847268
$end

$rem
   EXCHANGE          hf
   BASIS             cc-pvdz
   AUX_BASIS_CORR    rimp2-cc-pvdz
   THRESH            14
   INPUT_BOHR        true
   SCF_CONVERGENCE   8
   MAXSCF            1000
   SCF_GUESS         sad
   SYMMETRY          false
   GEN_SCFMAN_FINAL  true
   UNRESTRICTED      false    use restricted
   DO_BCC2           3        run BCC2
   REGULARIZED_O2    2        use kappa-regularizer
   REG_VARIABLE      1450     set kappa = 1.45
$end

View output

Example 6.17  Example of κ-BCC2 with the U orbital type applied to an OH radical

$molecule
0  2
O -2.766559046 0.187082886 0.566917837
H -3.6963043 1.179189102 -0.642506882
$end

$rem
   EXCHANGE          hf
   BASIS             cc-pvdz
   AUX_BASIS_CORR    rimp2-cc-pvdz
   THRESH            14
   INPUT_BOHR        true
   SCF_CONVERGENCE   8
   MAXSCF            1000
   SCF_GUESS         sad
   SYMMETRY          false
   GEN_SCFMAN_FINAL  true
   UNRESTRICTED      true     use unrestricted
   DO_BCC2           3        run BCC2
   REGULARIZED_O2    2        use kappa-regularizer
   REG_VARIABLE      1450     set kappa = 1.45
$end

View output