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

The statistical properties of ions in two-dimensional fully developed turbulence have been compared between two different numerical algorithms. In particular, we compare Hybrid Particle In Cell (hybrid PIC with fluid electrons) and full PIC simulations, focusing on particle diffusion and acceleration phenomena. To investigate several heliospheric plasma conditions, a series of numerical simulations has been performed by varying the plasma $\beta$ - the ratio between kinetic and magnetic pressure. These numerical studies allow the exploration of different scenarios, going from the solar corona (low $\beta$) to the solar wind ($\beta \sim 1$), as well as the Earth's magnetosheath (high $\beta$). It has been found that the two approaches compare pretty well, especially for the spectral properties of the magnetic field and the ion diffusion statistics. Small differences among the models have been found regarding the electric field behaviour at sub-ion scales and the acceleration statistics, due evidently to the more consistent treatment of the plasma in the full PIC approach.

A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.

A novel technique is presented for describing and visualizing the local topology of the magnetic field using single-spacecraft data in the solar wind. The approach merges two established techniques: the Grad-Shafranov (GS) reconstruction method, which provides a plausible regional two-dimensional magnetic field surrounding the spacecraft trajectory, and the Partial Variance of Increments (PVI) technique that identifies coherent magnetic structures, such as current sheets. When applied to one month of Wind magnetic field data at 1-minute resolution, we find that the quasi-two-dimensional turbulence emerges as a sea of magnetic islands and current sheets. Statistical analysis confirms that current sheets associated with high values of PVI are mostly located between and within the GS magnetic islands, corresponding to X-points and internal boundaries. The method shows great promise for visualizing and analyzing single-spacecraft data from missions such as Parker Solar Probe and Solar Orbiter, as well as 1 AU Space Weather monitors such as ACE, Wind and IMAP.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.

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.

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.

Understanding the nature of the turbulent fluctuations below the ion gyroradius in solar-wind turbulence is a great challenge. Recent studies have been mostly in favor of kinetic Alfv\'en wave (KAW) type of fluctuations, but other kinds of fluctuations with characteristics typical of magnetosonic, whistler and ion Bernstein modes, could also play a role depending on the plasma parameters. Here we investigate the properties of the sub-proton-scale cascade with high-resolution hybrid-kinetic simulations of freely-decaying turbulence in 3D3V phase space, including electron inertia effects. Two proton plasma beta are explored: the "intermediate" $\beta_p=1$ and "low" $\beta_p=0.2$ regimes, both typically observed in solar wind and corona. The magnetic energy spectum exhibits $k_\perp^{-8/3}$ and $k_\|^{-7/2}$ power laws at $\beta_p=1$, while they are slightly steeper at $\beta_p=0.2$. Nevertheless, both regimes develop a spectral anisotropy consistent with $k_\|\sim k_\perp^{2/3}$ at $k_\perp\rho_p>1$, and pronounced small-scale intermittency. In this context, we find that the kinetic-scale cascade is dominated by KAW-like fluctuations at $\beta_p=1$, whereas the low-$\beta$ case presents a more complex scenario suggesting the simultaneous presence of different types of fluctuations. In both regimes, however, a non-negligible role of ion Bernstein type of fluctuations at the smallest scales seems to emerge.

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.

The magnetic topology and field line random walk properties of a nanoflare-heated and magnetically confined corona are investigated in the reduced magnetohydrodynamic regime. Field lines originating from current sheets form coherent structures, called Current Sheet Connected (CSC) regions, extended around them. CSC field line random walk is strongly anisotropic, with preferential diffusion along the current sheets' in-plane length. CSC field line random walk properties remain similar to those of the entire ensemble but exhibit enhanced mean square displacements and separations due to the stronger magnetic field intensities in CSC regions. The implications for particle acceleration and heat transport in the solar corona and wind, and for solar moss formation are discussed.

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.

One view of the solar-wind turbulence is that the observed highly anisotropic fluctuations at spatial scales near the proton inertial length $d_p$ may be considered as Kinetic Alfv\'en waves (KAWs). In the present paper, we show how phase-mixing of large-scale parallel propagating Alfv\'en waves is an efficient mechanism for the production of KAWs at wavelengths close to $d_p$ and at large propagation angle with respect to the magnetic field. Magnetohydrodynamic (MHD), Hall-Magnetohydrodynamic (HMHD), and hybrid Vlasov-Maxwell (HVM) simulations modeling the propagation of Alfv\'en waves in inhomogeneous plasmas are performed. In linear regime, the role of dispersive effects is singled out by comparing MHD and HMHD results. Fluctuations produced by phase-mixing are identified as KAWs through a comparison of polarization of magnetic fluctuations and wave group velocity with analytical linear predictions. In the nonlinear regime, comparison of HMHD and HVM simulations allows to point out the role of kinetic effects in shaping the proton distribution function. We observe generation of temperature anisotropy with respect to the local magnetic field and production of field-aligned beams. The regions where the proton distribution function highly departs from thermal equilibrium are located inside the shear layers, where the KAWs are excited, this suggesting that the distortions of the proton distribution are driven by a resonant interaction of protons with KAW fluctuations. Our results are relevant in configurations where magnetic field inhomogeneities are present, as, for example, in the solar corona where the presence of Alfv\'en waves has been ascertained.

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.

Employing a simple ideal magnetohydrodynamic model in spherical geometry,we show that the presence of either rotation or finite magnetic helicity is sufficient to induce dynamical reversals of the magnetic dipole moment. The statistical character of the model is similar to that of terrestrial magnetic field reversals, with the similarity being stronger when rotation is present.The connection between long time correlations, $1/f$ noise, and statistics of reversals is supported, consistent with earlier suggestions.

Based on global conservation principles, magnetohydrodynamic (MHD) relaxation theory predicts the existence of several equilibria, such as the Taylor state or global dynamic alignment. These states are generally viewed as very long-time and large-scale equilibria, which emerge only after the termination of the turbulent cascade. As suggested by hydrodynamics and by recent MHD numerical simulations, relaxation processes can occur during the turbulent cascade that will manifest themselves as local patches of equilibrium-like configurations. Using multi-spacecraft analysis techniques in conjunction with Cluster data, we compute the current density and flow vorticity and for the first time demonstrate that these localized relaxation events are observed in the solar wind. Such events have important consequences for the statistics of plasma turbulence.

Magnetic reconnection at the interface between coronal holes and loops, so-called interchange reconnection, can release the hotter, denser plasma from magnetically confined regions into the heliosphere, contributing to the formation of the highly variable slow solar wind. The interchange process is often thought to develop at the apex of streamers or pseudo-streamers, near Y and X-type neutral points, but slow streams with loop composition have been recently observed along fanlike open field lines adjacent to closed regions, far from the apex. However, coronal heating models, with magnetic field lines shuffled by convective motions, show that reconnection can occur continuously in unipolar magnetic field regions with no neutral points: photospheric motions induce a magnetohydrodynamic turbulent cascade in the coronal field that creates the necessary small scales, where a sheared magnetic field component orthogonal to the strong axial field is created locally and can reconnect. We propose that a similar mechanism operates near and around boundaries between open and closed regions inducing a continual stochastic rearrangement of connectivity. We examine a reduced magnetohydrodynamic model of a simplified interface region between open and closed corona threaded by a strong unipolar magnetic field. This boundary is not stationary, becomes fractal, and field lines change connectivity continuously, becoming alternatively open and closed. This model suggests that slow wind may originate everywhere along loop-coronal hole boundary regions, and can account naturally and simply for outflows at and adjacent to such boundaries and for the observed diffusion of slow wind around the heliospheric current sheet.

Correlation anisotropy emerges dynamically in magnetohydrodynamics (MHD), producing stronger gradients across the large-scale mean magnetic field than along it. This occurs both globally and locally, and has significant implications in space and astrophysical plasmas, including particle scattering and transport, and theories of turbulence. Properties of local correlation anisotropy are further documented here by showing through numerical experiments that the effect is intensified in more localized estimates of the mean field. The mathematical formulation of this property shows that local anisotropy mixes second-order with higher order correlations. Sensitivity of local statistical estimates to higher order correlations can be understood in connection with the stochastic coordinate system inherent in such formulations. We demonstrate this in specific cases, illustrate the connection to higher order statistics by showing the sensitivity of local anisotropy to phase randomization, and thus establish that the local structure function is not a measure of the energy spectrum. Evidently the local enhancement of correlation anisotropy is of substantial fundamental interest, and this phenomenon must be understood in terms of higher order correlations, fourth-order and above.

This study shows the connection between three previously observed but seemingly unrelated phenomena in hydrodynamic (HD) and magnetohydrodynamic (MHD) turbulent flows, involving the emergence of fluctuations occurring on very long time scales: the low-frequency 1/f noise in the power frequency spectrum, the delayed ergodicity of complex valued amplitude fluctuations in wavenumber space, and the spontaneous flippings or reversals of large scale fields. Direct numerical simulations of ideal MHD and HD are employed in three space dimensions, at low resolution, for long periods of time, and with high accuracy to study several cases: Different geometries, presence of rotation and/or a uniform magnetic field, and different values of the associated conserved global quantities. It is conjectured that the origin of all these long-time phenomena is rooted in the interaction of the longest wavelength fluctuations available to the system with fluctuations at much smaller scales. The strength of this non-local interaction is controlled either by the existence of conserved global quantities with a back-transfer in Fourier space, or by the presence of a slow manifold in the dynamics.