Explore the vast range of titles available.

Turbulence plays a key role in the conversion of the energy of large-scale fields and flows to plasma heat, impacting the macroscopic evolution of the heliosphere and other astrophysical plasma systems. Although we have long been able to make direct spacecraft measurements of all aspects of the electromagnetic field and plasma fluctuations in near-Earth space, our understanding of the physical mechanisms responsible for the damping of the turbulent fluctuations in heliospheric plasmas remains incomplete. Here we propose an innovative field–particle correlation technique that can be used to measure directly the secular energy transfer from fields to particles associated with collisionless damping of the turbulent fluctuations. Furthermore, this novel procedure yields information about the collisionless energy transfer as a function of particle velocity, providing vital new information that can help to identify the dominant collisionless mechanism governing the damping of the turbulent fluctuations. Kinetic plasma theory is used to devise the appropriate correlation to diagnose Landau damping, and the field–particle correlation technique is thoroughly illustrated using the simplified case of the Landau damping of Langmuir waves in a 1D-1V (one dimension in physical space and one dimension in velocity space) Vlasov–Poisson plasma. Generalizations necessary to apply the field–particle correlation technique to diagnose the collisionless damping of turbulent fluctuations in the solar wind are discussed, highlighting several caveats. This novel field–particle correlation technique is intended to be used as a primary analysis tool for measurements from current, upcoming and proposed spacecraft missions that are focused on the kinetic microphysics of weakly collisional heliospheric plasmas, including the Magnetospheric Multiscale (MMS), Solar Probe Plus, Solar Orbiter and Turbulence Heating ObserveR (THOR) missions.

Turbulence is commonly observed in nearly collisionless heliospheric plasmas, including the solar wind and corona and the Earth’s magnetosphere. Understanding the collisionless mechanisms responsible for the energy transfer from the turbulent fluctuations to the particles is a frontier in kinetic turbulence research. Collisionless energy transfer from the turbulence to the particles can take place reversibly, resulting in non-thermal energy in the particle velocity distribution functions (VDFs) before eventual collisional thermalization is realized. Exploiting the information contained in the fluctuations in the VDFs is valuable. Here we apply a recently developed method based on VDFs, the field–particle correlation technique, to a $\\unicode[STIX]{x1D6FD}=1$ , solar-wind-like, low-frequency Alfvénic turbulence simulation with well-resolved phase space to identify the field–particle energy transfer in velocity space. The field–particle correlations reveal that the energy transfer, mediated by the parallel electric field, results in significant structuring of the VDF in the direction parallel to the magnetic field. Fourier modes representing the length scales between the ion and electron gyroradii show that energy transfer is resonant in nature, localized in velocity space to the Landau resonances for each Fourier mode. The energy transfer closely follows the Landau resonant velocities with varying perpendicular wavenumber $k_{\\bot }$ and plasma $\\unicode[STIX]{x1D6FD}$ . This resonant signature, consistent with Landau damping, is observed in all diagnosed Fourier modes that cover the dissipation range of the simulation.

Magnetic reconnection and plasma turbulence are ubiquitous processes important for laboratory, space and astrophysical plasmas. Reconnection has been suggested to play an important role in the energetics and dynamics of turbulence by observations, simulations and theory for two decades. The fundamental properties of reconnection at kinetic scales, essential to understanding the general problem of reconnection in magnetized turbulence, remain largely unknown at present. Here we present an application of the magnetic flux transport method that can accurately identify reconnection in turbulence to a three-dimensional simulation. Contrary to ideas that reconnection in turbulence would be patchy and unpredictable, highly extended reconnection X-lines, on the same order of magnitude as the system size, form at kinetic scales. Extended X-lines develop through bi-directional reconnection spreading. They satisfy critical balance characteristic of turbulence, which predicts the X-line extent at a given scale. These results present a picture of fundamentally extended reconnection in kinetic-scale turbulence.

Magnetic reconnection has been suggested to play an important role in the dynamics and energetics of plasma turbulence by spacecraft observations, simulations and theory over the past two decades, and recently, by magnetosheath observations of MMS. A new method based on magnetic flux transport (MFT) has been developed to identify reconnection activity in turbulent plasmas. This method is applied to a gyrokinetic simulation of two-dimensional (2D) plasma turbulence. Results on the identification of three active reconnection X-points are reported. The first two X-points have developed bi-directional electron outflow jets. Beyond the category of electron-only reconnection, the third X-point does not have bi-directional electron outflow jets because the flow is modified by turbulence. In all cases, this method successfully identifies active reconnection through clear inward and outward flux transport around the X-points. This transport pattern defines reconnection and produces a new quadrupolar structure in the divergence of MFT. This method is expected to be applicable to spacecraft missions such as MMS, Parker Solar Probe, and Solar Orbiter.

Magnetic reconnection plays an important role in converting energy while modifying field topology. This process takes place in varied plasma environments in which the transport of magnetic flux is intrinsic. Identifying active magnetic reconnection sites in in-situ observations is challenging. A new technique, Magnetic Flux Transport (MFT) analysis, has been developed recently and proven in numerical simulation for identifying active reconnection efficiently and accurately. In this study, we examine the MFT process in 37 previously reported electron diffusion region (EDR)/reconnection-line crossing events at the dayside magnetopause and in the magnetotail and turbulent magnetosheath using Magnetospheric Multiscale measurements. The coexisting inward and outward MFT flows at an X-point provides a signature that magnetic field lines become disconnected and reconnected. The application of MFT analysis to in-situ observations demonstrates that MFT can successfully identify active reconnection sites under complex varied conditions, including asymmetric and turbulent upstream conditions. It also provides a higher rate of identification than plasma outflow jets alone. MFT can be applied to in situ measurements from both single- and multi-spacecraft missions and laboratory experiments.

The spreading of the X-line out of the reconnection plane under a strong guide field is investigated using large-scale three-dimensional (3D) particle-in-cell (PIC) simulations in asymmetric magnetic reconnection. A simulation with a thick, ion-scale equilibrium current sheet (CS) reveals that the X-line spreads at the ambient ion/electron drift speeds, significantly slower than the Alfvén speed based on the guide field \\(V_{Ag}\\). Additional simulations with a thinner, sub-ion-scale CS show that the X-line spreads at \\(V_{Ag}\\) (Alfvénic spreading), much higher than the ambient species drifts. An Alfvénic signal consistent with kinetic Alfvén waves develops and propagates, leading to CS thinning and extending, which then results in reconnection onset. The continuous onset of reconnection in the signal propagation direction manifests as Alfvénic X-line spreading. The strong dependence on the CS thickness of the spreading speeds, and the X-line orientation are consistent with the collisionless tearing instability. Our simulations indicate that when the collisionless tearing growth is sufficiently strong in a thinner CS such that \\(\\gamma/\\Omega_{ci}\\gtrsim\\mathcal{O}(1)\\), Alfvénic X-line spreading can take place. Our results compare favorably with a number of numerical simulations and recent magnetopause observations. A key implications is that the magnetopause CS is typically too thick for Alfvénic X-line spreading to effectively take place.

Turbulence is ubiquitously observed in nearly collisionless heliospheric plasmas, including the solar wind and corona and the Earth's magnetosphere. Understanding the collisionless mechanisms responsible for the energy transfer from the turbulent fluctuations to the particles is a frontier in kinetic turbulence research. Collisionless energy transfer from the turbulence to the particles can take place reversibly, resulting in non-thermal energy in the particle velocity distribution functions (VDFs) before eventual collisional thermalization is realized. Exploiting the information contained in the fluctuations in the VDFs is valuable. Here we apply a recently developed method based on VDFs, the field-particle correlation technique, to a \\(\\beta=1\\), solar-wind-like, low-frequency Alfvénic turbulence simulation with well resolved phase space to identify the field-particle energy transfer in velocity space. The field-particle correlations reveal that the energy transfer, mediated by the parallel electric field, results in significant structuring of the ion and electron VDFs in the direction parallel to the magnetic field. Fourier modes representing the length scales between the ion and electron gyroradii show that energy transfer is resonant in nature, localized in velocity space to the Landau resonances for each Fourier mode. The energy transfer closely follows the Landau resonant velocities with varying perpendicular wavenumber \\(k_\\perp\\) and plasma \\(\\beta\\). This resonant signature, consistent with Landau damping, is observed in all diagnosed Fourier modes that cover the dissipation range of the simulation.

The orientation and stability of the reconnection x-line in asymmetric geometry is studied using three-dimensional (3D) particle-in-cell simulations. We initiate reconnection at the center of a large simulation domain to minimize the boundary effect. The resulting x-line has sufficient freedom to develop along an optimal orientation, and it remains laminar. Companion 2D simulations indicate that this x-line orientation maximizes the reconnection rate. The divergence of the non-gyrotropic pressure tensor breaks the frozen-in condition, consistent with its 2D counterpart. We then design 3D simulations with one dimension being short to fix the x-line orientation, but long enough to allow the growth of the fastest growing oblique tearing modes. This numerical experiment suggests that reconnection tends to radiate secondary oblique tearing modes if it is externally (globally) forced to proceed along an orientation not favored by the local physics. The development of oblique structure easily leads to turbulence inside small periodic systems.

Using 3D particle-in-cell (PIC) simulations, we study magnetic reconnection with the x-line being spatially confined in the current direction. We include thick current layers to prevent reconnection at two ends of a thin current sheet that has a thickness on an ion inertial (di) scale. The reconnection rate and outflow speed drop significantly when the extent of the thin current sheet in the current direction is < O(10 di). When the thin current sheet extent is long enough, we find it consists of two distinct regions; an inactive region (on the ion-drifting side) exists adjacent to the active region where reconnection proceeds normally as in a 2D case. The extent of this inactive region is ~ O(10 di), and it suppresses reconnection when the thin current sheet extent is comparable or shorter. The time-scale of current sheet thinning toward fast reconnection can be translated into the spatial-scale of this inactive region; because electron drifts inside the ion diffusion region transport the reconnected magnetic flux, that drives outflows and furthers the current sheet thinning, away from this region. This is a consequence of the Hall effect in 3D. While this inactive region may explain the shortest possible azimuthal extent of dipolarizing flux bundles at Earth, it may also explain the dawn-dusk asymmetry observed at the magnetotail of Mercury, that has a global dawn-dusk extent much shorter than that of Earth.

Plasma turbulence is ubiquitous in space and astrophysical plasmas, playing an important role in plasma energization, but the physical mechanisms leading to dissipation of the turbulent energy remain to be definitively identified. Kinetic simulations in two dimensions (2D) have been extensively used to study the dissipation process. How the limitation to 2D affects energy dissipation remains unclear. This work provides a model of comparison between two- and three-dimensional (3D) plasma turbulence using gyrokinetic simulations; it also explores the dynamics of distribution functions during the dissipation process. It is found that both 2D and 3D nonlinear gyrokinetic simulations of a low-beta plasma generate electron velocity-space structures with the same characteristics as that of linear Landau damping of Alfvén waves in a 3D linear simulation. The continual occurrence of the velocity-space structures throughout the turbulence simulations suggests that the action of Landau damping may be responsible for the turbulent energy transfer to electrons in both 2D and 3D, and makes possible the subsequent irreversible heating of the plasma through collisional smoothing of the velocity-space fluctuations. Although, in the 2D case where variation along the equilibrium magnetic field is absent, it may be expected that Landau damping is not possible, a common trigonometric factor appears in the 2D resonant denominator, leaving the resonance condition unchanged from the 3D case. The evolution of the 2D and 3D cases is qualitatively similar. However, quantitatively the nonlinear energy cascade and subsequent dissipation is significantly slower in the 2D case.