‘‘Ab initio’’ liquid water

Comparing machine learning potentials for water: Kernel-based regression and Behler-Parrinello neural networks

In this paper, we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE + D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1500 structures, as well as a smaller dataset, about half the size, obtained using only on-the-fly learning. This study reveals that despite minor differences between the MLPs, their agreement on observables such as the diffusion constant and pair-correlation functions is excellent, especially for the large training dataset. Variations in the predicted density isobars, albeit somewhat larger, are also acceptable, particularly given the errors inherent to approximate density functional theory. Overall, this study emphasizes the relevance of the database over the fitting method. Finally, this study underscores the limitations of root mean square errors and the need for comprehensive testing, advocating the use of multiple MLPs for enhanced certainty, particularly when simulating complex thermodynamic properties that may not be fully captured by simpler tests.

Consistent density functional theory-based description of ion hydration through density-corrected many-body representations

Delocalization error constrains the accuracy of density functional theory in describing molecular interactions in ion-water systems. Using Na+ and Cl- in water as model systems, we calculate the effects of delocalization error in the SCAN functional for describing ion-water and water-water interactions in hydrated ions, and demonstrate that density-corrected SCAN (DC-SCAN) predicts n-body and interaction energies with an accuracy approaching coupled cluster theory. The performance of DC-SCAN is size-consistent, maintaining an accurate description of molecular interactions well beyond the first solvation shell. Molecular dynamics simulations at ambient conditions with many-body MB-SCAN(DC) potentials, derived from the many-body expansion, predict the solvation structure of Na+ and Cl- in quantitative agreement with reference data, while simultaneously reproducing the structure of liquid water. Beyond rationalizing the accuracy of density-corrected models of ion hydration, our findings suggest that our unified density-corrected MB formalism holds great promise for efficient DFT-based simulations of condensed-phase systems with chemical accuracy.

Why Dissolving Salt in Water Decreases Its Dielectric Permittivity

The dielectric permittivity of salt water decreases on dissolving more salt. For nearly a century, this phenomenon has been explained by invoking saturation in the dielectric response of the solvent water molecules. Herein, we employ an advanced deep neural network (DNN), built using data from density functional theory, to study the dielectric permittivity of sodium chloride solutions. Notably, the decrease in the dielectric permittivity as a function of concentration, computed using the DNN approach, agrees well with experiments. Detailed analysis of the computations reveals that the dominant effect, caused by the intrusion of ionic hydration shells into the solvent hydrogen-bond network, is the disruption of dipolar correlations among water molecules. Accordingly, the observed decrease in the dielectric permittivity is mostly due to increasing suppression of the collective response of solvent waters.

Ab initio study of water dissociation on a charged Pd(111) surface

The interactions between molecules and electrode surfaces play a key role in electrochemical processes and are a subject of extensive research, both experimental and theoretical. In this paper, we address the water dissociation reaction on a Pd(111) electrode surface, modeled as a slab embedded in an external electric field. We aim at unraveling the relationship between surface charge and zero-point energy in aiding or hindering this reaction. We calculate the energy barriers with dispersion-corrected density-functional theory and an efficient parallel implementation of the nudged-elastic-band method. We show that the lowest dissociation barrier and consequently the highest reaction rate take place when the field reaches a strength where two different geometries of the water molecule in the reactant state are equally stable. The zero-point energy contributions to this reaction, on the other hand, remain nearly constant across a wide range of electric field strengths, despite significant changes in the reactant state. Interestingly, we show that the application of electric fields that induce a negative charge on the surface can make nuclear tunneling more significant for these reactions.

Ab Initio Cluster Approach for High Harmonic Generation in Liquids

High harmonic generation (HHG) takes place in all phases of matter. In gaseous atomic and molecular media, it has been extensively studied and is very well understood. In solids, research is ongoing, but a consensus is forming for the dominant microscopic HHG mechanisms. In liquids, on the other hand, no established theory yet exists, and approaches developed for gases and solids are generally inapplicable, hindering our current understanding. We develop here a powerful and reliable ab initio cluster-based approach for describing the nonlinear interactions between isotropic bulk liquids and intense laser pulses. The scheme is based on time-dependent density functional theory and utilizes several approximations that make it feasible yet accurate in realistic systems. We demonstrate our approach with HHG calculations in water, ammonia, and methane liquids and compare the characteristic response of polar and nonpolar liquids. We identify unique features in the HHG spectra of liquid methane that could be utilized for ultrafast spectroscopy of its chemical and physical properties, including a structural minimum at 15-17 eV that is associated solely with the liquid phase. Our results pave the way to accessible calculations of HHG in liquids and illustrate the unique nonlinear nature of liquid systems.

A rationale for non-linear responses to strong electric fields in molecular dynamics simulations

Many approaches for calculation of the field-dependent electric properties of water solutions rely on the Onsager and Kirkwood theories of polar dielectrics. Such basic theories implicitly consider the electric field intensity to fulfill the so-called 'weak field conditions', i.e. to produce a linear response in the system. In this work we made use of molecular dynamics simulations to investigate possible non-linear effects induced by high intensity electric fields, specifically continuous wave bursts with nanosecond duration, comparing them with the ones predicted by the theory. We found that field intensities above 0.15 V nm^-1 produce remarkable nonlinear responses in the whole 100 MHz-100 GHz frequency window considered, with the onset of higher order polarization signals, which are the clear fingerprint of harmonic distorsions. That non-linear response turned out to depend on the considered frequency. We finally show that MD outcomes are consistent with a modelization based on an extended formulation of the Langevin function including a frequency-dependent parameter.

Assessing the Interplay between Functional-Driven and Density-Driven Errors in DFT Models of Water

We investigate the interplay between functional-driven and density-driven errors in different density functional approximations within density functional theory (DFT) and the implications of these errors for simulations of water with DFT-based data-driven potentials. Specifically, we quantify density-driven errors in two widely used dispersion-corrected functionals derived within the generalized gradient approximation (GGA), namely BLYP-D3 and revPBE-D3, and two modern meta-GGA functionals, namely strongly constrained and appropriately normed (SCAN) and B97M-rV. The effects of functional-driven and density-driven errors on the interaction energies are first assessed for the water clusters of the BEGDB dataset. Further insights into the nature of functional-driven errors are gained from applying the absolutely localized molecular orbital energy decomposition analysis (ALMO-EDA) to the interaction energies, which demonstrates that functional-driven errors are strongly correlated with the nature of the interactions. We discuss cases where density-corrected DFT (DC-DFT) models display higher accuracy than the original DFT models and cases where reducing the density-driven errors leads to larger deviations from the reference energies due to the presence of large functional-driven errors. Finally, molecular dynamics simulations are performed with data-driven many-body potentials derived from DFT and DC-DFT data to determine the effect that minimizing density-driven errors has on the description of liquid water. Besides rationalizing the performance of widely used DFT models of water, we believe that our findings unveil fundamental relations between the shortcomings of some common DFT approximations and the requirements for accurate descriptions of molecular interactions, which will aid the development of a consistent, DFT-based framework for the development of data-driven and machine-learned potentials for simulations of condensed-phase systems.

Density Functional Theory of Water with the Machine-Learned DM21 Functional

The delicate interplay between functional-driven and density-driven errors in density functional theory (DFT) has hindered traditional density functional approximations (DFAs) from providing an accurate description of water for over 30 years. Recently, the deep-learned DeepMind 21 (DM21) functional has been shown to overcome the limitations of traditional DFAs as it is free of delocalization error. To determine if DM21 can enable a molecular-level description of the physical properties of aqueous systems within Kohn-Sham DFT, we assess the accuracy of the DM21 functional for neutral, protonated, and deprotonated water clusters. We find that the ability of DM21 to accurately predict the energetics of aqueous clusters varies significantly with cluster size. Additionally, we introduce the many-body MB-DM21 potential derived from DM21 data within the many-body expansion of the energy and use it in simulations of liquid water as a function of temperature at ambient pressure. We find that size-dependent functional-driven errors identified in the analysis of the energetics of small clusters calculated with the DM21 functional result in the MB-DM21 potential systematically overestimating the hydrogen-bond strength and, consequently, predicting a more ice-like local structure of water at room temperature.

Electrocatalysis in Alkaline Media and Alkaline Membrane-Based Energy Technologies

Hydrogen energy-based electrochemical energy conversion technologies offer the promise of enabling a transition of the global energy landscape from fossil fuels to renewable energy. Here, we present a comprehensive review of the fundamentals of electrocatalysis in alkaline media and applications in alkaline-based energy technologies, particularly alkaline fuel cells and water electrolyzers. Anion exchange (alkaline) membrane fuel cells (AEMFCs) enable the use of nonprecious electrocatalysts for the sluggish oxygen reduction reaction (ORR), relative to proton exchange membrane fuel cells (PEMFCs), which require Pt-based electrocatalysts. However, the hydrogen oxidation reaction (HOR) kinetics is significantly slower in alkaline media than in acidic media. Understanding these phenomena requires applying theoretical and experimental methods to unravel molecular-level thermodynamics and kinetics of hydrogen and oxygen electrocatalysis and, particularly, the proton-coupled electron transfer (PCET) process that takes place in a proton-deficient alkaline media. Extensive electrochemical and spectroscopic studies, on single-crystal Pt and metal oxides, have contributed to the development of activity descriptors, as well as the identification of the nature of active sites, and the rate-determining steps of the HOR and ORR. Among these, the structure and reactivity of interfacial water serve as key potential and pH-dependent kinetic factors that are helping elucidate the origins of the HOR and ORR activity differences in acids and bases. Additionally, deliberately modulating and controlling catalyst-support interactions have provided valuable insights for enhancing catalyst accessibility and durability during operation. The design and synthesis of highly conductive and durable alkaline membranes/ionomers have enabled AEMFCs to reach initial performance metrics equal to or higher than those of PEMFCs. We emphasize the importance of using membrane electrode assemblies (MEAs) to integrate the often separately pursued/optimized electrocatalyst/support and membranes/ionomer components. Operando/in situ methods, at multiscales, and ab initio simulations provide a mechanistic understanding of electron, ion, and mass transport at catalyst/ionomer/membrane interfaces and the necessary guidance to achieve fuel cell operation in air over thousands of hours. We hope that this Review will serve as a roadmap for advancing the scientific understanding of the fundamental factors governing electrochemical energy conversion in alkaline media with the ultimate goal of achieving ultralow Pt or precious-metal-free high-performance and durable alkaline fuel cells and related technologies.

Advances in the study of supercooled water

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

The Fermi level as an energy reference in liquid jet X-ray photoelectron spectroscopy studies of aqueous solutions

To advance the understanding of key electrochemical and photocatalytic processes that depend on the electronic structure of aqueous solutions, X-ray photoemission spectroscopy has become an invaluable tool, especially when practiced with liquid microjet setups. Determining vertical ionization energies referenced to the vacuum level, and binding energies referenced to the Fermi level, including the much-coveted reorganization energy of the oxidized species of a redox couple, requires that energy levels be properly defined. The present paper addresses specifically how the vacuum level "just outside the surface" can be known through the energy position of the rising edge of the secondary electrons, and how the Fermi level reference is uniquely determined via the introduction of a redox couple. Taking the case of the ferricyanide/ferrocyanide and ferric/ferrous couples, this study also tackles issues related to the electrokinetic effects inherent to the production of a liquid jet in a vacuum, which has become the standard water sample environment for photoemission experiments.

Nuclear Quantum Effect and Its Temperature Dependence in Liquid Water from Random Phase Approximation via Artificial Neural Network

We report structural and dynamical properties of liquid water described by the random phase approximation (RPA) correlation together with the exact exchange energy (EXX) within density functional theory. By utilizing thermostated ring polymer molecular dynamics, we examine the nuclear quantum effects and their temperature dependence. We circumvent the computational limitation of performing direct first-principles molecular dynamics simulation at this high level of electronic structure theory by adapting an artificial neural network model. We show that the EXX+RPA level of theory accurately describes liquid water in terms of both dynamical and structural properties.

Dielectric Constant of Liquid Water Determined with Neural Network Quantum Molecular Dynamics

The static dielectric constant ϵ_{0} and its temperature dependence for liquid water is investigated using neural network quantum molecular dynamics (NNQMD). We compute the exact dielectric constant in canonical ensemble from NNQMD trajectories using fluctuations in macroscopic polarization computed from maximally localized Wannier functions (MLWF). Two deep neural networks are constructed. The first, NNQMD, is trained on QMD configurations for liquid water under a variety of temperature and density conditions to learn potential energy surface and forces and then perform molecular dynamics simulations. The second network, NNMLWF, is trained to predict locations of MLWF of individual molecules using the atomic configurations from NNQMD. Training data for both the neural networks is produced using a highly accurate quantum-mechanical method, DFT-SCAN that yields an excellent description of liquid water. We produce 280×10^{6} configurations of water at 7 temperatures using NNQMD and predict MLWF centers using NNMLWF to compute the polarization fluctuations. The length of trajectories needed for a converged value of the dielectric constant at 0°C is found to be 20 ns (40×10^{6} configurations with 0.5 fs time step). The computed dielectric constants for 0, 15, 30, 45, 60, 75, and 90°C are in good agreement with experiments. Our scalable scheme to compute dielectric constants with quantum accuracy is also applicable to other polar molecular liquids.

Understanding the Dipole Moment of Liquid Water from a Self-Attractive Hartree Decomposition

The dipole moment of a single water molecule in liquid water has been a critical concept for understanding water's dielectric properties. In this work, we investigate the dipole moment of liquid water through a self-attractive Hartree (SAH) decomposition of total electron density computed by density functional theory, on water clusters sampled from ab initio molecular dynamics simulation of bulk water. By adjusting one parameter that controls the degree of density localization, we reveal two distinct pictures of water dipoles that are consistent with bulk dielectric properties: a localized picture with smaller and less polarizable monomer dipoles and a delocalized picture with larger and more polarizable monomer dipoles. We further uncover that the collective dipole-dipole correlation is stronger in the localized picture and is key to connecting individual dipoles with bulk dielectric properties. On the basis of these findings, we suggest considering both individual and collective dipole behaviors when studying the dipole moment of liquid water and propose new design strategies for developing water models.

Born–Oppenheimer molecular dynamics simulations on structures of high-density and low-density water: a comparison of the SCAN meta-GGA and PBE GGA functionals

Using Born–Oppenheimer ab initio molecular dynamics (BOAIMD) simulations, the high-density water (HDW) and low-density water (LDW) structures based on SCAN meta-GGA are compared with those based on PBE GGA.

Simulation of Liquids with the Tight-Binding Density-Functional Approach and Improved Atomic Charges

Theoretical description of liquids, especially liquid water, is an ongoing subject with important implications in various domains such as homogeneous catalysis; solvation of molecular, ionic, and biomolecular species; and reactivity. Various formalisms exist to describe liquids, each one displaying its own balance between accuracy and computational cost that defines its range of applications. The present article revisits the ability of the density-functional-based tight-binding (SCC-DFTB) approach to model liquids by focusing on liquid water and liquid benzene under ambient conditions. To do so, we benchmark a recent correction for the SCC-DFTB atomic charges that allows for a drastic improvement of the pair radial distribution functions of liquid water as compared to both experimental data and density-functional theory results performed in the generalized-gradient approximation. We also report the coupling of the deMonNano and i-PI codes to perform path-integral molecular dynamics. This allows us to rationalize the impact of nuclear quantum effects on the SCC-DFTB description of liquid water. This study evidences the rather good ability of SCC-DFTB to describe liquid water and liquid benzene. As the first example of application, we also present results for a benzene molecule solvated in water with the perspectives of further studies devoted to solvent/water interfaces.

Frontiers in molecular simulation of solvated ions, molecules and interfaces

This themed collection is a collection of articles on frontiers in molecular simulation of solvated ions, molecules and interfaces.

Hydrogen Bonding in the Liquid State of Linear Alcohols: Molecular Dynamics and Thermodynamics

Linear monohydroxy alcohols are strongly hydrogen-bonded liquids that are considered to be homologues of water. Here, we report ab initio molecular dynamics simulations of the liquid alcohols, methanol to pentanol, and from the combined radial-angular probability distribution of the intermolecular O···O distances and HO···O angles determine the geometrical parameters that define the hydrogen bonds in these systems. The key feature of hydrogen bonds in the liquid alcohols, irrespective of the size of the alkyl group, is the strong orientation dependence with the donor-acceptor HO···O angle being close to zero, similar to that observed in liquid water. Hydrogen bond formation is consequently considered to be the passage from a state where donor-acceptor pairs show no preferred orientation to one where they are almost linear. The potential of mean force, the reversible work associated with this process, is computed from the pair probability density distributions obtained from the simulations and that for a hypothetical state where donor-acceptor pairs are randomly oriented. We find that the magnitude of the free energy for hydrogen bond formation is maximum for ethanol and show that this arises from a larger electrostatic contribution to hydrogen bond formation in ethanol as compared to the other alcohols.

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.

A comparative study on Monte Carlo simulations of electron emission from liquid water

PURPOSE

Liquid water being the major constituent of the human body, is of fundamental importance in radiobiological research. Hence, the knowledge of electron-water interaction physics and particularly the secondary electron yield is essential. However, to date, only very little is known experimentally on the low energy electron interaction with liquid water because of certain practical limitations. The purpose of this study was to gain some useful information about electron emission from water using a Monte Carlo (MC) simulation technique that can numerically model electron transport trajectories in water.

METHODS

In this study, we have performed MC simulations of electron emission from liquid water in the primary energy range of 50 eV-30 keV by using two different codes, i.e., a classical trajectory MC (CMC) code developed in our laboratory and the Geant4-DNA (G4DNA) code. The calculated secondary electron yield and electron backscattering coefficient are compared with experimental results wherever applicable to verify the validity of physical models for the electron-water interaction.

RESULTS

The secondary electron yield vs. primary energy curves calculated using the two codes present the same generic curve shape as that of metals but in rather different absolute values. G4DNA underestimates the secondary electron yield due to the application of one step thermalization model by setting a cutoff energy at 10 eV so that the low energy losses due to phonon excitations are omitted. Our CMC code, using a full energy loss spectrum to model electron inelastic scattering, allows the simulation of individual phonon scattering events for very low energy losses down to 10 meV, which then enables the calculated secondary electron yields much closer to the experimental data and also gives quite reasonable energy distribution curve of secondary electrons.

CONCLUSIONS

It is concluded that full dielectric function data at low energy loss values below 10 eV are recommended for modeling of low energy electrons in liquid water.

First-Principles Study on the Peculiar Water Environment in a Hydrate-Melt Electrolyte

Aqueous electrolytes have great potential to improve the safety and production costs of Li-ion batteries. Our recent materials exploration led to the discovery of the Li-salt dihydrate melt Li(TFSI)_0.7(BETI)_0.3·2H₂O, which possesses an extremely wide potential window. To clarify the detailed liquid structure and electronic states of this unique aqueous system, a first-principles molecular dynamics study has been conducted. We found that water molecules in the hydrate melt exist as isolated monomers or clusters consisting of only a few (at most five) H₂O molecules. Both the monomers and the clusters have electronic structures largely deviating from that in bulk water, where the lowest unoccupied states are higher in energy than that of the Li-salt anions, which preferentially cause anion reduction leading to formation of an anion-derived stable solid-electrolyte interphase. This clearly shows the role of characteristic electronic structure inherent to the peculiar water environment for the extraordinary electrochemical stability of hydrate melts.

Car-Parrinello Monitor for More Robust Born-Oppenheimer Molecular Dynamics

Born-Oppenheimer molecular dynamics (BOMD) is a promising simulation method for exploring the possible reaction pathways of a chemical system, but one significant challenge is the increased difficulty of converging the self-consistent field (SCF) calculation that often accompanies the breaking and forming of chemical bonds. To address this challenge, we developed an enhancement to the BOMD simulation method called the Car-Parrinello monitor (CPMonitor) that uses Car-Parrinello molecular dynamics (CPMD) to recover from SCF convergence failures. CPMonitor works by detecting SCF convergence failures in BOMD and switching to a CPMD Hamiltonian to propagate through the region of configuration space where the SCF calculation is unable to converge, then switching back to BOMD when good convergence behavior is re-established. We present a series of simulation studies that use CPMonitor, including detailed studies of the thermodynamic and dynamical properties of simple systems, as well as ab initio nanoreactor simulations containing transition metal atoms that were previously not possible to simulate using standard BOMD methods. Our studies show that CPMonitor can make BOMD simulations robust to SCF convergence difficulties and improve simulation performance and stability in reaction discovery applications.

An ab initio molecular dynamics study of benzene in water at supercritical conditions: Structure, dynamics, and polarity of hydration shell water and the solute

Anisotropic structure and dynamics of the hydration shell of a benzene solute in supercritical water are investigated by means of ab initio molecular dynamics simulations. The polarity and structural distortion of the benzene solute in supercritical water are also investigated in this study. Calculations are done at 673 K for three different densities of the solvent. The simulations are carried out using the Becke-Lee-Yang-Parr (BLYP) and also the Becke-Lee-Yang-Parr functional including dispersion corrections of Grimme (BYLP-D). The structural anisotropy is found to exist even at supercritical conditions as elucidated by the radial distribution functions of different conical regions and also by angular and spatial distribution functions. The benzene-water πH-bond and also the water-water hydrogen bonds are found to exist even at the supercritical temperature of 673 K. However, the numbers of these hydrogen bonds are reduced substantially with a decrease in water density. The water molecules in the axial region of benzene are found to be preferably oriented with one OH vector pointing toward the benzene ring, whereas the water molecules located in the equatorial region are found to orient their dipoles mostly parallel to the ring plane. The orientational distributions, however, are found to be rather broad at the supercritical temperature due to thermal fluctuations. Although the water molecules have faster dynamics at these supercritical conditions, a slight difference is observed in the dynamics of the solvation shell and bulk molecules. The conformational flexibility of the ring is found to be enhanced which causes an increase in polarity of the benzene solute in water under supercritical conditions.

Hydrogen storage mechanism and diffusion in metal-organic frameworks

Diffusion and storage of hydrogen molecules in metal-organic frameworks are crucial for the development of next-generation energy storage devices. By resorting to the first principles modeling, we compute the diffusion coefficient of molecular hydrogen in these systems in a range of temperatures where MOF-based devices are expected to operate. The explicit inclusion of the electronic structure shows that diffusivities are one order of magnitude smaller than those reported by classical simulations, evidencing the insufficiency of the empirical force fields used so far. We show that hydrogen is mainly rolled up around the metal oxide nodes both in MOF-5 and IRMOF-6, and partly around the carbon atoms in the case of IRMOF-6, where charged linkers are present. Metal ions embedded in the junction sites exert an electrostatic attraction toward hydrogen and the resulting distribution shows some ordering around these same sites at low temperature, whereas this tendency vanishes at room temperature. The induced polarization of hydrogen molecules generates an electrostatic interaction with charged atoms inside these nano-scaffolds and this is a key factor for the enhancement in hydrogen storage both in MOF-5 and IRMOF-6. The mechanism discussed hereby provides a novel understanding of metal-organic frameworks and acts as a guide to tune their efficiency for hydrogen storage. Moreover it paves the way to a computer-aided design of effective MOFs indicating that a fine control of the distribution of electrostatic charges inside the hydrogen hosting structure is crucial.

Density anomaly of water at negative pressures from first principles

Using molecular dynamics simulations based on ab initio trained high-dimensional neural network potentials, we study the equation of state of liquid water at negative pressures. From density isobars computed for various pressures down to p  =  -230 MPa we determine the line of density maxima for two potentials based on the BLYP and the RPBE functionals, respectively. In both cases, dispersion corrections are included to account for non-local long-range correlations that give rise to van der Waals forces. We have followed the density maximum down to negative pressures close to the spinodal instability. For both functionals, the temperature of maximum density increases with decreasing pressure under moderate stretching, but changes slope at [Formula: see text] MPa and [Formula: see text] MPa for BLYP and RPBE, respectively. Our calculations confirm qualitatively the retracing shape of the line of density maxima found for empirical water models, indicating that the spinodal line maintains a positive slope even at strongly negative pressures.

Influence of the hydrogen bond quantum nature in liquid water and heavy water on stimulated Raman scattering

Stimulated Raman scattering (SRS) of liquid water and heavy water have been investigated using Nd:YAG laser. The SRS spectra of liquid heavy water indicate that ice-VII and ice-VIII structures are formed by shock-induced compression (SIC) in forward and backward directions, respectively. Simultaneously, the SRS spectra reveal of liquid water that only ice-VII structure is formed in the backward direction. The difference in ice structures formed by SIC in liquid water and heavy water could be attributed to the effect of the hydrogen bond quantum nature with H⁺. SRS spectra of 2M NaOH water solution with ice-VII and ice-VIII structures have been successfully obtained in forward and backward, respectively, as OH^- greatly reduce the quantum nature of hydrogen bonds by neutralizing H⁺ in water. The hydrogen bond quantum nature is important for understanding isotope calibration test structure and isotopic effect.

Equilibration and analysis of first-principles molecular dynamics simulations of water

First-principles molecular dynamics (FPMD) simulations based on density functional theory are becoming increasingly popular for the description of liquids. In view of the high computational cost of these simulations, the choice of an appropriate equilibration protocol is critical. We assess two methods of estimation of equilibration times using a large dataset of first-principles molecular dynamics simulations of water. The Gelman-Rubin potential scale reduction factor [A. Gelman and D. B. Rubin, Stat. Sci. 7, 457 (1992)] and the marginal standard error rule heuristic proposed by White [Simulation 69, 323 (1997)] are evaluated on a set of 32 independent 64-molecule simulations of 58 ps each, amounting to a combined cumulative time of 1.85 ns. The availability of multiple independent simulations also allows for an estimation of the variance of averaged quantities, both within MD runs and between runs. We analyze atomic trajectories, focusing on correlations of the Kohn-Sham energy, pair correlation functions, number of hydrogen bonds, and diffusion coefficient. The observed variability across samples provides a measure of the uncertainty associated with these quantities, thus facilitating meaningful comparisons of different approximations used in the simulations. We find that the computed diffusion coefficient and average number of hydrogen bonds are affected by a significant uncertainty in spite of the large size of the dataset used. A comparison with classical simulations using the TIP4P/2005 model confirms that the variability of the diffusivity is also observed after long equilibration times. Complete atomic trajectories and simulation output files are available online for further analysis.

First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional

We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen-oxygen partial correlation function at the experimental equilibrium density and 311 ± 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water.

Hydrogen-bond structure dynamics in bulk water: insights from ab initio simulations with coupled cluster theory

An accurate and efficient ab initio molecular dynamics (AIMD) simulation of liquid water was made possible using the fragment-based approach (J. F. Liu, X. He and J. Z. H. Zhang, Phys. Chem. Chem. Phys., 2017, 19, 11931-11936). In this study, we advance the AIMD simulations using the fragment-based coupled cluster (CC) theory, more accurately revealing the structural and dynamical properties of liquid water under ambient conditions. The results show that the double-donor hydrogen-bond configurations in liquid water are nearly in balance with the single-donor configurations, with a slight bias towards the former. Our observation is in contrast to the traditional tetrahedral water structure. The hydrogen-bond switching dynamics in liquid water are very fast, with a hydrogen-bond life time of around 0.78 picoseconds, determined using AIMD simulation at the CCD/aug-cc-pVDZ level. This time scale is remarkably shorter than the ∼3.0 picoseconds that is commonly obtained from traditional nonpolarized force fields and density functional theory (DFT) based first-principles simulations. Additionally, the obtained radial distribution functions, triplet oxygen angular distribution, diffusion coefficient, and the dipole moment of the water molecule are uniformly in good agreement with the experimental observations. The current high-level AIMD simulation sheds light on the understanding of the structural and dynamical properties of liquid water.

Structure of liquid water - a dynamical mixture of tetrahedral and 'ring-and-chain' like structures

The nature of the dynamical hydrogen-bond network of liquid water under ambient conditions has challenged both experimental and theoretical researchers for decades and remains a topic of intense debate. In this work, we addressed the structural issue of the hydrogen-bond network of liquid water based on an accurate ab initio molecular dynamics simulation. The present work showed clearly that liquid water is neither accurately described by a static picture of mostly tetrahedral water molecules nor dominated by "ring-and-chain" like structures. Instead, the structure of water is a dynamical mixture of tetrahedral and 'ring-and-chain' like structures with a slight bias toward the former. On average, each water molecule forms about three hydrogen bonds with the surrounding water molecules. The present accurate ab initio molecular dynamics simulation of liquid water was made possible by using a fragment-based second-order Møller-Plesset perturbation theory (MP2) with a large basis set to treat a large body of water molecules. This level of ab initio theory is sufficiently accurate for describing water interactions, and the simulated structural and dynamical properties of liquid water, including radial distribution functions, diffusion coefficient, dipole moment, etc., are uniformly in excellent agreement with experimental observations.