Hybrid EOS Framework for Neutron Stars

• Hybrid EOS Framework is a parameterized approach that unites lower-density hadronic matter with high-density deconfined quark matter through first-order transitions or crossovers.
• It employs a constant-speed-of-sound model where key parameters—transition pressure, energy density jump, and quark matter sound speed—govern stellar mass-radius relations.
• The framework links microscopic dense-QCD models with observable constraints from pulsars and gravitational waves, highlighting phenomena like mass twins and third-family branches.

Hybrid EOS Framework denotes a family of neutron-star equation-of-state constructions in which a hadronic sector at lower density is joined to a deconfined quark sector at higher density. In the standard compact-star setting, the matter is cold, catalyzed, and in beta equilibrium, and stellar configurations are obtained in general relativity from the Tolman–Oppenheimer–Volkoff equations. In its most widely used form, the framework is organized around a first-order hadron–quark transition with either a sharp Maxwell interface or a generalized mixed-phase or crossover construction; in the constant-speed-of-sound formulation, the high-density branch is controlled by the transition pressure $p_{\rm trans}$, the energy-density discontinuity $\Delta\epsilon$, and the quark-matter sound speed $c_s^2$, which together determine onset stability, mass–radius topology, and the existence of disconnected “third family” branches (Alford et al., 2013, Ji et al., 8 May 2025).

1. Thermodynamic setting and stellar structure

The framework is formulated for zero-temperature matter in beta equilibrium and charge neutrality. In the hadronic sector, the equilibrium conditions are

$\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$

with global neutrality $\sum_i q_i n_i = 0$, and for each baryon species $i$,

$\mu_i = b_i \mu_B - q_i \mu_e.$

At the thermodynamic level,

$$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$

with leptons $l=e,\mu$ where present (Ji et al., 8 May 2025).

For nonrotating, spherically symmetric stars, the EOS enters the Tolman–Oppenheimer–Volkoff system,

$$\frac{dP}{dr} = - \frac{G (\epsilon + P/c^2)\left(m + 4\pi r^3 P/c^2\right)}{r\left(r - 2Gm/c^2\right)},\qquad \frac{dm}{dr} = 4\pi r^2 \epsilon.$$

Choosing a central pressure or central energy density and integrating outward until $\Delta\epsilon$0 yields the radius $\Delta\epsilon$1 and gravitational mass $\Delta\epsilon$2. Varying the central condition generates the full mass–radius sequence $\Delta\epsilon$3 (Alford et al., 2013).

This setup fixes the hybrid-EOS problem as a coupled thermodynamic and stellar-structure problem. Microphysics specifies $\Delta\epsilon$4 or equivalently $\Delta\epsilon$5; the TOV equations translate that input into observables such as $\Delta\epsilon$6, $\Delta\epsilon$7, compactness, and—after linear tidal perturbations—Love numbers and dimensionless tidal deformabilities.

2. Two-phase construction and CSS parameterization

In the canonical two-phase realization, the hadronic EOS is used below a transition point and a quark EOS above it. For a sharp first-order transition, Maxwell matching imposes equality of pressure and baryon chemical potential at the phase boundary,

$\Delta\epsilon$8

while the energy density jumps by

$\Delta\epsilon$9

The simplest widely used quark parameterization is the constant-speed-of-sound form

$c_s^2$0

with $c_s^2$1. In the Maxwell-constructed hybrid EOS this becomes

$c_s^2$2

so the three control quantities are $c_s^2$3, $c_s^2$4, and $c_s^2$5 (Alford et al., 2013).

The same construction can be written in chemical-potential form. Introducing $c_s^2$6, one convenient CSS representation is

$c_s^2$7

which guarantees thermodynamic consistency and makes the matching to hadronic matter explicit (Alford et al., 2013).

Within the broader Hybrid EOS Framework, CSS is not the only option. Reviews and model studies also employ MIT bag, NJL, FRG quark–meson, V-QCD, multiquark mean-field, and Dyson–Schwinger quark sectors, and they extend the transition modeling from Maxwell to Gibbs mixed phases and quark–hadron crossover or interpolation schemes. Examples include pressure interpolation $c_s^2$8, speed-of-sound interpolation $c_s^2$9, and two-zone parabolic interpolation in $\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$0 over a finite matching window $\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$1 (Ji et al., 8 May 2025, Ayriyan et al., 2021, Ivanytskyi et al., 2022, Zhao et al., 2015).

3. Stability at deconfinement and branch topology

The onset of a quark core is controlled by the Seidov criterion. At $\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$2, the quark core is infinitesimal, so stability depends only on the energy-density jump: $\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$3 If this inequality holds, the hybrid branch is stably connected to the hadronic branch at onset; otherwise the sequence develops a cusp and becomes unstable immediately after deconfinement (Alford et al., 2013).

At fixed $\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$4, the CSS phase diagram is organized in the plane

$\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$5

The resulting branch topologies are conventionally classified as follows:

Region	Meaning
A	no hybrid branch
C	single connected hybrid branch
D	single disconnected hybrid branch
B	both connected and disconnected branches

The straight Seidov line

$\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$6

separates $\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$7 from $\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$8. A second, near-vertical boundary marks the appearance or disappearance of a disconnected branch through a stationary inflection in $\mu_n = \mu_p + \mu_e,\qquad \mu_\mu = \mu_e,$9, and its location depends strongly on $\sum_i q_i n_i = 0$0. A roughly horizontal boundary separates $\sum_i q_i n_i = 0$1 from $\sum_i q_i n_i = 0$2, where a stationary inflection within the connected branch turns into a maximum–minimum pair and creates an additional disconnected branch. A semi-analytic estimate for this split is

$\sum_i q_i n_i = 0$3

A notable feature of the CSS analysis is that the Seidov line is independent of $\sum_i q_i n_i = 0$4 and largely insensitive to the nuclear-matter EOS, whereas the existence and extent of disconnected branches depend strongly on $\sum_i q_i n_i = 0$5 (Alford et al., 2013).

Disconnected hybrid branches are the compact-star realization of the “third family.” When such a branch overlaps in mass with a hadronic branch, the framework admits mass twins: stars with the same gravitational mass but different radii and internal composition. Bayesian and excluded-volume studies identify these high-mass twins as a characteristic signature of a strong first-order deconfinement transition (Alvarez-Castillo et al., 2015, Alvarez-Castillo et al., 2016, Alvarez-Castillo et al., 2015).

4. Observable consequences and empirical constraints

The maximum hybrid-star mass is controlled by the stiffness of the high-density phase and the transition location. A useful scaling is

$\sum_i q_i n_i = 0$6

where $\sum_i q_i n_i = 0$7 is the central energy density of the maximum-mass star and $\sum_i q_i n_i = 0$8 is the effective sound speed in the core. For hybrid stars with large quark cores, $\sum_i q_i n_i = 0$9, so larger quark-matter sound speeds increase $i$0 (Alford et al., 2013).

Representative conditions that yield $i$1 in CSS studies are a transition density not too high, $i$2, a moderate jump $i$3, and stiff quark matter with $i$4. Values near or above the conformal limit $i$5 help, but often a larger $i$6 is required. For a soft nuclear EOS such as HLPS, hybrid stars can exceed the pure-hadronic maximum if the transition occurs early and the quark core is large and stiff; for a stiff hadronic EOS such as NL3, hybrid stars are typically lighter than the pure-hadronic maximum, although $i$7 hybrids remain readily obtainable for favorable quark parameters (Alford et al., 2013).

The framework connects directly to multimessenger observables. Hybridization generally reduces $i$8 and $i$9 at fixed mass when a softer quark core appears at moderate density, and a strong first-order transition can produce a kink or flattening in the mass–radius curve together with a reduction of $\mu_i = b_i \mu_B - q_i \mu_e.$0 into the GW170817-compatible range. Current benchmark constraints quoted in recent reviews are $\mu_i = b_i \mu_B - q_i \mu_e.$1, $\mu_i = b_i \mu_B - q_i \mu_e.$2, and $\mu_i = b_i \mu_B - q_i \mu_e.$3 (Ji et al., 8 May 2025).

Branch observability depends not only on stability but also on length in $\mu_i = b_i \mu_B - q_i \mu_e.$4 and $\mu_i = b_i \mu_B - q_i \mu_e.$5. Even when the Seidov criterion predicts a connected branch, the segment may be very short and difficult to detect if the transition lies close to the hadronic maximum or if a disconnected branch also exists. By contrast, a disconnected branch can be conspicuous through gaps in observed radii or by collapse-like transitions accompanied by neutrino, gamma-ray, or gravitational-wave emission (Alford et al., 2013).

5. Representative realizations and Bayesian inference

The Hybrid EOS Framework has been implemented in several distinct microphysical realizations. Two-phase models with excluded volume in both sectors combine a stiffened DD2 hadronic EOS and a confining quark density functional in which the effective string tension is reduced by an available-volume factor. In these constructions, increasing hadronic excluded volume stiffens the hadronic branch and lowers the deconfinement chemical potential, while quark-sector screening and vector repulsion control onset and high-density stiffness; the resulting Maxwell-constructed hybrid stars can satisfy the $\mu_i = b_i \mu_B - q_i \mu_e.$6 constraint and develop disconnected third-family branches with high-mass twins (Kaltenborn et al., 2017, Alvarez-Castillo et al., 2015).

Other realizations replace the CSS or bag-like quark phase by microscopic models. A Dyson–Schwinger construction compared smooth crossover interpolations in $\mu_i = b_i \mu_B - q_i \mu_e.$7 and $\mu_i = b_i \mu_B - q_i \mu_e.$8, and found that interpolation in the $\mu_i = b_i \mu_B - q_i \mu_e.$9–$$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$0 plane can stiffen the crossover region enough to produce $$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$1 even when the underlying quark EOS is relatively soft (Zhao et al., 2015). Finite-temperature extensions based on DD2 plus a stiff color-superconducting quark sector use a two-zone parabolic interpolation in $$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$2, permit crossover or first-order behavior, and encode two critical endpoints through a temperature-dependent density jump; in that framework the 2SC quark phase has $$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$3 over $$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$4, and entropy-conserving trajectories are driven to higher temperatures than for deconfinement to normal quark matter (Ivanytskyi et al., 2022).

Bayesian analyses use the framework as an inference engine. In an APR+CSS Maxwell model with parameter vector $$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$5, sampled over $$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$6 MeV/fm$$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$7, $$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$8, and $$n_B(\mu_B) = \frac{\partial P}{\partial \mu_B},\qquad \epsilon = -P + \mu_B n_B + \sum_l \mu_l n_l,$$9, current disjunct mass–radius constraints did not decisively distinguish strong first-order transitions from purely hadronic EOSs; if a transition is present, the posterior favored $l=e,\mu$0 MeV/fm$l=e,\mu$1 and $l=e,\mu$2 MeV/fm$l=e,\mu$3 (Alvarez-Castillo et al., 2014). Under large-radius priors, a two-parameter excluded-volume plus 8-quark-vector model favored $l=e,\mu$4 and $l=e,\mu$5, with an optimum near $l=e,\mu$6, $l=e,\mu$7, and corresponding high-mass twins with radius differences of order $l=e,\mu$8 km (Alvarez-Castillo et al., 2015). In multimessenger analyses based on APR plus a color-superconducting nonlocal NJL quark sector connected by two-zone interpolation, the favored solutions required $l=e,\mu$9 and $$\frac{dP}{dr} = - \frac{G (\epsilon + P/c^2)\left(m + 4\pi r^3 P/c^2\right)}{r\left(r - 2Gm/c^2\right)},\qquad \frac{dm}{dr} = 4\pi r^2 \epsilon.$$0, and with the actual NICER radius of PSR J0740+6620 they no longer supported $$\frac{dP}{dr} = - \frac{G (\epsilon + P/c^2)\left(m + 4\pi r^3 P/c^2\right)}{r\left(r - 2Gm/c^2\right)},\qquad \frac{dm}{dr} = 4\pi r^2 \epsilon.$$1 (Ayriyan et al., 2021).

6. Assumptions, limits, and extensions

The framework’s predictive power is inseparable from its assumptions. Maxwell constructions require sufficiently large surface tension and enforce local neutrality in each phase; Gibbs constructions instead allow a globally neutral mixed phase with smoothly varying pressure, and crossover models remove discontinuities altogether. Finite-size effects, pasta structures, and multiple first-order transitions can modify the effective $$\frac{dP}{dr} = - \frac{G (\epsilon + P/c^2)\left(m + 4\pi r^3 P/c^2\right)}{r\left(r - 2Gm/c^2\right)},\qquad \frac{dm}{dr} = 4\pi r^2 \epsilon.$$2, change onset stability, and generate more complex branch structures than the single-transition CSS picture (Ji et al., 8 May 2025, Alford et al., 2013).

Microphysical incompleteness remains significant at intermediate density. Studies differ in whether they include strange quarks, hyperons, color superconductivity, excluded volume, multiquark interactions, or confinement-inspired functionals. This suggests that Hybrid EOS Framework is best regarded not as a single EOS, but as a controlled parameterized map from uncertain dense-QCD microphysics to observable compact-star phenomenology.

Independent of the detailed modeling, causality imposes global limits. For a maximally compact causal EOS, the maximum-mass configuration satisfies

$$\frac{dP}{dr} = - \frac{G (\epsilon + P/c^2)\left(m + 4\pi r^3 P/c^2\right)}{r\left(r - 2Gm/c^2\right)},\qquad \frac{dm}{dr} = 4\pi r^2 \epsilon.$$3

and the corresponding compactness is $$\frac{dP}{dr} = - \frac{G (\epsilon + P/c^2)\left(m + 4\pi r^3 P/c^2\right)}{r\left(r - 2Gm/c^2\right)},\qquad \frac{dm}{dr} = 4\pi r^2 \epsilon.$$4. If the largest observed neutron-star mass is the true maximum mass, model-independent upper bounds on central thermodynamic quantities follow; for $$\frac{dP}{dr} = - \frac{G (\epsilon + P/c^2)\left(m + 4\pi r^3 P/c^2\right)}{r\left(r - 2Gm/c^2\right)},\qquad \frac{dm}{dr} = 4\pi r^2 \epsilon.$$5, the quoted limits are

$$\frac{dP}{dr} = - \frac{G (\epsilon + P/c^2)\left(m + 4\pi r^3 P/c^2\right)}{r\left(r - 2Gm/c^2\right)},\qquad \frac{dm}{dr} = 4\pi r^2 \epsilon.$$6

These bounds provide an external consistency filter for any hybrid construction (Prakash, 2013).

Rotation, finite temperature, and composition effects broaden the framework further. Uniform rotation shifts mass–radius sequences but does not erase the basic branch topology; finite-temperature extensions are needed for supernovae and mergers, where thermal pressure, trapped neutrinos, and entropy gradients move the phase boundaries and can imprint the transition in post-merger gravitational-wave spectra (Ji et al., 8 May 2025, Ivanytskyi et al., 2022).

In this sense, the Hybrid EOS Framework is a layered program: a thermodynamic ansatz for dense matter, a dynamical prescription for phase conversion, a TOV-based map to compact-star observables, and an inference machinery that confronts $$\frac{dP}{dr} = - \frac{G (\epsilon + P/c^2)\left(m + 4\pi r^3 P/c^2\right)}{r\left(r - 2Gm/c^2\right)},\qquad \frac{dm}{dr} = 4\pi r^2 \epsilon.$$7 masses, NICER radii, and gravitational-wave deformabilities. Its central technical contribution is to reduce the hadron–quark transition to a small set of control parameters—transition location, latent heat, and post-transition stiffness—while retaining enough structure to classify stable branches, quantify maximum masses, and test specific dense-QCD scenarios against data.