This dataset contains input files, training data and other files related to the machine learning models developed during the work by Muniz et al. In this work, we construct machine learning models based on the MB-pol many-body model. We find that the training set should include cluster configurations as well as liquid phase configurations in order to accurately represent both liquid and VLE properties. The results attest for the ability of machine learning models to accurately represent many-body potentials and provide an efficient avenue for water simulations.
Piaggi, Pablo M; Gartner, Thomas E; Car, Roberto; Debenedetti, Pablo G
The possible existence of a liquid-liquid critical point in deeply supercooled water has been a subject of debate in part due to the challenges associated with providing definitive experimental evidence. Pioneering work by Mishima and Stanley [Nature 392, 164 (1998) and Phys.~Rev.~Lett. 85, 334 (2000)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected location of the liquid-liquid transition and its associated critical point. Based on the continuous or discontinuous changes in slope of the melting curves, Mishima suggested that the liquid-liquid critical point lies between the melting curves of ice III and ice V. Here, we explore this conjecture using molecular dynamics simulations with a purely-predictive machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII using this model and find that the melting lines of all the studied ice polymorphs are supercritical and do not intersect the liquid-liquid transition locus. We also find a pronounced, yet continuous, change in slope of the melting lines upon crossing of the locus of maximum compressibility of the liquid. Finally, we analyze critically the literature in light of our findings, and conclude that the scenario in which melting curves are supercritical is favored by the most recent computational and experimental evidence. Thus, although the preponderance of experimental and computational evidence is consistent with the existence of a second critical point in water, the behavior of the melting lines of ice polymorphs does not provide strong evidence in support of this viewpoint, according to our calculations.
This dataset contains input and output files to reproduce the results of the manuscript "Homogeneous ice nucleation in an ab initio machine learning model" by Pablo M. Piaggi, Jack Weis, Athanassios Z. Panagiotopoulos, Pablo G. Debenedetti, and Roberto Car (arXiv preprint https://arxiv.org/abs/2203.01376). In this work, we studied the homogeneous nucleation of ice from supercooled liquid water using a machine learning model trained on ab initio energies and forces. Since nucleation takes place over times much longer than the simulation times that can be afforded using molecular dynamics simulations, we make use of the seeding technique that is based on simulating an ice cluster embedded in liquid water. 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 shrinking at the given supersaturation). Using data from the seeding simulations and the equations of classical nucleation theory we compute nucleation rates that can be compared with experiments.
These GROMACS trajectories show the existence of a critical point in deeply supercooled WAIL water. Also included is the code necessary to reproduce the figures in the corresponding paper from these trajectories. From this data the critical temperature, pressure, and density of the model can be found, and critical fluctuations in the deeply supercooled liquid can be directly observed (in a computer-simulation sense).
Muniz, Maria Carolina; Gartner III, Thomas E.; Riera, Marc; Knight, Christopher; Yue, Shuwen; Paesani, Francesco; Panagiotopoulos, Athanassios Z.
This dataset contains all data (including input files, simulation trajectories as well as other data files and analysis scripts) related to the publication "Vapor-liquid equilibrium of water with the MB-pol many-body potential" by Muniz et al. in preparation (2021). In this work, we assessed the performance of the MB-pol many-body potential with respect to water's vapor-liquid equilibrium properties. Through the use of direct coexistence molecular dynamics, we calculated properties such as coexistence densities, surface tension, vapor pressures and enthalpy of vaporization. We found that MB-pol is able to predict these properties in good agreement with experimental data. The results attest to the chemical accuracy of MB-pol and its large range of application across water's phase diagram.
Gartner, Thomas III; Zhang, Linfeng; Piaggi, Pablo; Car, Roberto; Panagiotopoulos, Athanassios; Debenedetti, Pablo
This dataset contains all data related to the publication "Signatures of a liquid-liquid transition in an ab initio deep neural network model for water", by Gartner et al., 2020. In this work, we used neural networks to generate a computational model for water using high-accuracy quantum chemistry calculations. Then, we used advanced molecular simulations to demonstrate evidence that suggests this model exhibits a liquid-liquid transition, a phenomenon that can explain many of water's anomalous properties. This dataset contains links to all software used, all data generated as part of this work, as well as scripts to generate and analyze all data and generate the plots reported in the publication.