Many-Body Basis Set Superposition Effect

Analysis of two overlapping fragmentation approaches in density matrix construction: GMBE-DM vs. ADMA

In this study, we conduct a comparative analysis of two density matrix construction methods: the generalized many-body expansion for building density matrices (GMBE-DM) based on the set-theoretical principle of inclusion/exclusion and the adjustable density matrix assembler (ADMA) based on the Mulliken-Mezey ansatz. We apply these methods to various noncovalent clusters, including water clusters, ion-water clusters, and ion-pair clusters, using both small 6-31G(d) and large def2-TZVPPD basis sets. Our findings reveal that the GMBE-DM method, particularly when combined with the purification scheme and truncation at the one-body level [GMBE(1)-DM-P], exhibits superior performance across all test systems and basis sets. In contrast, all ADMA set of methods show reasonable results only with small and compact basis sets. For example, GMBE(1)-DM-P outperforms the best ADMA method by at least 4 and 16 times with small and large basis sets, respectively, in the case of (H2O)N=6-55. This highlights the significance of the basis set choice for ADMA, which is even more critical than the fragmentation scheme, such as the size of subsystems, while GMBE-DM consistently produces accurate results irrespective of the chosen basis set. Consequently, the efficient and robust GMBE(1)-DM-P approach is recommended as a fragmentation method for generating accurate absolute and relative energies across different binding patterns and basis sets for noncovalent clusters.

Scalable generalized screening for high-order terms in the many-body expansion: Algorithm, open-source implementation, and demonstration

The many-body expansion lies at the heart of numerous fragment-based methods that are intended to sidestep the nonlinear scaling of ab initio quantum chemistry, making electronic structure calculations feasible in large systems. In principle, inclusion of higher-order n-body terms ought to improve the accuracy in a controllable way, but unfavorable combinatorics often defeats this in practice and applications with n ≥ 4 are rare. Here, we outline an algorithm to overcome this combinatorial bottleneck, based on a bottom-up approach to energy-based screening. This is implemented within a new open-source software application ("Fragme∩t"), which is integrated with a lightweight semi-empirical method that is used to cull subsystems, attenuating the combinatorial growth of higher-order terms in the graph that is used to manage the calculations. This facilitates applications of unprecedented size, and we report four-body calculations in (H2O)64 clusters that afford relative energies within 0.1 kcal/mol/monomer of the supersystem result using less than 10% of the unique subsystems. We also report n-body calculations in (H2O)20 clusters up to n = 8, at which point the expansion terminates naturally due to screening. These are the largest n-body calculations reported to date using ab initio electronic structure theory, and they confirm that high-order n-body terms are mostly artifacts of basis-set superposition error.

An Embedded Fragment Method for Molecules in Strong Magnetic Fields

An extension of the embedded fragment method for calculations on molecular clusters is presented, which includes strong external magnetic fields. The approach is flexible, allowing for calculations at the Hartree-Fock, current-density-functional theory, Møller-Plesset perturbation theory, and coupled-cluster levels using London atomic orbitals. For systems consisting of discrete molecular subunits, calculations using London atomic orbitals can be performed in a computationally tractable manner for systems beyond the reach of conventional calculations, even those accelerated by resolution-of-the-identity or Cholesky decomposition methods. To assess the applicability of the approach, applications to water clusters are presented, showing how strong magnetic fields enhance binding within the clusters. However, our calculations suggest that, contrary to previous suggestions in the literature, this enhanced binding may not be directly attributable to strengthening of hydrogen bonding. Instead, these results suggest that this arises for larger field strengths as a response of the system to the presence of the external field, which induces a charge density build up between the monomer units. The approach is embarrassingly parallel and its computational tractability is demonstrated for clusters of up to 103 water molecules in triple-ζ basis sets, which would correspond to conventional calculations with more than 12 000 basis functions.

From the gas phase to the solid state: The chemical bonding in the superheavy element flerovium

As early as 1975, Pitzer suggested that copernicium, flerovium and oganesson are volatile substances behaving noble-gas like because of their closed-shell configurations and accompanying relativistic effects. It is, however, precarious to predict the chemical bonding and physical behavior of a solid by knowledge of the atomic or molecular properties only. Copernicium and oganesson have been analyzed very recently by our group. Both are predicted to be semi-conductors and volatile substances with rather low melting and boiling points, which may justify a comparison with the noble gas elements. Here we study closed-shell flerovium in detail to predict solid-state properties including the melting point from a decomposition of the total energy into many-body forces derived from relativistic coupled-cluster and from density functional theory. The convergence of such a decomposition for flerovium is critically analyzed, and the problem of using density functional theory is highlighted. We predict that flerovium is in many ways not behaving like a typical noble gas element despite its closed-shell 7$p_{1/2}^2$ configuration and resulting weak interactions. Unlike for the noble gases, the many-body expansion in terms of the interaction energy is not converging smoothly. This makes the accurate prediction of phase transitions very difficult. Nevertheless, a first prediction by Monte-Carlo simulation estimates the melting point at $284\pm 50$ K. Furthermore, calculations for the electronic band gap suggests that flerovium is a semi-conductor similar to copernicium

Behavior of counterpoise correction in many-body molecular clusters of organic compounds: Hartree-Fock interaction energy perspective

The counterpoise (CP) correction by Boys and Bernardi has been well accepted as a reliable strategy to account for basis set superposition error (BSSE) in intermolecular complexes. The behavior of the CP correction was thoroughly studied in individual molecules of molecular complexes. This work studies the performance of the CP correction in many‐body clusters including three‐body clusters of organic compounds in the 3B‐69 dataset. Additionally, we used crystal structures of polymorphs of benzene, aspirin, and oxalyl dihydrazide (ODH) to construct a many‐body cluster dataset, abbreviated as the MBC‐36 dataset, consisting of two, four and eight molecules, and 16 molecules in the case of benzene. A series of Dunning's basis sets—cc‐pXZ and aug‐cc‐pXZ (X = D and T)—were used to predict CP‐corrected Hartree–Fock (HF) interaction energies of the 3B‐69 and MBC‐36 datasets. The CP‐corrected interaction energies were found to be basis‐set independent, whereas the non‐CP corrected interaction energies were found not to a follow a smooth exponential fitting as previously found for electronic energies of individual molecules. This observation was attributed to the presence of non‐additive induction forces in some clusters. Two 2 × 2 × 2 supercells of benzene polymorphs were constructed to explore the local nature of BSSE effects. A cut‐off radius of 10 Å was demonstrated to be sufficient to fully recover these effects. Although the behavior of CP correction was found to be non‐conventional in many‐body clusters of organic compounds, the use of a small basis set such as cc‐pVDZ showed excellent performance in the prediction of HF interaction energies.

Accelerating the Convergence of Self-Consistent Field Calculations Using the Many-Body Expansion

The balance between cost-effective and sufficiently accurate methods represents the proverbial "promised land" for quantum chemistry calculations. The burden thus falls upon theoretical and computational chemists to provide such alternatives to mitigate the issues that arise from the employ of finite computing resources. In this paper, we attempt to demonstrate the importance of the quality of the initial guess for the self-consistent field (SCF) calculation when considering cost reduction techniques. We broach this challenge by using the many body expansion (MBE) to yield high quality density matrices (DMs) which, in turn, are applied as an SCF initial guess. The MBE-DM approaches combined with purification schemes and distance-based cutoff schemes can serve as initial guesses to reduce the SCF cycles necessary for convergence or derive energy directly through one Fock build. To this end, four unique types of clusters including water clusters, fluoride anion water clusters, sodium cation water clusters, and ammonium-bisulfate salt clusters have been used to test the performance of MBE-DM where its truncation at three-body expansion, MBE(3)-DM, shows vast improvement for those four clusters with reductions in the number of SCF cycles up to 40% as compared with the traditional superposition of atomic densities (SAD) guess. Other types of typical initial guesses, superposition of atomic potentials (SAP) and basis set projection (BSP), perform much worse than MBE-DM and SAD. In addition, the MBE-DM shows consistency across an array of fragment types irrespective of charges, size, level of theory, and basis set selection. Through MBE(3)-DM with the distance cutoff and the average purification scheme, the energy can be obtained directly with a mere 3.2 mH of the mean absolute deviation (MAD) for (H₂O)_N=6-55 which is at least 73 times better than the energy prediction using the typical initial guesses (SAD, SAP, and BSP). The corresponding MAD per monomer is only 0.14 mH which reaches the threshold of the "dynamical accuracy". The promising results of the methods outlined in this paper not only indicate two direct routes for computational cost reduction but also lay the possible foundation for composite techniques (i.e., ab initio sampling) that make best use of near converged values as their starting point.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Density-Based Many-Body Expansion as an Efficient and Accurate Quantum-Chemical Fragmentation Method: Application to Water Clusters

Fragmentation methods based on the many-body expansion offer an attractive approach for the quantum-chemical treatment of large molecular systems, such as molecular clusters and crystals. Conventionally, the many-body expansion is performed for the total energy, but such an energy-based many-body expansion often suffers from a slow convergence with respect to the expansion order. For systems that show strong polarization effects such as water clusters, this can render the energy-based many-body expansion infeasible. Here, we establish a density-based many-body expansion as a promising alternative approach. By performing the many-body expansion for the electron density instead of the total energy and inserting the resulting total electron density into the total energy functional of density functional theory, one can derive a density-based energy correction, which in principle accounts for all higher-order polarization effects. Here, we systematically assess the accuracy of such a density-based many-body expansion for test sets of water clusters. We show that already a density-based two-body expansion is able to reproduce interaction energies per fragment within chemical accuracy and is able to accurately predict the energetic ordering as well as the relative interaction energies of different isomers of water clusters.

The Many-Body Expansion for Aqueous Systems Revisited: II. Alkali Metal and Halide Ion-Water Interactions

We present a detailed study of the many-body expansion (MBE) for alkali metal and halide ion-water interactions and quantify the effect of these ions on the strength of the surrounding aqueous hydrogen bonding environment. Building on our previous work on neutral water clusters [J. P. Heindel and S. S. Xantheas, J. Chem. Theor. Comput. 16 (11), 6843-6855 (2020)], we carry out the MBE for the alkali metal and halide ion-water clusters, Z^+/-(H₂O)₉, where Z = Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺, F^-, Cl^-, Br^-, and I^- and compare them with the results for a pure water cluster with the same number of "bodies", viz., (H₂O)₁₀. The 2-B ion-water (I-W) interaction accounts for a larger percentage of the total cluster binding energy compared to a pure water cluster of the same size, with the total 3-B term being smaller and of opposite sign (repulsive), whereas higher order terms are essentially negligible. The same oscillating behavior around zero for the MBE terms higher than the 5-B with a basis set that was reported for water clusters is also observed for the ion-water clusters considered here, with the basis set superposition error (BSSE) corrections amending this as in the water cluster case. A remarkable, linear anticorrelation between the total 2-B (I-W), the total 2-B (W-W), and also the 3-B (W-W-W) interactions is found, quantifying the effect of the different ions in disrupting and altering (weakening) the neighboring hydrogen bonded water network: stronger (I-W) interactions result in weaker (W-W) interactions. Additional linear correlations across the two series of alkali metals and halide ions were found between the 3-B (I-W-W) and the 2-B (I-W) as well as between the 3-B (I-W-W) and the 3-B (W-W-W) interactions, suggesting the existence of previously unrealized underlying physics governing these 2-B intermolecular and 3-B collective interactions. Our results further suggest a universal behavior of the two different families of ions (alkali metals and halides) for both the correlations of the various components of the total binding energies and the estimate of the 2-B BSSE correction, which is reported to follow a common profile for ion-water and water-water interactions when cast in terms of reduced distances and energies of the respective dimers. We expect the current results that quantify the interplay between ion-water and water-water interactions in aqueous clusters to impact the development of classical, ab initio-based force fields for monatomic ion solvation, whereas the insights into the nature of the BSSE to be critical in future ab initio-based, many-body molecular dynamics studies.

Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage

The many-body expansion (MBE) of energies of molecular clusters or solids offers a way to detect and analyze errors of theoretical methods that could go unnoticed if only the total energy of the system was considered. In this regard, the interaction between the methane molecule and its enclosing dodecahedral water cage, CH₄···(H₂O)₂₀, is a stringent test for approximate methods, including density functional theory (DFT) approximations. Hybrid and semilocal DFT approximations behave erratically for this system, with three- and four-body nonadditive terms having neither the correct sign nor magnitude. Here, we analyze to what extent these qualitative errors in different MBE contributions are conveyed to post-Kohn-Sham random-phase approximation (RPA), which uses approximate Kohn-Sham orbitals as its input. The results reveal a correlation between the quality of the DFT input states and the RPA results. Moreover, the renormalized singles energy (RSE) corrections play a crucial role in all orders of the many-body expansion. For dimers, RSE corrects the RPA underbinding for every tested Kohn-Sham model: generalized-gradient approximation (GGA), meta-GGA, (meta-)GGA hybrids, as well as the optimized effective potential at the correlated level. Remarkably, the inclusion of singles in RPA can also correct the wrong signs of three- and four-body nonadditive energies as well as mitigate the excessive higher-order contributions to the many-body expansion. The RPA errors are dominated by the contributions of compact clusters. As a workable method for large systems, we propose to replace those compact contributions with CCSD(T) energies and to sum up the remaining many-body contributions up to infinity with supermolecular or periodic RPA. As a demonstration of this approach, we show that for RPA(PBE0)+RSE it suffices to apply CCSD(T) to dimers and 30 compact, hydrogen-bonded trimers to get the methane-water cage interaction energy to within 1.6% of the reference value.

Energy-Screened Many-Body Expansion: A Practical Yet Accurate Fragmentation Method for Quantum Chemistry

We introduce an implementation of the truncated many-body expansion, MBE(n), in which the n-body corrections are screened using the effective fragment potential force field, and only those that exceed a specified energy threshold are computed at a quantum-mechanical level of theory. This energy-screened MBE(n) approach is tested at the n = 3 level for a sequence of water clusters, (H₂O)_N=6-34. A threshold of 0.25 kJ/mol eliminates more than 80% of the subsystem electronic structure calculations and is even more efficacious in that respect than is distance-based screening. Even so, the energy-screened MBE(3) method is faithful to a full-system quantum chemistry calculation to within 1-2 kJ/mol/monomer, even in good quality basis sets such as aug-cc-pVTZ. These errors can be reduced by means of a two-layer approach that involves a Hartree-Fock calculation for the entire cluster. Such a correction proves to be necessary in order to obtain accurate relative energies for conformational isomers of (H₂O)₂₀, but the cost of a full-system Hartree-Fock calculation remains smaller than the cost of three-body subsystem calculations at correlated levels of theory. At the level of second-order Møller-Plesset perturbation theory (MP2), a screened MBE(3) calculation plus a full-system Hartree-Fock calculation is less expensive than a full-system MP2 calculation starting at N = 12 water molecules. This is true even if all MBE(3) subsystem calculations are performed on a single 40-core compute node, i.e., without significant parallelization. Energy-screened MBE(n) thus provides a fragment-based method that is accurate, stable in large basis sets, and low in cost, even when the latter is measured in aggregate computer time.

Fantasy versus reality in fragment-based quantum chemistry

Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.

A Linear-Scaling Method for Noncovalent Interactions: An Efficient Combination of Absolutely Localized Molecular Orbitals and a Local Random Phase Approximation Approach

A novel method for the accurate and efficient calculation of interaction energies in weakly bound complexes composed of a large number of molecules is presented. The new ALMO+RPAd method circumvents the prohibitive scaling of coupled cluster singles and doubles while still providing similar accuracy across a diverse range of weakly bound chemical systems. Linear-scaling procedures for the Fock build are given utilizing absolutely localized molecular orbitals (ALMOs), resulting in the a priori exclusion of basis set superposition errors. A bespoke data structure and algorithm using density fitting are described, leading to linear scaling for the storage and computation of the two-electron integrals. Electron correlation is included through a new, linear-scaling pairwise local random phase approximation approach, including exchange interactions, and decomposed into purely dispersive excitations (RPAxd). Collectively, these allow meaningful decomposition of the interaction energy into physically distinct contributions: electrostatic, polarization, charge transfer, and dispersion. Comparison with symmetry-adapted perturbation theory shows good qualitative agreement. Tests on various dimers and the S66 benchmark set demonstrate results within 0.5 kcal mol^-1 of coupled cluster singles and doubles results. On a large cluster of water molecules, we achieve calculations involving over 3500 orbital and 12,000 auxiliary basis functions in under 10 min on a single CPU core.

Multi-Body Interactions in Molecular Docking Program Devised with Key Water Molecules in Protein Binding Sites

Water molecules play an important role in modeling protein-ligand interactions. However, traditional molecular docking methods often ignore the impact of the water molecules by removing them without any analysis or keeping them as a static part of the proteins or the ligands. Hence, the accuracy of the docking simulations will inevitably be damaged. Here, we introduce a multi-body docking program which incorporates the fixed or the variable number of the key water molecules in protein-ligand docking simulations. The program employed NSGA II, a multi-objective optimization algorithm, to identify the binding poses of the ligand and the key water molecules for a protein. To this end, a force-field-based hydration-specific scoring function was designed to favor estimate the binding affinity considering the key water molecules. The program was evaluated in aspects of the docking accuracy, cross-docking accuracy, and screening efficiency. When the numbers of the key water molecules were treated as fixed-length optimization variables, the docking accuracy of the multi-body docking program achieved a success rate of 80.58% for the best RMSD values for the recruit of the ligands smaller than 2.0 Å. The cross-docking accuracy was investigated on the presence and absence of the key water molecules by four protein targets. The screening efficiency was assessed against those protein targets. Results indicated that the proposed multi-body docking program was with good performance compared with the other programs. On the other side, when the numbers of the key water molecules were treated as variable-length optimization variables, the program obtained comparative performance under the same three evaluation criterions. These results indicated that the multi-body docking with the variable numbers of the water molecules was also efficient. Above all, the multi-body docking program developed in this study was capable of dealing with the problem of the water molecules that explicitly participating in protein-ligand binding.

Full QM Calculation of RNA Energy Using Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps Method

In this study, the electrostatically embedded generalized molecular fractionation with conjugate caps (concaps) method (EE-GMFCC) was employed for efficient linear-scaling quantum mechanical (QM) calculation of total energies of RNAs. In the EE-GMFCC approach, the total energy of RNA is calculated by taking a proper combination of the QM energy of each nucleotide-centric fragment with large caps or small caps (termed EE-GMFCC-LC and EE-GMFCC-SC, respectively) deducted by the energies of concaps. The two-body QM interaction energy between non-neighboring ribonucleotides which are spatially in close contact are also taken into account for the energy calculation. Numerical studies were carried out to calculate the total energies of a number of RNAs using the EE-GMFCC-LC and EE-GMFCC-SC methods at levels of the Hartree-Fock (HF) method, density functional theory (DFT), and second-order many-body perturbation theory (MP2), respectively. The results show that the efficiency of the EE-GMFCC-SC method is about 3 times faster than the EE-GMFCC-LC method with minimal accuracy sacrifice. The EE-GMFCC-SC method is also applied for relative energy calculations of 20 different conformers of two RNA systems using HF and DFT, respectively. Both single-point and relative energy calculations demonstrate that the EE-GMFCC method has deviations from the full system results of only a few kcal/mol.

Understanding the many-body expansion for large systems. II. Accuracy considerations

To complement our study of the role of finite precision in electronic structure calculations based on a truncated many-body expansion (MBE, or "n-body expansion"), we examine the accuracy of such methods in the present work. Accuracy may be defined either with respect to a supersystem calculation computed at the same level of theory as the n-body calculations, or alternatively with respect to high-quality benchmarks. Both metrics are considered here. In applications to a sequence of water clusters, (H2O)N=6-55 described at the B3LYP/cc-pVDZ level, we obtain mean absolute errors (MAEs) per H2O monomer of ∼1.0 kcal/mol for two-body expansions, where the benchmark is a B3LYP/cc-pVDZ calculation on the entire cluster. Three- and four-body expansions exhibit MAEs of 0.5 and 0.1 kcal/mol/monomer, respectively, without resort to charge embedding. A generalized many-body expansion truncated at two-body terms [GMBE(2)], using 3-4 H2O molecules per fragment, outperforms all of these methods and affords a MAE of ∼0.02 kcal/mol/monomer, also without charge embedding. GMBE(2) requires significantly fewer (although somewhat larger) subsystem calculations as compared to MBE(4), reducing problems associated with floating-point roundoff errors. When compared to high-quality benchmarks, we find that error cancellation often plays a critical role in the success of MBE(n) calculations, even at the four-body level, as basis-set superposition error can compensate for higher-order polarization interactions. A many-body counterpoise correction is introduced for the GMBE, and its two-body truncation [GMBCP(2)] is found to afford good results without error cancellation. Together with a method such as ωB97X-V/aug-cc-pVTZ that can describe both covalent and non-covalent interactions, the GMBE(2)+GMBCP(2) approach provides an accurate, stable, and tractable approach for large systems.

Molecular energies from an incremental fragmentation method

The systematic molecular fragmentation method by Collins and Deev [J. Chem. Phys. 125, 104104 (2006)] has been used to calculate total energies and relative conformational energies for a number of small and extended molecular systems. In contrast to the original approach by Collins, we have tested the accuracy of the fragmentation method by utilising an incremental scheme in which the energies at the lowest level of the fragmentation are calculated on an accurate quantum chemistry level while lower-cost methods are used to correct the low-level energies through a high-level fragmentation. In this work, the fragment energies at the lowest level of fragmentation were calculated using the random-phase approximation (RPA) and two recently developed extensions to the RPA while the incremental corrections at higher levels of the fragmentation were calculated using standard density functional theory (DFT) methods. The complete incremental fragmentation method has been shown to reproduce the supermolecule results with a very good accuracy, almost independent on the molecular type, size, or type of decomposition. The fragmentation method has also been used in conjunction with the DFT-SAPT (symmetry-adapted perturbation theory) method which enables a breakdown of the total nonbonding energy contributions into individual interaction energy terms. Finally, the potential problems of the method connected with the use of capping hydrogen atoms are analysed and two possible solutions are supplied.

Pair-Pair Approximation to the Generalized Many-Body Expansion: An Alternative to the Four-Body Expansion for ab Initio Prediction of Protein Energetics via Molecular Fragmentation

We introduce a "pair-pair" approximation to the generalized many-body expansion (pp-GMBE) as an approximation to a traditional four-body expansion, the latter of which is accurate but quickly becomes numerically unstable and ultimately intractable as the number of "bodies" (fragments) increases. The pp-GMBE method achieves a good balance between accuracy and efficiency by defining significant fragment pairs and then fragment quartets. An efficient fragmentation scheme is introduced for proteins such that the largest subsystems contain about 60 atoms. Application of the pp-GMBE method to proteins with as many as 70 amino acids (1142 atoms) reveals that pp-GMBE energies are quite faithful to those obtained when the same level of density functional theory is applied to the entire macromolecule. When combined with embedding charges obtained from natural population analysis, the pp-GMBE approach affords absolute energies that differ by 1-3 kcal/mol from full supersystem results, but it yields conformational energy profiles that are practically indistinguishable from the supersystem calculation at the same level of theory.

Are fragment-based quantum chemistry methods applicable to medium-sized water clusters?

The GEBF method is demonstrated to be more accurate than the EE-MB method for medium-sized water clusters.