Phase separation of binary fluids in porous media: Asymmetries in pore geometry and fluid composition

Surface-directed spinodal decomposition of fluids confined in a cylindrical pore

The surface-directed spinodal decomposition of a binary liquid confined inside a cylindrical pore is investigated using molecular dynamics simulations. One component of the liquid wets the pore surface while the other remains neutral. A variety of wetting conditions are studied. For the partial wetting case, after an initial period of phase separation, the domains organize themselves into pluglike structures and the system enters into a metastable state. Therefore, a complete phase separation is never achieved. Analysis of domain growth and the structure factor suggests a one-dimensional growth dynamics for the partial wetting case. As the wetting interaction is increased beyond a critical value, a transition from the pluglike to tubelike domain formation is observed, which corresponds to the full wetting morphology. Thus, a complete phase separation is achieved as the wetting species moves towards the pore surface and forms layers enclosing the nonwetting species residing around the axis of the cylinder. The coarsening dynamics of both the species are studied separately. The wetting species is found to follow a two-dimensional domain growth dynamics with a growth exponent 1/2 in the viscous hydrodynamic regime. This was substantiated by the Porod tail of the structure factor. On the other hand, the domain grows linearly with time for the nonwetting species. This suggests that the nonwetting species behaves akin to a three-dimensional bulk system. An appropriate reasoning is presented to justify the given observations.

Molecular Dynamics Simulation of Tolman Length and Interfacial Tension of Symmetric Binary Lennard–Jones Liquid

The Tolman length and interfacial tension of partially miscible symmetric binary Lennard–Jones (LJ) fluids (A, B) was revealed by performing a large-scale molecular dynamics (MD) simulation with a sufficient interfacial area and cutting distance. A unique phenomenon was observed in symmetric binary LJ fluids, where two surfaces of tension existed on both sides of an equimolar dividing surface. The range of interaction εAB between the different liquids and the temperature in which the two LJ fluids partially mixed was clarified, and the Tolman length exceeded 3 σ when εAB was strong at higher temperatures. The results show that as the temperature or εAB increases, the Tolman length increases and the interfacial tension decreases. This very long Tolman length indicates that one should be very careful when applying the concept of the liquid–liquid interface in the usual continuum approximation to nanoscale droplets and capillary phase separation in nanopores.

Phase diagrams for model heteronuclear dimers confined in narrow slitlike pores

Monte Carlo simulation is used to study the phase behavior of a model of heteronuclear dimers on a lattice. The calculations are carried out for two- and three-dimensional bulk systems and for dimers confined in slitlike pores of different widths and for different segment-wall interactions. The calculations recover a rich variety of phase behaviors that includes transitions between gaslike and liquidlike phases, the latter being ordered or disordered. The topology of the obtained diagrams depends on the energy of interactions between segments forming a dimer, on the segment-wall interaction, and on the pore width.

Phase separation of the Widom–Rowlinson mixture in confined geometry

The behavior of binary Widom-Rowlinson [J. Chem. Phys. 52, 1670 (1970)] mixture in confined geometry is investigated. The influence of confinement on phase separation is examined. Coexistence curves for the mixture in slitlike pores and pores of complex geometry are calculated in Monte Carlo simulations. Finite size scaling analysis is used to determine precisely the location of critical point for the bulk mixture.

Phase transitions and quantum effects in pore condensates: a path integral Monte Carlo study

Lennard-Jones condensates in cylindrical pores are studied by path integral Monte Carlo simulations with particular emphasis on phase transitions and quantum effects. The pore diameter effect and the influence of the interaction strength between the cylinder wall and the adsorbate particles on the structures and the location of the phase boundaries is studied and the quantum effect on the phase diagram is quantified by path integral Monte Carlo simulations. In case of strong wall-particle interactions good qualitative agreement with recent experimental results for the freezing of Ar-pore condensates is found. Meniscus structures in the solid phase are obtained as well as unexpected condensate structures for the system with the lighter Ne-particles due to quantum delocalization effects. The quantum effect on the freezing temperature can be as large as 10% in these systems.

Phase behavior of confined symmetric binary mixtures

We employ mean-field lattice density functional theory to investigate the phase behavior of a binary (A-B) mixture confined to nanoscopic slit pores with chemically homogeneous walls. We consider only nearest-neighbor interactions in symmetric mixtures, where epsilon(AA)=epsilon(BB) not equal epsilon(AB) and epsilon is a measure of attraction between molecules of like (subscripts AA and BB) and unlike species (subscript AB), respectively. In addition, molecules are exposed to short-range attraction by the substrates separated by z lattice planes where epsilon(W) is the relevant coupling parameter. Moreover, the chemical potentials of both components are the same, that is, mu(A)=mu(B)=mu. In thermodynamic equilibrium (for fixed temperature T and chemical potential mu) the grand-potential density omega[rho,m] (rho identical with [rho(1),...,rho(z)], m identical with [m(1),...,m(z)]) assumes a global minimum which we find by minimizing omega numerically with respect to the order parameters rho(l) identical with rho(A)(l)+rho(B)(l) (total local density) and m(l) identical with (rho(A)(l)-rho(B)(l))/rho(l) (local "miscibility") at lattice plane l parallel to the pore walls. By varying epsilon(AB) three generic types of bulk phase diagrams are observed. On account of confinement (i.e., by varying epsilon(W) as well as z) one may switch between these different types of phase diagrams. This may have profound practical repercussions for experimental nanophase separation since depending on pore width and chemical nature of its walls a bulk gas mixture may undergo capillary condensation and form either a stable mixed or demixed liquid phase.

Low-temperature phase separation of a binary liquid mixture in porous materials studied by cryoporometry and pulsed-field-gradient NMR

The low-temperature liquid-liquid phase separation of the partially miscible hexane-nitrobenzene mixture imbibed in porous glasses of different pore sizes from 7 to 130 nm has been studied using 1H NMR (nuclear magnetic resonance) cryoporometry and pulse field gradient NMR methods. The mixture was quenched below both its upper critical solution temperature (T(cr)) and the freezing point of nitrobenzene. The size distribution of frozen nitrobenzene domains was derived through their melting point suppression according to the Gibbs-Thompson relation. The obtained data reveal small initial droplets of nitrobenzene surrounded by hexane, which are created as the temperature is decreased below T(cr) and which thereafter coalesce by a droplet-diffusion mechanism. The inter-relation between the pore size and the found size distribution and shapes of nitrobenzene domains is discussed, as well as several aspects of molecular self-diffusion.