Equilibrium Phase Density Analysis

• Equilibrium phase density is the measure of steady-state mass or particle number per unit volume defined by minimizing thermodynamic potentials.
• It is determined through variational principles, common-tangent constructions, and solving Euler–Lagrange equations across diverse physical models.
• This concept underpins phase transition analysis, interfacial structure predictions, and fluctuation evaluations in systems from quantum gases to cosmological halos.

Equilibrium phase density is the fundamental quantity specifying the steady-state, spatially-resolved mass or particle number per unit volume in each phase of a system at thermodynamic equilibrium. The determination of equilibrium phase densities underpins classical and quantum descriptions of phase transitions, coexistence phenomena, and interfacial structure, spanning systems from condensed matter and statistical mechanics to quantum gases, complex fluids, and cosmological halos. The precise equilibrium densities for each phase are dictated by the minimization of appropriate thermodynamic potentials—typically free energy, grand potential, or Lyapunov-type entropy functionals—subject to the constraints imposed by conservation laws, system geometry, and microscopic interaction parameters associated with the model under study.

1. Thermodynamic Potentials and Variational Principles

Equilibrium densities for each phase arise as solutions to variational stationarity conditions for thermodynamic potentials under prescribed constraints. For example, in classical density functional theory, the equilibrium density field $\rho(\mathbf{r})$ minimizes the grand potential functional

$$\Omega[\rho] = F_{\text{id}}[\rho] + F_{\text{ex}}[\rho] + \int d^3r\,\rho(\mathbf{r})[V_{\text{ext}}(\mathbf{r}) - \mu]$$

where $F_{\text{id}}$ is the ideal-gas contribution, $F_{\text{ex}}$ captures interactions, $V_{\text{ext}}$ is an external field, and $\mu$ is the chemical potential. The stationarity condition $\delta\Omega[\rho]/\delta\rho(\mathbf{r}) = 0$ yields the Euler–Lagrange equation for $\rho(\mathbf{r})$ (Nold et al., 2014).

In phase-field and phase-field-crystal (PFC) models, the equilibrium phase densities are obtained by minimizing effective free-energy functionals $F[\{\phi_i\},\rho]$, where $\phi_i$ are phase indicators or density-order parameters and $\rho$ is the conserved density. The equilibrium conditions are a set of coupled Euler–Lagrange equations, always including $\delta F/\delta\rho = \mu_0$ (constant) and, for non-conserved variables, $\delta F/\delta\phi_i = 0$ (Toth et al., 2019, Kocher et al., 2014).

In systems with multiple conserved species as in binary alloys, equilibrium enforces constant chemical potentials and pressure across coexisting phases. This leads to the common-tangent construction or convex-hull conditions in the $(c,\mathcal{F}+\nu p)$ plane of concentration, molar free energy, and pressure, as implemented in the XPFC formalism (Frick et al., 2020).

2. Calculation of Equilibrium Densities Across Model Classes

Different physical models prescribe distinct functional forms and techniques for equilibrium density determination:

• Landau/Ising-like Mean-Field Theories: Near a critical point, the Helmholtz free energy density is expanded as a function of reduced variables $t = (T-T_c)/T_c$ and $\eta = (n-n_c)/n_c$:

$$f(n,T) = f_0(t) + f_1(t)\,\eta + f_2(t)\,\eta^2 + f_\sigma\,|\eta|^\sigma$$

The equilibrium density is the solution to $\partial f/\partial n = 0$, giving rise to scaling forms for phase coexistence densities:

$$n_{l,g}(T) = n_c[1 \pm A(-t)^{\beta}], \quad T < T_c$$

with critical exponents fixed by universality (e.g., 3D Ising) (Kapusta, 2010).

• Phase-Field Theories: These feature free-energy density functions $f_\rho(\phi,\rho)$ coupling the phase field $\phi$ to the density $\rho$. For instance, in "Model 2" formulations with $f_\rho = (\rho - q(\phi))^2/2$, the equilibrium densities in the bulk phases are set directly by $q(0)$ (liquid) and $q(1)$ (solid), e.g., $\rho_\ell^0 = 1-\epsilon$, $\rho_s^0 = 1+\epsilon$ (Toth et al., 2019).
• Microscopic DFT and Colloidal Microphases: Equilibrium densities in nanostructured and microphase-separated systems follow from minimizing the grand potential or Helmholtz free energy along constrained paths and utilizing thermodynamic integration. Coexistence densities are extracted by common-tangent constructions on the free-energy curves, permitting identification of equilibrium cluster, cylinder, and lamellar phases within precise density windows (Zhuang et al., 2017, Nold et al., 2014).
• Quantum Gases: In the Bose gas, the equilibrium momentum distribution is determined by a minimization of the grand potential in the density-matrix formalism, yielding the Bose–Einstein distribution:

$n_p = \frac{1}{\exp[\beta(E_k-\tilde\mu)] - 1}$

The total density is fixed by the self-consistency relation $\rho = \lambda_T^{-3} g_{3/2}(z)$, and at $T < T_c$ the condensate fraction and normal density $\rho_s$, $\rho_n$ can be written explicitly (Bondarev, 2013).

• Hydrodynamic and Multiphase Flows: For incompressible two-phase fluid systems, phase densities $\rho_1,\rho_2$ are constant in each phase and directly constrained by mass and energy conservation laws, with interface shape and spatial partition (e.g., volume of each phase) determined algebraically from global constraints, e.g.,

$$m\,\omega_n\,R^n = \frac{M_0 - \rho_2 |\Omega|}{\rho_1 - \rho_2}$$

(Pruess et al., 2013).

3. Phase Coexistence, Common-Tangent Constructions, and Criticality

At equilibrium, coexisting phases must have matching chemical potentials and pressure. This requirement manifests as the common-tangent construction on phase free-energy curves. In single-component PFC models, coexistence densities $(n_i, n_j)$ are solutions of: $$\frac{\partial f_i}{\partial n_i}(n_i) = \frac{\partial f_j}{\partial n_j}(n_j), \qquad -f_i(n_i) + \mu_i n_i = -f_j(n_j) + \mu_j n_j$$ with triple-phase coexistence requiring equal chemical potentials and pressures for all three phases (Kocher et al., 2014).

In binary systems and XPFC models, coexistence is enforced by a set of four equations for chemical potentials of both species and pressure across the two phases, together with the lever-rule constraint. Numerical mode-expansion and convex-hull methods yield pairs $(\rho^\alpha, c^\alpha), (\rho^\beta, c^\beta)$ at fixed $(T, p)$ (Frick et al., 2020).

Near a critical point, equilibrium densities satisfy non-analytic scaling with critical exponents, and the difference in densities between coexisting phases (the order parameter) exhibits a power-law vanishing as $T \nearrow T_c$ (Kapusta, 2010). In driven or complex systems (e.g., the ABC model), second-order transitions manifest as bifurcations from uniform to periodic density profiles when the critical temperature is reached, and the explicit form of the periodic order parameter can be traced to the underlying mean-field free energy (Barton et al., 2011).

4. Density Profiles and Interfaces

Beyond bulk phase densities, equilibrium phase density formalism determines interface structure and spatial density profiles in inhomogeneous systems. For 1D planar interfaces in phase-field models, the density profile is given by $\rho(x) = q(\phi_0(x))$ with $\phi_0(x)$ a tanh-profile, independent of coupling parameters (Toth et al., 2019).

In classical DFT, spatially resolved density maps $\rho(x,y)$ across contact lines or interfaces are obtained by solving the self-consistent Euler–Lagrange equations. The adsorption height $h(x)$ characterizes the spatially varying density at equilibrium and provides a coarse-grained link between microscopic DFT profiles and mesoscopic effective interfacial models (Nold et al., 2014).

Pseudo-phase-space density, an important construct in cosmological N-body systems, is defined as $Q(r) = \rho(r)/\sigma^3(r)$, where $\sigma(r)$ is the local velocity dispersion. For CDM halos, equilibrium $Q(r)$ profiles from both simulation and Jeans-based theory are nearly pure power laws, with the exponent set by the underlying density-profile shape parameter (Ludlow et al., 2011).

5. Fluctuations and Finite-Volume Effects

Equilibrium phase density in finite systems features mesoscopic fluctuations, whose variance and distribution are governed by the curvature of the free-energy landscape about the equilibrium. In the Landau expansion around $n_{\rm eq}$,

$$f(n, T) \approx f(n_{\rm eq}, T) + \tfrac12 a(T)(\delta n)^2 + \cdots$$

the probability to observe a fluctuation $\delta n$ in a system of finite volume $V$ is approximately Gaussian: $$\mathcal{P}(\eta) \approx \exp \left[- \frac{V}{2T} a(T) (\delta n)^2 \right]$$ Critical regions (large $\delta = \sigma - 1 \sim 4.8$) exhibit very flat free-energy wells, resulting in large fluctuations and enhanced variance—an effect of particular significance in relativistic heavy-ion collisions near the chiral critical point (Kapusta, 2010).

6. Role in Physical Systems and Applications

Equilibrium phase density models are essential for predicting coexistence compositions in metallurgy and alloy systems (Frick et al., 2020), analyzing phase-separating lattice models (Barton et al., 2011), and understanding hydrodynamic motion and phase boundaries in Navier–Stokes–Stefan problems (Pruess et al., 2013). In quantum fluids, explicit expressions for $\rho_n(T)$ and $\rho_s(T)$ as functions of temperature and system parameters explain superfluid transitions, while in cosmology the nearly universal form of equilibrium pseudo-phase-space density is used for predicting structure formation and related observables (Bondarev, 2013, Ludlow et al., 2011).

7. Methods of Determination and Computational Considerations

Equilibrium phase densities are computed via a range of analytical and numerical techniques:

• Analytic solution of coupled algebraic or integral equations (e.g., Bethe ansatz for 1D gases (Malatsetxebarria et al., 2013), mean-field stationarity).
• Common-tangent or convex-hull constructions on numerical free-energy curves (Kocher et al., 2014, Frick et al., 2020).
• Self-consistent solution of Euler–Lagrange equations in DFT (Nold et al., 2014).
• Mode-expansion and minimization procedures in Landau–Brazovskii field-theoretic models (Zhuang et al., 2017).
• Explicit use of symmetry (e.g., uniform phase versus periodic density waves) to reduce the stationarity problem's complexity (Barton et al., 2011).
• Finite-volume corrections and fluctuation analysis in systems with large but finite $V$ (Kapusta, 2010).

Each methodology matches the system's microphysics, order-parameter structure, and thermodynamic constraints, ensuring accurate phase equilibrium predictions in diverse physical contexts.