Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development.

Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born-Oppenheimer potential energy surface. In particular, we focus on ab initio-based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through (a) parameterization of distributed multipoles, (b) explicit fitting of the MBE, (c) inclusion of many-atom features in a neural network, and (d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.

Cn (n=2-4): current status.

The major aspects of the C2, C3 and C4 elemental carbon clusters are surveyed. For C2, a brief analysis of its current status is presented. Regarding C3, the most recent results obtained in our group are reviewed with emphasis on modelling its potential energy surface which is particularly complicated due to the presence of multiple conical intersections. As for C4, the most stable isomeric forms of both triplet and singlet spin states and their possible interconversion pathways are examined afresh by means of accurate ab initio calculations. The main strategies for modelling the ground triplet C4 potential are also discussed. Starting from a truncated cluster expansion and a previously reported DMBE form for C3, an approximate four-body term is calibrated from the ab initio energies. The final six-dimensional global DMBE form so obtained reproduces all known topographical aspects while providing an accurate description of the C4 linear-rhombic isomerization pathway. It is therefore commended for both spectroscopic and reaction dynamics studies.This article is part of the theme issue 'Modern theoretical chemistry'.

Towards J/mol Accuracy for the Cohesive Energy of Solid Argon.

The cohesive energies of argon in its cubic and hexagonal closed packed structures are computed with an unprecedented accuracy of about 5 J mol(-1) (corresponding to 0.05 % of the total cohesive energy). The same relative accuracy with respect to experimental data is also found for the face-centered cubic lattice constant deviating by ca. 0.003 Å. This level of accuracy was enabled by using high-level theoretical, wave-function-based methods within a many-body decomposition of the interaction energy. Static contributions of two-, three-, and four-body fragments of the crystal are all individually converged to sub-J mol(-1) accuracy and complemented by harmonic and anharmonic vibrational corrections. Computational chemistry is thus achieving or even surpassing experimental accuracy for the solid-state rare gases.

Automatic molecular fragmentation by evolutionary optimisation.

Molecular fragmentation is an effective suite of approaches to reduce the formal computational complexity of quantum chemistry calculations while enhancing their algorithmic parallelisability. However, the practical applicability of fragmentation techniques remains hindered by a dearth of automation and effective metrics to assess the quality of a fragmentation scheme. In this article, we present the Quick Fragmentation via Automated Genetic Search (QFRAGS), a novel automated fragmentation algorithm that uses a genetic optimisation procedure to generate molecular fragments that yield low energy errors when adopted in Many Body Expansions (MBEs). Benchmark testing of QFRAGS on protein systems with less than 500 atoms, using two-body (MBE2) and three-body (MBE3) MBE calculations at the HF/6-31G* level, reveals mean absolute energy errors (MAEE) of 20.6 and 2.2 kJ  mol - 1 , respectively. For larger protein systems exceeding 500 atoms, MAEEs are 181.5 kJ  mol - 1 for MBE2 and 24.3 kJ  mol - 1 for MBE3. Furthermore, when compared to three manual fragmentation schemes on a 40-protein dataset, using both MBE and Fragment Molecular Orbital techniques, QFRAGS achieves comparable or often lower MAEEs. When applied to a 10-lipoglycan/glycolipid dataset, MAEs of 7.9 and 0.3 kJ  mol - 1 were observed at the MBE2 and MBE3 levels, respectively.Scientific Contribution This Article presents the Quick Fragmentation via Automated Genetic Search (QFRAGS), an innovative molecular fragmentation algorithm that significantly improves upon existing molecular fragmentation approaches by specifically addressing their lack of automation and effective fragmentation quality metrics. With an evolutionary optimisation strategy, QFRAGS actively pursues high quality fragments, generating fragmentation schemes that exhibit minimal energy errors on systems with hundreds to thousands of atoms. The advent of QFRAGS represents a significant advancement in molecular fragmentation, greatly improving the accessibility and computational feasibility of accurate quantum chemistry calculations.

Adaptive Partitioning QM/MM Dynamics Simulations for Substrate Uptake, Product Release, and Solvent Exchange.

Combined quantum mechanics/molecular mechanics (QM/MM) plays an important role in multiscale simulations of biological systems including enzymes. The adaptive-partitioning (AP) schemes surpass the conventional QM/MM methods in that they allow the on-the-fly, smooth exchange of particles between QM and MM subsystems in molecular dynamics simulations, leading to a seamless and dynamic integration of the QM and MM realms. Originally developed for simulating ion solvation in bulk solutions, the AP schemes have recently been extended to the treatment of proteins, fostering applications in the simulations of enzymes. The present contribution provides a detailed account of the AP schemes. We delineate the background of the algorithms and their parallel implementation, as well as offer practical advice and examples for their applications in the simulations of biological systems.

Water 26-mers Drawn from Bulk Simulations: Benchmark Binding Energies for Unprecedentedly Large Water Clusters and Assessment of the Electrostatically Embedded Three-Body and Pairwise Additive Approximations.

It is important to test methods for simulating water, but small water clusters for which benchmarks are available are not very representative of the bulk. Here we present benchmark calculations, in particular CCSD(T) calculations at the complete basis set limit, for water 26-mers drawn from Monte Carlo simulations of bulk water. These clusters are large enough that each water molecule participates in 2.5 hydrogen bonds on average. The electrostatically embedded three-body approximation with CCSD(T) embedded dimers and trimers reproduces the relative binding energies of eight clusters with a mean unsigned error (MUE, kcal per mole of water molecules) of only 0.009 and 0.015 kcal for relative and absolute binding energies, respectively. Using only embedded dimers (electrostatically embedded pairwise approximation) raises these MUEs to 0.038 and 0.070 kcal, and computing the energies with the M11 exchange-correlation functional, which is very economical, yields errors of only 0.029 and 0.042 kcal.

Piezoelectric Effects of Applied Electric Fields on Hydrogen-Bond Interactions: First-Principles Electronic Structure Investigation of Weak Electrostatic Interactions.

The piezoelectric properties of 2-methyl-4-nitroaniline crystals were explored qualitatively and quantitatively using an electrostatically embedded many-body (EE-MB) expansion scheme for the correlation energies of a system of monomers within the crystal. The results demonstrate that hydrogen bonding is an inherently piezoelectric interaction, deforming in response to the electrostatic environment. We obtain piezo-coefficients in excellent agreement with the experimental values. This approach reduces computational cost and reproduces the total resolution of the identity (RI)-Møller-Plesset second-order perturbation theory (RI-MP2) energy for the system to within 1.3 × 10(-5)%. Furthermore, the results suggest novel ways to self-assemble piezoelectric solids and suggest that accurate treatment of hydrogen bonds requires precise electrostatic evaluation. Considering the ubiquity of hydrogen bonds across chemistry, materials, and biology, a new electromechanical view of these interactions is required.

Achieving the CCSD(T) Basis-Set Limit in Sizable Molecular Clusters: Counterpoise Corrections for the Many-Body Expansion.

An efficient procedure is introduced to obtain the basis-set limit in electronic structure calculations of large molecular and ionic clusters. This approach is based on a Boys-Bernardi-style counterpoise correction for clusters containing arbitrarily many monomer units, which is rendered computationally feasible by means of a truncated many-body expansion. This affords a tractable way to apply the sequence of correlation-consistent basis sets (aug-cc-pVXZ) to large systems and thereby obtain energies extrapolated to the complete basis set (CBS) limit. A three-body expansion with three-body counterpoise corrections is shown to afford errors of ≲0.1-0.2 kcal/mol with respect to traditional MP2/CBS results, even for challenging systems such as fluoride-water clusters. A triples correction, δCCSD(T) = ECCSD(T) - EMP2, can be estimated accurately and efficiently as well. Because the procedure is embarrassingly parallelizable and requires no electronic structure calculations in systems larger than trimers, it is extendible to very large clusters. As compared to traditional CBS extrapolations, computational time is dramatically reduced even without parallelization.