Many-Body Expansion

The total energy of an N-body system is expressed as a sum of 1-body (monomer), 2-body (dimer), 3-body (trimer), and higher-order contributions. Each term captures the incremental interaction energy beyond what lower-order terms already account for. The expansion is formally exact when carried to the N-body term, but in practice, truncation at 2-body or 3-body recovers >95% of the total interaction energy for most molecular systems due to the short-ranged nature of many-body effects in non-metallic systems.

The molecular system is partitioned into non-overlapping fragments (monomers) based on chemical intuition—typically along covalent bonds capped with hydrogen link atoms or via orbital-based localization methods. The quality of the expansion depends critically on the fragmentation scheme: small, electronically isolated fragments converge faster, while large conjugated systems require careful treatment. Common approaches include molecular fractionation with conjugate caps (MFCC), electrostatically embedded MBE, and generalized energy-based fragmentation (GEBF).

To accelerate convergence of the many-body expansion, electrostatic embedding places the charge distribution of all other fragments as a background potential when calculating each n-body term. This captures long-range polarization effects at the monomer level, dramatically reducing the need for higher-order terms. The embedding potential can be derived from atomic partial charges, multipole expansions, or self-consistently updated electron densities in polarizable embedding schemes.

While a full supersystem calculation scales as O(N³) to O(N⁷) depending on the quantum mechanical method, MBE reduces this to a series of calculations on small subsystems. For a system of M fragments truncated at 2-body, the cost scales as O(M²) but with a small constant prefactor. Critically, all n-body calculations are embarrassingly parallel—each dimer or trimer can be computed independently, enabling near-linear scaling on high-performance computing clusters.

Standard MBE partitions molecules into disjoint fragments, but overlapping fragmentation methods like the incremental scheme or cluster-in-molecule (CIM) approach include buffer regions around each fragment. These overlapping regions capture non-local correlation effects more efficiently, allowing truncation at lower orders. The overlap introduces double-counting corrections that must be carefully subtracted using inclusion-exclusion principles analogous to the MBE formalism itself.

Not all n-body terms contribute equally. Distance-based screening discards dimers or trimers where fragments are separated beyond a cutoff (typically 4-6 Å for neutral systems). Energy-based screening estimates contributions using a cheap method (e.g., semi-empirical) and only promotes significant terms to high-level theory. This ONIOM-style multi-layer treatment recovers >99% of the canonical energy while computing <10% of the possible n-body terms, making MBE practical for systems with thousands of atoms.

A class of computational techniques that divide a large molecular system into smaller, computationally tractable subsystems. The many-body expansion is a foundational fragmentation approach, alongside others like electrostatically embedded methods and systematic molecular fragmentation. These methods exploit the nearsightedness of electronic matter, where local chemical environments dominate interactions. Key distinctions include:

• Inclusion-exclusion principle: Ensures no overcounting when combining fragment energies
• Capping strategies: Hydrogen atoms or link atoms terminate severed covalent bonds
• Distance-based truncation: Higher-order n-body terms decay rapidly with separation

A systematic expansion where the total correlation energy is expressed as a sum of increments from individual orbital groups, pairs, triples, and higher tuples. Closely related to the many-body expansion, the method of increments applies the fragmentation logic at the orbital level rather than the atomic level. Developed for local correlation methods, it enables:

• Linear scaling with system size for insulators
• Recovery of 99%+ of the total correlation energy at the three-body level
• Natural synergy with localized orbital representations like Foster-Boys or Pipek-Mezey orbitals

A framework that places a fragment in the electrostatic field of its environment to improve convergence of the many-body expansion. Density functional embedding and projection-based embedding reduce the truncation error at low orders by capturing long-range polarization. Key approaches include:

• QM/QM embedding: Both fragment and environment treated quantum mechanically
• QM/MM embedding: Environment modeled with classical point charges
• Density matrix embedding theory (DMET): Matches fragment and environment density matrices at the boundary Embedding can reduce the required expansion order from 4-body to 2-body for polar systems.

The physical phenomenon where distant many-body interactions are damped by intervening charge distributions, accelerating the convergence of the many-body expansion. In polarizable environments, screening effects mean that 4-body and higher terms often contribute less than 0.1 kcal/mol to the total energy. Factors affecting screening:

• Dielectric constant of the medium: Higher dielectric = stronger screening
• Ionic strength: Mobile ions shield long-range Coulomb interactions
• Conjugation length: Delocalized electrons mediate longer-range n-body effects Understanding screening is critical for setting truncation thresholds in fragment-based calculations.

The application of gold-standard coupled cluster theory (CCSD(T)) to individual n-body terms, enabling near-exact total energies for systems with hundreds of atoms. This approach combines the accuracy of wavefunction methods with the scalability of fragmentation. Practical implementations include:

• Fragment molecular orbital (FMO) with CCSD(T) corrections
• Incremental CCSD(T) for molecular crystals
• Local natural orbital CCSD(T) as an alternative to explicit fragmentation These methods can achieve sub-kcal/mol accuracy for non-covalent interaction energies in supramolecular systems.

An artificial lowering of fragment interaction energies caused when fragments borrow basis functions from neighboring fragments in the supersystem calculation. The many-body expansion is particularly susceptible to BSSE because each n-body term uses a different basis set context. Mitigation strategies include:

• Counterpoise correction: Using the full supersystem basis for each fragment calculation
• Explicitly correlated methods (F12): Dramatically reduce basis set dependence
• Complete basis set extrapolation: Removes BSSE by converging to the basis set limit Uncorrected BSSE can artificially stabilize dimers by 0.5-2.0 kcal/mol with double-zeta basis sets.