11.6 Projector-Based Density Embedding

11.6.2 Job Control for DFT-in-DFT and WFT-in-DFT Calculations

(June 30, 2021)

The DFT-in-DFT and WFT-in-DFT calculations are implemented in Q-Chem under the framework of its default SCF engine GEN_SCFMAN, which is triggered by setting GEN_SCFMAN_EMBED = TRUE. The embedding job requires specification of fragments through the $molecule section (see Sec. 12.2). Unless otherwise specified, the first fragment is regarded as the embedded (“active”) fragment. The user is allowed to select the embedded fragment through the $embed_fragment section.

The high-level electronic structure method is given by $rem variable METHOD, and the low-level method can be specified via ENV_METHOD. Note that one should choose a mean-field method as the low-level theory, such as pure or hybrid density functionals. Unlike the previous implementation of this embedding theory based on “EmbedMan”, there is no specific limitation on supported high-level methods.

The calculation starts with an SCF calculation at the low-level theory, and then it performs an embedded SCF calculation with the higher-level DFT (for DFT-in-DFT cases) or Hartree-Fock (for WFT-in-DFT cases). To make post-SCF correlated methods directly applicable, the molecular orbitals are stored on disk in the “[Occ(environment) | Occ(active) | Vir (active) | Vir(inactive)]” order, and the environment occupied orbitals and the excluded inactive virtual orbitals are treated as frozen core and frozen virtual orbitals in the correlated WFT method, respectively. Note that the virtual orbitals are not truncated unless the concentric localization method is requested by setting CONCENTRIC_VIRTS = TRUE since that is the only virtual space truncation scheme for projector-based embedding theory that is available in Q-Chem to date.

In Q-Chem 5.3, the use of projector-based embedding theory is limited to overall closed-shell systems (“restricted-in-restricted” embedding). Extension to open-shell systems will be enabled in future releases.

GEN_SCFMAN_EMBED
       Run a projector-based embedding calculation using the implementation based onGEN_SCFMAN
TYPE:
       BOOLEAN
DEFAULT:
       FALSE
OPTIONS:
       TRUE Perform a projector-based embedding calculation FALSE Do not perform an embedding calculation
RECOMMENDATION:
       None

ENV_METHOD
       Specify the low-level theory in a projector-based embedding calculation
TYPE:
       STRING
DEFAULT:
       NONE
OPTIONS:
       Parsed in the same way as $rem variable “METHOD
RECOMMENDATION:
       A mean-field method (pure or hybrid density functional) should be chosen.

SPADE_PARTITION
       Use the SPADE approach to determine the initial set of embedded (active) orbitals
TYPE:
       BOOLEAN
DEFAULT:
       FALSE
OPTIONS:
       TRUE Use SPADE to partition the occupied space FALSE Use the Pipek-Mezey localization + Mulliken population to assign occupied orbitals
RECOMMENDATION:
       Use SPADE if a significant gap in the spectrum of singular values can be detected.

CONCENTRIC_VIRTS
       Use the concentric localization (CL) scheme to truncate the virtual space
TYPE:
       BOOLEAN
DEFAULT:
       FALSE
OPTIONS:
       TRUE Use the CL scheme to truncate the virtual space FALSE Leave the virtual space untruncated
RECOMMENDATION:
       Use CL truncation for WFT-in-DFT calculations.

CONCENTRIC_VIRTS_ZETA
       Specify the size of the truncated virtual space
TYPE:
       INTEGER
DEFAULT:
       2
OPTIONS:
       m The total number of the CL-truncated virtuals is m×noccactive
RECOMMENDATION:
       Use the default; set it to a larger value if higher accuracy is requested.

CONCENTRIC_REF_BASIS
       Specify the projection basis (PB) in the concentric localization procedure
TYPE:
       STRING
DEFAULT:
       NONE
OPTIONS:
       Parsed in the same way as BASIS; if unspecified, the working basis (WB) will be used as PB.
RECOMMENDATION:
       WB is usually a good choice; a smaller basis can chosen with caution to further reduce the computational cost.

Example 11.32  DFT-in-DFT embedding calculation: ωB97X-V (water 2) in the presence of water molecules 1 and 3 described by PBE. The embedded (“active”) fragment is specified by user through the $embed_fragment section.

$molecule
   0 1
   --
   0 1
   O    -0.00000   -0.00000   0.00000
   H     0.56704    0.41003   0.65488
   H     0.16427   -0.95671   0.05856
   --
   0 1
   O    -0.45701   -2.75857   0.16997
   H    -0.36918   -3.38806  -0.54777
   H    -1.36971   -2.42607   0.12617
   --
   0 1
   O    -2.62421   -0.99554  -0.00046
   H    -3.25789   -0.73194   0.66864
   H    -1.89698   -0.35254   0.04416
$end

$rem
   METHOD             wb97x-v
   BASIS              def2-svpd
   GEN_SCFMAN_EMBED   true
   THRESH             14
   SCF_CONVERGENCE    8
   ENV_METHOD         pbe
$end

$embed_fragment
   2
$end

View output

Example 11.33  RI-MP2 (fragment 1) embedded in the environment described by B3LYP. The SPADE scheme is employed to partition the occupied space and the concentric localization (CL) method to truncate virtual space. The working basis cc-pVDZ is used as the default projection basis for CL. The same setup works with other correlated methods such as MP2, CCSD, etc.

$molecule
   0 1
   ---
   -1 1
   O       -1.1867 -0.2472  0.0000
   H       -1.9237  0.3850  0.0000
   ---
   1 1
   H       -0.0227  1.1812  0.8852
   C        0.0000  0.5526  0.0000
   H       -0.0227  1.1812 -0.8852
   C        1.1879 -0.3829  0.0000
   H        2.0985  0.2306  0.0000
   H        1.1184 -1.0093  0.8869
   H        1.1184 -1.0093 -0.8869
$end

$rem
   METHOD             rimp2
   BASIS              cc-pvdz
   AUX_BASIS          rimp2-cc-pvdz
   GEN_SCFMAN_EMBED   true
   THRESH             14
   SCF_CONVERGENCE    8
   ENV_METHOD         b3lyp
   SPADE_PARTITION    true
   CONCENTRIC_VIRTS   true
   N_FROZEN_CORE      0   !no extra FC orbitals
   SYMMETRY           false
   SYM_IGNORE         true
$end

View output