Magnetized neutron stars: perturbative versus fully-numerical approaches
Pith reviewed 2026-05-20 04:17 UTC · model grok-4.3
The pith
Perturbative methods accurately model neutron star deformation by magnetic fields at all observed strengths but break down above a few times 10^16 G.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We compare the perturbative approach of Konno et al. (1999), which assumes small deformations from dipole currents, with a fully numerical solution based on the LORENE library. Both compute the magnetic field distribution and stellar deformation for varying polar field values, compactness, and equations of state, restricted to purely poloidal fields. The perturbative method breaks down for very high polar magnetic field values typically above a few times 10^16 G, yet yields very good results for observed values even in magnetars. In contrast, the numerical code exhibits resolution problems for relatively low magnetic field values typically around 10^10 G, leading to imprecise computation of
What carries the argument
Direct comparison of the Konno-99 perturbative expansion, which assumes small deformations induced by dipole currents, against the full numerical solution of the Einstein-Maxwell system using the LORENE library, applied to poloidal magnetic fields.
If this is right
- For polar fields up to 10^15 G typical of magnetars, the perturbative method supplies reliable values for stellar deformation and quadrupole moment.
- Above a few times 10^16 G the full numerical solution must be used because the perturbative expansion ceases to be valid.
- At field strengths near 10^10 G the numerical code requires substantially higher resolution to produce accurate deformation results.
- The preferred computational method for a given neutron star model depends on the expected polar field strength.
Where Pith is reading between the lines
- If real neutron stars contain substantial toroidal field components, the validity ranges reported for each method would shift and require separate verification.
- These findings can guide which technique to apply when estimating continuous gravitational wave signals from deformed magnetized stars.
- Extending the same comparison to rotating configurations or mixed poloidal-toroidal fields would test whether the current limits remain unchanged.
Load-bearing premise
The results apply only when the magnetic field inside the star is purely poloidal.
What would settle it
A side-by-side calculation of the mass quadrupole moment from both methods at a polar field of 10^15 G that shows close agreement, combined with clear disagreement at 5 times 10^16 G, would confirm the stated ranges of validity.
Figures
read the original abstract
(1) Background: for the study of highly magnetized neutron stars observed as magnetars, and to quantify the effect of this intense magnetic field onto the star's structure and shape which can be particularly relevant for the study of emission of continuous gravitational waves, both numerical and perturbative approaches have been developed. (2) Methods: we compare these two approaches in General Relativity with the limitation to the case where the magnetic field has a purely poloidal structure. The perturbative one (Konno-99) assumes that the deformation induced by the magnetic field is small and that this field arises only from dipole currents. The full numerical one is based on the library LORENE. (3) Results: we have used both approaches to compute the magnetic field distribution and the deformation of the star, varying the value of the magnetic field at the pole, the compactness of the star and its equation of state. (4) Conclusions: whereas the perturbative approach breaks down for very high polar magnetic field values (typically above a few times $10^{16}$ G), it gives very good results for observed values, even in magnetars. On the contrary, the numerical code exhibits resolution problems for relatively low magnetic field values (typically $10^{10}$ G), which translates into imprecise computation of the star's deformation and mass quadrupole moment.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper compares a perturbative approach (Konno et al. 1999) assuming small deformations from dipole currents and a fully numerical method based on the LORENE library for neutron stars with purely poloidal magnetic fields in general relativity. By systematically varying the polar magnetic field strength, stellar compactness, and equation of state, it concludes that the perturbative method remains accurate for observed magnetar fields but breaks down above a few times 10^16 G, while the numerical code shows resolution problems below ~10^10 G that impair deformation and quadrupole moment calculations.
Significance. If the validity ranges hold under the stated assumptions, the work is significant for providing practical guidance on selecting modeling approaches for highly magnetized neutron stars, particularly in contexts like continuous gravitational wave emission studies. The head-to-head comparison and exploration across compactness and EOS values represent a useful contribution to the field.
major comments (2)
- [Methods] Methods: The restriction to purely poloidal magnetic fields is stated explicitly, yet the Conclusions apply the reported thresholds and the statement of 'very good results... even in magnetars' without discussing or bounding the impact of possible toroidal components, which would alter current distributions, deformations, and quadrupole moments. This assumption is load-bearing for the central applicability claims.
- [Results] Results: The identification of specific breakdown thresholds (perturbative above a few times 10^16 G; numerical below 10^10 G) is presented without error bars, convergence tests, or tabulated data for the magnetic field distributions, deformations, and quadrupole moments across the parameter variations. This weakens the precision of the quantitative claims.
minor comments (1)
- [Abstract] Abstract and Conclusions: The phrases 'very good results' and 'imprecise computation' would benefit from explicit quantitative metrics or tolerance criteria used to define them.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable feedback on our manuscript. We address each major comment in detail below and outline the revisions we will make to improve the clarity and precision of our results.
read point-by-point responses
-
Referee: The restriction to purely poloidal magnetic fields is stated explicitly, yet the Conclusions apply the reported thresholds and the statement of 'very good results... even in magnetars' without discussing or bounding the impact of possible toroidal components, which would alter current distributions, deformations, and quadrupole moments. This assumption is load-bearing for the central applicability claims.
Authors: We thank the referee for highlighting this point. Our comparison is intentionally limited to purely poloidal magnetic fields because the perturbative method of Konno et al. (1999) is derived under the assumption of dipole currents producing poloidal fields, and the LORENE implementation we use is configured for the same. The conclusions regarding the validity for observed magnetar fields are therefore within this framework. However, to address the concern, we will add a paragraph in the Conclusions section explicitly noting that the reported thresholds apply specifically to purely poloidal configurations. We will also mention that the inclusion of toroidal components, which are expected in realistic models, could modify the deformation and quadrupole moment, and that extending the comparison to mixed fields is a natural direction for future work. This will better bound the applicability of our claims. revision: yes
-
Referee: The identification of specific breakdown thresholds (perturbative above a few times 10^16 G; numerical below 10^10 G) is presented without error bars, convergence tests, or tabulated data for the magnetic field distributions, deformations, and quadrupole moments across the parameter variations. This weakens the precision of the quantitative claims.
Authors: We agree that providing more quantitative support for the thresholds would strengthen the manuscript. In the revised version, we will include additional figures or a table presenting the relative differences between the two methods as a function of magnetic field strength, along with results from convergence tests for the LORENE code at various resolutions. Where feasible, we will report estimated uncertainties in the computed deformations and quadrupole moments. The thresholds were identified based on where the discrepancies exceed a certain percentage (e.g., 10%), but we will make this criterion explicit and provide the underlying data points for key cases across different compactnesses and equations of state. revision: yes
Circularity Check
No significant circularity: independent methods compared directly
full rationale
The paper performs a direct numerical comparison between the established perturbative formalism of Konno et al. (1999) and the independent LORENE library under a shared poloidal-field restriction. No parameters are fitted to subsets of the output data and then re-labeled as predictions; the reported breakdown thresholds (above a few times 10^16 G for perturbation, below 10^10 G for numerics) emerge from explicit computation of field distributions, deformations, and quadrupole moments rather than from any self-definitional loop or self-citation chain. Both methods are treated as external inputs whose results are cross-checked, satisfying the criterion of self-containment against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption General relativity governs the spacetime and hydrostatic equilibrium of the neutron star.
- ad hoc to paper The magnetic field is purely poloidal.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Methods: we compare these two approaches in General Relativity with the limitation to the case where the magnetic field has a purely poloidal structure. ... Conclusions: whereas the perturbative approach breaks down for very high polar magnetic field values (typically above a few times 10^{16} G)
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The perturbative one (Konno-99) assumes that the deformation induced by the magnetic field is small and that this field arises only from dipole currents.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Deformation of relativistic magnetized stars
K. Konno, T. Obata, and Y. Kojima, Astron. Astrophys.352, 211 (1999), arXiv:gr-qc/9910038 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 1999
-
[2]
LORENE: Spectral methods differential equations solver,
E. Gourgoulhon, P . Grandcl ´ement, J.-A. Marck, J. Novak, and K. Taniguchi, “LORENE: Spectral methods differential equations solver,” Astrophysics Source Code Library, record ascl:1608.018 (2016), ascl:1608.018. 24
work page 2016
-
[3]
Equations of state for supernovae and compact stars
M. Oertel, M. Hempel, T. Kl ¨ahn, and S. Typel, Rev. Mod. Phys.89, 015007 (2017), arXiv:1610.03361 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[4]
S. A. Olausen and V . M. Kaspi, Astrophys. J. Suppl.212, 6 (2014), arXiv:1309.4167 [astro- ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2014
- [5]
-
[6]
The illusion of neutron star magnetic field estimates
J. P ´etri, Month. Not. Roy. Astron. Soc.485, 4573 (2019), arXiv:1903.01528 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2019
- [7]
-
[8]
S. Bonazzola and E. Gourgoulhon, Astron. Astrophys.312, 675 (1996), arXiv:astro- ph/9602107 [astro-ph]
- [9]
- [10]
- [11]
- [12]
-
[13]
S. Dall’Osso, B. Giacomazzo, R. Perna, and L. Stella, Astrophys. J.798, 25 (2014)
work page 2014
- [14]
-
[15]
A. Mastrano, A. Melatos, A. Reisenegger, and T. Akg ¨un, Mon. Not. Roy. Astron. Soc.417, 2288–2299 (2011)
work page 2011
-
[16]
A. Mastrano, A. G. Suvorov, and A. Melatos, Mon. Not. Roy. Astron. Soc.447, 3475–3485 (2015)
work page 2015
- [17]
-
[18]
G. Pagliaro, M. A. Papa, J. Ming, and M. Muratore, Mon. Not. Roy. Astron. Soc.540, 1006 (2025), https://academic.oup.com/mnras/article-pdf/540/1/1006/63140705/staf774.pdf
work page 2025
-
[19]
Deformations of Accreting Neutron Star Crusts and Gravitational Wave Emission
G. Ushomirsky, C. Cutler, and L. Bildsten, Month. Not. Roy. Astron. Soc.319, 902 (2000), arXiv:astro-ph/0001136 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2000
-
[20]
C. J. Horowitz and K. Kadau, Phys. Rev. Lett.102, 191102 (2009), arXiv:0904.1986 [astro- ph.SR]
work page internal anchor Pith review Pith/arXiv arXiv 2009
- [21]
-
[22]
F. Gittins and N. Andersson, Mon. Not. Roy. Astron. Soc.507, 116 (2021), arXiv:2105.06493 [astro-ph.HE]
-
[23]
P . Adhikari, M. Ammon, S. S. Avancini, A. Ayala, A. Bandyopadhyay, D. Blaschke, F. L. Braghin, P . Buividovich, R. P . Cardoso, C. Cartwright, and et al., Prog. Part. Nuc. Phys.146, 104199 (2026), arXiv:2412.18632 [nucl-th]
-
[24]
ROTATING NEUTRON STAR MODELS WITH MAGNETIC FIELD
M. Bocquet, S. Bonazzola, E. Gourgoulhon, and J. Novak, Astron. Astrophys.301, 757 (1995), arXiv:gr-qc/9503044 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 1995
- [25]
-
[26]
Relativistic stars with purely toroidal magnetic fields
K. Kiuchi and S. Yoshida, Phys. Rev. D78, 044045 (2008), arXiv:0802.2983 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[27]
Structure and deformations of strongly magnetized neutron stars with twisted torus configurations
R. Ciolfi, V . Ferrari, and L. Gualtieri, Mon. Not. Roy. Astron. Soc.406, 2540 (2010), arXiv:1003.2148 [astro-ph.SR]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[28]
Twisted-torus configurations with large toroidal magnetic fields in relativistic stars
R. Ciolfi and L. Rezzolla, Mon. Not. Roy. Astron. Soc.435, L43 (2013), arXiv:1306.2803 [astro- ph.SR]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[29]
J. B. Hartle, Astrophys. J.150, 1005 (1967)
work page 1967
-
[30]
J. B. Hartle and K. S. Thorne, Astrophys. J.153, 807 (1968)
work page 1968
- [31]
-
[32]
S. Bonazzola, E. Gourgoulhon, M. Salgado, and J. A. Marck, Astron. Astrophys.278, 421 (1993)
work page 1993
-
[33]
Relativistic models of magnetars: structure and deformations
A. Colaiuda, V . Ferrari, L. Gualtieri, and J. A. Pons, Mon. Not. Roy. Astron. Soc.385, 2080 (2008), arXiv:0712.2162 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2080
-
[34]
Grad-Shafranov equation in noncircular stationary axisymmetric spacetimes
K. Ioka and M. Sasaki, Phys. Rev. D67, 124026 (2003), arXiv:gr-qc/0302106 [gr-qc]. 26
work page internal anchor Pith review Pith/arXiv arXiv 2003
-
[35]
Relativistic stars with poloidal and toroidal magnetic fields and meridional flow
K. Ioka and M. Sasaki, Astrophys. J.600, 296 (2004), arXiv:astro-ph/0305352 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2004
-
[36]
S. Chandrasekhar and J. C. Miller, Mon. Not. Roy. Astron. Soc.167, 63 (1974)
work page 1974
- [37]
-
[38]
Revising the multipole moments of numerical spacetimes, and its consequences
G. Pappas and T. A. Apostolatos, Phys. Rev. Lett.108, 231104 (2012), arXiv:1201.6067 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[39]
J. L. Friedman and N. Stergioulas, Rotating Relativistic Stars (Cambridge University Press, 2013)
work page 2013
-
[40]
R. Prix, J. Novak, and G. L. Comer, Phys. Rev. D71, 1 (2005)
work page 2005
-
[41]
Gourgoulhon, Lecture Notes in Physics, Lecture Notes in Physics, Vol
E. Gourgoulhon, Lecture Notes in Physics, Lecture Notes in Physics, Vol. 846 (Springer Verlag, Berlin, 2012)
work page 2012
-
[42]
P . Grandcl ´ement and J. Novak, Living Rev. Relativity12, 1 (2009), http://www.livingreviews.org/lrr-2009-1
work page 2009
-
[43]
P . Grandcl´ement, S. Bonazzola, E. Gourgoulhon, and J.-A. Marck, J. Comput. Phys.170, 231 (2001)
work page 2001
-
[44]
E. Gourgoulhon and S. Bonazzola, Class. Quantum Grav.11, 443 (1994)
work page 1994
- [45]
- [46]
-
[47]
Equation of state and thickness of the inner crust of neutron stars
F. Grill, H. Pais, C. Provid ˆencia, I. Vida˜na, and S. S. Avancini, Phys. Rev. C90, 045803 (2014), arXiv:1404.2753 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2014
- [48]
- [49]
-
[50]
Hyperons in neutron star matter within relativistic mean-field models
M. Oertel, C. Provid ˆencia, F. Gulminelli, and A. R. Raduta, J. Phys. G42, 075202 (2015), arXiv:1412.4545 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[51]
Composition and thermodynamics of nuclear matter with light clusters
S. Typel, G. Ropke, T. Klahn, D. Blaschke, and H. H. Wolter, Phys. Rev. C81, 015803 (2010), arXiv:0908.2344 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[52]
The equation of state for nucleon matter and neutron star structure
A. Akmal, V . R. Pandharipande, and D. G. Ravenhall, Phys. Rev. C58, 1804 (1998), arXiv:nucl-th/9804027
work page internal anchor Pith review Pith/arXiv arXiv 1998
- [53]
- [54]
- [55]
-
[56]
Ddfgos(apr) cold neutron star eos,
P . Davis et al., “Ddfgos(apr) cold neutron star eos,” (2025),https://compose.obspm.fr/ eos/327[Accessed: 2026-08-04]
work page 2025
-
[57]
Gppva(dd2) cold neutron star eos,
H. Pais et al., “Gppva(dd2) cold neutron star eos,” (2023),https://compose.obspm.fr/eos/ 217, [Accessed: 2026-08-04]
work page 2023
-
[58]
Oos(dd2-frg) cold neutron star eos,
K. Otto et al., “Oos(dd2-frg) cold neutron star eos,” (2021),https://compose.obspm.fr/ eos/197, [Accessed: 2026-08-04]
work page 2021
-
[59]
Opgr(gm1y5) cold neutron star eos,
M. Oertel et al., “Opgr(gm1y5) cold neutron star eos,” (2015),https://compose.obspm.fr/ eos/65, [Accessed: 2026-08-04]
work page 2015
- [60]
-
[61]
A. G. Pili, N. Bucciantini, and L. Del Zanna, Month. Not. Roy. Astron. Soc.439, 3541 (2014), arXiv:1401.4308 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2014
- [62]
-
[63]
Consistent neutron star models with magnetic field dependent equations of state
D. Chatterjee, T. Elghozi, J. Novak, and M. Oertel, Mon. Not. R. Astron. Soc.447, 3785 (2015), arXiv:1410.6332 [astro-ph.HE]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[64]
A. M. Raghoonundun, Exact Solutions for Compact Objects in General Relativity, Ph.D. thesis, University of Calgary (Canada) (2016), arXiv:arXiv:1604.08930 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[65]
B. Reina and R. Vera, Class. Quantum Grav.32, 155008 (2015), arXiv:1412.7083 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[66]
L. Rezzolla and B. J. Ahmedov, Mon. Not. R. Astron. Soc.352, 1161 (2004). 28
work page 2004
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.