A full-dimensional potential energy function (MB-pol) for simulations of water from the dimer to bulk phases is developed entirely from "first principles" by building upon the many-body expansion of the interaction energy. Specifically, the MB-pol potential is constructed by combining a highly accurate dimer potential energy surface [J. Chem. Theory Comput. 2013, 9, 5395] with explicit three-body and many-body polarization terms. The three-body contribution, expressed as a combination of permutationally invariant polynomials and classical polarizability, is derived from a fit to more than 12000 three-body energies calculated at the CCSD(T)/aug-cc-pVTZ level of theory, imposing the correct asymptotic behavior as predicted from "first principles". Here, the accuracy of MB-pol is demonstrated through comparison of the calculated third virial coefficient with the corresponding experimental data as well as through analysis of the relative energy differences of small clusters.

The MB-pol full-dimensional water potential introduced in the first two papers of this series [J. Chem. Theory Comput. 2013, 9, 5395 and J. Chem. Theory Comput. 2014, 10, 1599] is employed here in classical and quantum simulations of liquid water under ambient conditions. Comparisons with the available experimental data for several structural, thermodynamic, and dynamical properties indicate that MB-pol provides a highly accurate description of the liquid phase. Combined with previous analyses of the dimer vibration-rotation tunneling spectrum, second and third virial coefficients, and cluster structures and energies, the present results demonstrate that MB-pol represents a major step toward the long-sought "universal model" capable of describing the properties of water from the gas to the condensed phases.

Vibrational spectroscopy is a powerful technique to probe the structure and dynamics of water. However, deriving an unambiguous molecular-level interpretation of the experimental spectral features remains a challenge due to the complexity of the underlying hydrogen-bonding network. In this contribution, we present an integrated theoretical and computational framework (named many-body molecular dynamics or MB-MD) that, by systematically removing uncertainties associated with existing approaches, enables a rigorous modeling of vibrational spectra of water from quantum dynamical simulations. Specifically, we extend approaches used to model the many-body expansion of interaction energies to develop many-body representations of the dipole moment and polarizability of water. The combination of these "first-principles" representations with centroid molecular dynamics simulations enables the simulation of infrared and Raman spectra of liquid water under ambient conditions that, without relying on any ad hoc parameters, are in good agreement with the corresponding experimental results. Importantly, since the many-body energy, dipole, and polarizability surfaces employed in the simulations are derived independently from accurate fits to correlated electronic structure data, MB-MD allows for a systematic analysis of the calculated spectra in terms of both electronic and dynamical contributions. The present analysis suggests that, while MB-MD correctly reproduces both the shifts and the shapes of the main spectroscopic features, an improved description of quantum dynamical effects possibly combined with a dissociable water potential may be necessary for a quantitative representation of the OH stretch band.

The molecular characterization of the air/water interface is a key step in understanding fundamental multiphase phenomena ranging from heterogeneous chemical processes in the atmosphere to the hydration of biomolecules. The apparent simplicity of the air/water interface, however, masks an underlying complexity associated with the dynamic nature of the water hydrogen-bond network that has so far hindered an unambiguous characterization of its microscopic properties. Here, we demonstrate that the application of quantum many-body molecular dynamics, which enables spectroscopically accurate simulations of water from the gas to the condensed phase, leads to a definitive molecular-level picture of the interface region. For the first time, excellent agreement is obtained between the simulated vibrational sum-frequency generation spectrum and the most recent state-of-the-art measurements, without requiring any empirical frequency shift or ad hoc scaling of the spectral intensity. A systematic dissection of the spectral features demonstrates that a rigorous representation of nuclear quantum effects as well as of many-body energy and electrostatic contributions is necessary for a quantitative reproduction of the experimental data. The unprecedented accuracy of the simulations presented here indicates that quantum many-body molecular dynamics can enable predictive studies of aqueous interfaces, which by complementing analogous experimental measurements will provide unique molecular insights into multiphase and heterogeneous processes of relevance in chemistry, biology, materials science, and environmental research.

A full-dimensional model of water, HBB2-pol, derived entirely from "first-principles", is introduced and employed in computer simulations ranging from the dimer to the liquid. HBB2-pol provides excellent agreement with the measured second and third virial coefficients and, by construction, reproduces the dimer vibration-rotation-tunneling spectrum. The model also predicts the relative energy differences between isomers of small water clusters within the accuracy of highly correlated electronic structure methods. Importantly, when combined with simulation methods that explicitly include zero-point energy and quantum thermal motion, HBB2-pol accurately describes both structural and dynamical properties of the liquid phase.