An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2-25

Systematic improvement of empirical energy functions in the era of machine learning

The impact of targeted replacement of individual terms in empirical force fields is quantitatively assessed for pure water, dichloromethane (CH  2 Cl  2 ), and solvated K  + and Cl  - ions. For the electrostatic interactions, point charges (PCs) and machine learning (ML)-based minimally distributed charges (MDCM) fitted to the molecular electrostatic potential are evaluated together with electrostatics based on the Coulomb integral. The impact of explicitly including second-order terms is investigated by adding a fragment molecular orbital (FMO)-derived polarization energy to an existing force field, in this case CHARMM. It is demonstrated that anisotropic electrostatics reduce the RMSE for water (by 1.4 kcal/mol), CH  2 Cl  2 (by 0.8 kcal/mol) and for solvated Cl  - clusters (by 0.4 kcal/mol). An additional polarization term can be neglected for CH  2 Cl  2 but further improves the models for pure water (by ∼ 1.0 kcal/mol) and hydrated Cl  - (by 0.4 kcal/mol), and is key for solvated K  + , reducing the RMSE by 2.3 kcal/mol. A 12-6 Lennard-Jones functional form performs satisfactorily with PC and MDCM electrostatics, but is not appropriate for descriptions that account for the electrostatic penetration energy. The importance of many-body contributions is assessed by comparing a strictly 2-body approach with self-consistent reference data. Two-body interactions suffice for CH  2 Cl  2 whereas water and solvated K  + and Cl  - ions require explicit many-body corrections. Finally, a many-body-corrected dimer potential energy surface exceeds the accuracy attained using a conventional empirical force field, potentially reaching that of an FMO calculation. The present work systematically quantifies which terms improve the performance of an existing force field and what reference data to use for parametrizing these terms in a tractable fashion for ML fitting of pure and heterogeneous systems.

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. 1B-UCB for water is comparable in accuracy to those with much more complex functional forms, despite having drastically fewer parameters. The parametrization protocol has been implemented as part of the Q-Force automated workflow and requires only a quantum mechanical Hessian calculation as reference data, hence allowing it to be easily extended to a variety of molecular systems beyond water, which we demonstrate on a selection of small molecules with different symmetries.

The Near-Sightedness of Many-Body Interactions in Anharmonic Vibrational Couplings

Couplings between vibrational motions are driven by electronic interactions, and these couplings carry special significance in vibrational energy transfer, multidimensional spectroscopy experiments, and simulations of vibrational spectra. In this investigation, the many-body contributions to these couplings are analyzed computationally in the context of clathrate-like alkali metal cation hydrates, including Cs+(H2O)20, Rb+(H2O)20, and K+(H2O)20, using both analytic and quantum-chemistry potential energy surfaces. Although the harmonic spectra and one-dimensional anharmonic spectra depend strongly on these many-body interactions, the mode-pair couplings were, perhaps surprisingly, found to be dominated by one-body effects, even in cases of couplings to low-frequency modes that involved the motion of multiple water molecules. The origin of this effect was traced mainly to geometric distortion within water monomers and cancellation of many-body effects in differential couplings, and the effect was also shown to be agnostic to the identity of the ion. These outcomes provide new understanding of vibrational couplings and suggest the possibility of improved computational methods for the simulation of infrared and Raman spectra.

On the brink of self-hydration: the water heptadecamer

For pure, neutral, isolated molecular clusters, (H2O)17 marks the transition from structures with all water molecules on the cluster surface to water self-hydration, i.e., cluster structures around one central water molecule. Getting this right with water model potentials turns out to be challenging. Even the best water potentials currently available, which reproduce collective properties very well, still deliver contradicting results for (H2O)17, when different low-energy isomers from global structure optimizations are examined. Interestingly, ab initio quantum chemistry also struggles with the only seemingly simple question if (H2O)17 is all-surface or water-centered. Hence, although the long history of water potential development may be entering its final phase, it is not quite finished yet.

Attosecond-pump attosecond-probe x-ray spectroscopy of liquid water

Attosecond-pump/attosecond-probe experiments have long been sought as the most straightforward method for observing electron dynamics in real time. Although there has been much success with overlapped near-infrared femtosecond and extreme ultraviolet attosecond pulses combined with theory, true attosecond-pump/attosecond-probe experiments have been limited. We used a synchronized attosecond x-ray pulse pair from an x-ray free-electron laser to study the electronic response to valence ionization in liquid water through all x-ray attosecond transient absorption spectroscopy (AX-ATAS). Our analysis showed that the AX-ATAS response is confined to the subfemtosecond timescale, eliminating any hydrogen atom motion and demonstrating experimentally that the 1b1 splitting in the x-ray emission spectrum is related to dynamics and is not evidence of two structural motifs in ambient liquid water.

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. We mined a previously reported database of 4 948 959 water cluster minima for (H2O)n, n = 3-25 to analyze the evolution and correlation of these descriptors for the clusters within 5 kcal/mol of the putative minima. It was found that rOH and % HB correlated strongly with cluster size n, which was identified as the strongest predictor of energetic stability. Marked changes in the adjacencies and cycle count were observed, lending insight into changes in the hydrogen bond network upon aggregation. A Principal Component Analysis (PCA) was employed to identify descriptor dependencies and group clusters into specific structural patterns across different cluster sizes. The results of this study inform our understanding of how water clusters evolve in size and what appropriate descriptors of their structural and connectivity patterns are with respect to system size, stability, and similarity. The approach described in this study is general and can be easily extended to other hydrogen-bonded systems.

Ice structures assembled from cubic water clusters of D2d and S4 symmetry

The study of self-assembly processes is of key importance for fundamental science and modern technologies. Cubic water clusters of D2d and S4 symmetry show great potential as building blocks for self-assembly. The objective of this paper is to construct possible ice structures formed by hydrogen bonding of these very stable water clusters. A number of such structures are herein presented, including quasi-2D and 3D ices as well as spatial layered and tubular ices. The energetics and structure of many configurations differing in the arrangement of hydrogen atoms in hydrogen bonds have been studied. It was established that the proton disorder of all such ices is of island type. The residual entropy of these ices is equal to ln(3)/4 in dimensionless form. For layered structures formed by the stacking of multiple bilayers, the determining role of the van der Waals interactions is shown. Note that, for all considered ices, the lowest-energy configurations are formed only by clusters of D2d symmetry.

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.

Accurate Calculation of Many-Body Energies in Water Clusters Using a Classical Geometry-Dependent Induction Model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The moment expansion is carried out up to the hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA classical interaction potentials for water. The classical induction model with no adjustable parameters reproduces the ab initio 3-/4-body terms contained in these two databases with a root-mean-square error (RMSE) of 0.104/0.058 and a mean-absolute error (MAE) of 0.054/0.026 kcal/mol, respectively. These results are on par with the ones obtained by fitting the same data using over 14,000 (for the 3-body) and 200 (for the 4-body) parameters via Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy, included in the calculation of the 3-body energy, was found to be small but not quite negligible. The model represents a practical, efficient, and transferable approach for obtaining accurate nonadditive interactions for multicomponent systems without the need to perform tens of thousands of high level electronic structure calculations and fitting them with PIPs.

Anharmonic Assignment of the Water Octamer Spectrum in the OH Stretch Region

We interface the quasi-classical trajectory approach with an ab initio potential energy surface for water to assign the vibrational spectroscopical features of the OH stretch region of the water octamer cluster, which is considered to be a precursor of ice. An attempt by Li et al. to assign their recent reference experiment involved lower-level calculations based on an ad hoc scaled harmonic approach. Differently from the conclusions of this previous assignment, which invoked the contribution of 5 conformers and a solvated form of the water heptamer in the spectrum, we find out that the spectroscopic features can be related to the 4 conformers of the octamer lying lower in energy.