{"id":"31136851-d27e-4ec4-9d3e-8c75ea439129","arxiv_id":"2505.05299","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Using the Pappas neutron star metric in general relativistic hydrodynamics and electrodynamics simulations shows that near-surface energy and momentum densities differ from Kerr by tens of percent, while electric fields differ by about 10 percent.","lead":"This paper compares simulations of accretion and magnetospheres around neutron stars using a realistic neutron star metric versus the Kerr black hole metric. It finds large differences in fluid momentum near the surface, and smaller differences in electric field strength, which could affect simulations of gravitational collapse and jet formation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'black hole' baseline is not validated: setting α=β=γ=1 in the Pappas metric is asserted to give Kerr, but the metric is an approximate expansion, and the J=0 limit is not Schwarzschild, so the reported 30-70% discrepancies may partly be approximation error.","rationale":"The central claim depends on the comparison baseline being the actual Kerr spacetime. The paper's one-line assertion that α=β=γ=1 'corresponds to the Kerr metric' is load-bearing: if false, the 30-70% and 10-12% numbers measure differences between two members of the same approximate metric family, not between the Kerr metric and a physically realistic neutron-star metric. The End Matter formulas show the Pappas metric is a truncated multipolar expansion, and the J=0 limit is not Schwarzschild (it is not Ricci-flat), so the equivalence is far from established. This is the same weakest assumption the reader identified, and we agree with that assessment. The qualitative conclusion that neutron-star higher multipole moments matter near the surface is plausible and likely survives, but the quantitative claims in the abstract are not supported until the α=β=γ=1 baseline is validated against exact Kerr. Therefore the reader's CONDITIONAL verdict remains appropriate, pending the concrete test.","tokens_in":9060,"tokens_out":26439,"duration_ms":241178,"concrete_test":"Evaluate the Einstein tensor G_μν of the Pappas metric with α=β=γ=1, J=0.4M^2 at representative points on the r=3M sphere (e.g., pole and equator). If any component of G_μν normalized by M^{-2} is nonzero at machine precision, the metric is not Ricci-flat and therefore is not the exact Kerr solution, confirming that the reported discrepancies are at least partly an artifact of the approximate Pappas baseline. Additionally, compare the α=β=γ=1 metric components to the exact Kerr metric in the same Weyl-Lewis-Papapetrou coordinates at those points; a maximum relative deviation above a few percent would require the quantitative claims to be re-derived with an exact Kerr baseline.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative comparison is between a Pappas neutron-star metric with α=5 or 8 and a 'black hole' metric obtained by setting α=β=γ=1 in the same Pappas metric. The paper asserts that this degenerate case 'corresponds to the Kerr metric,' but no demonstration is provided, and the assertion is doubtful. The Pappas metric is an analytic approximation to numerical neutron-star spacetimes, and setting finitely many multipole moments to Kerr values does not uniquely determine a spacetime; the End Matter formulas are truncated multipolar expansions, not the exact Kerr solution. A concrete red flag is the nonrotating limit: with J=0, all higher moments vanish and the End Matter functions reduce to f=1-2M/r, ω=0, γ=0, giving ds^2=-(1-2M/r)dt^2+(1-2M/r)^{-1}(dρ^2+dz^2)+ρ^2 dφ^2. For a Weyl-Papapetrou metric, vacuum field equations require the function ψ=(1/2)ln f to satisfy Laplace's equation in cylindrical coordinates; here ∇²ψ = -2M^2/[r^2(r-2M)^2] ≠ 0, so this limit is not Ricci-flat and not Schwarzschild. Hence the α=β=γ=1 case is not the exact Kerr spacetime. The reported 30-70% fluid discrepancies and 10-12% electric-field discrepancies are therefore computed against an approximate Pappas baseline, not against the Kerr metric used in previous simulations, and the quantitative headline claims are unsubstantiated until the baseline is verified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":9460,"tokens_out":9947,"duration_ms":96509,"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":[{"comment":"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.","section":"Main text, paragraph following Eq. (3) ('Black hole spacetimes (for comparison purposes) are simulated by setting…"},{"comment":"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.","section":"Hydrodynamics section, text after Figure 2"},{"comment":"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.","section":"Numerical setup, hydrodynamics and electrodynamics simulations"}],"minor_comments":[{"comment":"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.","section":"Figure 1 caption and Figure 3 caption"},{"comment":"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.","section":"End Matter, metric functions"},{"comment":"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.","section":"End Matter, metric functions and main text"},{"comment":"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.","section":"Hydrodynamics section, description of Figure 2"},{"comment":"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.","section":"Main text, 'For the purposes of our numerical comparisons, we shall assume that the neutron star surface lies at r=3M…"},{"comment":"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.","section":"Main text, references"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of the journal and addresses a question of practical importance for numerical relativity and astrophysics. The referee's main concern is not the qualitative idea but the unvalidated 'Kerr' baseline, which is load-bearing for every quantitative claim. This is fixable by redoing the simulations with an exact Kerr metric or by providing a rigorous error estimate for the α=β=γ=1 Pappas metric, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper's central numbers are likely not sound because the 'Kerr' baseline is not actually Kerr. That said, the question it asks is worth asking, and the simulations are a concrete step.\n\nGorard et al. compare GRHD and GRED simulations around a Pappas neutron-star metric against the same metric with alpha=beta=gamma=1, which they claim reduces to Kerr. The qualitative message—that the no-hair metric can misrepresent near-surface accretion—is plausible and probably right, but the specific discrepancies (30–40% momentum, 50–70% energy, 10–12% electric field) are computed against an unvalidated stand-in. The End Matter formulas do not appear to reduce to Schwarzschild even for J=0: the metric function f retains terms like 2M^2/R^2 and higher, and a quick check of the Weyl-Papapetrou vacuum equation shows ∇²ψ≠0. So alpha=beta=gamma=1 is not exact Kerr. Without a check that this limit actually is Kerr (or that the deviations are small compared to the claimed effects), the numbers are not substantiated.\n\nWhat is genuinely useful: the paper picks two standard test problems (Bondi-Hoyle accretion and Wald magnetosphere), runs them in the Gkeyll tetrad-based solver, and makes a side-by-side comparison. That is the kind of concrete demonstration that could help the numerical relativity community decide when the Kerr approximation is safe. The authors also correctly note that magnetic fields suppress some of the fluid effects, so the practical impact for magnetized accretion is modest—that nuance is honest.\n\nOther soft spots: no convergence study, no resolution test, and the largest claimed effect (50–70% energy density) is stated but not shown in any figure. No code or data is released, which is not fatal for a Letter but would have made it easier to verify the baseline issue.\n\nMy sense: the qualitative conclusion is likely correct, but the paper overstates the quantitative case. A referee should insist on either using the actual Kerr metric in the same solver or demonstrating that the Pappas limit is Kerr to within negligible error. If that validation checks out, this becomes a useful practical caution. As it stands, I wouldn't cite the numbers.\n\nI'd still send this to peer review—the question is relevant, the simulations are serious, and the baseline issue is fixable. But it needs revision before publication.","headline":"The paper's quantitative claims likely rest on an invalid 'Kerr' baseline; the qualitative caution is plausible but needs a validated comparison.","tokens_in":9946,"tokens_out":3118,"would_cite":false,"duration_ms":29438,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Using a realistic neutron-star metric instead of Kerr changes near-surface accretion and fields by 10–70%.","keywords":["neutron star metric","Kerr metric","multipole moments","no-hair theorem","general relativistic hydrodynamics","general relativistic electrodynamics","accretion","magnetosphere"],"falsifier":"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.","tokens_in":8864,"feed_emoji":"🌀","tokens_out":7449,"duration_ms":67917,"temperature":0.7,"pith_summary":"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.","feed_headline":"Neutron star spacetimes alter accretion by up to 70 percent","feed_subtitle":"Ignoring higher multipole moments shifts simulated energy, momentum, and electric fields near the surface; collapse and jet models should…","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the approximate neutron-star metric used in all simulations.","marker":"[13]"},{"why":"Gives the parametrization of the higher multipole moments $M_2, S_3, M_4$ in terms of $(\\alpha,\\beta,\\gamma)$ and identifies the Kerr limit.","marker":"[15]"},{"why":"Provides the universal relations linking $\\beta$ and $\\gamma$ to the quadrupolar deformability $\\alpha$ for realistic equations of state.","marker":"[16]"},{"why":"Supplies the tetrad-based general relativistic hydrodynamic and electromagnetic solvers used to run both test problems.","marker":"[12]"},{"why":"Defines the non-axisymmetric Bondi-Hoyle wind accretion setup used for the hydrodynamics comparison.","marker":"[21]"},{"why":"Defines the uniform-field (Wald) magnetosphere setup used for the electromagnetic comparison.","marker":"[22]"},{"why":"Motivates the use of horizon-penetrating coordinates and explains why Kerr is the usual simulation choice.","marker":"[11]"},{"why":"Supports the conclusion that strongly magnetized neutron-star accretion is channeled along field lines, limiting the practical impact of the metric choice there.","marker":"[23]"}],"fun_headline_variants":["Kerr metric misses up to 70% of neutron star near-surface physics","Neutron star shape matters: Kerr metric error up to 70%","Beyond Kerr: multipole moments alter neutron star sims by 70%","Kerr metric fails near neutron stars, altering simulations by 70%","Neutron star multipole moments reshape near-surface physics by 70%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Kerr metric misses up to 70% of neutron star near-surface physics","Neutron star shape matters: Kerr metric error up to 70%","Beyond Kerr: multipole moments alter neutron star sims by 70%","Kerr metric fails near neutron stars, altering simulations by 70%","Neutron star multipole moments reshape near-surface physics by 70%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001362,"raw_usage":{"total_tokens":5540,"prompt_tokens":972,"completion_tokens":4568,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":588,"completion_tokens_details":{"reasoning_tokens":4468}},"tokens_in":588,"tokens_out":4568,"duration_ms":36213,"temperature":1.0,"reasoning_tokens":4468,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:07:25.265395+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Pappas, An accurate metric for the spacetime around rotating neutron stars., Monthly Notices of the Royal As- tronomical Society , stx019 (2017)","cited_arxiv_id":null,"evidence_quote":"Supplies the approximate neutron-star metric used in all simulations."},{"cited_title":"Pappas and T","cited_arxiv_id":null,"evidence_quote":"Gives the parametrization of the higher multipole moments $M_2, S_3, M_4$ in terms of $(\\alpha,\\beta,\\gamma)$ and identifies the Kerr limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the universal relations linking $\\beta$ and $\\gamma$ to the quadrupolar deformability $\\alpha$ for realistic equations of state."},{"cited_title":"A Tetrad-First Approach to Robust Numerical Algorithms in General Relativity","cited_arxiv_id":"2410.02549","evidence_quote":"Supplies the tetrad-based general relativistic hydrodynamic and electromagnetic solvers used to run both test problems."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the non-axisymmetric Bondi-Hoyle wind accretion setup used for the hydrodynamics comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the uniform-field (Wald) magnetosphere setup used for the electromagnetic comparison."},{"cited_title":"Horizon-adapted","cited_arxiv_id":null,"evidence_quote":"Motivates the use of horizon-penetrating coordinates and explains why Kerr is the usual simulation choice."},{"cited_title":"Philippov and M","cited_arxiv_id":null,"evidence_quote":"Supports the conclusion that strongly magnetized neutron-star accretion is channeled along field lines, limiting the practical impact of the metric choice there."}],"review_version":1}