Intrusion of fluids into nanogrooves: how geometry determines the shape of the gas-liquid interface

We study the shape of gas-liquid interfaces forming inside rectangular nanogrooves (i.e., slit-pores capped on one end). On account of purely repulsive fluid-substrate interactions the confining walls are dry (i.e., wet by vapor) and a liquid-vapor interface intrudes into the nanogrooves to a distance determined by the pressure (i.e., chemical potential). By means of Monte Carlo simulations in the grand-canonical ensemble (GCEMC) we obtain the density rho(z) along the midline (x = 0) of the nanogroove for various geometries (i.e., depths D and widths L) of the nanogroove. We analyze the density profiles with the aid of an analytic expression which we obtain through a transfer-matrix treatment of a one-dimensional effective interface Hamiltonian. Besides geometrical parameters such as D and L , the resulting analytic expression depends on temperature T , densities of coexisting gas and liquid phases in the bulk rho g,l(x) and the interfacial tension gamma. The latter three quantities are determined in independent molecular dynamics simulations of planar gas-liquid interfaces. Our results indicate that the analytic formula provides an excellent representation of rho(z) as long as L is sufficiently small. At larger L the meniscus of the intruding liquid flattens. Under these conditions the transfer-matrix analysis is no longer adequate and the agreement between GCEMC data and the analytic treatment is less satisfactory.