A globally accurate theory for a class of binary mixture models (original) (raw)
2001
https://doi.org/10.1080/00268970210124783
Sign up for access to the world's latest research
checkGet notified about relevant papers
checkSave papers to use in your research
checkJoin the discussion with peers
checkTrack your impact
Abstract
Using the self-consistent Ornstein-Zernike approximation (SCOZA) results for the 3D Ising model, we obtain phase diagrams for binary mixtures described by decorated models. We obtain the plait point, binodals, and closed-loop coexistence curves for the models proposed by Widom, Clark, Neece, and Wheeler. The results are in good agreement with series expansions and experiments.
Related papers
Equilibrium cluster distributions of the three-dimensional Ising model in the one phase region
1983
Abstract We analyse equilibrium cluster distributions obtained numerically from a ferromagnetic Ising model (simple cubic lattice, 125000 sites and periodic boundary conditions) along the coexistence line and in the one-phase region below T c. We find evidences that the distribution of sizes and energies scales with temperature and external magnetic field giving Binder's droplet exponent y≈ 4/9.
The phase behavior of two-dimensional symmetrical mixtures
The Journal of Chemical Physics, 2010
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 subjected to a weakly corrugated external field of a square symmetry. It is shown that the both vapor-liquid condensation and demixing transition in the liquid phase are not appreciably affected by a weakly corrugated external field. On the other hand, even a weakly corrugated external field considerably influences the structure of solid phases and the liquid-solid transition. In particular, the solid phases are found to exhibit uniaxially ordered distorted hexagonal structure. The triple point temperature increases with the corrugation of the external field, while the triple point density becomes lower when the surface corrugation increases. The changes in the location of the triple point are shown to lead to the changes of the phase diagram topology. It is also demonstrated that the solid phase undergoes a demixing transition, which is also very slightly affected by the weakly corrugated external potential. The demixing transition in the solid phase is shown to belong to the universality class of the Ising model.
Lattice Boltzmann model for binary mixtures
Physical Review E, 2002
An a priori derivation of the lattice Boltzmann equations for binary mixtures is provided by discretizing the Boltzmann equations that govern the evolution of binary mixtures. The present model leads to a set of two-fluid hydrodynamic equations for the mixture. In existing models, employing the single-relaxation-time approximation, the viscosity and diffusion coefficients are coupled through the relaxation parameter , thus limited to unity Prandtl number and Schmidt number. In the present model the viscosity and diffusion coefficient are independently controlled by two relaxation parameters, thus enabling the modeling of mixtures with an arbitrary Schmidt number. The theoretical framework developed here can be readily applied to multiplespecies mixing.
Zwanzig model of multi-component mixtures of biaxial particles:y3theory re-visited
Molecular Physics, 2006
The paper considers the thermodynamic and phase ordering properties of a multi-component Zwanzig mixture of hard rectangular biaxial parallelepipeds. An equation of state (EOS) is derived based on an estimate of the number of arrangements of the particles on a threedimensional cubic lattice. The methodology is a generalization of the Flory-DiMarzio counting scheme, but, unlike previous work, this treatment is thermodynamically consistent. The results are independent of the order in which particles are placed on the lattice. By taking the limit of zero lattice spacing, a translationally continuous variant of the model (the off-lattice variant) is obtained. The EOS is identical to that obtained previously by a wide variety of different approaches. In the off-lattice limit, it corresponds to a third-level y-expansion and, in the case of a binary mixture of square platelets, it also corresponds to the EOS obtained from fundamental measure theory. On the lattice it is identical to the EOS obtained by retaining only complete stars in the virial expansion. The off-lattice theory is used to study binary mixtures of rods (R 1 À R 2) and binary mixtures of platelets (P 1 À P 2). The particles were uniaxial, of length (thickness) L and width D. The aspect ratios À i ¼ L i /D i of the components were kept constant (À 1R ¼ 15, À 1P ¼ 1/15 and À 2R ¼ 150, À 2P ¼ 1/150), so the second virial coefficient of R 1 was identical to P 1 and similarly for R 2 and P 2. The volume ratio of particles 1 and 2, v 1 /v 2 , was then varied, with the constraints that v iR ¼ v iP and v 2R ¼ 150D 3 2R : Results on nematic-isotropic (N À I) phase coexistence at an infinite dilution of component 2, are qualitatively similar for rods and platelets. At small values of the ratio v 1 /v 2 , the addition of component 2 (i.e. a thin rod (e.g. a polymer) or a thin plate) results in the stabilization of the nematic phase. For larger values of v 1 /v 2 , however, this effect is reversed and the addition of component 2 destabilizes the nematic. For similar molecular volumes of the two components strong fractionation is observed: shorter rods and thicker platelets congregate in the isotropic phase. In general, the stabilization of the ordered phase and the fractionation between the phases are both weaker in the platelet mixtures. The calculated spinodal curves for isotropic-isotropic demixing are noticeably different between the R 1 À R 2 and the P 1 À P 2 systems. The platelet mixtures turn out to be stable with respect to de-mixing up to extremely high densities. The values of the consolute points for the R 1 À R 2 blends are remarkably similar to those obtained using the Parsons-Lee approximation for bi-disperse mixtures of freely rotating cylinders with similar aspect ratios [S. Varga. A. Galindo, G. Jackson, Mol. Phys., 101, 817 (2003)]. In a number of R 1 À R 2 mixtures, phase diagrams exhibiting both N À I equilibrium and I À I de-mixing were calculated. The latter is preempted by nematic ordering in all the cases studied. Calculations show the possible appearance of azeotropes in the N À I coexistence domain.
Effects oftheExchange Interaction Parameter onthePhase Diagrams fortheTransverse Ising Model
The microscopic effect of the exchange interaction parameter for the 2-Dimensional ising model with nearest neighbor interaction has been studied. By supposing simple temperature dependent relationship for the exchange parameter, graphs were straightforwardly obtained that show the reentrant closed looped phase diagrams symptomatic of some colloids and complex fluids and some binary liquids mixtures in particular. By parameter modifications, other phase diagrams were also obtained. Amongst which are the u-shapes and other exotic shapes of phase diagrams. Our results show that the exchange interaction parameter greatly influence the size of the ordered phase. Hence the larger the value of the constant, the larger the size of the ordered phase. This means that the higher values of the exchange parameter brings about phase transitions that straddle a wider range of polarizations and temperatures.
Theory and Simulation of the Ising Model
2021
We have provided a concise introduction to the Ising model as one of the most important models in statistical mechanics and in studying the phenomenon of phase transition. The required theoretical background and derivation of the Hamiltonian of the model have also been presented. We finally have discussed the computational method and details to numerically solve the twoand threedimensional Ising problems using Monte Carlo simulations. The related computer codes in both Python and Fortran, as well as a simulation trick to visualize the spin lattice, have also been provided.
Droplets in the coexistence region of the two-dimensional Ising model
Journal of Physics A: Mathematical and General, 2000
The two-dimensional Ising model with fixed magnetization is studied using Monte Carlo techniques. At the coexistence line, the macroscopic, extensive droplet of minority spins becomes thermally unstable by breaking up into microscopic clusters. Intriguing finite-size effects as well as singularities of thermal and cluster properties associated with the transition are discussed.
Projecting Ising Model Parameters for Fast Mixing
2014
Inference in general Ising models is difficult, due to high treewidth making treebased algorithms intractable. Moreover, when interactions are strong, Gibbs sampling may take exponential time to converge to the stationary distribution. We present an algorithm to project Ising model parameters onto a parameter set that is guaranteed to be fast mixing, under several divergences. We find that Gibbs sampling using the projected parameters is more accurate than with the original parameters when interaction strengths are strong and when limited time is available for sampling.
Phase behaviour of very asymmetric binary mixtures
Journal of Physics: Condensed Matter, 2000
The phase behaviour of very asymmetric binary mixtures can be understood in terms of the depletion interaction. For hard particles this yields a narrow deep attractive well surrounding the hard core. Colloids with similar interaction potentials are known to destabilize the liquid, causing it to show a wide fluid-solid coexistence, and in extreme cases they exhibit an exotic solid-solid condensation. For a mixture this means that phase separation is not fluid-fluid, as previously thought, but normally fluid-solid, and if the asymmetry is very large, even solid-solid. We present in this work the result of devising a density functional theory for an infinitely asymmetric mixture of parallel hard cubes. This model is singular and undergoes a collapse in a close-packed solid (an extreme fluid-solid demixing). We avoid this collapse by introducing a small amount of polydispersity in the large particles; the resulting phase diagram shows the fluid-solid and solid-solid demixing scenarios described above. † Current address:
Loading Preview
Sorry, preview is currently unavailable. You can download the paper by clicking the button above.
References (15)
- WHEELER, J. C., 1977, Ann. Rev. Phys. Chem. 28, 411.
- HØYE, J. S., and STELL, G., 1984, Mol. Phys. 52, 1071; 1985, Int. J. Thermophys. 6, 561; 1977, J. Chem. Phys. 67, 439.
- DICKMAN, R., and STELL, G., 1996, Phys. Rev. Lett. 77, 966 (1996);
- DICKMAN, R., STELL, G., and PINI, D., 1998, Phys. Rev. E. 57, 2862.
- WIDOM, B., 1967, J. Chem. Phys. 46, 3324.
- CLARK, R. K., 1968, J. Chem. Phys. 48, 741.
- NEECE, G., 1967, J. Chem. Phys. 47, 4112.
- WHEELER, J. C., 1972, Ber. Bunsenges. Phys. Chem. 76, 308.
- WHEELER, J. C., 1975, J. Chem. Phys. 62, 433.
- HØYE, J. S., PINI, D., and STELL, G., 2000, Physica A 279, 213.
- ZINN-JUSTIN, J., 2001, Phys. Rep. 344, 159.
- SCESNEY, P. E., 1970, Phys. Rev. B 1, 2274.
- McEWAN, B. C., 1923, J. Chem. Soc. (London) 123, 2284.
- PARVATIKER, R. R., and McEWAN, B. C., 1924, J. Chem. Soc. (London) 125, 1484.
- FISHER, M. E., 1968, Phys. Rev. 176, 257.
Related papers
Critical exponents in the Lin—Taylor model of asymmetrical associating binary mixtures
Molecular Physics, 1998
The critical behaviour of an associating¯uid mixture, the so called Lin± Taylor model, has been studied, with focus on the critical exponents which characterize di erent points on the critical manifold. In general the critical behaviour of the mixture is found to be Ising-like, but there are some special situations: when two ordinary critical points merge into a double critical point or three critical points coalesce at a critical in¯ection point. A check has been made on the validity of the relations between critical exponents and predictions from the scaling laws.
Computing the phase diagram of binary mixtures: A patchy particle case study
The Journal of Chemical Physics, 2013
We investigate the phase behaviour of 2D mixtures of bi-functional and three-functional patchy particles and 3D mixtures of bi-functional and tetra-functional patchy particles by means of Monte Carlo simulations and Wertheim theory. We start by computing the critical points of the pure systems and then we investigate how the critical parameters change upon lowering the temperature. We extend the Successive Umbrella Sampling method to mixtures to make it possible to extract information about the phase behaviour of the system at a fixed temperature for the whole range of densities and compositions of interest.
The Journal of Chemical Physics, 2002
In this work we present structure factors and triplet direct correlation functions extracted from extensive Monte Carlo simulations for a binary mixture of hard spheres. The results are compared with the predictions of two integral equation theories, namely, a recently proposed extension to mixtures of Attard's inhomogeneous integral equation approach, and Barrat, Hansen, and Pastore's factorization ansatz. In general, both theories yield acceptable estimates for the triplet structure functions, though, by construction, the inhomogeneous integral equation theory is more suited to furnish triplet distribution function results, whereas the factorization ansatz provides a more handy approach to triplet direct correlation functions.
Dynamics of discrete models of binary mixtures in two dimensions: Exact solution
Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics, 1994
The binary mixture of interacting disparate particles (heavy and light) in two dimensions has been modeled with the help of discrete-velocity Boltzmann-Broadwell models. These have been solved analytically for all integer and half-integer values of coupling constants. Depending on initial conditions the one-particle distribution functions and the entropy of light particles may exhibit a nonmonotonic behavior, as a function of time.
Physical Review E, 2021
A cluster weight Ising model is proposed by introducing an additional cluster weight in the partition function of the traditional Ising model. It is equivalent to the O(n) loop model or ncomponent face cubic loop model on the two-dimensional lattice, but on the three-dimensional lattice, it is still not very clear whether or not these models have the same universality. In order to simulate the cluster weight Ising model and search for new universality class, we apply a cluster algorithm, by combining the color-assignation and the Swendsen-Wang methods. The dynamical exponent for the absolute magnetization is estimated to be z = 0.45(3) at n = 1.5, consistent with that by the traditional Swendsen-Wang methods. The numerical estimation of the thermal exponent yt and magnetic exponent ym, show that the universalities of the two models on the three dimensional lattice are different. We obtain the global phase diagram containing paramagnetic and ferromagnetic phases. The phase transition between the two phases are second order at 1 ≤ n < nc and first order at n ≥ nc, where nc ≈ 2. The scaling dimension yt equals to the system dimension d when the first order transition occurs. Our results are helpful in the understanding of some traditional statistical mechanics models.
Entropy
Most undergraduate students who have followed a thermodynamics course would have been asked to evaluate the volume occupied by one mole of air under standard conditions of pressure and temperature. However, what is this task exactly referring to? If air is to be regarded as a mixture, under what circumstances can this mixture be considered as comprising only one component called “air” in classical statistical mechanics? Furthermore, following the paradigmatic Gibbs’ mixing thought experiment, if one mixes air from a container with air from another container, all other things being equal, should there be a change in entropy? The present paper addresses these questions by developing a prior-based statistical mechanics framework to characterise binary mixtures’ composition realisations and their effect on thermodynamic free energies and entropies. It is found that (a) there exist circumstances for which an ideal binary mixture is thermodynamically equivalent to a single component ideal...
Molecular Physics, 2001
A grand canonical Monte Carlo (GCMC) simulation method is presented for the determination of the phase equilibria of mixtures. The coexistence is derived by expanding the pressure into a Taylor series as a function of the temperature and the chemical potentials that are the independent intensive variables of the grand canonical ensemble. The coe cients of the Taylor series can be calculated from ensemble averages and¯uctuation formulae that are obtained from GCMC simulations in both phases. The method is able to produce the equilibrium data in a certain domain of the …T ; p † plane from two GCMC simulations. The vapour±liquid equilibrium results obtained for a Lennard-Jones mixture agree well with the corresponding Gibbs ensemble Monte Carlo data.
Binary mixtures in and beyond the critical region: thermodynamic properties
Fluid Phase Equilibria, 1998
An extended law of corresponding states, in combination with a parametric crossover equation of state, is applied to the prediction of thermodynamic properties and phase behavior of pure fluids and binary mixtures over a wide region around vapor-liquid critical points. For pure fluids, a generalized form of the model is given that requires only the acentric factor of the pure component. For mixtures, the critical locus is also required. A Ž . good representation of thermodynamic property data is achieved in the range of temperatures 0.8 T x F T F 1.5 c Ž .
Related topics
Cited by
Chemical Physics, 2005
The vapor-liquid phase behavior and the critical behavior of the square-well (SW) fluid are investigated as a function of the interaction range, λ∈ [1.25, 3], by means of the self-consistent Ornstein-Zernike approximation (SCOZA) and analytical equations of state based on a perturbation theory [A. L. Benavides and F. del Rio, Mol. Phys. 68, 983 (1989); A. Gil-Villegas, F. del Rio, and A. L. Benavides, Fluid Phase Equilib. 119, 97 (1996)]. For this purpose the SCOZA, which has been restricted up to now to a few model systems, has been generalized to hard-core systems with arbitrary interaction potentials requiring a fully numerical solution of an integro-partial differential equation. Both approaches, in general, describe well the liquid-vapor phase diagram of the square-well fluid when compared with simulation data. SCOZA yields very precise predictions for the coexistence curves in the case of long ranged SW interaction (λ>1.5), and the perturbation theory is able to predict the binodal curves and the saturated pressures, for all interaction ranges considered if one stays away from the critical region. In all cases, the SCOZA gives very good predictions for the critical temperatures and the critical pressures, while the perturbation theory approach tends to slightly overestimate these quantities. Furthermore, we propose analytical expressions for the critical temperatures and pressures as a function of the square-well range.