Two-Phase Coexistence Simulations

• Two-phase coexistence is a computational method that models stable coexistence between distinct phases to quantify phase boundaries and interfacial properties.
• These simulations employ techniques such as direct coexistence protocols, Gibbs ensemble methods, and field-theoretic frameworks to enforce equilibrium and nonequilibrium criteria.
• They provide actionable insights into phase transitions across systems from active matter to complex mixtures, guiding experimental validations and theoretical advances.

Two-phase coexistence simulations are computational protocols designed to directly realize, maintain, and interrogate the coexistence of two distinct thermodynamic phases—such as gas and liquid, fluid and crystal, or active and passive states—within a single simulation domain. These approaches are foundational for characterizing phase boundaries, measuring coexistence properties, studying interfacial phenomena, and validating theoretical coexistence criteria in both equilibrium and nonequilibrium systems across soft matter, active matter, porous media, and mixtures. In recent years, a suite of simulation techniques based on direct coexistence, advanced ensemble Monte Carlo, and mechanical/field-theoretic frameworks have been developed and applied to a wide array of physical and chemical systems.

1. Conceptual Framework for Two-Phase Coexistence

Two-phase coexistence refers to the stable or metastable realization of domains corresponding to distinct phases (e.g., liquid–gas, fluid–crystal, or banded phases) in a controlled simulation environment. The central requirement is the simultaneous satisfaction of mechanical, chemical, and (for equilibrium cases) thermal equilibrium between the phases, which are separated by an interface. In many modern studies, coexistence is identified and quantified by direct measurement of local densities, order parameters, pressure tensors, and associated structural observables across the interface.

The equilibrium foundation is the Gibbs phase rule, mechanical and chemical potential equality, and the equal-area (Maxwell) construction. Nonequilibrium systems (such as active matter) require generalized mechanical or dynamical criteria for phase coexistence, as detailed in recent theoretical advances (Omar et al., 2022). In quantum or multicomponent systems, further constraints—such as quantum statistics or species conservation—govern the coexistence conditions (Fantoni, 2015).

2. Methodologies for Simulating Two-Phase Coexistence

A variety of computational protocols enable the controlled realization of two-phase coexistence:

• Direct coexistence in periodic simulation cells: Classical approach where the simulation box is initially divided into slabs or regions corresponding to different phases. The interface relaxes and is maintained under periodic boundary conditions (Smallenburg et al., 2024, Wang et al., 7 Oct 2025, Algaba et al., 2024). Key steps include careful initialization, meticulous control of ensemble variables (NVT, NPT), and measurement of spatial profiles.
• Gibbs (Quantum) Ensemble Monte Carlo (GEMC/QGEMC): Utilizes two coupled simulation boxes that exchange particles and/or volume, enabling the phases to form and reach coexistence with rigorously imposed equality of pressure and chemical potentials, including path-integral variants for quantum systems (Fantoni, 2015).
• Slab geometry for inhomogeneous fluids or active matter: ABP and active matter studies employ elongated (slab) boxes to promote planar interfaces and minimize finite-size effects, directly measuring coexisting bulk and interfacial properties (Omar et al., 2022, Perez-Bastías et al., 17 Apr 2025).
• Lattice-based Monte Carlo for complex materials: Phase coexistence is directly realized as spatial domains of distinct states in lattice models of multicomponent systems, such as lipid membranes (Ehrig et al., 2010) or DNA (Fosado et al., 2017).
• Interface-capturing field methods: Phase-field, Allen-Cahn, and Cahn-Hilliard–Navier–Stokes frameworks use order parameters or concentration fields to capture coexistence and interface dynamics, especially in hydrodynamic two-phase flow regimes (Khanwale et al., 2019, Grave et al., 2022).

Each method is tailored to system specifics—e.g., molecular detail versus coarse-grained, equilibrium versus nonequilibrium, presence of hydrodynamics, quantum statistics, and computational scalability.

3. Mechanical and Thermodynamic Criteria for Coexistence

A central output of coexistence simulations is the quantitative determination of binodal points (coexistence densities, pressures, and compositions), interfacial free energies (line/surface tensions), and interface structure.

• Equilibrium systems: Coexistence points are found at the intersection of pressure and/or chemical potential curves for the two phases. Mechanical equilibrium (equal normal stress across the interface) and equality of chemical potential are enforced either intrinsically (GEMC, direct coexistence with interface relaxation) or via common-tangent constructions in free-energy space (Smallenburg et al., 2024, Wang et al., 7 Oct 2025).
• Nonequilibrium active matter: The mechanical theory replaces the equilibrium equal-area construction with a stress-balance framework. Here, the spatially constant normal stress and a generalized equal-area construction in the plane of effective pressure versus conservative pressure fully determine coexistence (Omar et al., 2022). For ABPs and related models, direct measurements of local conservative and swim pressures, polarization, and nematic order profiles inform the binodal and interface characteristics.
• Interfacial phenomena: The width, profile, and energetic cost of the interface are measured, often revealing rich physics—such as nonmonotonic interface broadening in ABP MIPS (Omar et al., 2022), thermal fluctuation spectra in lipid membranes (Ehrig et al., 2010), and topology-controlled scaling in DNA melting (Fosado et al., 2017).

4. Analysis Protocols and Observables

Across simulation methodologies, a suite of data analysis protocols is used to extract key quantities:

• Local fields: Density, polarization, nematic order (in active matter), order parameters (crystallinity, composition), and pressure tensors are resolved spatially, often by binning or coarse-graining (Smallenburg et al., 2024, Omar et al., 2022, Perez-Bastías et al., 17 Apr 2025).
• Order-parameter distributions: Probability density functions (PDFs) of local density, bond order, or other descriptors reveal bimodality and phase fraction in the coexistence window (Packard et al., 2024, Wang et al., 7 Oct 2025).
• Finite-size scaling and interface shape analysis: Binder cumulants, structure factors, and direct interface profile fitting (e.g., hyperbolic tangent for density) characterize the order of the transition, interface width, and correlation lengths (Packard et al., 2024, Ehrig et al., 2010).
• Kinetics and dynamical features: Mean-square displacement (MSD), domain growth exponents, and bubble/droplet coarsening laws provide information on the timescales and mechanisms underlying phase evolution (Roy et al., 2013, Fosado et al., 2017, Ehrig et al., 2010).
• Hydrodynamic and field-theoretic observables: In continuum field schemes (phase-field, level-set, Cahn-Hilliard), interface tracking, local curvature, and energy dissipation are fundamental for capturing two-phase dynamics with or without full Navier–Stokes coupling (Khanwale et al., 2019, Grave et al., 2022).

5. Case Studies Across Model Systems

1. Active Matter (ABP MIPS):

Slab-geometry ABP simulations underlie the nonequilibrium mechanical theory for MIPS. Direct measurement of local density and pressure profiles, with numerical integration of the mechanical equal-area construction, yields coexistence densities and captures unique nonequilibrium interfacial effects (e.g., anomalous interface widening) (Omar et al., 2022). Complementary two-field (density-polarization) theories reproduce both transient and steady-state features observed in simulations (Perez-Bastías et al., 17 Apr 2025).

2. Fluid–Crystal and Hard-Sphere Systems:

NVT direct-coexistence protocols, with pressure-matching and crystallinity analysis (via bond-order parameters), yield precise fluid-crystal coexistence curves. Finite seeding and tuning of particle "hardness" are used to facilitate nucleation and probe entropy-driven coexistence (Smallenburg et al., 2024, Wang et al., 7 Oct 2025).

3. Lattice-Based Models (Lipid Membranes):

MC simulations of binary lipid mixtures on square lattices, with dynamic exchange and single-site flips, capture fluid–gel coexistence, critical phenomena, line tension, and subdiffusive kinetics. The line tension scales as λ∼(1−T/T_c)^μ with μ≈1.17 (Ehrig et al., 2010).

4. Soft and Porous Media:

Two-field, second-gradient models describe phase coexistence in consolidating porous media. Coexistence profiles are computed via mechanical-analogy ODE reductions and are sensitive to second-gradient elastic moduli and external pressure, yielding nonmonotonic interface shapes for specific parameter regimes (Cirillo et al., 2011).

5. Complex and Quantum Mixtures:

The QGEMC algorithm for hydrogen–helium mixtures samples two boxes exchanging volume and particles, imposing equality of pressure and chemical potentials, with full quantum path-integral treatment (Fantoni, 2015).

6. Field-Theoretic Two-Phase Flow:

Cahn–Hilliard–Navier–Stokes and Allen–Cahn/level-set approaches on adaptively refined meshes afford accurate and energy-stable simulation of two-phase coexistence, including mass-conserving and interface-capturing strategies (Khanwale et al., 2019, Grave et al., 2022).

6. Nonequilibrium Coexistence and Emerging Phenomena

Simulation studies in nonequilibrium systems—such as active matter (ABPs), flocking models, and kinetic Ising systems—have revealed phenomena inaccessible in equilibrium:

• Nonmonotonic interface widths originating from self-organized polarization at the interface (Omar et al., 2022), and new coexistence mechanisms such as banded/flocking stripes in multi-species Vicsek models (Lardet et al., 22 Mar 2025, Packard et al., 2024).
• Topology-controlled broadening of DNA melting transitions due to competition between local order parameter and conserved torsional stress (Fosado et al., 2017).
• Dynamical scaling and anomalous coarsening, with growth exponents sensitive to topological constraints or nonequilibrium driving.

These findings motivate and validate generalized mechanical or dynamical coexistence criteria, reshaping conceptual understanding of phase separation in systems far from equilibrium.

7. Practical Guidelines and Limitations

Best practices for two-phase coexistence simulations emphasize:

• Large system sizes to minimize finite-size and interfacial effects, especially near critical points or in systems with giant fluctuations (Smallenburg et al., 2024, Wang et al., 7 Oct 2025, Packard et al., 2024).
• Careful ensemble and boundary condition choices (e.g., fixed box sizes spatially commensurate with equilibrium lattice, periodicity aligned with interface, NVT/NPT selection) (Smallenburg et al., 2024, Algaba et al., 2024).
• Diagnostic evaluation of equilibration (plateaued order parameters, pressure, and phase fractions), consistency across initialization protocols, and reproducibility over independent realizations (Wang et al., 7 Oct 2025).

Method-specific limitations include:

• Potential metastability or kinetic bottlenecks in nucleation, requiring explicit seeding or biasing (Wang et al., 7 Oct 2025).
• Interpretation of coexistence in finite time for glassy or extremely slow systems (pure hard spheres), which may require perturbations to facilitate observable phase separation within feasible computational windows (Wang et al., 7 Oct 2025).
• In field-theoretic approaches, interface thickness and grid/mesh refinement need to be systematically converged (Khanwale et al., 2019, Grave et al., 2022).

In summary, two-phase coexistence simulations, grounded in rigorous mechanical and thermodynamic frameworks, provide quantitative access to binodals, interfaces, and phase behaviors in both equilibrium and nonequilibrium systems. Continued developments in simulation methodologies, scaling analyses, and interfacial diagnostics are central to advancing the predictive capability and theoretical understanding of phase coexistence in complex fluids, soft and active matter, and strongly correlated mixtures.