- Award ID(s):
- 1654168
- NSF-PAR ID:
- 10328025
- 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
-
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
-
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
-
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 (
) is short and remains nearly unchanged because the electron gyrofrequency Ωe is independent of domain size. The quasi-steady stage of the EKHI (t MR) dominates the electron acceleration process and scales linearly with the size of the simulations asL /v A0, wherev A0is 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 ratiot MR/t G is sufficiently small, i.e., ≲10%. -
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 is , whereη rec≃ 0.06 is the inflow speed in units of the speed of light andω c=eB 0/mc is the gyrofrequency in the upstream magnetic field. They leave the region of active energization aftert esc, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measuret escin our simulations and find thatt esc∼t accforσ ≳ few. This leads to a universal (i.e.,σ -independent) power-law spectrum for the particles undergoing active acceleration, and for 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. -
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