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.