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 C₂, C₃ and C₄ elemental carbon clusters are surveyed. For C₂, a brief analysis of its current status is presented. Regarding C₃, 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 C₄, 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 C₄ potential are also discussed. Starting from a truncated cluster expansion and a previously reported DMBE form for C₃, 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 C₄ 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.

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.

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.

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.