# Electron-Only Reconnection in Plasma Turbulence

Francesco Califano, S. S. Cerri, Matteo Faganello, Dimitri Laveder, M. Sisti,

Matthew W. Kunz

► To cite this version:

Francesco Califano, S. S. Cerri, Matteo Faganello, Dimitri Laveder, M. Sisti, et al., Electron-Only Reconnection in Plasma Turbulence. Frontiers in Physics, 2020, 8, 10.3389/fphy.2020.00317. hal-03101917

HAL Id: hal-03101917

https://hal.science/hal-03101917

Submitted on 11 Jan 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Electron-Only Reconnection in Plasma Turbulence

Francesco Califano¹*, Silvio Sergio Cerri², Matteo Faganello³, Dimitri Laveder⁴,
Manuela Sisti¹,³ and Matthew W. Kunz²,⁵

¹ Dipartimento di Fisica 'E. Fermi', Università di Pisa, Pisa, Italy, ² Department of Astrophysical Sciences, Princeton University, Princeton, NJ, United States, ³ Aix-Marseille Université, CNRS, PIIM UMR 7345, Marseille, France, ⁴ Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Laboratoire J.L. Lagrange, Nice, France, ⁵ Princeton Plasma Physics Laboratory, Princeton, NJ, United States

Hybrid-Vlasov-Maxwell simulations of magnetized plasma turbulence including non-linear electron-inertia effects in a generalized Ohm's law are presented. When fluctuation energy is injected on scales sufficiently close to ion-kinetic scales, the ions efficiently become de-magnetized and electron-scale current sheets largely dominate the distribution of the emerging current structures, in contrast to the usual picture, where a full hierarchy of structure sizes is generally observed. These current sheets are shown to be the sites of electron-only reconnection (e-rec), in which the usual electron exhausts are unaccompanied by ion outflows and which are in qualitative agreement with those recently observed by MMS in the Earth's turbulent magnetosheath, downstream of the bow shock. Some features of the e-rec phenomenology are shown to be consistent with an electron magnetohydrodynamic description. Simulations suggest that this regime of collisionless reconnection may be found in turbulent systems where plasma processes, such as micro-instabilities and/or shocks, overpower the more customary turbulent cascade by directly injecting energy close to the ion-kinetic scales.

**Keywords:** magnetic fields, magnetic reconnection, plasma turbulence, solar wind, Earth magnetosheath, plasma simulations

## 1. INTRODUCTION

Highly accurate *in situ* satellite measurements of plasma fluctuations and particle distribution functions in the heliosphere are today the primary experimental tool for the study of plasma turbulence (e.g., [1-3]) and magnetic reconnection (e.g., [4, 5]). In particular, space missions such as Cluster and MMS provide unprecedented opportunities to investigate multi-scale plasma physics, from the ion-kinetic scales down to the electron-kinetic scales, and thereby to constrain theoretical models of kinetic turbulence and reconnection.

Recently, MMS has measured in the Earth's turbulent magnetosheath a series of electron-scale reconnection events in which super-Alfvénic electron jets are never accompanied by ion outflows [6]. Those authors dubbed these events "electron-only reconnection" (hereafter, e-rec), to differentiate them from the usual (fast) collisionless reconnection process (e.g., [7-10]) in which reconnection takes place in an "electron diffusion region" embedded within a larger "ion diffusion region" and in which both species are expelled in a collimated outflow named the "exhaust" (e.g., [11]). While e-rec has been recently shown to occur in kinetic simulations of magnetic reconnection when a single electron-scale CS is ad-hoc initialized [12] or after a quasi-parallel shock [13], it is difficult to explain e-rec in the context of a classical turbulent cascade, in which turbulent energy is

**OPEN ACCESS**

**Edited by:**
Luca Sorriso-Valvo,
Institute for Science and Technology
of Plasmas, National Research
Council, Italy

**Reviewed by:**
Stefano Markidis,
KTH Royal Institute of Technology,
Sweden
Keizo Fujimoto,
Beihang University, China

*Correspondence:
Francesco Califano
francesco.califano@unipi.it

**Specialty section:**
This article was submitted to
Space Physics,
a section of the journal
Frontiers in Physics

**Received:** 24 March 2020
**Accepted:** 10 July 2020
**Published:** 04 September 2020

**Citation:**
Califano F, Cerri SS, Faganello M, Laveder D, Sisti M and Kunz MW
(2020) Electron-Only Reconnection in
Plasma Turbulence.
Front. Phys. 8:317.
doi: 10.3389/fphy.2020.00317
transferred conservatively from magneto-fluid scales down to ion-kinetic scales, thereby coupling to the ion dynamics in the usual way [4, 14–24]. This difficulty also holds in scenarios in which the turbulent cascade of magnetic fluctuations is shown to be mediated by the reconnection of ion-scale current sheets (e.g., [20, 22, 25–27]). Indeed, to the best of our knowledge, simulations in which energy is injected at scales much larger than the ion-kinetic length scales ($d_i$, $\rho_i$ being of the same order) show reconnecting CSs with both electron and ion outflows, regardless of the nature of this injection—whether it be continuous forcing or an initial distribution of magnetic and/or velocity fluctuations. This is true also for satellite observations of heliospheric turbulence (now able to resolve the electron scales; e.g., [28]), except for those MMS observations cited above.

In this paper, we use hybrid-Vlasov-Maxwell (HVM) simulations of freely decaying, 2D-3V turbulence to show that e-rec events qualitatively similar to those observed by MMS can develop in an “MHD-scale” turbulent system if magnetic fluctuations are injected sufficiently close to the ion-kinetic scale. These simulations include the Hall and electron-inertia terms in a generalized Ohm’s law [29] that captures the decoupling of the magnetic field from the ion dynamics at ion-kinetic scales and allows the complete unfreezing of magnetic flux with respect to the electron fluid motion at the electron-inertial scale. We also demonstrate that simulations of plasmas under the same conditions but with turbulent energy injected farther away from the ion-kinetic scale do not show e-rec, and instead exhibit standard multi-scale reconnection with both ion and electron outflows. The transition from standard reconnection to e-rec is found to occur when the wavenumber range of the injected fluctuations gets close to the ion-kinetic scales, namely when the largest injected wavenumber is varied from ($k_\perp d_i$)$_{max}$ = 0.3 to ($k_\perp d_i$)$_{max}$ = 0.6. Finally, we show that the physics underlying some features of the e-rec phenomenology can be described by the equations of electron magnetohydrodynamics (EMHD; [30–32]).

# 2. METHOD OF SOLUTION

## 2.1. Basic Equations

We integrate the Vlasov equation for the ion distribution function $f_i(t, r, v)$ using an Eulerian approach [33], coupled with Faraday’s law of induction for the evolution of the magnetic field $B(t, r)$. In dimensionless units, these are, respectively,

$$ \frac{\partial f_i}{\partial t} + v \cdot \nabla f_i + (E + v \times B) \cdot \frac{\partial f_i}{\partial v} = 0, \quad (1) $$

$$ \frac{\partial B}{\partial t} = -\nabla \times E. \quad (2) $$

Within the hybrid-Vlasov-Maxwell (HVM) framework, these equations are coupled to an Ohm’s law for the electric field $E(t, r)$ that includes contributions from the pressure and inertia of an isothermal electron fluid [29]:

$$ (1-d_e^2 \nabla^2) E = -u_e \times B - T_e \nabla \ln n + d_e^2 \nabla \cdot [n(u_i u_i - u_e u_e)]. \quad (3) $$

All equations are normalized with respect to the ion quantities: the mass $m_i$, charge $e$, inertial length $d_i$, and cyclotron frequency $\Omega_i \doteq eB_0/m_i c$, where $B_0$ is the strength of the mean magnetic field. The electron fluid is characterized by a finite electron skin depth $d_e = \sqrt{m_e/m_i}$, a flow velocity $u_e = u_i - J/n$, and a constant electron temperature $T_e$. The number density $n$ and the ion momentum density $nu_i$ are computed as the zeroth and first velocity-space moments of $f_i$, respectively. Quasi-neutrality is assumed, $n_e \simeq n_i = n$, and the displacement current is neglected in Ampère’s law, so that the current density $J = \nabla \times B$. Electron inertia, see Equation (3), is a key ingredient in this work as it allows for collisionless reconnection to occur even in the framework of a fluid-electron model. At the same time, the Hall term (hidden in the $-u_e \times B$ term, equivalent to $-u_i \times B + J/n \times B$) is the key actor for speeding up the dynamics around the X-point and leading to “fast” magnetic reconnection (see [35] or Appendix B in [34] for further details). We mention that a fully kinetic approach would be more appropriate for the physical description of the diffusive effects since the agyrotropic electron pressure tensor plays a major role (see for instance [36, 37]). However, such an approach is out of the scope of this paper and left for future investigation.

## 2.2. Simulation Setup

We consider an initially uniform, Maxwellian, proton-electron plasma embedded in a homogeneous out-of-plane magnetic field, $B_0 = B_0 \hat{e}_z$. The plasma beta parameter for the ion species initially satisfies $\beta_i \doteq 8\pi n T_{i0}/B_0^2 = 1$, where $T_{i0}$ is the initial ion temperature. We set $T_e = T_{i0}$. A reduced mass ratio of $m_i/m_e = 144$ is employed, ensuring that $d_i$ and $d_e$ are well-separated. The HVM Equations (1)–(3) are then solved in a 2D-3V phase space with $1,024^2$ grid points spanning a real-space domain of size $L_x = L_y = 20\pi d_i$, corresponding to a uniform spatial resolution $\simeq 0.06 d_i \simeq 0.7 d_e$. In order to avoid spurious numerical effects at very small scales, we adopt high-order spectral filters [38] on the electromagnetic fields that act only on the high-k part of the spectrum for numerical stability, while reconnection is driven by electron inertia. The velocity space is sampled by $51^3$ uniformly distributed grid points spanning $[-5, 5]v_{th,i}$, where $v_{th,i} = \sqrt{\beta_i/2}$ is the initial ion thermal speed, corresponding to a resolution of $0.2 v_{th,i}$ in each velocity direction. We note that velocity fluctuations smaller than the grid resolution, $\Delta v = 0.2$ (in normalized units) in all directions, are well-recovered by the Eulerian approach thanks to the very low level of noise of the algorithm. Indeed, even when the mean flow is small (with respect to $\Delta v$) the asymmetries in the values that the distribution function assumes on the negative and positive part of the v-space grid are well captured by its first-order moment $\int v f(v) d^3 v$. We caution that this may not be true only if the width of the distribution is not well resolved, i.e., if $\Delta v / v_{th} > 1$.

We note that velocity fluctuations slower than the grid resolution, $dv = 0.2$ in all directions, are well recovered by the Eulerian approach where the full sampling of the distribution function allows to capture any mean-flow asymmetry except for the cold plasma limit.

Decaying turbulence is initialized by imposing random, isotropic magnetic-field perturbations, $\delta B = \delta B_{\parallel} \hat{e}_z + \delta B_{\perp}$
with $\nabla \cdot \delta\mathbf{B} = 0$. Such perturbations self-consistently excite "compressive" fluctuations in the velocity and density fields [23]. The choice of injecting compressive magnetic perturbations (rather than, for instance, purely Alfvénic fluctuations having $\delta B_\|$ = 0) is justified by the fact that inside the bow shock (i.e., in the magnetosheath) the majority of the fluctuations are observed to be compressive (see, e.g., [39, 40]). Furthermore, the injection scale in our simulations resides (purposefully) close to $d_i$, at which "Alfvénic" or "magneto-sonic" fluctuations are not clearly distinguishable. Accordingly, we let the (kinetic) response of the plasma to the magnetic-field perturbations to develop velocity and density fluctuations self-consistently. This particular choice has been indeed shown to not affect the development of the turbulent cascade at (and below) the ion scales [23]. However, a study about the precise dependence on a much wider variety of injection properties is beyond the scope of this work—which itself is meant to be a first proof-of-concept—and will be left for future investigation.

Two simulations, labeled sim.A and sim.B for convenience, have been performed. In sim.A, the wavenumbers $k$ of the injected magnetic-field perturbations lie in the range $0.1 \le (k_\perp d_i)_{\text{inj}} \le 0.6$ with a corresponding rms amplitude of the in-plane perturbations $|e_z \times \nabla \psi|^{\text{rms}} \simeq 0.15$. These values are chosen so that moderate-amplitude fluctuations are injected on scales close to the ion-kinetic scales, so as to separate the ion and the electron dynamics from the very beginning of the simulation (see the discussion in section 2.3). In sim.B we excite fewer modes whose wavenumbers lie farther away from the ion-kinetic scales, viz., $0.1 \le (k_\perp d_i)_{\text{inj}} \le 0.3$, corresponding to the regime studied in Franci et al. [20, 41], Cerri et al. [23, 42], and Cerri and Califano [22]. We adopt the same rms value of the magnetic fluctuations as in sim.A. A third simulation, sim.C, which is equivalent to sim.B in terms of injection scales (viz., $0.1 \le (k_\perp d_i)_{\text{inj}} \le 0.3$) but employs a larger rms value for the fluctuations (viz., $|e_z \times \nabla \psi|^{\text{rms}} \simeq 0.21$), has also been performed. It confirms that the qualitative results from sim.B do not depend strongly on the initial level of the fluctuations, although a more precise statement will require further studies that are left for future investigation. Therefore, since there is no more information in sim.C that is not already present in sim.B, in this paper we only present results from sim.A and sim.B.

## 2.3. The EMHD Limit at Sub-Ion Scales

As turbulent fluctuations cascade to smaller and smaller scales, their characteristic frequencies increase, eventually exceeding those described by the MHD model. At scales below the ion-kinetic scales but above the electron skin depth, the ions decouple from the magnetic field, which remains frozen in the electron fluid. In the limit where the spatial derivatives of the number density $n$ and of the ion bulk flow $u_i$ can be neglected with respect to the derivatives of the electron fluid velocity $u_e$, the plasma may be described by EMHD. In this case, Equations (2) and (3) reduce to Equation (7) of Bulanov et al. [31], in which the field $B' \doteq B - d_e^2 \nabla^2 B$ (rather than $B$) is frozen into an incompressible electron flow and the current is carried solely by the electron fluid, so that $u_e = -J/n$. Under these conditions, collisionless EMHD in two dimensions ensures the Lagrangian

conservation of $F \doteq \psi - d_e^2 \nabla^2 \psi$ [31, 43], where the flux function $\psi$ is related to the magnetic field via $\mathbf{B} = e_z \times \nabla \psi + B_z$. Note that while $B'$ and $F$ stay frozen in the electron fluid, $B$ and $\psi$ can reconnect. One advantage of the HVM model adopted here is that it includes EMHD as one of the limits of the generalized Ohm's law (3), while simultaneously accounting also for non-EMHD effects (e.g., compressibility) and the kinetic response of the ions. On this point, we also note that the electron inertia terms in Equation (3), aimed first at allowing reconnection without resistivity, should include a term proportional to $d_e^2 \nabla(\nabla \cdot \mathbf{E})$. Even if it is negligible in the EMHD limit, this term could influence the kinetic ion response since the electric field directly enters the ions' Vlasov equation. Our set of equations can be seen as the simplest model describing kinetic ions, fluid electrons, ion and electron decoupling at different scales without including dissipation, finally conserving the EMHD flux $F$ in the high-frequency limit. We nevertheless underline that in Equation (3) for the electric field we do not formally impose $\nabla \cdot \mathbf{E} = 0$. As a consequence the electrostatic feedback is anyway included in our hybrid model (see for instance the test on Landau damping on ion acoustic waves in [29]). Numerically, we have two points to solve the electron inertial length $d_e$ which is a somewhat borderline to completely separate the EMHD invariant from the flux function. This means that numerical filtering or intrinsic dissipation of the algorithm to advance the Vlasov equation (see [33]), could contribute to the reconnection process and partially dissipate $F$ and so $\psi$. We note (see section 3.1, in particular Figure 2) that $F$ and $\psi$ are quite well separated during their evolution suggesting that, even at a resolution that must menage at the same time the large-scale evolution and the small-scale reconnection events, electron inertia terms play an important role in allowing for reconnection to occur.

# 3. RESULTS

## 3.1. Sim.A: Signatures of e-rec in Kinetic Turbulence

We first describe the results from sim.A. Within roughly 1.5 eddy turnover times of the shortest initial wavelength ($\sim 30\Omega_i^{-1}$), CSs form and reconnection onsets. (Fully developed turbulence is realized only on the eddy turnover time of the outer-scale fluctuations, $\tau \sim 100\Omega_i^{-1}$.) In Figure 1, top left, we show a full-box view of the out-of-plane current density $J_z$ at $t = 42.5\Omega_i^{-1}$. For comparison, in the top right frame, we show the same quantity from sim.B (discussed in section 3.2). We observe the formation of several electron-scale CSs, e.g., CS1 (indicated by the rectangle box) at $(x, y) = (58, 11)$, CS2 at $(23, 11)$, and CS3 at $(21, 36)$. These CSs are characterized by widths of a few $d_e$. More interestingly, their lengths $L_{CS} \sim (2-5)d_i$ are shorter than those typically observed in other turbulence simulations (e.g., [19, 20, 22, 23, 44, 45]). On these scales, the electron fluid, and thus the magnetic field, is expected to decouple from the ions; it is only deep within these thin CSs that the magnetic field ultimately decouples from the electron fluid. These features are shown in Figure 1, bottom frame, which displays $J_z$ (black), the out-of-plane electric field in the frame of the electrons $E'_z \doteq$
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | (Left) From sim.A. (top) Shaded iso-contours of the out-of-plane current density $J_z$ at $t = 42.5\Omega_i^{-1}$. The rectangle indicates the position of CS1. (bottom) One-dimensional cut along $y$ at $x \simeq 58$ (marked by the vertical dashed line in the top panel) of $J_z$ (black line), $E'_z \doteq E_z + (u_e \times B)_z$ (red line), and $E_z + (u_i \times B)_z$ (blue line). (Right) From sim.B. (top) Shaded iso-contours of the out-of-plane current density $J_z$ at $t = 95\Omega_i^{-1}$. The square indicates the CS shown in **Figure 6** (the simulation plane has been slightly shifted to highlight the two magnetic islands forming this CS). (bottom) One-dimensional cut along $y$ at $x \simeq 18$ (marked by the vertical dashed line in the top panel) of $J_z$ (black line), $E'_z \doteq E_z + (u_e \times B)_z$ (red line), and $E_z + (u_i \times B)_z$ (blue line). (Note that the y-range of the bottom plot for $E'_z$ has just been shifted with respect to the corresponding plot for sim.A because the current peak in sim.B is negative. However, the length of the y-range shown, $|y_{max} - y_{min}|$, is the same to facilitate a meaningful comparison).

$E_z + (u_e \times B)_z$ (red), and the out-of-plane electric field in the frame of the ions $E_z + (u_i \times B)_z$ (blue), both vs. $y$ at $x \approx 58d_i$ (corresponding to the vertical dashed line in the accompanying top panel). While the electron fluid is seen to decouple from the magnetic field only deep within CS1 (at $y \approx 11d_i$, where $E'_z \not\approx 0$), the ions are approximately decoupled over much of the simulation domain (i.e., $E_z + (u_i \times B)_z \not\approx 0$ at all scales). The width of CS1 (either inferred by the condition $|z| > J_z^{rms}$ or by using the more accurate method discussed in section 3.3) is $\ell \approx 4d_e$, corresponding to the physical range on which e-rec has been observed to occur in the magnetosheath [6].

**Figure 2** presents a zoom-in of CS1 after the onset of e-rec and the formation of an active X-point structure at $t = 42.5\Omega_i^{-1}$. The top panel reveals a quadrupolar structure of $B_z$ near the CS, known to be caused by the Hall term allowing the electrons to locally decouple from the ion neutralizing background and to generate via their in-plane current the quadrupole [11, 32, 46, 47]. Quadrupole signature is today a kind of universal trademark of fast reconnection routinely used not only in simulations but also

in space observations to assert the presence of a reconnection event [48, 49]. Iso-contours of the flux function $\psi$ (dashed white lines) and of the corresponding EMHD invariant $F$ (black lines) coincide on scales larger than $d_e$, where electron inertia is unimportant. Both quantities are advected by the flow and well conserved everywhere except around the X-point, where $\psi$ locally breaks (and reconnects) on scales $\sim d_e$. This separation of the $\psi$ and $F$ iso-contours is a signature of the EMHD regime (see, e.g., [50]). Indeed, in the bottom panel of **Figure 2**, the x-component of the electron flow $u_{e,x}$ shows the typical electron jet structures coming out from the X-point (highlighted by $F$, the solid black lines). No evidence of corresponding ion outflows is found in the exhaust around the X-point.

In **Figure 3**, we show data taken from two virtual spacecraft passing through CS1 along the paths traced by the vertical dashed lines in **Figure 2**. These trajectories are chosen to be similar to those taken by MMS in Phan et al. [6]. In **Figures 3A, F**, we show the reversal of $B_x$ over a few $d_e$, highlighted by the vertical dashed lines corresponding to the CS1 boundaries
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | From sim.A. Zoom-in on CS1 at $t = 42.5\Omega_i^{-1}$. (Top) Shaded iso-contours of the out-of-plane magnetic field $B_z$ with superposed iso-contours of $\psi$ (dashed white) and $F$ (black). (Bottom) Shaded iso-contours of the x-component of the electron flow, $u_{\theta,x}$, and of $\psi$ (solid black lines). Dashed vertical lines trace the two virtual spacecraft trajectories shown in the two columns of **Figure 3**.</figcaption></figure>

(given by the condition $|J_z| > J_z^{\text{rms}}$). Within these boundaries, there are clear signatures of oppositely oriented super-Alfvénic electron jets, identified as the exhausts (Figures 3B,G), without any noticeable corresponding ion outflows (Figures 3C,H). This feature is in qualitative agreement with MMS measurements [6]. Note that ion jets do not appear even if we move the spacecraft trajectories further away from the X-point location, as would be the case in a “standard” reconnection event. The field-parallel electric field $E_\parallel$, $E'_z$, and the dissipation rate $J \cdot E'$ [51] all depart significantly from zero only across the CS (Figures 3D,E,I,J). Although $J \cdot E'$ oscillates, its integral across CS1 is positive, indicating the possibility of a region of net magnetic-to-particle energy conversion (our equations only allow conversion into electron bulk, rather than thermal, energy). All of these features seen in CS1 are observed also in the other CSs.

That the dynamics of these CSs is EMHD-like is supported by Figure 4, which shows spectra of the solenoidal $\mathbf{u}^{(\text{sol})}$ and irrotational $\mathbf{u}^{(\text{irr})}$ contribution to the in-plane ion and electron velocities and of the in-plane magnetic field during the onset of the first e-rec events (viz., $t = 42.5\Omega_i^{-1}$). These spectra demonstrate that the ion flow is almost incompressible in the range $d_i^{-1} < k_\perp < d_e^{-1}$, while the electron flow remains nearly incompressible across an even larger range. We have also verified that the solenoidal contribution to $\mathbf{u}_{e,\perp}$ around

the CSs largely dominates over its irrotational counterpart, and that the “EMHD terms” in our generalized Ohm’s law dominate the dynamics (viz., $|\nabla \cdot (nu_iu_{iz})| \ll |\nabla \cdot (nu_eu_{e,z})|$ and $|nu_{e,z}(\nabla \cdot \mathbf{u}_e)| \ll |\mathbf{u}_e \cdot \nabla (nu_{e,z})|$; see, e.g., [31]).

Crucially, the properties of the fluctuations shown in Figure 4 at $t \sim 42.5\Omega_i^{-1}$ do not change character once the fully developed turbulent stage is achieved at $t_p \sim 110\Omega_i^{-1}$ ($t_p$ corresponds to the time when the rms out-of-plane current $J_z^{(\text{rms})}$ reaches its peak; see, e.g., [44]). After nearly 3 eddy turnover times, corresponding to the fully turbulent regime and by which time the ions could have begun participating in the reconnection dynamics, no sign of ion outflows in sim.A is observed. Figure 5 (left panels) shows a comparison between the out-of-plane and in-plane electron and ion flow velocities at $t = 114\Omega_i^{-1}$ in sim.A; the iso-contours of $-u_{e,z} \approx J_z/n$ (top left), $u_{i,z}$ (top right), $u_{e,y}$ (bottom left), and $u_{i,y}$ (bottom right) are shown. All electron quantities exhibit many thin structures at the electron scale, which are well correlated with the CSs traced by $J_z$ (not shown here, as it is nearly identical to $-u_{e,z}$). On the other hand, all of the ion quantities exhibit much smoother variations and, in particular, are largely uncorrelated with the corresponding electron flows or CSs. As in Figure 2, electron jets are also visible (e.g., in $u_{e,y}$ at $(x,y) \approx (10,20)d_i$; note that here the jets are along $y$), while no corresponding ion outflows are apparent. Despite these e-rec processes, the fluctuation spectra in the ion-kinetic range of the fully developed turbulence are not significantly affected (when compared to sim.B). The magnetic and electric energy spectra continue to exhibit sub-ion-scale power laws close to $-2.8$ and $-1$, respectively, in the range $d_i^{-1} \lesssim k_\perp \lesssim d_e^{-1}$ (not shown here), similar to those found in previous gyro-kinetic, hybrid-kinetic, and fully kinetic simulations (e.g., [19, 41, 42, 52–61]) as well as satellite measurements (e.g., [2, 62–64]).

## 3.2. Sim.B: Standard Reconnection in Kinetic Turbulence

None of the e-rec characteristics highlighted in section 3.1 are observed in sim.B, in which the magnetic perturbations are injected farther away from $d_i$ (at least up to the time when the turbulence is fully developed, $t \approx 200\Omega_i^{-1}$). In particular, as soon as CSs form in sim.B (after 1.5 eddy-turnover times associated with the shortest initial wavelength, $\sim 60\Omega_i^{-1}$ here), we observe the development of “standard” reconnection with electron structures embedded in ion macro-layers. This is shown for sim.B in the top-right panel of Figure 1, which provides a full-box view of the out-of-plane current density $J_z$ at $t = 95\Omega_i^{-1}$. In particular, we see the formation of a CS located at $x = 18, y = 50$ (highlighted by the black square). As in sim.A, the CS width collapses down to the $d_e$ scale at the X-point. However, now its length along the in-plane magnetic field is much larger (more than $10d_i$) than those of the CSs observed in sim.A (left panels). As a consequence, at a sufficiently large distance from the X-point, ions couple to the magnetic-field/electron dynamics and ion outflows are generated. This is shown in the bottom-right panel (data vs. $y$ taken along the dashed line): by way of comparison with its sim.A counterpart (bottom-left panel), the ions are seen to be better coupled to the magnetic field away from
<figure><img src="image_5.png" /></figure>

the CS. This facilitates the ions' participation to the reconnection process; indeed, a zoom-in of this CS (see **Figure 6**) shows the X-point structure with both electron and ion outflows present. Even if smoother, the ion outflows point in opposite directions along the reconnecting in-plane magnetic field and are localized inside the separatrices.

To provide further evidence that the exhausts are composed of two coupled ion-electron jets, in **Figure 7** we show six successive plots of the electron and ion outflows measured to the right of the X-point of the CS shown in **Figure 6**. This sequence shows

a central, well-collimated, strong electron jet (continuous lines) that becomes less collimated farther away from the X-point. We also observe a smoother, less intense but clear ion jet superposed on the electron one (dashed lines).

Such an electron-ion structure is never observed in sim.A where the ions decouple from the magnetic-field and the electron dynamics over nearly the entire simulation domain (see **Figure 1**, bottom left frame, and **Figure 5**, left frame). This difference is observed not only when the first active CSs form, but also at later times when the turbulence is fully developed. We also note
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | From sim.A. Spectra at $t = 42.5\Omega_i^{-1}$ of the in-plane magnetic field (blue solid line), the solenoidal contributions to the in-plane ion and electron velocities $\mathbf{u}_{i,e}^{(\text{sol})}$ (orange and red solid lines, respectively), and the irrotational contributions of the same velocities $\mathbf{u}_{i,e}^{(\text{irr})}$ (green and black-dashed curves, respectively). Reference slopes of -1 and -3 slopes are provided.</figcaption></figure>

that, although the ion-electron coupling remains valid in the fully developed turbulence stage of sim.B, some electron-scale structures emerge, as shown in **Figure 5**, right frame, where finer structures develop inside an ion scale reconnecting structure. In summary, in sim.B the different CSs exhibit “standard” or even multi-scale reconnection with both well-developed ion and electron outflows. On the contrary, in sim.A only e-rec develops and “standard” reconnection is never observed, neither at the time when the first CSs activate nor when the turbulence is fully developed.

### 3.3. Statistical Analysis of Current Sheets' Properties

Quantitative differences between sim.A and sim.B can be further assessed by examining the statistical distribution of the CSs' characteristic widths and lengths. For this purpose we define a CS as a region where the magnitude of the current density is bigger than a given threshold $J_{thr}$. Following Zhdankin et al. [65], we define such a threshold as $J_{thr} = \sqrt{\langle J_z^2 \rangle + 3\sigma}$, where $\sigma = \sqrt{\langle J_z^4 \rangle - (\langle J_z^2 \rangle)^2}$. The widths and lengths of the CSs are then evaluated using two different techniques.

The first technique employs an automated procedure that calculates the width and length of each CS based on the shape of the local $|J_z|$. First, we define the (local) current peak as the maximum value of $|J_z|$. Then, given the eigenvector associated with the largest eigenvalue of the Hessian matrix of $|J_z|$ at the peak, we look at the interpolated profile along this direction and define the CS's width as the full width at half maximum of $|J_z|$. Performing a similar procedure to compute the length, namely, interpolating $|J_z|$ along the direction of the eigenvector associated with the smallest eigenvalue, could be misleading since the CSs are rarely “straight.” For this reason we infer the length as the maximal distance between two points belonging

the same CS (i.e., within the thresholded contour), following Zhdankin et al. [65].

We then evaluate the average CS width and length and their standard deviations by taking into account all of the current peaks present in a simulation between the time when the first CSs form ($t = 30\Omega_i^{-1}$ in sim.A and $t = 75\Omega_i^{-1}$ in sim.B) and the fully developed turbulent regime ($t = 128\Omega_i^{-1}$ in sim.A and $t = 203\Omega_i^{-1}$ in sim.B); in total, 7,182 (460) current peaks were identified in sim.A (sim.B). These quantities are shown in **Figure 8** for both sim.A (red markers) and sim.B (blue markers): the crosses indicate the average width/length of these runs, while the red (blue) shaded area indicates the parameter region within one standard deviation of these values for sim.A (sim.B). Note that the average width (a few $d_e$) and its standard deviation are very similar amongst sim.A and sim.B, suggesting that these characteristics are independent of the details of the energy injection. This is not particularly surprising, since the minimum width of a CS width is limited by the decoupling of the magnetic field from the electron dynamics. By contrast, the length distributions are noticeably different in sim.A and sim.B. The average CS length in sim.A is $\simeq 2.6d_i$ and the standard deviation is $\simeq 2.9d_i$, while in sim.B these values correspond to $\simeq 4.8d_i$ and $\simeq 6.2d_i$, respectively. Therefore, CSs in sim.B are not only longer on average, but they also show a more pronounced variability in length. The picture emerging from sim.B therefore conforms to the “standard” idea of plasma turbulence developing “a hierarchy of dissipative coherent structures extending down to electron scales” [66], whereas sim.A provides the first numerical evidence of the possibility to “shortcut” this hierarchy and develops only smaller-scale structures where e-rec occurs.

This first analysis procedure has the advantage of being automated, leading to a large statistical ensemble of CSs. That being said, this procedure does not distinguish between reconnecting CSs and non-reconnecting ones. For example, CSs across which the in-plane magnetic field (the only one that can reconnect in 2D) does not change sign but whose current densities are sufficiently large are counted by the automated procedure, despite their inability to reconnect. There are also CSs across which the in-plane magnetic field changes sign, but no clear signs of reconnection (e.g., electron and/or ion outflows emerging from an X point) are present. Our second technique is therefore to compute the widths and lengths of only those CSs that exhibit clear evidence of reconnection. At a fixed time, we visually inspect each CS (identified by our first technique) to assess whether the in-plane magnetic field changes sign and whether electron outflows are present. We then examine those CSs' time evolution to find evidence for the erosion of $\psi$ at the X point, indicating an ongoing topological modification of the magnetic field. The measured widths and lengths of these reconnecting CSs are marked by stars in **Figure 8**. We have verified, using the proxies $J_z$ and $E_z'$, that the estimated widths of these CSs represented by a star are correct.

Note that the widths of these reconnection sites are always smaller than the average CSs' width, in both sim.A and sim.B, corroborating the fact that active CSs collapse at the X point. By contrast, the lengths of reconnecting CSs exhibit quite different
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** | (Left) From sim.A. Shaded iso-contours of (minus) the out-of-plane electron flow, $-u_{e,z} \approx J_z/n$ (top left); the out-of-plane ion flow, $u_{i,z}$ (top right); and the y-component of the electron and ion flow, $u_{e,y}$ and $u_{i,y}$ (bottom left and bottom right, respectively) in a zoom-in region one-fourth the size of the simulation domain at $t = 114\Omega_i^{-1}$. (Right) From sim.B. Shaded iso-contours of (minus) the out-of-plane electron flow, $-u_{e,z} \approx J_z/n$ (top left); the out-of-plane ion flow, $u_{i,z}$ (top right); and the x-component of the electron and ion flow, $u_{e,x}$ and $u_{i,x}$ (bottom left and bottom right, respectively) in a zoom-in region one-fourth the size of the simulation domain at $t \approx 203\Omega_i^{-1}$.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**FIGURE 6** | From sim.B. In-plane velocity vector field for electrons (**left**) and ions (**right**) with superimposed the flux function $\psi$ iso-contours around the X-point at $x \approx 6d_i$, $y \approx 56d_i$, at $t = 105\Omega_i^{-1}$.</figcaption></figure>

behavior amongst these two simulations: in sim.A the majority of reconnecting CSs have a length in between the boundaries defined by the standard deviation, with only two cases (out of 8) of reconnecting CSs with a slightly larger length (although smaller than $9d_i$). On the other hand, all but one reconnecting CSs in sim.B are out of the boundaries defined by the standard deviation and have a length larger than $10d_i$. The only one with a

smaller length, at $\approx 2.9d_i$, is the CS found at $t = 203\Omega_i^{-1}$ where fine electron structures are visible (see **Figure 5**).

Even if the number of these reconnection sites is statistically smaller than the total number of CSs, we can identify a trend regarding CS width and length within the two simulations. In general, the ensemble of all the CSs in a simulation is not particularly representative of the reconnecting ones. In fact, all
<figure><img src="image_9.png" /><figcaption>**FIGURE 7** | From sim.B. Cuts vs. y of the rightward outflow from the X-point shown in **Figure 6**, at $t = 105\Omega_i^{-1}$. The x-component of the electron (solid) and ion (dash-dotted) velocities are plotted vs. y at different values of x starting close to the X-point and going as far away as $3d_i$.</figcaption></figure>

the reconnecting CSs have a width that is smaller than the average value of the corresponding distribution of CSs. This width is clustered around ~ (0.2–0.25)$d_i$ ≈ 3$d_e$ in both simulations, with a relatively small dispersion around that value (the width being always less than 0.35$d_i$). More interestingly, the lengths of the reconnecting CSs exhibit a completely different trend amongst the two simulations. In sim.A, the length of these reconnecting structures (which exhibit e-rec) is clustered around ~ 3.5$d_i$ ≈ 42$d_e$, and nearly all of them have a characteristic length which is within one standard deviation from the average value of the distribution of all the CSs in sim.A. In sim.B, on the other hand, nearly all of the reconnecting CSs (which exhibit ion outflows) have a length that is beyond one standard deviation from the average value of the distribution of all the CSs, i.e., larger than $10d_i = 120d_e$.

# 4. CONCLUSIONS AND DISCUSSION

Using numerical simulations of kinetic plasma turbulence in a hybrid Vlasov–Maxwell model that includes fully non-linear

electron-inertia effects (except for those related to the electron pressure tensor terms), we have demonstrated for the first time that short CSs in which electron-only reconnection takes place can naturally develop within kinetic plasma turbulence at $\beta \sim 1$ when fluctuations are injected on a range of wavenumbers near ion-kinetic scales (in our case, $k_\perp d_i = k_\perp \rho_i \le 0.6$). In this situation, all CSs showing clear evidence of ongoing magnetic reconnection exhibit electron exhausts that are unaccompanied by ion outflows, with the ions being decoupled from the magnetic-field and electron dynamics almost everywhere. The properties of these reconnecting CSs, including their statistical lengths and widths are in qualitative agreement with those “electron-scale” CSs recently measured in the turbulent magnetosheath by MMS [6], in which reconnection-driven electron outflows were unaccompanied by ion outflows. It is worth noticing that in this case where electron reconnection occurs associated with smaller CSs than for standard reconnection (sim.A), we do not observe the classical 2D island coalescence process [67], at least during the full duration of the simulation up to $t \simeq 140$.
<figure><img src="image_10.png" /><figcaption>**FIGURE 8** | A visual representation of the statistical distribution of width and length of all the CSs present in sim.A (red) and sim.B (blue): the cross indicates the average values and the shaded area the parameter region within one standard deviation. The widths and lengths of selected CSs that exhibit an X point and are clearly reconnecting are indicated by the stars.</figcaption></figure>

In our simplified model in which electrons are modeled as a fluid and decaying turbulence is initialized via compressive magnetic perturbations, this transition from “standard” to “electron-only” dynamics is observed by varying the maximum wave number of the initial fluctuations. In particular, if fluctuations are injected at wavelengths such that $k_{\perp}d_i \le 0.3$, a hierarchy of dissipative CSs emerge with a variety of lengths, whose evolution follows the “standard” reconnection dynamics in which both ion and electron exhausts are seen.

We have further shown that the wavelength range of energy injection has a measurable impact on the statistical properties of all the CSs: their average length and corresponding standard deviation are indeed different in sim.A and sim.B. In particular, sim.A shows a distribution of CSs' lengths that is centered around $\sim 2.6d_i$ with a standard deviation of their (asymmetric) distribution of $\sim 2.9d_i$, whereas the corresponding distribution in sim.B exhibits a larger average value ($\sim 4.8d_i$) and a much wider distribution, with a standard deviation of $\sim 6.2d_i$. This difference is even more apparent if we limit the analysis to only those CSs that can be clearly identified as reconnecting ones: in sim.A, these CSs are mostly clustered around a length of $\sim 3.5d_i$ and are never longer than $\approx 8-9d_i$. Moreover, their distribution is relatively thin and similar to the one obtained using all the CSs. By contrast, in sim.B we find both very long reconnecting CSs (up to $\sim 20d_i$) and short ones (a few $d_i$). These statistical properties directly relate to the physics of reconnection that is observed: short CSs in sim.A only produce e-rec, while long and short CSs in sim.B allow for the presence of ion outflows and

finer electron structures. Note that, differently to reconnection in a 1D initial Harris sheet [68, 69] where the length of the resulting CS is not limited in size and grows in time following the initial configuration, in the turbulent case CSs' characteristic length is limited. Indeed CSs develop in between magnetic flux ropes whose size sets the maximal length CSs can reach. Therefore, the energy injection scale plays an important role by setting the size of flux ropes, and so on the CSs' maximal length and on the possible kind of reconnection developing. Recent fully kinetic simulations of magnetic reconnection developing in an isolated current sheet [12] are in agreement with this picture.

Despite our initial conditions being oversimplified and not directly tied to any specific physical injection mechanism, we are encouraged by the close resemblance of sim.A's results to *MMS* data. That this resemblance holds not only during the onset of the first reconnection events in sim.B, but also during the fully developed turbulent state, is particularly convincing. We are then led to conjecture that energy injection occurring near ion-kinetic scales, perhaps due to velocity-space instabilities and/or shocks, can qualitatively alter the evolution of CSs and the dynamics of magnetic reconnection in plasma turbulence, potentially explaining the *MMS* results. In particular, if some energy-injection mechanism with properties similar to those employed in our simulations occurs past the bow shock, our results suggest that e-rec events could be detected relatively close to the bow shock and also further downstream in the magnetosheath. A more recent observational study using an extended *MMS* data set from the turbulent magnetosheath [70] seems to support this scenario. As the *MMS* mission goes on, the amount of collected data will enable a more detailed statistical analysis of these e-rec events. This kind of analysis should include also data collected outside of the magnetosheath, in the pristine solar wind, where different statistics for the reconnecting current sheets may be expected. However, due to the different level of fluctuations in the two regions, we stress that a high sensitivity of the spacecraft's instruments may be required in order to detect (the presumably smaller amount of) e-rec events in the pristine solar wind. Next-generation space missions would be most informative.

## DATA AVAILABILITY STATEMENT

The simulation dataset (UNIPI e-rec) is available at Cineca on the AIDA-DB. In order to access the meta-information and the link to the raw data, look at the tutorial at http://aida-space.eu/AIDAbd-iRODS.

## AUTHOR CONTRIBUTIONS

All authors have contributed to the data analysis, text editing, and to the discussions presented in the paper.

## FUNDING

This project has received funding from the European Union's Horizon 2020 research and innovation programme under Grant
Agreement No. 776262 (AIDA, www.aida-space.eu) SC and MK were supported by the National Aeronautics and Space Administration under Grant No. NNX16AK09G issued through the Heliophysics Supporting Research Program.

## ACKNOWLEDGMENTS

We acknowledge PRACE for awarding us access to the supercomputer Marconi, CINECA, Italy, where the calculations were performed under the grant no. 2017174107. FC and

SC thank Dr. C. Cavazzoni (CINECA, Italy) for his essential contribution to code parallelization and performance and Dr. M. Guarrasi (CINECA, Italy) for his contribution on code implementation on Marconi. This work was also granted access to the HPC resources of CINES/IDRIS under the allocations A0040407042 and A0060407042 made by GENCI. This paper originates from an unpublished arXiv letter at arXiv:1810.03957, but represents a significantly upgraded work. We acknowledge the Referee S. Markidis for his comments and discussion helping us to improve the content of the paper.

## REFERENCES

1. Alexandrova O, Saur J, Lacombe C, Mangeney A, Mitchell J, Schwartz SJ, et al. Universality of solar-wind turbulent spectrum from MHD to electron scales. *Phys Rev Lett.* (2009) **103**:165003. doi: 10.1103/PhysRevLett.103.165003

2. Sahraoui F, Goldstein ML, Robert P, Khotyaintsev YV. Evidence of a cascade and dissipation of solar-wind turbulence at the electron gyroscale. *Phys Rev Lett.* (2009) **102**:231102. doi: 10.1103/PhysRevLett.102.231102

3. Chen CHK, Boldyrev S. Nature of kinetic scale turbulence in the Earth's magnetosheath. *Astrophys J.* (2017) **842**:122. doi: 10.3847/1538-4357/aa74e0

4. Retinò A, Sundkvist D, Vaivads A, Mozer F, André M, Owen CJ. *In situ* evidence of magnetic reconnection in turbulent plasma. *Nat Phys.* (2007) **3**:236–8. doi: 10.1038/nphys574

5. Burch JL, Torbert RB, Phan TD, Chen LJ, Moore TE, Ergun RE, et al. Electron-scale measurements of magnetic reconnection in space. *Science*. (2016) **352**:aaf2939. doi: 10.1126/science.aaf2939

6. Phan TD, Eastwood JP, Shay MA, Drake JF, Sonnerup BUÖ, Fujimoto M, et al. Electron magnetic reconnection without ion coupling in Earth's turbulent magnetosheath. *Nature*. (2018) **557**:202–6. doi: 10.1038/s41586-018-0091-5

7. Coppi B. 'Inertial' instabilities in plasmas. *Phys Lett.* (1964) **11**:226. doi: 10.1016/0031-9163(64)90419-6

8. Drake JF, Kleva RG. Collisionless reconnection and the sawtooth crash. *Phys Rev Lett.* (1991) **66**:1458. doi: 10.1103/PhysRevLett.66.1458

9. Ottaviani M, Porcelli F. Nonlinear collisionless magnetic reconnection. *Phys Rev Lett.* (1993) **71**:3802. doi: 10.1103/PhysRevLett.71.3802

10. Drake JF, Biskamp D, Zeiler A. Breakup of the electron current layer during 3-D collisionless magnetic reconnection. *Geophys Res Lett.* (1997) **24**:2921. doi: 10.1029/97GL52961

11. Shay MA, Drake JF, Denton RE, Biskamp D. Structure of the dissipation region during collisionless magnetic reconnection. *J Geophys Res.* (1998) **103**:9165. doi: 10.1029/97JA03528

12. Sharma Pyakurel P, Shay MA, Phan TD, Matthaeus WH, Drake JF, TenBarge JM, et al. Transition from ion-coupled to electron-only reconnection: Basic physics and implications for plasma turbulence. *Phys Plasmas*. (2019) **26**:082307. doi: 10.1063/1.5090403

13. Bessho N, Chen LJ, Wang S, Hesse M, Wilson LB. Magnetic reconnection in a quasi-parallel shock: two-dimensional local particle-in-cell simulation. *Geophys Res Lett.* (2019) **46**:9352–61. doi: 10.1029/2019GL083397

14. Sundkvist D, Retinò A, Vaivads A, Bale SD. Dissipation in Turbulent Plasma due to Reconnection in Thin Current Sheets. *Phys Rev Lett.* (2007) **99**:025004. doi: 10.1103/PhysRevLett.99.025004

15. Servidio S, Valentini F, Califano F, Veltri P. Local kinetic effects in two-dimensional plasma turbulence. *Phys Rev Lett.* (2012) **108**:045001. doi: 10.1103/PhysRevLett.108.045001

16. Perrone D, Valentini F, Servidio S, Dalena S, Veltri P. Vlasov simulations of multi-ion plasma turbulence in the solar wind. *Astrophys J.* (2013) **762**:99. doi: 10.1088/0004-637X/762/2/99

17. Perrone D, Passot T, Laveder D, Valentini F, Sulem PL, Zouganelis I, et al. Fluid simulations of plasma turbulence at ion scales: comparison with Vlasov-Maxwell simulations. (2018) **25**:052302. doi: 10.1063/1.5026656

18. Karimabadi H, Roytershteyn V, Vu HX, Omelchenko YA, Scudder J, Daughton W, et al. The link between shocks, turbulence, and magnetic reconnection in collisionless plasmas. (2014) **21**:062308. doi: 10.1063/1.4882875

19. Franci L, Landi S, Matteini L, Verdini A, Hellinger P. Plasma beta dependence of the ion-scale spectral break of solar wind turbulence: high-resolution 2D hybrid simulations. *Astrophys J.* (2016) **833**:91. doi: 10.3847/1538-4357/833/1/91

20. Franci L, Cerri SS, Califano F, Landi S, Papini E, Verdini A, et al. Magnetic reconnection as a driver for a sub-ion-scale cascade in plasma turbulence. *Astrophys J Lett.* (2017) **850**:L16. doi: 10.3847/2041-8213/aa93fb

21. Wan M, Matthaeus WH, Roytershteyn V, Parashar TN, Wu P, Karimabadi H. Intermittency, coherent structures and dissipation in plasma turbulence. (2016) **23**:042307. doi: 10.1063/1.4945631

22. Cerri SS, Califano F. Reconnection and small-scale fields in 2D-3V hybrid-kinetic driven turbulence simulations. *New J Phys.* (2017) **19**:025007. doi: 10.1088/1367-2630/aa5c4a

23. Cerri SS, Franci L, Califano F, Landi S, Hellinger P. Plasma turbulence at ion scales: a comparison between particle in cell and Eulerian hybrid-kinetic approaches. *J Plasma Phys.* (2017) **83**:705830202. doi: 10.1017/S0022377817000265

24. Shay MA, Haggerty CC, Matthaeus WH, Parashar TN, Wan M, Wu P. Turbulent heating due to magnetic reconnection. (2018) **25**:012304. doi: 10.1063/1.4993423

25. Loureiro NL, Boldyrev S. Collisionless reconnection in magnetohydrodynamic and kinetic turbulence. *Astrophys J.* (2017) **850**:182. doi: 10.3847/1538-4357/aa9754

26. Mallet A, Schekochihin AA, Chandran BDG. Disruption of Alfvénic turbulence by magnetic reconnection in a collisionless plasma. *J Plasma Phys.* (2017) **83**:1–22. doi: 10.1017/S0022377817000812

27. Papini E, Franci L, Landi S, Verdini A, Matteini L, Hellinger P. Can hall magnetohydrodynamics explain plasma turbulence at sub-ion scales? *Astrophys J.* (2019) **870**:52. doi: 10.3847/1538-4357/aaf003

28. Vörös Z, Yordanova E, Varsani A, Genestreti KJ, Khotyaintsev YV, Li W, et al. MMS observation of magnetic reconnection in the turbulent magnetosheath. *J Geophys Res A.* (2017) **122**:11. doi: 10.1002/2017JA024535

29. Valentini F, Trávníček P, Califano F, Hellinger P, Mangeney A. A hybrid-Vlasov model based on the current advance method for the simulation of collisionless magnetized plasma. *J Comput Phys.* (2007) **225**:753. doi: 10.1016/j.jcp.2007.01.001

30. Kingsep AS, Chukbar KV, Ian'kov VV. Electron magnetohydrodynamics. *Voprosy Teorii Plazmy*. (1987) **16**:209.

31. Bulanov SV, Pegoraro F, Sakharov AS. Magnetic reconnection in electron magnetohydrodynamics. *Phys Fluids B*. (1992) **4**:2499. doi: 10.1063/1.860467

32. Mandt ME, Denton RE, Drake JF. Transition to whistler mediated magnetic reconnection. *Geophys Res Lett.* (1994) **21**:73. doi: 10.1029/93GL03382

33. Mangeney A, Califano F, Cavazzoni C, Travnicek P. A numerical scheme for the integration of the vlasov-maxwell system of equations. *J Comput Phys.* (2002) **179**:495. doi: 10.1006/jcp.2002.7071

34. Faganello M, Califano F, Pegoraro F. Being on time in magnetic reconnection. *New J Phys.* (2009) **11**:063008. doi: 10.1088/1367-2630/11/6/063008

35. Birn J, Drake JF, Shay MA, Rogers BN, Denton RE, Hesse M, et al. Geospace environmental modeling (gem) magnetic reconnection challenge. *J Geophys Res.* (2001) **106**:3715–9. doi: 10.1029/1999JA0900449

36. Cai HJ, Lee LC. The generalized Ohm's law in collisionless magnetic reconnection. *Phys Plasmas*. (1997) **4**:509–20. doi: 10.1063/1.872178
37. Hesse M, Schindler K, Birn J, Kuznetsova M. The diffusion region in collisionless magnetic reconnection. *Phys Plasmas*. **6**:1781–95. doi: 10.1063/1.873436

38. Lele SK. Compact finite difference schemes with spectral-like resolution. *J Comput Phys* (1992) **103**:16. doi: 10.1016/0021-9991(92)90324-R

39. Huang SY, Hadid LZ, Sahraoui F, Yuan ZG, Deng XH. On the existence of the Kolmogorov inertial range in the terrestrial magnetosheath turbulence. *Astrophys J Lett*. (2017) **836**:L10. doi: 10.3847/2041-8213/836/1/L10

40. Hadid LZ, Sahraoui F, Galtier S, Huang SY. Compressible Magnetohydrodynamic Turbulence in the Earth's Magnetosheath: Estimation of the Energy Cascade Rate Using in situ Spacecraft Data. *Phys Rev Lett*. (2018) **120**:055102. doi: 10.1103/PhysRevLett.120.055102

41. Franci L, Landi S, Matteini L, Verdini A, Hellinger P. High-resolution hybrid simulations of kinetic plasma turbulence at proton scales. *Astrophys J*. (2015) **812**:21. doi: 10.1088/0004-637X/812/1/

42. Cerri SS, Califano F, Jenko F, Told D, Rincon F. Subproton-scale cascades in solar wind turbulence: driven hybrid-kinetic simulations. *Astrophys J Lett*. (2016) **822**:L12. doi: 10.3847/2041-8205/822/1/L12

43. Morrison PJ, Greene JM. Noncanonical Hamiltonian density formulation of hydrodynamics and ideal magnetohydrodynamics. *Phys Rev Lett*. (1980) **45**:790. doi: 10.1103/PhysRevLett.45.790

44. Servidio S, Valentini F, Perrone D, Greco A, Califano F, Matthaeus WH, et al. A kinetic model of plasma turbulence. *J Plasma Phys*. (2015) **81**:325810107. doi: 10.1017/S0022377814000841

45. Valentini F, Servidio S, Perrone D, Califano F, Matthaeus WH, Veltri P. Hybrid Vlasov-Maxwell simulations of two-dimensional turbulence in plasmas. (2014) **21**:082307. doi: 10.1063/1.4893301

46. Biskamp D, Schwarz E, Drake JF. Two-fluid theory of collisionless magnetic reconnection. *Phys Plasmas*. (1997) **4**:1002–9. doi: 10.1063/1.872211

47. Uzdensky DA, Kulsrud RM. Physical origin of the quadrupole out-of-plane magnetic field in Hall-magnetohydrodynamic reconnection. (2006) **13**:062305. doi: 10.1063/1.2209627

48. Öïeriset M, Phan TD, Fujimoto M, Lin RP, Lepping RP. *In situ* detection of collisionless reconnection in the Earth's magnetotail. *Nature*. (2001) **412**:414–7. doi: 10.1038/35086520

49. Vaivads A, Khotyaintsev Y, André M, Retinó A, Buchert SC, Rogers BN, et al. Structure of the magnetic reconnection diffusion region from four-spacecraft observations. *Phys Rev Lett*. **93**:105001. doi: 10.1103/PhysRevLett.93.105001

50. Attico N, Califano F, Pegoraro F. Kinetic regimes of high frequency magnetic reconnection in a neutral sheet configuration. (2002) **9**:458. doi: 10.1063/1.1435002

51. Zenitani S, Hesse M, Klimas A, Kuznetsova M. New measure of the dissipation region in collisionless magnetic reconnection. *Phys Rev Lett*. (2011) **106**:195003. doi: 10.1103/PhysRevLett.106195003

52. Howes GG, Tenbarge JM, Dorland W, Quataert E, Schekochihin AA, Numata R, et al. Gyrokinetic simulations of solar wind turbulence from ion to electron scales. *Phys Rev Lett*. (2011) **107**:035004. doi: 10.1103/PhysRevLett.107.035004

53. Tenbarge JM, Podesta JJ, Klein KG, Howes GG. Interpreting magnetic variance anisotropy measurements in the solar wind. *Astrophys J*. (2012) **753**:107. doi: 10.1088/0004-637X/753/2/107

54. Told D, Jenko F, Tenbarge JM, Howes GG, Hammett GW. Multiscale nature of the dissipation range in gyrokinetic simulations of Alfvénic turbulence. *Phys Rev Lett*. (2015) **115**:025003. doi: 10.1103/PhysRevLett.115.025003

55. Franci L, Landi S, Verdini A, Matteini L, Hellinger P. Solar wind turbulent cascade from MHD to sub-ion scales: large-size 3D hybrid particle-in-cell simulations. *Astrophys J*. (2018) **853**:26. doi: 10.3847/1538-4357/aaa3e8

56. Cerri SS, Servidio S, Califano F. Kinetic cascade in solar-wind turbulence: 3D3V hybrid-kinetic simulations with electron inertia. *Astrophys J Lett*. (2017) **846**:L18. doi: 10.3847/2041-8213/aa87b0

57. Cerri SS, Kunz MW, Califano F. Dual phase-space cascades in 3D hybrid-vlasov-maxwell turbulence. *Astrophys J Lett*. (2018) **856**:L13. doi: 10.3847/2041-8213/aab557

58. Grošelj D, Cerri SS, Bañón Navarro A, Willmott C, Told D, Loureiro NF, et al. Fully kinetic versus reduced-kinetic modeling of collisionless plasma turbulence. *Astrophys J*. (2017) **847**:28. doi: 10.3847/1538-4357/aa894d

59. Grošelj D, Mallet A, Loureiro NF, Jenko F. Fully kinetic simulation of 3D kinetic Alfvén turbulence. *Phys Rev Lett*. (2018) **120**:105101. doi: 10.1103/PhysRevLett.120.105101

60. Arzamasskiy L, Kunz MW, Chandran BDG, Quataert E. Hybrid-kinetic simulations of ion heating in Alfvénic turbulence. *Astrophys J*. (2019) **879**:1. doi: 10.3847/1538-4357/ab20cc

61. Cerri SS, GroÅjelj D, Franci L. Kinetic plasma turbulence: recent insights and open questions from 3D3V simulations. *Front Astron Space Sci*. (2019) **6**:64. doi: 10.3389/fspas.2019.00064

62. Alexandrova O, Carbone V, Veltri P, Sorriso-Valvo L. Small-scale energy cascade of the solar wind turbulence. *Astrophys J*. (2008) **674**:1153–7. doi: 10.1086/524056

63. Chen CHK, Horbury TS, Schekochihin AA, Wicks RT, Alexandrova O, Mitchell J. Anisotropy of solar wind turbulence between ion and electron scales. *Phys Rev Lett*. (2010) **104**:255002. doi: 10.1103/PhysRevLett.104.255002

64. Lacombe C, Alexandrova O, Matteini L. Anisotropies of the magnetic field fluctuations at kinetic scales in the solar wind: cluster observations. *Astrophys J*. (2017) **848**:45. doi: 10.3847/1538-4357/aa8c06

65. Zhdankin V, Uzdensky DA, Perez JC, Boldyrev S. Statistical analysis of current sheets in three-dimensional magnetohydrodynamic turbulence. *Astrophys J*. (2013) **771**:124. doi: 10.1088/0004-637X/771/2/124

66. Karimabadi H, Roytershteyn V, Wan M, Matthaeus WH, Daughton W, Wu P, et al. Coherent structures, intermittent turbulence, and dissipation in high-temperature plasmas. (2013) **20**:012303. doi: 10.1063/1.4773205

67. Malara F, Veltri P, Carbone V. Competition among nonlinear effects in tearing instability saturation. *Phys Fluids B*. **4**:3070–86. doi: 10.1063/1.860477

68. Fujimoto K. Time evolution of the electron diffusion region and the reconnection rate in fully kinetic and large system. *Phys Plasmas*. (2006) **13**:072904. doi: 10.1063/1.2220534

69. Daughton W, Scudder J, Karimabadi H. Fully kinetic simulations of undriven magnetic reconnection with open boundary conditions. *Phys Plasmas*. (2006) **13**:072101. doi: 10.1063/1.2218817

70. Stawarcy JE, Eastwood JP, Phan TD, Gingell IL, Shay MA, Burch JL, et al. Properties of the turbulence associated with electron-only magnetic reconnection in earth's magnetosheath. *Astrophys J*. (2019) **877**:L37. doi: 10.3847/2041-8213/ab21c8

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Califano, Cerri, Faganello, Laveder, Sisti and Kunz. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.