Magnetized neutron stars: perturbative versus fully-numerical approaches

This paper compares perturbative and fully numerical approaches for modeling magnetized neutron stars with purely poloidal fields, finding that while the perturbative method (Konno-99) is accurate for observed magnetar fields but fails at extreme strengths ($>10^{16}$ G), the fully numerical LORENE code struggles with resolution issues at lower field strengths ($\sim10^{10}$ G).

The Gist

Imagine a neutron star as a cosmic city, incredibly dense and heavy, where the laws of physics are pushed to their absolute limits. Now, imagine this city is being squeezed and stretched by an invisible, super-powerful magnetic force field. This is the world of magnetars, a type of neutron star with magnetic fields so strong they could wipe a credit card from halfway across the galaxy.

Scientists want to understand exactly how these magnetic fields warp the shape of the star. Why? Because if a star is perfectly round, it spins silently. But if the magnetic field squishes it into an egg shape, it might wobble as it spins, sending out ripples in space-time called gravitational waves. Detecting these waves is like listening for a whisper in a hurricane; we need to know exactly what the "whisper" should sound like to find it.

To figure this out, scientists have developed two different ways of doing the math: a simplified shortcut (the perturbative approach) and a brute-force super-computation (the fully numerical approach). This paper is like a referee stepping in to see which method is better and when.

The Two Methods: A Map vs. A 3D Scan

1. The Perturbative Approach (The "Small Stretch" Map)
Think of this method as drawing a map of a slightly bumpy road. It starts with a perfect, smooth sphere (the star without a magnetic field) and then asks, "What happens if we add a tiny amount of magnetic stretch?"

• The Assumption: It assumes the magnetic field is simple (like a bar magnet) and that the star doesn't change shape very much.
• The Analogy: It's like calculating how much a trampoline sags when you place a single bowling ball on it. It works great for small weights because the math stays simple and linear.

2. The Fully Numerical Approach (The "Full 3D Scan")
This method doesn't assume the star is round to begin with. It builds the star from scratch, calculating every single point of pressure and magnetic force simultaneously, allowing the star to twist, squash, and deform as much as it wants.

• The Assumption: It lets the physics speak for itself without forcing the star to stay round.
• The Analogy: This is like using a high-end 3D scanner to model a trampoline with a giant boulder on it. It captures every wrinkle and dip, but it requires a massive amount of computing power and is very sensitive to tiny errors in the calculation.

The Showdown: Who Wins?

The authors ran both methods side-by-side, testing them with different star sizes and different types of "star soup" (equations of state). Here is what they found:

Scenario A: The "Normal" Magnetar (Low to Medium Magnetic Fields)

• The Result: Both methods agree perfectly.
• The Takeaway: For the magnetic fields we actually observe in the universe (even in the strongest magnetars), the "Small Stretch" map is just as accurate as the "Full 3D Scan." The shortcut works! You don't need a supercomputer to get the right answer for the stars we know today.

Scenario B: The "Super-Magnetar" (Extremely High Magnetic Fields)

• The Result: The "Small Stretch" map breaks down.
• The Takeaway: If the magnetic field gets crazy strong (above a few times $10^{16}$ Gauss), the star deforms so much that the "tiny stretch" assumption is no longer true. The shortcut fails, and you must use the heavy-duty 3D scan to get the right answer.

Scenario C: The "Ghost" Problem (Very Low Magnetic Fields)

• The Result: Surprisingly, the "Full 3D Scan" struggles here.
• The Takeaway: When the magnetic field is weak, the star is almost perfectly round. The 3D scanner tries to calculate the difference between "perfectly round" and "almost perfectly round." Because these numbers are so close, the computer gets confused by tiny rounding errors (like trying to measure the thickness of a hair by subtracting two huge numbers). The "Small Stretch" map, which was built to handle these small changes, is actually more accurate for weak fields.

The Verdict

The paper concludes with a clear rule of thumb for astronomers hunting for gravitational waves:

1. For the stars we see today: The simple, fast, "perturbative" method is sufficient. It gives accurate results for the magnetic fields we actually measure, making it much easier to model these stars and predict the gravitational waves they might emit.
2. For the extreme edge cases: If we ever encounter a star with a magnetic field far stronger than anything we've seen yet, we will need the complex, numerical method.
3. For the very weak fields: If you are looking at very subtle deformations, the simple method is actually more precise because the complex method gets tripped up by computer math errors.

In short, for the current "cosmic city" we are observing, the shortcut is not just a good guess—it's the right tool for the job. The heavy machinery is only needed if we discover a monster star that breaks the rules of our current observations.

Technical Summary

Technical Summary: Magnetized Neutron Stars – Perturbative versus Fully-Numerical Approaches

Problem Statement
The study of highly magnetized neutron stars (magnetars) requires accurate modeling of how intense magnetic fields ($B \sim 10^{14} - 10^{15}$ G, potentially up to $10^{18}$ G theoretically) deform stellar structure. These deformations are critical for predicting the emission of continuous gravitational waves, a target for current and future detectors (e.g., LIGO-Virgo-Kagra, DECIGO, Einstein Telescope). While both perturbative and fully numerical methods exist to model these systems, their respective domains of validity and accuracy had not been systematically compared. The central question is whether the simpler, computationally cheaper perturbative approach is sufficient for observed magnetars, or if the more complex fully numerical approach is strictly necessary to capture non-linear general-relativistic effects and higher-order multipoles.

Methodology
The authors compare two distinct theoretical frameworks within General Relativity, restricting the analysis to non-rotating stars with purely poloidal magnetic fields and axial symmetry:

1. Perturbative Approach: Based on the work of Konno et al. [1], this method treats the magnetic field as a small perturbation on a spherically symmetric, non-magnetized background (Tolman-Oppenheimer-Volkoff solution). It assumes the magnetic field arises solely from dipole currents ($\ell=1$) and that the induced deformation is small. The equations are solved using a second-order Runge-Kutta integrator, with careful treatment of the origin and surface matching conditions.
2. Fully-Numerical Approach: Utilizing the magstar code and the LORENE library [2], this method solves the coupled Einstein-Maxwell and hydrostatic equilibrium equations without assuming small deformations. It employs spectral methods (Chebyshev polynomials in radius, Fourier series in angle) and multi-domain decomposition to handle the infinite domain. The magnetic field is assumed purely poloidal, but higher multipoles are naturally allowed.

The comparison focuses on gauge-independent physical quantities: the magnetic field distribution (central, equatorial, and polar values), the mass quadrupole moment ($Q$), and the surface ellipticity ($\epsilon_{surf}$). The study varies the polar magnetic field strength ($B_{pole}$), the stellar compactness ($C = M/R$), and the Equation of State (EoS), testing four models: DDFGOS(APR), GPPVA(DD2), OOS(DD2-FRG), and OPGR(GM1Y5).

Key Results

• Magnetic Field Distribution: For polar magnetic fields up to approximately $B_{pole} \sim 10^{16}$ G, the perturbative and numerical approaches yield nearly identical results for the magnetic field strength at the center and equator. Significant deviations (up to 100%) only appear near the maximum magnetic field a star can support ($B_{pole} > 10^{17}$ G), where the assumptions of a spherical background and pure dipole currents in the perturbative model break down.
• Deformation Metrics (Quadrupole and Ellipticity): The behavior of deformation indicators differs significantly from the magnetic field results.
  • Low Field Regime: For relatively low magnetic fields ($B_{pole} \lesssim 10^{15}$ G for quadrupole, $\lesssim 10^{10}$ G for ellipticity), the perturbative approach is more accurate than the fully numerical one. The numerical code suffers from resolution issues when computing these quantities, as they are derived from the difference between nearly equal numbers (e.g., asymptotic metric coefficients), leading to numerical noise.
  • High Field Regime: As $B_{pole}$ increases beyond $\sim 5 \times 10^{16}$ G, the perturbative approach fails due to non-linear general-relativistic effects and the emergence of higher multipoles.
• Thresholds and Dependencies: The authors define threshold values ($B^5_{pole}$ and $B^{50}_{pole}$) where the relative difference in ellipticity exceeds 5% and 50%, respectively. These thresholds increase with stellar compactness (more compact stars require stronger fields to deform significantly). The dependence on the EoS is present but weak; softer EoSs allow the perturbative approach to remain valid up to slightly higher magnetic fields than stiffer ones.

Significance and Conclusions
The paper concludes that for the modeling of all currently observed neutron stars and magnetars (where $B_{pole}$ is typically below $10^{16}$ G), the perturbative approach is sufficiently accurate and computationally efficient. This finding is particularly relevant for gravitational wave searches, as it suggests that simpler models can reliably estimate the deformations induced by magnetic fields without the computational cost of full numerical relativity.

The authors note that the numerical approach currently exhibits limitations in precision at low magnetic fields for deformation metrics, a technical issue that requires improvement in how quadrupole moments and ellipticities are extracted from the code. Furthermore, the study is limited to non-rotating, purely poloidal fields; future work must address mixed poloidal-toroidal configurations and the inclusion of magnetic field effects within the EoS itself to fully capture realistic magnetar physics.

Caveats
The authors emphasize that this is the first systematic comparison of these two specific approaches. The results are modest in scope, applying strictly to non-rotating stars with poloidal fields. The validity of the perturbative approach is confirmed only for the range of magnetic fields currently observed in nature, not necessarily for the theoretical maximums ($10^{18}$ G) where non-linear effects dominate.