Simulations of water nano-confined between corrugated planes

Water confined to nanoscale widths in two dimensions between ideal planar walls has been the subject of ample study, aiming at understanding the intrinsic response of water to confinement, avoiding the consideration of the chemistry of actual confining materials. In this work, we study the response of such nanoconfined water to the imposition of a periodicity in the confinement by means of computer simulations, both using empirical potentials and from first-principles. For that we propose a periodic confining potential emulating the atomistic oscillation of the confining walls, which allows varying the lattice parameter and amplitude of the oscillation. We do it for a triangular lattice, with several values of the lattice parameter: one which is ideal for commensuration with layers of Ih ice and other values that would correspond to more realistic substrates. For the former, the phase diagram shows an overall rise of the melting temperature. The liquid maintains a bi-layer triangular structure, however, despite the fact that it is not favoured by the external periodicity. The first-principles liquid is significantly affected by the modulation in its layering and stacking even at relatively small amplitudes of the confinement modulation. Beyond some critical modulation amplitude, the hexatic phase present in flat confinement is replaced by a trilayer crystalline phase unlike any of the phases encountered for flat confinement. For more realistic lattice parameters, the liquid does not display higher tendency to freeze, but it clearly shows inhomogeneous behaviour as the strength of the rugosity increases. In spite of this expected inhomogeneity, the structural and dynamical response of the liquid is surprisingly insensitive to the external modulation. Although the first-principles calculations give a more triangular liquid than the one observed with empirical potentials (TIP4P/2005), both agree remarkably well for the main conclusions of the study.