Signature properties of water: Their molecular electronic origins

Universal Pairwise Interatomic van der Waals Potentials Based on Quantum Drude Oscillators

Repulsive short-range and attractive long-range van der Waals (vdW) forces play an appreciable role in the behavior of extended molecular systems. When using empirical force fields, the most popular computational methods applied to such systems, vdW forces are typically described by Lennard-Jones-like potentials, which unfortunately have a limited predictive power. Here, we present a universal parameterization of a quantum-mechanical vdW potential, which requires only two free-atom properties─the static dipole polarizability α1 and the dipole-dipole C6 dispersion coefficient. This is achieved by deriving the functional form of the potential from the quantum Drude oscillator (QDO) model, employing scaling laws for the equilibrium distance and the binding energy, and applying the microscopic law of corresponding states. The vdW-QDO potential is shown to be accurate for vdW binding energy curves, as demonstrated by comparing to the ab initio binding curves of 21 noble-gas dimers. The functional form of the vdW-QDO potential has the correct asymptotic behavior at both zero and infinite distances. In addition, it is shown that the damped vdW-QDO potential can accurately describe vdW interactions in dimers consisting of group II elements. Finally, we demonstrate the applicability of the atom-in-molecule vdW-QDO model for predicting accurate dispersion energies for molecular systems. The present work makes an important step toward constructing universal vdW potentials, which could benefit (bio)molecular computational studies.

Optimized Quantum Drude Oscillators for Atomic and Molecular Response Properties

The quantum Drude oscillator (QDO) is an efficient yet accurate coarse-grained approach that has been widely used to model electronic and optical response properties of atoms and molecules as well as polarization and dispersion interactions between them. Three effective parameters (frequency, mass, and charge) fully characterize the QDO Hamiltonian and are adjusted to reproduce response properties. However, the soaring success of coupled QDOs for many-atom systems remains fundamentally unexplained, and the optimal mapping between atoms/molecules and oscillators has not been established. Here we present an optimized parametrization (OQDO) where the parameters are fixed by using only dipolar properties. For the periodic table of elements as well as small molecules, our model accurately reproduces atomic (spatial) polarization potentials and multipolar dispersion coefficients, elucidating the high promise of the presented model in the development of next-generation quantum-mechanical force fields for (bio)molecular simulations.

Supercritical Water is not Hydrogen Bonded

Thinking about water is inextricably linked to hydrogen bonds, which are highly directional in character and determine the unique structure of water, in particular its tetrahedral H-bond network. Here, we assess if this common connotation also holds for supercritical water. We employ extensive ab initio molecular dynamics simulations to systematically monitor the evolution of the H-bond network mode of water from room temperature, where it is the hallmark of its fluctuating three-dimensional network structure, to supercritical conditions. Our simulations reveal that the oscillation period required for H-bond vibrations to occur exceeds the lifetime of H-bonds in supercritical water by far. Instead, the corresponding low-frequency intermolecular vibrations of water pairs as seen in supercritical water are found to be well represented by isotropic van-der-Waals interactions only. Based on these findings, we conclude that water in its supercritical phase is not a H-bonded fluid.

Assessing the properties of supercritical water in terms of structural dynamics and electronic polarization effects

Evolution of water's structural dynamics from ambient liquid to supercritical dense liquid-like and dilute gas-like conditions.

Quantum mechanics of proteins in explicit water: The role of plasmon-like solute-solvent interactions

Quantum-mechanical van der Waals dispersion interactions play an essential role in intraprotein and protein-water interactions-the two main factors affecting the structure and dynamics of proteins in water. Typically, these interactions are only treated phenomenologically, via pairwise potential terms in classical force fields. Here, we use an explicit quantum-mechanical approach of density-functional tight-binding combined with the many-body dispersion formalism and demonstrate the relevance of many-body van der Waals forces both to protein energetics and to protein-water interactions. In contrast to commonly used pairwise approaches, many-body effects substantially decrease the relative stability of native states in the absence of water. Upon solvation, the protein-water dispersion interaction counteracts this effect and stabilizes native conformations and transition states. These observations arise from the highly delocalized and collective character of the interactions, suggesting a remarkable persistence of electron correlation through aqueous environments and providing the basis for long-range interaction mechanisms in biomolecular systems.

A FFLUX Water Model: Flexible, Polarizable and with a Multipolar Description of Electrostatics

Key to progress in molecular simulation is the development of advanced models that go beyond the limitations of traditional force fields that employ a fixed, point charge‐based description of electrostatics. Taking water as an example system, the FFLUX framework is shown capable of producing models that are flexible, polarizable and have a multipolar description of the electrostatics. The kriging machine‐learning methods used in FFLUX are able to reproduce the intramolecular potential energy surface and multipole moments of a single water molecule with chemical accuracy using as few as 50 training configurations. Molecular dynamics simulations of water clusters (25–216 molecules) using the new FFLUX model reveal that incorporating charge‐quadrupole, dipole–dipole, and quadrupole–charge interactions into the description of the electrostatics results in significant changes to the intermolecular structuring of the water molecules. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Theory and practice of modeling van der Waals interactions in electronic-structure calculations

The accurate description of long-range electron correlation, most prominently including van der Waals (vdW) dispersion interactions, represents a particularly challenging task in the modeling of molecules and materials. vdW forces arise from the interaction of quantum-mechanical fluctuations in the electronic charge density. Within (semi-)local density functional approximations or Hartree-Fock theory such interactions are neglected altogether. Non-covalent vdW interactions, however, are ubiquitous in nature and play a key role for the understanding and accurate description of the stability, dynamics, structure, and response properties in a plethora of systems. During the last decade, many promising methods have been developed for modeling vdW interactions in electronic-structure calculations. These methods include vdW-inclusive Density Functional Theory and correlated post-Hartree-Fock approaches. Here, we focus on the methods within the framework of Density Functional Theory, including non-local van der Waals density functionals, interatomic dispersion models within many-body and pairwise formulation, and random phase approximation-based approaches. This review aims to guide the reader through the theoretical foundations of these methods in a tutorial-style manner and, in particular, highlight practical aspects such as the applicability and the advantages and shortcomings of current vdW-inclusive approaches. In addition, we give an overview of complementary experimental approaches, and discuss tools for the qualitative understanding of non-covalent interactions as well as energy decomposition techniques. Besides representing a reference for the current state-of-the-art, this work is thus also designed as a concise and detailed introduction to vdW-inclusive electronic structure calculations for a general and broad audience.

Observation of the thermal influenced quantum behaviour of water near a solid interface

Water is essential for life. However, the structure and properties of water are still not well understood. It has been introduced that anomalies are in vicinal water near solid interfaces. We performed capillary flow experiments on water with a silica colloid sample using a high-performance liquid chromatography (HPLC) system by accurately varying the temperature and analysed the peak shape rigorously. We obtained a novel anomalous temperature spectrum from the peak-shape analysis. Here we report the observed distinct specific anomalous temperature (SAT) behaviour in vicinal water at silica interface. The anomaly appeared in the viscous force that was derived from a relationship between the shape of the HPLC peak and the velocity profile for the capillary flow. The observations were highly reproducible, and we conclude that the SAT is related to the quantum mechanical behaviour of water, in agreement of the characteristic acceptance of thermal displacement according to the Franck-Condon principle. We performed the same experiments using heavy water and water mixed with a small amount of methanol, and the results support the quantum phenomenological origin.

Building better water models using the shape of the charge distribution of a water molecule

The unique properties of liquid water apparently arise from more than just the tetrahedral bond angle between the nuclei of a water molecule since simple three-site models of water are poor at mimicking these properties in computer simulations. Four- and five-site models add partial charges on dummy sites and are better at modeling these properties, which suggests that the shape of charge distribution is important. Since a multipole expansion of the electrostatic potential describes a charge distribution in an orthogonal basis set that is exact in the limit of infinite order, multipoles may be an even better way to model the charge distribution. In particular, molecular multipoles up to the octupole centered on the oxygen appear to describe the electrostatic potential from electronic structure calculations better than four- and five-site models, and molecular multipole models give better agreement with the temperature and pressure dependence of many liquid state properties of water while retaining the computational efficiency of three-site models. Here, the influence of the shape of the molecular charge distribution on liquid state properties is examined by correlating multipoles of non-polarizable water models with their liquid state properties in computer simulations. This will aid in the development of accurate water models for classical simulations as well as in determining the accuracy needed in quantum mechanical/molecular mechanical studies and ab initio molecular dynamics simulations of water. More fundamentally, this will lead to a greater understanding of how the charge distribution of a water molecule leads to the unique properties of liquid water. In particular, these studies indicate that p-orbital charge out of the molecular plane is important.

Structure and hydrogen bonding at the limits of liquid water stability

Liquid water exhibits unconventional behaviour across its wide range of stability - from its unusually high liquid-vapour critical point down to its melting point and below where it reaches a density maximum and exhibits negative thermal expansion allowing ice to float. Understanding the molecular underpinnings of these anomalies presents a challenge motivating the study of water for well over a century. Here we examine the molecular structure of liquid water across its range of stability, from mild supercooling to the negative pressure and high temperature regimes. We use a recently-developed, electronically-responsive model of water, constructed from gas-phase molecular properties and incorporating many-body, long-range interactions to all orders; as a result the model has been shown to have high transferability from ice to the supercritical regime. We report a link between the anomalous thermal expansion of water and the behaviour of its second coordination shell and an anomaly in hydrogen bonding, which persists throughout liquid water's range of stability - from the high temperature limit of liquid water to its supercooled regime.

Long-Range Repulsion Between Spatially Confined van der Waals Dimers

It is an undisputed textbook fact that nonretarded van der Waals (vdW) interactions between isotropic dimers are attractive, regardless of the polarizability of the interacting systems or spatial dimensionality. The universality of vdW attraction is attributed to the dipolar coupling between fluctuating electron charge densities. Here, we demonstrate that the long-range interaction between spatially confined vdW dimers becomes repulsive when accounting for the full Coulomb interaction between charge fluctuations. Our analytic results are obtained by using the Coulomb potential as a perturbation over dipole-correlated states for two quantum harmonic oscillators embedded in spaces with reduced dimensionality; however, the long-range repulsion is expected to be a general phenomenon for spatially confined quantum systems. We suggest optical experiments to test our predictions, analyze their relevance in the context of intermolecular interactions in nanoscale environments, and rationalize the recent observation of anomalously strong screening of the lateral vdW interactions between aromatic hydrocarbons adsorbed on metal surfaces.

First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications

Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions

Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.

How van der Waals interactions determine the unique properties of water

Whereas the interactions between water molecules are dominated by strongly directional hydrogen bonds (HBs), it was recently proposed that relatively weak, isotropic van der Waals (vdW) forces are essential for understanding the properties of liquid water and ice. This insight was derived from ab initio computer simulations, which provide an unbiased description of water at the atomic level and yield information on the underlying molecular forces. However, the high computational cost of such simulations prevents the systematic investigation of the influence of vdW forces on the thermodynamic anomalies of water. Here, we develop efficient ab initio-quality neural network potentials and use them to demonstrate that vdW interactions are crucial for the formation of water's density maximum and its negative volume of melting. Both phenomena can be explained by the flexibility of the HB network, which is the result of a delicate balance of weak vdW forces, causing, e.g., a pronounced expansion of the second solvation shell upon cooling that induces the density maximum.

A Stochastic, Resonance-Free Multiple Time-Step Algorithm for Polarizable Models That Permits Very Large Time Steps

Molecular dynamics remains one of the most widely used computational tools in the theoretical molecular sciences to sample an equilibrium ensemble distribution and/or to study the dynamical properties of a system. The efficiency of a molecular dynamics calculation is limited by the size of the time step that can be employed, which is dictated by the highest frequencies in the system. However, many properties of interest are connected to low-frequency, long time-scale phenomena, requiring many small time steps to capture. This ubiquitous problem can be ameliorated by employing multiple time-step algorithms, which assign different time steps to forces acting on different time scales. In such a scheme, fast forces are evaluated more frequently than slow forces, and as the former are often computationally much cheaper to evaluate, the savings can be significant. Standard multiple time-step approaches are limited, however, by resonance phenomena, wherein motion on the fastest time scales limits the step sizes that can be chosen for the slower time scales. In atomistic models of biomolecular systems, for example, the largest time step is typically limited to around 5 fs. Previously, we introduced an isokinetic extended phase-space algorithm (Minary et al. Phys. Rev. Lett. 2004, 93, 150201) and its stochastic analog (Leimkuhler et al. Mol. Phys. 2013, 111, 3579) that eliminate resonance phenomena through a set of kinetic energy constraints. In simulations of a fixed-charge flexible model of liquid water, for example, the time step that could be assigned to the slow forces approached 100 fs. In this paper, we develop a stochastic isokinetic algorithm for multiple time-step molecular dynamics calculations using a polarizable model based on fluctuating dipoles. The scheme developed here employs two sets of induced dipole moments, specifically, those associated with short-range interactions and those associated with a full set of interactions. The scheme is demonstrated on the polarizable AMOEBA water model. As was seen with fixed-charge models, we are able to obtain large time steps exceeding 100 fs, allowing calculations to be performed 10 to 20 times faster than standard thermostated molecular dynamics.

Molecular-scale remnants of the liquid-gas transition in supercritical polar fluids

An electronically coarse-grained model for water reveals a persistent vestige of the liquid-gas transition deep into the supercritical region. A crossover in the density dependence of the molecular dipole arises from the onset of nonpercolating hydrogen bonds. The crossover points coincide with the Widom line in the scaling region but extend farther, tracking the heat capacity maxima, offering evidence for liquidlike and gaslike state points in a "one-phase" fluid. The effect is present even in dipole-limit models, suggesting that it is common for all molecular liquids exhibiting dipole enhancement in the liquid phase.

Surface-active organic matter induces salt morphology transitions during new atmospheric particle formation and growth

The salt within an aerosol nucleus assumes a brine morphology in increasing presence of organic matter on the surface. This affects, in turn, the water uptake dynamics.

Low polarity water, a novel transition species at the polyethylene–water interface

Polyethylene sandwich or a single window cell: dark bars represent polyethylene (PE) windows (double or single) of a cell. Wavy lines are water (W) and low polarity water (LPW). Subtraction of the single window spectrum from the double window spectrum leaves the LPW spectrum as illustrated in the figure.