Energy conversion via reconnecting current sheets is common in space and astrophysical plasmas. Frequently, current sheets disrupt at multiple reconnection sites, leading to the formation of plasmoid structures between sites, which might affect energy conversion. We present in situ evidence of the firehose instability in multiple reconnection in the Earth's magnetotail. The observed proton beams accelerated in the direction parallel to magnetic field and ion-scale fluctuations of whistler type imply the development of firehose instability between two active reconnection sites. The linear wave dispersion relation, estimated for the measured plasma parameters, indicates a positive growth rate of firehose-related electromagnetic fluctuations. Simulations of temporal evolution of the observed multiple reconnection by using a 2.5D implicit particle-in-cell code show that, as the plasmoid formed between two reconnection sites evolves, the plasma at its edge becomes anisotropic and overcomes the firehose marginal stability threshold, leading to the generation of magnetic field fluctuations. The combined results of observations and simulations suggest that the firehose instability, operating between reconnection sites, converts plasma kinetic energy into energy of magnetic field fluctuations, counteracting the conversion of magnetic energy into plasma energy occurring at reconnection sites. This suggests that magnetic energy conversion in multiple reconnection can be less efficient than in the case of the single-site reconnection.

We present a Vlasov-DArwin numerical code (ViDA) specifically designed to address plasma physics problems, where small-scale high accuracy is requested even during the non linear regime to guarantee a clean description of the plasma dynamics at fine spatial scales. The algorithm provides a low-noise description of proton and electron kinetic dynamics, by splitting in time the multi-advection Vlasov equation in phase space. Maxwell equations for the electric and magnetic fields are reorganized according to Darwin approximation to remove light waves. Several numerical tests show that ViDA successfully reproduces the propagation of linear and nonlinear waves and captures the physics of magnetic reconnection. We also discuss preliminary tests of the parallelization algorithm efficiency, performed at CINECA on the Marconi-KNL cluster. ViDA will allow to run Eulerian simulations of a non-relativistic fully-kinetic collisionless plasma and it is expected to provide relevant insights on important problems of plasma astrophysics such as, for instance, the development of the turbulent cascade at electron scales and the structure and dynamics of electron-scale magnetic reconnection, such as the electron diffusion region.

In weakly collisional space plasmas, the turbulent cascade provides most of the energy that is dissipated at small scales by various kinetic processes. Understanding the characteristics of such dissipative mechanisms requires the accurate knowledge of the fluctuations that make energy available for conversion at small scales, as different dissipation processes are triggered by fluctuations of a different nature. The scaling properties of different energy channels are estimated here using a proxy of the local energy transfer, based on the third-order moment scaling law for magnetohydrodynamic turbulence. In particular, the sign-singularity analysis was used to explore the scaling properties of the alternating positive-negative energy fluxes, thus providing information on the structure and topology of such fluxes for each of the different type of fluctuations. The results show the highly complex geometrical nature of the flux, and that the local contributions associated with energy and cross-helicity nonlinear transfer have similar scaling properties. Consequently, the fractal properties of current and vorticity structures are similar to those of the Alfv\'enic fluctuations.

Particles are accelerated to very high, non-thermal energies in solar and space plasma environments. While energy spectra of accelerated electrons often exhibit a power law, it remains unclear how electrons are accelerated to high energies and what processes determine the power-law index $\delta$. Here, we review previous observations of the power-law index $\delta$ in a variety of different plasma environments with a particular focus on sub-relativistic electrons. It appears that in regions more closely related to magnetic reconnection (such as the `above-the-looptop' solar hard X-ray source and the plasma sheet in Earth's magnetotail), the spectra are typically soft ($\delta \gtrsim$ 4). This is in contrast to the typically hard spectra ($\delta \lesssim$ 4) that are observed in coincidence with shocks. The difference implies that shocks are more efficient in producing a larger non-thermal fraction of electron energies when compared to magnetic reconnection. A caveat is that during active times in Earth's magnetotail, $\delta$ values seem spatially uniform in the plasma sheet, while power-law distributions still exist even in quiet times. The role of magnetotail reconnection in the electron power-law formation could therefore be confounded with these background conditions. Because different regions have been studied with different instrumentations and methodologies, we point out a need for more systematic and coordinated studies of power-law distributions for a better understanding of possible scaling laws in particle acceleration as well as their universality.

Low frequency turbulence in Saturn's magnetosheath is investigated using in-situ measurements of the Cassini spacecraft. Focus is put on the magnetic energy spectra computed in the frequency range $\sim[10^{-4}, 1]$Hz. A set of 42 time intervals in the magnetosheath were analyzed and three main results that contrast with known features of solar wind turbulence are reported: 1) The magnetic energy spectra showed a $\sim f^{-1}$ scaling at MHD scales followed by an $\sim f^{-2.6}$ scaling at the sub-ion scales without forming the so-called inertial range; 2) The magnetic compressibility and the cross-correlation between the parallel component of the magnetic field and density fluctuations $ C(\delta n,\delta B_{||}) $ indicates the dominance of the compressible magnetosonic slow-like modes at MHD scales rather than the Alfv\'en mode; 3) Higher order statistics revealed a monofractal (resp. multifractal) behaviour of the turbulent flow downstream of a quasi-perpendicular (resp. quasi-parallel) shock at the sub-ion scales. Implications of these results on theoretical modeling of space plasma turbulence are discussed.

To investigate the internal structure of the magnetopause with spacecraft data, it is crucial to be able to determine its normal direction and to convert the measured time series into spatial profiles. We propose here a new single-spacecraft method, called the BV method, to reach these two objectives. Its name indicates that the method uses a combination of the magnetic field (B) and velocity (V) data. The method is tested on simulation and Cluster data, and a short overview of the possible products is given. We discuss its assumptions and show that it can bring a valuable improvement with respect to previous methods.

The magnetopause boundary implies two kinds of variations: a density/ temperature gradient and a magnetic field rotation. These two kinds are always observed in a close vicinity of each other, if not inseparably mixed. We present a case study from the Cluster data where the two are clearly separated and investigate the natures of both layers. We evidence that the first one is a slow shock while the second is a rotational discontinuity. The interaction between these two kinds of discontinuities is then studied with the help of 1,5-D magnetohydrodynamics simulations. The comparison with the data is quite positive and leads to think that most of the generic properties of the magnetopause may be interpreted in this sense.

Electrons can be accelerated to ultrarelativistic energies at strong (high-Mach number) collisionless shock waves that form when stellar debris rapidly expands after a supernova. Collisionless shock waves also form in the flow of particles from the Sun (the solar wind), and extensive spacecraft observations have established that electron acceleration at these shocks is effectively absent whenever the upstream magnetic field is roughly parallel to the shock surface normal (quasi-parallel conditions). However, it is unclear whether this magnetic dependence of electron acceleration also applies to the far stronger shocks around young supernova remnants, where local magnetic conditions are poorly understood. Here we present Cassini spacecraft observations of an unusually strong solar system shock wave (Saturn's bow shock) where significant local electron acceleration has been confirmed under quasi-parallel magnetic conditions for the first time, contradicting the established magnetic dependence of electron acceleration at solar system shocks. Furthermore, the acceleration led to electrons at relativistic energies (~MeV), comparable to the highest energies ever attributed to shock-acceleration in the solar wind. These observations demonstrate that at high-Mach numbers, like those of young supernova remnant shocks, quasi-parallel shocks become considerably more effective electron accelerators.