skip to main content

Title: Transition from ion-coupled to electron-only reconnection: Basic physics and implications for plasma turbulence
Using 2.5 dimensional kinetic particle-in-cell (PIC) simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion inertial length, the ions do not respond to the reconnection dynamics leading to “electron-only” reconnection with very large quasi-steady reconnection rates. Note that in these simulations the ion Larmor radius is comparable to the ion inertial length. The transition to more traditional “ion-coupled” reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasi-steady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection more » Alfven speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths. « less
Award ID(s):
Publication Date:
Journal Name:
Physics of plasmas
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. Abstract One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfvénic magnetic field reversals termed switchbacks . These δ B R / B ∼  ( 1 ) fluctuations occur over a range of timescales and in patches separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate that patches of switchbacks are localized within the extensions of plasma structures originating at the base of the corona. These structures are characterized by an increase in alpha particle abundance, Mach number, plasma β and pressure, and by depletions in the magnetic field magnitude and electron temperature. These intervals are in pressure balance, implying stationary spatial structure, and the field depressions are consistent with overexpanded flux tubes. The structures are asymmetric in Carrington longitude with a steeper leading edge and a small (∼1°) edge of hotter plasma and enhanced magnetic field fluctuations. Some structures contain suprathermal ions to ∼85 keV that we argue are the energetic tail of the solar wind alpha population. The structures are separated in longitude by angular scales associated with supergranulation. This suggests that these switchbacks originate near the leadingmore »edge of the diverging magnetic field funnels associated with the network magnetic field—the primary wind sources. We propose an origin of the magnetic field switchbacks, hot plasma and suprathermals, alpha particles in interchange reconnection events just above the solar transition region and our measurements represent the extended regions of a turbulent outflow exhaust.« less
  3. Observations in Earth’s turbulent magnetosheath downstream of a quasiparallel bow shock reveal a prevalence of electron-scale current sheets favorable for electron-only reconnection where ions are not coupled to the reconnecting magnetic fields. In small-scale turbulence, magnetic structures associated with intense current sheets are limited in all dimensions. And since the coupling of ions are constrained by a minimum length scale, the dynamics of electron reconnection is likely to be 3D. Here, both 2D and 3D kinetic particle-in-cell simulations are used to investigate electron-only reconnection, focusing on the reconnection rate and associated electron flows. A new form of 3D electron-only reconnection spontaneously develops where the magnetic X-line is localized in the out-of-plane (z) direction. The consequence is an enhancement of the reconnection rate compared with two dimensions, which results from differential mass flux out of the diffusion region along z, enabling a faster inflow velocity and thus a larger reconnection rate. This outflow along z is due to the magnetic tension force in z just as the conventional exhaust tension force, allowing particles to leave the diffusion region efficiently along z unlike the 2D configuration.
  4. We study the time-dependent formation and evolution of a current sheet (CS) in a magnetised, collisionless, high-beta plasma using hybrid-kinetic particle-in-cell simulations. An initially tearing-stable Harris sheet is frozen into a persistently driven incompressible flow so that its characteristic thickness gradually decreases in time. As the CS thins, the strength of the reconnecting field increases, and adiabatic invariance in the inflowing fluid elements produces a field-biased pressure anisotropy with excess perpendicular pressure. At large plasma beta, this anisotropy excites the mirror instability, which deforms the reconnecting field on ion-Larmor scales and dramatically reduces the effective thickness of the CS. Tearing modes whose wavelengths are comparable to that of the mirrors then become unstable, triggering reconnection on smaller scales and at earlier times than would have occurred if the thinning CS were to have retained its Harris profile. A novel method for identifying and tracking X-points is introduced, yielding X-point separations that are initially intermediate between the perpendicular and parallel mirror wavelengths in the upstream plasma. These mirror-stimulated tearing modes ultimately grow and merge to produce island widths comparable to the CS thickness, an outcome we verify across a range of CS formation timescales and initial CS widths. Our results maymore »find their most immediate application in the tearing disruption of magnetic folds generated by turbulent dynamo in weakly collisional, high-beta, astrophysical plasmas.« less
  5. The Hall Magnetohydrodynamic (MHD) equations are an extension of the standard MHD equations that include the “Hall” term from the general Ohm’s law. The Hall term decouples ion and electron motion physically on the ion inertial length scales. Implementing the Hall MHD equations in a numerical solver allows more physical simulations for plasma dynamics on length scales less than the ion inertial scale length but greater than the electron inertial length. The present effort is an important step towards producing physically correct results to important problems, such as the Geospace Environmental Modeling (GEM) Magnetic Reconnection problem. The solver that is being modified is currently capable of solving the resistive MHD equations on unstructured grids using the spectral difference scheme which is an arbitrarily high-order method that is relatively simple to parallelize. The GEM Magnetic Reconnection problem is used to evaluate whether the Hall MHD equations have been correctly implemented in the solver using the spectral difference method with divergence cleaning (SDDC) algorithm by comparing against the reconnection rates reported in the literature.