Pith. sign in

REVIEW 3 major objections 6 minor 32 references

Hydrodynamic and Electromagnetic Discrepancies between Neutron Star and Black Hole Spacetimes

T0 review · 3 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Using a realistic neutron-star metric instead of Kerr changes near-surface accretion and fields by 10–70%.

desk verdict The paper's quantitative claims likely rest on an invalid 'Kerr' baseline; the qualitative caution is plausible but needs a validated comparison. read the letter →

arxiv 2505.05299 v1 pith:RNWCMYNN submitted 2025-05-08 gr-qc astro-ph.HE

Jonathan Gorard , James Juno , Ammar Hakim This is my paper · ORCID
classification gr-qcastro-ph.HE
keywords neutronstarmetricKerrmultipolemomentsno-hairtheoremgeneralrelativistichydrodynamicselectrodynamicsaccretionmagnetosphere
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper argues that replacing a black-hole spacetime with a physically realistic neutron-star spacetime changes simulated accretion and magnetospheres enough to matter. Using the Kerr metric, which encodes only mass and spin, the authors find fluid momentum density discrepancies of 30–40%, energy density discrepancies of 50–70%, and electric field strength discrepancies of 10–12% near the object's surface compared with a neutron-star metric carrying higher multipole moments. The message is that the no-hair simplification is not a safe approximation for hydrodynamic and electromagnetic processes close to a neutron star, even though it remains acceptable at large distances and for strongly magnetized accretion. If correct, this matters for modeling gravitational collapse, particle acceleration, and jet launching around neutron stars.

What carries the argument

The load-bearing object is an approximate analytic neutron-star spacetime (the Pappas metric) written in Weyl-Lewis-Papapetrou coordinates, with metric functions depending on mass $M$, angular momentum $J$, and the higher multipoles $M_2$, $S_3$, $M_4$. These higher moments are parameterized by $(\alpha,\beta,\gamma)$, with $\alpha$ the quadrupolar deformability, and the degenerate choice $\alpha=\beta=\gamma=1$ is taken to reproduce the Kerr metric. The multipole relations $\beta(\alpha)$ and $\gamma(\alpha)$ come from universal relations for realistic neutron-star equations of state. The comparison mechanism is simply the difference between running identical tetrad-based general relativistic hydrodynamics and electrodynamics solvers on this metric with $\alpha=5$ or $8$ versus $\alpha=1$, holding mass and spin fixed.

What would settle it

Compute the metric functions of the approximate neutron-star metric at $\alpha=\beta=\gamma=1$ for $J=0.4M^2$ and compare them, or the resulting geodesics and curvature invariants, with the exact Kerr solution across $2M\le r\le 5M$; if the differences are comparable to the 10–70% discrepancies reported in the fluid quantities, the comparison does not isolate the physics of neutron-star multipoles. Equivalently, rerun the same two test problems on exact Kerr in Kerr-Schild coordinates and on the $\alpha=\beta=\gamma=1$ case of the approximate metric; any large mismatch would invalidate the claimed metric equivalence.

Watch

Extended reading notes

Core claim

The central claim is that the exterior multipole structure of a rotating neutron star, not just its mass and spin, measurably changes the outcome of general relativistic hydrodynamic and electromagnetic simulations near the surface. The authors quantify this by running the same two test problems, Bondi-Hoyle-type wind accretion and a Wald-type magnetosphere, on two background spacetimes: the Kerr metric and an approximate neutron-star metric whose mass quadrupole, spin octupole, and mass hexadecapole are set to values appropriate for realistic equations of state at spin $j=0.4$. Across the neutron-star surface at $r\simeq 3M$, the momentum density changes by 30–40% on the counter-rotating side, the relativistic energy density by 50–70%, and the electric field strength by 10–12% in the axial direction. The paper concludes that using the Kerr metric for neutron-star simulations is a reasonable approximation for quantitative accretion rates in highly magnetized systems, but not for collapse or particle-acceleration physics near the surface.

Load-bearing premise

That the special case of the neutron-star metric used as the stand-in for the Kerr spacetime is a faithful stand-in, close enough that the differences in the simulations come from the neutron star's extra multipole moments rather than from the approximation.

Editorial extensions

If this is right

  • The Kerr metric remains a workable approximation for quantitative accretion-rate simulations around highly magnetized neutron stars, since the magnetic field dominates and plasma is forced to co-rotate.
  • Simulations of gravitational collapse, which are sensitive to accretion geometry, should use a neutron-star metric rather than Kerr to capture near-surface energy and momentum densities.
  • The 10–12% electric-field discrepancy in the axial direction is large enough to affect Blandford-Znajek-type jet launching and particle acceleration around neutron stars.
  • Differences seen in simple unmagnetized hydrodynamics and vacuum electromagnetic fields are expected to be amplified when more realistic plasma physics, such as instabilities, is included.
  • At large distances the Kerr approximation is safe, since realistic axisymmetric metrics converge to it in the weak-field limit.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the discrepancies scale with the quadrupolar deformability $\alpha$, the same comparison could serve as a diagnostic of the neutron-star equation of state: measuring near-surface flow or electric-field structure might constrain $\alpha$.
  • The same multipole sensitivity is likely to appear in neutron-star merger remnants and post-merger accretion disks, where matter is routinely modeled with black-hole-like metrics; testing that setting is a natural extension.
  • A sharper test would run the same simulations on a numerically constructed exact spacetime from a specific equation of state, rather than an approximate analytic metric, to separate physical multipole effects from approximation error.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. This paper uses the Gkeyll tetrad-based general-relativistic solvers to compare hydrodynamics and electrodynamics around a rotating neutron star and around a black hole. The neutron-star spacetime is represented by the approximate Pappas metric with quadrupolar deformability parameters α=5 and α=8 at spin j=0.4, while the black-hole baseline is generated by setting α=β=γ=1 inside the same Pappas metric, which the authors assert reproduces the Kerr metric. The reported results are large qualitative differences in fluid momentum density contours near r=3M, quantitative differences in momentum density of 15–40%, energy-density differences of 50–70%, and electric-field strength differences of 5–12%. The paper concludes that using the Kerr metric for neutron-star simulations may be inadequate for collapse dynamics and jet physics.

Significance. If the results are correct, they would provide a practically important caveat for the common practice of modeling neutron-star accretion and magnetospheres with Kerr or Schwarzschild backgrounds, and they would showcase a capable non-Kerr metric implementation. The manuscript is clearly written and the simulations are concrete. However, the central quantitative claims are currently not fully supported: the black-hole baseline is asserted to be Kerr without validation, the headline energy-density discrepancy is not displayed in any figure, and no convergence study is presented. These gaps are load-bearing because the specific percentages are the main output of the paper.

major comments (3)
  1. [Main text, paragraph following Eq. (3) ('Black hole spacetimes (for comparison purposes) are simulated by setting…] The equivalence of the α=β=γ=1 Pappas metric to the Kerr metric is asserted but not demonstrated. This is the foundation of the entire comparison, so it must be verified. In fact, the nonrotating limit is not Schwarzschild: with J=0 all higher multipole moments vanish, but the End Matter expression for f(ρ,z) still contains a 2M²/(ρ²+z²) term (and γ=0), so the resulting Weyl-Papapetrou metric is not Ricci-flat; (1/2)ln f does not satisfy the vacuum Laplace equation. Hence the α=β=γ=1 metric is only an approximate truncated expansion, not exact Kerr. The authors should either rerun the comparisons against the exact Kerr metric in Kerr-Schild coordinates, or quantitatively demonstrate that the α=β=γ=1 Pappas metric agrees with Kerr to an accuracy well below the 5–70% discrepancies they report, for example by computing the metric deviation or the Ricci scalar near r=3M for j=0.4.
  2. [Hydrodynamics section, text after Figure 2] The abstract and the main text claim that the peak relativistic energy density √γτ is modified by up to 50–70% near the neutron star surface, but no figure or table in the manuscript displays √γτ or its discrepancy. Only √γS_y profiles are shown (Figure 2). Since the 50–70% energy-density discrepancy is one of the headline quantitative results, the authors should show the corresponding energy-density profiles or provide the numerical values used to compute this range.
  3. [Numerical setup, hydrodynamics and electrodynamics simulations] All simulations are run at a single resolution of 1024×1024 cells with no convergence study. Because the reported discrepancies are quantitative percentages (30–40%, 50–70%, 10–12%), it is necessary to rule out numerical resolution as a significant contributor. Please provide at least one higher-resolution run (or a Richardson-extrapolation estimate) for one hydrodynamics case and one electrodynamics case, and state the resulting change in the quoted discrepancies.
minor comments (6)
  1. [Figure 1 caption and Figure 3 caption] Both captions describe the neutron-star metric as 'α=0.8', while the text and Figures 2 and 4 consistently use α=8. This is likely a typographical error, but it must be corrected because α=0.8 is not among the values simulated.
  2. [End Matter, metric functions] The notation is ambiguous: 'M2' is used for the mass quadrupole moment, but the expression for f(ρ,z) also contains a term that reads as '2M 2' (apparently 2M²), and 'M3' appears both as a multipole moment and as a power of M. Please use unambiguous symbols, e.g., M (mass), \mathcal{M}_2 (quadrupole), S_3 (spin octupole), and M_4 (hexadecapole), and distinguish powers of M with explicit exponents.
  3. [End Matter, metric functions and main text] The line element is written with '-f(dt - ω dφ)²' in the main text but with '-f(dt - Ω dφ)²' in the End Matter. Please standardize the symbol for the frame-dragging function.
  4. [Hydrodynamics section, description of Figure 2] The text states that the peak value of |√γS_y| is modified by 15–20% on the co-rotating side and 30–40% on the counter-rotating side, but it does not specify how these percentages are computed from the profiles (e.g., peak of the absolute value over the whole profile, or at a specific radius). Please define the estimator used.
  5. [Main text, 'For the purposes of our numerical comparisons, we shall assume that the neutron star surface lies at r=3M…] The choice of r=3M as the neutron-star surface is stated without justification or sensitivity analysis. Since all discrepancy estimates are quoted near this radius, please cite a source for the 3M–4M radius range and state whether the reported percentages are sensitive to the exact location of the surface within this range.
  6. [Main text, references] Reference [12] is an arXiv preprint (arXiv:2410.02549) on which the numerical methods rely; if a published version exists, it should be cited, otherwise the reliance on an unreviewed code description should be noted.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the two spacetimes being compared are externally specified metric inputs, and no parameter is fitted to the simulation output.

full rationale

The paper's central comparison is between the Pappas neutron-star metric (with α=5 and α=8, and β(α), γ(α) taken from the independent universal relations of Pappas & Apostolatos and Yagi et al.) and a 'black hole' baseline obtained by setting α=β=γ=1 in the same metric family. These are externally specified inputs, not quantities fit to the simulations' output. The reported discrepancy percentages are diagnostics computed after evolving the same initial data in two different fixed spacetimes; they are not fitted parameters renamed as predictions. The β(α) and γ(α) relations are imported from independent fits to numerical equation-of-state models, not derived from the present simulations. The self-citation to the Gkeyll solver paper [12] supports the numerical method but is not load-bearing for the physical conclusion: the hydrodynamics and electrodynamics results follow from direct simulation in the prescribed metrics. Whether α=β=γ=1 accurately reproduces the Kerr spacetime is a question of baseline validity and correctness, not circularity: if the identification is imperfect, the quantitative claims are weakened, but the derivation does not reduce to its own inputs. No load-bearing step is equivalent to its inputs by construction, and the paper is self-contained against external benchmarks for the metric and equation-of-state relations.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The central comparison depends on the Pappas metric and the assumption that setting its multipole deviations to Kerr values recovers Kerr. The key free parameter is the quadrupolar deformability alpha.

free parameters (5)
  • quadrupolar deformability alpha = 5 and 8
    Selected to represent moderate and low stiffness equations of state; governs the mass quadrupole M2 = -alpha j^2 M^3.
  • dimensionless spin j = 0.4
    Chosen as astrophysically plausible, below breakup for standard equations of state.
  • mass M = 0.3 (accretion), 0.5 (Wald)
    Set separately for the two test problems; not justified by physical scale.
  • spin multipole coefficient beta = beta = (-0.36 + 1.48 sqrt(alpha)^0.65)^3
    From empirical EOS fits by Pappas and Apostolatos (2014); not fitted here.
  • mass hexadecapole coefficient gamma = gamma = (-4.749 + 0.27613 sqrt(alpha)^1.5146 + 5.5168 sqrt(alpha)^0.22229)^4
    From empirical EOS fits by Yagi et al. (2014); not fitted here.
assumptions (5)
  • domain assumption The Pappas metric accurately describes the exterior spacetime of rotating neutron stars for alpha <= 8 and j <= 0.5.
    The paper relies on Pappas (2017) without validating against numerical spacetimes in the specific setup.
  • ad hoc to paper The Pappas metric with alpha=beta=gamma=1 is equivalent to the Kerr metric for the purposes of comparison.
    This is assumed for the baseline comparison; the paper provides no test that the parameter reduction reproduces Kerr.
  • domain assumption The empirical relations beta(alpha) and gamma(alpha) are valid for the chosen alpha values.
    These universal relations come from prior fits to numerical EOS models.
  • domain assumption The region 2M <= r <= 3M treated as unphysical does not contaminate the exterior flow.
    Fluids crossing the assumed surface at r=3M may experience unphysical metric behavior inside; the simulation excises at r=2M.
  • standard math The ideal gas equation of state with Gamma=5/3 and vacuum Maxwell equations are appropriate test problems.
    Standard choices for Bondi-Hoyle and Wald-type problems.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Hydrodynamic and Electromagnetic Discrepancies between Neutron Star and Black Hole Spacetimes." pith.science (2026). https://pith.science/paper/RNWCMYNN

@misc{pith2026250505299,
  author       = {Pith},
  title        = {Pith review of: Hydrodynamic and Electromagnetic Discrepancies between Neutron Star and Black Hole Spacetimes},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RNWCMYNN}},
  note         = {Machine review of arXiv:2505.05299}
}
read the original abstract

The exterior spacetime geometry surrounding an uncharged, spinning black hole in general relativity depends only upon its mass and spin. However, the exterior geometry surrounding any other rotating compact object, for example a neutron star, will generally depend upon higher moments in its multipole expansion, which will in turn be dependent upon the object's equation of state. Using general relativistic hydrodynamics and electrodynamics simulations, we illustrate that the presence or absence of these higher moments (assuming a physically realistic neutron star equation of state) has a significant qualitative effect near the surface of the compact object on the dynamics of unmagnetized accretion, and a smaller quantitative effect on the electromagnetic field configuration of its magnetosphere. In some places, the discrepancies in energy-momentum density are found to reach or exceed 50%, with electric field strength discrepancies in excess of 10%. We argue that many of these differences are likely to be amplified by the inclusion of more sophisticated plasma physics models, and are therefore likely to be relevant for the dynamics of gravitational collapse, and potentially also for particle acceleration and jet launching. These discrepancies suggest important limitations regarding the use of the Kerr metric when performing numerical simulations around neutron stars.

Figures

Figures reproduced from arXiv: 2505.05299 by the authors.

Figure 1
Figure 1. FIG. 1. The configuration of momentum density contours (as perceived by an Eulerian observer) at time [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Cross-sectional profiles of the [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. The configuration of electric flux surfaces (as perceived by an Eulerian observer) at time [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Cross-sectional profiles of the electric field strength [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 31 canonical work pages

  1. [1]

    G. D. Birkhoff and R. E. Langer,Relativity and Modern Physics:(Harvard University Press, 1927)

    work page 1927

  2. [2]

    J. B. Hartle, Slowly Rotating Relativistic Stars. I. Equa- tions of Structure, The Astrophysical Journal150, 1005 (1967)

    work page 1967

  3. [3]

    J. B. Hartle and K. S. Thorne, Slowly Rotating Relativis- tic Stars. II. Models for Neutron Stars and Supermassive Stars, The Astrophysical Journal153, 807 (1968)

    work page 1968

  4. [4]

    Heusler,Black hole uniqueness theorems, transferred to digital print ed., Cambridge lecture notes in physics No

    M. Heusler,Black hole uniqueness theorems, transferred to digital print ed., Cambridge lecture notes in physics No. 6 (Cambridge University Press, Cambridge, 2003)

    work page 2003

  5. [5]

    Chandrasekhar and J

    S. Chandrasekhar and J. C. Miller, On Slowly Rotating Homogeneous Masses in General Relativity, Monthly No- tices of the Royal Astronomical Society167, 63 (1974)

    work page 1974

  6. [6]

    Bradley and G

    M. Bradley and G. Fodor, Quadrupole moment of slowly rotating fluid balls, Physical Review D79, 044018 (2009)

    work page 2009

  7. [7]

    Paschalidis and N

    V. Paschalidis and N. Stergioulas, Rotating stars in rel- ativity, Living Reviews in Relativity20, 7 (2017)

    2017

  8. [8]

    C ¸ ıkınto˘ glu, K

    S. C ¸ ıkınto˘ glu, K. Y. Ek¸ si, and L. Rezzolla, GRMHD simulations of accreting neutron stars I: Non-rotating dipoles, Monthly Notices of the Royal Astronomical So- ciety517, 3212 (2022)

    work page 2022

Show all 32 references
  1. [9]

    Torres, F

    R. Torres, F. Cruz, T. Grismayer, R. Fonseca, and L. Silva, General relativistic particle-in-cell simulations of compact neutron star magnetospheres, inAPS Divi- sion of Plasma Physics Meeting Abstracts, APS Meeting Abstracts, Vol. 2021 (2021) p. JP11.008

    2021

  2. [10]

    Parfrey and A

    K. Parfrey and A. Tchekhovskoy, Accreting Neutron Stars in 3D General-relativistic Magnetohydrodynamic Simulations: Jets, Magnetic Polarity, and the Inter- change Slingshot, The Astrophysical Journal975, 57 (2024)

    2024

  3. [11]

    Horizon-adapted

    J. A. Font, J. M. Ib´ a˜ nez, and P. Papadopoulos, A “Horizon-adapted” Approach to the Study of Relativistic Accretion Flows onto Rotating Black Holes, The Astro- physical Journal507, L67 (1998)

    1998

  4. [12]

    Gorard, A

    J. Gorard, A. Hakim, J. Juno, and J. M. TenBarge, A tetrad-first approach to robust numerical algorithms in general relativity (2025), arXiv:2410.02549 [gr-qc]

    2025 arXiv

  5. [13]

    Pappas, An accurate metric for the spacetime around rotating neutron stars., Monthly Notices of the Royal As- tronomical Society , stx019 (2017)

    G. Pappas, An accurate metric for the spacetime around rotating neutron stars., Monthly Notices of the Royal As- tronomical Society , stx019 (2017)

    2017

  6. [14]

    Arnowitt, S

    R. Arnowitt, S. Deser, and C. W. Misner, Dynamical Structure and Definition of Energy in General Relativity, Physical Review116, 1322 (1959)

    1959

  7. [15]

    Pappas and T

    G. Pappas and T. A. Apostolatos, Effectively Universal Behavior of Rotating Neutron Stars in General Relativity Makes Them Even Simpler than Their Newtonian Coun- terparts, Physical Review Letters112, 121101 (2014)

    2014

  8. [16]

    K. Yagi, K. Kyutoku, G. Pappas, N. Yunes, and T. A. 6 Apostolatos, Effective no-hair relations for neutron stars and quark stars: Relativistic results, Physical Review D 89, 124013 (2014)

    2014

  9. [17]

    G. B. Cook, S. L. Shapiro, and S. A. Teukolsky, Rapidly rotating neutron stars in general relativity: Realistic equations of state, The Astrophysical Journal424, 823 (1994)

    1994

  10. [18]

    J. W. T. Hessels, S. M. Ransom, I. H. Stairs, P. C. C. Freire, V. M. Kaspi, and F. Camilo, A Radio Pulsar Spin- ning at 716 Hz, Science311, 1901 (2006)

    2006

  11. [19]

    Wang and A

    Y.-F. Wang and A. H. Nitz, Targeted search for gravita- tional waves from highly spinning light compact binaries, Monthly Notices of the Royal Astronomical Society528, 3891 (2024)

    2024

  12. [20]

    Haensel, Apparent radii of neutron stars and equation of state of dense matter, Astronomy & Astrophysics380, 186 (2001)

    P. Haensel, Apparent radii of neutron stars and equation of state of dense matter, Astronomy & Astrophysics380, 186 (2001)

    2001

  13. [21]

    J. A. Font, J. M. Ib´ a˜ nez, and P. Papadopoulos, Non- axisymmetric relativistic Bondi-Hoyle accretion on to a Kerr black hole, Monthly Notices of the Royal Astronom- ical Society305, 920 (1999)

    1999

  14. [22]

    R. M. Wald, Black hole in a uniform magnetic field, Phys- ical Review D10, 1680 (1974)

    1974

  15. [23]

    Philippov and M

    A. Philippov and M. Kramer, Pulsar magnetospheres and their radiation, Annual Review of Astronomy and Astro- physics60, 495 (2022)

    2022

  16. [24]

    Sekiguchi and M

    Y.-i. Sekiguchi and M. Shibata, Axisymmetric collapse simulations of rotating massive stellar cores in full gen- eral relativity: Numerical study for prompt black hole formation, Physical Review D71, 084013 (2005)

    2005

  17. [25]

    Parfrey and A

    K. Parfrey and A. Tchekhovskoy, Accreting neutron stars in 3d general-relativistic magnetohydrodynamic simula- tions: Jets, magnetic polarity, and the interchange sling- shot, The Astrophysical Journal975, 57 (2024)

    2024

  18. [26]

    A. M. Beloborodov, Electron–positron flows around mag- netars, The Astrophysical Journal777, 114 (2013)

    2013

  19. [27]

    R. D. Blandford and R. L. Znajek, Electromagnetic ex- traction of energy from Kerr black holes, Monthly Notices of the Royal Astronomical Society179, 433 (1977)

    1977

  20. [28]

    Parfrey, A

    K. Parfrey, A. Philippov, and B. Cerutti, First-Principles Plasma Simulations of Black-Hole Jet Launching, Phys- ical Review Letters122, 035101 (2019)

    2019

  21. [29]

    Parfrey and A

    K. Parfrey and A. Tchekhovskoy, General-relativistic simulations of four states of accretion onto millisecond pulsars, The Astrophysical Journal Letters851, L34 (2017)

    2017

  22. [30]

    Das and O

    P. Das and O. Porth, Three-dimensional GRMHD Simu- lations of Neutron Star Jets, The Astrophysical Journal Letters960, L12 (2024)

    2024

  23. [31]

    Banyuls, J

    F. Banyuls, J. A. Font, J. M. Ibanez, J. M. Marti, and J. A. Miralles, Numerical{3 + 1}General Relativistic Hydrodynamics: A Local Characteristic Approach, The Astrophysical Journal476, 221 (1997)

    1997

  24. [32]

    Gottlieb, C.-W

    S. Gottlieb, C.-W. Shu, and E. Tadmor, Strong Stability- Preserving High-Order Time Discretization Methods, SIAM Review43, 89 (2001). END MA TTER Metric F unctions for the Pappas Metric The particular forms of the metric functionsf(ρ,z),ω(ρ,z), andγ(ρ,z) appearing within the g...

    2001

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.