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.

The shape of the electron velocity distribution function plays an important role in the dynamics of the solar wind acceleration. Electrons are normally modelled with three components, the core, the halo, and the strahl. We investigate how well the fast strahl electrons in the inner heliosphere preserve the information about the coronal electron temperature at their origin. We analysed the data obtained by two missions, Helios spanning the distances between 65 and 215 R$_S$, and Parker Solar Probe (PSP) reaching down to 35 R$_S$ during its first two orbits around the Sun. The electron strahl was characterised with two parameters, pitch-angle width (PAW), and the strahl parallel temperature (T$_{s\parallel}$). PSP observations confirm the already reported dependence of strahl PAW on core parallel plasma beta ($\beta_{ec\parallel}$)\citep{Bercic2019}. Most of the strahl measured by PSP appear narrow with PAW reaching down to 30$^o$. The portion of the strahl velocity distribution function aligned with the magnetic field is for the measured energy range well described by a Maxwellian distribution function. T$_{s\parallel}$ was found to be anti-correlated with the solar wind velocity, and independent of radial distance. These observations imply that T$_{s\parallel}$ carries the information about the coronal electron temperature. The obtained values are in agreement with coronal temperatures measured using spectroscopy (David et al. 2998), and the inferred solar wind source regions during the first orbit of PSP agree with the predictions using a PFSS model (Bale et al. 2019, Badman et al. 2019).

Although the origins of slow solar wind are unclear, there is increasing evidence that at least some of it is released in a steady state on over-expanded coronal hole magnetic field lines. This type of slow wind has similar properties to the fast solar wind, including a high degree of Alfv\'enicity. In this study a combination of proton, alpha particle, and electron measurements are used to investigate the kinetic properties of a single interval of slow Alfv\'enic wind at 0.35 AU. It is shown that this slow Alfv\'enic interval is characterised by high alpha particle abundances, pronounced alpha-proton differential streaming, strong proton beams, and large alpha to proton temperature ratios. These are all features observed consistently in the fast solar wind, adding evidence that at least some Alfv\'enic slow solar wind also originates in coronal holes. Observed differences between speed, mass flux, and electron temperature between slow Alfv\'enic and fast winds are explained by differing magnetic field geometry in the lower corona.

Parker Solar Probe (PSP) observations during its first encounter at 35.7 $R_\odot$ have shown the presence of magnetic field lines which are strongly perturbed to the point that they produce local inversions of the radial magnetic field, known as switchbacks. Their counterparts in the solar wind velocity field are local enhancements in the radial speed, or jets, displaying (in all components) the velocity-magnetic field correlation typical of large amplitude Alfv\'en waves propagating away from the Sun. Switchbacks and radial jets have previously been observed over a wide range of heliocentric distances by Helios, WIND and Ulysses, although they were prevalent in significantly faster streams than seen at PSP. Here we study via numerical MHD simulations the evolution of such large amplitude Alfv\'enic fluctuations by including, in agreement with observations, both a radial magnetic field inversion and an initially constant total magnetic pressure. Despite the extremely large excursion of magnetic and velocity fields, switchbacks are seen to persist for up to hundreds of Alfv\'en crossing times before eventually decaying due to the parametric decay instability. Our results suggest that such switchback/jet configurations might indeed originate in the lower corona and survive out to PSP distances, provided the background solar wind is sufficiently calm, in the sense of not being pervaded by strong density fluctuations or other gradients, such as stream or magnetic field shears, that might destabilize or destroy them over shorter timescales.

Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfv\'enic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and electron betas, turbulence strength, and injection scale) exhibits a quantitative agreement on the spectral, intermittency, and cascade-rate properties with in situ observations, despite the different driving mechanisms. Such agreement demonstrates a certain universality of the turbulent cascade from magnetohydrodynamic (MHD) to sub-ion scales, whose properties are mainly determined by the selected parameters, also indicating that the KH instability-driven turbulence has a quasi-2D nature. The validity of the Taylor hypothesis in the sub-ion spatial range suggests that the fluctuations at sub-ion scales have predominantly low frequencies, consistent with a kinetic Alfv\'en wave-like nature or with quasi-static structures. Finally, the third-order structure function analysis indicates that the cascade rate of the turbulence generated by a KH event in the magnetopause is an order of magnitude larger than in the ambient magnetosheath.

The relationship between a decaying plasma turbulence and proton fire hose instabilities in a slowly expanding plasma is investigated using three-dimensional (3-D) hybrid expanding box simulations. We impose an initial ambient magnetic field along the radial direction, and we start with an isotropic spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations with zero cross-helicity. A turbulent cascade rapidly develops and leads to a weak proton heating that is not sufficient to overcome the expansion-driven perpendicular cooling. The plasma system eventually drives the parallel and oblique fire hose instabilities that generate quasi-monochromatic wave packets that reduce the proton temperature anisotropy. The fire hose wave activity has a low amplitude with wave vectors quasi-parallel/oblique with respect to the ambient magnetic field outside of the region dominated by the turbulent cascade and is discernible in one-dimensional power spectra taken only in the direction quasi-parallel/oblique with respect to the ambient magnetic field; at quasi-perpendicular angles the wave activity is hidden by the turbulent background. These waves are partly reabsorbed by protons and partly couple to and participate in the turbulent cascade. Their presence reduces kurtosis, a measure of intermittency, and the Shannon entropy but increases the Jensen-Shannon complexity of magnetic fluctuations; these changes are weak and anisotropic with respect to the ambient magnetic field and it's not clear if they can be used to indirectly discern the presence of instability-driven waves.

Electron velocity distribution functions in the solar wind according to standard models consist of 4 components, of which 3 are symmetric - the core, the halo, and the superhalo, and one is magnetic field-aligned, beam-like population, referred to as the strahl. We analysed in-situ measurements provided by the two Helios spacecrafts to study the behaviour of the last, the strahl electron population, in the inner Solar system between 0.3 and 1 au. The strahl is characterised with a pitch-angle width (PAW) depending on electron energy and evolving with radial distance. We find different behaviour of the strahl electrons for solar wind separated into types by the core electron beta parallel value ($\beta_{ec\parallel}$). For the low-$\beta_{ec\parallel}$ solar wind the strahl component is more pronounced, and the variation of PAW is electron energy dependent. At low energies a slight focusing over distance is observed, and the strahl PAW measured at 0.34 au agrees with the width predicted by a collisionless focusing model. The broadening observed for higher-energy strahl electrons during expansion can be described by an exponential relation, which points toward an energy dependent scattering mechanism. In the high-$\beta_{ec\parallel}$ solar wind the strahl appears broader in consistence with the high-$\beta_{ec\parallel}$ plasma being more unstable with respect to kinetic instabilities. Finally we extrapolate our observations to the distance of 0.16 au, predicting the strahl PAWs in the low-$\beta_{ec\parallel}$ solar wind to be $\sim$ 29$^o$ for all energies, and in the high-$\beta_{ec\parallel}$ solar wind a bit broader, ranging between 37$^o$ and 65$^o$.

By means of three dimensional high-resolution hybrid simulations we study the properties of the magnetic field spectral anisotropies near and beyond ion kinetic scales. By using both a Fourier analysis and a local analysis based on multi-point 2nd-order structure function techniques, we show that the anisotropy observed is less than what expected by standard wave normal modes turbulence theories although the non linear energy transfer is still in the perpendicular direction, only advected in the parallel direction as expected balancing the non-linear energy transfer time and the decorrelation time. Such result can be explained by a phenomenological model based on the formation of strong intermittent two-dimensional structures in the plane perpendicular to the local mean field that have some prescribed aspect ratio eventually depending on the scale. This model support the idea that small scales structures, such as reconnecting current sheets, contribute significantly to the formation of the turbulent cascade at kinetic scales.

We analyse velocity fluctuations in the solar wind at magneto-fluid scales in two datasets, extracted from Wind data in the period 2005-2015, that are characterised by strong or weak expansion. Expansion affects measurements of anisotropy because it breaks axisymmetry around the mean magnetic field. Indeed, the small-scale three-dimensional local anisotropy of magnetic fluctuations ({\delta}B) as measured by structure functions (SF_B) is consistent with tube-like structures for strong expansion. When passing to weak expansion, structures become ribbon-like because of the flattening of SFB along one of the two perpendicular directions. The power-law index that is consistent with a spectral slope -5/3 for strong expansion now becomes closer to -3/2. This index is also characteristic of velocity fluctuations in the solar wind. We study velocity fluctuations ({\delta}V) to understand if the anisotropy of their structure functions (SF_V ) also changes with the strength of expansion and if the difference with the magnetic spectral index is washed out once anisotropy is accounted for. We find that SF_V is generally flatter than SF_B. When expansion passes from strong to weak, a further flattening of the perpendicular SF_V occurs and the small-scale anisotropy switches from tube-like to ribbon-like structures. These two types of anisotropy, common to SF_V and SF_B, are associated to distinct large-scale variance anisotropies of {\delta}B in the strong- and weak-expansion datasets. We conclude that SF_V shows anisotropic three-dimensional scaling similar to SF_B, with however systematic flatter scalings, reflecting the difference between global spectral slopes.

Plasma processes occurring in the corona and solar wind can be probed by studying the thermodynamic properties of different ion species. However, most in-situ observations of positive ions in the solar wind are taken at 1 AU, where information on their solar source properties may have been irreversibly erased. In this study we aimed to use the properties of alpha particles at heliocentric distances between 0.3 and 1 AU to study plasma processes occurring at the points of observation, and to infer processes occurring inside 0.3 AU by comparing our results to previous remote sensing observations of the plasma closer to the Sun. We reprocessed the original Helios positive ion distribution functions, isolated the alpha particle population, and computed the alpha particle number density, velocity, and magnetic field perpendicular and parallel temperatures. We then investigated the radial variation of alpha particle temperatures in fast solar wind observed between 0.3 and 1 AU. Between 0.3 and 1 AU alpha particles are heated in the magnetic field perpendicular direction, and cooled in the magnetic field parallel direction. Alpha particle evolution is bounded by the alpha firehose instability threshold, which provides one possible mechanism to explain the observed parallel cooling and perpendicular heating. Closer to the Sun our observations suggest that the alpha particles undergo heating the perpendicular direction, whilst the large magnetic field parallel temperatures observed at 0.3 AU may be due the combined effect of double adiabatic expansion and alpha particle deceleration inside 0.3 AU.

We discuss properties of Alfv\'enic fluctuations with large amplitude in plasmas characterised by low magnetic field compression. We note that in such systems power laws can not develop with arbitrarily steep slopes at large scales, i.e. when $|\delta \bf{B}|$ becomes of the order of the background field $|\bf{B}|$. In such systems there is a scale $l_0$ at which the spectrum has to break due to the condition of weak compressibility. A very good example of this dynamics is offered by solar wind fluctuations in Alfv\'enic fast streams, characterised by the property of constant field magnitude. We show here that the distribution of $\delta B=|\delta \bf{B}|$ in the fast wind displays a strong cut-off at $\delta B/|{\bf B}|\lesssim2$, as expected for fluctuations bounded on a sphere of radius $B=|{\bf B}|$. This is also associated with a saturation of the rms of the fluctuations at large scales and introduces a specific length $l_0$ above which the amplitude of the fluctuations becomes independent on the scale $l$. Consistent with that, the power spectrum at $l>l_0$ is characterised by a -1 spectral slope, as expected for fluctuations that are scale-independent. Moreover, we show that the spectral break between the 1/f and inertial range in solar wind spectra indeed corresponds to the scale $l_0$ at which $\left<\delta B/B\right>\sim1$. Such a simple model provides a possible alternative explanation of magnetic spectra observed in interplanetary space, also pointing out the inconsistency for a plasma to simultaneously maintain $|\bf{B}|\sim$const. at arbitrarily large scales and satisfy a Kolmogorov scaling.

Alfv\'enic fluctuations are very common features in the solar wind and are found especially within the main portion of fast wind streams while the slow wind usually is less Alfv\'enic and more variable. In general, fast and slow wind show many differences which span from the large scale structure to small scale phenomena including also a different turbulent behaviour. Recent studies, however, have shown that even slow wind can be sometimes highly Alfv\'enic with fluctuations as large as those of the fast wind. The present study is devoted to present many facets of this Alfv\'enic slow solar wind including for example the study of the source regions and their connection to coronal structures, large-scale properties and micro-scale phenomena and also impact on the spectral features. This study will be conducted performing a comparative analysis with the typical slow wind and with the fast wind. It has been found that the fast wind and the Alfv\'enic slow wind share common characteristics, probably attributable to their similar solar origin, i.e. coronal-hole solar wind. Given these similarities, it is suggested that in the Alfv\'enic slow wind a major role is played by the super-radial expansion responsible for the lower velocity. Relevant implications of these new findings for the upcoming Solar Orbiter and Solar Probe Plus missions, and more in general for turbulence measurements close to the Sun, will be discussed.

Robustly identifying the solar sources of individual packets of solar wind measured in interplanetary space remains an open problem. We set out to see if this problem is easier to tackle using solar wind measurements closer to the Sun than 1 AU, where the mixing and dynamical interaction of different solar wind streams is reduced. Using measurements from the Helios mission, we examined how the proton core temperature anisotropy and cross helicity varied with distance. At 0.3 AU there are two clearly separated anisotropic and isotropic populations of solar wind, that are not distinguishable at 1 AU. The anisotropic population is always Alfv\'enic and spans a wide range of speeds. In contrast the isotropic population has slow speeds, and contains a mix of Alfv\'enic wind with constant mass fluxes, and non-Alfv\'enic wind with large and highly varying mass fluxes. We split the in-situ measurements into three categories according these observations, and suggest that these categories correspond to wind that originated in the core of coronal holes, in or near active regions or the edges of coronal holes, and as small transients form streamers or pseudostreamers. Although our method by itself is simplistic, it provides a new tool that can be used in combination with other methods for identifying the sources of solar wind measured by Parker Solar Probe and Solar Orbiter.

In the near future, Parker Solar Probe and Solar Orbiter will provide the first comprehensive in-situ measurements of the solar wind in the inner heliosphere since the Helios mission in the 1970s. We describe a reprocessing of the original Helios ion distribution functions to provide reliable and reproducible data to characterise the proton core population of the solar wind in the inner heliosphere. A systematic fitting of bi-Maxwellian distribution functions was performed to the raw Helios ion distribution function data to extract the proton core number density, velocity, and temperatures parallel and perpendicular to the magnetic field. We present radial trends of these derived proton parameters, forming a benchmark from which new measurements in the inner heliosphere will be compared to. The new dataset has been made openly available for other researchers to use, along with the source code used to generate it.

A dynamical vectorial equation for homogeneous incompressible Hall-MHD turbulence together with the exact scaling law for third-order correlation tensors, analogous to that for the incompressible MHD, is rederived and applied to the results of two-dimensional hybrid simulations of plasma turbulence. At large (MHD) scales the simulations exhibits a clear inertial range where the MHD dynamic law is valid. In the sub-ion range the cascade continues via the Hall term but the dynamic law derived in the framework of incompressible Hall MHD equations is obtained only in a low plasma beta simulation. For a higher beta plasma the cascade rate decreases in the sub-ion range and the change becomes more pronounced as the plasma beta increases. This break in the cascade flux can be ascribed to non thermal (kinetic) features or to others terms in the dynamical equation that are not included in the Hall-MHD incompressible approximation.

We investigate the spectral properties of plasma turbulence from fluid to sub-ion scales by means of high-resolution three-dimensional (3D) numerical simulations performed with the hybrid particle-in-cell (HPIC) code CAMELIA. We produce extended turbulent spectra with well-defined power laws for the magnetic, ion bulk velocity, density, and electric fluctuations. The present results are in good agreement with previous two-dimensional (2D) HPIC simulations, especially in the kinetic range of scales, and reproduce several features observed in solar wind spectra. By providing scaling tests on many different architectures and convergence studies, we prove CAMELIA to represent a very efficient, accurate and reliable tool for investigating the develpoment of the turbulent cascade in the solar wind, being able to cover simultaneously several decades in wavenumber, also in 3D.

Spectral properties of the turbulent cascade from fluid to kinetic scales in collisionless plasmas are investigated by means of large-size three-dimensional (3D) hybrid (fluid electrons, kinetic protons) particle-in-cell simulations. Initially isotropic Alfv\`enic fluctuations rapidly develop a strongly anisotropic turbulent cascade, mainly in the direction perpendicular to the ambient magnetic field. The omnidirectional magnetic field spectrum shows a double power-law behavior over almost two decades in wavenumber, with a Kolmogorov-like index at large scales, a spectral break around ion scales, and a steepening at sub-ion scales. Power laws are also observed in the spectra of the ion bulk velocity, density, and electric field, both at magnetohydrodynamic (MHD) and at kinetic scales. Despite the complex structure, the omnidirectional spectra of all fields at ion and sub-ion scales are in remarkable quantitative agreement with those of a two-dimensional (2D) simulation with similar physical parameters. This provides a partial, a-posteriori validation of the 2D approximation at kinetic scales. Conversely, at MHD scales, the spectra of the density and of the velocity (and, consequently, of the electric field) exhibit differences between the 2D and 3D cases. Although they can be partly ascribed to the lower spatial resolution, the main reason is likely the larger importance of compressible effects in a full geometry. Our findings are also in remarkable quantitative agreement with solar wind observations.

A new path for the generation of a sub-ion scale cascade in collisionless plasma turbulence, triggered by magnetic reconnection, is uncovered by means of high-resolution two-dimensional hybrid-kinetic simulations employing two complementary approaches, Lagrangian and Eulerian, and different driving mechanisms. The simulation results provide clear numerical evidences that the development of power-law energy spectra below the so-called ion break occurs as soon as the first magnetic reconnection events take place, regardless of the actual state of the turbulent cascade at MHD scales. In both simulations, the reconnection-mediated small-scale energy spectrum of parallel magnetic fluctuations exhibits a very stable spectral slope of ~-2.8, whether or not a large-scale turbulent cascade has already fully developed. Once a quasi-stationary turbulent state is achieved, the spectrum of the total magnetic fluctuations settles towards a spectral index of -5/3 in the MHD range and of ~-3 at sub-ion scales.

We review the presence and signatures of the non-equilibrium processes, both non-Maxwellian distributions and non-equilibrium ionization, in the solar transition region, corona, solar wind, and flares. Basic properties of the non-Maxwellian distributions are described together with their influence on the heat flux as well as on the rates of individual collisional processes and the resulting optically thin synthetic spectra. Constraints on the presence of high-energy electrons from observations are reviewed, including positive detection of non-Maxwellian distributions in the solar corona, transition region, flares, and wind. Occurrence of non-equilibrium ionization is reviewed as well, especially in connection to hydrodynamic and generalized collisional-radiative modelling. Predicted spectroscopic signatures of non-equilibrium ionization depending on the assumed plasma conditions are summarized. Finally, we discuss the future remote-sensing instrumentation that can be used for detection of these non-equilibrium phenomena in various spectral ranges.

The relationship between a decaying strong turbulence and the mirror instability in a slowly expanding plasma is investigated using two-dimensional hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfvenic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at sub-ion scales. The imposed expansion (taking a strictly transverse ambient magnetic field) leads to generation of an important perpendicular proton temperature anisotropy that eventually drives the mirror instability. This instability generates large-amplitude, nonpropagating, compressible, pressure-balanced magnetic structures in a form of magnetic enhancements/humps that reduce the perpendicular temperature anisotropy.

We investigate properties of the ion-scale spectral break of solar wind turbulence by means of two-dimensional high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field perpendicular to the simulation box and add a spectrum of in-plane, large-scale, magnetic and kinetic fluctuations. We perform a set of simulations with different values of the plasma beta, distributed over three orders of magnitude, from 0.01 to 10. In all the cases, once turbulence is fully developed, we observe a power-law spectrum of the fluctuating magnetic field on large scales (in the inertial range) with a spectral index close to -5/3, while in the sub-ion range we observe another power-law spectrum with a spectral index systematically varying with $\beta$ (from around -3.6 for small values to around -2.9 for large ones). The two ranges are separated by a spectral break around ion scales. The length scale at which this transition occurs is found to be proportional to the ion inertial length, $d_i$, for $\beta \ll 1$ and to the ion gyroradius, $\rho_i = d_i\sqrt{\beta}$, for $\beta \gg 1$, i.e., to the larger between the two scales in both the extreme regimes. For intermediate cases, i.e., $\beta \sim 1$, a combination of the two scales is involved. We infer an empiric relation for the dependency of the spectral break on $\beta$ that provides a good fit over the whole range of values. We compare our results with in situ observations in the solar wind and suggest possible explanations for such a behavior.

The firehose and mirror instabilities are thought to arise in a variety of space and astrophysical plasmas, constraining the pressure anisotropies and drifts between particle species. The plasma stability depends on all species simultaneously, meaning that a combined analysis is required. Here, we present the first such analysis in the solar wind, using the long-wavelength stability parameters to combine the anisotropies and drifts of all major species (core and beam protons, alphas, and electrons). At the threshold, the firehose parameter was found to be dominated by protons (67%), but also to have significant contributions from electrons (18%) and alphas (15%). Drifts were also found to be important, contributing 57% in the presence of a proton beam. A similar situation was found for the mirror, with contributions of 61%, 28%, and 11% for protons, electrons, and alphas, respectively. The parallel electric field contribution, however, was found to be small at 9%. Overall, the long-wavelength thresholds constrain the data well (<1% unstable), and the implications of this are discussed.

Proton temperature anisotropies between the directions parallel and perpendicular to the mean magnetic field are usually observed in the solar wind plasma. Here, we employ a high-resolution hybrid particle-in-cell simulation in order to investigate the relation between spatial properties of the proton temperature and the peaks in the current density and in the flow vorticity. Our results indicate that, although regions where the proton temperature is enhanced and temperature anisotropies are larger correspond approximately to regions where many thin current sheets form, no firm quantitative evidence supports the idea of a direct causality between the two phenomena. On the other hand, quite a clear correlation between the behavior of the proton temperature and the out-of-plane vorticity is obtained.

In the outer heliosphere, beyond the solar wind termination shock, it is expected that the warped heliospheric current sheet forms a region of closely-packed, multiple, thin current sheets. Such a system may be subject to the ion-kinetic tearing instability, and hence generate magnetic islands and hot populations of ions associated with magnetic reconnection. Reconnection processes in this environment have important implications for local particle transport, and for particle acceleration at reconnection sites and in turbulence. We study this complex environment by means of three-dimensional hybrid simulations over long time scales, in order to capture the evolution from linear growth of the tearing instability to a fully developed turbulent state at late times. The final state develops from the highly ordered initial state via both forward and inverse cascades. Component and spectral anisotropy in the magnetic fluctuations is present when a guide field is included. The inclusion of a population of new-born interstellar pickup protons does not strongly affect these results. Finally, we conclude that reconnection between multiple current sheets can act as an important source of turbulence in the outer heliosphere, with implications for energetic particle acceleration and propagation.

The relationship between a decaying strong turbulence and kinetic instabilities in a slowly expanding plasma is investigated using two-dimensional (2-D) hybrid expanding box simulations. We impose an initial ambient magnetic field perpendicular to the simulation box, and we start with a spectrum of large-scale, linearly-polarized, random-phase Alfv\'enic fluctuations which have energy equipartition between kinetic and magnetic fluctuations and vanishing correlation between the two fields. A turbulent cascade rapidly develops, magnetic field fluctuations exhibit a Kolmogorov-like power-law spectrum at large scales and a steeper spectrum at ion scales. The turbulent cascade leads to an overall anisotropic proton heating, protons are heated in the perpendicular direction, and, initially, also in the parallel direction. The imposed expansion leads to generation of a large parallel proton temperature anisotropy which is at later stages partly reduced by turbulence. The turbulent heating is not sufficient to overcome the expansion-driven perpendicular cooling and the system eventually drives the oblique firehose instability in a form of localized nonlinear wave packets which efficiently reduce the parallel temperature anisotropy. This work demonstrates that kinetic instabilities may coexist with strong plasma turbulence even in a constrained 2-D regime.

The solar wind magnetic field contains rotations at a broad range of scales, which have been extensively studied in the MHD range. Here we present an extension of this analysis to the range between ion and electron kinetic scales. The distribution of rotation angles was found to be approximately log-normal, shifting to smaller angles at smaller scales almost self-similarly, but with small, statistically significant changes of shape. The fraction of energy in fluctuations with angles larger than $\alpha$ was found to drop approximately exponentially with $\alpha$, with e-folding angle $9.8^\circ$ at ion scales and $0.66^\circ$ at electron scales, showing that large angles ($\alpha > 30^\circ$) do not contain a significant amount of energy at kinetic scales. Implications for kinetic turbulence theory and the dissipation of solar wind turbulence are discussed.

We investigate properties of a solar wind-like plasma including a secondary alpha particle population exhibiting a parallel temperature anisotropy with respect to the background magnetic field, using linear and quasi-linear predictions and by means of one-dimensional hybrid simulations. We show that anisotropic alpha particles can drive a parallel fire hose instability analogous to that generated by protons, but that, remarkably, the instability can be triggered also when the parallel plasma beta of alpha particles is below unity. The wave activity generated by the alpha anisotropy affects the evolution of the more abundant protons, leading to their anisotropic heating. When both ion species have sufficient parallel anisotropies both of them can drive the instability, and we observe generation of two distinct peaks in the spectra of the fluctuations, with longer wavelengths associated to alphas and shorter ones to protons. If a non-zero relative drift is present, the unstable modes propagate preferentially in the direction of the drift associated with the unstable species. The generated waves scatter particles and reduce their temperature anisotropy to marginally stable state, and, moreover, they significantly reduce the relative drift between the two ion populations. The coexistence of modes excited by both species leads to saturation of the plasma in distinct regions of the beta/anisotropy parameter space for protons and alpha particles, in good agreement with in situ solar wind observations. Our results confirm that fire hose instabilities are likely at work in the solar wind and limit the anisotropy of different ion species in the plasma.

We investigate properties of plasma turbulence from magneto-hydrodynamic (MHD) to sub-ion scales by means of two-dimensional, high-resolution hybrid particle-in-cell simulations. We impose an initial ambient magnetic field, perpendicular to the simulation box, and we add a spectrum of large-scale magnetic and kinetic fluctuations, with energy equipartition and vanishing correlation. Once the turbulence is fully developed, we observe a MHD inertial range, where the spectra of the perpendicular magnetic field and the perpendicular proton bulk velocity fluctuations exhibit power-law scaling with spectral indices of -5/3 and -3/2, respectively. This behavior is extended over a full decade in wavevectors and is very stable in time. A transition is observed around proton scales. At sub-ion scales, both spectra steepen, with the former still following a power law with a spectral index of ~-3. A -2.8 slope is observed in the density and parallel magnetic fluctuations, highlighting the presence of compressive effects at kinetic scales. The spectrum of the perpendicular electric fluctuations follows that of the proton bulk velocity at MHD scales, and flattens at small scales. All these features, which we carefully tested against variations of many parameters, are in good agreement with solar wind observations. The turbulent cascade leads to on overall proton energization with similar heating rates in the parallel and perpendicular directions. While the parallel proton heating is found to be independent on the resistivity, the number of particles per cell and the resolution employed, the perpendicular proton temperature strongly depends on these parameters.

We present results from a high-resolution and large-scale hybrid (fluid electrons and particle-in-cell protons) two-dimensional numerical simulation of decaying turbulence. Two distinct spectral regions (separated by a smooth break at proton scales) develop with clear power-law scaling, each one occupying about a decade in wave numbers. The simulation results exhibit simultaneously several properties of the observed solar wind fluctuations: spectral indices of the magnetic, kinetic, and residual energy spectra in the magneto-hydrodynamic (MHD) inertial range along with a flattening of the electric field spectrum, an increase in magnetic compressibility, and a strong coupling of the cascade with the density and the parallel component of the magnetic fluctuations at sub-proton scales. Our findings support the interpretation that in the solar wind large-scale MHD fluctuations naturally evolve beyond proton scales into a turbulent regime that is governed by the generalized Ohm's law.

We investigate properties of the plasma fluid motion in the large amplitude low frequency fluctuations of highly Alfv\'enic fast solar wind. We show that protons locally conserve total kinetic energy when observed from an effective frame of reference comoving with the fluctuations. For typical properties of the fast wind, this frame can be reasonably identified by alpha particles, which, owing to their drift with respect to protons at about the Alfv\'en speed along the magnetic field, do not partake in the fluid low frequency fluctuations. Using their velocity to transform proton velocity into the frame of Alfv\'enic turbulence, we demonstrate that the resulting plasma motion is characterized by a constant absolute value of the velocity, zero electric fields, and aligned velocity and magnetic field vectors as expected for unidirectional Alfv\'enic fluctuations in equilibrium. We propose that this constraint, via the correlation between velocity and magnetic field in Alfv\'enic turbulence, is at the origin of the observed constancy of the magnetic field: while the constant velocity corresponding to constant energy can be only observed in the frame of the fluctuations, the correspondingly constant total magnetic field, invariant for Galilean transformations, remains the observational signature, in the spacecraft frame, of the constant total energy in the Alfv\'en turbulence frame.

The nature of the magnetic field fluctuations in the solar wind between the ion and electron scales is still under debate. Using the Cluster/STAFF instrument, we make a survey of the power spectral density and of the polarization of these fluctuations at frequencies $f\in[1,400]$ Hz, during five years (2001-2005), when Cluster was in the free solar wind. In $\sim 10\%$ of the selected data, we observe narrow-band, right-handed, circularly polarized fluctuations, with wave vectors quasi-parallel to the mean magnetic field, superimposed on the spectrum of the permanent background turbulence. We interpret these coherent fluctuations as whistler mode waves. The life time of these waves varies between a few seconds and several hours. Here we present, for the first time, an analysis of long-lived whistler waves, i.e. lasting more than five minutes. We find several necessary (but not sufficient) conditions for the observation of whistler waves, mainly a low level of the background turbulence, a slow wind, a relatively large electron heat flux and a low electron collision frequency. When the electron parallel beta factor $\beta_{e\parallel}$ is larger than 3, the whistler waves are seen along the heat flux threshold of the whistler heat flux instability. The presence of such whistler waves confirms that the whistler heat flux instability contributes to the regulation of the solar wind heat flux, at least for $\beta_{e\parallel} \ge$ 3, in the slow wind, at 1 AU.

The long-term evolution of large-amplitude Alfven waves propagating in the solar wind is investigated by performing two-dimensional MHD simulations within the expanding box model. The linear and nonlinear phases of the parametric decay instability are studied for both circularly polarized waves in parallel propagation and for arc-polarized waves in oblique propagation. The non-monochromatic case is also considered. In the oblique case, the direct excitation of daughter modes transverse to the local background field is found for the first time in an expanding environment, and this transverse cascade seems to be favored for monochromatic mother waves. The expansion effect reduces the instability growth rate, and it can even suppress its onset for the lowest frequency modes considered here, possibly explaining the persistence of these outgoing waves in the solar wind.

We investigate the role of kinetic instabilities driven by a proton anisotropy on the onset of magnetic reconnection by means of 2-D hybrid simulations. The collisionless tearing of a current sheet is studied in the presence of a proton temperature anisotropy in the surrounding plasma. Our results confirm that anisotropic protons within the current sheet region can significantly enhance/stabilize the tearing instability of the current. Moreover, fluctuations associated to linear instabilities excited by large proton temperature anisotropies can significantly influence the stability of the plasma and perturb the current sheets, triggering the tearing instability. We find that such a complex coupling leads to a faster tearing evolution in a regime with larger perpendicular temperature when an ion-cyclotron instability is generated by the anisotropic proton distribution functions. On the contrary, in the presence of the opposite anisotropy, fire hose fluctuations excited by the unstable background protons with larger parallel temperature are not able to efficiently destabilize the current sheets, which remain stable for a long time after fire hose saturation. We discuss possible influences of this novel coupling on the solar wind and heliospheric plasma dynamics.