Gas-liquid critical point in ionic fluids

Based on the method of collective variables we develop the statistical field theory for the study of a simple charge-asymmetric $1:z$ primitive model (SPM). It is shown that the well-known approximations for the free energy, in particular DHLL and ORPA, can be obtained within the framework of this theory. In order to study the gas-liquid critical point of SPM we propose the method for the calculation of chemical potential conjugate to the total number density which allows us to take into account the higher order fluctuation effects. As a result, the gas-liquid phase diagrams are calculated for $z=2-4$. The results demonstrate the qualitative agreement with MC simulation data: critical temperature decreases when $z$ increases and critical density increases rapidly with $z$.Comment: 18 pages, 1 figur

The density functional theory of classical fluids revisited

We reconsider the density functional theory of nonuniform classical fluids from the point of view of convex analysis. From the observation that the logarithm of the grand-partition function $\log \Xi [\phi]$ is a convex functional of the external potential $\phi$ it is shown that the Kohn-Sham free energy ${\cal A}[\rho]$ is a convex functional of the density $\rho$. $\log \Xi [\phi]$ and ${\cal A}[\rho]$ constitute a pair of Legendre transforms and each of these functionals can therefore be obtained as the solution of a variational principle. The convexity ensures the unicity of the solution in both cases. The variational principle which gives $\log \Xi [\phi]$ as the maximum of a functional of $\rho$ is precisely that considered in the density functional theory while the dual principle, which gives ${\cal A}[\rho]$ as the maximum of a functional of $\phi$ seems to be a new result.Comment: 10 page

A mesoscopic field theory of ionic systems versus a collective variable approach

We establish a link between the two functional approaches: a mesoscopic field theory developed recently by A.Ciach and G.Stell [A. Ciach and G. Stell, J. Mol. Liq. 87 (2000) 253] for the study of ionic models and an exact statistical field theory based on the method of collective variables.Comment: 7 page

Thermo-statistical description of gas mixtures from space partitions

The new mathematical framework based on the free energy of pure classical fluids presented in [R. D. Rohrmann, Physica A 347, 221 (2005)] is extended to multi-component systems to determine thermodynamic and structural properties of chemically complex fluids. Presently, the theory focuses on $D$-dimensional mixtures in the low-density limit (packing factor $\eta < 0.01$). The formalism combines the free-energy minimization technique with space partitions that assign an available volume $v$ to each particle. $v$ is related to the closeness of the nearest neighbor and provides an useful tool to evaluate the perturbations experimented by particles in a fluid. The theory shows a close relationship between statistical geometry and statistical mechanics. New, unconventional thermodynamic variables and mathematical identities are derived as a result of the space division. Thermodynamic potentials $\mu_{il}$, conjugate variable of the populations $N_{il}$ of particles class $i$ with the nearest neighbors of class $l$ are defined and their relationships with the usual chemical potentials $\mu_i$ are established. Systems of hard spheres are treated as illustrative examples and their thermodynamics functions are derived analytically. The low-density expressions obtained agree nicely with those of scaled-particle theory and Percus-Yevick approximation. Several pair distribution functions are introduced and evaluated. Analytical expressions are also presented for hard spheres with attractive forces due to K\^ac-tails and square-well potentials. Finally, we derive general chemical equilibrium conditions.Comment: 14 pages, 8 figures. Accepted for publication in Physical Review

Implementation of the Hierarchical Reference Theory for simple one-component fluids

Combining renormalization group theoretical ideas with the integral equation approach to fluid structure and thermodynamics, the Hierarchical Reference Theory is known to be successful even in the vicinity of the critical point and for sub-critical temperatures. We here present a software package independent of earlier programs for the application of this theory to simple fluids composed of particles interacting via spherically symmetrical pair potentials, restricting ourselves to hard sphere reference systems. Using the hard-core Yukawa potential with z=1.8/sigma for illustration, we discuss our implementation and the results it yields, paying special attention to the core condition and emphasizing the decoupling assumption's role.Comment: RevTeX, 16 pages, 2 figures. Minor changes, published versio

Modelling colloids with Baxter's adhesive hard sphere model

The structure of the Baxter adhesive hard sphere fluid is examined using computer simulation. The radial distribution function (which exhibits unusual discontinuities due to the particle adhesion) and static structure factor are calculated with high accuracy over a range of conditions and compared with the predictions of Percus--Yevick theory. We comment on rigidity in percolating clusters and discuss the role of the model in the context of experiments on colloidal systems with short-range attractive forces.Comment: 14 pages, 7 figures. (For proceedings of "Structural arrest in colloidal systems with short-range attractive forces", Messina, December 2003

A Monte Carlo Approach for Studying Microphases Applied to the Axial Next-Nearest-Neighbor Ising and the Ising-Coulomb Models

The equilibrium phase behavior of microphase-forming systems is notoriously difficult to obtain because of the extended metastability of their modulated phases. In this paper we present a systematic simulation methodology for studying layered microphases and apply the approach to two prototypical lattice-based systems: the three-dimensional axial next-nearest-neighbor Ising (ANNNI) and Ising-Coulomb (IC) models. The method involves thermodynamically integrating along a reversible path established between a reference system of free spins under an ordering field and the system of interest. The resulting free energy calculations unambiguously locate the phase boundaries. The simple phases are not observed to play a particularly significant role in the devil's flowers. With the help of generalized order parameters, the paramagnetic-modulated critical transition of the ANNNI model is also studied. We confirm the XY universality of the paramagnetic-modulated transition and its isotropic nature. Interfacial roughening is found to play at most a small role in the ANNNI layered regime.Comment: 15 pages, 11 figures, 2 table

Field-theoretic description of ionic crystallization in the restricted primitive model

Effects of charge-density fluctuations on a phase behavior of the restricted primitive model (RPM) are studied within a field-theoretic formalism. We focus on a $\lambda$-line of continuous transitions between charge-ordered and charge-disordered phases that is observed in several mean-field (MF) theories, but is absent in simulation results. In our study the RPM is reduced to a $\phi^6$ theory, and a fluctuation contribution to a grand thermodynamic potential is obtained by generalizing the Brazovskii approach. We find that in a presence of fluctuations the $\lambda$-line disappears. Instead, a fluctuation-induced first-order transition to a charge-ordered phase appears in the same region of a phase diagram, where the liquid -- ionic-crystal transition is obtained in simulations. Our results indicate that the charge-ordered phase should be identified with an ionic crystal.Comment: 31 pages, 10 figure

Integral equations for simple fluids in a general reference functional approach

The integral equations for the correlation functions of an inhomogeneous fluid mixture are derived using a functional Taylor expansion of the free energy around an inhomogeneous equilibrium distribution. The system of equations is closed by the introduction of a reference functional for the correlations beyond second order in the density difference from the equilibrium distribution. Explicit expressions are obtained for energies required to insert particles of the fluid mixture into the inhomogeneous system. The approach is illustrated by the determination of the equation of state of a simple, truncated Lennard--Jones fluid and the analysis of the behavior of this fluid near a hard wall. The wall--fluid integral equation exhibits complete drying and the corresponding coexisting densities are in good agreement with those obtained from the standard (Maxwell) construction applied to the bulk fluid. Self--consistency of the approach is examined by analyzing the virial/compressibility routes to the equation of state and the Gibbs--Duhem relation for the bulk fluid, and the contact density sum rule and the Gibbs adsorption equation for the hard wall problem. For the bulk fluid, we find good self--consistency for stable states outside the critical region. For the hard wall problem, the Gibbs adsorption equation is fulfilled very well near phase coexistence where the adsorption is large.For the contact density sum rule, we find some deviationsnear coexistence due to a slight disagreement between the coexisting density for the gas phase obtained from the Maxwell construction and from complete drying at the hard wall.Comment: 29 page