Understanding wetting of immiscible liquids near a solid surface using molecular simulation

We introduce Monte Carlo simulation methods for determining interfacial properties of fluids that exhibit bulk liquid-liquid immiscibility. 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, is utilized to deduce the wetting characteristics of these systems. We present a framework for implementing this general method within both the grand canonical and semigrand isobaric-isothermal ensembles. Tracking the evolution of interfacial properties along various thermodynamic paths is also examined. This task is accomplished by implementing variants of the expanded ensemble technique, which enables one to obtain components of the interface potential along a path of interest. We also discuss how these concepts are employed to calculate bulk liquid-liquid coexistence properties in an efficient manner. The computational strategies introduced here are applied to three model Lennard-Jones systems. For each system, we compile the evolution of the liquid-liquid surface tension and contact angle with temperature or pressure. For one of the model systems we compare our results with literature data. We also examine how interfacial properties evolve upon variation of the relative affinity of the fluid components for the substrate. Overall, we find that the approach pursued here is generally applicable and provides an efficient and precise means to calculate the bulk and interfacial properties of fluids that exhibit liquid-liquid immiscibility.