2D Relativistic MHD simulations of the Kruskal–Schwarzschild instability in a relativistic striped wind. Issue 3 (23rd November 2017)

• Record Type: Journal Article
• Authors: Gill, Ramandeep; Granot, Jonathan; Lyubarsky, Yuri
• Is Part Of: Monthly notices of the Royal Astronomical Society. Volume 474:Issue 3(2018)
• Page Start: 3535
• Page End: 3546
• Publication Date: 2017-11-23
• Subjects: instabilities -- magnetic reconnection -- MHD -- plasmas -- stars: jets; Astronomy -- Periodicals; Periodicals; 520.5
• DOI: 10.1093/mnras/stx3000
• Languages: English
• ISSNs: 0035-8711

Abstract

Abstract: We study the linear and non-linear development of the Kruskal–Schwarzchild instability in a relativisitically expanding striped wind. This instability is the generalization of Rayleigh–Taylor instability in the presence of a magnetic field. It has been suggested to produce a self-sustained acceleration mechanism in strongly magnetized outflows found in active galactic nuclei, gamma-ray bursts, and micro-quasars. The instability leads to magnetic reconnection, but in contrast with steady-state Sweet–Parker reconnection, the dissipation rate is not limited by the current layer's small aspect ratio. We performed two-dimensional (2D) relativistic magnetohydrodynamic (RMHD) simulations featuring two cold and highly magnetized (1 ≤ σ ≤ 10 3 ) plasma layers with an anti-parallel magnetic field separated by a thin layer of relativistically hot plasma with a local effective gravity induced by the outflow's acceleration. Our simulations show how the heavier relativistically hot plasma in the reconnecting layer drips out and allows oppositely oriented magnetic field lines to reconnect. The instability's growth rate in the linear regime matches the predictions of linear stability analysis. We find turbulence rather than an ordered bulk flow near the reconnection region, with turbulent velocities up to ∼0.1c, largely independent of model parameters. However, the magnetic energy dissipation rate is found to be much slower, corresponding to an effective ordered bulk velocity inflow into the reconnection region v in = βin c of 10 −3 ≲ βin ≲ 5 × 10 −3 . This occurs due to the slow evacuation of hot plasma from the current layer, largely because of the Kelvin–Helmholtz instability experienced by the dripping plasma. 3D RMHD simulations are needed to further investigate the non-linear regime.