Fluid-fluid Transition Research Articles

High-pressure phase equilibrium data of carbon dioxide/food-relevant systems (2011-2022): Experimental methods, multiphase behavior, thermodynamic modeling, and applications

High-pressure phase equilibrium data involving food-type substances and pure or cosolvent-modified pressurized carbon dioxide (CO2) published in the last decade are reviewed. Experimental data regarding binary, ternary, and multicomponent systems studied at 293–373 K and pressures up to 55 MPa are compiled. CO2/lipid-like substances (fatty acids, esters, triacylglycerols), essential oils, and their components are the most studied food systems, followed by systems containing biopolymers. The Vapor-Liquid Equilibrium is the most observed fluid-fluid transition at experimental conditions. The cosolvent effect of organic solvents in removing immiscibility gaps like the Vapor-Liquid-Liquid and the Liquid-Liquid boundaries, the antisolvent role of CO2 to some solutes, and its plasticizing effect on polymers and solid lipids are also highlighted in the context of solid-fluid equilibrium. The Peng and Robinson Equation of State coupled with the Van der Waals type mixing rules is still the main thermodynamic model applied for modeling phase equilibrium of these systems, despite some correlative limitations regarding multicomponent mixtures containing polar molecules. Combining classical Equations of State with the associative term and predictive models as Group Contribution methods is a reliable tool to improve the modeling of systems with complex phase behaviors. Fluid phase immiscibility regions and experimental conditions that favor the presence of a solid phase are mainly useful to set the appropriate operating point based on phase diagrams and the observed equilibria. For many applications, phase equilibrium is relevant in processing food-type substances destined for food or food-related applications and human consumption by processes employing CO2 as an extraction solvent, a fractionating, dispersant, precipitative, or atomizing (antisolvent) agent, or as a reaction medium for biocatalysis.

Condensation of Methane in the Metal-Organic Framework IRMOF-1: Evidence for Two Critical Points.

Extensive grand canonical Monte Carlo simulations in combination with successive umbrella sampling are used to investigate the condensation of methane in the nanoporous crystalline material IRMOF-1. Two different types of novel condensation transitions are found, each of them ending in a critical point: (i) a fluid-fluid transition at higher densities (the analog of the liquid-gas transition in the bulk) and (ii) a phase transition at low densities on the surface of the IRMOF-1 structure. The nature of these transitions is different from the usual capillary condensation in thin films and cylindrical pores where the coexisting phases are confined in one or two of the three spatial dimensions. In contrast to that, in IRMOF-1 the different phases can be described as bulk phases that are inhomogeneous due to the presence of the metal-organic framework. As a consequence, the condensation transitions in IRMOF-1 belong to the three-dimensional (3D) Ising universality class.

Novel structure formation of a phase separating colloidal fluid in a ratchet potential

Based on Dynamical Density Functional Theory (DDFT) we investigate a binary mixture of interacting Brownian particles driven over a substrate via a one-dimensional ratchet potential. The particles are modeled as soft spheres where one component carries a classical Heisenberg spin. In the absence of a substrate field, the system undergoes a first-order fluid-fluid demixing transition driven by the spin-spin interaction. We demonstrate that the interplay between the intrinsic spinodal decomposition and time-dependent external forces leads to a novel dynamical instability where stripes against the symmetry of the external potential form. This structural transition is observed for a broad range of parameters related to the ratchet potential. Moreover, we find intriguing effects for the particle transport.

Demixing and surface-induced population inversion of a binary Gaussian mixture in nanopores

We investigate the demixing and the surface-induced population inversion of a binary Gaussian mixture, which exhibits the fluid-fluid demixing transition of the bulk system, in a nanopore. The result shows that the population inversion is the typical shift of the first-order fluid-fluid demixing transition due to the confinement effect. The population inversion line is affected by the crossinteraction between unlike species σ12, the mole fraction of species x2, and the slit width H. At a low cross-interaction between unlike species, the population inversion line for x2 xc. At a high cross-interaction, for x2 xc, it is shifted toward a lower density.

Communication: Virial coefficients and demixing in highly asymmetric binary additive hard-sphere mixtures

The problem of demixing in a binary fluid mixture of highly asymmetric additive hard spheres is revisited. A comparison is presented between the results derived previously using truncated virial expansions for three finite size ratios with those that one obtains with the same approach in the extreme case in which one of the components consists of point particles. Since this latter system is known not to exhibit fluid-fluid segregation, the similarity observed for the behavior of the critical constants arising in the truncated series in all instances, while not being conclusive, may cast serious doubts as to the actual existence of a demixing fluid-fluid transition in disparate-sized binary additive hard-sphere mixtures.

Nonadditive Penetrable Mixtures in Nanopores: Surface-Induced Population Inversion

We investigate the surface-induced population inversion of the nonadditive penetrable mixtures which exhibits the fluid-fluid demixing transition of the bulk system due to the confinement effect. The result shows that the population inversions are strongly affected by the extra repulsion between unlike species, the mole fraction of species, the width of nanopores, and the nonadditive walls. The extra repulsion between unlike species in a confined system increases the contact density of both species at the wall and promotes the population inversion in nanopores. The population inversion is the typical shift first-order fluid-fluid demixing transition due to the confinement effect in nanopores. The population inversions are only observed in nanopores with finite widths. The population inversion line is shifted toward a higher fluid density with decreasing width of the nanopores and lies slightly in lower density compared with the coexistence curves of the bulk system. The nonadditive wall for the big particles leads to the population inversion in lower density compared with that of the nonadditive wall for the small particles. The population inversion line is terminated at a lower mole fraction.

Self-assembly mechanism in colloids: perspectives from statistical physics

AbstractMotivated by recent experimental findings in chemical synthesis of colloidal particles, we draw an analogy between self-assembly processes occurring in biological systems (e.g. protein folding) and a new exciting possibility in the field of material science. We consider a self-assembly process whose elementary building blocks are decorated patchy colloids of various types, that spontaneously drive the system toward a unique and predetermined targeted macroscopic structure. To this aim, we discuss a simple theoretical model — the Kern-Frenkel model — describing a fluid of colloidal spherical particles with a pre-defined number and distribution of solvophobic and solvophilic regions on their surface. The solvophobic and solvophilic regions are described via a short-range square-well and a hard-sphere potentials, respectively. Integral equation and perturbation theories are presented to discuss structural and thermodynamical properties, with particular emphasis on the computation of the fluid-fluid (or gas-liquid) transition in the temperaturedensity plane. The model allows the description of both one and two attractive caps, as a function of the fraction of covered attractive surface, thus interpolating between a square-well and a hard-sphere fluid, upon changing the coverage. By comparison with Monte Carlo simulations, we assess the pros and the cons of both integral equation and perturbation theories in the present context of patchy colloids, where the computational effort for numerical simulations is rather demanding.

A numerical test of a high-penetrability approximation for the one-dimensional penetrable-square-well model

The one-dimensional penetrable-square-well fluid is studied using both analytical tools and specialized Monte Carlo simulations. The model consists of a penetrable core characterized by a finite repulsive energy combined with a short-range attractive well. This is a many-body one-dimensional problem, lacking an exact analytical solution, for which the usual van Hove theorem on the absence of phase transition does not apply. We determine a high-penetrability approximation complementing a similar low-penetrability approximation presented in previous work. This is shown to be equivalent to the usual Debye-Hückel theory for simple charged fluids for which the virial and energy routes are identical. The internal thermodynamic consistency with the compressibility route and the validity of the approximation in describing the radial distribution function is assessed by a comparison against numerical simulations. The Fisher-Widom line separating the oscillatory and monotonic large-distance behaviors of the radial distribution function is computed within the high-penetrability approximation and compared with the opposite regime, thus providing a strong indication of the location of the line in all possible regimes. The high-penetrability approximation predicts the existence of a critical point and a spinodal line, but this occurs outside the applicability domain of the theory. We investigate the possibility of a fluid-fluid transition by the Gibbs ensemble Monte Carlo techniques, not finding any evidence of such a transition. Additional analytical arguments are given to support this claim. Finally, we find a clustering transition when Ruelle's stability criterion is not fulfilled. The consequences of these findings on the three-dimensional phase diagrams are also discussed.

Virial coefficients, thermodynamic properties, and fluid-fluid transition of nonadditive hard-sphere mixtures

Different theoretical approaches for the thermodynamic properties and the equation of state for multicomponent mixtures of nonadditive hard spheres in d dimensions are presented in a unified way. These include the theory by Hamad, our previous formulation, the original MIX1 theory, a recently proposed modified MIX1 theory, as well as a nonlinear extension of the MIX1 theory proposed in this paper. Explicit expressions for the compressibility factor, Helmholtz free energy, and second, third, and fourth virial coefficients are provided. A comparison is carried out with recent Monte Carlo data for the virial coefficients of asymmetric mixtures and with available simulation data for the compressibility factor, the critical consolute point, and the liquid-liquid coexistence curves. The merits and limitations of each theory are pointed out.

Gelation and phase coexistence in colloidal suspensions with short-range forces: Generic behavior versus specificity

The interplay between physical gelation and equilibrium phase transitions in asymmetric binary mixtures is analyzed from the effective fluid approach, in which the big particles interact via a short-range effective attraction beyond the core due to the depletion mechanism. The question of the universality of the scenario for dynamical arrest is then addressed. The comparison of the phase diagrams of the hard-sphere mixture and the Asakura-Oosawa models at various size ratios shows that strong specificity is observed for nonideal depletants. In particular, equilibrium gelation, without the competition with fluid-fluid transition is possible in mixtures of hard-sphere colloids. This is interpreted from the specificities of the effective potential, such as its oscillatory behavior and its complex variation with the physical parameters. The consequences on the dynamical arrest and the fluid-fluid transition are then investigated by considering in particular the role of the well at contact and the first repulsive barrier. This is done for the actual effective potential in the hard-sphere mixture and for a square well and shoulder model, which allows a separate discussion of the role of the different parameters, in particular on the localization length and the escape time. This study is next extended to mixtures of "hard-sphere-like" colloids with residual interactions. It confirms the trends relative to equilibrium gelation and illustrates a diversity of the phase behavior well beyond the scenarios expected from simple models.

Penetrable-square-well fluids: Analytical study and Monte Carlo simulations

We study structural and thermophysical properties of a one-dimensional classical fluid made of penetrable spheres interacting via an attractive square-well potential. Penetrability of the spheres is enforced by reducing from infinite to finite the repulsive energy barrier in the pair potentials As a consequence, an exact analytical solution is lacking even in one dimension. Building upon previous exact analytical work in the low-density limit [A. Santos, R. Fantoni, and A. Giacometti, Phys. Rev. E 77, 051206 (2008)], we propose an approximate theory valid at any density and in the low-penetrable regime. By comparison with specialized Monte Carlo simulations and integral equation theories, we assess the regime of validity of the theory. We investigate the degree of inconsistency among the various routes to thermodynamics and explore the possibility of a fluid-fluid transition. Finally we locate the dependence of the Fisher-Widom line on the degree of penetrability. Our results constitute the first systematic study of penetrable spheres with attractions as a prototype model for soft systems.

Solid-solid transition of the size-polydisperse hard sphere system

The solid-solid coexistence of a polydisperse hard sphere system is studied by using the Monte Carlo simulation. The results show that for large enough polydispersity the solid-solid coexistence state is more stable than the single-phase solid. The two coexisting solids have different composition distributions but the same crystal structure. Moreover, there is evidence that the solid-solid transition terminates in a critical point as in the case of the fluid-fluid transition.

Triple Point on the Melting Curve and Polymorphism of Nitrogen at High Pressure

Raman spectra of solid and fluid nitrogen to pressures up to 120 GPa and temperatures up to 2500 K reveal that the melting line exhibits a maximum near 70 GPa, followed by a triple point near 87 GPa, after which the melting temperature rises again. Fluid nitrogen remains molecular over the entire pressure range studied, and there is no sign of a fluid-fluid transition. Solid phases obtained on quenching from the melt above 48 GPa are identical to the recently discovered iota and zeta' phases. We find that kinetics plays a major role in the experimentally observed phase changes and account for the metastability of various crystalline molecular phases and the existence of an amorphous single bonded eta-N.

Patchy sticky hard spheres: Analytical study and Monte Carlo simulations

We consider a fluid of hard spheres bearing one or two uniform circular adhesive patches, distributed so as not to overlap. Two spheres interact via a "sticky" Baxter potential if the line joining the centers of the two spheres intersects a patch on each sphere, and via a hard sphere potential otherwise. We analyze the location of the fluid-fluid transition and of the percolation line as a function of the size of the patch (the fractional coverage of the sphere's surface) and of the number of patches within a virial expansion up to third order and within the first two terms (C0 and C1) of a class of closures Cn hinging on a density expansion of the direct correlation function. We find that the locations of the two lines depend sensitively on both the total adhesive coverage and its distribution. The treatment is almost fully analytical within the chosen approximate theory. We test our findings by means of specialized Monte Carlo simulations and find the main qualitative features of the critical behavior to be well captured in spite of the low density perturbative nature of the closure. The introduction of anisotropic attractions into a model suspension of spherical particles is a first step toward a more realistic description of globular proteins in solution.

Large attractive depletion interactions in soft repulsive–sphere binary mixtures

We consider binary mixtures of soft repulsive spherical particles and calculate the depletion interaction between two big spheres mediated by the fluid of small spheres, using different theoretical and simulation methods. The validity of the theoretical approach, a virial expansion in terms of the density of the small spheres, is checked against simulation results. Attention is given to the approach toward the hard-sphere limit and to the effect of density and temperature on the strength of the depletion potential. Our results indicate, surprisingly, that even a modest degree of softness in the pair potential governing the direct interactions between the particles may lead to a significantly more attractive total effective potential for the big spheres than in the hard-sphere case. This might lead to significant differences in phase behavior, structure, and dynamics of a binary mixture of soft repulsive spheres. In particular, a perturbative scheme is applied to predict the phase diagram of an effective system of big spheres interacting via depletion forces for a size ratio of small and big spheres of 0.2; this diagram includes the usual fluid-solid transition but, in the soft-sphere case, the metastable fluid-fluid transition, which is probably absent in hard-sphere mixtures, is close to being stable with respect to direct fluid-solid coexistence. From these results, the interesting possibility arises that, for sufficiently soft repulsive particles, this phase transition could become stable. Possible implications for the phase behavior of real colloidal dispersions are discussed.

Theoretical prediction of multiple fluid-fluid transitions in monocomponent fluids

The authors use the analytical equation of state obtained by the discrete perturbation theory [A. L. Benavides and A. Gil-Villegas, Mol. Phys. 97, 1225 (1999)] to study the phase diagram of fluids with discrete spherical potentials formed by a repulsive square-shoulder plus an attractive square-well interaction (SS+SW). This interaction is characterized by the usual energy and size parameters plus three dimensionless parameters: two of them measuring the widths of the SS and the SW and the third the relative height of the SS. The matter of interest is that, for certain values of the interaction parameters, the SS+SW systems exhibit more than one first-order fluid-fluid transition. The evidence that several real substances (such as water, phosphorus, carbon, and silica, among others) exhibit an extra liquid-liquid transition has drawn interest into the study of interactions responsible for this behavior. The simple SS+SW fluid is one of the systems that, in spite of being spherically symmetric, shows multiple fluid-fluid transitions. In this work the authors investigate systematically the effect on the phase diagram of varying the interaction parameters. The use of an analytical free-energy equation gives a clear thermodynamic picture of the emergence of different types of critical points, throwing new light on the phase behavior of these fluids and thus clarifying previous results obtained by other techniques. The interplay of attractive and repulsive forces with several scale lengths produces very rich phase diagrams, including cases with three critical points. The region of the interaction-parameter space where multiple critical points appear is mapped for various families of interactions.

Phase behavior and structure of model colloid-polymer mixtures confined between two parallel planar walls

Using Gibbs ensemble Monte Carlo simulations and density functional theory we investigate the fluid-fluid demixing transition in inhomogeneous colloid-polymer mixtures confined between two parallel plates with separation distances between one and ten colloid diameters covering the complete range from quasi-two-dimensional to bulklike behavior. We use the Asakura-Oosawa-Vrij model in which colloid-colloid and colloid-polymer interactions are hard-sphere like, while the pair potential between polymers vanishes. Two different types of confinement induced by a pair of parallel walls are considered--namely, either through two hard walls or through two semipermeable walls that repel colloids but allow polymers to freely penetrate. For hard (semipermeable) walls we find that the capillary binodal is shifted towards higher (lower) polymer fugacities and lower (higher) colloid fugacities as compared to the bulk binodal; this implies capillary condensation (evaporation) of the colloidal liquid phase in the slit. A macroscopic treatment is provided by a symmetric Kelvin equation for general binary mixtures based on the proximity in chemical potentials of statepoints at capillary coexistence and the reference bulk coexistence. Results for capillary binodals compare well with those obtained from the classic version of the Kelvin equation due to [Evans and Marini Bettolo Marconi, J. Chem. Phys. 86, 7138 (1987)] and are quantitatively accurate away from the fluid-fluid critical point, even at small wall separations. However, the significant shift of the critical polymer fugacity towards higher values upon increasing confinement, as found in simulations, is not reproduced. For hard walls the density profiles of polymers and colloids inside the slit display oscillations due to packing effects for all statepoints. For semipermeable walls either similar structuring or flat profiles are found, depending on the statepoint considered.

Miscibility of small colloidal spheres with large polymers in good solvent

Nearly athermal colloid-polymer mixtures were studied in the "protein limit." A fluid-fluid transition was observed in mixtures of stearyl-alcohol-coated silica particles and large polystyrene coils in toluene. The ratios of the polymer radius of gyration to the particle radii were q=4.1 and q=5.2. The binodal curves and the critical points were determined. Turbidity measurements and analysis for one set of particles allowed the systems to be mapped onto hard sphere-polymer mixtures. A comparison with recent predictions for the miscibility of model mixtures shows that the experimental binodals lie between the two extreme results for ideal and interacting polymers. The critical colloid volume fraction is also found to decrease with increasing size ratios.

Decay of correlation functions in hard-sphere mixtures: structural crossover.

We investigate the decay of pair correlation functions in a homogeneous (bulk) binary mixture of hard spheres. At a given state point the asymptotic decay r-->infinity of all three correlation functions is governed by a common exponential decay length and a common wavelength of oscillations. Provided the mixture is sufficiently asymmetric, size ratios q less than or approximately 0.7, we find that the common wavelength reflects either the size of the small or that of the big spheres. By analyzing the (complex) poles of the partial structure factors we find a sharp structural crossover line in the phase diagram. On one side of this line the common wavelength is approximately the diameter of the smaller sized spheres whereas on the other side it is approximately the diameter of the bigger ones; the wavelength of the longest ranged oscillations changes discontinuously at the structural crossover line. Using density functional theory and Monte Carlo simulations we show that structural crossover also manifests itself in the intermediate range behavior of the pair correlation functions and we comment on the relevance of this observation for real (colloidal) mixtures. In highly asymmetric mixtures, q< or =0.1, where there is metastable fluid-fluid transition, we find a Fisher-Widom line with two branches. This line separates a region of the phase diagram where the decay of pair correlations is oscillatory from one in which it is monotonic.

Phase behavior of aqueous solutions containing dipolar proteins from second-order perturbation theory.

Due to the interplay of Coulombic repulsion and attractive dipolar and van der Waals interactions, solutions of globular proteins display a rich variety of phase behavior featuring fluid-fluid and fluid-solid transitions that strongly depend on solution pH and salt concentration. Using a simple model for charge, dispersion and dipole-related contributions to the interprotein potential, we calculate phase diagrams for protein solutions within the framework of second-order perturbation theory. For each phase, we determine the Helmholtz energy as the sum of a hard-sphere reference term and a perturbation term that reflects both the electrostatic and dispersion interactions. Dipolar effects can induce fluid-fluid phase separation or crystallization even in the absence of any significant dispersion attraction. Because dissolved electrolytes screen the charge-charge repulsion more strongly than the dipolar attraction, the ionic strength dependence of the potential of mean force can feature a minimum at intermediate ionic strengths offering an explanation for the observed nonmonotonic dependence of the phase behavior on salt concentration. Inclusion of correlations between charge-dipole and dipole-dipole interactions is essential for a reliable calculation of phase diagrams for systems containing charged dipolar proteins and colloids.