REVIEW 3 major objections 4 minor 1 cited by
Fully general relativistic description of rapidly-rotating axially-symmetric neutron stars for constraining nuclear matter equations of state
T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Using fully general-relativistic rotating-star models, this paper claims that neutron stars can spin at 716 Hz in mass–radius regions that the standard empirical Kepler-limit formula marks as forbidden, so the 716 Hz constraint on…
desk verdict Careful KEH implementation with a genuinely interesting but overstated challenge to the 716 Hz constraint; the stability language goes beyond what is computed. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the KEH method, a fully general-relativistic self-consistent field scheme that constructs rigidly rotating, axially symmetric perfect-fluid stars. It solves elliptic integral equations for the metric potentials $\gamma$, $\rho$, $\alpha$ and the frame-dragging frequency $\omega$ in a quasi-isotropic coordinate system, together with the first integral of hydrostatic equilibrium that ties the fluid to the equation of state; a compactified radial coordinate brings spatial infinity onto the grid. The Kepler (mass-shedding) limit is evaluated from the geodesic condition for equatorial circular orbits, and the paper compares that general-relativistic limit against the Newtonian empirical formula. The other essential ingredient is a tabulated npeμ equation of state built from five Skyrme-type parameter sets plus standard crust equations of state, which lets the code output not only masses and radii but internal density and composition profiles.
What would settle it
Evolve one of the reported low-mass 716 Hz equilibria (for example the SLy4 configuration with central density $n_c \simeq 0.54\,n_0$ and mass about $0.75\,M_\odot$) with a fully general-relativistic hydrodynamics code: if it sheds mass, collapses, or develops an r-mode instability within a short dynamical time, the claim that this region is stable is falsified. A less expensive check is a linear radial-oscillation or r-mode analysis of the same equilibria looking for unstable modes.
Extended reading notes
Core claim
On its own terms, the paper claims that a fully general-relativistic treatment of axisymmetric rotating stars changes how the 716 Hz spin constraint should be read. Using the KEH method, a fully general-relativistic self-consistent field scheme, with core matter described by five Skyrme-type equations of state in npeμ (neutron–proton–electron–muon) composition and standard crust equations of state, it computes mass–radius curves at fixed spin frequencies of 0, 205, 346, and 716 Hz. At 716 Hz, a significant portion of the gray 'forbidden' region defined by the empirical formula $\nu_K = 1045\,(M/1.4M_\odot)^{1/2}(10\,\mathrm{km}/R)^{3/2}\,\mathrm{Hz}$ is populated by convergent solutions whose angular frequency stays below the general-relativistic Kepler limit; the paper concludes that the parameter space accessible to stable rotating configurations is broader than previously estimated and that the 716 Hz constraint should be interpreted more conservatively. It also finds that the rotational increase in maximum mass at 346 Hz can push an EoS like SLy4 over the $2\,M_\odot$ boundary, and that radius enhancement at a given spin scales with the slope of the symmetry energy, $L$.
Load-bearing premise
The word 'stable' here only means the computed spin is below the mass-shedding limit; the paper does not check whether the 716 Hz configurations survive perturbations, and if any of them are dynamically unstable the claim of a broader allowed region is weakened.
Editorial extensions
If this is right
- The empirical 716 Hz limit curve should be treated as approximate: fully general-relativistic equilibria occupy a substantial part of the mass–radius region the curve labels forbidden, so the constraint is less restrictive than usually stated.
- For spin frequencies above about 400 Hz, rotation changes mass and radius enough that comparing TOV models with observations can bias equation-of-state constraints; rotating models become necessary.
- At 346 Hz the rotational rise in maximum mass, though small, can matter for marginal cases: for SLy4 the maximum mass increases from $2.055$ to $2.062\,M_\odot$, moving it from barely consistent to more comfortably consistent with the observed mass of PSR J0740+6620.
- The rotational increase in radius at fixed mass correlates with the symmetry-energy slope $L$, so equations of state with similar incompressibility but different $L$ respond differently to the same spin.
- Low-mass stars deform more at a given frequency and develop thicker equatorial crusts, making low-mass rapidly rotating stars the most sensitive probes of rotation-induced structural changes.
Reading between the lines
- Editorial inference: if the 716 Hz equilibria pass a full dynamical-stability test, pulsar searches should not dismiss low-mass candidates that fall inside the empirically 'forbidden' band, since such objects could be genuinely stable millisecond pulsars.
- Editorial inference: a measured mass and radius for any pulsar near 716 Hz would turn the spin frequency from a veto into a direct EoS probe, because only equations of state stiff enough in the right density range can support that spin in full general relativity.
- Editorial inference: extending the calculation to differential rotation, which the paper notes can support larger maximum masses, would likely shift the allowed region further outward, but it would also make the missing dynamical-stability analysis more urgent.
- Editorial inference: a controlled scan that varies a single saturation parameter such as $L$ at a time, as the paper suggests for future work, could convert the observed correlations into a quantitative map from nuclear parameters to maximum spin frequency.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a newly developed fully general relativistic code based on the KEH method for computing equilibrium structures of rapidly rotating, axially symmetric neutron stars. Using five Skyrme-type EoS parameter sets with npeμ core matter and BPS/BBP crusts, the authors compute mass-radius relations at fixed spin frequencies, study how rotation shifts these relations, correlate the shifts with nuclear matter parameters (especially the symmetry-energy slope L), and compare with the spin frequencies of PSR J0740+6620, PSR J0030+0451, and PSR J1748-2446ad. The central claim, stated in Sec. IV.C.3 and Fig. 7, is that fully general relativistic rotating configurations occupy a significant portion of the parameter space that the empirical 716 Hz constraint of Eq. (44) marks as forbidden, so that this empirical constraint should be interpreted more conservatively.
Significance. If the central claim holds, the paper provides a useful, fully relativistic reassessment of a widely used observational constraint on neutron star equations of state, and it offers concrete M-R curves for several modern Skyrme EoSs at realistic spin frequencies. The numerical work is self-contained: the KEH inputs are central density and axis ratio, the spin frequency is an output, the EoS parameters come from independent prior fits, and no parameter is fitted to the 716 Hz constraint. The convergence study in Table II shows mass and radius stable to about the third decimal across resolutions from 257×257 to 2049×2049, which supports the internal consistency of the code. However, the paper's main conclusion is weakened by the fact that 'stability' is identified only with being below the Kepler mass-shedding limit, and by the absence of any quantitative external validation of the code in the manuscript.
major comments (3)
- [Sec. IV.C.3, Fig. 7, Eq. (42)] The claim that a significant portion of the empirically forbidden region 'can actually correspond to stable rotating neutron star configurations' is not supported by the calculations shown. The only stability criterion applied is that the computed angular frequency lies below the general relativistic Kepler (mass-shedding) limit of Eq. (42). The KEH method solves for equilibrium configurations, not for dynamical, secular, or non-axisymmetric stability; radial/turning-point stability, the bar-mode instability, and the CFS r-mode instability are not examined. The manuscript itself states in Sec. V that 'stability analysis becomes crucial, particularly in connection with phenomena such as r-mode instabilities' and defers it to future work. Because the paper's central conclusion about the 716 Hz constraint depends on these equilibria being physically realizable and not merely formal solutions, the word 'stable' should either be replaced with 'equilibrium' or the statement should be accompanied by an explicit stability analysis or a clear caveat that stability has not been established.
- [Sec. III] The claimed validation of the new code against Stergioulas and Friedman (Ref. [18]) and Morsink and Stella (Ref. [41]) is only asserted verbally, without any quantitative comparison. Since the central conclusion rests on the accuracy of the KEH implementation and the code is not released, please include a table comparing masses, equatorial radii, and Kepler frequencies for identical input EoS, central density, and axis ratio against the published results, or at least provide a quantitative convergence statement for the Kepler frequency itself. The current convergence table (Table II) is for one configuration and does not demonstrate that the code reproduces established rotating-star results to the accuracy needed for the 716 Hz claim.
- [Sec. V] The discussion in Sec. V acknowledges that 'in order to systematically explore the 716 Hz constraint across a wider range of EoSs—especially to identify the lowest-mass neutron star that can stably rotate at this frequency—it will be necessary to develop a more numerically stable and automated framework.' This admission undercuts the strength of the Sec. IV.C.3 conclusion, because the paper does not actually identify the boundary of the allowed region or demonstrate that the configurations shown at 716 Hz are stable. Please either temper the central claim accordingly or provide the additional analysis needed to support it.
minor comments (4)
- [Abstract and Sec. II.D.2] The abstract describes the KEH method as providing 'stable solutions' for rotating equilibrium configurations; this wording conflates equilibrium with stability. Please change it to 'equilibrium configurations' or explicitly qualify that stability is meant only with respect to mass shedding.
- [Sec. V] There are several typographical and grammatical errors, e.g., 'relays on' should be 'relies on' and 'metioned' should be 'mentioned'. A careful proofread is recommended.
- [Table I and Ref. [25]] The SkT5a parameter set is cited only as 'M. Duta, private communication'. Please provide a citable published source for this parameter set, since it is used in several quantitative conclusions.
- [Sec. IV.C.1] The statement that rotation at 346 Hz makes the SLy4 EoS 'more satisfactory' relative to the PSR J0740+6620 mass constraint is based on a single EoS and a small mass increase (ΔM ≈ 0.007 M⊙). It would be helpful to state explicitly that the shift is only marginally larger than the numerical uncertainty implied by Table II, and that this is not a decisive EoS test.
Circularity Check
No significant circularity: the central results are self-contained and validated against independent external benchmarks.
full rationale
The paper's derivation chain is self-contained. Rotating neutron star structures are computed with a newly developed code implementing the KEH method, where the inputs are the central density and the axis ratio, and the angular frequency, mass, and radius are self-consistently obtained outputs. The five Skyrme EoS parameter sets (BSk24, SkI4, SkM*, SkT5a, SLy4) are taken from independent prior fits to nuclear data, and the paper verifies that their saturation properties match the original references. The code is validated against published external results: stellar structures agree with Stergioulas and Friedman and Kepler frequencies are consistent with Morsink and Stella. The central claim about the 716 Hz constraint compares the computed Kepler-limited configurations against the empirical Lattimer-Prakash formula, which is an external benchmark and not a quantity fitted in this work. No parameter is fitted to the target result, and the maximum angular frequencies are not defined in terms of the empirical constraint. The paper's use of the word 'stable' relies on the general relativistic Kepler (mass-shedding) limit rather than a full dynamical stability analysis; this is a scientific caveat, not a circularity, because the equilibrium solutions are obtained independently and the stability criterion is not an input to the calculation. Self-citations are limited to methodological references (KEH, CST, RNS) that are independently established and externally validated by the community. Therefore, no load-bearing step reduces to its own input or to a self-citation chain.
Assumptions & free parameters
free parameters (1)
- A (differential rotation parameter) =
∞ (rigid rotation)
assumptions (5)
- domain assumption Rigid-body rotation with A tending to infinity
- domain assumption npeμ composition throughout the core, no exotic degrees of freedom
- domain assumption BPS and BBP crust EoS adopted
- standard math Validity of the KEH formulation
- domain assumption Skyrme EoS parameters from published fits
Cite this review
Pith. "Pith review of Fully general relativistic description of rapidly-rotating axially-symmetric neutron stars for constraining nuclear matter equations of state." pith.science (2026). https://pith.science/paper/QQ4FC77F
@misc{pith2026250520990,
author = {Pith},
title = {Pith review of: Fully general relativistic description of rapidly-rotating axially-symmetric neutron stars for constraining nuclear matter equations of state},
year = {2026},
howpublished = {\url{https://pith.science/paper/QQ4FC77F}},
note = {Machine review of arXiv:2505.20990}
}
abstract
Background: Constraining the nuclear matter equation of state (EoS) from neutron star observations is one of the main subjects in nuclear physics today. In general, neutron stars rotate rapidly and structure of neutron stars can be affected, especially in millisecond pulsars. To better constrain the nuclear EoS, it is important to describe neutron star structure taking into account the effects of rotation in a fully relativistic manner. Purpose: In this study, we investigate the internal structure of neutron stars under the influence of rotation. We explore correlations between rotational effects and EoS parameters, based on realistic calculations of rapidly rotating neutron stars based on the KEH method, which provides stable solutions for axially symmetric rotating equilibrium configurations. Results: Using 5 different Skyrme EoS parameter sets, we find that the maximum angular frequency achievable by rotating neutron stars, as calculated via the KEH method, varies depending on the stiffness of the equation of state. We confirm that an increase in the rotating frequency leads to an overall increase in both the mass and radius along the M-R curve. By performing calculations at two frequently referenced neutron stars, we further examine how the changes in mass and radius correlate with the nuclear matter properties at saturation density. Our results suggest that the 716Hz rotational constraint may require a more conservative interpretation when accounting for realistic stellar deformation effects. Conclusions: To place stringent constraints on the nuclear EoS based on observational data, it is sometimes essential to account for the effects of rotation in neutron star models. In particular, the influence of rotation becomes increasingly significant at higher spin frequencies and cannot be neglected in rapidly rotating systems with $>$400Hz.
Figures
Figures from the paper (4 more)
Forward citations
Cited by 1 Pith paper
-
Accuracy and Applicability of the Hartle-Thorne and Komatsu-Eriguchi-Hachisu Methods for Modeling Rotating Neutron Stars
For OMEG equations of state, Hartle-Thorne and fully relativistic KEH rotation models agree closely below ~200 Hz but differ by up to ~27% in deformation at 800 Hz, with larger discrepancies for stiffer symmetry-energ...
Reference graph
Works this paper leans on
- [18]
-
[41]
G. Baym, C. Pethick, and P. Sutherland, ApJ170, 299 (1971)
work page 1971
-
[1]
Static case To calculate the structure of a neutron star in hydrostatic equilibrium within general relativity, we begin with the Ein- stein field equations, which relate the curvature of spacetime to the energy-momentum content of the star: Gµν =R µν − 1 2 gµνR= 8πG c4 Tµν,(26) whereG µν is the Einstein tensor,R µν is the Ricci tensor,g µν is the metric t...
-
[2]
Rotating case To handle such rotating configurations, it is convenient to adopt a quasi-isotropic coordinate system. In this system, the line element takes the following form: ds2 =−e γ+ϱdt2+e2α(dr2+r2dθ2)+eγ−ϱr2 sin2 θ(dϕ−ωdt)2, (31) whereγ(r, θ),ϱ(r, θ), andα(r, θ)are metric potentials, and ω(r, θ)is the angular velocity of frame dragging. In the KEH me...
-
[3]
This pulsar is also known to be rotating at a spin frequency of 346 Hz
PSR J0740+6620 The mass constraint from PSR J0740+6620 is one of the most stringent observational constraints currently known, in- dicating a neutron star mass exceeding2M ⊙. This pulsar is also known to be rotating at a spin frequency of 346 Hz. In Fig. 5, we present theM-Rcurves for rotating neutron stars computed with each EoS parameter set at this ang...
-
[4]
Observational anal- ysis of PSR J0030+0451 is based on two Bayesian studies TABLE V
PSR J0030+0451 We have also performed calculations at a spin frequency of 205 Hz, corresponding to PSR J0030+0451, which is of- ten referenced in connection with radius constraints near the canonical neutron star mass of1.4M ⊙. Observational anal- ysis of PSR J0030+0451 is based on two Bayesian studies TABLE V . Comparison of neutron star properties from ...
-
[5]
PSR J1748–2446ad Lastly, we consider the so-called Keplerian constraint de- rived from PSR J1748–2446ad, which has the highest known spin frequency of 716 Hz. Lattimeret al.[45] suggested an empirical upper limit on the rotation frequency for rigid New- tonian spheres, which is given by νK ≈1045 M 1.4M ⊙ 1/2 10 km R 3/2 Hz.(44) Using this expression, one ...
- [6]
Show all 50 references
-
[7]
Dutra, O
M. Dutra, O. Lourenço, S. S. Avancini, B. V . Carlson, A. Delfino, D. P. Menezes, C. Providência, S. Typel, and J. R. Stone, Phys. Rev. C90, 055203 (2014)
2014
-
[8]
B. Sun, S. Bhattiprolu, and J. M. Lattimer, Phys. Rev. C109, 055801 (2024)
2024
-
[9]
J. R. Oppenheimer and G. M. V olkoff, Phys. Rev.55, 374 (1939)
1939
-
[10]
Einstein, Annalen der Physik354, 769 (1916), https://onlinelibrary.wiley.com/doi/pdf/10.1002/andp.19163540702
A. Einstein, Annalen der Physik354, 769 (1916), https://onlinelibrary.wiley.com/doi/pdf/10.1002/andp.19163540702
1916 doi
-
[11]
M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Ar- zoumanian, K. C. Gendreau, S. Guillot, A. K. Harding, W. C. G. Ho, J. M. Lattimer, R. M. Ludlam, S. Mahmoodifar, S. M. Morsink, P. S. Ray, T. E. Strohmayer, K. S. Wood, T. Enoto, R. Foster, T. Okajima, G. Prigozhin,...
2019
-
[12]
M. T. Wolff, S. Guillot, S. Bogdanov, P. S. Ray, M. Kerr, Z. Ar- zoumanian, K. C. Gendreau, M. C. Miller, A. J. Dittmann, W. C. G. Ho, L. Guillemot, I. Cognard, G. Theureau, and K. S. Wood, The Astrophysical Journal Letters918, L26 (2021)
2021
-
[13]
J. W. T. Hessels, S. M. Ransom, I. H. Stairs, P. C. C. Freire, V . M. Kaspi, and F. Camilo, Science311, 1901 (2006), https://www.science.org/doi/pdf/10.1126/science.1123430
2006 doi
-
[14]
J. B. Hartle, ApJ150, 1005 (1967)
1967
-
[15]
J. B. Hartle and K. S. Thorne, ApJ153, 807 (1968)
1968
-
[16]
Kacskovics, D
B. Kacskovics, D. Barta, and M. Vasúth, As- tronomische Nachrichten344, e220109 (2023), https://onlinelibrary.wiley.com/doi/pdf/10.1002/asna.20220109
2023 doi
-
[17]
R. A. James, ApJ140, 552 (1964)
1964
-
[19]
Eriguchi and E
Y . Eriguchi and E. Mueller, A&A147, 161 (1985)
1985
-
[20]
Hachisu, ApJS61, 479 (1986)
I. Hachisu, ApJS61, 479 (1986)
1986
-
[21]
Komatsu, Y
H. Komatsu, Y . Eriguchi, and I. Hachisu, Monthly Notices of the Royal Astronomical Society237, 355 (1989), https://academic.oup.com/mnras/article- pdf/237/2/355/2983266/mnras237-0355.pdf
1989
-
[22]
G. B. Cook, S. L. Shapiro, and S. A. Teukolsky, The Astrophys- ical Journal398, 203 (1992)
1992
-
[23]
Stergioulas and J
N. Stergioulas and J. L. Friedman, The Astrophysical Journal 444, 306 (1995)
1995
-
[24]
X. Qu, S. Wang, and H. Tong, The Astrophysical Journal980, 3 (2025). 14
2025
-
[25]
Konstantinou, The Astrophysical Journal968, 83 (2024)
A. Konstantinou, The Astrophysical Journal968, 83 (2024)
2024
-
[26]
I. A. Rather, U. Rahaman, M. Imran, H. C. Das, A. A. Usmani, and S. K. Patra, Phys. Rev. C103, 055814 (2021)
2021
-
[27]
Goriely, N
S. Goriely, N. Chamel, and J. M. Pearson, Phys. Rev. C88, 024308 (2013)
2013
-
[28]
Reinhard and H
P.-G. Reinhard and H. Flocard, Nuclear Physics A584, 467 (1995)
1995
-
[29]
Bartel, P
J. Bartel, P. Quentin, M. Brack, C. Guet, and H.-B. Håkansson, Nuclear Physics A386, 79 (1982)
1982
-
[30]
Duta, private communication
M. Duta, private communication
-
[31]
Chabanat, P
E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, Nucl. Phys. A635, 231 (1998)
1998
-
[32]
J. M. Lattimer, Particles6, 30 (2023)
2023
-
[33]
Margueron, R
J. Margueron, R. Hoffmann Casali, and F. Gulminelli, Phys. Rev. C97, 025805 (2018)
2018
-
[34]
Shlomo, V
S. Shlomo, V . Kolomietz, and G. Colò, Eur. Phys. J. A20, 23 (2006)
2006
-
[35]
Piekarewicz, Journal of Physics G: Nuclear and Particle Physics37, 064038 (2010)
J. Piekarewicz, Journal of Physics G: Nuclear and Particle Physics37, 064038 (2010)
2010
-
[36]
Oertel, M
M. Oertel, M. Hempel, T. Klähn, and S. Typel, Rev. Mod. Phys. 89, 015007 (2017)
2017
-
[37]
Li and X
B.-A. Li and X. Han, Physics Letters B727, 276 (2013)
2013
-
[38]
Skyrme, Nuclear Physics9, 615 (1958)
T. Skyrme, Nuclear Physics9, 615 (1958)
1958
-
[39]
Vautherin and D
D. Vautherin and D. M. Brink, Phys. Rev. C5, 626 (1972)
1972
-
[40]
Chamel, S
N. Chamel, S. Goriely, and J. M. Pearson, Phys. Rev. C80, 065804 (2009)
2009
-
[42]
G. Baym, H. A. Bethe, and C. J. Pethick, Nucl. Phys. A175, 225 (1971)
1971
-
[43]
E. M. Butterworth and J. R. Ipser, ApJ204, 200 (1976)
1976
-
[44]
J. M. Bardeen, ApJ162, 71 (1970)
1970
-
[45]
G. B. Cook, S. L. Shapiro, and S. A. Teukolsky, ApJ422, 227 (1994)
1994
-
[46]
S. M. Morsink and L. Stella, ApJ513, 827 (1999), arXiv:astro- ph/9808227 [astro-ph]
1999
-
[47]
Astrophys
Nozawa, T., Stergioulas, N., Gourgoulhon, E., and Eriguchi, Y ., Astron. Astrophys. Suppl. Ser.132, 431 (1998)
1998
-
[48]
H. T. Cromartie, E. Fonseca, S. M. Ransom, P. B. Demor- est, Z. Arzoumanian, H. Blumer, P. R. Brook, M. E. DeCesar, T. Dolch, J. A. Ellis, R. D. Ferdman, E. C. Ferrara, N. Garver- Daniels, P. A. Gentile, M. L. Jones, M. T. Lam, D. R. Lorimer, R. S. Lynch, M. A. McLaughlin, C. ...
2020 arXiv
-
[49]
T. E. Riley, A. L. Watts, S. Bogdanov, P. S. Ray, R. M. Lud- lam, S. Guillot, Z. Arzoumanian, C. L. Baker, A. V . Bilous, D. Chakrabarty, K. C. Gendreau, A. K. Harding, W. C. G. Ho, J. M. Lattimer, S. M. Morsink, and T. E. Strohmayer, The As- trophysical Journal Letters887, L21 (2019)
2019
-
[50]
J. M. Lattimer and M. Prakash, Science304, 536 (2004), https://www.science.org/doi/pdf/10.1126/science.1090720
2004 doi
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.