Phase boundaries near critical end points. II. General spherical models

Emergence of a bicritical end point in the random-crystal-field Blume-Capel model

We obtain the phase diagram for the Blume-Capel model with the bimodal distribution for random crystal fields, in the space of three fields: temperature (T), crystal field (Δ), and magnetic field (H) on a fully connected graph. We find three different topologies for the phase diagram, depending on the strength of disorder. Three critical lines meet at a tricritical point only for weak disorder. As disorder strength increases there is no tricritical point in the phase diagram. We instead find a bicritical end point, where only two of the critical lines meet on a first-order surface in the H=0 plane. For intermediate strengths of disorder, the phase diagram has critical end points along with the bicritical end point. One needs to look at the phase diagram in the space of three fields to identify various such multicritical points.

Critical charge and density coupling in ionic spherical models

We investigate ionic criticality on the basis of a specially devised spherical model that accounts both for Coulomb and nonionic forces in binary systems. We show in detail here the consequences of the entanglement of density and charge correlation functions G_{NN} and G_{ZZ} on criticality and screening. We also show on this soluble model how, because of electroneutrality, the long-range Coulomb interactions do not change the universality class of criticality in the model driven primarily by sufficiently attractive nonionic interactions. Near criticality, G_{NN} and G_{ZZ} are fully decoupled in charge symmetric systems. However, in more realistic nonsymmetric models, charge and density fluctuations couple in leading order so that the charge and density correlation lengths diverge asymptotically in a similar way. Similarly, the Stillinger-Lovett sum rule, which characterizes a conducting fluid, is violated at criticality in nonsymmetric models when the critical-point density-decay exponent η vanishes. In addition, if quantum effects are accounted for semiclassically by incorporating algebraically decaying interactions, G_{ZZ} decays only as a power law in the whole phase space, contrary to the usually expected exponential Debye screening. We expect these results on this soluble toy model to be general and to reveal general mechanisms ruling ionic criticality.

Criticality in multicomponent spherical models: results and cautions

To enable the study of criticality in multicomponent fluids, the standard spherical model is generalized to describe an S -species hard-core lattice gas. On introducing S spherical constraints, the free energy may be expressed generally in terms of an SxS matrix describing the species interactions. For binary systems, thermodynamic properties have simple expressions, while all the pair correlation functions are combinations of just two eigenmodes. When only hard-core and short-range overall attractive interactions are present, a choice of variables relates the behavior to that of one-component systems. Criticality occurs on a locus terminating a coexistence surface; however, except at some special points, an unexpected "demagnetization effect" suppresses the normal divergence of susceptibilities at criticality and distorts two-phase coexistence. This effect, unphysical for fluids, arises from a general lack of symmetry and from the vectorial and multicomponent character of the spherical model. Its origin can be understood via a mean-field treatment of an XY spin system below criticality.

Three-dimensional Ising model with nearest- and next-nearest-neighbor interactions

The phase diagram of the Ising model in the presence of nearest- and next-nearest-neighbor interactions on a simple cubic lattice is studied within the framework of the differential operator technique. The Hamiltonian is solved by employing an effective-field theory with finite clusters consisting of a pair of spins. A functional form is also proposed for the free energy, similar to the Landau expansion, in order to obtain the phase diagram of the model. The transition from the ferromagnetic (or antiferromagnetic) phase to the disordered paramagnetic phase is of second order. On the other hand, a first-order transition is obtained from the lamellar phase to the paramagnetic phase, as well as from the lamellar phase to the ferromagnetic (or antiferromagnetic) phase, with the presence of a critical end point. An expected singular behavior of the first-order line at the critical end point is also obtained.

Critical endpoint behavior in an asymmetric Ising model: application of Wang-Landau sampling to calculate the density of states

Using the Wang-Landau sampling method with a two-dimensional random walk we determine the density of states for an asymmetric Ising model with two- and three-body interactions on a triangular lattice, in the presence of an external field. With an accurate density of states we were able to map out the phase diagram accurately and perform quantitative finite-size analyses at, and away from, the critical endpoint. We observe a clear divergence of the curvature of the spectator phase boundary and of the magnetization coexistence diameter derivative at the critical endpoint, and the exponents for both divergences agree well with previous theoretical predictions.

Ionic criticality: an exactly soluble model

Gas-liquid criticality in ionic fluids is studied in exactly soluble spherical models that use interlaced sublattices to represent hard-core multicomponent systems. Short-range attractions in the uncharged fluid drive criticality, but charged ions do not alter the universality class. Debye screening remains exponential at criticality in fully ion-symmetric 1:1 models. However, ionic asymmetry couples charge and density fluctuations in a direct manner: The charge correlation length then diverges precisely as the density correlation length and the Stillinger-Lovett rule is violated at criticality.

Phase diagram of symmetric binary fluid mixtures: first-order or second-order demixing

Binary fluid mixtures of 1:1 concentration can demix in a phase transition of first order or of second order. We analyze the two scenarios in density-concentration space and relate them to the structure of the line at which the demixing coexistence surface cuts the liquid-vapor coexistence surface. These scenarios help us to decide between first and second order for a model of a symmetric Lennard-Jones mixture. An optimized reference hypernetted chain integral equation method is employed for calculating the correlation functions and from there the pressure and chemical potentials. We conclude that demixing of a 1:1 mixture is of first order in the whole range of parameters that we have investigated. We did not find a critical point in the 1:1 concentration plane.