pith. sign in

arxiv: 2607.00070 · v1 · pith:JBDKAWP6new · submitted 2026-06-30 · 🌌 astro-ph.HE · gr-qc· hep-ph

An Analytical Toy Equation of State for Neutron Stars Consistent with Current Observations

Tian-Shun Chen , Xiao-Ding Zhou , Kilar Zhang This is my paper

Pith reviewed 2026-07-02 18:00 UTC · model grok-4.3

classification 🌌 astro-ph.HE gr-qchep-ph
keywords neutron starequation of statepolytropicGW170817NICERtidal deformabilitymaximum masscausality
0
0 comments X

The pith

A smooth double-polytropic equation of state matches current neutron star observations while remaining causal.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The authors construct an analytic equation of state by scanning a double polytropic form for energy density in terms of pressure. They apply filters from massive pulsar masses, GW170817 tidal deformability, and NICER radius measurements to find viable parameter regions. A curve-integral score against public data sets selects benchmark models. These models achieve maximum masses of 2.44 to 2.49 solar masses, radii of about 11.3 km for 1.4 solar mass stars, and tidal deformabilities of 485 to 512. The family offers convenient reference cases for relativistic stellar structure calculations at current observational precision.

Core claim

Scanning the smooth double polytropic relation epsilon-hat(p-hat)=a1 p-hat^Gamma1 + a2 p-hat^Gamma2 identifies a viable region of parameters that passes all observational filters, with selected representatives supporting M_max = 2.44--2.49 solar masses, R_1.4 approx 11.3 km and Lambda_1.4=485--512 while remaining causal on the stable branch.

What carries the argument

The smooth double polytropic relation for energy density as a function of pressure, epsilon-hat(p-hat)=a1 p-hat^Gamma1 + a2 p-hat^Gamma2, which is scanned to find parameters compatible with observations.

Load-bearing premise

The double polytropic form is flexible enough to produce sequences compatible with all current constraints without violating causality or stability.

What would settle it

Detection of a neutron star with mass exceeding 2.5 solar masses or a 1.4 solar mass star with radius far from 11.3 km would rule out the benchmark models.

Figures

Figures reproduced from arXiv: 2607.00070 by Kilar Zhang, Tian-Shun Chen, Xiao-Ding Zhou.

Figure 1
Figure 1. Figure 1: One-term baseline feasibility scan in observable space. Panel (a) shows the R1.4-Λ1.4 projection; green and orange points enter the canonical radius–tidal window. Panel (b) projects the same canonical-window candidates into the Mmax-R2.08 plane; orange points also satisfy the adopted high-mass filters. Panel (c) shows the corresponding Mmax-max c 2 s/c2 projection. All orange points lie above the causal bo… view at source ↗
Figure 2
Figure 2. Figure 2: Double-polytrope feasibility scan in parameter space. Grey points are numerically converged filter rejections and green points pass all filters. Blue curves mark 50%, 68%, and 90% Gaussian KDE highest-density regions of the passing points. The second exchange-symmetric solution region is obtained by interchanging the two polytropic terms. guides. The representative curves pass through both the low-mass and… view at source ↗
Figure 3
Figure 3. Figure 3: Observable-space projection of the double-polytrope scan. Colors identify the filter status of each numerically converged model. Green points mark models that satisfy the adopted R1.4, Λ1.4, Mmax, R2.08, and stable-branch causality requirements simultaneously. 4×10 4 1×10 3 3×10 3 / 10 6 10 5 10 4 10 3 p/p H4 SLy MAP mean median [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Representative EOS pressure as a function of total energy density for the maximum-score, score-weighted median, and score-weighted mean models. The grey curves show the SLy and H4 Read et al. piecewise-polytrope fits. 4×10 4 1×10 3 3×10 3 / 0.0 0.2 0.4 0.6 0.8 1.0 1.2 c 2 s /c 2 H4 SLy MAP mean median [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Representative stable-branch sound speeds with the SLy and H4 Read et al. piecewise-polytrope fits shown in grey. The dashed line marks c 2 s/c2 = 1. imations, and sensitivity tests in which a tabulated microphysical EOS would make the EOS-dependence harder to isolate. The analytic form also gives direct access to dϵ/dp and the sound speed, making ther￾9 10 11 12 13 14 15 R [km] 0.8 1.0 1.2 1.4 1.6 1.8 2.0… view at source ↗
Figure 7
Figure 7. Figure 7: Representative dimensionless tidal deformability along the stable branch. The shaded region is the GW170817 component M-Λ highest-posterior-density projection from the public parametrized-EOS posterior samples. The ver￾tical bracket at 1.4 M⊙ marks the adopted Λ1.4 screening window; thin grey patterned curves are the same Read et al. reference fits as in [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
read the original abstract

Fast analytic and semi-analytic studies of neutron stars often require an equation of state that is convenient to evaluate while producing relativistic stellar sequences compatible with current multimessenger constraints. We construct such a benchmark by scanning a smooth double polytropic relation for the energy density as a function of pressure, $\hat\epsilon (\hat p)=a_1\hat p^{\Gamma_1}+a_2\hat p^{\Gamma_2}$. The parameters are selected with filters based on massive pulsars, tidal deformability from the binary-neutron-star event GW170817, and NICER mass-radius measurements. A single polytropic baseline scan finds no model passing all filters, whereas a double-polytrope scan identifies a viable region. A curve-integral score, evaluated against public NICER and GW170817 posterior data sets, is then used to choose benchmark equations of state within this region. The selected representatives support $M_{\max}=2.44$--$2.49\,M_\odot$, with $R_{1.4}\simeq 11.3$ km and $\Lambda_{1.4}=485$--$512$, and remain causal on the stable branch. This compact analytic family provides reference cases for relativistic stellar-structure tests at current observational scales.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript constructs an analytical toy equation of state for neutron stars via a smooth double-polytropic form ε̂(p̂)=a1 p̂^Γ1 + a2 p̂^Γ2. A parameter scan over a1, a2, Γ1, Γ2 applies observational filters drawn from massive pulsars, GW170817 tidal deformability, and NICER mass-radius data; a curve-integral score against public posteriors then selects benchmark models. These yield M_max=2.44–2.49 M_⊙, R_1.4≃11.3 km, Λ_1.4=485–512 while remaining causal on the stable branch. A single-polytrope scan produces no viable models.

Significance. If the scan is exhaustive and causality is enforced at every point on the stable branch, the resulting compact analytic family supplies convenient reference EOS for relativistic stellar-structure calculations at current observational scales. The explicit demonstration that a single polytrope is insufficient while the double-polytrope form succeeds, together with the use of public data sets for scoring, constitutes a reproducible and transparent construction method.

minor comments (3)
  1. Abstract: the statement that the selected representatives 'support' the quoted M_max, R_1.4 and Λ_1.4 values should be rephrased to make explicit that these properties are direct consequences of the observational filters and scoring procedure used for selection, rather than independent outcomes.
  2. The normalized variables ε̂ and p̂ and the precise normalization constants should be defined at first use, together with the pressure range over which the double-polytropic relation is applied and how causality is verified throughout the stable branch.
  3. The manuscript should state whether the curve-integral score is computed on the full posterior samples or on summary statistics, and whether any hold-out subset of the NICER or GW170817 data is reserved for validation.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the accurate summary of the manuscript, the positive significance assessment, and the recommendation for minor revision. No specific major comments were listed in the report.

Circularity Check

0 steps flagged

No significant circularity; explicit observational selection presented as such

full rationale

The paper's workflow consists of an explicit parameter scan of the double-polytropic form, application of filters drawn from external observational datasets (pulsars, GW170817, NICER), and selection of benchmark models that pass those filters. The reported M_max, R_1.4, and Lambda_1.4 values are direct outputs of this selection and are presented as reference cases consistent with the inputs, without any claim of independent first-principles derivation or prediction. No self-definitional equations, fitted inputs relabeled as predictions, or load-bearing self-citations appear in the described chain. The single-polytrope versus double-polytrope contrast is reported as an empirical scan result.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim depends on the choice of the double-polytrope functional form and the specific observational filters used to select viable parameters from the scan.

free parameters (1)
  • a1, a2, Gamma1, Gamma2
    Coefficients and exponents in the double polytrope ε̂(p̂)=a1 p̂^Γ1 + a2 p̂^Γ2, selected via observational filters.
axioms (2)
  • domain assumption The equation of state can be adequately represented by a smooth double polytropic function of pressure.
    Invoked in the construction of the benchmark EOS.
  • standard math General relativity governs the stellar structure and stability.
    Standard assumption for relativistic stellar models.

pith-pipeline@v0.9.1-grok · 5761 in / 1535 out tokens · 41442 ms · 2026-07-02T18:00:13.179831+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

  1. [1]

    GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral,

    Abbott, B. P., et al. 2017, Phys. Rev. Lett., 119, 161101, doi: 10.1103/PhysRevLett.119.161101 —. 2018, Phys. Rev. Lett., 121, 161101, doi: 10.1103/PhysRevLett.121.161101

    work page doi:10.1103/physrevlett.119.161101 2017

  2. [2]

    2020, Nat

    Vuorinen, A. 2020, Nat. Phys., 16, 907, doi: 10.1038/s41567-020-0914-9

    work page doi:10.1038/s41567-020-0914-9 2020

  3. [3]

    2013, Science, 340, 1233232, doi: 10.1126/science.1233232

    Antoniadis, J., et al. 2013, Science, 340, 1233232, doi: 10.1126/science.1233232

    work page doi:10.1126/science.1233232 2013

  4. [4]

    M., Thorne, K

    Bardeen, J. M., Thorne, K. S., & Meltzer, D. W. 1966, Astrophys. J., 145, 505, doi: 10.1086/148791

    work page doi:10.1086/148791 1966

  5. [5]

    Bertulani, C. A. 2023, Particles, 6, 198, doi: 10.3390/particles6010011

    work page doi:10.3390/particles6010011 2023

  6. [6]

    2024, Astrophys

    Choudhury, D., et al. 2024, Astrophys. J. Lett., 971, L20, doi: 10.3847/2041-8213/ad5a6f

    work page doi:10.3847/2041-8213/ad5a6f 2024

  7. [7]

    2009, Phys

    Damour, T., & Nagar, A. 2009, Phys. Rev. D, 80, 084035, doi: 10.1103/PhysRevD.80.084035

    work page doi:10.1103/physrevd.80.084035 2009

  8. [8]

    B., Pennucci, T., Ransom, S

    Demorest, P. B., Pennucci, T., Ransom, S. M., Roberts, M. S. E., & Hessels, J. W. T. 2010, Nature, 467, 1081, doi: 10.1038/nature09466

    work page doi:10.1038/nature09466 2010

  9. [9]

    J., Miller, M

    Dittmann, A. J., Miller, M. C., Lamb, F. K., et al. 2024, Astrophys. J., 974, 295, doi: 10.3847/1538-4357/ad5f1e

    work page doi:10.3847/1538-4357/ad5f1e 2024

  10. [10]

    T., Pennucci, T

    Fonseca, E., Cromartie, H. T., Pennucci, T. T., et al. 2021, Astrophys. J. Lett., 915, L12, doi: 10.3847/2041-8213/ac03b8

    work page doi:10.3847/2041-8213/ac03b8 2021

  11. [11]

    A., Gamba, R., Radice, D., & Bernuzzi, S

    Godzieba, D. A., Gamba, R., Radice, D., & Bernuzzi, S. 2021, Phys. Rev. D, 103, 063036, doi: 10.1103/PhysRevD.103.063036

    work page doi:10.1103/physrevd.103.063036 2021

  12. [12]

    Watts, A. L. 2019, Mon. Not. R. Astron. Soc., 485, 5363, doi: 10.1093/mnras/stz654

    work page doi:10.1093/mnras/stz654 2019

  13. [13]

    2001, Bernoulli, 7, 223, doi: 10.2307/3318737

    Haario, H., Saksman, E., & Tamminen, J. 2001, Bernoulli, 7, 223, doi: 10.2307/3318737

    work page doi:10.2307/3318737 2001

  14. [14]

    M., Pethick, C

    Hebeler, K., Lattimer, J. M., Pethick, C. J., & Schwenk, A. 2013, Astrophys. J., 773, 11, doi: 10.1088/0004-637X/773/1/11

    work page doi:10.1088/0004-637x/773/1/11 2013

  15. [15]

    2008, Astrophys

    Hinderer, T. 2008, Astrophys. J., 677, 1216, doi: 10.1086/533487

    work page doi:10.1086/533487 2008

  16. [16]

    2026, arXiv e-prints, arXiv:2603.12331, doi: 10.48550/arXiv.2603.12331

    Jarequi, G., Mitra, S., & Vaidya, V. 2026, arXiv e-prints, arXiv:2603.12331, doi: 10.48550/arXiv.2603.12331

    work page doi:10.48550/arxiv.2603.12331 2026

  17. [17]

    2019, Phys

    Kumar, B., & Landry, P. 2019, Phys. Rev. D, 99, 123026, doi: 10.1103/PhysRevD.99.123026

    work page doi:10.1103/physrevd.99.123026 2019

  18. [18]

    2019, Phys

    Landry, P., & Essick, R. 2019, Phys. Rev. D, 99, 084049, doi: 10.1103/PhysRevD.99.084049

    work page doi:10.1103/physrevd.99.084049 2019

  19. [19]

    Lattimer, J. M. 2012, Annu. Rev. Nucl. Part. Sci., 62, 485, doi: 10.1146/annurev-nucl-102711-095018

    work page doi:10.1146/annurev-nucl-102711-095018 2012

  20. [20]

    2022, Phys

    Landry, P. 2022, Phys. Rev. D, 105, 043016, doi: 10.1103/PhysRevD.105.043016

    work page doi:10.1103/physrevd.105.043016 2022

  21. [21]

    2010, Phys

    Lindblom, L. 2010, Phys. Rev. D, 82, 103011, doi: 10.1103/PhysRevD.82.103011 Analytical Toy EOS9 —. 2018, Phys. Rev. D, 97, 123019, doi: 10.1103/PhysRevD.97.123019

    work page doi:10.1103/physrevd.82.103011 2010

  22. [22]

    C., et al

    Miller, M. C., et al. 2019, Astrophys. J. Lett., 887, L24, doi: 10.3847/2041-8213/ab50c5 —. 2021, Astrophys. J. Lett., 918, L28, doi: 10.3847/2041-8213/ac089b O’Boyle, M. F., Markakis, C., Stergioulas, N., & Read, J. S. 2020, Phys. Rev. D, 102, 083027, doi: 10.1103/PhysRevD.102.083027

    work page doi:10.3847/2041-8213/ab50c5 2019

  23. [23]

    R., & Volkoff, G

    Oppenheimer, J. R., & Volkoff, G. M. 1939, Phys. Rev., 55, 374, doi: 10.1103/PhysRev.55.374 ¨Ozel, F., & Freire, P. 2016, Annu. Rev. Astron. Astrophys., 54, 401, doi: 10.1146/annurev-astro-081915-023322

    work page doi:10.1103/physrev.55.374 1939

  24. [24]

    Y., Fantina, A

    Potekhin, A. Y., Fantina, A. F., Chamel, N., Pearson, J. M., & Goriely, S. 2013, Astron. Astrophys., 560, A48, doi: 10.1051/0004-6361/201321697

    work page doi:10.1051/0004-6361/201321697 2013

  25. [25]

    E., & Watts, A

    Raaijmakers, G., Riley, T. E., & Watts, A. L. 2018, Mon. Not. R. Astron. Soc., 478, 2177, doi: 10.1093/mnras/sty1052

    work page doi:10.1093/mnras/sty1052 2018

  26. [26]

    S., Lackey, B

    Read, J. S., Lackey, B. D., Owen, B. J., & Friedman, J. L. 2009, Phys. Rev. D, 79, 124032, doi: 10.1103/PhysRevD.79.124032

    work page doi:10.1103/physrevd.79.124032 2009

  27. [27]

    E., et al

    Riley, T. E., et al. 2019, Astrophys. J. Lett., 887, L21, doi: 10.3847/2041-8213/ab481c —. 2021, Astrophys. J. Lett., 918, L27, doi: 10.3847/2041-8213/ac0a81

    work page doi:10.3847/2041-8213/ab481c 2019

  28. [28]

    J., Novak, J., Oertel, M., & Pons, J

    Servignat, G., Davis, P. J., Novak, J., Oertel, M., & Pons, J. A. 2024, Phys. Rev. D, 109, 103022, doi: 10.1103/PhysRevD.109.103022

    work page doi:10.1103/physrevd.109.103022 2024

  29. [29]

    W., Lattimer, J

    Steiner, A. W., Lattimer, J. M., & Brown, E. F. 2010, Astrophys. J., 722, 33, doi: 10.1088/0004-637X/722/1/33

    work page doi:10.1088/0004-637x/722/1/33 2010

  30. [30]

    L., & Providˆ encia, C

    Suleiman, L., Fortin, M., Zdunik, J. L., & Providˆ encia, C. 2022, Phys. Rev. C, 106, 035805, doi: 10.1103/PhysRevC.106.035805 ter Braak, C. J. F. 2006, Stat. Comput., 16, 239, doi: 10.1007/s11222-006-8769-1

    work page doi:10.1103/physrevc.106.035805 2022

  31. [31]

    Tolman, R. C. 1939, Phys. Rev., 55, 364, doi: 10.1103/PhysRev.55.364

    work page doi:10.1103/physrev.55.364 1939

  32. [32]

    Zhao, T., & Lattimer, J. M. 2018, Phys. Rev. D, 98, 063020, doi: 10.1103/PhysRevD.98.063020 10Chen, Zhou, and Zhang APPENDIX A.NUMERICAL SETTINGS The one-term baseline uses a uniform Sobol scan in Γ∈[0.20,1.40] and log 10 a∈[−3.00,2.00]. The double-polytrope scan used for the final feasible-region figures is uniform in the ordered box Γ1 ∈[0.10,0.60], Γ2 ...

    work page doi:10.1103/physrevd.98.063020 2018