Monte Carlo simulation methods for computing the wetting and drying properties of model systems

Abstract

We introduce general Monte Carlo simulation methods for determining the wetting and drying properties of model systems. We employ an interface-potential-based approach in which the interfacial properties of a system are related to the surface excess free energy of a thin fluid film in contact with a surface. Two versions of this approach are explored: a "spreading" method focused on the growth of a thin liquid film from a surface in a mother vapor and a "drying" method focused on the growth of a thin vapor film from a surface in a mother liquid. The former provides a direct measure of the spreading coefficient while the latter provides an analogous drying coefficient. When coupled with an independent measure of the liquid-vapor surface tension, these coefficients enable one to compute the contact angle. We also show how one can combine information gathered from application of the spreading and drying methods at a common state point to obtain direct measures of the contact angle and liquid-vapor surface tension. The computational strategies introduced here are applied to two model systems. One includes a monatomic Lennard-Jones fluid that interacts with a structureless substrate via a long-ranged substrate potential. The second model contains a monatomic Lennard-Jones fluid that interacts with an atomistically detailed substrate via a short-ranged potential. Expanded ensemble techniques are coupled with the interface potential approach to compile the temperature- and substrate strength-dependence of various interfacial properties for these systems. Overall, we find that the approach pursued here provides an efficient and precise means to calculate the wetting and drying properties of model systems.