Heat capacity of water: A signature of nuclear quantum effects

Accurate nuclear quantum statistics on machine-learned classical effective potentials

The contribution of nuclear quantum effects (NQEs) to the properties of various hydrogen-bound systems, including biomolecules, is increasingly recognized. Despite the development of many acceleration techniques, the computational overhead of incorporating NQEs in complex systems is sizable, particularly at low temperatures. In this work, we leverage deep learning and multiscale coarse-graining techniques to mitigate the computational burden of path integral molecular dynamics (PIMD). In particular, we employ a machine-learned potential to accurately represent corrections to classical potentials, thereby significantly reducing the computational cost of simulating NQEs. We validate our approach using four distinct systems: Morse potential, Zundel cation, single water molecule, and bulk water. Our framework allows us to accurately compute position-dependent static properties, as demonstrated by the excellent agreement obtained between the machine-learned potential and computationally intensive PIMD calculations, even in the presence of strong NQEs. This approach opens the way to the development of transferable machine-learned potentials capable of accurately reproducing NQEs in a wide range of molecular systems.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

Unmasking quantum effects in the surface thermodynamics of fluid nanodrops

The focus of our study is an in-depth investigation of the quantum effects associated with the surface tension and other thermodynamic properties of nanoscopic liquid drops. The behavior of drops of quantum Lennard-Jones fluids is investigated with path-integral Monte Carlo simulations, and the test-area method is used to determine the surface tension of the spherical vapor-liquid interface. As the thermal de Broglie wavelength, λB, becomes more significant, the average density of the liquid drop decreases, with the drop becoming mechanically unstable at large wavelengths. As a consequence, the surface tension is found to decrease monotonically with λB, vanishing altogether for dominant quantum interactions. Quantum effects can be significant, leading to values that are notably lower than the classical thermodynamic limit, particularly for smaller drops. For planar interfaces (with infinite periodicity in the direction parallel to the interface), quantum effects are much less significant with the same values of λB but are, nevertheless, consequential for values representative of hydrogen or helium-4 at low temperatures corresponding to vapor-liquid coexistence. Large quantum effects are found for small drops of molecules with quantum interactions corresponding to water, ethane, methanol, and carbon dioxide, even at ambient conditions. The notable decrease in the density and tension has important consequences in reducing the Gibbs free-energy barrier of a nucleating cluster, enhancing the nucleation kinetics of liquid drops and of bubble formation. This implies that drops would form at a much greater rate than is predicted by classical nucleation theory.

Collective Proton Transfers in Cyclic Water-Ammonia Tetramers: A Path Integral Machine-Learning Study

We present results from machine-learning-based path integral molecular dynamics simulations that describe isomerization paths articulated via collective proton transfers along mixed, cyclic tetramers combining water and ammonia at cryogenic conditions. The net result of such isomerizations is a reverse of the chirality of the global hydrogen-bonding architecture along the different cyclic moieties. In monocomponent tetramers, the classical free energy profiles associated with these isomerizations present the usual symmetric double-well characteristics whereas the reactive paths exhibit full concertedness among the different intermolecular transfer processes. Contrastingly, in mixed water/ammonia tetramers, the incorporation of a second component introduces imbalances in the strengths of the different hydrogen bonds leading to a partial loss of concertedness, most notably at the vicinity of the transition state. As such, the highest and lowest degrees of progression are registered along OH···N and O···HN coordinations, respectively. These characteristics lead to polarized transition state scenarios akin to solvent-separated ion-pair configurations. The explicit incorporation of nuclear quantum effects promotes drastic depletions in the activation free energies and modifications in the overall shape of the profiles which include central plateau-like stages, indicating the prevalence of deep tunneling regimes. On the other hand, the quantum treatment of the nuclei partially restores the degree of concertedness among the evolutions of the individual transfers.

Coarse-Graining of Imaginary Time Feynman Path Integrals: Inclusion of Intramolecular Interactions and Bottom-up Force-Matching

Feynman's imaginary time path integral formalism of quantum statistical mechanics and the corresponding quantum-classical isomorphism provide a tangible way of incorporating nuclear quantum effect (NQE) in the simulation of condensed matter systems using well-developed classical simulation techniques. Our previous work has presented the many-body coarse-graining of path integral (CG-PI) theory that builds an isomorphism between the quantum partition function of N distinguishable particles and the classical partition function of 2N pseudoparticles. In this present work, we develop a generalized version of the many-body CG-PI theory that incorporates many-body interactions in the force field. Based on the new derivation, we provide a numerical CG-PI (n-CG-PI) modeling strategy parametrized from the underlying path integral molecular dynamics (PIMD) trajectories using force matching and Boltzmann inversion. The n-CG-PI models for two liquid systems are shown to capture well both the intramolecular and intermolecular structural correlations of the reference PIMD simulations. The generalized derivation of the many-body CG-PI theory and the n-CG-PI model presented in this work extend the scope of the CG-PI formalism by generalizing the previously limited theory to incorporate force fields of realistic molecular systems.

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

Evidence of a liquid-liquid phase transition in H[Formula: see text]O and D[Formula: see text]O from path-integral molecular dynamics simulations

We perform path-integral molecular dynamics (PIMD), ring-polymer MD (RPMD), and classical MD simulations of H[Formula: see text]O and D[Formula: see text]O using the q-TIP4P/F water model over a wide range of temperatures and pressures. The density [Formula: see text], isothermal compressibility [Formula: see text], and self-diffusion coefficients D(T) of H[Formula: see text]O and D[Formula: see text]O are in excellent agreement with available experimental data; the isobaric heat capacity [Formula: see text] obtained from PIMD and MD simulations agree qualitatively well with the experiments. Some of these thermodynamic properties exhibit anomalous maxima upon isobaric cooling, consistent with recent experiments and with the possibility that H[Formula: see text]O and D[Formula: see text]O exhibit a liquid-liquid critical point (LLCP) at low temperatures and positive pressures. The data from PIMD/MD for H[Formula: see text]O and D[Formula: see text]O can be fitted remarkably well using the Two-State-Equation-of-State (TSEOS). Using the TSEOS, we estimate that the LLCP for q-TIP4P/F H[Formula: see text]O, from PIMD simulations, is located at [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Isotope substitution effects are important; the LLCP location in q-TIP4P/F D[Formula: see text]O is estimated to be [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]. Interestingly, for the water model studied, differences in the LLCP location from PIMD and MD simulations suggest that nuclear quantum effects (i.e., atoms delocalization) play an important role in the thermodynamics of water around the LLCP (from the MD simulations of q-TIP4P/F water, [Formula: see text] MPa, [Formula: see text] K, and [Formula: see text] g/cm[Formula: see text]). Overall, our results strongly support the LLPT scenario to explain water anomalous behavior, independently of the fundamental differences between classical MD and PIMD techniques. The reported values of [Formula: see text] for D[Formula: see text]O and, particularly, H[Formula: see text]O suggest that improved water models are needed for the study of supercooled water.

Freezing point depression of salt aqueous solutions using the Madrid-2019 model

Salt aqueous solutions are relevant in many fields, ranging from biological systems to seawater. Thus, the availability of a force-field that is able to reproduce the thermodynamic and dynamic behavior of salt aqueous solutions would be of great interest. Unfortunately, this has been proven challenging, and most of the existing force-fields fail to reproduce much of their behavior. In particular, the diffusion of water or the salt solubility are often not well reproduced by most of the existing force-fields. Recently, the Madrid-2019 model was proposed, and it was shown that this force-field, which uses the TIP4P/2005 model for water and non-integer charges for the ions, provides a good description of a large number of properties, including the solution densities, viscosities, and the diffusion of water. In this work, we assess the performance of this force-field on the evaluation of the freezing point depression. Although the freezing point depression is a colligative property that at low salt concentrations depends solely on properties of pure water, a good model for the electrolytes is needed to accurately predict the freezing point depression at moderate and high salt concentrations. The coexistence line between ice and several salt aqueous solutions (NaCl, KCl, LiCl, MgCl₂, and Li₂SO₄) up to the eutectic point is estimated from direct coexistence molecular dynamics simulations. Our results show that this force-field reproduces fairly well the experimentally measured freezing point depression with respect to pure water freezing for all the salts and at all the compositions considered.

A tiny charge-scaling in the OPLS-AA + L-OPLS force field delivers the realistic dynamics and structure of liquid primary alcohols

We carry out molecular dynamics simulations for pure liquid primary alcohols ranging from methanol to 1-decanol under ambient conditions. Based on the OPLS-AA force field with the L-OPLS correction, we demonstrate that a few % increases in the partial charges deliver the realistic dynamics (self-diffusion coefficient and shear viscosity) and structure (density and X-ray scattering intensity) as well as enthalpy of vaporization and isothermal compressibility. The validity against thermal expansion coefficient, isobaric heat capacity, and static dielectric constant are also discussed.

Thermophysical properties of water using reactive force fields

The widescale importance and rich phenomenology of water continue to motivate the development of computational models. ReaxFF force fields incorporate many characteristics desirable for modeling aqueous systems: molecular flexibility, polarization, and chemical reactivity (bond formation and breaking). However, their ability to model the general properties of water has not been evaluated in detail. We present comprehensive benchmarks of the thermophysical properties of water for two ReaxFF models, the water-2017 and CHON-2017_weak force fields. These include structural, electrostatic, vibrational, thermodynamic, coexistence, and transport properties at ambient conditions (300 K and 0.997 g cm^-3) and along the standard pressure (1 bar) isobar. Overall, CHON-2017_weak predicts more accurate thermophysical properties than the water-2017 force field. Based on our results, we recommend potential avenues for improvement: the dipole moment to quadrupole moment ratio, the self-diffusion coefficient, especially for water-2017, and the gas phase vibrational frequencies with the aim to improve the vibrational properties of liquid water.

Classical nucleation theory of ice nucleation: Second-order corrections to thermodynamic parameters

Accurately estimating the nucleation rate is crucial in studying ice nucleation and ice-promoting and anti-freeze strategies. In classical nucleation theory, estimates of the ice nucleation rate are very sensitive to thermodynamic parameters, such as the chemical potential difference between water and ice Δμ and the ice-water interfacial free energy γ. However, even today, there are still many contradictions and approximations when estimating these thermodynamic parameters, introducing a large uncertainty in any estimate of the ice nucleation rate. Starting from basic concepts for a general solid-liquid crystallization system, we expand the Gibbs-Thomson equation to second order and derive second-order analytical formulas for Δμ, γ, and the nucleation barrier ΔG*, which are used in molecular dynamics simulations. These formulas describe well the temperature dependence of these thermodynamic parameters. This may be a new method of estimating Δμ, γ, and ΔG*.

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

Efficacy and Safety of a Nasopharyngeal Catheter for Selective Brain Cooling in Patients with Traumatic Brain Injury: A Prospective, Non-randomized Pilot Study

BACKGROUND

The efficacy objective was to determine whether a novel nasopharyngeal catheter could be used to cool the human brain after traumatic brain injury, and the safety objective was to assess the local and systemic effects of this therapeutic strategy.

METHODS

This was a prospective, non-randomized, interventional clinical trial that involved five patients with severe traumatic brain injury. The intervention consisted of inducing and maintaining selective brain cooling for 24 h by positioning a catheter in the nasopharynx and circulating cold water inside the catheter in a closed-loop arrangement. Core temperature was maintained at ≥ 35 °C using counter-warming.

RESULTS

In all study participants, a brain temperature reduction of ≥ 2 °C was achieved. The mean brain temperature reduction from baseline was 2.5 ± 0.9 °C (P = .04, 95% confidence interval). The mean systemic temperature was 37.3 ± 1.1 °C at baseline and 36.0 ± 0.8 °C during the intervention. The mean difference between the brain temperature and the systemic temperature during intervention was - 1.2 ± 0.8 °C (P = .04). The intervention was well tolerated with no significant changes observed in the hemodynamic parameters. No relevant variations in intracranial pressure and transcranial Doppler were observed. The laboratory results underwent no major changes, aside from the K+ levels and blood counts. The K+ levels significantly varied (P = .04); however, the variation was within the normal range. Only one patient experienced an event of mild localized and superficial nasal discoloration, which was re-evaluated on the seventh day and indicated complete recovery.

CONCLUSION

The results suggest that our noninvasive method for selective brain cooling, using a novel nasopharyngeal catheter, was effective and safe for use in humans.

Nuclear quantum effects on the thermodynamic, structural, and dynamical properties of water

We perform path-integral molecular dynamics (PIMD) simulations of H2O and D2O using the q-TIP4P/F model. Simulations are performed at P = 1 bar and over a wide range of temperatures that include the equilibrium (T≥ 273 K) and supercooled (210 ≤T < 273 K) liquid states of water. The densities of both H2O and D2O calculated from PIMD simulations are in excellent agreement with experiments in the equilibrium and supercooled regimes. We also evaluate important thermodynamic response functions, specifically, the thermal expansion coefficient αP(T), isothermal compressibility κT(T), isobaric heat capacity CP(T), and static dielectric constant ε(T). While these properties are in excellent [αP(T) and κT(T)] or semi-quantitative agreement [CP(T) and ε(T)] with experiments in the equilibrium regime, they are increasingly underestimated upon further cooling. It follows that the inclusion of nuclear quantum effects in PIMD simulations of (q-TIP4P/F) water is not sufficient to reproduce the anomalous large fluctuations in density, entropy, and electric dipole moment characteristic of supercooled water. It has been hypothesized that water may exhibit a liquid-liquid critical point (LLCP) in the supercooled regime at P > 1 bar and that such a LLCP generates a maximum in CP(T) and κT(T) at 1 bar. Consistent with this hypothesis and in particular, with experiments, we find a maximum in the κT(T) of q-TIP4P/F light and heavy water at T≈ 230-235 K. No maximum in CP(T) could be detected down to T≥ 210 K. We also calculate the diffusion coefficient D(T) of H2O and D2O using the ring-polymer molecular dynamics (RPMD) technique and find that computer simulations are in remarkable good agreement with experiments at all temperatures studied. The results from RPMD/PIMD simulations are also compared with the corresponding results obtained from classical MD simulations of q-TIP4P/F water where atoms are represented by single interacting sites. Surprisingly, we find minor differences in most of the properties studied, with CP(T), D(T), and structural properties being the only (expected) exceptions.

Bulk supercooled water versus adsorbed films on silica surfaces: specific heat by Monte Carlo simulation

Between 150 and 230.6 K, bulk supercooled water freezes upon cooling, and amorphous ice crystallizes upon heating: bulk water thus exists only in its stable ice form. To circumvent this problem, experiments are generally performed on water adsorbed in SiO2 based porous systems. In this work, we take advantage of Monte Carlo simulations to explore this metastable supercooled region inaccessible to experiments. Using three rigid, non-polarizable water models, namely SPC, TIP4P and TIP4P/2005, we investigate the isobaric specific heat capacity (Cp), between 100 and 300 K, of bulk water and water films of few monolayers adsorbed on different SiO2 surfaces: a smooth surface, a non-hydroxylated (0001) surface of quartz, and a fully hydroxylated (001) surface of cristobalite. As Cp is directly related to the entropy fluctuations and we focus on low temperatures, the convergence of the Monte Carlo simulations is a critical point of this work. Also, due to the small mass of the hydrogen atoms, quantum corrections are taken into account, and lead to an excellent agreement of the simulated and experimental Cp values at low temperature (100 K region). Altogether, we conclude that, in bulk, Cp is shown to exhibit a broad peak around 225 K for the SPC and TIP4P models, and around 250 K for the TIP4P/2005 model, in qualitative agreement with the experimentally observed features in Cp measurements. For interfacial water, in all cases, the broad Cp peak disappears. This result, at odds with experimental observations, suggests that disorder and hydrogen bonding at the interface (not yet taken into account) have a fundamental role in confined water transitions.

Real Space Triplets in Quantum Condensed Matter: Numerical Experiments Using Path Integrals, Closures, and Hard Spheres

Path integral Monte Carlo and closure computations are utilized to study real space triplet correlations in the quantum hard-sphere system. The conditions cover from the normal fluid phase to the solid phases face-centered cubic (FCC) and cI16 (de Broglie wavelengths , densities ). The focus is on the equilateral and isosceles features of the path-integral centroid and instantaneous structures. Complementary calculations of the associated pair structures are also carried out to strengthen structural identifications and facilitate closure evaluations. The three closures employed are Kirkwood superposition, Jackson-Feenberg convolution, and their average (AV3). A large quantity of new data are reported, and conclusions are drawn regarding (i) the remarkable performance of AV3 for the centroid and instantaneous correlations, (ii) the correspondences between the fluid and FCC salient features on the coexistence line, and (iii) the most conspicuous differences between FCC and cI16 at the pair and the triplet levels at moderately high densities (. This research is expected to provide low-temperature insights useful for the future related studies of properties of real systems (e.g., helium, alkali metals, and general colloidal systems).

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.

Accurate Description of Nuclear Quantum Effects with High-Order Perturbed Path Integrals (HOPPI)

Imaginary time path-integral (PI) simulations that account for nuclear quantum effects (NQE) beyond the harmonic approximation are increasingly employed together with modern electronic-structure calculations. Existing PI methods are applicable to molecules, liquids, and solids; however, the computational cost of such simulations increases dramatically with decreasing temperature. To address this challenge, here, we propose to combine high-order PI factorization with perturbation theory (PT). Already for conventional second-order PI simulations, the PT ansatz increases the accuracy 2-fold compared to fourth-order schemes with the same settings. In turn, applying PT to high-order path integrals (HOPI) further improves the efficiency of simulations for molecular and condensed matter systems especially at low temperatures. We present results for bulk liquid water, the aspirin molecule, and the CH₅⁺ molecule. Perturbed HOPI simulations remain both efficient and accurate down to 20 K and provide a convenient method to estimate the convergence of quantum-mechanical observables.

Implementation of Geometry-Dependent Charge Flux into the Polarizable AMOEBA+ Potential

Molecular dynamics (MD) simulations employing classical force fields (FFs) have been widely used to model molecular systems. The important ingredient of the current FFs, atomic charge, remains fixed during MD simulations despite the atomic environment or local geometry changes. This approximation hinders the transferability of the potential being used in multiple phases. Here we implement a geometry-dependent charge flux (GDCF) model into the multipole-based AMOEBA+ polarizable potential. The CF in the current work explicitly depends on the local geometry (bond and angle) of the molecule. To our knowledge, this is the first study that derives energy and force expressions due to GDCF in a multipole-based polarizable FF framework. Due to the inclusion of GDCF, the AMOEBA+ water model is noticeably improved in terms of describing the monomer properties, cluster binding/interaction energy, and a variety of liquid properties, including the infrared spectra that previous flexible water models were not able to capture.

Coarse-graining of many-body path integrals: Theory and numerical approximations

Feynman's imaginary time path integral approach to quantum statistical mechanics provides a theoretical formalism for including nuclear quantum effects (NQEs) in simulation of condensed matter systems. Sinitskiy and Voth [J. Chem. Phys. 143, 094104 (2015)] have presented the coarse-grained path integral (CG-PI) theory, which provides a reductionist coarse-grained representation of the imaginary time path integral based on the quantum-classical isomorphism. In this paper, the many-body generalization of the CG-PI theory is presented. It is shown that the N interacting particles obeying quantum Boltzmann statistics can be represented as a system of N pairs of classical-like pseudoparticles coupled to each other analogous to the pseudoparticle pair of the one-body theory. Moreover, we present a numerical CG-PI (n-CG-PI) method applying a simple approximation to the coupling scheme between the pseudoparticles due to numerical challenges of directly implementing the full many-body CG-PI theory. Structural correlations of two liquid systems are investigated to demonstrate the performance of the n-CG-PI method. Both the many-body CG-PI theory and the n-CG-PI method not only present reductionist views of the many-body quantum Boltzmann statistics but also provide theoretical and numerical insight into how to explicitly incorporate NQEs in the representation of condensed matter systems with minimal additional degrees of freedom.

Water-retaining, tough and self-healing hydrogels and their uses as fire-resistant materials

Hydrogels as fire-resistant materials have attracted great attention due to their high water content and tailored shapes that can cover various surfaces.

A new intermolecular potential for simulations of methanol: The OPLS/2016 model

In this work, a new rigid-nonpolarizable model of methanol is proposed. The model has three sites, located at the same positions as those used in the OPLS model previously proposed by Jorgensen [J. Phys. Chem. 90, 1276 (1986)]. However, partial charges and the values of the Lennard-Jones parameters were modified by fitting to an adequately selected set of target properties including solid-fluid experimental data. The new model was denoted as OPLS/2016. The overall performance of this model was evaluated and compared to that obtained with other popular models of methanol using a similar test to that recently proposed for water models. In the test, a certain numerical score is given to each model. It was found that the OPLS/2016 obtained the highest score (7.4 of a maximum of 10) followed by L1 (6.6), L2 (6.4), OPLS (5.8), and H1 (3.5) models. The improvement of OPLS/2016 with respect to L1 and L2 is mainly due to an improvement in the description of fluid-solid equilibria (the melting point is only 14 K higher than the experimental value). In addition, it was found that no methanol model was able to reproduce the static dielectric constant and the isobaric heat capacity, whereas the better global performance was found for models that reproduce the vaporization enthalpy once the so-called polarization term is included. Similar conclusions were suggested previously in the analysis of water models and are confirmed here for methanol.

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.

Dielectric Anomaly in Ice near 20 K: Evidence of Macroscopic Quantum Phenomena

H2O is one of the most important substances needed in sustaining life; but not much is known about its ground state. Here a previously unidentified anomaly is identified in the form of a minimum in the imaginary part of the dielectric constant with respect to temperature near 20 K, while the real part remains monotonic. Isothermal dispersion and absorption measurements show coinciding results. For the case of heavy ice (D2O), no anomaly was identified, confirming an apparent isotope effect. Concerted quantum tunneling of protons is believed to be the main cause behind the reported anomaly. Our findings identify another system that exhibits macroscopic quantum phenomena that rarely occur in nature.

Thermophysical properties of supercritical water and bond flexibility

Molecular dynamics results are reported for the thermodynamic properties of supercritical water using examples of both rigid (TIP4P/2005) and flexible (TIP4P/2005f) transferable interaction potentials. Data are reported for pressure, isochoric and isobaric heat capacities, the thermal expansion coefficient, isothermal and adiabatic compressibilities, Joule-Thomson coefficient, speed of sound, self-diffusion coefficient, viscosities, and thermal conductivity. Many of these properties have unusual behavior in the supercritical phase such as maximum and minimum values. The effectiveness of bond flexibility on predicting these properties is determined by comparing the results to experimental data. The influence of the intermolecular potential on these properties is both variable and state point dependent. In the vicinity of the critical density, the rigid and flexible potentials yield very different values for the compressibilities, heat capacities, and thermal expansion coefficient, whereas the self-diffusion coefficient, viscosities, and thermal conductivities are much less potential dependent. Although the introduction of bond flexibility is a computationally expedient way to improve the accuracy of an intermolecular potential, it can be counterproductive in some cases and it is not an adequate replacement for incorporating the effects of polarization.

Understanding the rates and molecular mechanism of water-exchange around aqueous ions using molecular simulations

Solvation processes occurring around aqueous ions are of fundamental importance in physics, chemistry, and biology. Over the past few decades, several experimental and theoretical studies were devoted to understanding ion solvation and the processes involved in it. In this article, we present a summary of our recent efforts that, through computer simulations, focused on providing a comprehensive understanding of solvent-exchange processes around aqueous ions. To accomplish these activities, we have looked at the mechanistic properties associated with the water-exchange process, such as potentials of mean force, time-dependent transmission coefficients, and the corresponding rate constants using transition state theory, the reactive flux method, and Grote-Hynes treatments of the dynamic response of the solvent.

Path-integral simulation of solids

The path-integral formulation of the statistical mechanics of quantum many-body systems is described, with the purpose of introducing practical techniques for the simulation of solids. Monte Carlo and molecular dynamics methods for distinguishable quantum particles are presented, with particular attention to the isothermal-isobaric ensemble. Applications of these computational techniques to different types of solids are reviewed, including noble-gas solids (helium and heavier elements), group-IV materials (diamond and elemental semiconductors), and molecular solids (with emphasis on hydrogen and ice). Structural, vibrational, and thermodynamic properties of these materials are discussed. Applications also include point defects in solids (structure and diffusion), as well as nuclear quantum effects in solid surfaces and adsorbates. Different phenomena are discussed, as solid-to-solid and orientational phase transitions, rates of quantum processes, classical-to-quantum crossover, and various finite-temperature anharmonic effects (thermal expansion, isotopic effects, electron-phonon interactions). Nuclear quantum effects are most remarkable in the presence of light atoms, so that especial emphasis is laid on solids containing hydrogen as a constituent element or as an impurity.

A systematic development of a polarizable potential of water

Based on extensive studies of existing potentials we propose a new molecular model for water. The new model is rigid and contains three Gaussian charges. Contrary to other models, all charges take part in the polarization of the molecule. They are connected by harmonic springs to their gas-phase positions: the negative charge to a prescribed point on the main axis of the molecule; the positive charges to the hydrogens. The mechanical equilibrium between the electrostatic forces and the spring forces determines the polarization of the molecule which is established by iteration at every timestep. The model gives excellent estimates for ambient liquid properties and reasonably good results from high-pressure solids to gas-phase clusters. We present a detailed description of the development of this model and a large number of calculated properties compared to the estimates of the nonpolarizable TIP4P∕2005 [J. L. F. Abascal and C. Vega, J. Chem. Phys. 123, 234505 (2005)], the polarizable GCPM [P. Paricaud, M. Predota, A. A. Chialvo, and P. T. Cummings, J. Chem. Phys. 122, 244511 (2005)], and our earlier BKd3 model [P. T. Kiss and A. Baranyai, J. Chem. Phys. 137, 084506 (2012)]. The best overall performance is shown by the new model.

Quantum mechanical force field for water with explicit electronic polarization

A quantum mechanical force field (QMFF) for water is described. Unlike traditional approaches that use quantum mechanical results and experimental data to parameterize empirical potential energy functions, the present QMFF uses a quantum mechanical framework to represent intramolecular and intermolecular interactions in an entire condensed-phase system. In particular, the internal energy terms used in molecular mechanics are replaced by a quantum mechanical formalism that naturally includes electronic polarization due to intermolecular interactions and its effects on the force constants of the intramolecular force field. As a quantum mechanical force field, both intermolecular interactions and the Hamiltonian describing the individual molecular fragments can be parameterized to strive for accuracy and computational efficiency. In this work, we introduce a polarizable molecular orbital model Hamiltonian for water and for oxygen- and hydrogen-containing compounds, whereas the electrostatic potential responsible for intermolecular interactions in the liquid and in solution is modeled by a three-point charge representation that realistically reproduces the total molecular dipole moment and the local hybridization contributions. The present QMFF for water, which is called the XP3P (explicit polarization with three-point-charge potential) model, is suitable for modeling both gas-phase clusters and liquid water. The paper demonstrates the performance of the XP3P model for water and proton clusters and the properties of the pure liquid from about 900 × 10(6) self-consistent-field calculations on a periodic system consisting of 267 water molecules. The unusual dipole derivative behavior of water, which is incorrectly modeled in molecular mechanics, is naturally reproduced as a result of an electronic structural treatment of chemical bonding by XP3P. We anticipate that the XP3P model will be useful for studying proton transport in solution and solid phases as well as across biological ion channels through membranes.