Piecewise Polytropes in Neutron-Star EoS

• Piecewise polytropes are equation-of-state models that approximate dense matter by dividing it into segments each described by a polytropic relation.
• They support both forward-modeling of tabulated data and inverse reconstruction from observational data, enhancing neutron-star property estimation.
• Generalized formulations address discontinuities by enforcing sound speed continuity, improving simulation fidelity and astrophysical inference.

Searching arXiv for recent and foundational papers on piecewise polytropes to support the article. Piecewise polytropes (PP) are a class of equation-of-state (EoS) representations in which matter is modeled by adjoining density intervals, with each interval assigned a polytropic relation of the form $P = K_i \rho^{\Gamma_i}$. In neutron-star applications, PP representations are used to approximate tabulated cold dense-matter EoS, to generate broad ensembles of candidate EoS under uncertainty, and to support forward and inverse problems involving mass, radius, tidal deformability, moment of inertia, and mode spectra (Benitez et al., 2020). The formalism is valued because it is computationally tractable and flexible, but standard PP constructions do not, in general, preserve continuity of the sound speed across segment boundaries; this limitation has motivated several generalizations and revisions (O'Boyle et al., 2020).

1. Definition and basic formalism

In the standard PP construction, the high-density EoS is divided into $N$ intervals using prescribed dividing densities $\rho_1, \rho_2, \ldots, \rho_N$, and in each interval the pressure is written as

$$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$

The constant $K_i$ is chosen for continuity in pressure at the dividing densities, while the adiabatic index $\Gamma_i$ controls the stiffness of the segment (O'Boyle et al., 2020).

For cold neutron stars, the corresponding energy density for a given segment can be written as

$$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$

with the integration constant $a_i$ determined by demanding continuity at each segment boundary and with the initial value fixed by $\epsilon=\rho$ as $\rho\to 0$ (Suleiman et al., 2022).

A standard use case is to represent a tabulated EoS by a small number of polytropic regions. One formulation divides each unified EoS into seven polytropic segments—four for the crust and three for the core—with adaptive placement of transition densities to minimize the total fitting error (Suleiman et al., 2022). Another formulation, used for broad uncertainty studies, parameterizes the dense-matter EoS above a known low-density regime with five contiguous segments whose transition densities are set at $N$0, $N$1, $N$2, $N$3, and $N$4, with $N$5 (Benitez et al., 2020).

The PP formalism is therefore both a fitting device for known EoS tables and a generative parameterization for candidate EoS spaces. This dual role is central to its prominence in neutron-star phenomenology.

2. Segmentation strategies and parameter choices

Segment placement is a defining design choice in PP models. In one broad-sampling construction, the low-density regime $N$6 is taken to be known from laboratory experiments, with the QHC19 EoS used up to half nuclear saturation density, while the higher-density regime is represented by five contiguous polytropic segments (Benitez et al., 2020). Within each segment $N$7, the pressure is given by

$N$8

the polytropic index $N$9 is drawn uniformly from $\rho_1, \rho_2, \ldots, \rho_N$0, and $\rho_1, \rho_2, \ldots, \rho_N$1 is fixed to ensure pressure continuity at segment boundaries (Benitez et al., 2020).

In fitting modern unified neutron-star EoS, the segmentation can instead be chosen adaptively. A set of 52 unified equations of state was divided into seven polytropes via an adaptive segmentation, and two parameters per polytrope were fitted to the tabulated EoS (Suleiman et al., 2022). The pressure-density relation $\rho_1, \rho_2, \ldots, \rho_N$2 was interpolated from the tabulated EoS at 1500 points, with $\rho_1, \rho_2, \ldots, \rho_N$3 in the crust and $\rho_1, \rho_2, \ldots, \rho_N$4 in the core, logarithmically distributed for optimal accuracy across relevant density regimes (Suleiman et al., 2022). The fitting algorithm iteratively proposed initial guesses for transition densities, fitted each segment with $\rho_1, \rho_2, \ldots, \rho_N$5 and $\rho_1, \rho_2, \ldots, \rho_N$6 using nonlinear least squares, enforced pressure continuity, refined transition densities, and repeated until convergence (Suleiman et al., 2022).

A related development is the search for preferred dividing densities. In the generalized piecewise-polytropic framework, preferred dividing densities are determined by minimizing an error norm consisting of integral astrophysical observables rather than only pointwise thermodynamic quantities (O'Boyle et al., 2020). For the core-candidate EoS used by LIGO/Virgo, single universal values for the dividing densities were found: $\rho_1, \rho_2, \ldots, \rho_N$7 This suggests that segmentation can be optimized relative to observable fidelity rather than imposed a priori (O'Boyle et al., 2020).

3. Continuity, differentiability, and generalized piecewise polytropes

A major limitation of the standard PP parameterization is that the sound speed $\rho_1, \rho_2, \ldots, \rho_N$8 is not generally continuous at dividing densities, because the adiabatic index jumps from one segment to the next (O'Boyle et al., 2020). The stated consequences include unphysical behavior in quantities sensitive to second derivatives, artifacts in simulations near junctions, and non-differentiable observable mappings (O'Boyle et al., 2020).

The generalized piecewise polytrope (GPP) addresses this by imposing continuity in pressure, energy density, and sound speed. In each segment, the EoS may be written as

$\rho_1, \rho_2, \ldots, \rho_N$9

$$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$0

where the offset $$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$1 is the additional integration constant that is absent in standard PP constructions (O'Boyle et al., 2020). At each dividing density $$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$2, the GPP imposes simultaneous continuity of pressure, pressure derivative, energy density, and energy-density derivative. As a result, both pressure and its first derivative—directly related to the sound speed,

$$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$3

—are continuous across every dividing density (O'Boyle et al., 2020).

The practical implications described for GPPs are specific. Generalized piecewise polytropes accurately reproduce astrophysical observables such as mass, radius, tidal deformability and mode frequencies, as well as thermodynamic quantities such as the adiabatic index (O'Boyle et al., 2020). The formalism is algebraically simple and fully differentiable, can improve pointwise convergence in numerical relativity simulations of neutron stars, and can improve Bayesian parameter estimation from gravitational waveforms (O'Boyle et al., 2020). Existing implementations of piecewise polytropes can easily accommodate this generalization, and the overall number of free parameters is unchanged (O'Boyle et al., 2020).

A further point is that the GPP structure allows a controlled discontinuity in sound speed to be introduced, if desired, to model phase transitions by relaxing continuity at a chosen dividing density (O'Boyle et al., 2020). This makes continuity an adjustable modeling choice rather than an unavoidable structural defect.

4. Piecewise polytropes in neutron-star universal relations

PP ensembles have been used extensively to study approximate EoS-independence in neutron-star observables. One such application modeled dense matter with piecewise polytropic representations with five segments and used them as EoS input for solving the stellar structure equations and calculating the moment of inertia $$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$4, tidal deformability $$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$5, rotational mass quadrupole moment $$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$6, and axial $$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$7-mode frequencies (Benitez et al., 2020).

The associated dimensionless quantities were defined as

$$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$8

with $$p = K_i \rho^{\Gamma_i}, \qquad \rho_{i-1} < \rho < \rho_i .$$9 and $K_i$0 (Benitez et al., 2020). Relations among $K_i$1, $K_i$2, $K_i$3, and the real and imaginary parts of the scaled $K_i$4-mode frequency were fitted by a fourth-order polynomial in $K_i$5,

$K_i$6

Here $K_i$7 denotes $K_i$8, $K_i$9, $\Gamma_i$0, or $\Gamma_i$1 (Benitez et al., 2020).

The reported results are two-tiered. On the first stable branch, the I-Love-Q relations are tight, with root-mean-square residuals at the $\Gamma_i$2–$\Gamma_i$3 level, and the $\Gamma_i$4-$\Gamma_i$5 relation is tighter than the $\Gamma_i$6-$\Gamma_i$7 relation (Benitez et al., 2020). On the second branch associated with twin stars, the relations are less tight by a factor of $\Gamma_i$8, and branch-specific fitting improves modeling compared with applying first-branch fits to second-branch stars (Benitez et al., 2020). The paper also proposes new empirical relations between $\Gamma_i$9, $$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$0, $$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$1, and the complex frequency $$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$2 of the fundamental axial $$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$3-mode and finds that they are comparably tight to the I-Love-Q correlations (Benitez et al., 2020).

These results are presented as evidence that the PP formalism is flexible yet robust for EoS uncertainty quantification. The largest outliers are associated with extreme choices of $$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$4, including a maximum in the low-density segment followed by sharp drops, indicating that some polytropic sequences create more diverse structural features (Benitez et al., 2020).

5. Reconstruction, fitting accuracy, and observational constraints

PP parameterizations are also used in inverse problems. A piecewise polytropic meshing and refinement method reconstructs the neutron-star EoS from experimental data of the mass and the tidal Love parameter using a 4-dimensional parameter space

$$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$5

with an initial mesh of $$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$6 equations of state (Mena-Fernández et al., 2019). The initial mesh is defined by $$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$7 equally spaced points in each parameter direction: $$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$8 For each EOS, the method solves the Tolman-Oppenheimer-Volkoff equations and the tidal Love number equation, computes the $$\epsilon(\rho) = (1 + a_i)\rho + \frac{\kappa_i}{\Gamma_i - 1}\rho^{\Gamma_i},$$9 curve, evaluates an error metric against the input data, and refines the mesh around the best candidates until a target tolerance is met (Mena-Fernández et al., 2019).

In a test using six $a_i$0 configurations from APR4, the best candidate after two local refinements reached $a_i$1, and the relative difference for all the parameters was smaller than $a_i$2 (Mena-Fernández et al., 2019). Reported maximum errors over 20 sequences between the original and reconstructed APR4 EOS include radius $a_i$3, mass $a_i$4, moment of inertia and quadrupole $a_i$5, tidal deformability $a_i$6, quasinormal mode frequencies $a_i$7, quasinormal mode damping times $a_i$8, and scaled $a_i$9-mode frequencies $\epsilon=\rho$0 (Mena-Fernández et al., 2019).

PP fits to tabulated unified EoS provide a complementary forward-modeling perspective. For 52 unified EoS, the total mass, radius, tidal deformability and moment of inertia of neutron stars were modeled from seven-segment PP fits and compared with the quantities calculated from the original tables (Suleiman et al., 2022). The reported fit errors are small: maximum mass is typically reproduced to within $\epsilon=\rho$1, radii at $\epsilon=\rho$2 and $\epsilon=\rho$3 are reproduced with $\epsilon=\rho$4 error for all 51 EoS with $\epsilon=\rho$5, the moment of inertia error is $\epsilon=\rho$6, and the tidal deformability is usually $\epsilon=\rho$7 for the vast majority, peaking up to $\epsilon=\rho$8 for a few outliers (Suleiman et al., 2022).

The observational role of PP parameter spaces is illustrated by GW170817. Using the TaylorF2 waveform model for the low-spin scenario, EOSs with $\epsilon=\rho$9 are marked as excluded at $\rho\to 0$0 confidence in the equal-mass case, and the excluded region in parameter space is described as largely independent of $\rho\to 0$1 (Mena-Fernández et al., 2019). A plausible implication is that, within this specific parameterization, current gravitational-wave constraints predominantly act on lower-to-intermediate density segments.

6. Extensions, physical interpretation, and limitations

Several works place PP methods within broader polytropic modeling. One development identifies the polytrope index with the pressure derivative of the bulk modulus,

$\rho\to 0$2

where $\rho\to 0$3 is the bulk modulus and $\rho\to 0$4 is the compressibility (Weppner et al., 2014). In that framework, traditional PP corresponds to assigning piecewise-constant $\rho\to 0$5, while a variable-index formulation allows $\rho\to 0$6 to vary smoothly with $\rho\to 0$7 or $\rho\to 0$8, recovering PP as a limiting case (Weppner et al., 2014). The generalized Lane-Emden equation with variable index then provides a route from PP-style discretizations to smoothly varying material models (Weppner et al., 2014). This suggests that PP can be interpreted as a coarse-grained approximation to a more continuous elastic response.

Another extension concerns multi-component, rigidly rotating polytropes. The theory has been improved to allow subsystems with nonzero density on the boundary and subsystems with intersecting boundaries, with equilibrium configurations independent of whether the fluid is collisional or collisionless provided the polytropic index lies within $\rho\to 0$9 (Caimmi, 2016). The formalism is restricted to two subsystems for simplicity but can, in principle, be extended to $N$00 subsystems (Caimmi, 2016). The paper’s summary explicitly notes that this is directly applicable to piecewise polytropes where segments meet at interfaces with generally nonzero density (Caimmi, 2016).

The principal limitations of standard PP are also explicit in the literature summarized here. First, the sound speed is not continuously defined at segment boundaries, which may limit use in studies where high-order derivatives of the EoS are needed (Suleiman et al., 2022). Second, seven segments suffice for current needs in one comparative study, but exotic EoS with sharp phase transitions or highly nontrivial structure may require more (Suleiman et al., 2022). Third, using nonunified EoS by gluing crusts and cores from different microphysics can create significant artificial errors in macroscopic neutron-star properties, whereas unified fits reduce these errors (Suleiman et al., 2022). Finally, broad PP parameter surveys may include extreme choices of $N$01 that generate physically less likely or more diverse structural features (Benitez et al., 2020).

Taken together, these developments present piecewise polytropes as a hierarchy rather than a single recipe: standard PP for efficient parametrization, generalized PP for continuity and differentiability, adaptive unified fits for high-accuracy reproduction of tabulated EoS, and inverse PP meshes for observational reconstruction. Across these uses, the method’s core premise remains unchanged: approximate the EoS by adjoining polytropic segments while enforcing the continuity conditions appropriate to the intended application (O'Boyle et al., 2020).