Critical-point and coexistence curve properties of a symmetric mixture of nonadditive hard spheres: A finite size scaling study

Investigation of Fluid-Fluid and Solid-Solid Phase Separation of Symmetric Nonadditive Hard Spheres at High Density

We have calculated the values of the critical packing fractions for the mixtures of symmetric nonadditive hard spheres at high densities for small values of the nonadditivity parameter. Calculations have been performed for solid-solid and fluid-fluid demixing transitions. A cluster algorithm for Monte Carlo simulations in a semigrand ensemble was used, and the waste recycling method was applied to improve the accuracy of the calculations. The finite size scaling analysis was employed to compute the critical packing fractions for infinite systems with high accuracy.

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).

Monte Carlo simulation of the nonadditive restricted primitive model of ionic fluids: phase diagram and clustering

We report an accurate Monte Carlo calculation of the phase diagram and clustering properties of the restricted primitive model with nonadditive hard-sphere diameters. At high density the positively nonadditive fluid shows more clustering than in the additive model and the negatively nonadditive fluid shows less clustering than in the additive model; at low density the reverse scenario appears. A negative nonadditivity tends to favor the formation of neutrally charged clusters starting from the dipole. A positive nonadditivity favors the pairing of like ions at high density. The critical point of the gas-liquid phase transition moves at higher temperatures and higher densities for a negative nonadditivity and at lower temperatures and lower densities for a positive nonadditivity. The law of corresponding states does not seem to hold strictly. Our results can be used to interpret recent experimental works on room temperature ionic liquids.

Nonadditive hard-sphere fluid mixtures: a simple analytical theory

We construct a nonperturbative fully analytical approximation for the thermodynamics and the structure of nonadditive hard-sphere fluid mixtures. The method essentially lies in a heuristic extension of the Percus-Yevick solution for additive hard spheres. Extensive comparison with Monte Carlo simulation data shows a generally good agreement, especially in the case of like-like radial distribution functions.

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.

Binary non-additive hard sphere mixtures: fluid demixing, asymptotic decay of correlations and free fluid interfaces

Using a fundamental measure density functional theory we investigate both bulk and inhomogeneous systems of the binary non-additive hard sphere model. For sufficiently large (positive) non-additivity the mixture phase separates into two fluid phases with different compositions. We calculate bulk fluid-fluid coexistence curves for a range of size ratios and non-additivity parameters and find that they compare well to simulation results from the literature. Using the Ornstein-Zernike equation, we investigate the asymptotic, [Formula: see text], decay of the partial pair correlation functions, g(ij)(r). At low densities a structural crossover occurs in the asymptotic decay between two different damped oscillatory modes with different wavelengths corresponding to the two intra-species hard-core diameters. On approaching the fluid-fluid critical point there is a Fisher-Widom crossover from exponentially damped oscillatory to monotonic asymptotic decay. Using the density functional we calculate the density profiles for the planar free fluid-fluid interface between coexisting fluid phases. We show that the type of asymptotic decay of g(ij)(r) not only determines the asymptotic decay of the interface profiles, but is also relevant for intermediate and even short-ranged behaviour. We also determine the surface tension of the free fluid interface, finding that it increases with non-additivity, and that on approaching the critical point mean-field scaling holds.

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.

Surface tension of the Widom-Rowlinson model

We consider the computation of the surface tension of the fluid-fluid interface for the Widom-Rowlinson [J. Chem. Phys. 52, 1670 (1970)] binary mixture from direct simulation of the inhomogeneous system. We make use of the standard mechanical route, in which the surface tension follows from the computation of the normal and tangential components of the pressure tensor of the system. In addition to the usual approach, which involves simulations of the inhomogeneous system in the canonical ensemble, we also consider the computation of the surface tension in an ensemble where the pressure perpendicular (normal) to the planar interface is kept fixed. Both approaches are seen to provide consistent values of the interfacial tension. The issue of the system-size dependence of the surface tension is addressed. In addition, simulations of the fluid-fluid coexistence properties of the mixture are performed in the semigrand canonical ensemble. Our results are compared with existing data of the Widom-Rowlinson mixture and are also examined in the light of the vapor-liquid equilibrium of the thermodynamically equivalent one-component penetrable sphere model.

Fundamental measure density functional theory for nonadditive hard-core mixtures: the one-dimensional case

We treat a one-dimensional binary mixture of hard-core particles that possess nonadditive diameters. For this model, a density functional theory is constructed following similar principles as an earlier extension of Rosenfeld's fundamental measure theory to three-dimensional nonadditive hard-sphere mixtures. The theory applies to arbitrary positive and moderate negative nonadditivity and reduces to Percus' exact functional in the additive case. Bulk direct correlation functions are obtained as functional derivatives of the excess free energy functional. Results for the partial pair correlation functions in bulk, as calculated via the Ornstein-Zernike route and using the direct correlation functions as input, show very good agreement with results from our Monte Carlo computer simulations of the mixture.

Virial Coefficients and Demixing of Athermal Nonadditive Mixtures

We compute the fourth virial coefficient of a binary nonadditive, hard-sphere mixture over a wide range of deviations from diameter additivity and size ratios. Hinging on this knowledge, we build up a y expansion (Barboy, B.; Gelbart, W. N. J. Chem. Phys. 1979, 71, 3053) in order to trace the fluid-fluid coexistence lines, which we then compare with the available Gibbs-ensemble Monte Carlo data and with the estimates obtained through two refined integral-equation theories of the fluid state. We find that in a regime of moderately negative nonadditivity and largely asymmetric diameters, relevant to the modeling of sterically and electrostatically stabilized colloidal mixtures, the fluid-fluid critical point is unstable with respect to crystallization.

Structure and phase behavior of Widom-Rowlinson model calculated from a nonuniform Ornstein-Zernike equation

The second-order integral-equation formalism of [Attard J. Chem. Phys. 91, 3072 (1989); 95, 4471 (1991)], applied previously to one-component hard spheres and Lennard-Jones fluids, as well as to their mixtures, is used to binary Widom-Rowlinson mixtures. Comparison with Monte Carlo simulations of the pair correlation functions and of the demixing phase diagram shows that this method is also quite accurate in the case of highly nonadditive mixtures. Moreover, the results of the second-order theory are compared with previous theoretical predictions. Our interest is also in the calculation of the bridge functions, i.e., parts of the radial distribution functions either not included or simply approximated in the usual theories.

Isotropic-nematic interfacial tension of hard and soft rods: Application of advanced grand canonical biased-sampling techniques

Coexistence between the isotropic and the nematic phase in suspensions of rods is studied using grand canonical Monte Carlo simulations with a bias on the nematic order parameter. The biasing scheme makes it possible to estimate the interfacial tension gamma(IN) in systems of hard and soft rods. For hard rods with LD=15, we obtain gammaIN approximately 1.4kBT/L2, with L the rod length, D the rod diameter, T the temperature, and kB the Boltzmann constant. This estimate is in good agreement with theoretical predictions, and the order of magnitude is consistent with experiments.

On the separation of nonadditive symmetric mixtures in nanoscopic slitlike pores: A simple model for racemic fluids

A grand canonical ensemble Monte Carlo simulation method is used to study the adsorption of nonadditive symmetric mixtures of Lennard-Jones spherical particles in nanoscopic slitlike pores. The walls of the pore are assumed to be formed by the parallel (100) planes of the model face centered cubic crystal of adjustable corrugation potential. It is demonstrated that depending on the nonadditivity effects in the mixture and the pore width the condensed phases formed inside the pore may have different structures. In particular, it is shown that the mixture may separate into layers containing only one component each and the stacking may depend on the pore width and properties of the mixture.

Interfacial tension of the isotropic-nematic interface in suspensions of soft spherocylinders

The isotropic to nematic transition in a system of soft spherocylinders is studied by means of grand canonical Monte Carlo simulations. The probability distribution of the particle density is used to determine the coexistence densities of the isotropic and the nematic phases. The distributions are also used to compute the interfacial tension of the isotropic-nematic interface, including an analysis of finite size effects. Our results confirm that the Onsager limit is not recovered until for very large elongation, exceeding at least L/D=40 , with L the spherocylinder length and D the diameter. For smaller elongation, we find that the interfacial tension increases with increasing L/D , in agreement with theoretical predictions.

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.

Equation of state of nonadditive d-dimensional hard-sphere mixtures

An equation of state for a multicomponent mixture of nonadditive hard spheres in d dimensions is proposed. It yields a rather simple density dependence and constitutes a natural extension of the equation of state for additive hard spheres proposed by us [A. Santos, S. B. Yuste, and M. Lopez de Haro, Mol. Phys. 96, 1 (1999)]. The proposal relies on the known exact second and third virial coefficients and requires as input the compressibility factor of the one-component system. A comparison is carried out both with another recent theoretical proposal based on a similar philosophy and with the available exact results and simulation data in d=1, 2, and 3. Good general agreement with the reported values of the virial coefficients and of the compressibility factor of binary mixtures is observed, especially for high asymmetries and/or positive nonadditivities.

Cluster algorithm for nonadditive hard-core mixtures

In this paper, we present a cluster algorithm for the numerical simulations of nonadditive hard-core mixtures. This algorithm allows one to simulate and equilibrate systems with a number of particles two orders of magnitude larger than previous simulations. The phase separation for symmetric binary mixtures is studied for different nonadditivities as well as for the Widom-Rowlinson model [B. Widom and J. S. Rowlinson, J. Chem. Phys. 52, 1670 (1970)] in two and three dimensions. The critical densities are determined from finite size scaling. The critical exponents for all the nonadditivities are consistent with the Ising universality class.

Integral Equation Theory for Symmetric Nonadditive Hard Sphere Mixtures

An integral equation theory is presented for the pair correlation functions and phase behavior of symmetric nonadditive hard sphere mixtures with hard sphere diameters given by sigma(A)(A)() = sigma(BB) = lambdad and sigma(AB) = d. This mixture exhibits a fluid-fluid phase separation into an A-rich phase and a B-rich phase at high densities. The theory incorporates, into the closure approximation, all terms that can be calculated exactly in the density expansion of the direct correlation functions. We find that the closure approximation developed in this work is accurate for the structure and phase behavior over the entire range of lambda, when compared to computer simulations, and is significantly more accurate than the previous theories.

Grand canonical Monte Carlo simulation of a model colloid–polymer mixture: Coexistence line, critical behavior, and interfacial tension

Grand canonical Monte Carlo simulations are used to study phase separation in a simple colloid-polymer model, the so-called Asakura-Oosawa model. To overcome the problem of small acceptance rates of the grand-canonical moves, cluster moves are introduced. Successive umbrella sampling, recently introduced by Virnau and Muller [J. Chem. Phys. 120, 10925 (2004)], is used to access the phase-separated regime. The unmixing binodal and the interfacial tension are measured and compared to theoretical predictions. By means of finite-size scaling, the behavior close to the critical point is also investigated. Close to criticality, we observe substantial deviations from mean-field behavior.