Connection between water's dynamical and structural properties: Insights from ab initio simulations

Unraveling the Hydration Shell Structure and Dynamics of Group 10 Aqua Ions

We present ab initio simulations based on subsystem DFT of group 10 aqua ions accurately compared against experimental data on hydration structure. Our simulations provide insights into the molecular structures and dynamics of hydration shells, offering recalibrated interpretations of experimental results. We observe a soft, but distinct second hydration shell in Palladium (Pd) due to a balance between thermal fluctuations, metal-water interactions, and hydrogen bonding. Nickel (Ni) and platinum (Pt) exhibit more rigid hydration shells. Notably, our simulations align with experimental findings for Pd, showing axial hydration marked by a broad peak at about 3 Å in the Pd-O radial distribution function, revising the previously sharp "mesoshell" prediction. We introduce the "hydrogen bond dome" concept to describe a resilient network of hydrogen-bonded water molecules around the metal, which plays a critical role in the axial hydration dynamics.

Understanding X-ray absorption in liquid water using triple excitations in multilevel coupled cluster theory

X-ray absorption (XA) spectroscopy is an essential experimental tool to investigate the local structure of liquid water. Interpretation of the experiment poses a significant challenge and requires a quantitative theoretical description. High-quality theoretical XA spectra require reliable molecular dynamics simulations and accurate electronic structure calculations. Here, we present the first successful application of coupled cluster theory to model the XA spectrum of liquid water. We overcome the computational limitations on system size by employing a multilevel coupled cluster framework for large molecular systems. Excellent agreement with the experimental spectrum is achieved by including triple excitations in the wave function and using molecular structures from state-of-the-art path-integral molecular dynamics. We demonstrate that an accurate description of the electronic structure within the first solvation shell is sufficient to successfully model the XA spectrum of liquid water within the multilevel framework. Furthermore, we present a rigorous charge transfer analysis of the XA spectrum, which is reliable due to the accuracy and robustness of the electronic structure methodology. This analysis aligns with previous studies regarding the character of the prominent features of the XA spectrum of liquid water.

Structure and dynamics of liquid water from ab initio simulations: adding Minnesota density functionals to Jacob's ladder

The accurate representation of the structural and dynamical properties of water is essential for simulating the unique behavior of this ubiquitous solvent. Here we assess the current status of describing liquid water using ab initio molecular dynamics, with a special focus on the performance of all the later generation Minnesota functionals. Findings are contextualized within the current knowledge on DFT for describing bulk water under ambient conditions and compared to experimental data. We find that, contrary to the prevalent idea that local and semilocal functionals overstructure water and underestimate dynamical properties, M06-L, revM06-L, and M11-L understructure water, while MN12-L and MN15-L overdistance water molecules due to weak cohesive effects. This can be attributed to a weakening of the hydrogen bond network, which leads to dynamical fingerprints that are over fast. While most of the hybrid Minnesota functionals (M06, M08-HX, M08-SO, M11, MN12-SX, and MN15) also yield understructured water, their dynamical properties generally improve over their semilocal counterparts. It emerges that exact exchange is a crucial component for accurately describing hydrogen bonds, which ultimately leads to corrections in both the dynamical and structural properties. However, an excessive amount of exact exchange strengthens hydrogen bonds and causes overstructuring and slow dynamics (M06-HF). As a compromise, M06-2X is the best performing Minnesota functional for water, and its D3 corrected variant shows very good structural agreement. From previous studies considering nuclear quantum effects (NQEs), the hybrid revPBE0-D3, and the rung-5 RPA (RPA@PBE) have been identified as the only two approximations that closely agree with experiments. Our results suggest that the M06-2X(-D3) functionals have the potential to further improve the reproduction of experimental properties when incorporating NQEs through path integral approaches. This work provides further proof that accurate modeling of water interactions requires the inclusion of both exact exchange and balanced (non-local) correlation, highlighting the need for higher rungs on Jacob's ladder to achieve predictive simulations of complex biological systems in aqueous environments.

Structural effects of water clusters on viscosity at high shear rates

In this study, we use molecular dynamics simulations of liquid water to investigate how shear thinning affects the viscosity of liquid water by structural changes of the hydrogen bond network. The effect of shear on viscosity can be divided into two parts: shear-induced destruction of the hydrogen bond network and the influence of the water structure on shear viscosity. First, strong shear destroys tetrahedral structures and thus reduces the connectivity of the hydrogen bond network. It is mainly because shear deformation, characterized by compression and expansion axes, respectively, triggers the destruction and formation of hydrogen bonds, resulting in anisotropic effects on water structures. At the same time, shear destroys large clusters and enhances the formation of small ones, resulting in a decrease in average cluster sizes. Second, the change of viscosity obeys a power law relationship with the change of hydrogen bond structures, highlighting a one-to-one correspondence between structure and property. Meanwhile, in order to explain why the structure affects viscosity, we define hydrogen-bond viscosity and find that the cooperative motion of the water structures can promote momentum transfer in the form of aggregations. Hydrogen-bond viscosity accounts for 5%-50% of the total viscosity. Our results elucidate that water structures are the important structural units to explain the change of water properties.

Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations

Molecular dynamics (MD) simulations are widely used to investigate molecular systems at atomic resolution including biomolecular structures, drug-receptor interactions, and novel materials. Frequently, MD simulations are performed in an aqueous solution with explicit models of water molecules. Commonly, such models are parameterized to reproduce the liquid phase of water under ambient conditions. However, often, simulations at significantly higher temperatures are also of interest. Hence, it is important to investigate the equilibrium of the liquid and vapor phases of molecular models of water at elevated temperatures. Here, we evaluate the behavior of 11 common rigid three-point water models over a wide range of temperatures. From liquid-vapor coexistence simulations, we estimated the critical points and studied the spontaneous evaporation of these water models. Moreover, we investigated the influence of the system size, choice of the pressure-coupling algorithm, and rate of heating on the process and compared them with the experimental data. We found that modern rigid three-point water models reproduce the critical point surprisingly well. Furthermore, we discovered that the critical temperature correlates with the quadrupole moment of the respective water model. This indicates that the spatial arrangement of the partial charges is important for reproducing the liquid-vapor phase transition. Our findings may guide the selection of water models for simulations conducted at high temperatures.

Transport properties in liquids from first-principles: The case of liquid water and liquid argon

Shear and bulk viscosities of liquid water and argon are evaluated from first-principles in the density functional theory (DFT) framework, by performing molecular dynamics simulations in the NVE ensemble and using the Kubo–Greenwood equilibrium approach. The standard DFT functional is corrected in such a way to allow for a reasonable description of van der Waals effects. For liquid argon, the thermal conductivity has been also calculated. Concerning liquid water, to our knowledge, this is the first estimate of the bulk viscosity and of the shear-viscosity/bulk-viscosity ratio from first-principles. By analyzing our results, we can conclude that our first-principles simulations, performed at a nominal average temperature of 366 to guarantee that the systems are liquid-like, actually describe the basic dynamical properties of liquid water at about 330 K. In comparison with liquid water, the normal, monatomic liquid Ar is characterized by a much smaller bulk-viscosity/shear-viscosity ratio (close to unity) and this feature is well reproduced by our first-principles approach, which predicts a value of the ratio in better agreement with experimental reference data than that obtained using the empirical Lennard-Jones potential. The computed thermal conductivity of liquid argon is also in good agreement with the experimental value.