skip to main content

Title: Spectral signatures of recursive magnetic field reconnection

We use 2.5D magnetohydrodynamic simulations to investigate the spectral signatures of the non-linear disruption of a tearing unstable current sheet via the generation of multiple secondary current sheets and magnetic islands. During the non-linear phase of tearing mode evolution, there develops a regime in which the magnetic energy density shows a spectrum with a power law close to B(k)2 ∼ k−0.8. Such an energy spectrum is found in correspondence of the neutral line, within the diffusion region of the primary current sheet, where energy is conveyed towards smaller scales via a ‘recursive’ process of fast tearing-type instabilities. Far from the neutral line, we find that magnetic energy spectra evolve towards slopes compatible with the ‘standard’ Kolmogorov spectrum. Starting from a self-similar description of the non-linear stage at the neutral line, we provide a model that predicts a reconnecting magnetic field energy spectrum scaling as k−4/5, in good agreement with numerical results. An extension of the predicted power law to generic current sheet profiles is also given and possible implications for turbulence phenomenology are discussed. These results provide a step forward to understand the ‘recursive’ generation of magnetic islands (plasmoids), which has been proposed as a possible explanation for the more » energy release during flares, but which, more in general, can have an impact on the subsequent turbulent evolution of unstable sheets that naturally form in the high Lundquist number and collisionless plasmas found in most of the astrophysical environments.

« less
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 4267-4276
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. We investigate particle acceleration in an MHD-scale system of multiple current sheets by performing 2D and 3D MHD simulations combined with a test particle simulation. The system is unstable for the tearing-mode instability, and magnetic islands are produced by magnetic reconnection. Due to the interaction of magnetic islands, the system relaxes to a turbulent state. The 2D (3D) case both yield −5/3 (− 11/3 and −7/3) power-law spectra for magnetic and velocity fluctuations. Particles are efficiently energized by the generated turbulence, and form a power-law tail with an index of −2.2 and −4.2 in the energy distribution function for the 2D and 3D case, respectively. We find more energetic particles outside magnetic islands than inside. We observe super-diffusion in the 2D (∼ t 2.27 ) and 3D (∼ t 1.2 ) case in the energy space of energetic particles.
  2. Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially two-dimensional) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in three dimensions, we conduct a comprehensive particle-in-cell simulation study comparing two- and three-dimensional evolution of long, thin current sheets in moderately magnetized, collisionless, relativistically hot electron–positron plasma, and find dramatic differences. We first systematically characterize this process in two dimensions, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths and system sizes. We then show that three-dimensional (3-D) simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3-D current sheet evolution is not always fundamentally classical reconnection with perturbing 3-D effects but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details and even stochastic events. It often yields slower but longer-lasting andmore »ultimately greater magnetic energy release than in two dimensions. Intriguingly, non-thermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in two and three dimensions. Although the variety of underlying current sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.« less
  3. Abstract Particle acceleration behind a shock wave due to interactions between magnetic islands in the heliosphere has attracted attention in recent years. The downstream acceleration may yield a continuous increase of particle flux downstream of the shock wave. Although it is not obvious how the downstream magnetic islands are produced, it has been suggested that current sheets are involved in the generation of magnetic islands due to their interaction with a shock wave. We perform 2D hybrid kinetic simulations to investigate the interaction between multiple current sheets and a shock wave. In the simulation, current sheets are compressed by the shock wave and a tearing instability develops at the compressed current sheets downstream of the shock. As the result of this instability, the electromagnetic fields become turbulent and magnetic islands form well downstream of the shock wave. We find a “post-cursor” region in which the downstream flow speed normal to the shock wave in the downstream rest frame is decelerated to ∼ 1 V A immediately behind the shock wave, where V A is the upstream Alfvén speed. The flow speed then gradually decelerates to 0 accompanied by the development of the tearing instability. We also observe an efficient productionmore »of energetic particles above 100 E 0 during the development of the instability some distance downstream of the shock wave, where E 0 = m p V A 2 and m p is the proton mass. This feature corresponds to Voyager observations showing that the anomalous cosmic-ray intensity increase begins some distance downstream of the heliospheric termination shock.« less
  4. The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde {t}^{-1}$ , the number of islands decreases as $\tilde {t}^{-1}$ and the averaged area of islands grows as $\tilde {t}$ , where $\tilde {t}$ is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decaymore »of magnetohydrodynamic turbulence.« less
  5. Abstract We investigate the interaction of turbulence with shock waves by performing 2D hybrid kinetic simulations. We inject force-free magnetic fields upstream that are unstable to the tearing-mode instability. The magnetic fields evolve into turbulence and interact with a shock wave whose sonic Mach number is 2.4. Turbulence properties, the total and normalized residual energy and the normalized cross helicity, change across the shock wave. While the energy of velocity and magnetic fluctuations is mostly distributed equally upstream, the velocity fluctuations are amplified dominantly downstream of the shock wave. The amplitude of turbulence spectra for magnetic, velocity, and density fluctuations are also increased at the shock wave while their spectral index remains unchanged. We compare our results with the Zank et al. model of turbulence transmission across a shock, and find that it provides a reasonable explanation for the spectral change across the shock wave. We find that particles are efficiently accelerated at the shock front, and a power-law spectrum forms downstream. This can be explained by diffusive shock acceleration, in which particles gain energy by being scattered upstream and downstream of a shock wave. The trajectory of an accelerated particle suggests that upstream turbulence plays a role scattering ofmore »particles.« less