Generation of Large-scale Magnetic Fields through Kelvin-Helmholtz Instability

Magnetic field generation via large-scale unstable shear flows. (a) The large-scale component of the magnetic field quasi-cyclically reverses, analogous to how the Sun’s magnetic field switches polarity every 11 years. (b) The magnetic fields are generated via the non-traditional Upsilon-effect—a mean-field dynamo process postulated in 1990 and validated here for the first time. (c) The Upsilon-effect is driven by nonlinear jets induced by large-scale unstable flow.

Magnetic field generation via large-scale unstable shear flows. (a) The large-scale component of the magnetic field quasi-cyclically reverses, analogous to how the Sun’s magnetic field switches polarity every 11 years. (b) The magnetic fields are generated via the non-traditional Upsilon-effect—a mean-field dynamo process postulated in 1990 and validated here for the first time. (c) The Upsilon-effect is driven by nonlinear jets induced by large-scale unstable flow.
Proton Synchrotron Emission in M87 Flares

Structure of the reconnection layer in the three-dimensional PIC simulation. The blue rendering shows the proton density distribution, while the red rendering corresponds to pairs. Protons are accelerated by bouncing between converging, reconnection-driven inflows, leading to a “puffy” distribution around the layer, while pairs are strongly cooled and remain localized near the center of the layer. Colored lines represent magnetic field lines.

Structure of the reconnection layer in the three-dimensional PIC simulation. The blue rendering shows the proton density distribution, while the red rendering corresponds to pairs. Protons are accelerated by bouncing between converging, reconnection-driven inflows, leading to a “puffy” distribution around the layer, while pairs are strongly cooled and remain localized near the center of the layer. Colored lines represent magnetic field lines.
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.
Kinetically Informed Reconnection in Black Hole Magnetospheres

Left: The reconnection electric field, indicating a reconnection rate of 0.1 for the new resistivity model (top half), and instead 0.03 for uniform resistivity (bottom half). Right: A zoom into the resistive electric field in the X-point where reconnection occurs (top for the new resistivity model, bottom for uniform) as indicated by the silver rectangle on the left. White lines show magnetic field lines.

Left: The reconnection electric field, indicating a reconnection rate of 0.1 for the new resistivity model (top half), and instead 0.03 for uniform resistivity (bottom half). Right: A zoom into the resistive electric field in the X-point where reconnection occurs (top for the new resistivity model, bottom for uniform) as indicated by the silver rectangle on the left. White lines show magnetic field lines.
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.
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.
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.
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.
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.