Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the energization of transrelativistic electrons in turbulence using hybrid particle-in-cell, which provide a realistic model of Alfv\'{e}nic turbulence from MHD to sub-ion scales, and test particle simulations for electrons. We find that, depending on the electron initial energy and turbulence strength, electrons may undergo a fast and efficient phase of energization due to the magnetic curvature drift during the time they are trapped in dynamic magnetic structures. In addition, electrons are accelerated stochastically which is a slower process that yields lower maximum energies. The combined effect of these two processes determines the overall electron acceleration. With appropriate turbulence parameters, we find that superthermal electrons can be accelerated up to relativistic energies. For example, with heliospheric parameters and a relatively high turbulence level, rapid energization to MeV energies is possible.

Turbulence develops in any stressed flow when the scales of the forcing are much larger than those of the dissipation. In neutral fluids, it consists of chaotic motions in physical space but with a universal energy spectrum in Fourier space. Intermittency (non-Gaussian statistics of fluctuations) is another general property and it is related to the presence of coherent structures. Space plasmas are turbulent as well. Here, we focus on the kinetic plasma scales, which are not yet well understood. We address the following fundamental questions: (1) Do the turbulent fluctuations at kinetic scales form a universal spectrum? and (2) What is the nature of the fluctuations? Using measurements in the solar wind we show that the magnetic spectra of kinetic turbulence at 0.3, 0.6 and 0.9 AU from the Sun have the same shape as the ones close to the Earth orbit at 1 AU, indicating universality of the phenomenon. The fluctuations, which form this spectrum, are typically non-linearly interacting eddies that tend to generate magnetic filaments.

Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfv\'enic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and electron betas, turbulence strength, and injection scale) exhibits a quantitative agreement on the spectral, intermittency, and cascade-rate properties with in situ observations, despite the different driving mechanisms. Such agreement demonstrates a certain universality of the turbulent cascade from magnetohydrodynamic (MHD) to sub-ion scales, whose properties are mainly determined by the selected parameters, also indicating that the KH instability-driven turbulence has a quasi-2D nature. The validity of the Taylor hypothesis in the sub-ion spatial range suggests that the fluctuations at sub-ion scales have predominantly low frequencies, consistent with a kinetic Alfv\'en wave-like nature or with quasi-static structures. Finally, the third-order structure function analysis indicates that the cascade rate of the turbulence generated by a KH event in the magnetopause is an order of magnitude larger than in the ambient magnetosheath.

The relationship between a decaying plasma turbulence and proton fire hose instabilities in a slowly expanding plasma is investigated using three-dimensional (3-D) hybrid expanding box simulations. We impose an initial ambient magnetic field along the radial direction, and we start with an isotropic spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations with zero cross-helicity. A turbulent cascade rapidly develops and leads to a weak proton heating that is not sufficient to overcome the expansion-driven perpendicular cooling. The plasma system eventually drives the parallel and oblique fire hose instabilities that generate quasi-monochromatic wave packets that reduce the proton temperature anisotropy. The fire hose wave activity has a low amplitude with wave vectors quasi-parallel/oblique with respect to the ambient magnetic field outside of the region dominated by the turbulent cascade and is discernible in one-dimensional power spectra taken only in the direction quasi-parallel/oblique with respect to the ambient magnetic field; at quasi-perpendicular angles the wave activity is hidden by the turbulent background. These waves are partly reabsorbed by protons and partly couple to and participate in the turbulent cascade. Their presence reduces kurtosis, a measure of intermittency, and the Shannon entropy but increases the Jensen-Shannon complexity of magnetic fluctuations; these changes are weak and anisotropic with respect to the ambient magnetic field and it's not clear if they can be used to indirectly discern the presence of instability-driven waves.

By means of three dimensional high-resolution hybrid simulations we study the properties of the magnetic field spectral anisotropies near and beyond ion kinetic scales. By using both a Fourier analysis and a local analysis based on multi-point 2nd-order structure function techniques, we show that the anisotropy observed is less than what expected by standard wave normal modes turbulence theories although the non linear energy transfer is still in the perpendicular direction, only advected in the parallel direction as expected balancing the non-linear energy transfer time and the decorrelation time. Such result can be explained by a phenomenological model based on the formation of strong intermittent two-dimensional structures in the plane perpendicular to the local mean field that have some prescribed aspect ratio eventually depending on the scale. This model support the idea that small scales structures, such as reconnecting current sheets, contribute significantly to the formation of the turbulent cascade at kinetic scales.

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.

Measurements and simulations of inertial compressive turbulence in the solar wind are characterized by anti-correlated magnetic fluctuations parallel to the mean field and density structures. This signature has been interpreted as observational evidence for non-propagating pressure balanced structures (PBS), kinetic ion acoustic waves, as well as the MHD slow-mode. Given the high damping rates of parallel propagating compressive fluctuations, their ubiquity in satellite observations is surprising, and suggestive of a local driving process. One possible candidate for the generation of compressive fluctuations in the solar wind is Alfv\'en wave parametric instability. Here we test the parametric decay process as a source of compressive waves in the solar wind by comparing the collisionless damping rates of compressive fluctuations with the growth rates of the parametric decay instability daughter waves. Our results suggest that generation of compressive waves through parametric decay is overdamped at 1 AU, but that the presence of slow-mode like density fluctuations is correlated with the parametric decay of Alfv\'en waves.

We investigate the spectral properties of plasma turbulence from fluid to sub-ion scales by means of high-resolution three-dimensional (3D) numerical simulations performed with the hybrid particle-in-cell (HPIC) code CAMELIA. We produce extended turbulent spectra with well-defined power laws for the magnetic, ion bulk velocity, density, and electric fluctuations. The present results are in good agreement with previous two-dimensional (2D) HPIC simulations, especially in the kinetic range of scales, and reproduce several features observed in solar wind spectra. By providing scaling tests on many different architectures and convergence studies, we prove CAMELIA to represent a very efficient, accurate and reliable tool for investigating the develpoment of the turbulent cascade in the solar wind, being able to cover simultaneously several decades in wavenumber, also in 3D.

Spectral properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size three-dimensional (3D) hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfv\`enic 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, both at magnetohydrodynamic (MHD) and at 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 two-dimensional (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 a full geometry. Our findings are also in remarkable quantitative agreement with solar wind observations.

A new path for the generation of a sub-ion scale cascade in collisionless 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 power-law 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.

The relationship between a decaying strong turbulence and the mirror instability in a slowly expanding plasma is investigated using two-dimensional hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at sub-ion scales. The imposed expansion (taking a strictly transverse ambient magnetic field) leads to generation of an important perpendicular proton temperature anisotropy that eventually drives the mirror instability. This instability generates large-amplitude, nonpropagating, compressible, pressure-balanced magnetic structures in a form of magnetic enhancements/humps that reduce the perpendicular temperature anisotropy.

Kinetic-range turbulence in magnetized plasmas and, in particular, in the context of solar-wind turbulence has been extensively investigated over the past decades via numerical simulations. Among others, one of the widely adopted reduced plasma model is the so-called hybrid-kinetic model, where the ions are fully kinetic and the electrons are treated as a neutralizing (inertial or massless) fluid. Within the same model, different numerical methods and/or approaches to turbulence development have been employed. In the present work, we present a comparison between two-dimensional hybrid-kinetic simulations of plasma turbulence obtained with two complementary approaches spanning about two decades in wavenumber - from MHD inertial range to scales well below the ion gyroradius - with a state-of-the-art accuracy. One approach employs hybrid particle-in-cell (HPIC) simulations of freely-decaying Alfv\'enic turbulence, whereas the other consists of Eulerian hybrid Vlasov-Maxwell (HVM) simulations of turbulence continuously driven with partially-compressible large-scale fluctuations. Despite the completely different initialization and injection/drive at large scales, the same properties of turbulent fluctuations at $k_\perp\rho_i\gtrsim1$ are observed. The system indeed self-consistently "reprocesses" the turbulent fluctuations while they are cascading towards smaller and smaller scales, in a way which actually depends on the plasma beta parameter. Small-scale turbulence has been found to be mainly populated by kinetic Alfv\'en wave (KAW) fluctuations for $\beta\geq1$, whereas KAW fluctuations are only sub-dominant for low-$\beta$.

Using two-dimensional hybrid expanding box simulations we study the competition between the continuously driven parallel proton temperature anisotropy and fire hose instabilities in collisionless homogeneous plasmas. For quasi radial ambient magnetic field the expansion drives $T_{\mathrm{p}\|}>T_{\mathrm{p}\perp}$ and the system becomes eventually unstable with respect to the dominant parallel fire hose instability. This instability is generally unable to counteract the induced anisotropization and the system typically becomes unstable with respect to the oblique fire hose instability later on. The oblique instability efficiently reduces the anisotropy and the system rapidly stabilizes while a significant part of the generated electromagnetic fluctuations is damped to protons. As long as the magnetic field is in the quasi radial direction, this evolution repeats itself and the electromagnetic fluctuations accumulate. For sufficiently oblique magnetic field the expansion drives $T_{\mathrm{p}\perp}>T_{\mathrm{p}\|}$ and brings the system to the stable region with respect to the fire hose instabilities.

We investigate properties of the ion-scale spectral break of solar wind turbulence 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 add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma beta, 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 $\beta$ (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 $\beta \ll 1$ and to the ion gyroradius, $\rho_i = d_i\sqrt{\beta}$, for $\beta \gg 1$, i.e., to the larger between the two scales in both the extreme regimes. For intermediate cases, i.e., $\beta \sim 1$, a combination of the two scales is involved. We infer an empiric relation for the dependency of the spectral break on $\beta$ 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.

Using one-dimensional hybrid expanding box model we investigate properties of the solar wind in the outer heliosphere. We assume a proton-electron plasma with a strictly transverse ambient magnetic field and, beside the expansion, we take into account influence of a continuous injection of cold pick-up protons through the charge-exchange process between the solar wind protons and hydrogen of interstellar origin. The injected cold pick-up protons form a ring distribution function that rapidly becomes unstable and generate Alfv\'en cyclotron waves. The Alfv\'en cyclotron waves scatter pick-up protons to a spherical shell distribution function that thickens over that time owing to the expansion-driven cooling. The Alf\'ven cyclotron waves heat solar wind protons in the perpendicular direction (with respect to the ambient magnetic field) through the cyclotron resonance. At later times, the Alfv\'en cyclotron waves become parametrically unstable and the generated ion acoustic waves heat protons in the parallel direction through the Landau resonance. The resulting heating of the solar wind protons is efficient on the expansion time scale.

Proton temperature anisotropies between the directions parallel and perpendicular to the mean magnetic field are usually observed in the solar wind plasma. Here, we employ a high-resolution hybrid particle-in-cell simulation in order to investigate the relation between spatial properties of the proton temperature and the peaks in the current density and in the flow vorticity. Our results indicate that, although regions where the proton temperature is enhanced and temperature anisotropies are larger correspond approximately to regions where many thin current sheets form, no firm quantitative evidence supports the idea of a direct causality between the two phenomena. On the other hand, quite a clear correlation between the behavior of the proton temperature and the out-of-plane vorticity is obtained.

Supercritical collisionless perpendicular shocks have an average macrostructure determined primarily by the dynamics of ions specularly reflected at the magnetic ramp. Within the overall macrostructure, instabilities, both linear and nonlinear, generate fluctuations and microstructure. To identify the sources of such microstructure, high-resolution two- and three-dimensional simulations have been carried out using the hybrid method, wherein the ions are treated as particles and the electron response is modelled as a massless fluid. We confirm the results of earlier 2-D simulations showing both field-parallel aligned propagating fluctuations and fluctuations carried by the reflected-gyrating ions. In addition, it is shown that, for 2-D simulations of the shock coplanarity plane, the presence of short-wavelength fluctuations in all magnetic components is associated with the ion Weibel instability driven at the upstream edge of the foot by the reflected-gyrating ions. In 3-D simulations we show for the first time that the dominant microstructure is due to a coupling between field-parallel propagating fluctuations in the ramp and the motion of the reflected ions. This results in a pattern of fluctuations counter-propagating across the surface of the shock at an angle inclined to the magnetic field direction, due to a combination of field-parallel motion at the Alfv\'en speed of the ramp, and motion in the sense of gyration of the reflected ions.

The relationship between a decaying strong turbulence and kinetic instabilities in a slowly expanding plasma is investigated using two-dimensional (2-D) hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfv\'enic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at ion scales. The turbulent cascade leads to an overall anisotropic proton heating, protons are heated in the perpendicular direction, and, initially, also in the parallel direction. The imposed expansion leads to generation of a large parallel proton temperature anisotropy which is at later stages partly reduced by turbulence. The turbulent heating is not sufficient to overcome the expansion-driven perpendicular cooling and the system eventually drives the oblique firehose instability in a form of localized nonlinear wave packets which efficiently reduce the parallel temperature anisotropy. This work demonstrates that kinetic instabilities may coexist with strong plasma turbulence even in a constrained 2-D regime.

We investigate properties of a solar wind-like plasma including a secondary alpha particle population exhibiting a parallel temperature anisotropy with respect to the background magnetic field, using linear and quasi-linear predictions and by means of one-dimensional hybrid simulations. We show that anisotropic alpha particles can drive a parallel fire hose instability analogous to that generated by protons, but that, remarkably, the instability can be triggered also when the parallel plasma beta of alpha particles is below unity. The wave activity generated by the alpha anisotropy affects the evolution of the more abundant protons, leading to their anisotropic heating. When both ion species have sufficient parallel anisotropies both of them can drive the instability, and we observe generation of two distinct peaks in the spectra of the fluctuations, with longer wavelengths associated to alphas and shorter ones to protons. If a non-zero relative drift is present, the unstable modes propagate preferentially in the direction of the drift associated with the unstable species. The generated waves scatter particles and reduce their temperature anisotropy to marginally stable state, and, moreover, they significantly reduce the relative drift between the two ion populations. The coexistence of modes excited by both species leads to saturation of the plasma in distinct regions of the beta/anisotropy parameter space for protons and alpha particles, in good agreement with in situ solar wind observations. Our results confirm that fire hose instabilities are likely at work in the solar wind and limit the anisotropy of different ion species in the plasma.

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 large-scale 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.

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.

The measure of the third-order structure function, Y, is employed in the solar wind to compute the cascade rate of turbulence. In the absence of a mean field B0=0, Y is expected to be isotropic (radial) and independent of the direction of increments, so its measure yields directly the cascade rate. For turbulence with mean field, as in the solar wind, Y is expected to become more two dimensional (2D), that is, to have larger perpendicular components, loosing the above simple symmetry. To get the cascade rate one should compute the flux of Y, which is not feasible with single-spacecraft data, thus measurements rely upon assumptions about the unknown symmetry. We use direct numerical simulations (DNS) of magneto-hydrodynamic (MHD) turbulence to characterize the anisotropy of Y. We find that for strong guide field B0=5 the degree of two-dimensionalization depends on the relative importance of shear and pseudo polarizations (the two components of an Alfv\'en mode in incompressible MHD). The anisotropy also shows up in the inertial range. The more Y is 2D, the more the inertial range extent differs along parallel and perpendicular directions. We finally test the two methods employed in observations and find that the so-obtained cascade rate may depend on the angle between B0 and the direction of increments. Both methods yield a vanishing cascade rate along the parallel direction, contrary to observations, suggesting a weaker anisotropy of solar wind turbulence compared to our DNS. This could be due to a weaker mean field and/or to solar wind expansion.

The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime.