Magneto-elastic oscillations of neutron stars with dipolar magnetic fields

By means of two dimensional, general-relativistic, magneto-hydrodynamical simulations we investigate the oscillations of magnetized neutron star models (magnetars) including the description of an extended solid crust. The aim of this study is to understand the origin of the QPOs observed in the giant flares of SGRs. We confirm the existence of three different regimes: (a) a weak magnetic field regime B<5 x 10^13 G, where crustal shear modes dominate the evolution; (b) a regime of intermediate magnetic fields 5 x 10^13 G<B< 10^15 G, where Alfv'en QPOs are mainly confined to the core of the neutron star and the crustal shear modes are damped very efficiently; and (c) a strong field regime B>10^15 G, where magneto-elastic oscillations reach the surface and approach the behavior of purely Alfv'en QPOs. When the Alfv'en QPOs are confined to the core of the neutron star, we find qualitatively similar QPOs as in the absence of a crust. The lower QPOs associated with the closed field lines of the dipolar magnetic field configuration are reproduced as in our previous simulations without crust, while the upper QPOs connected to the open field lines are displaced from the polar axis. Additionally, we observe a family of edge QPOs. Our results do not leave much room for a crustal-mode interpretation of observed QPOs in SGR giant flares, but can accommodate an interpretation of these observations as originating from Alfv'en-like, global, turning-point QPOs in models with dipolar magnetic field strengths in the narrow range of 5 x 10^15 G < B < 1.4 x 10^16 G. This range is somewhat larger than estimates for magnetic field strengths in known magnetars. The discrepancy may be resolved in models including a more complicated magnetic field structure or with models taking superfluidity of the neutrons and superconductivity of the protons in the core into account.

💡 Research Summary

This paper presents a comprehensive study of magneto‑elastic oscillations in magnetized neutron stars (magnetars) using two‑dimensional, general‑relativistic magnetohydrodynamic (GRMHD) simulations that incorporate an extended solid crust. The authors aim to explain the quasi‑periodic oscillations (QPOs) observed in the tails of giant flares from soft gamma‑ray repeaters (SGRs).

The theoretical framework starts from the 3+1 split of general relativity, combining the ideal‑MHD stress‑energy tensor with an elastic contribution based on Carter‑Samuelsson’s formalism. In the linear, small‑amplitude limit appropriate for torsional shear oscillations, the elastic terms reduce to a shear modulus μS and a shear tensor Σμν, which are coupled to the conserved GRMHD variables. The resulting hyperbolic system is solved with a high‑resolution shock‑capturing (HRSC) Riemann solver; an alternative linearized wave‑equation approach for the crust is also implemented as a consistency check.

Equilibrium background models are constructed for several realistic equations of state (EoS) and stellar masses, each endowed with a purely dipolar magnetic field. The magnetic field strength at the surface is used to define three distinct regimes:

1. Weak‑field regime (B < 5 × 10¹³ G) – The crustal shear modes dominate. The simulations recover the familiar torsional shear spectrum ranging from tens of Hz up to a few kHz, consistent with earlier crust‑only studies.

2. Intermediate‑field regime (5 × 10¹³ G < B < 10¹⁵ G) – The crustal shear modes are rapidly damped by resonant absorption into the Alfvén continuum in the fluid core. Damping times are ≤ 0.2 s, far shorter than the observed QPO lifetimes, indicating that pure shear modes cannot explain the SGR QPOs in this regime.

3. Strong‑field regime (B > 10¹⁵ G) – The Alfvén oscillations are no longer confined to the core; they couple efficiently through the crust and reach the stellar surface. Two families of global Alfvén QPOs appear: “lower” QPOs associated with turning points of closed field‑line regions near the equator, and “upper” QPOs linked to turning points of open field lines close to the magnetic pole. The upper QPOs have their maximum amplitudes at the surface, making them viable candidates for the observed QPOs. Additionally, a set of “edge” QPOs emerges at the last open field line, and a new upper QPO not seen in crust‑less models is identified.

The authors extend the semi‑analytic short‑wavelength model of Cerdá‑Durán et al. (2009) to compute the Alfvén continuum for the dipolar configuration. The semi‑analytic frequencies match the numerical results with high accuracy, confirming that the QPOs correspond to turning points or edges of the continuum. Frequency ratios of overtones are integer multiples, naturally reproducing observed QPO clusters such as 30 Hz, 92 Hz, and 150 Hz in SGR 1806‑20.

A key outcome is that the observed SGR QPOs can be interpreted as global Alfvén‑like turning‑point QPOs only if the surface magnetic field lies in a relatively narrow window of ≈ 5 × 10¹⁵ G – 1.4 × 10¹⁶ G (for the stiff EoS models examined). This range exceeds typical estimates for known magnetars (≈ 10¹⁴ – 10¹⁵ G). The authors suggest that more complex magnetic topologies (e.g., tangled or toroidal components) or the inclusion of superfluid neutrons and superconducting protons in the core could shift the required field strength, potentially reconciling theory with observations.

In summary, the paper demonstrates that crustal shear modes are quickly damped in realistic magnetar magnetic fields and cannot account for the low‑frequency QPOs. Instead, strong‑field magneto‑elastic oscillations that propagate to the surface and form global Alfvén QPOs provide a plausible explanation for the observed QPO spectrum, provided that the magnetic field strength and internal composition are appropriately modeled.