Phase stability in the three-dimensional open-source code for the chiral mean-field model

In this paper we explore independently for the first time three chemical potentials (baryon ${\ensuremath{\mu}}_{B}$, charged ${\ensuremath{\mu}}_{Q}$, and strange ${\ensuremath{\mu}}_{S}$) in the chiral mean field (CMF) model. We designed and implemented cmf++, a new version of the CMF model rewritten in c++ that is optimized, modular, and well documented. cmf++ has been integrated into the MUSES Calculation Engine as a free and open-source software module. The runtime improved in more than 4 orders of magnitude across all 3 chemical potentials, when compared to the legacy code. Here we focus on the zero temperature case and study stable, as well as metastable and unstable, vacuum, hadronic, and quark phases, showing how phase boundaries vary with the different chemical potentials. Due to the significant numerical improvements in cmf++, we can now for the first time sweep the entire ${\ensuremath{\mu}}_{B}$, ${\ensuremath{\mu}}_{S}$, ${\ensuremath{\mu}}_{Q}$ phase space, investigate metastable phases, and calculate high-order susceptibilities within the CMF framework. This allows us to find phases of matter that include a light hadronic phase, a strangeness-dominated hadronic phase, and a quark phase. The numerical improvements also allow us to identify the order of the transitions among these phases, finding a first-order chiral symmetry restoration phase transition among the hadronic phases (favored for negative ${\ensuremath{\mu}}_{Q}$ and/or ${\ensuremath{\mu}}_{S}$ for some coupling schemes), in addition to third-order phase transitions. In particular, we identify for the first time triple points in the CMF model, where both chiral symmetry restoration and deconfinement phase transitions meet in the chemical potential phase space. Such points could potentially be identified in low-energy heavy-ion collisions.