---
title: First-Order Hadron–Quark Transition
url: https://www.emergentmind.com/topics/first-order-hadron-quark-phase-transition-fopt
type: topic
---

# First-Order Hadron–Quark Transition

Searching arXiv for recent and relevant papers on first-order hadron–quark phase transitions to ground the article in the cited literature.
First-order hadron–quark phase transition (FOPT) denotes a non-analytic change between hadronic matter and deconfined quark matter characterized by coexistence, a discontinuity in density or energy density, and, depending on the thermodynamic ensemble and conserved charges, either a constant-pressure plateau or a non-constant-pressure mixed phase. In neutron-star phenomenology, heavy-ion matter, and early-Universe cosmology, the FOPT is typically formulated through matched hadronic and quark equations of state (EOSs), with the transition specified by quantities such as the onset density, the latent heat or energy-density jump, and the quark-matter stiffness. Recent work has emphasized both direct inversion of neutron-star observables within constant-speed-of-sound parameterizations and more microscopic or phenomenological constructions based on RMF, NJL, Dyson–Schwinger, Friedberg–Lee, and relativistic density-functional frameworks [2304.07381], [2009.13653], [1903.12336], [2006.08298].

## 1. Definition and thermodynamic structure

A first-order hadron–quark phase transition is conventionally identified by equality of thermodynamic potentials across two phases together with a discontinuous change in density-like observables. In one-component matter, the standard Maxwell construction imposes continuity of pressure and chemical potential at a single transition point, producing a constant-pressure mixed region. In the CSS hybrid-star parameterization, the EOS is written piecewise in energy density,
\[
P(\varepsilon) =
\begin{cases}
P_{\rm HM}(\varepsilon), & \varepsilon<\varepsilon_t,\\
P_t + c_{\rm QM}^2(\varepsilon-\varepsilon_t), & \varepsilon\ge \varepsilon_t,
\end{cases}
\]
with pressure continuity at \(\varepsilon_t\) and chemical-potential equality implied by the Maxwell construction when only one conserved charge is considered [2304.07381].

The principal parameters in this setting are the transition density \(\rho_t\), defined by \(\varepsilon_t=\varepsilon(\rho_t)\) on the hadronic branch, the energy-density jump
\[
\Delta\varepsilon \equiv \varepsilon_t^+ - \varepsilon_t^-,
\]
and the quark-matter squared sound speed \(c_{\rm QM}^2\) [2304.07381]. Closely related CSS formulations instead express the quark branch as
\[
\varepsilon(p)=
\begin{cases}
\varepsilon_{\rm HM}(p),&p<p_t,\\
\varepsilon_{\rm HM}(p_t)+\Delta\varepsilon+\dfrac{p-p_t}{c_{\rm QM}^2},&p>p_t,
\end{cases}
\]
again making the latent heat and stiffness explicit [2009.13653].

In systems with two conserved charges, the mixed phase is not generally a constant-pressure plateau. Under Gibbs conditions one requires equality of temperature, baryon chemical potential, electron or charge chemical potential, and pressure, while charge neutrality is imposed globally rather than locally. In this case pressure varies continuously with the quark volume fraction \(\chi\), and the mixed phase exhibits intrinsically non-linear structure [1401.0865], [1102.4964]. This distinction is central in comparing Maxwell-type neutron-star constructions to heavy-ion or multicomponent stellar matter calculations.

The latent heat may be represented either as an energy-density discontinuity \(\Delta\varepsilon\) in compact-star EOS models or as a coexistence-region discontinuity in \(\varepsilon\), entropy density, or baryon density in finite-temperature QCD-motivated descriptions. In an analytic mean-field Ising mapping to QCD, discontinuities \(\Delta n_B\) and \(\Delta s\) across the coexistence line imply the Clapeyron relation
\[
\frac{dT_{\rm coex}}{d\mu_B}=\frac{\Delta n_B}{\Delta s},
\]
making the FOPT part of a broader coexistence and spinodal structure [2511.14891].

## 2. EOS constructions and phase-equilibrium prescriptions

A large fraction of the literature models the FOPT through two-phase EOS matching. In compact-star applications, hadronic matter has been described by RMF or meta-modeling approaches, while quark matter is represented by CSS, Dyson–Schwinger, NJL, string-flip or confinement density-dependent mass models [2009.13653], [1904.01978], [2009.10846], [1503.02739]. The transition is then implemented by Maxwell or Gibbs conditions depending on the conserved-charge content.

In CSS-based neutron-star studies, the hadronic EOS is fixed or sampled, then matched at a transition pressure \(p_t\) to a quark branch with constant \(c_{\rm QM}^2\). The three independent FOPT parameters are commonly taken as \(\rho_t/\rho_0\), \(\Delta\varepsilon/\varepsilon_t\) or \(\Delta\varepsilon\), and \(c_{\rm QM}^2\) [2009.13653], [2304.07381], [2006.00839]. In one implementation, the priors are
\[
\rho_t/\rho_0\in[1,6],\quad
\Delta\varepsilon/\varepsilon_t\in[0.2,1],\quad
c_{\rm QM}^2\in[0,1],
\]
with Maxwell matching through continuity of \(p\) and \(\mu\) [2009.13653].

More microscopic quark-sector descriptions alter both the coexistence region and the existence of critical endpoints. In the RMF + lattice-motivated QGP treatment, the hadronic sector includes the full baryon octet, pions, and kaons interacting through \(\{\sigma,\sigma^*,\omega,\rho,\phi\}\), while the quark–gluon sector is an MIT-bag model with \({\cal O}(\alpha_s)\) corrections. Solving \(P_{\rm H}=P_{\rm Q}\) yields a first-order line ending at a critical endpoint near \((\mu_B^{\rm CEP},T^{\rm CEP})\simeq(360\,\mathrm{MeV},154\,\mathrm{MeV})\) [1201.3710].

In asymmetric matter relevant to heavy-ion collisions, the two-EOS RMF + NJL approach gives Gibbs equilibrium under fixed global asymmetry. The mixed phase is characterized by total densities
\[
\rho_B=(1-\chi)\rho_B^H+\chi\rho_B^Q,\qquad
\rho_3=(1-\chi)\rho_3^H+\chi\rho_3^Q,
\]
and the coexistence region terminates around
\[
T_{\rm CEP}\approx 80\,\mathrm{MeV},\qquad \rho_{B,\rm CEP}\approx 2\,\rho_0,
\]
with the FOPT occurring for \(T\simeq 50\mbox{--}80\) MeV and \(\rho_B\simeq 2\mbox{--}4\,\rho_0\) [1102.4964].

Alternative constructions replace explicit mixed-phase microphysics with interpolation. In the Dyson–Schwinger/RMF framework, the sound-speed interpolation scheme parametrizes
\[
c_s^2(\mu)=A\mu^2+B\mu+C,\qquad \mu\in[\mu_0,\mu_1],
\]
and reconstructs the coexistence-window EOS through
\[
\frac{d\rho_B}{d\mu}=\frac{\rho_B}{\mu c_s^2(\mu)},\qquad \frac{dp}{d\mu}=\rho_B(\mu),
\]
matched to hadronic and quark branches at \(\mu_0,\mu_1\) [1904.01978], [1903.12336].

Finite-size and surface-tension effects motivate mixed-phase smoothing. In a twin-star study based on the ACB4 piecewise-polytropic EOS, the sharp Maxwell plateau is replaced by a quadratic interpolation in chemical potential,
\[
P_M(\mu)=\alpha_2(\mu-\mu_c)^2+\alpha_1(\mu-\mu_c)+(1+\Delta p)P_c,
\]
where \(\Delta p\) controls the degree of pasta-induced smoothing; \(\Delta p=0\) reproduces the sharp transition [2309.08775].

## 3. Neutron-star constraints and inversion of transition parameters

Recent neutron-star studies have shifted from forward modeling to inverse constraints on FOPT parameters. A CSS inversion in a three-dimensional \((c_{\rm QM}^2,\rho_t,\Delta\varepsilon)\) parameter space fixes all hadronic EOS parameters at their best-known values, constructs the hybrid EOS for each parameter triplet, solves the TOV equations, and computes \(M_{\rm max}\), \(R_{2.01}\), and \(\Lambda_{1.4}\). Imposing
\[
M_{\rm max}>2.01\,M_\odot,\qquad
R_{2.01}>11.41~\mathrm{km}\ \text{or}\ 12.2~\mathrm{km},\qquad
\Lambda_{1.4}<427
\]
delimits an allowed “tube” in CSS parameter space [2304.07381].

The most striking result of that inversion is a lower limit on quark-matter stiffness above the conformal bound:
\[
c_{\rm QM}^2 \ge 0.35 \quad \text{for } R_{2.01}>11.41~\mathrm{km},
\]
\[
c_{\rm QM}^2 \ge 0.43 \quad \text{for } R_{2.01}>12.2~\mathrm{km}.
\]
No meaningful upper limit on \(c_{\rm QM}^2\) is obtained from current data [2304.07381]. This indicates that, within that framework, quark matter in neutron stars must be stiffer than a conformal gas.

A Bayesian analysis with nine EOS parameters—six hadronic meta-model parameters and three CSS parameters—finds posterior peaks at
\[
\rho_t/\rho_0=1.6^{+1.2}_{-0.4},\qquad
\Delta\varepsilon/\varepsilon_t=0.40^{+0.20}_{-0.15},\qquad
c_{\rm QM}^2=0.95^{+0.05}_{-0.35}
\]
at \(68\%\) confidence level [2009.13653]. The same study reports that the total probability of \(c_{\rm QM}^2\le 1/3\) is very small, and that the posterior distributions of hadronic and quark EOS parameters are very weakly correlated.

Within the quark-mean-field + CSS framework, systematic exploration of \((n_{\rm trans}/n_0,\Delta\varepsilon/\varepsilon_{\rm trans},c_{\rm QM}^2)\) together with the symmetry-energy slope \(L\) shows that increasing either \(n_{\rm trans}/n_0\) or \(\Delta\varepsilon/\varepsilon_{\rm trans}\) softens the core and reduces \(M_{\rm max}\), while a stiffer quark phase raises \(M_{\rm max}\) and yields more compact stars [2006.00839]. That analysis quotes robust limits \(M_{\rm max}<3.6\,M_\odot\) and \(R_{1.4}\gtrsim 9.6\) km, and finds that a too-weak transition,
\[
\Delta\varepsilon/\varepsilon_{\rm trans}\lesssim 0.2,
\]
taking place at low densities \(\lesssim 1.3\mbox{--}1.5\,n_0\) is strongly disfavored [2006.00839].

Sound-speed interpolation combined with astronomical calibration yields a complementary description in terms of transition chemical potentials rather than direct CSS parameters. In one such analysis, the most probable values are
\[
\langle\mu_0\rangle=(1.00\pm0.05)\,\mathrm{GeV},\qquad
\langle\mu_1\rangle=(1.51\pm0.07)\,\mathrm{GeV},
\]
corresponding to \(\langle n_0\rangle=1.50\,n_s\) and \(\langle n_1\rangle=7.9\,n_s\) [1904.01978]. A closely related study using extensive sampling finds peak values near
\[
\mu_0\simeq 1.00\pm0.01~\mathrm{GeV},\qquad
\mu_1\simeq 1.53\pm0.01~\mathrm{GeV},
\]
with modest sensitivity to interpolation order and hyperon content [1903.12336].

These results collectively favor early onset and moderate latent heat together with a very stiff quark branch. A plausible implication is that present multimessenger data are more compatible with moderate first-order transitions than with either vanishingly weak crossovers in the same parameterizations or very strong latent-heat jumps that destabilize the stellar sequence.

## 4. Correlations, critical criteria, and the twin-star problem

The structure of the FOPT in neutron stars is strongly constrained by stability criteria. In CSS-type models, a necessary condition for a connected stable hybrid branch is the Seidov inequality
\[
\frac{\Delta\varepsilon_{\rm crit}}{\varepsilon_{\rm trans}}
=
\frac12+\frac32\frac{p_{\rm trans}}{\varepsilon_{\rm trans}},
\]
with larger jumps allowing a disconnected third family [2006.00839]. The Bayesian study of canonical stars used the Seidov condition in its default likelihood filter, although removing it changed the posterior only slightly [2009.13653].

An important empirical finding is the approximately linear anti-correlation between transition density and latent heat at the causality boundary \(c_{\rm QM}^2=1\). The upper envelope of allowed points in \((\rho_t/\rho_0,\Delta\varepsilon)\) is fitted by
\[
\Delta\varepsilon~[\mathrm{MeV\,fm^{-3}}] = -112.91(\rho_t/\rho_0)+345.86
\]
for \(R_{2.01}=11.41\) km, and
\[
\Delta\varepsilon~[\mathrm{MeV\,fm^{-3}}] = -119.89(\rho_t/\rho_0)+296.33
\]
for \(R_{2.01}=12.2\) km, both with Pearson \(r\simeq 0.9994\) [2304.07381]. This suggests a narrow tradeoff: later transitions permit smaller latent heat if the EOS is to satisfy current mass-radius-tidal constraints.

Twin-star scenarios remain a focal controversy. A commonly cited heuristic is that twins require a strong first-order transition, typically \(\Delta\varepsilon\gtrsim350\,\mathrm{MeV\,fm^{-3}}\) [2304.07381]. Yet the inversion study finds upper limits
\[
\Delta\varepsilon\lesssim231\,\mathrm{MeV\,fm^{-3}}
\]
for \(R_{2.01}>11.41\) km and
\[
\Delta\varepsilon\lesssim175\,\mathrm{MeV\,fm^{-3}}
\]
for \(R_{2.01}>12.2\) km; a survey of 137 hybrid EOS mass–radius curves in the allowed region showed no twin-star branches [2304.07381].

By contrast, models purposely engineered for large discontinuities or explicit mixed-phase control can generate third-family branches. In the ACB4 plus pasta interpolation framework, \(\Delta p=0\) produces a vertical jump in \(f(M)\) and \(\tau(M)\) at branch onset, whereas increasing \(\Delta p\) up to approximately \(5\mbox{--}8\%\) smooths the transition and eventually removes the third family [2309.08775]. In the excluded-volume/chiral hybrid model, an abrupt first-order-like jump to pure quark matter at \(a=0.8\) produces a plateau in \(P(\varepsilon)\) at \(P_c\simeq 0.3\,\mathrm{GeV\,fm^{-3}}\) and a latent heat \(\Delta\varepsilon\simeq 0.5\,\mathrm{GeV\,fm^{-3}}\), which may produce a “mass twin” phenomenon if the latent heat is large enough [2310.09738].

The apparent tension is model-dependent rather than contradictory. Present observational inversion in constrained CSS spaces disfavors twins [2304.07381], while broader phenomenological constructions can still realize them if the transition is sufficiently abrupt or if mixed-phase smoothing is weak [2309.08775], [2310.09738].

## 5. Bubble nucleation, spinodals, and finite-temperature dynamics

Finite-temperature FOPT dynamics are often formulated in terms of nucleation and metastability. In the Friedberg–Lee model, the tree-level potential
\[
U(\sigma)=\frac12 a\sigma^2+\frac16 b\sigma^3+\frac1{24}c\sigma^4
\]
has two minima when \(b^2>3ac\), and the cubic term \(b\sigma^3\) generates a barrier even at \(T=0\), producing a strong FOPT [2006.08298]. Including one-loop thermal and chemical-potential corrections yields an effective potential \(V_{\rm eff}(\sigma;T,\mu)\) with degenerate minima at the coexistence point.

The critical bubble is the \(O(3)\)-symmetric bounce solving
\[
\frac{d^2\sigma}{dr^2}+\frac{2}{r}\frac{d\sigma}{dr}=\frac{\partial V_{\rm eff}}{\partial\sigma},
\]
with \(\sigma(r\to\infty)=0\) and \(d\sigma/dr|_{r=0}=0\) [2006.08298]. In the thin-wall limit, the bubble free energy is
\[
F(R)\simeq 4\pi R^2\Sigma-\frac{4\pi}{3}R^3\varepsilon,
\]
leading to the critical radius
\[
R_c(T,\mu)=\frac{2\Sigma(T,\mu)}{\varepsilon(T,\mu)}.
\]
The surface tension and bubble free-energy shift are
\[
\Sigma(T,\mu)=\int_0^\infty dr\left[\frac12\left(\frac{d\sigma_b}{dr}\right)^2+V_{\rm eff}(\sigma_b;T,\mu)\right],
\]
\[
\Delta F_b(T,\mu)=4\pi\int_0^\infty dr\,r^2\left[\frac12\left(\frac{d\sigma_b}{dr}\right)^2+V_{\rm eff}(\sigma_b;T,\mu)\right],
\]
and the thermal nucleation rate is approximated by
\[
\Gamma \simeq T^4 \exp[-\Delta F_b(T,\mu)/T].
\]
For the strong transition considered there, \(\Sigma(T)\) does not go to zero at low \(T\), \(R_c\) remains finite as \(T\to0\), and \(\Delta F_b/T\) is non-monotonic, unlike in weak FOPT scenarios approaching spinodal decomposition [2006.08298].

Spinodal physics has re-emerged in analytic QCD EOS constructions. A mean-field Ising mapping to the QCD phase diagram reconstructs both coexistence and spinodal regions by solving the cubic Landau equation
\[
h = a r M + b M^3
\]
and locating stability loss through
\[
\frac{\partial h}{\partial M}=a r + 3 b M^2 =0.
\]
This yields explicit spinodal curves \(T_{\rm sp}^{(\pm)}(\mu_B)\), providing access to super-heated hadronic matter and super-cooled quark–gluon plasma [2511.14891]. Such formulations are intended for hydrodynamic simulations at low collision energy where metastable trajectories matter.

This finite-temperature perspective clarifies a frequent misconception: the coexistence line does not exhaust FOPT physics. Metastable and unstable domains, nucleation barriers, and spinodal instabilities are essential for dynamical realizations in heavy-ion collisions and, by analogy, in astrophysical transient environments [2006.08298], [2511.14891].

## 6. Physical settings and observational signatures

The FOPT has been investigated in three principal arenas: neutron stars, heavy-ion collisions, and the early Universe. In neutron stars, the main observables are maximum mass, radii, tidal deformabilities, and oscillation spectra [2304.07381], [2009.13653], [2309.08775]. In heavy-ion matter, emphasis falls on coexistence regions, spinodals, isospin distillation, and critical endpoints [1102.4964], [2511.14891]. In cosmology, the focus is on transition timing, bubble growth, and the modification of thermodynamic evolution by non-standard gravity [1401.5029], [1201.0372].

Gravitational-wave asteroseismology has been proposed as a probe of the transition character. In the full general relativistic \(f\)-mode analysis of sharp and smooth transitions, the sharp-Maxwell case produces a gap in the \(f(M)\) curve around onset mass: two stars of the same mass but with \(\Delta R\sim1\) km have \(f\) differing by \(\sim5\%\). With Einstein Telescope-level errors of \(\Delta f_f/f_f\lesssim1\%\) and \(\Delta\tau_f/\tau_f\lesssim5\%\), twin branches could be distinguished if \(\Delta R\gtrsim0.4\) km at \(2\sigma\); smooth transitions with \(\Delta p\gtrsim5\%\) erase these jumps [2309.08775].

Heavy-ion studies of asymmetric matter predict isospin-sensitive signals of the mixed phase, including enhanced \(\pi^-/\pi^+\) and \(K^0/K^+\) ratios, increased \(\Delta^{-}/\Delta^{++}\) and hyperon yields, and event-by-event neutron/proton fluctuations due to isospin distillation [1102.4964]. In the analytic Ising-based EOS, the spinodal region permits super-heating and super-cooling and may seed nucleation and cavitation in hydrodynamic evolution at \(\sqrt{s_{NN}\lesssim20}\) GeV [2511.14891].

Astrophysical table-based EOS models with explicit first-order transitions have also been deployed in core-collapse supernovae and binary neutron-star mergers. In the phenomenological quark–hadron EOS, the FOPT can trigger a second shock wave in core-collapse supernova simulations and produce a distinct neutrino burst when the quark core forms; in binary mergers it can shift the post-merger gravitational-wave frequency \(f_2\) and alter collapse outcomes [2009.10846].

In early-Universe applications, first-order quark–hadron transitions have been studied in DGP + Brans–Dicke and chameleon Brans–Dicke brane cosmologies. In the DGP + BD model, using quark and hadron EOSs matched at \(T_c\simeq125\) MeV for \(m_s=200\) MeV and \(B^{1/4}=200\) MeV, the transition occurs on microsecond scales, with faster completion on the self-accelerating branch [1401.5029]. In the chameleon Brans–Dicke brane model, the modified Friedmann law \(H^2\propto\rho^2\) pushes the QCD transition into the nanosecond regime, with \(t_c\sim0.05\mbox{--}0.25\) ns and completion around \(t\sim1\) ns [1201.0372]. These cosmological results are model-specific and should not be conflated with standard FRW timing.

## 7. Open issues and recurrent points of contention

Several unresolved issues organize current debate. The first concerns the stiffness of deconfined matter. Multiple neutron-star inferences place the favored \(c_{\rm QM}^2\) well above the conformal limit \(1/3\), with lower bounds \(0.35\) or \(0.43\) in one inversion and a posterior peak at \(0.95^{+0.05}_{-0.35}\) in another [2304.07381], [2009.13653]. This does not establish a unique microscopic mechanism, but it does constrain viable quark-matter descriptions within the adopted hybrid-star frameworks.

A second issue is whether the mixed phase should be sharp or smooth. Maxwell constructions enforce abrupt jumps at a single pressure, while Gibbs constructions and pasta-inspired interpolations yield extended non-constant-pressure coexistence. Surface tension, finite-size effects, and the number of conserved charges determine which description is appropriate in a given setting [1401.0865], [2309.08775]. Therefore, statements about “the” hadron–quark transition must be qualified by thermodynamic context.

A third issue concerns the critical endpoint. Some models locate a CEP near \((T,\mu_B)\simeq(154\,\mathrm{MeV},360\,\mathrm{MeV})\) [1201.3710], while NJL-based asymmetric-matter calculations place it near \(T\approx80\) MeV and \(\rho_B\approx2\rho_0\) [1102.4964]. A plausible implication is that CEP coordinates are highly model-dependent because they arise in a regime where chiral symmetry breaking, deconfinement, and dense-matter interactions compete strongly.

A fourth issue is whether twin stars are generic signatures of FOPT. Current observational inversion under tightly constrained CSS assumptions makes them improbable [2304.07381], whereas broader EOS classes still permit them under sufficiently large latent heat or sufficiently sharp transitions [2309.08775], [2310.09738]. The more precise statement is that twin stars are possible in some FOPT realizations but not generically implied by the existence of a hadron–quark transition.

Finally, there is a methodological divide between direct microphysical constructions and inference-driven phenomenology. Microphysical models specify the origin of the barrier, the role of chiral restoration, or the behavior of metastable branches [2006.08298], [1003.2843], [2310.09738]. Inference approaches instead ask which transition parameters are compatible with current astrophysical data [2304.07381], [2009.13653], [1903.12336]. Together they indicate that, while the existence of a first-order hadron–quark transition remains unconfirmed, a consistent research program now links coexistence thermodynamics, compact-star observables, heavy-ion signatures, and finite-temperature nucleation dynamics within a shared parameter language of onset, latent heat, stiffness, and phase stability.

Source: https://www.emergentmind.com/topics/first-order-hadron-quark-phase-transition-fopt