On the Quasi-Periodic Oscillations of Magnetars

We study torsional Alfv'en oscillations of magnetars, i.e., neutron stars with a strong magnetic field. We consider the poloidal and toroidal components of the magnetic field and a wide range of equilibrium stellar models. We use a new coordinate system (X,Y), where $X=\sqrt{a_1} \sin \theta$, $Y=\sqrt{a_1}\cos \theta$ and $a_1$ is the radial component of the magnetic field. In this coordinate system, the 1+2-dimensional evolution equation describing the quasi-periodic oscillations, QPOs, see Sotani et al. (2007), is reduced to a 1+1-dimensional equation, where the perturbations propagate only along the Y-axis. We solve the 1+1-dimensional equation for different boundary conditions and open magnetic field lines, i.e., magnetic field lines that reach the surface and there match up with the exterior dipole magnetic field, as well as closed magnetic lines, i.e., magnetic lines that never reach the stellar surface. For the open field lines, we find two families of QPOs frequencies; a family of “lower” QPOs frequencies which is located near the X-axis and a family of “upper” frequencies located near the Y-axis. According to Levin (2007), the fundamental frequencies of these two families can be interpreted as the turning points of a continuous spectrum. We find that the upper frequencies are constant multiples of the lower frequencies with a constant equaling 2n+1. For the closed lines, the corresponding factor is n+1 . By these relations, we can explain both the lower and the higher observed frequencies in SGR 1806-20 and SGR 1900+14.

💡 Research Summary

The paper investigates torsional Alfvén oscillations in magnetars—highly magnetized neutron stars—by extending previous work on quasi‑periodic oscillations (QPOs) to include both poloidal and toroidal magnetic field components across a broad set of equilibrium stellar models. The authors introduce a novel coordinate transformation defined by X = √a₁ sin θ and Y = √a₁ cos θ, where a₁ is the radial component of the magnetic field. In this (X, Y) system the original 1 + 2‑dimensional evolution equation governing the QPOs collapses to a 1 + 1‑dimensional wave equation in which perturbations propagate solely along the Y‑axis. This reduction dramatically simplifies the numerical treatment while preserving the essential physics of Alfvén wave propagation along magnetic field lines.

Two distinct families of magnetic field lines are considered: open lines that intersect the stellar surface and match onto an external dipole field, and closed lines that remain confined within the star. For each family the authors solve the reduced wave equation under appropriate boundary conditions (fixed, free, or reflective at the stellar surface and at the magnetic axis). The solutions reveal two characteristic QPO frequency families for the open lines. “Lower” QPOs are localized near the X‑axis (i.e., near the magnetic equator), while “upper” QPOs concentrate near the Y‑axis (i.e., near the magnetic pole). According to Levin (2007), the fundamental frequencies of these families correspond to turning points of a continuous Alfvén spectrum. The key result is that the upper‑frequency family is an integer multiple of the lower‑frequency family, with the multiplier given by 2 n + 1 (n = 0, 1, 2, …). For the closed‑field‑line case the multiplier reduces to n + 1. These simple proportionalities emerge directly from the eigenvalue structure of the 1 + 1‑dimensional equation and are confirmed by extensive numerical experiments.

The authors then confront their theoretical frequency ladders with the observed QPOs from the soft‑gamma‑repeaters SGR 1806‑20 and SGR 1900+14. Both sources exhibit a rich set of frequencies ranging from a few tens of hertz up to several hundred hertz. By assigning the lowest observed peaks to the lower‑frequency family and applying the 2 n + 1 scaling, the higher observed peaks fall naturally into the upper‑frequency family. Likewise, the closed‑line scaling (n + 1) accounts for additional frequencies that do not fit the open‑line pattern. This unified interpretation demonstrates that a single magnetar model, equipped with realistic mixed poloidal‑toroidal fields, can reproduce the full spectrum of observed QPOs without invoking ad‑hoc mode selection.

Beyond the immediate astrophysical application, the paper makes several methodological contributions. The coordinate transformation effectively decouples the angular dependence of the magnetic field, turning a computationally intensive 2‑D problem into a tractable 1‑D wave propagation problem. This approach could be extended to include rotation, non‑linear coupling, or crust‑core interaction in future work. Moreover, the explicit identification of turning‑point frequencies as the origin of discrete QPO families provides a clear physical picture linking the continuous Alfvén spectrum to the discrete signals detected by X‑ray telescopes.

In summary, the study presents a robust, analytically transparent, and numerically efficient framework for modeling magnetar torsional Alfvén QPOs. By demonstrating that the observed QPO frequencies in SGR 1806‑20 and SGR 1900+14 follow simple integer‑multiple relations derived from the underlying magnetic geometry, the authors bridge the gap between theoretical magnetohydrodynamic spectra and real astrophysical data, offering a compelling explanation for the rich phenomenology of magnetar oscillations.