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

This White Paper outlines the importance of addressing the fundamental science theme <<How are charged particles energized in space plasmas>> through a future ESA mission. The White Paper presents five compelling science questions related to particle energization by shocks, reconnection,waves and turbulence, jets and their combinations. Answering these questions requires resolving scale coupling, nonlinearity and nonstationarity, which cannot be done with existing multi-point observations. In situ measurements from a multi-point, multi-scale L-class plasma observatory consisting of at least 7 spacecraft covering fluid, ion and electron scales are needed. The plasma observatory will enable a paradigm shift in our comprehension of particle energization and space plasma physics in general, with very important impact on solar and astrophysical plasmas. It will be the next logical step following Cluster, THEMIS and MMS for the very large and active European space plasmas community. Being one of the cornerstone missions of the future ESA Voyage 2035-2050 science program, it would further strengthen the European scientific and technical leadership in this important field.

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.

Fluctuations of solar wind magnetic field and plasma parameters exhibit a typical turbulence power spectrum with a spectral index ranging between $\sim -5/3$ and $\sim -3/2$. In particular, at $1$ AU, the magnetic field spectrum, observed within fast corotating streams, also shows a clear steepening for frequencies higher than the typical proton scales, of the order of $\sim 3\times10^{-1}$ Hz, and a flattening towards $1/f$ at frequencies lower than $\sim 10^{-3}$ Hz. However, the current literature reports observations of the low-frequency break only for fast streams. Slow streams, as observed to date, have not shown a clear break, and this has commonly been attributed to slow wind intervals not being long enough. Actually, because of the longer transit time from the Sun, slow wind turbulence would be older and the frequency break would be shifted to lower frequencies with respect to fast wind. Based on this hypothesis, we performed a careful search for long-lasting slow wind intervals throughout $12$ years of Wind satellite measurements. Our search, based on stringent requirements not only on wind speed but also on the level of magnetic compressibility and Alfv\'enicity of the turbulent fluctuations, yielded $48$ slow wind streams lasting longer than $7$ days. This result allowed us to extend our study to frequencies sufficiently low and, for the first time in the literature, we are able to show that the $1/f$ magnetic spectral scaling is also present in the slow solar wind, provided the interval is long enough. However, this is not the case for the slow wind velocity spectrum, which keeps the typical Kolmogorov scaling throughout the analysed frequency range. After ruling out the possible role of compressibility and Alfv\'enicity for the 1/f scaling, a possible explanation in terms of magnetic amplitude saturation, as recently proposed in the literature, is suggested.

The degree of deviation from the thermodynamic equilibrium in the ion velocity distribution functions (VDFs), measured by the Magnetospheric Multiscale (MMS) mission in the Earth's turbulent magnetosheath, is quantitatively investigated. Taking advantage of MMS ion data, having a resolution never reached before in space missions, and of the comparison with Vlasov-Maxwell simulations, this analysis aims at relating any deviation from Maxwellian equilibrium to typical plasma parameters. Correlations of the non-Maxwellian features with plasma quantities such as electric fields, ion temperature, current density and ion vorticity are very similar in both magnetosheath data and numerical experiments, and suggest that distortions in the ion VDFs occur close to (but not exactly at) peaks in current density and ion temperature. Similar results have also been found during a magnetopause crossing by MMS. This work could help clarifying the origin of distortion of the ion VDFs in space plasmas.

The properties of the turbulence which develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. Here we provide a direct comparison of our simulations with reported observations of reconnection events in the magnetotail investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show the development of a turbulent cascade consistent with spacecraft observations, statistics of the the dissipation mechanisms in the turbulent outflows similar to the one observed in reconnection jets in the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.

The collision of magnetic reconnection jets is studied by means of a three dimensional numerical simulation at kinetic scale, in the presence of a strong guide field. We show that turbulence develops due to the jets collision producing several current sheets in reconnection outflows, aligned with the guide field direction. The turbulence is mainly two-dimensional, with stronger gradients in the plane perpendicular to the guide field and a low wave-like activity in the parallel direction. First, we provide a numerical method to isolate the central turbulent region. Second, we analyze spatial second-order structure function and prove that turbulence is confined in this region. Finally, we compute local magnetic and electric frequency spectra, finding a trend in the sub-ion range that differs from typical cases for which the Taylor hypothesis is valid, as well as wave activity in the range between ion and electron cyclotron frequencies. Our results are relevant to understand observations of reconnection jets collisions in space plasmas.

In the context of space and astrophysical plasma turbulence and particle heating, several vocabularies emerge for estimating turbulent energy dissipation rate, including Kolmogorov-Yaglom third-order law and, in its various forms, $\boldsymbol{j}\cdot\boldsymbol{E}$ (work done by the electromagnetic field on particles), and $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ (pressure-strain interaction), to name a couple. It is now understood that these energy transfer channels, to some extent, are correlated with coherent structures. In particular, we find that different energy dissipation proxies, although not point-wise correlated, are concentrated in proximity to each other, for which they decorrelate in a few $d_i$(s). However, the energy dissipation proxies dominate at different scales. For example, there is an inertial range over which the third-order law is meaningful. Contributions from scale bands stemming from scale-dependent spatial filtering show that, the energy exchange through $\boldsymbol{j}\cdot\boldsymbol{E}$ mainly results from large scales, while the energy conversion from fluid flow to internal through $-\left( \boldsymbol{P} \cdot \nabla \right) \cdot \boldsymbol{u}$ dominates at small scales.

Reconnection outflows are regions of intense recent scrutiny, from in situ observations and from simulations. These regions are host to a variety of instabilities and intense energy exchanges, often even superior to the main reconnection site. We report here a number of results drawn from investigation of simulations. First, the outflows are observed to become unstable to drift instabilities. Second, these instabilities lead to the formation of secondary reconnection sites. Third, the secondary processes are responsible for large energy exchanges and particle energization. Finally, the particle distribution function are modified to become non-Maxwellian and include multiple interpenetrating populations.

The properties of inertial and kinetic range solar wind turbulence have been investigated with the arbitrary-order Hilbert spectral analysis method, applied to high-resolution density measurements. Due to the small sample size, and to the presence of strong non-stationary behavior and large-scale structures, the classical structure function analysis fails to detect power law behavior in the inertial range, and may underestimate the scaling exponents. However, the Hilbert spectral method provides an optimal estimation of the scaling exponents, which have been found to be close to those for velocity fluctuations in fully developed hydrodynamic turbulence. At smaller scales, below the proton gyroscale, the system loses its intermittent multiscaling properties, and converges to a monofractal process. The resulting scaling exponents, obtained at small scales, are in good agreement with those of classical fractional Brownian motion, indicating a long-term memory in the process, and the absence of correlations around the spectral break scale. These results provide important constraints on models of kinetic range turbulence in the solar wind.

The transfer of energy from large to small scales in solar wind turbulence is an important ingredient of the longstanding question about the mechanism of the interplanetary plasma heating. Previous studies have shown that magnetohydrodynamic (MHD) turbulence is statistically compatible with the observed solar wind heating as it expands in the heliosphere. However, in order to understand which processes contribute to the plasma heating, it is necessary to have a local description of the energy flux across scales. To this aim, it is customary to use indicators such as the magnetic field partial variance of increments (PVI), which is associated with the local, relative, scale-dependent magnetic energy. A more complete evaluation of the energy transfer should also include other terms, related to velocity and cross-helicity. This is achieved here by introducing a proxy for the local, scale dependent turbulent energy transfer rate epsilon_Dt(t), based on the third-order moment scaling law for MHD turbulence. Data from Helios 2 are used to determine the statistical properties of such a proxy in comparison with the magnetic and velocity fields PVI, and the correlation with local solar wind heating is computed. PVI and epsilon_Dt(t) are generally well correlated, however epsilon_Dt(t) is a very sensitive proxy that can exhibit large amplitude values, both positive and negative, even for low amplitude peaks in the PVI. Furthermore, epsilon_Dt(t) is very well correlated with local increases of temperature when large amplitude bursts of energy transfer are localized, thus suggesting an important role played by this proxy in the study of plasma energy dissipation.

Thanks to multi-spacecraft mission, it has recently been possible to directly estimate the current density in space plasmas, by using magnetic field time series from four satellites flying in a quasi perfect tetrahedron configuration. The technique developed, commonly called 'curlometer' permits a good estimation of the current density when the magnetic field time series vary linearly in space. This approximation is generally valid for small spacecraft separation. The recent space missions Cluster and Magnetospheric Multiscale (MMS) have provided high resolution measurements with inter-spacecraft separation up to 100 km and 10 km, respectively. The former scale corresponds to the proton gyroradius/ion skin depth in 'typical' solar wind conditions, while the latter to sub-proton scale. However, some works have highlighted an underestimation of the current density via the curlometer technique with respect to the current computed directly from the velocity distribution functions, measured at sub-proton scales resolution with MMS. In this paper we explore the limit of the curlometer technique studying synthetic data sets associated to a cluster of four artificial satellites allowed to fly in a static turbulent field, spanning a wide range of relative separation. This study tries to address the relative importance of measuring plasma moments at very high resolution from a single spacecraft with respect to the multi-spacecraft missions in the current density evaluation.

Solar Active Region NOAA 11158 has hosted a number of strong flares, including one X2.2 event. The complexity of current density and current helicity are studied through cancellation analysis of their sign-singular measure, which features power-law scaling. Spectral analysis is also performed, revealing the presence of two separate scaling ranges with different spectral index. The time evolution of parameters is discussed. Sudden changes of the cancellation exponents at the time of large flares, and the presence of correlation with EUV and X-ray flux, suggest that eruption of large flares can be linked to the small scale properties of the current structures.

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.

The intermittency of density fluctuations in the solar wind at kinetic scales has been examined using high time resolution Faraday cup measurements from the Spektr-R spacecraft. It was found that the probability density functions (PDFs) of the fluctuations are highly non-Gaussian over this range, but do not show large changes in shape with scale. These properties are statistically similar to those of the magnetic fluctuations and are important to understanding the dynamics of small scale turbulence in the solar wind. Possible explanations for the behavior of the density and magnetic fluctuations are discussed.

In situ observations of the solar wind frequently show the temperature of $\alpha$-particles (fully ionized helium), $T_\alpha$, to significantly differ from that of protons (ionized hydrogen), $T_p$. Many heating processes in the plasma act preferentially on $\alpha$-particles, even as collisions among ions act to gradually establish thermal equilibrium. Measurements from the $\textit{Wind}$ spacecraft's Faraday cups reveal that, at $r=1.0\ \textrm{AU}$ from the Sun, the observed values of the $\alpha$-proton temperature ratio, $\theta_{\alpha p} \equiv T_\alpha\,/\,T_p$ has a complex, bimodal distribution. This study applied a simple model for the radial evolution of $\theta_{\alpha p}$ to these data to compute expected values of $\theta_{\alpha p}$ at $r=0.1\ \textrm{AU}$. These inferred $\theta_{\alpha p}$-values have no trace of the bimodality seen in the $\theta_{\alpha p}$-values measured at $r=1.0\ \textrm{AU}$ but are instead consistent with the actions of the known mechanisms for $\alpha$-particle preferential heating. This result underscores the importance of collisional processes in the dynamics of the solar wind and suggests that similar mechanisms may lead to preferential $\alpha$-particle heating in both slow and fast wind.

Solar wind is probably the best laboratory to study turbulence in astrophysical plasmas. In addition to the presence of magnetic field, the differences with neutral fluid isotropic turbulence are: weakness of collisional dissipation and presence of several characteristic space and time scales. In this paper we discuss observational properties of solar wind turbulence in a large range from the MHD to the electron scales. At MHD scales, within the inertial range, turbulence cascade of magnetic fluctuations develops mostly in the plane perpendicular to the mean field. Solar wind turbulence is compressible in nature. The spectrum of velocity fluctuations do not follow magnetic field one. Probability distribution functions of different plasma parameters are not Gaussian, indicating presence of intermittency. At the moment there is no global model taking into account all these observed properties of the inertial range. At ion scales, turbulent spectra have a break, compressibility increases and the density fluctuation spectrum has a local flattening. Around ion scales, magnetic spectra are variable and ion instabilities occur as a function of the local plasma parameters. Between ion and electron scales, a small scale turbulent cascade seems to be established. Approaching electron scales, the fluctuations are no more self-similar: an exponential cut-off is usually observed indicating an onset of dissipation. The nature of the small scale cascade and a possible dissipation mechanism are still under debate.

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.

Interplanetary magnetic field magnitude fluctuations are notoriously more intermittent than velocity fluctuations in both fast and slow wind. This behaviour has been interpreted in terms of the anomalous scaling observed in passive scalars in fully developed hydrodynamic turbulence. In this paper, the strong intermittent nature of the interplanetary magnetic field is briefly discussed comparing results performed during different phases of the solar cycle. The scaling properties of the interplanetary magnetic field magnitude show solar cycle variation that can be distinguished in the scaling exponents revealed by structure functions. The scaling exponents observed around solar maximum coincide, within the errors, to those measured for passive scalars in hydrodynamic turbulence. However, it is also found that the values are not universal in the sense that the solar cycle variation may be reflected in dependence on the structure of the velocity field.

The asymmetric shape of reversals of the Earth's magnetic field indicates a possible connection with relaxation oscillations as they were early discussed by van der Pol. A simple mean-field dynamo model with a spherically symmetric $\alpha$ coefficient is analysed with view on this similarity, and a comparison of the time series and the phase space trajectories with those of paleomagnetic measurements is carried out. For highly supercritical dynamos a very good agreement with the data is achieved. Deviations of numerical reversal sequences from Poisson statistics are analysed and compared with paleomagnetic data. The role of the inner core is discussed in a spectral theoretical context and arguments and numerical evidence is compiled that the growth of the inner core might be important for the long term changes of the reversal rate and the occurrence of superchrons.

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.

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 interplanetary induced electric field e=vxb is studied, using solar wind time series. The probability distribution functions (PDFs) of the electric field components are measured from the data and their non-gaussianity is discussed. Moreover, for the first time we show that the electric field turbulence is characterized by intermittency. This point is addressed by studying, as usual, the scaling of the PDFs of field increments, which allows a quantitative characterization of intermittency.

We analyzed intermittency in the solar wind, as observed on the ecliptic plane, looking at magnetic field and velocity fluctuations between 0.3 and 1 AU, for both fast and slow wind and for compressive and directional fluctuations. Our analysis focused on the property that probability distribution functions of a fluctuating field affected by intermittency become more and more peaked at smaller and smaller scales. Since the peakedness of a distribution is measured by its flatness factor we studied the behavior of this parameter for different scales to estimate the degree of intermittency of our time series. We confirmed that both magnetic field and velocity fluctuations are rather intermittent and that compressive magnetic fluctuations are generally more intermittent than the corresponding velocity fluctuations. In addition, we observed that compressive fluctuations are always more intermittent than directional fluctuations and that while slow wind intermittency does not depend on the radial distance from the sun, fast wind intermittency of both magnetic field and velocity fluctuations clearly increases with the heliocentric distance....

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