skip to main content


Title: Statistical description of coalescing magnetic islands via magnetic reconnection
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) decay of magnetohydrodynamic turbulence.  more » « less
Award ID(s):
1654168
NSF-PAR ID:
10328025
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Plasma Physics
Volume:
87
Issue:
6
ISSN:
0022-3778
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract It has been recently shown numerically that there exists an inverse transfer of magnetic energy in decaying, nonhelical, magnetically dominated, magnetohydrodynamic turbulence in 3-dimensions (3D). We suggest that magnetic reconnection is the underlying physical mechanism responsible for this inverse transfer. In the two-dimensional (2D) case, the inverse transfer is easily inferred to be due to smaller magnetic islands merging to form larger ones via reconnection. We find that the scaling behaviour is similar between the 2D and the 3D cases, i.e., the magnetic energy evolves as t−1, and the magnetic power spectrum follows a slope of k−2. We show that on normalizing time by the magnetic reconnection timescale, the evolution curves of the magnetic field in systems with different Lundquist numbers collapse onto one another. Furthermore, transfer function plots show signatures of magnetic reconnection driving the inverse transfer. We also discuss the conserved quantities in the system and show that the behaviour of these quantities is similar between the 2D and 3D simulations, thus making the case that the dynamics in 3D could be approximately explained by what we understand in 2D. Lastly, we also conduct simulations where the magnetic field is subdominant to the flow. Here, too, we find an inverse transfer of magnetic energy in 3D. In these simulations, the magnetic energy evolves as t−1.4 and, interestingly, a dynamo effect is observed. 
    more » « less
  2. null (Ed.)
    ABSTRACT The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell vpic, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov–Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed. 
    more » « less
  3. Abstract

    To investigate how magnetic reconnection (MR) accelerates electrons to a power-law energy spectrum in solar flares, we explore the scaling of a kinetic model proposed by Che & Zank (CZ) and compare it to observations. Focusing on thin current sheet MR particle-in-cell (PIC) simulations, we analyze the impact of domain size on the evolution of the electron Kelvin–Helmholtz instability (EKHI). We find that the duration of the growth stage of the EKHI (tGΩe1) is short and remains nearly unchanged because the electron gyrofrequency Ωeis independent of domain size. The quasi-steady stage of the EKHI (tMR) dominates the electron acceleration process and scales linearly with the size of the simulations asL/vA0, wherevA0is the Alfvén speed. We use the analytical results obtained by CZ to calculate the continuous temporal evolution of the electron energy spectra from PIC simulations and linearly scale them to solar flare observational scales. For the first time, an electron acceleration model predicts the sharp two-stage transition observed in typical soft–hard–harder electron energy spectra, implying that the electron acceleration model must be efficient with an acceleration timescale that is a small fraction of the duration of solar flares. Our results suggest that we can use PIC MR simulations to investigate the observational electron energy spectral evolution of solar flares if the ratiotMR/tGis sufficiently small, i.e., ≲10%.

     
    more » « less
  4. Abstract

    Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model—supported and benchmarked with large-scale three-dimensional kinetic particle-in-cell simulations in electron–positron plasmas—that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factorγ≳ 3σ(here,σis the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time istaccγηrec1ωc120γωc1, whereηrec≃ 0.06 is the inflow speed in units of the speed of light andωc=eB0/mcis the gyrofrequency in the upstream magnetic field. They leave the region of active energization aftertesc, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measuretescin our simulations and find thattesctaccforσ≳ few. This leads to a universal (i.e.,σ-independent) power-law spectrumdNfree/dγγ1for the particles undergoing active acceleration, anddN/dγγ2for the overall particle population. Our results help to shed light on the ubiquitous presence of power-law particle and photon spectra in astrophysical nonthermal sources.

     
    more » « less
  5. 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. 
    more » « less