Recent Research on Neutron Stars
Magnetohydrodynamic simulation of Magnetar Giant Flare

Onset and evolution of the eruption. The magnetosphere initially inflates as the rotating stellar surface twists the magnetic field lines. Following magnetic reconnection in the equatorial plane, the magnetosphere relaxes toward a nearly dipolar configuration, retaining only a small residual toroidal field.

Time evolution of the plasma internal energy in the inner magnetosphere, illustrating the formation of a hot, magnetized fireball.

Onset and evolution of the eruption. The magnetosphere initially inflates as the rotating stellar surface twists the magnetic field lines. Following magnetic reconnection in the equatorial plane, the magnetosphere relaxes toward a nearly dipolar configuration, retaining only a small residual toroidal field.
Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts, K. Chatterjee, A. Philippov, A. M. Beloborodov, K. Parfrey, B. Ripperda, and E. R. Most, 2026
This work, led by postdoctoral fellow Koushik Chatterjee in collaboration with Alexander Philippov (UMd), Andrei Beloborodov (Columbia), Bart Ripperda (CITA), and other collaborators, presents the first relativistic magnetohydrodynamic simulation of a surface shear–driven eruption in a magnetar’s magnetosphere. The study shows how gradual crustal deformation can trigger giant flares with energies exceeding 10⁴⁴ erg.
The simulations show that once the magnetic twist exceeds a critical threshold, the magnetosphere inflates and reconnection occurs in its tail. A large fraction of the dipole energy is released in a giant plasmoid, whose trailing hot plasma may produce the main spike of the gamma-ray flare. Hot plasma trapped in the closed magnetosphere forms a fireball that powers the extended X-ray tail. The ejecta also launches a strong compressive wave—a fast magnetosonic pulse—carrying enough energy to potentially produce fast radio bursts.
Although the simulated flare timescale is shorter than observed, uncertainties in reconnection physics—especially in the highly radiative environment near the magnetar—may substantially influence both the reconnection rate and the flare duration, highlighting an important challenge for future studies.
Radiatively driven two-stream instability as a possible mechanism for radio emission in magnetars

(Left panel): Snapshot of the phase-space distribution function, f(x,px), in the steady state of the one-dimensional simulation. (Right panel): Development of turbulence in the current-carrying region, along with the generation and propagation of electromagnetic waves that can escape through the simulation boundaries.

(Left panel): Snapshot of the phase-space distribution function, f(x,px), in the steady state of the one-dimensional simulation. (Right panel): Development of turbulence in the current-carrying region, along with the generation and propagation of electromagnetic waves that can escape through the simulation boundaries.
Origin of Pulsed Radio Emission from Magnetars, S. Zeng, A. Philippov, J. Juno , A. M. Beloborodov, and E. Popova, 2025
This study proposes a plasma mechanism to explain the unusual radio emission from quiescent magnetars. Led by Shuzhe Zeng, a graduate student at the University of Maryland, together with Alexander Philippov (Maryland) and Andrei Beloborodov (Columbia), the team shows that the emission originates from a closed, twisted bundle of magnetic field lines near the magnetic pole, stretching roughly 15–50 stellar radii. In this region, an electron–positron plasma must carry electric current while also experiencing strong drag from the magnetar’s intense X-ray radiation. This drag locks the plasma into a state that becomes unstable to a two-stream plasma instability. Kinetic simulations using a novel computational code show that the turbulence from this instability generates electromagnetic waves that can escape as radio emission.
The model naturally explains several observed features of magnetar radio signals, including their brightness (~10³⁰ erg/s), wide pulses, strong linear polarization, delayed appearance following X-ray outbursts, and a broad spectrum from GHz frequencies up to nearly 100 GHz. Researchers are currently investigating the plasma mechanisms responsible for producing superluminal waves and their spectral properties, with the goal of comparing them to magnetar observations.
Monster Shock and GHz Precursor Waves in Magnetars

Zoom-in on the shock structure in the steepening zone, with the shock centered on the leading soliton-like structure.

Zoom-in on the shock structure in the steepening zone, with the shock centered on the leading soliton-like structure.
Fast radio bursts as precursor radio emission from monster shocks, A. Vanthieghem and A. Levinson, 2025
Using large-scale particle-in-cell simulations, Amir Levinson and Arno Vanthieghem studied how a fast magnetosonic wave, produced near a magnetar, propagates through the declining magnetic field of the inner magnetosphere, steepens, and eventually breaks. The wave forms a ‘monster shock’ when the electric and magnetic fields reach comparable strengths, B² − E² → 0. After forming the shock, the wave dissipates approximately half of its initial energy.
For the first time, the simulations show that the shock generates a high-frequency precursor wave carrying roughly ~10⁻³ of the dissipated energy. Its spectrum features sharp harmonic peaks in the GHz band, consistent with the frequencies observed in fast radio bursts. Under magnetar conditions, the precursor from the monster shock is strong enough to account for the fast radio bursts observed from the galactic magnetar SGR 1935+2154.
These results provide a self-consistent mechanism linking magnetic wave breaking in magnetars to observable radio transients and open new directions for understanding fast radio bursts from magnetars.
Superconducting Phase Transition in Pulsars Powers Giant Hall Waves

Evolution of the magnetic field at the core–crust interface, shown at 10 kyr (left panel) and 200 kyr (right panel) after the phase transition. Green lines trace the poloidal field, while color indicates the toroidal magnetic-field component. The core field evolves due to the flux tube tension, whereas the crustal field is shaped by Hall drift, Ohmic dissipation, and the elastic back-reaction of the crust.

Evolution of the magnetic field at the core–crust interface, shown at 10 kyr (left panel) and 200 kyr (right panel) after the phase transition. Green lines trace the poloidal field, while color indicates the toroidal magnetic-field component. The core field evolves due to the flux tube tension, whereas the crustal field is shaped by Hall drift, Ohmic dissipation, and the elastic back-reaction of the crust.
Giant Hall Waves Launched by Superconducting Phase Transition in Pulsars, A. Bransgrove, Y. Levin, A. M. Beloborodov, 2025
Using specialized numerical simulations, Ashley Bransgrove (Princeton), Yuri Levin, and Andrei Beloborodov (Columbia) studied how the transition to proton superconductivity in pulsar cores affects magnetic field evolution in the crust. The transition to type-II superconductivity in the outer core occurs soon after the pulsar is born and is associated with the quantization of the magnetic field into discrete superconducting flux tubes that possess enormous tension.
The simulations show how the enhanced tension of the flux tubes creates sharp magnetic gradients at the crust-core boundary, which launch giant magnetic Hall waves into the crust. The Maxwell stress associated with Hall waves may be strong enough to fracture the crust and could trigger starquakes, leading to rotation glitches and changes in the radio pulse profile. The Hall waves also couple to slow magnetospheric changes that cause measurable spin-down irregularities of the pulsar.
This work provides a robust mechanism to activate magnetic activity in pulsar crusts, suggesting a new direction to probe the exotic physics of neutron star interiors using astrophysical observables.