The mechanism of heating for hot, dilute, and turbulent plasmas represents a long-standing problem in space physics, whose implications concern both near-Earth environments and astrophysical systems. In order to explore the possible role of interparticle collisions, simulations of plasma turbulence -- in both collisionless and weakly collisional regimes -- have been compared by adopting Eulerian Hybrid Boltzmann-Maxwell simulations, being proton-proton collisions explicitly introduced through the nonlinear Dougherty operator. Although collisions do not significantly influence the statistical characteristics of the turbulence, they dissipate nonthermal features in the proton distribution function and suppress the enstrophy/entropy cascade in the velocity space, damping the spectral transfer toward large Hermite modes. This enstrophy dissipation is particularly effective in regions where the plasma distribution function is strongly distorted, suggesting that collisional effects are enhanced by fine velocity-space structures. A qualitative connection between the turbulent energy cascade in fluids and the enstrophy cascade in plasmas has been established, opening a new path to the understanding of astrophysical plasma turbulence

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.

Plasma turbulence is investigated using high-resolution ion velocity distributions measured by the Magnetospheric Multiscale Mission (MMS) in the Earth's magnetosheath. The particle distribution is highly structured, suggesting a cascade-like process in velocity space. This complex velocity space structure is investigated using a three-dimensional Hermite transform that reveals a power law distribution of moments. In analogy to hydrodynamics, a Kolmogorov approach leads directly to a range of predictions for this phase-space cascade. The scaling theory is in agreement with observations, suggesting a new path for the study of plasma turbulence in weakly collisional space and astrophysical plasmas.

Using high resolution Cluster satellite observations, we show that the turbulent solar wind is populated by magnetic discontinuities at different scales, going from proton down to electron scales. The structure of these layers resembles the Harris equilibrium profile in plasmas. Using a multi-dimensional intermittency technique, we show that these structures are connected through the scales. Supported by numerical simulations of magnetic reconnection, we show that observations are consistent with a scenario where many current layers develop in turbulence, and where the outflow of these reconnection events are characterized by complex sub-proton networks of secondary islands, in a self-similar way. The present work establishes that the picture of "reconnection in turbulence" and "turbulent reconnection", separately invoked as ubiquitous, coexist in space plasmas.

We study the radial evolution of intermittency of density fluctuations in the fast solar wind. The study is performed analyzing the plasma density measurements provided by Helios 2 in the inner heliosphere between $0.3$ and $0.9$ AU. The analysis is carried out by means of a complete set of diagnostic tools, including the flatness factor at different time scales to estimate intermittency, the Kolmogorov-Smirnov test to estimate the degree of intermittency, and the Fourier transform to estimate the power spectral densities of these fluctuations. Density fluctuations within fast wind are rather intermittent and their level of intermittency, together with the amplitude of intermittent events, decreases with distance from the Sun, at odds with intermittency of both magnetic field and all the other plasma parameters. Furthermore, the intermittent events are strongly correlated, exhibiting temporal clustering. This indicates that the mechanism underlying their generation departs from a time-varying Poisson process. A remarkable, qualitative similarity with the behavior of plasma density fluctuations obtained from a numerical study of the nonlinear evolution of parametric instability in the solar wind supports the idea that this mechanism has an important role in governing density fluctuations in the inner heliosphere.

Using the magnetohydrodynamic (MHD) description, we develop a nonlinear dynamo model that couples the evolution of the large scale magnetic field with turbulent dynamics of the plasma at small scale by electromotive force (e.m.f.) in the induction equation at large scale. The nonlinear behavior of the plasma at small scale is described by using a MHD shell model for velocity field and magnetic field fluctuations.The shell model allow to study this problem in a large parameter regime which characterizes the dynamo phenomenon in many natural systems and which is beyond the power of supercomputers at today. Under specific conditions of the plasma turbulent state, the field fluctuations at small scales are able to trigger the dynamo instability. We study this transition considering the stability curve which shows a strong decrease in the critical magnetic Reynolds number for increasing inverse magnetic Prandlt number $\textrm{Pm}^{-1}$ in the range $[10^{-6},1]$ and slows an increase in the range $[1,10^{8}]$. We also obtain hysteretic behavior across the dynamo boundary reveling the subcritical nature of this transition. The system, undergoing this transition, can reach different dynamo regimes, depending on Reynolds numbers of the plasma flow. This shows the critical role that the turbulence plays in the dynamo phenomenon. In particular the model is able to reproduce the dynamical situation in which the large-scale magnetic field jumps between two states which represent the opposite polarities of the magnetic field, reproducing the magnetic reversals as observed in geomagnetic dynamo and in the VKS experiments.

Cluster observations in the near-Earth magnetotail have shown that sometimes the current sheet is bifurcated, i.e. it is divided in two layers. The influence of magnetic turbulence on ion motion in this region is investigated by numerical simulation, taking into account the presence of both protons and oxygen ions. The magnetotail current sheet is modeled as a magnetic field reversal with a normal magnetic field component $B_n$, plus a three-dimensional spectrum of magnetic fluctuations $\delta {\bf B}$, which represents the observed magnetic turbulence. The dawn-dusk electric field E$_y$ is also included. A test particle simulation is performed using different values of $\delta {\bf B}$, E$_y$ and injecting two different species of particles, O$^+$ ions and protons. O$^+$ ions can support the formation of a double current layer both in the absence and for large values of magnetic fluctuations ($\delta B/B_0 = 0.0$ and $\delta B/B_0 \geq 0.4$, where B$_0$ is the constant magnetic field in the magnetospheric lobes).

Magnetic fluctuations in the solar wind are distributed according to Kolmogorov's power law $f^{-5/3}$ below the ion cyclotron frequency $f_{ci}$. Above this frequency, the observed steeper power law is usually interpreted in two different ways: a dissipative range of the solar wind turbulence or another turbulent cascade, the nature of which is still an open question. Using the Cluster magnetic data we show that after the spectral break the intermittency increases toward higher frequencies, indicating the presence of non-linear interactions inherent to a new inertial range and not to the dissipative range. At the same time the level of compressible fluctuations raises. We show that the energy transfer rate and intermittency are sensitive to the level of compressibility of the magnetic fluctuations within the small scale inertial range. We conjecture that the time needed to establish this inertial range is shorter than the eddy-turnover time, and is related to dispersive effects. A simple phenomenological model, based on the compressible Hall MHD, predicts the magnetic spectrum $\sim k^{-7/3+2\alpha}$, which depends on the degree of plasma compression $\alpha$.

We show in this article direct evidence for the presence of an inertial energy cascade, the most characteristic signature of hydromagnetic turbulence (MHD), in the solar wind as observed by the Ulysses spacecraft. After a brief rederivation of the equivalent of Yaglom's law for MHD turbulence, we show that a linear relation is indeed observed for the scaling of mixed third order structure functions involving Els\"asser variables. This experimental result, confirming the prescription stemming from a theorem for MHD turbulence, firmly establishes the turbulent character of low-frequency velocity and magnetic field fluctuations in the solar wind plasma.

Often in nature the temporal distribution of inhomogeneous stochastic point processes can be modeled as a realization of renewal Poisson processes with a variable rate. Here we investigate one of the classical examples, namely the temporal distribution of polarity reversals of the geomagnetic field. In spite of the commonly used underlying hypothesis, we show that this process strongly departs from a Poisson statistics, the origin of this failure stemming from the presence of temporal clustering. We find that a Levy statistics is able to reproduce paleomagnetic data, thus suggesting the presence of long-range correlations in the underlying dynamo process.

We model the hydrodynamic evolution of the plasma confined in a coronal loop, 30000 km long, subject to the heating of nanoflares due to intermittent magnetic dissipative events in the MHD turbulence produced by loop footpoint motions. We use the time-dependent distribution of energy dissipation along the loop obtained from a hybrid shell model, occurring for a magnetic field of about 10 G in corona; the relevant heating per unit volume along the loop is used in the Palermo-Harvard loop plasma hydrodynamic model. We describe the results focussing on the effects produced by the most intense heat pulses, which lead to loop temperatures between 1 and 1.5 MK.

The spatio-temporal dynamics of the solar photosphere is studied by performing a Proper Orthogonal Decomposition (POD) of line of sight velocity fields computed from high resolution data coming from the MDI/SOHO instrument. Using this technique, we are able to identify and characterize the different dynamical regimes acting in the system. Low frequency oscillations, with frequencies in the range 20-130 microHz, dominate the most energetic POD modes (excluding solar rotation), and are characterized by spatial patterns with typical scales of about 3 Mm. Patterns with larger typical scales of 10 Mm, are associated to p-modes oscillations at frequencies of about 3000 microHz.

In spite of a large number of papers dedicated to study MHD turbulence in the solar wind there are still some simple questions which have never been sufficiently addressed like: a)do we really know how the magnetic field vector orientation fluctuates in space? b) what is the statistics followed by the orientation of the vector itself? c) does the statistics change as the wind expands into the interplanetary space? A better understanding of these points can help us to better characterize the nature of interplanetary fluctuations and can provide useful hints to investigators who try to numerically simulate MHD turbulence. This work follows a recent paper presented by the same authors. This work follows a recent paper presented by some of the authors which shows that these fluctuations might resemble a sort of random walk governed by a Truncated Leevy Flight statistics. However, the limited statistics used in that paper did not allow final conclusions but only speculative hypotheses. In this work we aim to address the same problem using a more robust statistics which on one hand forces us not to consider velocity fluctuations but, on the other hand allows us to establish the nature of the governing statistics of magnetic fluctuations with more confidence. In addition, we show how features similar to those found in the present statistical analysis for the fast speed streams of solar wind, are qualitatively recovered in numerical simulations of the parametric instability. This might offer an alternative viewpoint for interpreting the questions raised above.

The observations of magnetic field variations as a signature of flaring activity is one of the main goal in solar physics. Some efforts in the past give apparently no unambiguous observations of changes. We observed that the scaling laws of the current helicity inside a given flaring active region change clearly and abruptly in correspondence with the eruption of big flares at the top of that active region. Comparison with numerical simulations of MHD equations, indicates that the change of scaling behavior in the current helicity, seems to be associated to a topological reorganization of the footpoint of the magnetic field loop, namely to dissipation of small scales structures in turbulence. It is evident that the possibility of forecasting in real time high energy flares, even if partially, has a wide practical interest to prevent the effects of big flares on Earth and its environment.

Intermittency in fluid turbulence can be emphasized through the analysis of Probability Distribution Functions (PDF) for velocity fluctuations, which display a strong non-gaussian behavior at small scales. Castaing et al. (1990) have introduced the idea that this behavior can be represented, in the framework of a multiplicative cascade model, by a convolution of gaussians whose variances is distributed according to a log-normal distribution. In this letter we have tried to test this conjecture on the MHD solar wind turbulence by performing a fit of the PDF of the bulk speed and magnetic field intensity fluctuations calculated in the solar wind, with the model. This fit allows us to calculate a parameter depending on the scale, which represents the width of the log-normal distribution of the variances of the gaussians. The physical implications of the obtained values of the parameter as well as of its scaling law are finally discussed