The hierarchical reference theory as applied to square well fluids of variable range

Microcanonical ensemble simulation method applied to discrete potential fluids

In this work we extend the applicability of the microcanonical ensemble simulation method, originally proposed to study the Ising model [A. Hüller and M. Pleimling, Int. J. Mod. Phys. C 13, 947 (2002)0129-183110.1142/S0129183102003693], to the case of simple fluids. An algorithm is developed by measuring the transition rates probabilities between macroscopic states, that has as advantage with respect to conventional Monte Carlo NVT (MC-NVT) simulations that a continuous range of temperatures are covered in a single run. For a given density, this new algorithm provides the inverse temperature, that can be parametrized as a function of the internal energy, and the isochoric heat capacity is then evaluated through a numerical derivative. As an illustrative example we consider a fluid composed of particles interacting via a square-well (SW) pair potential of variable range. Equilibrium internal energies and isochoric heat capacities are obtained with very high accuracy compared with data obtained from MC-NVT simulations. These results are important in the context of the application of the Hüller-Pleimling method to discrete-potential systems, that are based on a generalization of the SW and square-shoulder fluids properties.

Improved renormalization group theory for critical asymmetry of fluids

We develop an improved renormalization group (RG) approach incorporating the critical vapor-liquid equilibrium asymmetry. In order to treat the critical asymmetry of vapor-liquid equilibrium, the integral measure is introduced in the Landau-Ginzbug partition function to achieve a crossover between the local order parameter in Ising model and the density of fluid systems. In the implementation of the improved RG approach, we relate the integral measure with the inhomogeneous density distribution of a fluid system and combine the developed method with SAFT-VR (statistical associating fluid theory of variable range) equation of state. The method is applied to various fluid systems including square-well fluid, square-well dimer fluid and real fluids such as methane (CH4), ethane (C2H6), trifluorotrichloroethane (C2F3Cl3), and sulfur hexafluoride (SF6). The descriptions of vapor-liquid equilibria provided by the developed method are in excellent agreement with simulation and experimental data. Furthermore, the improved method predicts accurate and qualitatively correct behavior of coexistence diameter near the critical point and produces the non-classical 3D Ising criticality.

Critical asymmetry in renormalization group theory for fluids

The renormalization-group (RG) approaches for fluids are employed to investigate critical asymmetry of vapour-liquid equilibrium (VLE) of fluids. Three different approaches based on RG theory for fluids are reviewed and compared. RG approaches are applied to various fluid systems: hard-core square-well fluids of variable ranges, hard-core Yukawa fluids, and square-well dimer fluids and modelling VLE of n-alkane molecules. Phase diagrams of simple model fluids and alkanes described by RG approaches are analyzed to assess the capability of describing the VLE critical asymmetry which is suggested in complete scaling theory. Results of thermodynamic properties obtained by RG theory for fluids agree with the simulation and experimental data. Coexistence diameters, which are smaller than the critical densities, are found in the RG descriptions of critical asymmetries of several fluids. Our calculation and analysis show that the approach coupling local free energy with White's RG iteration which aims to incorporate density fluctuations into free energy is not adequate for VLE critical asymmetry due to the inadequate order parameter and the local free energy functional used in the partition function.

Reformulation of liquid perturbation theory for low temperatures

A perturbation strategy is proposed which allows liquid perturbation methodology to be applied to extremely low temperature situations for which the available liquid integral equation theory and traditional thermodynamic perturbation theory fail. The possibility of avoiding the low temperature problem of the thermodynamic perturbation theory not only is of relevance to the investigation of complex fluids, but also may be useful for reformulation of other liquid theories to achieve higher accuracy and avoid the respective low temperature problems.

Liquid-vapor transition of systems with mean field universality class

We have considered a system where the interaction, v(r)=v(IS)(r)+xi(2)v(MF)(r), is given as a linear combination of two potentials, each of which being characterized with a well-defined critical behavior: for v(IS)(r) we have chosen the potential of the restricted primitive model which is known to belong to the three-dimensional Ising universality class, while for v(MF)(r) we have considered a long-range interaction in the Kac [J. Math. Phys. 4, 216 (1963)] limit, displaying mean field (MF) behavior. We study the performance of two theoretical approaches and of computer simulations in the critical region for this particular system and give a detailed comparison between theories and simulation of the critical region and the location of the critical point. Having shown by theoretical arguments that the system belongs to the MF universality class for any positive value of xi and shows nonclassical behavior only for xi=0, we examine to which extent theoretical approximations and simulation can reproduce this behavior. While in this limiting case theoretical approaches are known to fail, we find good agreement for the critical properties between the theoretical approaches and the simulations for xi(2) larger than 0.05.

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.

Thermodynamic properties of short-range square well fluid

The interfacial properties of short-range square well fluid with lambda=1.15, 1.25, and 1.375 were determined by using single canonical Monte Carlo simulations. Simulations were carried out in the vapor-liquid region. The coexistence curves of these models were calculated and compared to those previously reported in the literature and good agreement was found among them. We found that the surface tension curves for any potential model of short range form a single master curve when we plot gamma* vs TT(c). It is demonstrated that the critical reduced second virial coefficient B(2)* as a function of interaction range or T(c)* is not constant.

Thermodynamically self-consistent liquid state theories for systems with bounded potentials

The mean spherical approximation (MSA) can be solved semianalytically for the Gaussian core model (GCM) and yields exactly the same expressions for the energy and the virial equations. Taking advantage of this semianalytical framework, we apply the concept of the self-consistent Ornstein-Zernike approximation (SCOZA) to the GCM: a state-dependent function K is introduced in the MSA closure relation which is determined to enforce thermodynamic consistency between the compressibility route and either the energy or virial route. Utilizing standard thermodynamic relations this leads to two differential equations for the function K that have to be solved numerically. Generalizing our concept we propose an integrodifferential-equation-based formulation of the SCOZA which, although requiring a fully numerical solution, has the advantage that it is no longer restricted to the availability of an analytic solution for a particular system. Rather it can be used for an arbitrary potential and even in combination with other closure relations, such as a modification of the hypernetted chain approximation.

Vapor-liquid equilibrium and critical behavior of the square-well fluid of variable range: A theoretical study

The vapor-liquid phase behavior and the critical behavior of the square-well (SW) fluid are investigated as a function of the interaction range, lambdain [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 (lambda>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.

Towards a unification of the hierarchical reference theory and the self-consistent Ornstein-Zernike approximation

The hierarchical reference theory and the self-consistent Ornstein-Zernike approximation are two liquid state theories that both furnish a largely satisfactory description of the critical region as well as phase coexistence and the equation of state in general. Furthermore, there are a number of similarities that suggest the possibility of a unification of both theories. As a first step towards this goal, we consider the problem of combining the lowest order gamma expansion result for the incorporation of a Fourier component of the interaction with the requirement of consistency between internal and free energies, leaving aside the compressibility relation. For simplicity, we restrict ourselves to a simplified lattice gas that is expected to display the same qualitative behavior as more elaborate models. It turns out that the analytically tractable mean spherical approximation is a solution to this problem, as are several of its generalizations. Analysis of the characteristic equations shows the potential for a practical scheme and yields necessary conditions that any closure to the Ornstein-Zernike relation must fulfill for the consistency problem to be well posed and to have a unique differentiable solution. These criteria are expected to remain valid for more general discrete and continuous systems, even if consistency with the compressibility route is also enforced where possible explicit solutions will require numerical evaluations.

Pair correlation function of short-ranged square-well fluids

We have performed extensive Monte Carlo simulations in the canonical (NVT) ensemble of the pair correlation function for square-well fluids with well widths lambda-1 ranging from 0.1 to 1.0, in units of the diameter sigma of the particles. For each one of these widths, several densities rho and temperatures T in the ranges 0.1< or =rhosigma(3)< or =0.8 and T(c)(lambda) less or approximately T less or approximately 3T(c)(lambda), where T(c)(lambda) is the critical temperature, have been considered. The simulation data are used to examine the performance of two analytical theories in predicting the structure of these fluids: the perturbation theory proposed by Tang and Lu [Y. Tang and B. C.-Y. Lu, J. Chem. Phys. 100, 3079 (1994); 100, 6665 (1994)] and the nonperturbative model proposed by two of us [S. B. Yuste and A. Santos, J. Chem. Phys. 101 2355 (1994)]. It is observed that both theories complement each other, as the latter theory works well for short ranges and/or moderate densities, while the former theory works for long ranges and high densities.

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.