Experimental and theoretical bulk phase diagram and interfacial tension of ouzo.

Ouzo is a well-known drink in Mediterranean countries, with ingredients water, alcohol and trans-anethole oil. The oil is insoluble in water, but completely soluble in alcohol, so when water is added to the spirit, the available alcohol is depleted and the mixture exhibits spontaneous emulsification. This process is commonly known as the louche or ouzo effect. Although the phase boundaries of this archetypal ternary mixture are well known, the properties of coexisting phases have not previously been studied. Here, we present a detailed experimental investigation into the phase behaviour, including tie-lines connecting coexisting phases, determination of the critical point (also called the plait point in ternary systems) and measurements of the surface tension and density for varying alcohol concentrations. Additionally, we present a theory for the thermodynamics and phase diagram of the system. With suitable selection of the interaction parameters, the theory captures nearly all features of the experimental work. This simple model can be used to determine both bulk and non-uniform (e.g. interfacial) properties, paving the way for a wide range of future applications of the model to ternary mixtures in general. We show how our accurate equilibrium phase diagram can be used to provide improved understanding of non-equilibrium phenomena.

How ice grows from premelting films and water droplets.

Close to the triple point, the surface of ice is covered by a thin liquid layer (so-called quasi-liquid layer) which crucially impacts growth and melting rates. Experimental probes cannot observe the growth processes below this layer, and classical models of growth by vapor deposition do not account for the formation of premelting films. Here, we develop a mesoscopic model of liquid-film mediated ice growth, and identify the various resulting growth regimes. At low saturation, freezing proceeds by terrace spreading, but the motion of the buried solid is conveyed through the liquid to the outer liquid-vapor interface. At higher saturations water droplets condense, a large crater forms below, and freezing proceeds undetectably beneath the droplet. Our approach is a general framework that naturally models freezing close to three phase coexistence and provides a first principle theory of ice growth and melting which may prove useful in the geosciences.

Rounded Layering Transitions on the Surface of Ice.

Understanding the wetting properties of premelting films requires knowledge of the film's equation of state, which is not usually available. Here we calculate the disjoining pressure curve of premelting films and perform a detailed thermodynamic characterization of premelting behavior on ice. Analysis of the density profiles reveals the signature of weak layering phenomena, from one to two and from two to three water molecular layers. However, disjoining pressure curves, which closely follow expectations from a renormalized mean field liquid state theory, show that there are no layering phase transitions in the thermodynamic sense along the sublimation line. Instead, we find that transitions at mean field level are rounded due to capillary wave fluctuations. We see signatures that true first order layering transitions could arise at low temperatures, for pressures between the metastable line of water-vapor coexistence and the sublimation line. The extrapolation of the disjoining pressure curve above water-vapor saturation displays a true first order phase transition from a thin to a thick film consistent with experimental observations.

Binding potentials for vapour nanobubbles on surfaces using density functional theory.

We calculate density profiles of a simple model fluid in contact with a planar surface using density functional theory (DFT), in particular for the case where there is a vapour layer intruding between the wall and the bulk liquid. We apply the method of Hughes et al (2015 J. Chem. Phys. 142 074702) to calculate the density profiles for varying (specified) amounts of the vapour adsorbed at the wall. This is equivalent to varying the thickness h of the vapour at the surface. From the resulting sequence of density profiles we calculate the thermodynamic grand potential as h is varied and thereby determine the binding potential as a function of h. The binding potential obtained via this coarse-graining approach allows us to determine the disjoining pressure in the film and also to predict the shape of vapour nano-bubbles on the surface. Our microscopic DFT based approach captures information from length scales much smaller than some commonly used models in continuum mechanics.

The pressure tensor across a liquid-vapour interface.

Inhomogeneous fluids exhibit physical properties that are neither uniform nor isotropic. The pressure tensor is a case in point, key to the mechanical description of the interfacial region. Kirkwood and Buff and, later, Irving and Kirkwood, obtained a formal treatment based on the analysis of the pressure across a planar surface [J. G. Kirkwood and F. P. Buff, J. Chem. Phys. 17(3), 338 (1949); J. H. Irving and J. G. Kirkwood, J. Chem. Phys. 18, 817 (1950)]. We propose a generalisation of Irving and Kirkwood's argument to fluctuating, non-planar surfaces and obtain an expression for the pressure tensor that is not smeared by thermal fluctuations at the molecular scale and corresponding capillary waves [F. P. Buff et al., Phys. Rev. Lett. 15, 621-623 (1965)]. We observe the emergence of surface tension, defined as an excess tangential stress, acting exactly across the dividing surface at the sharpest molecular resolution. The new statistical mechanical expressions extend current treatments to fluctuating inhomogeneous systems far from equilibrium.

Films, layers, and droplets: The effect of near-wall fluid structure on spreading dynamics.

We present a study of the spreading of liquid droplets on a solid substrate at very small scales. We focus on the regime where effective wetting energy (binding potential) and surface tension effects significantly influence steady and spreading droplets. In particular, we focus on strong packing and layering effects in the liquid near the substrate due to underlying density oscillations in the fluid caused by attractive substrate-liquid interactions. We show that such phenomena can be described by a thin-film (or long-wave or lubrication) model including an oscillatory Derjaguin (or disjoining or conjoining) pressure and explore the effects it has on steady droplet shapes and the spreading dynamics of droplets on both an adsorption (or precursor) layer and completely dry substrates. At the molecular scale, commonly used two-term binding potentials with a single preferred minimum controlling the adsorption layer height are inadequate to capture the rich behavior caused by the near-wall layered molecular packing. The adsorption layer is often submonolayer in thickness, i.e., the dynamics along the layer consists of single-particle hopping, leading to a diffusive dynamics, rather than the collective hydrodynamic motion implicit in standard thin-film models. We therefore modify the model in such a way that for thicker films the standard hydrodynamic theory is realized, but for very thin layers a diffusion equation is recovered.

On the moving contact line singularity: asymptotics of a diffuse-interface model.

The behaviour of a solid-liquid-gas system near the three-phase contact line is considered using a diffuse-interface model with no-slip at the solid and where the fluid phase is specified by a continuous density field. Relaxation of the classical approach of a sharp liquid-gas interface and careful examination of the asymptotic behaviour as the contact line is approached is shown to resolve the stress and pressure singularities associated with the moving contact line problem. Various features of the model are scrutinised, alongside extensions to incorporate slip, finite-time relaxation of the chemical potential, or a precursor film at the wall.