New Computational Approach to Determine Liquid-Solid Phase Equilibria of Water Confined to Slit Nanopores

We devise a new computational approach to compute solid-liquid phase equilibria of confined fluids. Specifically, we extend the multibaric-multithermal ensemble method with an anisotropic pressure control to achieve the solid-liquid phase equilibrium for confined water inside slit nanopores (with slit width h ranging from 5.4 angstrom to 7.2 angstrom). A unique feature of this multibaric-multithermal ensemble is that the freezing points of confined water can be determined from the heat-capacity peaks. The new approach has been applied to compute the freezing point of two monolayer ices, namely, a high-density flat rhombic monolayer ice (HD-fRMI) and a high-density puckered rhombic monolayer ice (HD-pRMI) observed in our simulation. We find that the liquid-to-solid transition temperature (or the freezing point) of HD-pRMI is dependent on the slit width h, whereas that of HD-fRMI is nearly independent of the h.