# Simulations of water nano-confined between corrugated planes

## TL;DR

This study uses computer simulations to explore how water confined between corrugated planes responds structurally and dynamically, revealing phase changes and inhomogeneous behavior influenced by wall oscillations.

## Contribution

It introduces a novel periodic confining potential to simulate atomistic wall oscillations, analyzing their effects on water's phase behavior with both empirical and first-principles methods.

## Key findings

- Melting temperature increases with commensurate lattice parameters.
- Hexatic phase is replaced by trilayer crystalline phase at high modulation amplitudes.
- Structural and dynamical responses are surprisingly insensitive to wall rugosity.

## Abstract

Two-dimensionally nanoconfined water 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 under a periodic confining potential by means of computer simulations, both using empirical potentials and from first-principles. 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, an overall rise of the melting temperature is observed. 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, both agree remarkably well for the main conclusions of the study.