How Neutron Star Radii Encode the Dense-Matter Equation of State and Hadron-Quark Transition

Technical Summary: How Neutron Star Radii Encode the Dense-Matter Equation of State and Hadron-Quark Transition

Problem Statement
A central goal of modern neutron star (NS) physics is to infer the microscopic properties of supradense matter from macroscopic observations. While Bayesian inference has become the standard framework for extracting Equation of State (EOS) constraints from NS data, existing methods often provide limited insight into the physical mechanisms linking microscopic EOS parameters to observable stellar radii. Specifically, it remains unclear why improving the precision of radius measurements dramatically shifts the inferred values of certain EOS parameters while leaving others largely unaffected. This question is particularly critical when considering first-order hadron–quark phase transitions, which can produce "twin stars" (stars with identical masses but different radii) and result in distinct mass–radius topologies (Connected, Disconnected, Both, or No-Quark-Matter). The paper addresses the open question of how microscopic dense-matter physics is encoded into measurable NS radii and how this encoding depends on the global topology of the mass–radius sequence.

Methodology
The authors employ a Bayesian framework using a meta-model EOS that encapsulates a first-order hadron–quark phase transition via the Constant Speed of Sound (CSS) model. The hadronic sector is parameterized by nine microscopic parameters: the incompressibility ($K_0$), skewness ($J_0$), symmetry energy slope ($L$), curvature ($K_{sym}$), skewness ($J_{sym}$), and symmetry energy at saturation ($E_{sym}(\rho_0)$), alongside transition parameters: transition density ($\rho_t$), normalized energy-density jump ($\Delta\epsilon/\epsilon_t$), and quark-matter sound speed ($c_s^2/c^2$).

The study utilizes mock measurements of the canonical NS radius ($R_{1.4}$) with a central value of 11.9 km and two uncertainty levels: $\sigma_R = 0.9$ km (representing current precision) and $\sigma_R = 0.1$ km (representing expected future precision). Accepted EOSs are classified into four topological categories based on the structure of their mass–radius sequences:

1. Connected: The hybrid branch connects directly to the hadronic branch.
2. Disconnected: The hybrid branch is separated from the hadronic branch by an unstable region.
3. Both: Both connected and disconnected branches coexist.
4. No-Quark-Matter: Purely hadronic stars.

The core methodological innovation is the introduction of inverse EOS–radius mappings. Instead of solely analyzing posterior probability distributions (PDFs), the authors calculate the posterior mean of each EOS parameter $\langle \theta_i \rangle$ as a function of the inferred radius $R_{1.4}$, marginalizing over all other parameters. These mappings are analyzed in conjunction with the Jensen expansion, which relates the curvature of the inverse mapping to the systematic shift in the posterior mean as observational precision improves.

Key Contributions

1. Inverse Mappings as Physical Probes: The paper introduces inverse mappings as a tool to visualize how information in radius measurements is transmitted back to the underlying EOS. The slope of the mapping quantifies the sensitivity of a parameter to radius variations, while the curvature determines the precision-dependent shift in the posterior mean.
2. Topology-Resolved Analysis: By separating the analysis into the four mass–radius topologies, the authors reveal a hierarchy of information that is obscured when all EOSs are treated as a single ensemble.
3. Decoupling Precision and Topology Dependence: The study distinguishes between parameters that are primarily constrained by the value of the radius (precision-dependent) and those that are constrained by the global structure of the mass–radius sequence (topology-dependent).

Key Results

• Hadronic Parameters:

  • Symmetry Energy ($L$ and $K_{sym}$): These parameters are strongly encoded in $R_{1.4}$. Their inverse mappings exhibit substantial slopes, indicating high sensitivity to the radius. Consequently, their posterior means shift appreciably as radius precision improves from 0.9 km to 0.1 km. These parameters are largely topology-independent; different topologies produce similar $L$–$R_{1.4}$ and $K_{sym}$–$R_{1.4}$ relations because $R_{1.4}$ is primarily determined by the pressure around $2\rho_0$.
  • High-Density Parameters ($J_0$ and $J_{sym}$): These parameters show weaker dependence on $R_{1.4}$ and exhibit pronounced topology dependence. Their inverse mappings vary significantly across topological categories, reflecting the fact that different high-density continuations of the hadronic EOS can yield the same canonical radius.

• Transition Parameters:

  • Transition Density ($\rho_t$): This is the most strongly encoded transition parameter in $R_{1.4}$. The inverse mapping shows a non-monotonic but significant dependence on radius. Improving radius precision leads to better localization of $\rho_t$ and systematic shifts in its posterior mean, particularly for Connected and Both categories.
  • Energy-Density Jump ($\Delta\epsilon/\epsilon_t$) and Sound Speed ($c_s^2$): These parameters show weak dependence on radius precision but strong separation among topological classes. Large energy-density jumps and low sound speeds favor the Disconnected topology, while smaller jumps and higher sound speeds favor the Connected topology. The posterior means of these parameters change little with improved radius precision within a given topology, indicating that $R_{1.4}$ alone cannot uniquely determine the strength of the transition or the stiffness of the quark phase.

• Topology Identification: Despite the strong separation in parameter space for different topologies, the posterior distributions of $R_{1.4}$ for the four categories overlap significantly. Therefore, even precise radius measurements ($\sigma_R = 0.1$ km) cannot by themselves identify the mass–radius topology or uniquely determine the high-density transition properties.

Significance and Claims
The paper claims that the scientific return of future high-precision radius measurements is intrinsically parameter-dependent.

• What is Achievable: High-precision $R_{1.4}$ measurements are highly effective for constraining the symmetry-energy sector ($L, K_{sym}$) governing the pressure around $2\rho_0$ and for providing useful constraints on the transition density ($\rho_t$).
• What is Not Achievable: Radius measurements alone cannot uniquely determine the nature of the hadron–quark transition (specifically the energy-density jump and sound speed) or identify the global mass–radius topology (Connected vs. Disconnected).
• Conclusion: The authors conclude that a complete understanding of the dense-matter EOS requires complementary probes. While precise radii constrain the EOS leading up to the transition, additional observations sensitive to the global high-density stellar structure are necessary to distinguish topologies and constrain the physical nature of the phase transition itself. The inverse mappings provide a unified physical interpretation for why different parameters respond differently to observational improvements.