Completely Multipolar Model as a General Framework for Many-Body Interactions as Illustrated for Water.

We introduce a general framework for many-body force fields, the Completely Multipolar Model (CMM), that utilizes multipolar electrical moments modulated by exponential decay of electron density as a common functional form for all terms of an energy decomposition analysis of intermolecular interactions. With this common functional form, the CMM model establishes well-formulated damped tensors that reach the correct asymptotes at both long- and short-range while formally ensuring no short-range catastrophes. CMM describes the separable EDA terms of dispersion, exchange polarization, and Pauli repulsion with short-ranged anisotropy, polarization as intramolecular charge fluctuations and induced dipoles, while charge transfer describes explicit movement of charge between molecules, and naturally describes many-body charge transfer by coupling into the polarization equations.

Simple and Accurate One-Body Energy and Dipole Moment Surfaces for Water and Beyond.

Water is often the testing ground for new, advanced force fields. While advanced functional forms for intermolecular interactions have been integral to the development of accurate water models, less attention has been paid to a transferable model for intramolecular valence terms. In this work, we present a one-body energy and dipole moment surface model, named 1B-UCB, that is simple yet accurate and can be feasibly adapted for both standard and advanced potentials.

The role of charge in microdroplet redox chemistry.

In charged water microdroplets, which occur in nature or in the lab upon ultrasonication or in electrospray processes, the thermodynamics for reactive chemistry can be dramatically altered relative to the bulk phase. Here, we provide a theoretical basis for the observation of accelerated chemistry by simulating water droplets of increasing charge imbalance to create redox agents such as hydroxyl and hydrogen radicals and solvated electrons. We compute the hydration enthalpy of OH and H that controls the electron transfer process, and the corresponding changes in vertical ionization energy and vertical electron affinity of the ions, to create OH and H reactive species.

Descriptors of water aggregation.

We rely on a total of 23 (cluster size, 8 structural, and 14 connectivity) descriptors to investigate structural patterns and connectivity motifs associated with water cluster aggregation. In addition to the cluster size n (number of molecules), the 8 structural descriptors can be further categorized into (i) one-body (intramolecular): covalent OH bond length (rOH) and HOH bond angle (θHOH), (ii) two-body: OO distance (rOO), OHO angle (θOHO), and HOOX dihedral angle (ϕHOOX), where X lies on the bisector of the HOH angle, (iii) three-body: OOO angle (θOOO), and (iv) many-body: modified tetrahedral order parameter (q) to account for two-, three-, four-, five-coordinated molecules (qm, m = 2, 3, 4, 5) and radius of gyration (Rg). The 14 connectivity descriptors are all many-body in nature and consist of the AD, AAD, ADD, AADD, AAAD, AAADD adjacencies [number of hydrogen bonds accepted (A) and donated (D) by each water molecule], Wiener index, Average Shortest Path Length, hydrogen bond saturation (% HB), and number of non-short-circuited three-membered cycles, four-membered cycles, five-membered cycles, six-membered cycles, and seven-membered cycles.

Using machine learning to go beyond potential energy surface benchmarking for chemical reactivity.

We train an equivariant machine learning (ML) model to predict energies and forces for hydrogen combustion under conditions of finite temperature and pressure. This challenging case for reactive chemistry illustrates that ML potential energy surfaces are difficult to make complete, due to overreliance on chemical intuition of what data are important for training. Instead, a 'negative design' data acquisition strategy using metadynamics as part of an active learning workflow helps to create a ML model that avoids unforeseen high-energy or unphysical energy configurations.

Predicting the Raman Spectra of Liquid Water with a Monomer-Field Model.

The Raman spectrum of liquid water is quite complex, reflecting its strong sensitivity to the local environment of the individual waters. The OH-stretch region of the spectrum, which captures the influence of hydrogen bonding, has only just begun to be unraveled. Here we develop a model for predicting the Raman spectra of the OH-stretch region by considering how local electric fields distort the energy surface of each water monomer.

Intra-cluster Charge Migration upon Hydration of Protonated Formic Acid Revealed by Anharmonic Analysis of Cold Ion Vibrational Spectra.

The rates of many chemical reactions are accelerated when carried out in micron-sized droplets, but the molecular origin of the rate acceleration remains unclear. One example is the condensation reaction of 1,2-diaminobenzene with formic acid to yield benzimidazole. The observed rate enhancements have been rationalized by invoking enhanced acidity at the surface of methanol solvent droplets with low water content to enable protonation of formic acid to generate a cationic species (protonated formic acid or PFA) formed by attachment of a proton to the neutral acid.

Using Diffusion Maps to Analyze Reaction Dynamics for a Hydrogen Combustion Benchmark Dataset.

We use local diffusion maps to assess the quality of two types of collective variables (CVs) for a recently published hydrogen combustion benchmark dataset that contains ab initio molecular dynamics (MD) trajectories and normal modes along minimum energy paths. This approach was recently advocated in for assessing CVs and analyzing reactions modeled by classical MD simulations. We report the effectiveness of this approach to molecular systems modeled by quantum ab initio MD.

M-Chem: a Modular Software Package for Molecular Simulation that Spans Scientific Domains.

We present a new software package called M-Chem that is designed from scratch in C++ and parallelized on shared-memory multi-core architectures to facilitate efficient molecular simulations. Currently, M-Chem is a fast molecular dynamics (MD) engine that supports the evaluation of energies and forces from two-body to many-body all-atom potentials, reactive force fields, coarse-grained models, combined quantum mechanics molecular mechanics (QM/MM) models, and external force drivers from machine learning, augmented by algorithms that are focused on gains in computational simulation times. M-Chem also includes a range of standard simulation capabilities including thermostats, barostats, multi-timestepping, and periodic cells, as well as newer methods such as fast extended Lagrangians and high quality electrostatic potential generation.

Isotope Effects in the Zundel-Eigen Isomerization of H(HO).

The isomerization pathway between the energetically low-lying Zundel and Eigen isomers of the protonated water hexamer was investigated using high-level calculations including a treatment of zero-point corrections. On the basis of these calculations, the Zundel-Eigen isomerization was found to proceed through a stable intermediate isomer, which consists of a four-membered ring with two single acceptor water molecules. The inclusion of vibrational zero-point energy is shown to be important for accurately establishing the relative energies of the three relevant isomers involved in the Zundel-Eigen isomerization.

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.

Spontaneous Formation of Hydrogen Peroxide in Water Microdroplets.

There is accumulating evidence that many chemical reactions are accelerated by several orders of magnitude in micrometer-sized aqueous or organic liquid droplets compared to their corresponding bulk liquid phase. However, the molecular origin of the enhanced rates remains unclear as in the case of spontaneous appearance of 1 μM hydrogen peroxide in water microdroplets. In this Letter, we consider the range of ionization energies and whether interfacial electric fields of a microdroplet can feasibly overcome the high energy step from hydroxide ions (OH) to hydroxyl radicals (OH) in a primary HO mechanism.

Hydrogen bond arrangements in (HO) clathrate hydrate cages: Optimization and many-body analysis.

We provide a detailed study of hydrogen bonding arrangements, relative stability, residual entropy, and an analysis of the many-body effects in the (HO) (D-cage), (HO) (T-cage), and (HO) (H-cage) hollow cages making up structures I (sI) and II (sII) of clathrate hydrate lattices. Based on the enumeration of the possible hydrogen bonding networks for a fixed oxygen atom scaffold, the residual entropy (S) of these three gas phase cages was estimated at 0.754 82, 0.

Molecular Dynamics Driven by the Many-Body Expansion (MBE-MD).

We present a protocol for classical and nuclear quantum dynamics, in which the energies and forces are generated by the many-body expansion (MBE), and apply it to water clusters using the TTM2.1-F and MB-Pol interaction potentials at various temperatures. We carry out MBE-molecular dynamics (MD) classical and nuclear quantum dynamical simulations, in which the energies and forces of the full system are approximated by the two-, three-, and four-body terms of the MBE, and compare the average potential and the vibrational density of states with the full simulation, i.

Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH, CO, and HS in (HO) Cages.

We present benchmark binding energies of naturally occurring gas molecules CH, CO, and HS in the small cage, namely, the pentagonal dodecahedron (5) (HO), which is one of the constituent cages of the 3 major lattices (structures I, II, and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with an increasing basis set to the complete basis set (CBS) limit, necessitating the use of large basis sets up to the aug-cc-pV5Z and subsequent correction for basis set superposition error (BSSE). For the host hollow (HO) cages, we have identified a most stable isomer with binding energy of -200.

The many-body expansion for aqueous systems revisited: III. Hofmeister ion-water interactions.

We report a Many Body Energy (MBE) analysis of aqueous ionic clusters containing anions and cations at the two opposite ends of the Hofmeister series, viz. the kosmotropes Ca and SO and the chaotropes NH and ClO, with 9 water molecules to quantify how these ions alter the interaction between the water molecules in their immediate surroundings. We specifically aim at quantifying how various ions (depending on their position in the Hofmeister series) affect the interaction between the surrounding water molecules and probe whether there is a qualitatively different behavior between kosmotropic vs.

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.

The Many-Body Expansion for Aqueous Systems Revisited: I. Water-Water Interactions.

We revisit the many-body expansion (MBE) for water-water interactions by examining the effects of the basis set, including those resulting from the basis set superposition error (BSSE) correction, and electron correlation on the various terms for selected sizes of water clusters up to = 21. The analysis is performed at the second-order Møller-Plesset (MP2) perturbation theory with the family of augmented correlation consistent basis sets up to five zeta quality (aug-cc-pVZ, = D, T, Q, 5) for the (HO), = 7, 10, 13, 16, and 21, clusters for which we report either the complete MBE (for = 7, 10) or the ones through the 6-body (for = 13) and the 5-body terms (for = 16, 21). For the = 3 and 7 clusters, we also report the analysis at the coupled cluster with single, double, and perturbative triple replacements in order to assess the effects of a higher correlation on the magnitude and percentage of the various MBE terms.

Mapping the temperature-dependent and network site-specific onset of spectral diffusion at the surface of a water cluster cage.

We explore the kinetic processes that sustain equilibrium in a microscopic, finite system. This is accomplished by monitoring the spontaneous, time-dependent frequency evolution (the frequency autocorrelation) of a single OH oscillator, embedded in a water cluster held in a temperature-controlled ion trap. The measurements are carried out by applying two-color, infrared-infrared photodissociation mass spectrometry to the DO·(HDO)(DO) isotopologue of the "magic number" protonated water cluster, H·(HO) The OH group can occupy any one of the five spectroscopically distinct sites in the distorted pentagonal dodecahedron cage structure.

A look inside the black box: Using graph-theoretical descriptors to interpret a Continuous-Filter Convolutional Neural Network (CF-CNN) trained on the global and local minimum energy structures of neutral water clusters.

We describe a method for the post-hoc interpretation of a neural network (NN) trained on the global and local minima of neutral water clusters. We use the structures recently reported in a newly published database containing over 5 × 10 unique water cluster networks (HO) of size N = 3-30. The structural properties were first characterized using chemical descriptors derived from graph theory, identifying important trends in topology, connectivity, and polygon structure of the networks associated with the various minima.

Atlas of putative minima and low-lying energy networks of water clusters n = 3-25.

We report a database consisting of the putative minima and ∼3.2 × 10 local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3-25 using an improved version of the Monte Carlo temperature basin paving (MCTBP) global optimization procedure in conjunction with the ab initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.

Beyond Badger's Rule: The Origins and Generality of the Structure-Spectra Relationship of Aqueous Hydrogen Bonds.

The structure of hydrogen bonded networks is intimately intertwined with their dynamics. Despite the incredibly wide range of hydrogen bond strengths encountered in water clusters, ion-water clusters, and liquid water, we demonstrate that the previously reported correlation between the change in the equilibrium bond length of the hydrogen bonded OH covalent bond and the corresponding shift in its harmonic frequency in water clusters is much more broadly applicable. Surprisingly, this correlation describes the ratios for both the equilibrium OH bond length/harmonic frequency and the vibrationally averaged bond length/anharmonic frequency in water, hydronium water, and halide water clusters.

Benchmark Electronic Structure Calculations for HO(HO) , n = 0-5, Clusters and Tests of an Existing 1,2,3-Body Potential Energy Surface with a New 4-Body Correction.

We report extensive benchmark CCSD(T) Complete Basis Set (CBS) estimates of the binding energies, structures, and harmonic frequencies of HO(HO) clusters, n = 0-5, including all currently known low-lying energy isomers. These are used to test a previously reported many-body (up to 3-body interactions) CCSD(T)-based potential energy surface (PES) for the hydrated proton. A new 4-body term for the hydronium-water-water-water interactions is introduced.

Origin of Many-Body Vibrational Frequency Shifts in Water Clusters.

We have demonstrated the application of many-body expansions to calculations of the anharmonic, local-mode, OH-stretching vibrational frequencies of water clusters. We focused on five low-lying isomers of the water hexamer and the DD*(20,1) isomer of (HO). Our approach provides accurate OH-stretching vibrational frequencies when treating one- and two-body interactions with the CCSD(T)-F12 level of theory and the three- and four-body interactions with the DF-MP2-F12 level.