Nuclear quantum effects in liquid water from path-integral simulations using anab initioforce-matching approach

Nuclear quantum effects on the vibrational dynamics of the water-air interface

We have applied path-integral molecular dynamics simulations to investigate the impact of nuclear quantum effects on the vibrational dynamics of water molecules at the water-air interface. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time series analysis, while the time scales of vibrational spectral diffusion are determined from frequency-time correlation functions and joint probability distributions. We find that the inclusion of nuclear quantum effects leads not only to a redshift in the vibrational frequency distribution by about 120 cm-1 for both the bulk and interfacial water molecules but also to an acceleration of the vibrational dynamics at the water-air interface by as much as 35%. In addition, a blueshift of about 45 cm-1 is seen in the vibrational frequency distribution of interfacial water molecules compared to that of the bulk. Furthermore, the dynamics of water molecules beyond the topmost molecular layer was found to be rather similar to that of bulk water.

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. The effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.

Accurate SCC-DFTB Parametrization of Liquid Water with Improved Atomic Charges and Iterative Boltzmann Inversion

This work presents improvements of the description of liquid water within the self-consistent-charge density-functional based tight-binding scheme combining the use of Weighted Mulliken (WMull) charges and optimized O-H repulsive potential through the iterative Boltzmann inversion (IBI) process. The quality of the newly developed models is validated considering pair radial distribution functions (RDFs), as well as other structural, energetic, thermodynamic, and dynamic properties. The use of WMull charges certainly improves the agreement with experimental data, however leading to over-structured RDFs at short distance, that can be further improved by considering an optimized O-H repulsive potential obtained by the IBI process. Three different schemes were used to optimize this potential: (i) optimization including short O-H distances. This led to accurate RDFs as well as improved self-diffusion coefficient and heat of vaporization, while the proton transfer energy barrier is severely deteriorated; (ii) optimization starting at long distance. The proton transfer energy barrier is recovered while the heat of vaporization is deteriorated and the O-H RDF is less accurate at short distance; (iii) optimization within the path-integral molecular dynamics scheme which allows us to exclude nuclear quantum effects from the repulsive potential. The latter potential, in conjunction with the WMull improved atomic charges, provides similar results as (i) for structural, dynamic, and thermodynamic properties while recovering a large part of the proton transfer energy barrier. It therefore offers a good compromise to study both dynamic properties and chemistry within liquid water at a quantum chemical level.

Dynamics of Aqueous Electrolyte Solutions: Challenges for Simulations

This Perspective article focuses on recent simulation work on the dynamics of aqueous electrolytes. It is well-established that full-charge, nonpolarizable models for water and ions generally predict solution dynamics that are too slow in comparison to experiments. Models with reduced (scaled) charges do better for solution diffusivities and viscosities but encounter issues describing other dynamic phenomena such as nucleation rates of crystals from solution. Polarizable models show promise, especially when appropriately parametrized, but may still miss important physical effects such as charge transfer. First-principles calculations are starting to emerge for these properties that are in principle able to capture polarization, charge transfer, and chemical transformations in solution. While direct ab initio simulations are still too slow for simulations of large systems over long time scales, machine-learning models trained on appropriate first-principles data show significant promise for accurate and transferable modeling of electrolyte solution dynamics.

Hydrogen bond dynamics of interfacial water molecules revealed from two-dimensional vibrational sum-frequency generation spectroscopy

Vibrational sum-frequency generation (vSFG) spectroscopy allows the study of the structure and dynamics of interfacial systems. In the present work, we provide a simple recipe, based on a narrowband IR pump and broadband vSFG probe technique, to computationally obtain the two-dimensional vSFG spectrum of water molecules at the air-water interface. Using this technique, to study the time-dependent spectral evolution of hydrogen-bonded and free water molecules, we demonstrate that at the interface, the vibrational spectral dynamics of the free OH bond is faster than that of the bonded OH mode.

Nuclear Quantum Effects from the Analysis of Smoothed Trajectories: Pilot Study for Water

Nuclear quantum effects have significant contributions to thermodynamic quantities and structural properties; furthermore, very expensive methods are necessary for their accurate computation. In most calculations, these effects, for instance, zero-point energies, are simply neglected or only taken into account within the quantum harmonic oscillator approximation. Herein, we present a new method, Generalized Smoothed Trajectory Analysis, to determine nuclear quantum effects from molecular dynamics simulations. The broad applicability is demonstrated with the examples of a harmonic oscillator and different states of water. Ab initio molecular dynamics simulations have been performed for ideal gas up to the temperature of 5000 K. Classical molecular dynamics have been carried out for hexagonal ice, liquid water, and vapor at atmospheric pressure. With respect to the experimental heat capacity, our method outperforms previous calculations in the literature in a wide temperature range at lower computational cost than other alternatives. Dynamic and structural nuclear quantum effects of water are also discussed.

Opposing Electronic and Nuclear Quantum Effects on Hydrogen Bonds in H2 O and D2 O

The effect of extending the O-H bond length(s) in water on the hydrogen-bonding strength has been investigated using static ab initio molecular orbital calculations. The "polar flattening" effect that causes a slight σ-hole to form on hydrogen atoms is strengthened when the bond is stretched, so that the σ-hole becomes more positive and hydrogen bonding stronger. In opposition to this electronic effect, path-integral ab initio molecular-dynamics simulations show that the nuclear quantum effect weakens the hydrogen bond in the water dimer. Thus, static electronic effects strengthen the hydrogen bond in H2 O relative to D2 O, whereas nuclear quantum effects weaken it. These quantum fluctuations are stronger for the water dimer than in bulk water.

Ab initio thermodynamics of liquid and solid water

Thermodynamic properties of liquid water as well as hexagonal (Ih) and cubic (Ic) ice are predicted based on density functional theory at the hybrid-functional level, rigorously taking into account quantum nuclear motion, anharmonic fluctuations, and proton disorder. This is made possible by combining advanced free-energy methods and state-of-the-art machine-learning techniques. The ab initio description leads to structural properties in excellent agreement with experiments and reliable estimates of the melting points of light and heavy water. We observe that nuclear-quantum effects contribute a crucial [Formula: see text] to the stability of ice Ih, making it more stable than ice Ic. Our computational approach is general and transferable, providing a comprehensive framework for quantitative predictions of ab initio thermodynamic properties using machine-learning potentials as an intermediate step.

Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost

It is now established that nuclear quantum motion plays an important role in determining water's hydrogen bonding, structure, and dynamics. Such effects are important to include in density functional theory (DFT) based molecular dynamics simulation of water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads required. In this work we introduce a method whereby PIMD can be incorporated into a DFT simulation with little extra cost and little loss in accuracy. The method is based on the many body expansion of the energy and has the benefit of including a monomer level correction to the DFT energy. Our method calculates intramolecular forces using the highly accurate monomer potential energy surface developed by Partridge-Schwenke, which is cheap to evaluate. Intermolecular forces and energies are calculated with DFT only once per timestep using the centroid positions. We show how our method may be used in conjunction with a multiple time step algorithm for an additional speedup and how it relates to ring polymer contraction and other schemes that have been introduced recently to speed up PIMD simulations. We show that our method, which we call "monomer PIMD", correctly captures changes in the structure of water found in a full PIMD simulation but at much lower computational cost.

Computational Methods for Ab Initio Molecular Dynamics

Ab initio molecular dynamics is an irreplaceable technique for the realistic simulation of complex molecular systems and processes from first principles. This paper proposes a comprehensive and self-contained review of ab initio molecular dynamics from a computational perspective and from first principles. Quantum mechanics is presented from a molecular dynamics perspective. Various approximations and formulations are proposed, including the Ehrenfest, Born–Oppenheimer, and Hartree–Fock molecular dynamics. Subsequently, the Kohn–Sham formulation of molecular dynamics is introduced as well as the afferent concept of density functional. As a result, Car–Parrinello molecular dynamics is discussed, together with its extension to isothermal and isobaric processes. Car–Parrinello molecular dynamics is then reformulated in terms of path integrals. Finally, some implementation issues are analysed, namely, the pseudopotential, the orbital functional basis, and hybrid molecular dynamics.

Harmonic-phase path-integral approximation of thermal quantum correlation functions

We present an approximation to the thermal symmetric form of the quantum time-correlation function in the standard position path-integral representation. By transforming to a sum-and-difference position representation and then Taylor-expanding the potential energy surface of the system to second order, the resulting expression provides a harmonic weighting function that approximately recovers the contribution of the phase to the time-correlation function. This method is readily implemented in a Monte Carlo sampling scheme and provides exact results for harmonic potentials (for both linear and non-linear operators) and near-quantitative results for anharmonic systems for low temperatures and times that are likely to be relevant to condensed phase experiments. This article focuses on one-dimensional examples to provide insights into convergence and sampling properties, and we also discuss how this approximation method may be extended to many-dimensional systems.

Nuclear quantum effects on the vibrational dynamics of liquid water

Based on quantum-mechanical path-integral molecular dynamics simulations, the impact of nuclear quantum effects on the vibrational and hydrogen bond dynamics in liquid water is investigated. The instantaneous fluctuations in the frequencies of the O-H stretch modes are calculated using the wavelet method of time-series analysis, while the time scales of the vibrational spectral diffusion are determined from frequency-time correlation functions, joint probability distributions, and the slope of three-pulse photon echo. We find that the inclusion of nuclear quantum effects leads not only to a redshift of the vibrational frequency distribution by around 130 cm^-1 but also to an acceleration of the vibrational dynamics by as much as 30%. In addition, quantum fluctuations also entail a significantly faster decay of correlation in the initial diffusive regime, which is in agreement with recent vibrational echo experiments.

Self-consistent molecular dynamics calculation of diffusion in higher n-alkanes

Diffusion is one of the key subjects of molecular modeling and simulation studies. However, there is an unresolved lack of consistency between Einstein-Smoluchowski (E-S) and Green-Kubo (G-K) methods for diffusion coefficient calculations in systems of complex molecules. In this paper, we analyze this problem for the case of liquid n-triacontane. The non-conventional long-time tails of the velocity autocorrelation function (VACF) are found for this system. Temperature dependence of the VACF tail decay exponent is defined. The proper inclusion of the long-time tail contributions to the diffusion coefficient calculation results in the consistency between G-K and E-S methods. Having considered the major factors influencing the precision of the diffusion rate calculations in comparison with experimental data (system size effects and force field parameters), we point to hydrogen nuclear quantum effects as, presumably, the last obstacle to fully consistent n-alkane description.

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.

Quantum ring-polymer contraction method: Including nuclear quantum effects at no additional computational cost in comparison to ab initio molecular dynamics

We present a simple and accurate computational method which facilitates ab initio path-integral molecular dynamics simulations, where the quantum-mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions using density functional theory. This development will enable routine inclusion of nuclear quantum effects in ab initio molecular dynamics simulations of condensed-phase systems.