Study of phase separation of a binary fluid mixture in confined geometry

Confined microemulsions: pore diameter induced change of the phase behavior

In the present work, the temperature-dependent phase behavior of a C10E4 based microemulsion is studied in different meso-macroporous glasses, as a function of their pore diameter. The phase behavior in these pores is investigated by small-angle X-ray scattering (SAXS). The crucial parameter we discuss based on the SAXS results is the domain size of the bicontinuous phase. Using a simplified model to fit the scattering data, we can observe the microemulsion inside the pores. These experiments reveal a temperature-dependent change in domain sizes of the bicontinuous microemulsion only for large pores.

Machine-Learned Free Energy Surfaces for Capillary Condensation and Evaporation in Mesopores

Using molecular simulations, we study the processes of capillary condensation and capillary evaporation in model mesopores. To determine the phase transition pathway, as well as the corresponding free energy profile, we carry out enhanced sampling molecular simulations using entropy as a reaction coordinate to map the onset of order during the condensation process and of disorder during the evaporation process. The structural analysis shows the role played by intermediate states, characterized by the onset of capillary liquid bridges and bubbles. We also analyze the dependence of the free energy barrier on the pore width. Furthermore, we propose a method to build a machine learning model for the prediction of the free energy surfaces underlying capillary phase transition processes in mesopores.

Water/PEG Mixtures: Phase Behavior, Dynamics and Soft Confinement

Polyethylene glycol is water soluble and forms an eutectic system with water. The eutectic temperature is −19 °C for M=1500 g mol−1 and increases with molecular weight. The dielectric relaxation spectrum of the mixtures exhibits a strong loss maximum in ϵ″ (ω) similar to pure water. Relaxation time increases with the addition of PEG. Activation energies exhibit a maximum of 0.35 eV at molar fraction χp ≈0.2. This compares well with results on ethanol water mixtures. Adding PEG molecules to nanoscopic water droplets of inverse microemulsions has only small impact on the bending modulus κ of a non-ionic microemulsion. In AOT based microemulsions an increase or decrease of κ is found in dependence on the size of the droplets. This is in accordance with the variation of the dynamic percolation transition in the same systems.

Impact of complex topology of porous media on phase separation of binary mixtures

Porous materials, which are characterized by the large surface area and percolated nature crucial for transport, play an important role in many technological applications including battery, ion exchange, catalysis, microelectronics, medical diagnosis, and oil recovery. Phase separation of a mixture in such a porous structure should be strongly influenced by both surface wetting and strong geometrical confinement effects. Despite its fundamental and technological importance, however, this problem has remained elusive for a long time because of the difficulty associated with the complex geometry of pore structures. We overcome this by developing a novel phase-field model of two coupled order parameters, the composition field of a binary mixture and the density field of a porous structure. We find that demixing behavior in complex pore structures is severely affected by the topological characteristics of porous materials, contrary to the conventional belief that it can be inferred from the behavior in a simple cylindrical pore. Our finding not only reveals the physical mechanism of demixing in random porous structures but also has an impact on technological applications.

Phase separation of triethylamine and water in native and organically modified silica nanopores

A mixture of triethylamine and water is a lower critical solution temperature system that demixes (separates into individual phases) on heating. Differential scanning calorimetry has been applied to study the process of demixing in native and organically modified silica nanopores whose size varied from 4 to 30 nm. It has been found that in both types of nanopores, the temperature and enthalpy of demixing decrease significantly with decreasing the pore size. Isoconversional kinetic analysis has been utilized to determine the activation energy and pre-exponential factor of the process. It has been demonstrated that the depression of the transition temperature upon nanoconfinement is associated with acceleration of the process due to lowering of the activation energy. Nanoconfinement has also been found to lower the pre-exponential factor of the process that has been linked to a decrease in the molecular mobility.

Phase separation of binary nonadditive hard sphere fluid mixture confined in random porous media

I analyze the fluid-fluid phase separation of nonadditive hard sphere fluid mixture absorbed in random porous media. An equation of state is derived by using the perturbation theory to this complex system with quenched disorders. The results of this theory are in good agreement with those obtained from semi-grand canonical ensemble Monte Carlo simulations. The contact value of the fluid-fluid radial distribution functions of the reference which is the key point of the perturbation process is derived as well, the comparison against Monte Carlo simulations shows that it has an excellent accuracy.

Three-dimensional reconstruction of liquid phases in disordered mesopores usingin situsmall-angle scattering

Small-angle scattering of X-rays (SAXS) or neutrons is one of the few experimental methods currently available for thein situanalysis of phenomena in mesoporous materials at the mesoscopic scale. In the case of disordered mesoporous materials, however, the main difficulty of the method lies in the data analysis. A stochastic model is presented, which enables one to reconstruct the three-dimensional nanostructure of liquids confined in disordered mesopores starting from small-angle scattering data. This so-called plurigaussian model is a multi-phase generalization of clipped Gaussian random field models. Its potential is illustrated through the synchrotron SAXS analysis of a gel permeated with a critical nitrobenzene/hexane solution that is progressively cooled below its consolute temperature. The reconstruction brings to light a wetting transition whereby the nanostructure of the pore-filling liquids passes from wetting layers that uniformly cover the solid phase of the gel to plugs that locally occlude the pores. Using the plurigaussian model, the dewetting phenomenon is analyzed quantitatively at the nanometre scale in terms of changing specific interface areas, contact angle and specific length of the triple line.

Confinement-driven phase separation of quantum liquid mixtures

We report small-angle neutron scattering studies of liquid helium mixtures confined in Mobil Crystalline Material-41 (MCM-41), a porous silica glass with narrow cylindrical nanopores (d=3.4 nm). MCM-41 is an ideal model adsorbent for fundamental studies of gas sorption in porous media because its monodisperse pores are arranged in a 2D triangular lattice. The small-angle scattering consists of a series of diffraction peaks whose intensities are determined by how the imbibed liquid fills the pores. Pure (4)He adsorbed in the pores show classic, layer-by-layer film growth as a function of pore filling, leaving the long range symmetry of the system intact. In contrast, the adsorption of (3)He-(4)He mixtures produces a structure incommensurate with the pore lattice. Neither capillary condensation nor preferential adsorption of one helium isotope to the pore walls can provide the symmetry-breaking mechanism. The scattering is consistent with the formation of randomly distributed liquid-liquid microdomains ∼2.3 nm in size, providing evidence that confinement in a nanometer scale capillary can drive local phase separation in quantum liquid mixtures.

Phase separation of a binary liquid in anodic aluminium oxide templates: a structural study by small angle neutron scattering

The radial nanostructure of the binary liquid triethylamine/water confined in 60 nm diameter independent cylindrical pores of anodic aluminium oxide membranes is studied by small angle neutron scattering. It is shown that composition inhomogeneities are present in the confined mixtures well below the bulk critical point. An analysis of the neutron scattering form factor reveals the existence of an adsorbed water layer of a few nanometers at the liquid/alumina interface, coexisting with a TEA-rich phase in the core.

Capillary condensation in cylindrical pores: Monte Carlo study of the interplay of surface and finite size effects

When a fluid that undergoes a vapor to liquid transition in the bulk is confined to a long cylindrical pore, the phase transition is shifted (mostly due to surface effects at the walls of the pore) and rounded (due to finite size effects). The nature of the phase coexistence at the transition depends on the length of the pore: for very long pores, the system is axially homogeneous at low temperatures. At the chemical potential where the transition takes place, fluctuations occur between vapor- and liquidlike states of the cylinder as a whole. At somewhat higher temperatures (but still far below bulk criticality), the system at phase coexistence is in an axially inhomogeneous multidomain state, where long cylindrical liquid- and vaporlike domains alternate. Using Monte Carlo simulations for the Ising/lattice gas model and the Asakura-Oosawa model of colloid-polymer mixtures, the transition between these two different scenarios is characterized. It is shown that the density distribution changes gradually from a double-peak structure to a triple-peak shape, and the correlation length in the axial direction (measuring the equilibrium domain length) becomes much smaller than the cylinder length. The (rounded) transition to the disordered phase of the fluid occurs when the axial correlation length has decreased to a value comparable to the cylinder diameter. It is also suggested that adsorption hysteresis vanishes when the transition from the simple domain state to the multidomain state of the cylindrical pore occurs. We predict that the difference between the pore critical temperature and the hysteresis critical temperature should increase logarithmically with the length of the pore.

Rounding of phase transitions in cylindrical pores

Phase transitions of systems confined in long cylindrical pores (capillary condensation, wetting, crystallization, etc.) are intrinsically not sharply defined but rounded. The finite size of the cross section causes destruction of long range order along the pore axis by spontaneous nucleation of domain walls. This rounding is analyzed for two models (Ising or lattice gas and Asakura-Oosawa model for colloid-polymer mixtures) by Monte Carlo simulations and interpreted by a phenomenological theory. We show that characteristic differences between the behavior of pores of finite length and infinitely long pores occur. In pores of finite length a rounded transition occurs first, from phase coexistence between two states towards a multidomain configuration. A second transition to the axially homogeneous phase follows near pore criticality.

Critical behavior of symmetrical fluid mixtures in random pores

We study the liquid-liquid demixing of a binary mixture with a symmetrical coupling to the quenched disorder by means of computer simulation. The critical point in the thermodynamic limit is estimated both by assuming the knowledge of the critical exponents and independently of them. The finite-size scaling analysis of the susceptibilities and the values of the critical amplitudes show that the universality class of the fluid mixture is compatible with the diluted quenched-Ising model. Our findings extend the class of systems exhibiting the same critical behavior of diluted antiferromagnets.

Condensation phenomena in nanopores: A Monte Carlo study

The nonequilibrium dynamics of condensation phenomena in nanopores is studied via Monte Carlo simulations of a lattice-gas model. Hysteretic behavior of the particle density as a function of the density of a reservoir is obtained for various pore geometries in two and three dimensions. The shape of the hysteresis loops depend on the characteristics of the pore geometry. The evaporation of particles from a pore can be fitted to a stretched exponential decay of the particle density. Phase-separation dynamics inside the pore is effectively described by a random walk of the non-wetting phases. Domain evolution is significantly slowed down in the presence of a random wall-particle potential and gives rise to a temperature-dependent growth exponent. A geometric roughness of the pore wall only delays the onset of a pure domain growth.

Local structure of a phase-separating binary mixture in a mesoporous glass matrix studied by small-angle neutron scattering

The mesoscopic structure of the binary system isobutyric acid + heavy water (D(2)O) confined in a porous glass (controlled-pore silica glass, mean pore width ca. 10 nm) was studied by small-angle neutron scattering at off-critical compositions in a temperature range above and below the upper critical solution point. The scattering data were analyzed in terms of a structure factor model similar to that proposed by Formisano and Teixeira [Eur. Phys. J. E 1, 1 (2000)], but allowing for both Ornstein-Zernike-type composition fluctuations and domainlike structures in the microphase-separated state of the pore liquid. The results indicate that the phase separation in the pores is shifted by ca. 10 K and spread out in temperature. Microphase separation is pictured as a transition from partial segregation at high temperature, due to the strong preferential adsorption of water at the pore wall, to a tube or capsule configuration of the two phases at low temperatures, depending on the overall composition of the pore liquid. Results for samples in which the composition of the pore liquid can vary with temperature due to equilibration with extra-pore liquid are consistent with this picture.

Computer simulation of the phase diagram for a fluid confined in a fractal and disordered porous material

We present a grand canonical Monte Carlo simulation study of the phase diagram of a Lennard-Jones fluid adsorbed in a fractal and highly porous aerogel. The gel environment is generated from an off-lattice diffusion limited cluster-cluster aggregation process. Simulations have been performed with the multicanonical ensemble sampling technique. The biased sampling function has been obtained by histogram reweighting calculations. Comparing the confined and the bulk system liquid-vapor coexistence curves we observe a decrease of both the critical temperature and density in qualitative agreement with experiments and other Monte Carlo studies on Lennard-Jones fluids confined in random matrices of spheres. At variance with these numerical studies we do not observe upon confinement a peak on the liquid side of the coexistence curve associated with a liquid-liquid phase coexistence. In our case only a shouldering of the coexistence curve appears upon confinement. This shoulder can be associated with high density fluctuations in the liquid phase. The coexisting vapor and liquid phases in our system show a high degree of spatial disorder and inhomogeneity.

Density fluctuations near the liquid-gas critical point of a confined fluid

We report the results of an experimental study of the effect of a dilute silica network on liquid-gas critical phenomena in carbon dioxide (CO2). Using small-angle neutron scattering, we measured the correlation length of the density fluctuations in bulk (xi(bulk)) and confined CO2 (xi(conf)) as a function of temperature and average fluid density. We find that quenched disorder induced by an aerogel suppresses density fluctuations: xi(conf) loses the Ising model divergence characteristic of xi(bulk) and does not exceed the size of pores in the homogeneous region.

De-mixing dynamics of a binary liquid system in a controlled-pore glass. A neutron spin-echo spectroscopy and small angle neutron scattering study

The temperature-induced microphase separation of the binary liquid system iso-butyric acid+heavy water (iBA + D(2)O) in a mesoporous silica glass (CPG-10-75) of nominal pore width 7.5 nm was investigated by neutron spin-echo spectroscopy (NSE) and neutron small-angle scattering (SANS). Two mixtures of different composition were studied at different scattering angles at temperatures above and below the bulk phase transition temperature. The phase separation in the pore space is found to occur at a lower temperature than the bulk transition and extends over a significant temperature range. The effective diffusion coefficient derived from NSE at low scattering angles is found to decrease by one order of magnitude from 70 degrees C to 20 degrees C. This observation is attributed to the growing size of concentration fluctuations having a cut-off at ca. 8 nm, which corresponds to the mean pore size. The dynamics of the concentration fluctuations appears to be strongly influenced by the confinement in the pores, as it differs strongly from the bulk behaviour. These results are consistent with the preliminary results of the SANS study.

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.

Effects of pore walls and randomness on phase transitions in porous media

We study spin models within the mean field approximation to elucidate the topology of the phase diagrams of systems modeling the liquid-vapor transition and the separation of 3He-4He mixtures in periodic porous media. These topologies are found to be identical to those of the corresponding random field and random anisotropy spin systems with a bimodal distribution of the randomness. Our results suggest that the presence of walls (periodic or otherwise) are a key factor determining the nature of the phase diagram in porous media.

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.

Capillary condensation as a morphological transition

The process of capillary condensation/evaporation in cylindrical pores is considered within the idea of symmetry breaking. Capillary condensation/evaporation is treated as a morphological transition between the wetting film configurations of different symmetry. We considered two models: (i) the classical Laplace theory of capillarity and (ii) the Derjaguin model which takes into account the surface forces expressed in terms of the disjoining pressure. Following the idea of Everett and Haynes, the problem of condensation/evaporation is considered as a transition from bumps/undulations to lenses. Using the method of phase portraits, we discuss the mathematical mechanisms of this transition hidden in the Laplace and Derjaguin equations. Analyzing the energetic barriers of the bump and lens formation, it is shown that the bump formation is a prerogative of capillary condensation: for the vapor-liquid transition in a pore, the bump plays the same role as the spherical nucleus in a bulk fluid. We show also that the Derjaguin model admits a variety of interfacial configurations responsible for film patterning at specific conditions.

Effect of templated quenched disorder on fluid phase equilibrium

Templating offers a means to direct the structure of quenched disorder. We show here that changes in phase equilibrium due to the presence of quenched disorder can themselves be altered by templating. We calculate the phase diagram of a fluid in a collection of template-directed, quenched particles by solving a set of replica Ornstein-Zernike equations within the mean spherical approximation and show templating to enhance phase behavior, that is, shift the phase envelope upward from its location for a nontemplated system of identical available volume. This enhancement is due to an augmented number of fluid-fluid interactions.

Phase diagram of a symmetric binary fluid in a porous matrix

The phase behavior of a binary symmetric fluid in thermal equilibrium with a porous matrix has been studied with the optimized random phase approximation and grand canonical Monte Carlo simulations. Depending on the matrix properties and the matrix-fluid and fluid-fluid interactions we find three types of phase diagram characterized by a tricritical point, a tricritical point with a triple point, or a critical end point. Small changes in the properties of the matrix or in the interactions are demonstrated to lead to drastic modifications of the phase diagram of the fluid, in qualitative agreement with observations in experimental studies. We show, in particular, that the change between the different types of phase diagram is triggered not only by the fluid-fluid interactions (internal parameters) but also by the properties of the matrix and of the matrix-fluid potentials (external parameters).

Simple views on critical binary liquid mixtures in porous glass

A simple scenario, different from previous attempts, is proposed to resolve the problem of the slow phase separation dynamics of binary liquid mixtures confined in porous Vycor glass. We demonstrate that simply mutual diffusion, renormalized by critical composition fluctuations and geometrical hindrance of the porous glass, accounts for the slow phase separation kinetics. Capillary invasion studies of porous Vycor glass by the critical isobutyric acid-water mixture, close to the consolute solution temperature, corroborate our analysis.

Monte Carlo approach to the gas-liquid transition in porous materials

The gas-liquid transition of a "quenched-annealed"(QA) system is studied by grand-canonical Monte Carlo simulation. The "quenched" particles are hard spheres within configurations chosen randomly from those of an equilibrium hard-sphere system at given density. The fluid particles interact with the matrix particles by a hard-core potential and with each other by a hard-core potential and an additional potential of a Lennard-Jones type. Our results are in good qualitative agreement with various theoretical approaches. With increasing matrix density the critical temperature is lowered compared to that of the bulk system and the gap between the gas and liquid densities narrowed. Our simulations confirm, for this QA system, the possibility of two fluid-fluid transitions substituting for the unique gas-liquid transition of the bulk system. They demonstrate the necessity to average over a significant number of matrix realizations in order to obtain a quantitative location of the phase coexistence lines.