Phase diagram of a symmetric binary fluid in a porous matrix

Deformable hard particles confined in a disordered porous matrix

With suitably designed Monte Carlo simulations, we have investigated the properties of mobile, impenetrable, yet deformable particles that are immersed into a porous matrix, the latter one realized by a frozen configuration of spherical particles. By virtue of a model put forward by Batista and Miller [Phys. Rev. Lett. 105, 088305 (2010)], the fluid particles can change in their surroundings, formed by other fluid particles or the matrix particles, their shape within the class of ellipsoids of revolution; such a change in shape is related to a change in energy, which is fed into suitably defined selection rules in the deformation "moves" of the Monte Carlo simulations. This concept represents a simple yet powerful model of realistic, deformable molecules with complex internal structures (such as dendrimers or polymers). For the evaluation of the properties of the system, we have used the well-known quenched-annealed protocol (with its characteristic double average prescription) and have analyzed the simulation data in terms of static properties (the radial distribution function and aspect ratio distribution of the ellipsoids) and dynamic features (notably the mean squared displacement). Our data provide evidence that the degree of deformability of the fluid particles has a distinct impact on the aforementioned properties of the system.

Connect the Thermodynamics of Bulk and Confined Fluids: Confinement-Adsorption Scaling

A fluid (a gas or a liquid) adsorbed in a porous material can behave very differently from its bulk counterpart. The advent of various synthesized materials with nanopores and their wide applications have provided strong impetus for studying fluids in confinement because our current understanding is still incomplete. From a large number of Monte Carlo simulations, we found a scaling relation that allows for connecting some thermodynamic properties (chemical potential, free energy per particle, and grand potential per particle) of a confined fluid to the bulk ones. Upon rescaling the adsorbed fluid density, the adsorption isotherms for many different confining environments collapse to the corresponding bulk curve. We also reveal the intimate connection of the reported scaling relation to Gibbs theory of inhomogeneous fluids and morphological thermodynamics. The advance in our understanding of confined fluids, gained from this study, also opens attractive perspectives for circumventing experimental difficulty for directly measuring some fluid thermodynamic properties in nanoporous materials.

Fluids in porous media. IV. Quench effect on chemical potential

It appears to be a common sense to measure the crowdedness of a fluid system by the densities of the species constituting it. In the present work, we show that this ceases to be valid for confined fluids under some conditions. A quite thorough investigation is made for a hard sphere (HS) fluid adsorbed in a hard sphere matrix (a quench-annealed system) and its corresponding equilibrium binary mixture. When fluid particles are larger than matrix particles, the quench-annealed system can appear much more crowded than its corresponding equilibrium binary mixture, i.e., having a much higher fluid chemical potential, even when the density of each species is strictly the same in both systems, respectively. We believe that the insight gained from this study should be useful for the design of functionalized porous materials.

Vapour-liquid phase diagram for an ionic fluid in a random porous medium

We study the vapour-liquid phase behaviour of an ionic fluid confined in a random porous matrix formed by uncharged hard sphere particles. The ionic fluid is modelled as an equimolar binary mixture of oppositely charged equisized hard spheres, the so-called restricted primitive model (RPM). Considering the matrix-fluid system as a partly-quenched model, we develop a theoretical approach which combines the method of collective variables with the extension of the scaled-particle theory (SPT) for a hard-sphere fluid confined in a disordered hard-sphere matrix. The approach allows us to formulate the perturbation theory using the SPT for the description of the thermodynamics of the reference system. The phase diagrams of the RPM in matrices of different porosities and for different size ratios of matrix and fluid particles are calculated in the random-phase approximation and also when the effects of higher-order correlations between ions are taken into account. Both approximations correctly reproduce the basic effects of porous media on the vapour-liquid phase diagram, i.e. with a decrease of porosity the critical point shifts towards lower fluid densities and lower temperatures and the coexistence region gets narrower. For the fixed matrix porosity, both the critical temperature and the critical density increase with an increase of size of matrix particles and tend to the critical values of the bulk RPM.

Fluids in porous media: the case of neutral walls

The bulk phase behavior of a fluid is typically altered when the fluid is brought into confinement by the walls of a random porous medium. Inside the porous medium, phase-transition points are shifted, or may disappear altogether. A crucial determinant is how the walls interact with the fluid particles. In this work, we consider the situation whereby the walls are neutral with respect to the liquid and vapor phases. In order to realize the condition of strict neutrality, we use a symmetric binary mixture inside a porous medium that interacts identically with mixture species. Monte Carlo simulations are then used to obtain the phase behavior. Our main finding is that, in the presence of the porous medium, a liquid-vapor critical point still exists. At the critical point, the distribution of the order parameter remains scale invariant, but self-averaging is violated. These findings provide further evidence that random confinement by neutral walls induces critical behavior of the random Ising model (i.e., Ising models with dilution type disorder, where the disorder couples to the energy).

Reentrant miscibility in two-dimensional symmetrical mixtures

The Monte Carlo simulation method in the grand canonical ensemble is used to study the phase behavior of two-dimensional symmetrical binary mixtures of Lennard-Jones particles with negative nonadditivity and the weaker interaction between the pairs of unlike than between the pairs of like particles. We have determined the evolution of the phase diagram topology when the parameters describing the interaction between unlike particles vary. It has been found that such systems may exhibit reentrant miscibility in the liquid and the solid phases.

The phase behavior of two-dimensional symmetrical mixtures

Using Monte Carlo simulation methods in the grand canonical and semigrand canonical ensembles, we study the phase behavior of two-dimensional symmetrical binary mixtures of Lennard-Jones particles. We discuss the interplay between the demixing transition in a liquid and the freezing in detail. Phase diagrams for several systems characterized by different parameters describing interactions in the system are presented. It is explicitly demonstrated that different scenarios involving demixing and freezing transitions, described in our earlier paper [A. Patrykiejew and S. Sokołowski, Phys. Rev. E, 81, 012501 (2010)], are possible. In one class of systems, the λ-line representing a continuous demixing transition in a liquid phase starts at the liquid side of either the vapor-liquid or liquid-solid coexistence. The second class involves the systems in which the λ-line begins at the liquid side of the vapor-liquid coexistence, in the lower critical end point, and then terminates at the liquid side of the liquid-solid coexistence, in the upper critical end point. It is also shown that in such systems the solid phase may undergo a demixing transition at the temperature above the upper critical end point.

Impact of random obstacles on the dynamics of a dense colloidal fluid

Using molecular dynamics simulations we study the slow dynamics of a colloidal fluid annealed within a matrix of obstacles quenched from an equilibrated colloidal fluid. We choose all particles to be of the same size and to interact as hard spheres, thus retaining all features of the porous confinement while limiting the control parameters to the packing fraction of the matrix, φ(m), and that of the fluid, φ(f). We conduct detailed investigations on several dynamic properties, including the tagged-particle and collective intermediate scattering functions, the mean-squared displacement, and the van Hove function. We show the confining obstacles to profoundly impact the relaxation pattern of various quantifiers pertinent to the fluid. Varying the type of quantifier (tagged-particle or collective) as well as φ(m) and φ(f), we unveil both discontinuous and continuous arrest scenarios. Furthermore, we discover subdiffusive behavior and demonstrate its close connection to the matrix structure. Our findings partly confirm the various predictions of a recent extension of mode-coupling theory to the quenched-annealed protocol.

Interplay between demixing and freezing in two-dimensional symmetrical mixtures

The interplay between demixing and freezing in two-dimensional symmetrical binary mixtures of Lennard-Jones particles is studied using Monte Carlo simulation. It is demonstrated that different scenarios are possible. For example, the line of continuous liquid demixing transition can start at the liquid side of the vapor-liquid coexistence at the lower critical end point and then it can terminate at the liquid side of the liquid-demixed solid coexistence at the upper critical end point. Other situations are also possible. We distinguish four different scenarios depending on the interactions between unlike particles.

Colloid-polymer mixtures in the presence of quenched disorder: a theoretical and computer simulation study

We use theory and computer simulation to study the structure and phase behavior of colloid-polymer mixtures in the presence of quenched disorder. The Asakura-Oosawa model (AO) (Asakura and Oosawa 1954 J. Chem. Phys. 22 1255) is used to describe the colloid-colloid, colloid-polymer, and polymer-polymer pair interactions. We then investigate the behavior of this model in the presence of frozen-in (quenched) obstacles. The obstacles will be placed according to two different scenarios, both of which are experimentally feasible. In the first scenario, polymers are distributed at positions drawn from an ideal gas configuration. In the second scenario, colloidal particles are distributed at positions drawn from an equilibrium hard sphere configuration. We investigate how the unmixing transition of the AO model is affected by the type of quenched disorder. The theoretical formalism is based on the replica method of Given and Stell (1994 Physica A 209 495). Our foremost aim is to test the accuracy of three common closures to the replica Ornstein-Zernike equations, namely the hypernetted chain, the Percus-Yevick, and the Martinov-Sarkisov equations. The accuracy is determined by comparison with grand canonical Monte Carlo simulations. We find that, for quenched polymer disorder, all three closures perform remarkably well. However, when quenched colloid disorder is considered, i.e. the second mentioned scenario, the predictions of all three closures worsen dramatically.

Phase separation of model adsorbates in random matrices

The liquid-liquid phase behavior of binary mixtures in random pores is investigated with non-additive hard spheres using both ROZ (Replica Ornstein-Zernike) integral equations and cavity biased grand canonical Monte Carlo simulations. The critical densities of the coexistence phase envelopes are determined as function of the non-additivity parameter Delta, varying from Delta = 0.2, 0.4, 0.6, up to 0.8. The matrix is made of quenched hard spheres. Its porosity is varied to ascertain the effects of confinement, with packing densities rho(0) ranging from 0.1, 0.3, to 0.5. To obtain fiduciary results from ROZ, we use the accurate ZSEP closure relation proposed earlier with and without thermodynamic consistency. The ZSEP closure is known to enforce the zero-separation theorems via special adjustable parameters in the bridge function. Two versions of this closure are used to assess their accuracies (vis-à-vis the Monte Carlo data): first ZSEP-T, namely, the ZSEP closure with added thermodynamic consistency (the Gibbs-Duhem relation); and second purely ZSEP without adding thermodynamic consistency. It is found that both closures give correct qualitative trend, with errors of ZSEP falling within 8-9%, while ZSEP-T, being more accurate, to within 1-2%. As non-additivity is increased, both versions become more accurate. The critical density rho(c) is found to decrease with decreasing porosity. In addition, rho(c) also decreases with increasing Delta, in a non-monotone fashion.

Fluids in porous media. I. A hard sponge model

The morphology of many porous materials is spongelike. Despite the abundance of such materials, simple models which allow for a theoretical description of these materials are still lacking. Here, we propose a hard sponge model which is made by digging spherical cavities in a solid continuum. We found an analytical expression for describing the interaction potential between fluid particles and the spongelike porous matrix. The diagrammatic expansions of different correlation functions are derived as well as that of grand potential. We derived also the Ornstein-Zernike (OZ) equations for this model. In contrast to Madden-Glandt model of random porous media [W. G. Madden and E. D. Glandt, J. Stat. Phys. 51, 537 (1988)], the OZ equations for a fluid confined in our hard sponge model have some similarity to the OZ equations of a three-component fluid mixture. We show also how the replica method can be extended to study our sponge model and that the same OZ equations can be derived also from the extended replica method.

Phase behavior of a symmetrical binary fluid mixture

We have investigated the phase behavior of a symmetrical binary fluid mixture for the situation where the chemical potentials mu(1) and mu(2) of the two species differ. Attention is focused on the set of interparticle interaction strengths for which, when mu(1)=mu(2), the phase diagram exhibits both a liquid-vapor critical point and a tricritical point. The corresponding phase behavior for the case mu(1) not equalmu(2) is investigated via integral-equation theory calculations within the mean spherical approximation and grand canonical Monte Carlo (GCMC) simulations. We find that two possible subtypes of phase behavior can occur, these being distinguished by the relationship between the triple lines in the full phase diagram in the space of temperature, density, and concentration. We present the detailed form of the phase diagram for both subtypes and compare with the results from GCMC simulations, finding good overall agreement. The scenario via which one subtype evolves into the other is also studied, revealing interesting features.

Critical behavior of colloid-polymer mixtures in random porous media

We show that the critical behavior of a colloid-polymer mixture inside a random porous matrix of quenched hard spheres belongs to the universality class of the random-field Ising model. We also demonstrate that random-field effects in colloid-polymer mixtures are surprisingly strong. This makes these systems attractive candidates to study random-field behavior experimentally.

A modified version of a self-consistent Ornstein-Zernike approximation for a fluid with a one-Yukawa pair potential

We present a modified version of a thermodynamically self-consistent Ornstein-Zernike approximation (SCOZA) for a fluid of spherical particles with a pair potential given by a hard core repulsion and a Yukawa tail [Formula: see text]. We take advantage of the known analytical properties of the solution of the Ornstein-Zernike equation for the case in which the direct correlation function outside the repulsive core is given by the multi-screened Coulomb plus power series (multi-SCPPS) tails [Formula: see text] and the radial distribution function g(r) satisfies the exact core condition g(r) = 0 for r<1. The SCOZA is known to provide very good overall thermodynamics and a remarkably accurate critical point and coexistence curve. However, the SCOZA presented so far for continuum fluids has the deficiency that the solution behaves singularly at a density ρ where the screening length z(1)(ρ) of the hard sphere fluid nearly coincides with the Yukawa-tail screening length z(2) (>3.8). This is by no means a rare case in the studies of real fluids and colloidal suspensions. We show that the deficiency is resolved in the modified version of the SCOZA with multi-SCPPS tails. As a demonstration, we present some numerical results for z(2) = 8.0.

Type-IV phase behavior in fluids with an internal degree of freedom

We have identified a fourth archetype of phase diagram in binary symmetrical mixtures, which is encountered when the ratio of the interaction between the unlike and the like particles is sufficiently small. This type of phase diagram is characterized by the fact that the lambda line (i.e., the line of the second-order demixing transition) intersects the first-order liquid-vapor curve at densities smaller than the liquid-vapor critical density.

XY-spin fluids in an external magnetic field: an integral equation approach

We develop an integral equation approach to study anisotropic fluids with planar spins in the presence of an external field. As a result, the integral equation calculations for these systems appear to be no more difficult than those for ordinary isotropic liquids. The method presented is applied to the investigation of phase coexistence properties of ferromagnetic XY-spin fluids in a magnetic field. The soft mean spherical approximation is used for the closure relation connecting the orientationally dependent two-particle direct and total correlation functions. The Lovett-Mou-Buff-Wertheim and Born-Green-Yvon equations are employed to describe the one-particle orientational distribution. The phase diagrams are obtained in the whole range of varying the external field for a wide class of XY-spin fluid models with various ratios of the strengths of magnetic to nonmagnetic Yukawa-like interactions. The influence of changing the screening radii of the interaction potentials is also considered. Different types of the phase diagram topology are identified. They are characterized by the existence of critical, tricritical, critical end, and triple points related to transitions between gas, liquid, and para- and ferromagnetic states, accompanied by different external field dependencies of critical temperatures and densities corresponding to the gas-liquid and liquid-liquid transitions. As is demonstrated, the integral equation approach leads to accurate predictions of the complicated phase diagram behavior which coincide well with those evaluated by the cumbersome Gibbs ensemble simulation and multiple-histogram reweighting techniques.

Phase behavior and local structure of a binary mixture in pores: Mean-field lattice model calculations for analyzing neutron scattering data

We investigate the phase behavior of an asymmetric binary liquid A-W mixture confined between two planar homogenous substrates (slit pore). Molecules of species W interact preferentially with the solid walls via a long-range potential. Assuming nearest-neighbor attractions between the liquid molecules, we employ a lattice-gas model and a mean-field approximation for the grand potential. Minimization of this potential yields the density profiles of thermodynamically stable phases for fixed temperature, chemical potentials of both species, pore width and strengths of attraction. This model is used to analyze experimental small-angle neutron-scattering (SANS) data on the microscopic structure of the binary system isobutyric acid (iBA)+heavy water (D2O) inside a mesoscopic porous matrix (controlled-pore glass of about 10 nm mean pore width). Confinement-independent model parameters are adjusted so that the theoretical liquid-liquid coexistence curve in the bulk matches its experimental counterpart. By choosing appropriate values of the pore width and the attraction strength between substrates and water we analyze the effect of confinement on the phase diagram. In addition to a depression of the liquid-liquid critical point we observe surface induced phase transitions as well as water-film adsorption near the walls. The temperature dependence of the structure of water-rich and iBA-rich phases of constant composition are discussed in detail. The theoretical predictions are consistent with results of the SANS study and assist their interpretation.

Phase diagram of a binary symmetric hard-core Yukawa mixture

We assess the accuracy of the self-consistent Ornstein-Zernike approximation for a binary symmetric hard-core Yukawa mixture by comparison with Monte Carlo simulations of the phase diagrams obtained for different choices of the ratio alpha of the unlike-to-like interactions. In particular, from the results obtained at alpha=0.75 we find evidence for a critical endpoint in contrast to recent studies based on integral equation and hierarchical reference theories. The variation of the phase diagrams with range of the Yukawa potential is investigated.

Vapor-liquid transitions of dipolar fluids in disordered porous media: Performance of angle-averaged potentials

Using replica integral equations in the reference hypernetted-chain (RHNC) approximation we calculate vapor-liquid spinodals, chemical potentials, and compressibilities of fluids with angle-averaged dipolar interactions adsorbed to various disordered porous media. Comparison with previous RHNC results for systems with true angle-dependent Stockmayer (dipolar plus Lennard-Jones) interactions indicate that, for a dilute hard sphere matrix, the angle-averaged fluid-fluid (ff) potential is a reasonable alternative for reduced fluid dipole moments m( *2)=mu(2)/(epsilon(0)sigma(3))< or =2.0. This range is comparable to that estimated in bulk fluids, for which RHNC results are presented as well. Finally, results for weakly polar matrices suggest that angle-averaged fluid-matrix (fm) interactions can reproduce main features observed for true dipolar (fm) interactions such as the shift of the vapor-liquid spinodals towards lower temperatures and higher densities. However, the effective attraction induced by dipolar (fm) interaction is underestimated rather than overestimated as in the case of angle-averaged ff interactions.

Phase behavior of a binary symmetric mixture in slitlike pores with opposing walls: application of density functional approach

We study adsorption of a symmetric binary Lennard-Jones mixture, which exhibits partial mixing in a bulk phase, in slitlike pores formed by the walls having antisymmetric properties with respect to the components. The calculations are carried out by means of a density functional approach. We show that under suitable conditions the pore filling may occur as a sequence of two first-order transitions. The capillary condensation may lead to an "antisymmetric" liquidlike film, the symmetry of which follows the symmetry of the adsorbing potential, or to a "demixed" film, the symmetry of which is only weakly associated with the symmetry of the adsorption potential. The additional first-order antisymmetric-demixed film transition begins at the triple point temperature and ends at the critical end point temperature.

Replica Ornstein-Zernike self-consistent theory for mixtures in random pores

We present a self-consistent integral equation theory for a binary liquid in equilibrium with a disordered medium, based on the formalism of the replica Ornstein-Zernike (ROZ) equations. Specifically, we derive direct formulas for the chemical potentials and the zero-separation theorems (the latter provide a connection between the chemical potentials and the fluid cavity distribution functions). Next we solve a modified-Verlet closure to ROZ equations, which has built-in parameters that can be adjusted to satisfy the zero-separation theorems. The degree of thermodynamic consistency of the theory is also kept under control. We model the binary fluid in random pores as a symmetrical binary mixture of nonadditive hard spheres in a disordered hard-sphere matrix and consider two different values of the nonadditivity parameter and of the quenched matrix packing fraction, at different mixture concentrations. We compare the theoretical structural properties as obtained through the present approach with Percus-Yevick and Martinov-Sarkisov integral equation theories, and assess both structural and thermodynamic properties by performing canonical standard and biased grand canonical Monte Carlo simulations. Our theory appears superior to the other integral equation schemes here examined and provides reliable estimates of the chemical potentials. This feature should be useful in studying the fluid phase behavior of model adsorbates in random pores in general.

Ising fluids in an external magnetic field: an integral equation approach

The phase behavior of Ising spin fluids is studied in the presence of an external magnetic field with the integral equation method. The calculations are performed on the basis of a soft mean spherical approximation using an efficient algorithm for solving the coupled set of the Ornstein-Zernike equations, the closure relations, and the external field constraint. The phase diagrams are obtained in the whole thermodynamic space including the magnetic field H for a wide class of Ising fluid models with various ratios R of the strengths of magnetic to nonmagnetic Yukawa-like interactions. The influence of varying the inverse screening lengths z(1) and z(2), corresponding to the magnetic and nonmagnetic Yukawa parts of the potential, is investigated too. It is shown that changes in R as well as in z(1) and z(2) can lead to different topologies of the phase diagrams. In particular, depending on the value of R, the critical temperature of the liquid-gas transition either decreases monotonically, behaves nonmonotonically, or increases monotonically with increasing H. The para-ferro magnetic transition is also affected by changes in R and the screening lengths. At H=0, the Ising fluid maps onto a simple model of a symmetric nonmagnetic binary mixture. For H--> infinity, it reduces to a pure nonmagnetic fluid. The results are compared with available simulations and the predictions of other theoretical methods. It is demonstrated that the mean spherical approximation appears to be more accurate compared with mean field theory, especially for systems with short ranged attraction potentials (when z(1) and z(2) are large). In the Kac limit z(1), z(2) -->+0, both approaches tend to nearly the same results.

Integral equation study of a Stockmayer fluid adsorbed in polar disordered matrices

Based on replica integral equations in the (reference) hypernetted chain approximation we investigate the structural features and phase properties of a dipolar Stockmayer fluid confined to a disordered dipolar matrix. The integral equations are applied to the homogeneous high-temperature phase where the system is globally isotropic. At low densities we find the influence of dipolar interactions between fluid (f) and matrix (m) particles to be surprisingly similar to the previously investigated effect of attractive isotropic (fm) interactions: the critical temperature of the vapor-liquid transition decreases with increasing (fm) coupling, while the critical density increases. The anisotropic nature of the dipolar (fm) interactions turns out to play a more dominant role at high fluid densities where we observe a pronounced sensitivity in the dielectric constant and a strong degree of local orientational ordering of the fluid particles along the local fields generated by the matrix. Moreover, an instability of the dielectric constant, which is a precursor of ferroelectric ordering occurring both in bulk Stockmayer fluids and in fluids in nonpolar matrices, is observed only for very small dipolar (fm) couplings.

Continuous demixing at liquid-vapor coexistence in a symmetrical binary fluid mixture

We report a Monte Carlo finite-size scaling study of the demixing transition of a symmetrical Lennard-Jones binary fluid mixture. For equal concentration of species, and for a choice of the unlike-to-like interaction ratio delta=0.7, this transition is found to be continuous at liquid-vapor coexistence. The associated critical end point exhibits an Ising-like universality. These findings confirm those of earlier smaller scale simulation studies of the same model, but contradict the findings of recent integral equation and hierarchical reference theory investigations.

Phase diagram of symmetric binary mixtures at equimolar and nonequimolar concentrations: a systematic investigation

We consider symmetric binary mixtures consisting of spherical particles with equal diameters interacting via a hard-core plus attractive tail potential with strengths epsilon(ij), i,j=1,2, such that epsilon(11)=epsilon(22)>epsilon(12). The phase diagram of the system at all densities and concentrations is investigated as a function of the unlike-to-like interaction ratio delta=epsilon(12)/epsilon(11) by means of the hierarchical reference theory. The results are related to those of previous investigations performed at equimolar concentration, as well as to the topology of the mean-field critical lines. As delta is increased in the interval 0<delta<1, we find first a regime where the phase diagram at equal species concentration displays a tricritical point, then one where both a tricritical and a liquid-vapor critical point are present. We did not find any clear evidence of the critical end point topology predicted by mean-field theory as delta approaches 1, at least up to delta=0.8, which is the largest value of delta investigated here. Particular attention was paid to the description of the critical-plus-tricritical point regime in the whole density-concentration plane. In this situation, the phase diagram shows, in a certain temperature interval, a coexistence region that encloses an island of homogeneous, one-phase fluid.

Capillary condensation of a binary mixture in slit-like pores

We investigate the capillary condensation of two model fluid mixtures in slit-like pores, which exhibit different demixing properties in the bulk phase. The interactions between adsorbate particles are modeled by using Lennard-Jones (12,6) potentials and the adsorbing potentials are of the Lennard-Jones (9,3) type. The calculations are performed for different pore widths and at different concentrations of the bulk gas, by means of density functional theory. We evaluate the capillary phase diagrams and discuss their dependence on the parameters of the model. Our calculations indicate that a binary mixture confined to a slit-like pore may exhibit rich phase behavior.

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.

Binary mixtures of magnetic fluids

We study a binary mixture of a van der Waals fluid and a ferromagnetic fluid at zero magnetic field on the basis of the mean field Ising fluid model and the van der Waals theory with quadratic mixing rules. Depending on three reduced parameters, the phase diagram shows a surface of magnetic phase transitions and lines of tricritical points, critical end points, and magnetic consolute points. First-order phase transition surfaces and critical lines are calculated numerically. For the line of tricritical points, which can occur in two different topologies, an analytic expression is derived. All higher-order lines and coexistence surfaces are visualized in three-dimensional x, T, p and xi, T, p diagrams, where xi is a mapping of delta, the conjugated field of the mole fraction x, on the unit interval.