Merging multidimensional equations of state of strongly interacting matter via a statistical mixture
Yumu Yang, Prachi Garella, Musa R. Khan, Tulio E. Restrepo, Joaquin Grefa, Johannes Jahan, Mauricio Hippert, Jorge Noronha, Claudia Ratti, Romulo Rougemont
TL;DR
This work introduces a thermodynamically consistent two-fluid mixing framework to merge distinct equations of state (EoSs) for strongly interacting matter by minimizing a grand potential density $ω(T,μ_B)$ with respect to an internal order-parameter $p$. The method replaces naive switching with an entropy- and interaction-driven mixing term, enabling crossover, a critical point with mean-field Ising universality, and a first-order line, while preserving convexity and stability. As a concrete application, the authors merge a van der Waals hadron-resonance-gas EoS with a holographic Einstein–Maxwell–Dilaton EoS, yielding a single global EoS that spans hadronic and deconfined matter up to $μ_B\sim 1$ GeV and $T\sim 600$ MeV, in good agreement with lattice QCD at $μ_B=0$ and with $T'$ expansions at finite density. The resulting EoS supports realistic hydrodynamic simulations and provides a robust, extensible framework for exploring the QCD phase diagram, including the possible observability of phase-transition signatures in heavy-ion collisions and neutron-star phenomena.
Abstract
We introduce a general method to merge multidimensional equations of state (EoSs) by combining them in a two-fluid equilibrium statistical mixture in the grand canonical ensemble. The merged grand potential density $ω$ is built directly from the input EoSs and the fluid fractions are fixed by minimizing $ω$ at fixed temperature $T$ and baryon chemical potential $μ_B$. Thermodynamic consistency and stability are guaranteed as all thermodynamic quantities are consistently derived from a single merged grand potential $ω(T,μ_B)$ with the correct convexity properties. Our method can accommodate a first-order phase transition and a critical endpoint with mean-field critical exponents. We use this method to merge a van der Waals Hadron-Resonance-Gas EoS with a holographic Einstein-Maxwell-Dilaton EoS that has a critical point and a first-order line. The result is a single EoS, spanning hadronic and deconfined matter over a broad range in $(T,μ_B)$, which can be readily used in heavy-ion hydrodynamic simulations. Our merging method can be generalized to consider a higher dimensional phase diagram (e.g., by considering more chemical potentials) and more than two input EoSs.
