Multi-dimensional EoS Merger via a Statistical Mixture

We present results from a new 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 $\omega$ is built directly from the input EoSs. The weight of each EoS is determined by minimizing $\omega$ at fixed temperature $T$ and baryon chemical potential $\mu_B$. Thermodynamic consistency and stability are ensured via the single merged grand potential $\omega(T,\mu_B)$ with the correct convexity properties. The code is specifically designed to merge EoS with two independent variables, though the method is not restricted to two independent variables. Here, we include the final results for our merged EoS, along with the two components---a van der Waals Hadron--Resonance--Gas EoS and a holographic Einstein--Maxwell--Dilaton EoS.