Publication date: May 03, 2019
We present a method for performing multithermal-multibaric molecular dynamics simulations that sample entire regions of the temperature-pressure (TP) phase diagram. The method uses a variational principle [Valsson and Parrinello, Phys. Rev. Lett. 113, 090601 (2014)] in order to construct a bias that leads to a uniform sampling in energy and volume. The intervals of temperature and pressure are taken as inputs and the relevant energy and volume regions are determined on the fly. In this way the method guarantees adequate statistics for the chosen TP region. We show that our multithermal-multibaric simulations can be used to calculate all static physical quantities for all temperatures and pressures in the targeted region of the TP plane. We illustrate our approach by studying the density anomaly of TIP4P/Ice water. This record includes input and output files, and Jupyter Notebooks describing the analysis of the simulations and the creation of the figures for the paper.
No Explore or Discover sections associated with this archive record.
|15.6 KiB||Source files that should be added to Plumed v2.4 in order to perform and analyze the simulations.|
|9.8 MiB||A simulation of liquid Na in the multicanonical ensemble. References computed in the canonical ensemble at different temperatures. A Jupyter Notebook with the code used to generate the results shown in the paper.|
|3.1 GiB||A simulation of liquid TIP4P/Ice water in the multithermal-multibaric ensemble. References computed in the isothermal-isobaric ensemble. A Jupyter Notebook with the code used to generate the results shown in the paper.|
|1.4 KiB||README file with description of the other files|