Magnetic field dynamics in isolated neutron stars with an external dipole field

Through long-term numerical-relativity simulations, this study demonstrates that isolated neutron stars with initial mixed magnetic fields dynamically relax into stable configurations where the toroidal component contributes less than 10% of the total magnetic energy, a process driven by Tayler instabilities and gravitational-wave emission that constrains the long-term evolution of pulsar and magnetar magnetic fields.

The Gist

Imagine a neutron star as a cosmic dead star, incredibly dense and small, but spinning with a magnetic force so strong it could rip a credit card apart from a million miles away. For a long time, scientists have wondered: How does the magnetic field inside this star actually hold together?

Is it a simple bar magnet? A twisted knot? Or something else entirely?

This paper by Capobianco, Cook, and their team uses supercomputer simulations to answer that question. They treated the neutron star like a giant, invisible ball of fluid and watched how its magnetic field behaved over time. Here is what they found, explained simply:

1. The Setup: A Tangled Mess

The scientists started their simulation with a neutron star that had a strong, simple magnetic field on the outside (like a standard bar magnet), but a messy, complex mix of fields inside. They specifically tested what happens if the inside is dominated by a "twisted" magnetic component (called the toroidal field), which is like a rubber band wrapped tightly around the star's equator.

They tested different scenarios, some where the twist was weak and some where it was extremely strong (up to 80% of the total magnetic energy).

2. The Chaos: The "Sausage" and the "Kink"

As soon as they turned on the simulation, the magnetic field didn't stay calm. It started to wiggle and break apart. The paper describes two main ways the field tried to tear itself apart:

• The "Sausage" Instability: Imagine a long, thin tube of magnetic force. Suddenly, it pinches in the middle and bulges out at the ends, looking like a string of sausages.
• The "Kink" Instability: Imagine twisting a rubber band until it snaps and kinks over on itself.

These instabilities caused the magnetic field lines to tangle, twist, and churn violently, creating a chaotic storm inside the star.

3. The Calm After the Storm: Finding a Stable Shape

Here is the most important discovery: The chaos didn't last forever.

After about 150 milliseconds (a blink of an eye in cosmic time), the magnetic field stopped fighting itself. It settled down into a new, stable shape.

• The Result: The star didn't keep the massive, twisted "rubber band" it started with. Instead, it relaxed into a mixed configuration.
• The Ratio: In this final, stable state, the "twisted" part of the magnetic field shrank dramatically. It ended up contributing only about 0.5% to 10% of the total magnetic energy. The rest was a more standard, flowing field.

Think of it like a child playing with a tangled ball of yarn. At first, they pull and twist it into a huge, messy knot. But eventually, they let go, and the yarn settles into a neat, manageable ball. The neutron star's magnetic field does the same thing: it untangles itself until it finds a stable, mixed shape that won't fall apart.

4. The "Leak" and the Wave

During this process, two other things happened:

• The Leak: Because the magnetic field was so strong, some of the "twisted" energy actually leaked out of the star's surface and into the space around it, like steam escaping a pressure cooker. This helped the star calm down faster.
• The Rumble: As the magnetic field rearranged itself, it made the star vibrate. These vibrations sent out ripples in space and time called gravitational waves. The paper detected these waves, noting that the specific "song" the star sang changed as the magnetic field settled down.

5. Why This Matters

The paper concludes that no matter how messy or twisted the magnetic field starts out inside a neutron star, it naturally evolves toward a specific, stable "sweet spot." It won't stay a chaotic mess, and it won't stay a purely twisted knot. It will always settle into a mixed state where the twisted part is small but necessary for stability.

This finding helps astronomers understand:

• How long these magnetic fields can last.
• Why pulsars (spinning neutron stars) emit light the way they do.
• What kind of "ripples" in space (gravitational waves) we should expect to detect from these stars.

In short: The universe seems to have a rule for neutron stars: if you twist their magnetic fields too much, they will eventually untangle themselves just enough to find a comfortable, stable balance.

Technical Summary

Technical Summary: Magnetic Field Dynamics in Isolated Neutron Stars with an External Dipole Field

Problem Statement
While neutron stars (NSs) possess the strongest magnetic fields in the universe, the structure and stability of their internal magnetic configurations remain poorly constrained by observation. Observational data suggests a predominantly dipolar external field, yet the internal topology is theoretically expected to be a complex, mixed poloidal-toroidal geometry to ensure stability against kink and varicose instabilities. Previous numerical studies have often focused on purely interior fields or employed Newtonian approximations, leaving a gap in understanding how an external dipole field interacts with mixed internal configurations over long timescales in a fully relativistic framework. Specifically, it remains unclear whether these systems relax into stable equilibrium states and how the presence of an external field influences the toroidal-to-poloidal energy ratio in the final configuration.

Methodology
The authors perform long-term, nonlinear General-Relativistic Magnetohydrodynamic (GRMHD) simulations of isolated, non-rotating neutron stars using the AthenaK code. The simulations couple the ideal GRMHD equations with the Z4c formulation of the 3+1 Einstein equations to evolve the spacetime dynamically. For comparison, identical setups are also run in the Cowling approximation, where the spacetime metric is fixed to that of a static Tolman-Oppenheimer-Volkoff (TOV) star.

• Initial Conditions: The models feature a static TOV star with a Gamma-law equation of state ($\Gamma=2$). The magnetic field is initialized with a continuous, axisymmetric configuration consisting of an external dipole field ($B_p = 10^{15}$ G) and a mixed poloidal-toroidal field confined to the interior.
• Parameter Space: Six models are simulated: a non-magnetized reference and five magnetized models where the initial toroidal-to-total magnetic energy ratio ranges from 0% to 80% ($B_{t0}$ to $B_{t80}$).
• Numerical Setup: The computational domain extends to $\sim 470$ km (approx. $40R$) to minimize boundary effects, with the innermost mesh resolution at $\Delta x = 230$ m. Simulations evolve for 150 ms (up to 200 ms for Cowling runs).
• Analysis: The study tracks the evolution of magnetic field magnitude, energy decomposition (poloidal vs. toroidal, interior vs. exterior), mode analysis via spherical harmonic decomposition of density and radial magnetic field, and gravitational wave (GW) emission via the Newman-Penrose scalar $\Psi_4$.

Key Contributions and Results

1. Relaxation to Stable Mixed Configurations:
   The simulations demonstrate that initially unstable magnetic configurations, regardless of the initial toroidal energy fraction (up to 80%), evolve through Tayler instabilities (varicose and kink modes) and relax toward a dynamically stable, mixed poloidal-toroidal geometry.

   • Final State: In the final quasi-stationary state, the toroidal component contributes only $\lesssim 10\%$ of the total magnetic energy, both in the interior and the exterior. This is significantly lower than the initial 80% in the $B_{t80}$ model.
   • Timescale: This relaxation occurs within approximately one Alfvén crossing time ($\tau_A \sim 240$ ms) following the saturation of instabilities.

2. Role of External Field and Flux Leakage:
   A critical finding is the role of the external dipole field and the stellar atmosphere. Unlike previous studies where toroidal fields were confined to the interior, the presence of an external field allows toroidal magnetic flux to leak into the atmosphere.

   • Energy Dissipation: This leakage, combined with gravitational-wave emission, drives a more efficient decay of magnetic energy. The final toroidal energy fraction (0.5%–10%) is lower than in simulations where the field is strictly interior-confined (where fractions of $\sim 15\%$ were observed).
   • Comparison with Cowling: Fully dynamical GR simulations show a delayed onset of instabilities and energy decay compared to Cowling runs, highlighting the importance of spacetime dynamics in energy redistribution.

3. Energetics and Electromagnetic Emission:

   • Magnetic Energy: The magnetic energy initially increases by 20–50% during the rearrangement phase before decaying. By 150 ms, the internal magnetic energy drops to $\sim 23\%$ of its initial value in GR runs.
   • Luminosity: The electromagnetic luminosity ($L_{EM}$) exhibits an initial burst, followed by a decrease and a secondary small increase around 40 ms (coinciding with toroidal flux leakage). It settles into a quasi-stationary phase with $L_{EM} \sim 10^{46}$ erg/s.

4. Gravitational Wave Signatures:
   The magnetically driven oscillations produce GW signals dominated by the $(\ell, m) = (2,0)$ and $(2,2)$ modes.

   • Spectral Features: Fourier analysis reveals that magnetized runs exhibit an additional peak at early times (0–25 ms) corresponding to a nonlinear coupling between the fundamental pressure quadrupole mode ($2f$) and the first quasiradial overtone ($H_1$). This coupling is more prominent in models with higher initial toroidal fields.
   • Decay: The GW signal amplitude tracks the evolution of magnetic and kinetic energies, showing peaks during instability saturation and subsequent decay.

Significance and Claims
The paper claims that these results provide the first long-term GRMHD investigation of isolated neutron stars with external dipole fields and mixed internal configurations. The primary significance lies in the demonstration that long-lived neutron star magnetic fields are strongly constrained toward stable mixed configurations where the toroidal component is a minor fraction of the total energy ($\lesssim 10\%$).

The authors assert that this convergence to a specific equilibrium state, largely insensitive to initial conditions within the explored class, suggests that strongly magnetized NSs naturally evolve toward a limited set of equilibria under ideal GRMHD evolution. These findings have direct implications for:

• Pulsar Emission Models: The final field geometry influences the magnetosphere structure.
• Magnetar Evolution: The decay of magnetic energy and the specific equilibrium state inform models of magnetar activity.
• Gravitational Wave Interpretation: The specific spectral imprints (nonlinear mode couplings) and the timescales of magnetic rearrangement provide templates for interpreting GW signals from magnetized remnants or magnetar flares.

The authors note limitations, including the absence of resistive MHD effects (reconnection is purely numerical) and the use of a simplified atmosphere, which may affect the quantitative fraction of surviving toroidal fields. However, the qualitative conclusion regarding the stability of mixed configurations and the suppression of the toroidal component remains robust across the tested scenarios.