REVIEW 3 major objections 5 minor 1 cited by
Stability of the spacetime of a magnetized compact object
T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Magnetized compact objects are stable, and stronger magnetization speeds up the relaxation of scalar perturbations.
desk verdict A credible first QNM study of a magnetized compact object, but the magnetization-driven damping claim needs convergence evidence and a more honest scope statement. 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 central object is the Gutsunaev-Manko metric in prolate spheroidal coordinates $(x,y,\phi)$: an exact axisymmetric electrovacuum solution describing a body with a dipole magnetic moment, regularized by a perfect-reflection (mirror) boundary at $x=x_{\mathrm{min}}$ that hides the spacetime singularities. Scalar perturbations are decomposed into spherical harmonics, reducing the Klein-Gordon equation to a coupled system of one-dimensional wave equations with coupling matrices $A_{\ell'}^{\ell}$ and $B_{\ell'}^{\ell}$; these are evolved in the time domain with a fourth-order-in-space, second-order-in-time finite-difference scheme. In the $\alpha=0$ limit the system decouples and the quasinormal-mode problem reduces to a transcendental equation involving confluent Heun functions, which provides the analytic benchmark for the numerics.
What would settle it
Recompute the fundamental mode at small magnetization, for example $\alpha=0.005$ at $\tilde{x}_{\mathrm{min}}=-15$ where the reported shift in $\mathrm{Im}\,\omega$ is about $10^{-6}$, using $\ell_{\mathrm{max}}=7$ or $9$ and a finer grid; if $\mathrm{Im}\,\omega$ no longer increases monotonically with $\alpha$, or if the tiny shift vanishes at higher resolution, the claimed magnetization-driven stabilization is a numerical artifact.
Extended reading notes
Core claim
The paper argues that a spheroidal, non-rotating magnetized body built by imposing a Dirichlet reflection condition on a central region of the Gutsunaev-Manko spacetime has a stable exterior for scalar perturbations: every tested combination of the magnetization parameter $\alpha$, the mirror position $\tilde{x}_{\mathrm{min}}$, and the multipole number $\ell$ shows decaying quasinormal oscillations. The central quantitative finding is that the fundamental mode shifts as $\omega(\alpha)/\omega(0) \approx 1 + b\alpha^2$ with $b>0$, and the fitted coefficient for $\mathrm{Im}\,\omega$ is systematically larger than the coefficient $4$ expected from the $\alpha$-dependence of the mass alone. Because of this excess, the paper concludes that the faster decay is not merely a mass effect but a genuine magnetization-induced stabilization of the spacetime. The paper also reports that the power-law tail exponent $r$ increases linearly with the size of the star and decreases quadratically with $\alpha$, and it validates the numerical scheme against an analytic confluent-Heun solution in the $\alpha=0$ Schwarzschild limit.
Load-bearing premise
The extracted quasinormal frequencies are numerically converged: the multipole cutoff $\ell_{\mathrm{max}}=5$ and the chosen grid resolve the magnetization-induced shifts, some of which are smaller than the last quoted decimal.
Editorial extensions
If this is right
- For this model, mirror-reflected scalar perturbations never grow: every scanned value of $\alpha$, $\tilde{x}_{\mathrm{min}}$, and $\ell$ shows decaying ringdown, so the exterior region is stable throughout the allowed parameter space.
- Stronger magnetization shortens the ringdown: because $\mathrm{Im}\,\omega$ rises faster with $\alpha$ than the mass rescaling predicts, magnetized stars should relax to equilibrium more quickly than their mass alone would indicate.
- The late-time signal remains Schwarzschild-like: the tail exponent clusters near $r\approx 6$, increasing linearly with the star's size and decreasing quadratically with magnetization.
- Because the model reproduces magnetar-scale field strengths for small $\alpha$, the result gives a qualitative expectation that highly magnetized neutron-star exteriors are stable and faster-relaxing in their scalar channel.
- The confluent-Heun solution in the Schwarzschild limit matches the time-domain numerics, providing a cross-check of the numerical method and a starting point for analytic approximations at small magnetization.
Reading between the lines
- If the stabilizing effect carries over to gravitational and electromagnetic perturbations, which the paper does not treat, then highly magnetized neutron stars would radiate their ringdown more quickly than unmagnetized ones at the same mass, a difference that could show up in the damping of post-merger or flare oscillations.
- The mirror boundary is a crude stand-in for a stellar surface; replacing it with an absorbing or radiating boundary, or with a fluid interior, could either erase or amplify the $\alpha$-effect, so repeating the analysis with a more physical surface is a direct testable extension.
- The tail exponent's slight decrease with $\alpha$ means the ringdown is faster while the very late-time signal decays a little more slowly; for long-lived magnetar afterglows, this competition between the two regimes could be observationally relevant.
- A charged scalar field would couple directly to the dipole vector potential $A_\phi$; the paper mentions this as future work, and it would provide a sharper test of whether the stabilization is tied to the geometry itself rather than to the neutrality of the perturbing field.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies massless scalar perturbations on the exterior of a magnetized compact object described by the Gutsunaev-Manko spacetime, with a perfectly reflecting (mirror) boundary at x=xmin. It derives a coupled system of wave equations from a multipole decomposition, integrates them in the time domain with a fourth-order spatial/second-order temporal finite-difference scheme, and cross-checks the alpha=0 limit against an analytic confluent-Heun calculation. The authors report that the exterior is stable in all cases considered, that the imaginary part of the fundamental quasinormal frequency grows with the magnetization parameter alpha beyond the simple mass-rescaling M(alpha), and that the late-time power-law tail depends on both the mirror position and alpha. The qualitative stability result is plausible, but the quantitative claim of magnetization-induced stabilization is currently not established because the reported alpha-dependent shifts are close to or smaller than the apparent numerical precision and no convergence or error analysis is provided.
Significance. If the central effect is real, the paper would provide a useful model of how strong magnetic fields affect QNM stability of compact objects, with a concrete analytic benchmark in the Schwarzschild limit and a reproducible public code. The alpha=0 comparison between time-domain data and Heun-function roots gives genuine support to the numerical scheme, and the public availability of the code is a strength. However, the significance of the main new finding depends on showing that the alpha-induced shifts in Im omega are not numerical artifacts and that the explored parameter region supports the 'entire parameter space' claim; neither point is demonstrated in the current manuscript.
major comments (3)
- [Section III.B, Table I, Eqs. (30)-(31), Figure 4] The central claim that larger magnetization alpha increases Im omega beyond the M(alpha) rescaling rests on differences that are at or below the reported precision, with no convergence or uncertainty quantification. For xmin=-15, Table I lists Im omega = 0.000749 for every alpha from 0 to 0.01, while the Figure 4 fit with b_im approximately 14.8 would imply a change in the sixth decimal at alpha=0.01; for xmin=-10 the sequence is non-monotonic, with Im omega = 0.008041 at alpha=0.0075 and 0.008039 at alpha=0.01. The finite-difference scheme is described, but the grid spacing, time step, outer-boundary placement and treatment, and lmax truncation error are not reported, and no Richardson or lmax-convergence test is shown. In particular, the statement in Section III.B that the QNM spectrum is effectively independent of the cutoff is not demonstrated by any numerical data for the reported runs, all of which use lmax=5. Please provide convergence data and error bars for the quoted frequencies; without them, the magnetization-induced shift cannot be distinguished from discretization or truncation error.
- [Abstract, Section IV, Figure 6] The abstract and final remarks claim stability 'in the entire parameter space', but the numerical scans cover only three mirror positions and alpha up to 0.01, 0.03, and 0.11 respectively for the QNM analysis (Table I). The tail analysis in Figure 6 extends alpha up to 2, which is far beyond the physical bound alpha^2 < 1/3 from Eq. (9). No results are presented for other values of xmin or alpha, and no discussion is given of how the claim extends to the full allowed range of parameters or to regions near the spacetime singularities. Please either restrict the conclusion to the explored region or extend the scans and explicitly address the allowed parameter space.
- [Section III.A, Figures 1 and 2, Table I] The safety of the computational domain with respect to the singularities of the coupling matrices is not established for the parameter values actually used. Figures 1 and 2 show divergences at x~sing approximately -9.07 for alpha=0.057, but the text does not report where these divergences occur for the values alpha=0.01, 0.03, and 0.11 used in Table I. Since the mirror at xmin=-10 is used with alpha up to 0.03, and the singular location may depend on alpha, it is not clear that the singularities remain outside the integration domain for all reported runs. This is a load-bearing issue because a divergence inside the domain would invalidate the fourth-order finite-difference stencil. Please provide the singular-location curve as a function of alpha and verify that xmin is outside the singular region for every reported case.
minor comments (5)
- [Section III.B] There is a typo in the text: 'a par of integers' should be 'a pair of integers'.
- [Figure 4] The caption refers to '(Top)' and '(Bottom)' panels, but the figure itself does not appear to label the panels explicitly; please add (a) and (b) labels or otherwise make the panel correspondence unambiguous.
- [Table I] The number of quoted decimals varies between columns and rows, and the xmin=-15 column shows no change in Im omega for all alpha values. Please use uniform significant figures and state the numerical precision so that unchanged entries are not misinterpreted as a null result.
- [Section III.B, Figure 6] The tail-exponent fit r0(1+a xmin)(1+b alpha^2) is presented without residuals or error bars, and the text notes that the data points are clustered up to the graphics resolution. Please quantify the fit uncertainty, especially for the alpha dependence, since this is a secondary but still quantitative claim.
- [Section II] The statement that x=1 and -1<=y<=1 represents 'the event horizon' in the alpha=0 limit should be phrased more carefully when the object is later interpreted as a reflecting star, since the same surface becomes the mirror boundary x=xmin for the compact-object model.
Circularity Check
No significant circularity: the stability claim is derived from direct time-domain integration with an independent analytic α=0 benchmark; the sole same-group citation is a numerical-method reference whose discretization is reproduced in the paper and whose code is public.
full rationale
The paper's central claim — that the imaginary part of the fundamental quasinormal frequency grows with magnetization beyond the mass rescaling M(α) — is obtained by numerically integrating the coupled wave equations (22) with finite differences (30)-(31), starting from the Gutsunaev-Manko metric (5), the Klein-Gordon equation (16), and a Dirichlet mirror boundary (39). No parameter is fitted to the QNM output and then renamed a prediction: the b coefficients in Fig. 4 are post-hoc quadratic fits summarizing the same extracted frequencies, not inputs to the extraction. The α=0 sector is independently benchmarked: Eq. (48) is a transcendental equation in confluent Heun functions, solved with Newton's method, and Table II agrees with Table I, validating the numerical scheme. The only same-group citation is Ref. [24] for the finite-difference discretization, but the paper explicitly writes the discretized derivatives (30)-(31) and provides the public code at [25], so the method is self-contained and the citation is not load-bearing. The magnetization-dependent stabilization is therefore an honest numerical result, not an identity or a self-citation chain. Any concerns about truncation error or the extent of the parameter scan are numerical/correctness issues, not circularity.
Assumptions & free parameters
free parameters (4)
- alpha =
scanned up to 0.01 (xmin=-15), 0.03 (xmin=-10), 0.11 (xmin=-5); up to 2 in tail analysis
- xmin =
-15, -10, -5 for QNM runs; [-3,1] for tail runs
- b_re, b_im =
b_re approx 4.9, 5.4, 5.4; b_im approx 9.7, 13.7, 14.8 for xmin=-5, -10, -15
- r0, a, b_tail =
r0 approx 5.97, a approx 7e-4, b_tail approx -3.3e-2
assumptions (5)
- standard math The Gutsunaev-Manko metric (Eqs. 5-7) is an exact electrovacuum solution of Einstein's equations with a dipole magnetic moment.
- domain assumption The curvature singularities of the GM spacetime are hidden by the reflecting boundary at x=xmin and do not affect the exterior dynamics.
- domain assumption Massless scalar perturbations obey Eq. (16), and the multipole decomposition (18) with cutoff lmax=5 is sufficient.
- domain assumption The finite-difference scheme (fourth-order space, second-order time) is stable and converged with the reported step sizes.
- standard math Physically admissible stars require 3 alpha^2 < 1 so that the mass in Eq. (9) is positive.
Cite this review
Pith. "Pith review of Stability of the spacetime of a magnetized compact object." pith.science (2026). https://pith.science/paper/HNKJ4QQ7
@misc{pith2026241111117,
author = {Pith},
title = {Pith review of: Stability of the spacetime of a magnetized compact object},
year = {2026},
howpublished = {\url{https://pith.science/paper/HNKJ4QQ7}},
note = {Machine review of arXiv:2411.11117}
}
read the original abstract
We investigate the stability of scalar perturbations around a magnetized stationary compact object in General Relativity. The considered object is one of the simplest exact solutions of Einstein electrovacuum equations corresponding to a spheroidal body endowed with a dipole magnetic moment. It is effectively constructed by imposing a perfect reflection (mirror) boundary condition on a central region of the Gutsunaev-Manko spacetime. A time-domain analysis of the perturbations reveals a quasinormal phase followed by a power-law decaying tail. Our findings suggest that the exterior region of the magnetized compact object is stable in the entire parameter space. Moreover, the system tends to become generically more stable the stronger the magnetization of the central object is. Such findings can be useful for the qualitative understanding of more realistic astrophysical situations involving highly magnetized sources.
Figures
Figures from the paper (2 more)
Forward citations
Cited by 1 Pith paper
-
Universality in quasinormal modes of a magnetized black hole
For charged scalar perturbations of an Ernst-Schwarzschild black hole, the quasinormal-mode frequency scales as |q - q_c|^{1/2} near a critical charge q_c, with a mode-independent exponent.
Reference graph
Works this paper leans on
- [1]
-
[2]
The Event Horizon Telescope Collaboration, ApJL 930 L17 (2022). [arXiv:2311.09484]
arXiv 2022
-
[3]
H. Stephani, D. Kramer, M. MacCallum, C.Hoenselaers, and E. Herlt, Exact solutions of Einstein ’s field equations, Cambridge University Press (2009)
work page 2009
- [4]
- [5]
- [6]
- [7]
-
[8]
Stix, The Sun: an introduction , 2nd Ed., Springer (2002)
M. Stix, The Sun: an introduction , 2nd Ed., Springer (2002)
work page 2002
Show all 31 references
-
[9]
Kaspi and A.M
V.M. Kaspi and A.M. Beloborodov. Annu. Rev. Astron. Astrophys. 55, 261 (2017). [arXiv:1703.00068]
2017 arXiv
-
[10]
de Freitas and A
V.P. de Freitas and A. Saa, Phys. Rev. D 95, 124040 (2017). [arXiv:1703.10883]
2017 arXiv
-
[11]
Polanco, J
J.D. Polanco, J. Ayala, and M. Ujevic, Class. Quantum Grav. 41, 065017 (2024). [arXiv:2310.15979]
2024 arXiv
-
[12]
Kirsten, M
F. Kirsten, M. P. Snelders, M. Jenkins, and et al, Nature Astronomy 5, 414 (2020). [arXiv:2007.05101]
2020
- [13]
- [14]
-
[15]
Abdalla, E.G.M
E. Abdalla, E.G.M. Ferreira, R.G. Landim, and et al, Astron. Astrophys. 664, A14 (2022). [arXiv:2107.01633]
2022 arXiv
-
[16]
dos Santos, R.G
M.V. dos Santos, R.G. Landim, G.A. Hoerning, and et al, Astron. Astrophys. 681, A120 (2023). [arXiv:2308.06805]
2023 arXiv
-
[17]
Abdalla, B
E. Abdalla, B. Cuadros-Melgar, R.D.B. Fontana, J. Oliveira, E. Papantonopoulos, and A.B. Pavan, Phys. Rev. D 99, 104065 (2019). [arXiv:1903.10850]
2019 arXiv
-
[18]
Zhu, S.J
Z. Zhu, S.J. Zhang, C.E. Pellicer, B. Wang, and E. Abdalla, Phys. Rev. D, 90, 044042 (2014). [arXiv:1405.4931] [Addendum: Phys. Rev. D 90, 049904 (2014)]
2014 arXiv
-
[19]
Mora-Chaverri, E
A. Mora-Chaverri, E. Santiago-Leandro, and F. Frutos- Alfaro, Classical general relativity effects by magne- tars with slow rotation and a magnetic dipole (2024). [arXiv:2405.03503]
2024 arXiv
-
[20]
Kokkotas and B.G
K.D. Kokkotas and B.G. Schmidt, Living Rev. Relat. 2, 1 (1999). [arXiv:gr-qc/9909058]
1999 arXiv
-
[21]
Konoplya and A
R.A. Konoplya and A. Zhidenko, Rev. Mod. Phys. 83, 793 (2011). [arXiv:1102.4014]
2011 arXiv
- [22]
-
[23]
Manko, Class
V.S. Manko, Class. Quantum Grav. 10, L239 (1993)
1993
-
[24]
Molina, D
C. Molina, D. Giugno, E. Abdalla, A. Saa, Phys. Rev. D 69, 104013 (2004). [arXiv:gr-qc/0309079]
2004 arXiv
-
[25]
https://github.com/bingotcollab/Finite-difference-GM
-
[26]
Thus, the gen- eral solution UX (z) around the point X = 0 or X = ∞ can be written as UX (z) = C + X H + X (z) + C − X H − X (z) , (38) where C ± X are constants
such that the linearly independent solutions take the forms H + ∞ (z) ∼ e−β(z+ln z) z X ν≥0 a+ ν zν , a + 0 = 1, (37a) H − ∞ (z) ∼ 1 z X ν≥0 a− ν zν , a − 0 = 1, (37b) with the coefficients a± ν determined by the known recur- rence relations of the Heun functions [26]. Thus, t...
2014
-
[27]
Slavyanov and W
S.Y. Slavyanov and W. Lay, Special functions: a uni- fied theory based on singularities, Oxford University Press (2000)
2000
-
[28]
Horta¸ csu, Adv.High Energy Phys
M. Horta¸ csu, Adv.High Energy Phys. 2018, 8621573 (2018). [arXiv:1101.0471]
2018 arXiv
-
[29]
Cavalcante, B
J.P. Cavalcante, B. Carneiro da Cunha, Phys. Rev. D 104, 124040 (2021). [arXiv:2109.06929]
2021 arXiv
-
[30]
Amado, B
J.B. Amado, B. Carneiro da Cunha, and E. Pallante, Phys. Rev. D 105, 044028 (2022). [arXiv:2110.08349]
2022 arXiv
-
[31]
London, D
L. London, D. Shoemaker, and J. Healy, Phys. Rev. D, 90, 124032 (2014). [arXiv:1404.3197]
2014 arXiv
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.