Search

We investigate a class of ion-scale magnetic solitary structures in the solar wind, characterized by distinct magnetic field enhancements and bipolar rotations over spatial scales of several proton inertial lengths. Previously tentatively identified as Alfv\'enic solitons, these structures are revisited using high-resolution data from the Solar Orbiter and Parker Solar Probe missions. Using a machine learning-based method, we identified nearly a thousand such structures, providing new insights into their evolution and physical properties. Statistical analysis shows that these structures are more abundant closer to the Sun, with occurrence rates peaking around 30-40 solar radii and declining at greater distances, suggesting that they decay. High-cadence measurements reveal that these structures are predominantly found in low-beta environments, with consistent fluctuations in density, velocity, and magnetic field. Magnetic field enhancements are often accompanied by plasma density drops, which, under near pressure balance, limit field increases. This leads to small fractional field enhancements near the Sun (approximately 0.01 at 20 solar radii), making detection challenging. Magnetic field variance analysis indicates that these structures are primarily oblique to the local magnetic field. Alfv\'enic velocity-magnetic field correlations suggest that most of these structures propagate sunward in the plasma frame, distinguishing them from typical solar wind fluctuations. We compare these findings with previous studies, discussing possible generation mechanisms and their implications for the turbulent cascade in the near-Sun Alfv\'enic solar wind. Further high-resolution observations and simulations are needed to fully understand their origins and impacts., Comment: 11 pages, 4 figures. Submitted to The Astrophysical Journal Letters and currently under review

Where and under what conditions the transfer of energy between electromagnetic fields and particles takes place in the solar wind remains an open question. We investigate the conditions that promote the growth of kinetic instabilities predicted by linear theory, to infer how turbulence and temperature-anisotropy-driven instabilities are interrelated. Using a large dataset from Solar Orbiter, we introduce the radial rate of strain, a novel measure computed from single-spacecraft data, that we interpret as a proxy for the double-adiabatic strain rate. The solar wind exhibits high absolute values of the radial rate of strain at locations with large temperature anisotropy. We measure the kurtosis and skewness of the radial rate of strain from the statistical moments to show that it is non-Gaussian for unstable intervals and increasingly intermittent at smaller scales with a power-law scaling. We conclude that the velocity field fluctuations in the solar wind contribute to the presence of temperature anisotropy sufficient to create potentially unstable conditions., Comment: Accepted for publication in the Journal of Plasma Physics

Plasma in the terrestrial magnetosheath is characterised by very weak particle-particle collisions, so kinetic microinstabilities are thought to be responsible for regulating the thermodynamics of the plasma. By exciting electromagnetic waves, these instabilities redistribute free energy in velocity space, moulding the velocity distribution function (VDF) into a lower energy state. In the high-beta magnetosheath, relatively small perturbations to the VDF can easily excite instabilities compared to in the low-beta inner heliosphere. Since magnetic fields cannot do work on the particles, electric fields mediate energy exchange between the electromagnetic field and the bulk fluid properties of the plasma. We investigate signatures of non-ideal energy conversion associated with turbulent fluctuations in the context of electron and ion temperature anisotropy-beta instabilities, utilising over 24 hours of data spread over 163 distinct intervals of in situ magnetosheath observations from Magnetospheric Multiscale (MMS). We find that average energy conversion into fluid flow is enhanced along instability boundaries, suggesting that turbulence is playing a role in how free energy is redistributed in the plasma. The work enables a quantification of the energetics which are associated with the role of kinetic microinstabilities in regulating collisionless plasma thermodynamics. This work provides insight into the open question of how specific plasma processes couple into the turbulent dynamics and ultimately lead to energy dissipation and particle energisation in collisionless plasmas., Comment: 10 pages, 3 figures, 1 table. Submitted for peer review at Astrophysical Journal Letters (ApJL), July 2024

One of the primary science objectives of Parker Solar Probe (PSP) is to determine the structures and dynamics of the plasma and magnetic fields at the sources of the solar wind. However, establishing the connection between {\it in situ} measurements and structures and dynamics in the solar atmosphere is challenging: most of the magnetic footpoint mapping techniques have significant uncertainties in the source localization of a plasma parcel observed {\it in situ}, and the PSP plasma measurements suffer from a limited field of view. Therefore it is of interest to investigate whether {\it in situ} measurements can be used on their own to identify streams originating from the same structures in the corona more finely than the well known fast wind-coronal hole, slow wind-elsewhere distinction. Here we develop a novel time series visualization method \textcolor{red}{(time-frequency representation or TFR)} named Gaussianity Scalogram. Utilizing this method, by analyzing the magnetic magnitude data from both PSP and Ulysses, we successfully identify {\it in situ} structures that are possible remnants of solar atmospheric and magnetic structures spanning more than seven orders of magnitude, from years to seconds, including polar and mid-latitude coronal holes, as well as structures compatible with super-granulation , ``jetlets'' and ``picoflares''. \textcolor{red}{Furthermore, computer simulations of Alfv\'enic turbulence successfully reproduce the Gaussianization of the magnetic magnitude for locally homogeneous structures.} Building upon these discoveries, the Gaussianity Scalogram can help future studies to reveal the fractal-like fine structures in the solar wind time series from both PSP and decades-old data archive., Comment: Accepted by ApJL

We investigate properties of solar wind-like plasma turbulence using direct numerical simulations. We analyze the transition from large, magnetohydrodynamic (MHD) scales to the ion characteristic ones using two-dimensional hybrid (fluid electrons, kinetic ions) simulations. To capture and quantify turbulence properties, we apply the Karman-Howarth-Monin (KHM) equation for compressible Hall MHD (extended by considering the plasma pressure as a tensor quantity) to the numerical results. The KHM analysis indicates that the transition from MHD to ion scales (the so called ion break in the power spectrum) results from a combination of an onset of Hall physics and of an effective dissipation owing to the pressure-strain energy-exchange channel and resistivity. We discuss the simulation results in the context of the solar wind., Comment: 7 pages, 3 figures

Turbulent magnetic field fluctuations observed in the solar wind often maintain a constant magnitude condition accompanied by spherically polarized velocity fluctuations; these signatures are characteristic of large-amplitude Alfv\'{e}n waves. Nonlinear energy transfer in Alfv\'{e}nic turbulence is typically considered in the small-amplitude limit where the constant magnitude condition may be neglected; in contrast, nonlinear energy transfer in the large-amplitude limit remains relatively unstudied. We develop a method to analyze finite-amplitude turbulence through studying fluctuations as constant magnitude rotations in the stationary wave (de Hoffmann-Teller) frame, which reveals that signatures of finite-amplitude effects exist deep into the MHD range. While the dominant fluctuations are consistent with spherically-polarized large-amplitude Alfv\'{e}n waves, the subdominant mode is relatively compressible. Signatures of nonlinear interaction between the finite-amplitude spherically polarized mode with the subdominant population reveal highly aligned transverse components. In theoretical models of Alfv\'{e}nic turbulence, alignment is thought to reduce nonlinearity; our observations require that alignment is sufficient to either reduce shear nonlinearity such that non-Alfv\'{e}nic interactions may be responsible for energy transfer in spherically polarized states, or that counter-propagating fluctuations maintain anomalous coherence, which is a predicted signature of reflection-driven turbulence.

We analyze magnetic field data from the first six encounters of PSP, three Helios fast streams and two Ulysses south polar passes covering heliocentric distances $0.1\lesssim R\lesssim 3$ au. We use this data set to statistically determine the evolution of switchbacks of different periods and amplitudes with distance from the Sun. We compare the radial evolution of magnetic field variances with that of the mean square amplitudes of switchbacks, and quantify the radial evolution of the cumulative counts of switchbacks per km. We find that the amplitudes of switchbacks decrease faster than the overall turbulent fluctuations, in a way consistent with the radial decrease of the mean magnetic field. This could be the result of a saturation of amplitudes and may be a signature of decay processes of large amplitude Alfv\'enic fluctuations in the solar wind. We find that the evolution of switchback occurrence in the solar wind is scale-dependent: the fraction of longer duration switchbacks increases with radial distance whereas it decreases for shorter switchbacks. This implies that switchback dynamics is a complex process involving both decay and in-situ generation in the inner heliosphere. We confirm that switchbacks can be generated by the expansion although other type of switchbacks generated closer to the sun cannot be ruled out., Comment: Accepted for publication on ApJL

The Parker Solar Probe (PSP) mission presents a unique opportunity to study the near-Sun solar wind closer than any previous spacecraft. During its fourth and fifth solar encounters, PSP had the same orbital trajectory, meaning that solar wind was measured at the same latitudes and radial distances. We identify two streams measured at the same heliocentric distance ($\sim$0.13au) and latitude ($\sim$-3.5$^{\circ}$) across these encounters to reduce spatial evolution effects. By comparing the plasma of each stream, we confirm that they are not dominated by variable transient events, despite PSP's proximity to the heliospheric current sheet. Both streams are consistent with a previous slow Alfv\'enic solar wind study once radial effects are considered, and appear to originate at the Southern polar coronal hole boundary. We also show that the switchback properties are not distinctly different between these two streams. Low $\alpha$-particle abundance ($\sim$ 0.6 %) is observed in the encounter 5 stream, suggesting that some physical mechanism must act on coronal hole boundary wind to cause $\alpha$-particle depletion. Possible explanations for our observations are discussed, but it remains unclear whether the depletion occurs during the release or the acceleration of the wind. Using a flux tube argument, we note that an $\alpha$-particle abundance of $\sim$ 0.6 % in this low velocity wind could correspond to an abundance of $\sim$ 0.9 % at 1 au. Finally, as the two streams roughly correspond to the spatial extent of a switchback patch, we suggest that patches are distinct features of coronal hole wind., Comment: 9 pages, 5 figures

We derive two new forms of the K\'arm\'an-Howarth-Monin equation for decaying compressible Hall magnetohydrodynamic (MHD) turbulence. We test them on results of a weakly-compressible, two-dimensional, moderate-Reynolds-number Hall MHD simulation and compare them with an isotropic spectral transfer (ST) equation. The KHM and ST equations are automatically satisfied during the whole simulation owing to the periodic boundary conditions and have complementary cumulative behavior. They are used here to analyze the onset of turbulence and its properties when it is fully developed. These approaches give equivalent results characterizing: the decay of the kinetic + magnetic energy at large scales, the MHD and Hall cross-scale energy transfer/cascade, the pressure dilatation, and the dissipation. The Hall cascade appears when the MHD one brings the energy close to the ion inertial range and is related to the formation of reconnecting current sheets. At later times, the pressure-dilation energy-exchange rate oscillates around zero with no net effect on the cross-scale energy transfer when averaged over a period of its oscillations. A reduced one-dimensional analysis suggests that all three methods may be useful to estimate the energy cascade rate from in situ observations., Comment: 11 pages, 7 figues

Parker Solar Probe (PSP) routinely observes magnetic field deflections in the solar wind at distances less than 0.3 au from the Sun. These deflections are related to structures commonly called 'switchbacks' (SBs), whose origins and characteristic properties are currently debated. Here, we use a database of visually selected SB intervals - and regions of solar wind plasma measured just before and after each SB - to examine plasma parameters, turbulent spectra from inertial to dissipation scales, and intermittency effects in these intervals. We find that many features, such as perpendicular stochastic heating rates and turbulence spectral slopes are fairly similar inside and outside of SBs. However, important kinetic properties, such as the characteristic break scale between the inertial to dissipation ranges differ inside and outside these intervals, as does the level of intermittency, which is notably enhanced inside SBs and in their close proximity, most likely due to magnetic field and velocity shears observed at the edges. We conclude that the plasma inside and outside of a SB, in most of the observed cases, belongs to the same stream, and that the evolution of these structures is most likely regulated by kinetic processes, which dominate small scale structures at the SB edges.

One of the main discoveries from the first two orbits of Parker Solar Probe (PSP) was the presence of magnetic switchbacks, whose deflections dominated the magnetic field measurements. Determining their shape and size could provide evidence of their origin, which is still unclear. Previous work with a single solar wind stream has indicated that these are long, thin structures although the direction of their major axis could not be determined. We investigate if this long, thin nature extends to other solar wind streams, while determining the direction along which the switchbacks within a stream were aligned. We try to understand how the size and orientation of the switchbacks, along with the flow velocity and spacecraft trajectory, combine to produce the observed structure durations for past and future orbits. We searched for the alignment direction that produced a combination of a spacecraft cutting direction and switchback duration that was most consistent with long, thin structures. The expected form of a long, thin structure was fitted to the results of the best alignment direction, which determined the width and aspect ratio of the switchbacks for that stream. The switchbacks had a mean width of $50,000 \, \rm{km}$, with an aspect ratio of the order of $10$. We find that switchbacks are not aligned along the background flow direction, but instead aligned along the local Parker spiral, perhaps suggesting that they propagate along the magnetic field. Since the observed switchback duration depends on how the spacecraft cuts through the structure, the duration alone cannot be used to determine the size or influence of an individual event. For future PSP orbits, a larger spacecraft transverse component combined with more radially aligned switchbacks will lead to long duration switchbacks becoming less common.

During Parker Solar Probe's first two orbits there are widespread observations of rapid magnetic field reversals known as switchbacks. These switchbacks are extensively found in the near-Sun solar wind, appear to occur in patches, and have possible links to various phenomena such as magnetic reconnection near the solar surface. As switchbacks are associated with faster plasma flows, we questioned whether they are hotter than the background plasma and whether the microphysics inside a switchback is different to its surroundings. We have studied the reduced distribution functions from the Solar Probe Cup instrument and considered time periods with markedly large angular deflections, to compare parallel temperatures inside and outside switchbacks. We have shown that the reduced distribution functions inside switchbacks are consistent with a rigid phase space rotation of the background plasma. As such, we conclude that the proton core parallel temperature is the same inside and outside of switchbacks, implying that a T-V relationship does not hold for the proton core parallel temperature inside magnetic field switchbacks. We further conclude that switchbacks are consistent with Alfv\'enic pulses travelling along open magnetic field lines. The origin of these pulses, however, remains unknown. We also found that there is no obvious link between radial Poynting flux and kinetic energy enhancements suggesting that the radial Poynting flux is not important for the dynamics of switchbacks., Comment: Submitted to Monthly Notices of the Royal Astronomical Society

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.

Properties of the turbulent cascade of kinetic energy are studied using direct numerical simulations of three-dimensional hydrodynamic decaying turbulence with a moderate Reynolds number and the initial Mach number $M=1$. Compressible and incompressible versions of the Karman-Howarth-Monin (KHM) and low-pass filtering/coarse-graining approaches are compared. In the simulation the total energy is well conserved; the scale dependent KHM and coarse-grained energy equations are also well conserved; the two approaches show similar results, the system does not have an inertial range for the cascade of kinetic energy, the region where this cascade dominates also have a non-negligible contribution of the kinetic-energy decay, dissipation, and pressure-dilatation effects. While the two approaches give semi-quantitatively similar results for the kinetic energy cascade, dissipation and pressure-dilatation rates, they differ in the increment separation and filtering scales; these scales are not simply related. The two approaches may be used to find the inertial range and to determine the cascade/dissipation rate of the kinetic energy., Comment: 13 pages, 5 figures

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).

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., Comment: Accepted for publication on APJ PSP special issue

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., Comment: 11 pages, 6 figures, submitted to The Astrophysical Journal

We present estimates of the turbulent energy cascade rate, derived from a Hall-MHD third-order law. We compute the contribution from the Hall term and the MHD term to the energy flux. We use MMS data accumulated in the magnetosheath and the solar wind, and compare the results with previously established simulation results. We find that in observation, the MHD contribution is dominant at inertial scales, as in the simulations, but the Hall term becomes significant in observations at larger scales than in the simulations. Possible reasons are offered for this unanticipated result., Comment: Accepted for Publication in Physical Review Letters

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., Comment: accepted in MNRAS

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., Comment: 9 pages, 4 figures

We present a local approach to the study of scale-to-scale energy transfers in magnetohydrodynamic (MHD) turbulence. This approach is based on performing local averages of the physical fields, which amounts to filtering scales smaller than some parameter $\ell$. A key step in this work is the derivation of a local K\'arm\'an-Howarth-Monin relation which provides a local form of Politano and Pouquet's 4/3-law, without any assumption of homogeneity or isotropy. Our approach is exact, non-random, and we show its connection to the usual statistical results of turbulence. Its implementation on data obtained via a three dimensional direct numerical simulation of the forced, incompressible MHD equations from the John Hopkins turbulence database constitutes the main part of our study. First, we show that the local K\'arm\'an-Howarth-Monin relation holds well. The space statistics of local cross-scale transfers is studied next, their means and standard deviations being maximum at inertial scales, and their probability density functions (PDFs) displaying very wide tails. Events constituting the tails of the PDFs are shown to form structures of strong transfers, either positive or negative, and which can be observed over the whole available range of scales. As $\ell$ is decreased, these structures become more and more localized in space while contributing to an increasing fraction of the mean energy cascade rate. Finally, we highlight their quasi 1D (filament-like) or quasi 2D (sheet/ribbon-lke) nature, and show that they appear in areas of strong vorticity or electric current density.

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., Comment: 6 pages, 5 figures, Accepted by The Astrophysical Journal Letters

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.

We investigate the properties of plasma turbulence by means of two-dimensional Hall-magnetohydrodynamic (HMHD) and hybrid particle-in-cell (HPIC) numerical simulations. We find that HMHD simulations exhibit spectral properties that are in most cases in agreement with the results of the HPIC simulations and with solar wind observations. The energy spectra of magnetic fluctuations exhibit a double power-law with spectral index $-5/3$ at MHD scales and $-3$ at kinetic scales, while for velocity fluctuations the spectral index is $-3/2$ at MHD scales. The break between the MHD and the kinetic scales occurs at the same scale in both simulations. In the MHD range the slopes of the total energy and residual energy spectra satisfy a fast Alfv\'en-dynamo balance. The development of a turbulent cascade is concurrently characterized by magnetic reconnection events taking place in thin current sheets that form between large eddies. A statistical analysis reveals that reconnection is qualitatively the same and fast in both the HMHD and HPIC models, characterized by inverse reconnection rates much smaller than the characteristic large-eddy nonlinear time. The agreement extends to other statistical properties, such us the kurtosis of the magnetic field. Moreover, the observation of a direct energy transfer from the large vortices to the small sub-ion scales triggered by magnetic reconnection, further supports the existence of a reconnection-mediated turbulent regime at kinetic scales. We conclude that the Hall-MHD fluid description captures to a large extent the transition of the turbulent cascade between the large MHD scales and the sub-ion scales., Comment: 12 pages, 8 figures

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., Comment: 3 Figures, 5 pages, ApJL, submitted

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., Comment: 11 pages, 8 figures, submitted to The Astrophysical Journal

We present the first statistical study of the anisotropy of the magnetic field turbulence in the solar wind between 1 and 200Hz, i.e. from proton to sub-electron scales. We consider 93 10-minute intervals of Cluster/STAFF measurements. We find that the fluctuations $\delta B_{\perp}^2$ are not gyrotropic at a given frequency $f$, a property already observed at larger scales ($\parallel$/$\perp$ mean parallel/perpendicular to the average magnetic ${\bf B_0}$). This non-gyrotropy gives indications on the angular distribution of the wave vectors ${\bf k}$: at $f <$ 10Hz, we find that $k_{\perp}\gg k_{\parallel}$, mainly in the fast wind; at $f >$ 10Hz, fluctuations with a non-negligible $k_{\parallel}$ are also present. We then consider the anisotropy ratio $\delta B_{\parallel}^2/\delta B_{\perp}^2$, which is a measure of the magnetic compressibility of the fluctuations. This ratio, always smaller than 1, increases with $f.$ It reaches a value showing that the fluctuations are more or less isotropic at electron scales, for $f \geq 50$Hz. From 1 to 15-20Hz, there is a strong correlation between the observed compressibility and the one expected for the kinetic Alfv\'en waves (KAW), which only depends on the total plasma $\beta$. For $f > $ 15-20Hz, the observed compressibility is larger than expected for KAW; and it is stronger in the slow wind: this could be an indication of the presence of a slow-ion acoustic mode of fluctuations, which is more compressive and is favored by the larger values of the electron to proton temperature ratio generally observed in the slow wind., Comment: 13 pages, 12 figures, accepted for publication in ApJ

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., Comment: Solar Physics, accepted

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., Comment: 8 pages, 4 figures, accepted for publication on The Astrophysical Journal

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., Comment: 4 pages, 2 figures, Proceedings of the Fourteenth International Solar Wind Conference

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., Comment: Accepted for publication in The Astrophysical Journal, March 21, 2016

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., Comment: 10 pages, 9 figures, Accepted for publication in The Astrophysical Journal

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., Comment: 6 pages, 5 figures, in press in ApJL

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., Comment: 15 pages, 13 figures, submitted to ApJ

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., Comment: 5 pages, 5 figures; introduction and conclusions changed, references updated, accepted for publication in ApJL

We present the first three-dimensional hybrid simulations of the evolution of ion-scale current sheets, with an investigation of the role of temperature anisotropy and associated kinetic instabilities on the growth of the tearing instability and particle heating. We confirm the ability of the ion cyclotron and firehose instabilities to enhance or suppress reconnection, respectively. The simulations demonstrate the emergence of persistent three-dimensional structures, including patchy reconnection sites and the fast growth of a narrow-band drift-kink instability, which suppresses reconnection for thin current sheets with weak guide fields. Potential observational signatures of the three-dimensional evolution of solar wind current sheets are also discussed. We conclude that kinetic instabilities, arising from non-Maxwellian ion populations, are significant to the evolution of three-dimensional current sheets, and two-dimensional studies of heating rates by reconnection may therefore over-estimate the ability of thin, ion-scale current sheets to heat the solar wind by reconnection.

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., Comment: The Astrophysical Journal, 2014, in press

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., Comment: 6 pages, 6 figures, Accepted for publication in the Astrophysical Journal