We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of largescale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbulence is fully developed, we observe a MHD inertial range, where the spectra of the perpendicular magnetic field and the perpendicular proton bulk velocity fluctuations exhibit power-law scaling with spectral indices of −5/3 and −3/2, respectively. This behavior is extended over a full decade in wavevectors and is very stable in time. A transition is observed around proton scales. At sub-ion scales, both spectra steepen, with the former still following a power law with a spectral index of ∼ −3. A −2.8 slope is observed in the density and parallel magnetic fluctuations, highlighting the presence of compressive effects at kinetic scales. The spectrum of the perpendicular electric fluctuations follows that of the proton bulk velocity at MHD scales, and flattens at small scales. All these features, which we carefully tested against variations of many parameters, are in good agreement with solar wind observations. The turbulent cascade leads to on overall proton energization with similar heating rates in the parallel and perpendicular directions. While the parallel proton heating is found to be independent on the resistivity, the number of particles per cell and the resolution employed, the perpendicular proton temperature strongly depends on these parameters.

A new path for the generation of a sub-ion scale cascade in collisionless space and astrophysical plasma turbulence, triggered by magnetic reconnection, is uncovered by means of high-resolution two-dimensional hybrid-kinetic simulations employing two complementary approaches, Lagrangian and Eulerian, and different driving mechanisms. The simulation results provide clear numerical evidences that the development of powerlaw energy spectra below the so-called ion break occurs as soon as the first magnetic reconnection events take place, regardless of the actual state of the turbulent cascade at MHD scales. In both simulations, the reconnection-mediated small-scale energy spectrum of parallel magnetic fluctuations exhibits a very stable spectral slope of ∼ −2.8, whether or not a large-scale turbulent cascade has already fully developed. Once a quasi-stationary turbulent state is achieved, the spectrum of the total magnetic fluctuations settles towards a spectral index of −5/3 in the MHD range and of ∼ −3 at sub-ion scales.

21 Corresponding author: J. E. Stawarz j.stawarz@imperial.ac.uk 2 Turbulent plasmas generate intense current structures, which have long been sug-22 gested as magnetic reconnection sites. Recent Magnetospheric Multiscale (MMS) ob-23 servations in Earth's magnetosheath revealed a novel form of reconnection where the 24 dynamics only couple to electrons, without ion involvement. It was suggested that 25 such dynamics were driven by magnetosheath turbulence. In this study, the fluctua-26 tions are examined to determine the properties of the turbulence and if a signature of 27 reconnection is present in the turbulence statistics. The study reveals statistical prop-28 erties consistent with plasma turbulence with a correlation length of ∼ 10 ion inertial 29 lengths. When reconnection is more prevalent, a steepening of the magnetic spectrum 30 occurs at the length scale of the reconnecting current sheets. The statistics of intense 31 currents suggest the prevalence of electron-scale current sheets favorable for electron-32 reconnection. The results support the hypothesis that electron-reconnection is driven 33 by turbulence and highlight diagnostics that may provide insight into reconnection in 34 other turbulent plasmas.35

A dynamical vectorial equation for homogeneous incompressible Hall-MHD turbulence together with the exact scaling law for third-order correlation tensors, analogous to that for the incompressible MHD, is rederived and applied to the results of two-dimensional hybrid simulations of plasma turbulence. At large (MHD) scales the simulations exhibits a clear inertial range where the MHD dynamic law is valid. In the sub-ion range the cascade continues via the Hall term but the dynamic law derived in the framework of incompressible Hall MHD equations is obtained only in a low plasma beta simulation. For a higher beta plasma the cascade rate decreases in the sub-ion range and the change becomes more pronounced as the plasma beta increases. This break in the cascade flux can be ascribed to non thermal (kinetic) features or to others terms in the dynamical equation that are not included in the Hall-MHD incompressible approximation.

We investigate properties of the ion-scale spectral break of solar wind turbulence by means of twodimensional high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box and add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma β, distributed over three orders of magnitude, from 0.01 to 10. In all the cases, once turbulence is fully developed, we observe a power-law spectrum of the fluctuating magnetic field on large scales (in the inertial range) with a spectral index close to −5/3, while in the sub-ion range we observe another power-law spectrum with a spectral index systematically varying with β (from around −3.6 for small values to around −2.9 for large ones). The two ranges are separated by a spectral break around ion scales. The length scale at which this transition occurs is found to be proportional to the ion inertial length, d i , for β 1 and to the ion gyroradius, ρ i = d i √ β, for β 1, i.e., to the larger between the two scales in both the extreme regimes. For intermediate cases, i.e., β ∼ 1, a combination of the two scales is involved. We infer an empiric relation for the dependency of the spectral break on β that provides a good fit over the whole range of values. We compare our results with in situ observations in the solar wind and suggest possible explanations for such a behavior.

We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several properties of the observed solar wind fluctuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magneto-hydrodynamic (MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic compressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind large-scale MHD fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized Ohm's law.

Properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size 3D hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfvénic fluctuations rapidly develop a strongly anisotropic turbulent cascade, mainly in the direction perpendicular to the ambient magnetic field. The omnidirectional magnetic field spectrum shows a double power-law behavior over almost two decades in wavenumber, with a Kolmogorov-like index at large scales, a spectral break around ion scales, and a steepening at sub-ion scales. Power laws are also observed in the spectra of the ion bulk velocity, density, and electric field, at both magnetohydrodynamic (MHD) and kinetic scales. Despite the complex structure, the omnidirectional spectra of all fields at ion and sub-ion scales are in remarkable quantitative agreement with those of a 2D simulation with similar physical parameters. This provides a partial, a posteriori validation of the 2D approximation at kinetic scales. Conversely, at MHD scales, the spectra of the density and of the velocity (and, consequently, of the electric field) exhibit differences between the 2D and 3D cases. Although they can be partly ascribed to the lower spatial resolution, the main reason is likely the larger importance of compressible effects in the full 3D geometry. Our findings are also in remarkable quantitative agreement with solar wind observations.

We investigate the properties of plasma turbulence by means of two-dimensional Hallmagnetohydrodynamic (HMHD) and hybrid particle-in-cell (HPIC) numerical simulations. We find that HMHD simulations exhibit spectral properties that are in most cases in agreement with the results of the HPIC simulations and with solar wind observations. The energy spectra of magnetic fluctuations exhibit a double power-law with spectral index −5/3 at MHD scales and −3 at kinetic scales, while for velocity fluctuations the spectral index is −3/2 at MHD scales. The break between the MHD and the kinetic scales occurs at the same scale in both simulations. In the MHD range the slopes of the total energy and residual energy spectra satisfy a fast Alfvén-dynamo balance. The development of a turbulent cascade is concurrently characterized by magnetic reconnection events taking place in thin current sheets that form between large eddies. A statistical analysis reveals that reconnection is qualitatively the same and fast in both the HMHD and HPIC models, characterized by inverse reconnection rates much smaller than the characteristic large-eddy nonlinear time. The agreement extends to other statistical properties, such us the kurtosis of the magnetic field. Moreover, the observation of a direct energy transfer from the large vortices to the small sub-ion scales triggered by magnetic reconnection, further supports the existence of a reconnection-mediated turbulent regime at kinetic scales. We conclude that the Hall-MHD fluid description captures to a large extent the transition of the turbulent cascade between the large MHD scales and the sub-ion scales.