Effect of confinement on water properties in super-hydrophilic pores using MD simulations with the mW model

Context

We explore the influence of strongly hydrophilic confinement on various properties of water, such as density, enthalpy, potential energy, radial distribution function, entropy, specific heat capacity, structural dynamics, and transition temperatures (freezing and melting temperatures), using monatomic water (mW) model. The properties of water are found to be dependent on confinement and the wall-fluid surface interaction. Hysteresis loops are observed for density, enthalpy, potential energy, and entropy around the transition temperatures, while the size of hysteresis loops varies with confinement and surface interaction. In smaller pore sizes (H ≤ 20), the solid phase displays a higher density compared to the liquid phase, which is unconventional behavior compared to bulk water systems due to the pronounced hydrophilic properties of the confinement surface. Specific heat capacity exhibits more oscillations in the confined system compared to bulk water, stemming from uneven enthalpy differences across equal temperature intervals. During phase transformation in both heating and quenching processes, there is an abrupt change observed in specific heat capacity. Confinement exerts a notable impact on entropy in the solid phase, but its influence is negligible in the liquid phase. At lower pore sizes (H < 25 Å), there is more fluctuation in freezing temperature for all wall-fluid interactions, which diminishes beyond pore sizes of H > 25 Å. Similarly, more oscillatory behavior is observed in melting temperatures at lower pore sizes (H < 40 Å), which diminishes at higher pore sizes (H > 40 Å). During the quenching process, a sudden jump in the in-plane orientational and tetrahedral order parameters indicates the formation of an ordered phase, specifically a diamond crystalline structure. The percentages of different crystalline structures (cubic diamond, hexagonal diamond, and 2D hexagonal) vary with both the confinement size and the wall-fluid interaction strength.

Methods

Cooling and heating simulations are conducted with the mW water model using LAMMPS for different nanoscale confinement separation sizes ranging from 8.5 to 70 Å within the temperature range of 100–350 K. The water is modeled using two-body and three-body interaction potential (Stillinger–Weber potential) and the confinement is introduced using LJ 9–3 water-wall interaction potential. Entropy is calculated using RDF data obtained from the simulation experiments for each temperature point with increments or decrements of 2.5 K. The transition temperatures are estimated using the specific heat capacity analysis.