X

Search Results

Searching....

12.13 The Generalized Many-Body Expansion Method

12.13.1 Introduction

(May 21, 2025)

The generalized many-body expansion (GMBE) method approximates the energy of a supersystem using the energies of its fragments or subsystems, determined by a distance-based threshold. 1120 Richard R. M., Herbert J. M.
J. Chem. Phys.
(2012), 137, pp. 064113.
Link
, 893 Mayhall N. J., Raghavachari K.
J. Chem. Theory Comput.
(2012), 8, pp. 2669.
Link
, 1121 Richard R. M., Herbert J. M.
J. Chem. Theory Comput.
(2013), 9, pp. 1408.
Link
This threshold can be defined by the minimum atomic distance between “unit fragments” of the supersystem. Alternatively, one may use a criterion based on the heavy-atom distance, which excludes hydrogen atoms. The resulting subsystems, identified based on this threshold, are referred to as primitive fragments (monomers). These primitive fragments serve as the foundation for generating a set of fragment calculations used to approximate the supersystem’s energy. To prevent redundancy, the principle of inclusion and exclusion (PIE) is enforced.

In the n-body GMBE, denoted as GMBE(n), the energy of the supersystem is expressed as

E=I(Nn)I(n), (12.74)

where N is the number of primitive fragments.

The GMBE implementation in Q-Chem is currently limited to first-order, GMBE(1), and employs a novel binning algorithm for efficiently generating the set of primitive fragments. 69 Ballesteros F., Tan J. A., Lao K. U.
J. Chem. Phys.
(2023), 159, pp. 074107.
Link
The supersystem’s energy is approximated using GMBE(1) as:

E InI(1). (12.75)

The intersection-corrected energy for fragment I is

I(1)=EI(1)-J(J>I)EIJ(1)+J,K(K>J>I)EIJK(1)-J,K,L(L>K>J>I)EIJKL(1)+, (12.76)

where EI(1), EIJ(1), EIJK(1), and EIJKL(1) denote energies respectively of the fragment I, the intersection of fragments I and J, the intersection of fragments I, J, and K, and the intersection of fragments I, J, K, and L.

Similarly, the GMBE(1) approximation for the supersystem’s density matrix is given by:

P In𝒫I(1) (12.77)

where the intersection-corrected density matrix sub-blocks for subset I are defined as:

𝒫I(1)=PI(1)-J(J>I)PIJ(1)+J,K(K>J>I)PIJK(1)-J,K,L(L>K>J>I)PIJKL(1)+ (12.78)

The GMBE(1) density matrix can be used to predict the supersystem’s energy through a single Fock build. Alternatively, it can serve as an initial guess in a SCF calculation for the supersystem.