skip to main content

Title: Reconnection and particle acceleration in three-dimensional current sheet evolution in moderately magnetized astrophysical pair plasma
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 and more » 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
Authors:
;
Award ID(s):
1903335 1806084
Publication Date:
NSF-PAR ID:
10329489
Journal Name:
Journal of Plasma Physics
Volume:
87
Issue:
6
ISSN:
0022-3778
Sponsoring Org:
National Science Foundation
More Like this
  1. Over three decades of in-situ observations illustrate that the Kelvin–Helmholtz (KH) instability driven by the sheared flow between the magnetosheath and magnetospheric plasma often occurs on the magnetopause of Earth and other planets under various interplanetary magnetic field (IMF) conditions. It has been well demonstrated that the KH instability plays an important role for energy, momentum, and mass transport during the solar-wind-magnetosphere coupling process. Particularly, the KH instability is an important mechanism to trigger secondary small scale (i.e., often kinetic-scale) physical processes, such as magnetic reconnection, kinetic Alfvén waves, ion-acoustic waves, and turbulence, providing the bridge for the coupling of cross scale physical processes. From the simulation perspective, to fully investigate the role of the KH instability on the cross-scale process requires a numerical modeling that can describe the physical scales from a few Earth radii to a few ion (even electron) inertial lengths in three dimensions, which is often computationally expensive. Thus, different simulation methods are required to explore physical processes on different length scales, and cross validate the physical processes which occur on the overlapping length scales. Test particle simulation provides such a bridge to connect the MHD scale to the kinetic scale. This study applies different testmore »particle approaches and cross validates the different results against one another to investigate the behavior of different ion species (i.e., H+ and O+), which include particle distributions, mixing and heating. It shows that the ion transport rate is about 10 25  particles/s, and mixing diffusion coefficient is about 10 10  m 2  s −1 regardless of the ion species. Magnetic field lines change their topology via the magnetic reconnection process driven by the three-dimensional KH instability, connecting two flux tubes with different temperature, which eventually causes anisotropic temperature in the newly reconnected flux.« less
  2. In a magnetized, collisionless plasma, the magnetic moment of the constituent particles is an adiabatic invariant. An increase in the magnetic-field strength in such a plasma thus leads to an increase in the thermal pressure perpendicular to the field lines. Above a $\unicode[STIX]{x1D6FD}$ -dependent threshold (where $\unicode[STIX]{x1D6FD}$ is the ratio of thermal to magnetic pressure), this pressure anisotropy drives the mirror instability, producing strong distortions in the field lines on ion-Larmor scales. The impact of this instability on magnetic reconnection is investigated using a simple analytical model for the formation of a current sheet (CS) and the associated production of pressure anisotropy. The difficulty in maintaining an isotropic, Maxwellian particle distribution during the formation and subsequent thinning of a CS in a collisionless plasma, coupled with the low threshold for the mirror instability in a high- $\unicode[STIX]{x1D6FD}$ plasma, imply that the geometry of reconnecting magnetic fields can differ radically from the standard Harris-sheet profile often used in simulations of collisionless reconnection. As a result, depending on the rate of CS formation and the initial CS thickness, tearing modes whose growth rates and wavenumbers are boosted by this difference may disrupt the mirror-infested CS before standard tearing modes can develop. Amore »quantitative theory is developed to illustrate this process, which may find application in the tearing-mediated disruption of kinetic magnetorotational ‘channel’ modes.« less
  3. Context. Magnetic reconnection plays a fundamental role in plasma dynamics under many different conditions, from space and astrophysical environments to laboratory devices. High-resolution in situ measurements from space missions allow naturally occurring reconnection processes to be studied in great detail. Alongside direct measurements, numerical simulations play a key role in the investigation of the fundamental physics underlying magnetic reconnection, also providing a testing ground for current models and theory. The choice of an adequate plasma model to be employed in numerical simulations, while also compromising with computational cost, is crucial for efficiently addressing the problem under study. Aims. We consider a new plasma model that includes a refined electron response within the “hybrid-kinetic framework” (fully kinetic protons and fluid electrons). The extent to which this new model can reproduce a full-kinetic description of 2D reconnection, with particular focus on its robustness during the nonlinear stage, is evaluated. Methods. We perform 2D simulations of magnetic reconnection with moderate guide field by means of three different plasma models: (i) a hybrid-Vlasov-Maxwell model with isotropic, isothermal electrons, (ii) a hybrid-Vlasov-Landau-fluid (HVLF) model where an anisotropic electron fluid is equipped with a Landau-fluid closure, and (iii) a full-kinetic model. Results. When compared to themore »full-kinetic case, the HVLF model effectively reproduces the main features of magnetic reconnection, as well as several aspects of the associated electron microphysics and its feedback onto proton dynamics. This includes the global evolution of magnetic reconnection and the local physics occurring within the so-called electron-diffusion region, as well as the evolution of species’ pressure anisotropy. In particular, anisotropy-driven instabilities (such as fire-hose, mirror, and cyclotron instabilities) play a relevant role in regulating electrons’ anisotropy during the nonlinear stage of magnetic reconnection. As expected, the HVLF model captures all these features, except for the electron-cyclotron instability.« less
  4. Context. Accreting black holes (BHs) may be surrounded by a highly magnetized plasma threaded by an organized poloidal magnetic field. Nonthermal flares and power-law spectral components at high energy could originate from a hot, collisionless, and nearly force-free corona. The jets we often observe from these systems are believed to be rotation-powered and magnetically driven. Aims. We study axisymmetric BH magnetospheres, where a fraction of the magnetic field lines anchored in a surrounding disk are connected to the event horizon of a rotating BH. For different BH spins, we identify the conditions and sites of magnetic reconnection within 30 gravitational radii. Methods. With the fully general relativistic particle-in-cell code GRZeltron , we solve the time-dependent dynamics of the electron–positron pair plasma and of the electromagnetic fields around the BH. The aligned disk is represented by a steady and perfectly conducting plasma in Keplerian rotation, threaded by a dipolar magnetic field. Results. For prograde disks around Kerr BHs, the topology of the magnetosphere is hybrid. Twisted open magnetic field lines crossing the horizon power a Blandford-Znajek jet, while open field lines with their footpoint beyond a critical distance on the disk could launch a magneto-centrifugal wind. In the innermost regions, couplingmore »magnetic field lines ensure the transfer of significant amounts of angular momentum and energy between the BH and the disk. From the Y point at the intersection of these three regions, a current sheet forms where vivid particle acceleration via magnetic reconnection takes place. We compute the synchrotron images of the current sheet emission. Conclusions. Our estimates for jet power and BH–disk exchanges match those derived from purely force-free models. Particles are accelerated at the Y point, which acts as a heat source for the so-called corona. It provides a physically motivated ring-shaped source of hard X-rays above the disk for reflection models. Episodic plasmoid ejection might explain millisecond flares observed in Cygnus X-1 in the high-soft state, but are too fast to account for daily nonthermal flares from Sgr A * . Particles flowing from the Y point down to the disk could produce a hot spot at the footpoint of the outermost closed magnetic field line.« less
  5. Magnetic reconnection is a fundamental process providing topological changes of the magnetic field, reconfiguration of space plasmas and release of energy in key space weather phenomena, solar flares, coronal mass ejections and magnetospheric substorms. Its multiscale nature is difficult to study in observations because of their sparsity. Here we show how the lazy learning method, known as K nearest neighbors, helps mine data in historical space magnetometer records to provide empirical reconstructions of reconnection in the Earth’s magnetotail where the energy of solar wind-magnetosphere interaction is stored and released during substorms. Data mining reveals two reconnection regions (X-lines) with different properties. In the mid tail ( ∼ 30 R E from Earth, where R E is the Earth’s radius) reconnection is steady, whereas closer to Earth ( ∼ 20 R E ) it is transient. It is found that a similar combination of the steady and transient reconnection processes can be reproduced in kinetic particle-in-cell simulations of the magnetotail current sheet.