Abstract:Conspectus The past 15 years have witnessed an explosion of activity in the field of fragment-based quantum chemistry, whereby ab initio electronic structure calculations are performed on very large systems by decomposing them into a large number of relatively small subsystem calculations and then reassembling the subsystem data in order to approximate supersystem properties. Most of these methods are based, at some level, on the so-called many-body (or "n-body") expansion, which ultimately requires calculations on monomers, dimers, ..., n-mers of fragments. To the extent that a low-order n-body expansion can reproduce supersystem properties, such methods replace an intractable supersystem calculation with a large number of easily distributable subsystem calculations. This holds great promise for performing, for example, "gold standard" CCSD(T) calculations on large molecules, clusters, and condensed-phase systems. The literature is awash in a litany of fragment-based methods, each with their own working equations and terminology, which presents a formidable language barrier to the uninitiated reader. We have sought to unify these methods under a common formalism, by means of a generalized many-body expansion that provides a universal energy formula encompassing not only traditional n-body cluster expansions but also methods designed for macromolecules, in which the supersystem is decomposed into overlapping fragments. This formalism allows various fragment-based methods to be systematically classified, primarily according to how the fragments are constructed and how higher-order n-body interactions are approximated. This classification furthermore suggests systematic ways to improve the accuracy. Whereas n-body approaches have been thoroughly tested at low levels of theory in small noncovalent clusters, we have begun to explore the efficacy of these methods for large systems, with the goal of reproducing benchmark-quality calculations, ideally meaning complete-basis CCSD(T). For high accuracy, it is necessary to deal with basis-set superposition error, and this necessitates the use of many-body counterpoise corrections and electrostatic embedding methods that are stable in large basis sets. Tests on small noncovalent clusters suggest that total energies of complete-basis CCSD(T) quality can indeed be obtained, with dramatic reductions in aggregate computing time. On the other hand, naive applications of low-order n-body expansions may benefit from significant error cancellation, wherein basis-set superposition error partially offsets the effects of higher-order n-body terms, affording fortuitously good results in some cases. Basis sets that afford reasonable results in small clusters behave erratically in larger systems and when high-order n-body expansions are employed. For large systems, and (H2O)N≳30 is large enough, the combinatorial nature of the many-body expansion presents the possibility of serious loss-of-precision problems that are not widely appreciated. Tight thresholds are required in the subsystem calculations in order to stave off size-dependent errors, and high-order expansions may be inherently numerically ill-posed. Moreover, commonplace script- or driver-based implementations of the n-body expansion may be especially susceptible to loss-of-precision problems in large systems. These results suggest that the many-body expansion is not yet ready to be treated as a "black-box" quantum chemistry method.

Analytical Structure Matching and Very Precise Approach to the Coulombic Quantum Three-Body Problem

Kinetic energy operator approach to the quantum three-body problem with Coulomb interactions

Hyperspherical Calculations of Ultralow-Energy Collisions in Coulomb Three-Body Systems

Wavefunction matching for solving quantum many-body problems

Hyperspherical explicitly correlated Gaussian approach for few-body systems with finite angular momentum

Molecular structure calculations: a unified quantum mechanical description of electrons and nuclei using explicitly correlated Gaussian functions and the global vector representation

The Formalism and Matrix Elements of a Complete Potential-Harmonic Scheme for Directly Solving the Schrödinger Equation of the Helium Atom

Computational method for general multicenter electronic structure calculations

Addressing Strong Correlation by Approximate Coupled-Pair Methods with Active-Space and Full Treatments of Three-Body Clusters

Application of the complex scaling method in solving three-body Coulomb scattering problem

Solving Schrodinger Equation for Three-Electron Quantum Systems by the Use of The Hyperspherical Function Method

Calculating the many-potential vacuum polarization density of the Dirac equation in the finite-basis approximation

Combining symmetry breaking and restoration with configuration interaction: a highly accurate many-body scheme applied to the pairing Hamiltonian

Structure and Decays of Nuclear Three-Body Systems: the Gamow Coupled-Channel Method in Jacobi Coordinates

Double Exponential Transformation For Computing Three-Centre Nuclear Attraction Integrals

Precision of Finite-Difference Representation in 3D Coordinate-Space Hartree-Fock-Bogoliubov Calculations

High-precision Calculation of the Geometric Quantities of Two-Electron Atoms Based on the Hylleraas Configuration-Interaction Basis

Computationally Efficient Molecular Integrals of Solid Harmonic Gaussian Orbitals Using Quantum Entanglement of Angular Momentum

Finite-Basis-Set Approach to the Two-Center Heteronuclear \\ Dirac Problem

Natural Orbitals and Targeted Non-Orthogonal Orbital Sets for Atomic Hyperfine Structure Multiconfiguration Calculations

Aiming for benchmark accuracy with the many-body expansion