Search

We perform a high-resolution two-dimensional fully-kinetic numerical simulation of a turbulent plasma system with observation-driven conditions, in order to investigate the interplay between turbulence, magnetic reconnection, and particle heating from ion to sub-electron scales in the near-Sun solar wind. We find that the power spectra of the turbulent plasma and electromagnetic fluctuations show multiple power-law intervals down to scales smaller than the electron gyroradius. Magnetic reconnection is observed to occur in correspondence of current sheets with a thickness of the order of the electron inertial length, which form and shrink due to interacting ion-scale vortexes. In some cases, both ion and electron outflows are observed (the classic reconnection scenario), while in others -- typically for the shortest current sheets -- only electron jets are presents ("electron-only reconnection"). At the onset of reconnection, the electron temperature starts to increase and a strong parallel temperature anisotropy develops. This suggests that in strong turbulence electron-scale coherent structures may play a significant role for electron heating, as impulsive and localized phenomena such as magnetic reconnection may transfer energy from the electromagnetic fields to particles more efficiently than damping mechanisms related to interactions with wave-like fluctuations., Comment: 18 pages, 6 figures

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

We investigate spectral properties of turbulence in the solar wind that is a weakly collisional astrophysical plasma, accessible to in-situ observations. Using the Helios search coil magnetometer measurements in the fast solar wind, in the inner heliosphere, we focus on properties of the turbulent magnetic fluctuations at scales smaller than the ion characteristic scales, the so-called kinetic plasma turbulence. At such small scales, we show that the magnetic power spectra between 0.3 and 0.9 AU from the Sun have a generic shape $\sim f^{-8/3}\exp{(-f/f_d)}$ where the dissipation frequency $f_d$ is correlated with the Doppler shifted frequency $f_{\rho e}$ of the electron Larmor radius. This behavior is statistically significant: all the observed kinetic spectra are well described by this model, with $f_d = f_{\rho e}/1.8$. Our results indicate that the electron gyroradius plays the role of the dissipation scale and marks the end of the electromagnetic cascade in the solar wind., Comment: Published in PhysRevE in June 2021 https://ui.adsabs.harvard.edu/abs/2021PhRvE.103f3202A/abstract. arXiv admin note: substantial text overlap with arXiv:2004.01102

Spatiotemporal properties of two-dimensional (2D) Hall-magnetohydrodynamic turbulence at intermediate plasma $\beta=2$ are studied by means of Fast Iterative Filtering, a new technique for the decomposition of nonstationary nonlinear signals. Results show that the magnetic energy at sub-ion scales is concentrated in perturbations with frequencies smaller than the ion-cyclotron (IC) frequency and with polarization properties that are incompatible with both kinetic Alfv\'en waves (KAWs) and IC waves. At higher frequencies, we clearly identify signatures of both whistler waves and KAWs, however their energetic contribution to the magnetic power spectrum is negligible. We conclude that the dynamics of 2D Hall-MHD turbulence at sub-ion scales is mainly driven by localized intermittent structures, with no significant contribution of wavelike fluctuations., Comment: 10 pages, 5 figures. Accepted for publication on The Astrophysical Journal Letters

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

How the turbulent energy cascade develops below the magnetohydrodynamic scales in space and astrophysical plasmas is a major open question. Here, we measure the power spectrum of magnetic fluctuations in Parker Solar Probe's observations close to the Sun and in state-of-the-art numerical simulations of plasma turbulence. Both reveal a power-law behavior with a slope compatible with $-11/3$ at scales smaller than the ion characteristic scales, steeper than what is typically observed in the solar wind and in the Earth's magnetosheath. We explain such behavior by developing a simple two-fluid model which does not require any kinetic processes nor electron-inertia effects. This is characterized by a significant contribution of the ion kinetic energy to the total turbulent energy cascade at sub-ion scales, although the dynamics is driven by the magnetic field through the current density. We expect that this regime may be relevant for a broad class of low-beta plasmas, e.g. the solar corona, non-relativistic magnetized jets and disks, and laboratory plasmas.

Turbulent space and astrophysical plasmas exhibit a complex dynamics, which involves nonlinear coupling across different temporal and spatial scales. There is growing evidence that impulsive events, such as magnetic reconnection instabilities, lead to a spatially localized enhancement of energy dissipation, thus speeding up the energy transfer at small scales. Capturing such a diverse dynamics is challenging. Here, we employ the Multidimensional Iterative Filtering (MIF) method, a novel technique for the analysis of nonstationary multidimensional signals. Unlike other traditional methods (e.g., based on Fourier or wavelet decomposition), MIF does not require any previous assumption on the functional form of the signal to be identified. Using MIF, we carry out a multiscale analysis of Hall-magnetohydrodynamic (HMHD) and hybrid particle-in-cell (HPIC) numerical simulations of decaying plasma turbulence. The results assess the ability of MIF to spatially identify and separate the different scales (the MHD inertial range, the sub-ion kinetic, and the dissipation scales) of the plasma dynamics. Furthermore, MIF decomposition allows to detect localized current structures and to characterize their contribution to the statistical and spectral properties of turbulence. Overall, MIF arises as a very promising technique for the study of turbulent plasma environments., Comment: 19 pages, 9 figures, submitted to Journal of Plasma Physics

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

Turbulence develops in any stressed flow when the scales of the forcing are much larger than those of the dissipation. In neutral fluids, it consists of chaotic motions in physical space but with a universal energy spectrum in Fourier space. Intermittency (non-Gaussian statistics of fluctuations) is another general property and it is related to the presence of coherent structures. Space plasmas are turbulent as well. Here, we focus on the kinetic plasma scales, which are not yet well understood. We address the following fundamental questions: (1) Do the turbulent fluctuations at kinetic scales form a universal spectrum? and (2) What is the nature of the fluctuations? Using measurements in the solar wind we show that the magnetic spectra of kinetic turbulence at 0.3, 0.6 and 0.9 AU from the Sun have the same shape as the ones close to the Earth orbit at 1 AU, indicating universality of the phenomenon. The fluctuations, which form this spectrum, are typically non-linearly interacting eddies that tend to generate magnetic filaments., Comment: 12 pages, 7 figures, submitted to Nature Communications on March 2, 2020

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

Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the energization of transrelativistic electrons in turbulence using hybrid particle-in-cell, which provide a realistic model of Alfv\'{e}nic turbulence from MHD to sub-ion scales, and test particle simulations for electrons. We find that, depending on the electron initial energy and turbulence strength, electrons may undergo a fast and efficient phase of energization due to the magnetic curvature drift during the time they are trapped in dynamic magnetic structures. In addition, electrons are accelerated stochastically which is a slower process that yields lower maximum energies. The combined effect of these two processes determines the overall electron acceleration. With appropriate turbulence parameters, we find that superthermal electrons can be accelerated up to relativistic energies. For example, with heliospheric parameters and a relatively high turbulence level, rapid energization to MeV energies is possible., Comment: Accepted for publication in 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

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

Properties of the electron mirror instability and its competition with the usually dominant whistler (electron cyclotron) instability driven by the electron perpendicular temperature anisotropy are investigated on the linear level using a Vlasov linear solver and on the nonlinear level using a two-dimensional full particle code. The simulation results show that the linearly subdominant electron mirror instability may compete on the nonlinear level with the whistler instability and may even become eventually the dominant mode that generates robust non-propagating sub-ion-scale coherent structures in the form of magnetic peaks., Comment: 12 pages, 7 figures, J. Plasma Phys

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

Measurements and simulations of inertial compressive turbulence in the solar wind are characterized by anti-correlated magnetic fluctuations parallel to the mean field and density structures. This signature has been interpreted as observational evidence for non-propagating pressure balanced structures (PBS), kinetic ion acoustic waves, as well as the MHD slow-mode. Given the high damping rates of parallel propagating compressive fluctuations, their ubiquity in satellite observations is surprising, and suggestive of a local driving process. One possible candidate for the generation of compressive fluctuations in the solar wind is Alfv\'en wave parametric instability. Here we test the parametric decay process as a source of compressive waves in the solar wind by comparing the collisionless damping rates of compressive fluctuations with the growth rates of the parametric decay instability daughter waves. Our results suggest that generation of compressive waves through parametric decay is overdamped at 1 AU, but that the presence of slow-mode like density fluctuations is correlated with the parametric decay of Alfv\'en waves.

We study the stability of large amplitude, circularly polarized Alfv\'en waves in an anisotropic plasma described by the double-adiabatic/CGL closure, and in particular the effect of a background thermal pressure anisotropy on the well-known properties of Alfv\'en wave parametric decay in Magnetohydrodynamics (MHD). Anisotropy allows instability over a much wider range of values of parallel plasma beta ($\beta_\parallel$) when $\xi=p_{0\bot}/p_{0\parallel} > 1$. When the pressure anisotropy exceeds a critical value, $\xi\geq\xi^*$ with $\xi^*\simeq2.7$, there is a new regime in which the parametric instability is no longer quenched at high $\beta_\parallel$ and in the limit $\beta_\parallel\gg1$ the growth rate becomes independent of $\beta_\parallel$. In the opposite case of $\xi< \xi^*$, the instability is strongly suppressed with increasing parallel plasma beta, similarly to the MHD case. We analyze marginal stability conditions for parametric decay in the $(\xi,\,\beta_\parallel)$ parameter space, and discuss possible implications for Alfv\'enic turbulence in the solar wind., Comment: 9 pages, 9 figures. Under review

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

Using two-dimensional hybrid expanding box simulations we study the competition between the continuously driven parallel proton temperature anisotropy and fire hose instabilities in collisionless homogeneous plasmas. For quasi radial ambient magnetic field the expansion drives $T_{\mathrm{p}\|}>T_{\mathrm{p}\perp}$ and the system becomes eventually unstable with respect to the dominant parallel fire hose instability. This instability is generally unable to counteract the induced anisotropization and the system typically becomes unstable with respect to the oblique fire hose instability later on. The oblique instability efficiently reduces the anisotropy and the system rapidly stabilizes while a significant part of the generated electromagnetic fluctuations is damped to protons. As long as the magnetic field is in the quasi radial direction, this evolution repeats itself and the electromagnetic fluctuations accumulate. For sufficiently oblique magnetic field the expansion drives $T_{\mathrm{p}\perp}>T_{\mathrm{p}\|}$ and brings the system to the stable region with respect to the fire hose instabilities., Comment: Journal of Plasma Physics, 14 pages, 9 figures

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

Using one-dimensional hybrid expanding box model we investigate properties of the solar wind in the outer heliosphere. We assume a proton-electron plasma with a strictly transverse ambient magnetic field and, beside the expansion, we take into account influence of a continuous injection of cold pick-up protons through the charge-exchange process between the solar wind protons and hydrogen of interstellar origin. The injected cold pick-up protons form a ring distribution function that rapidly becomes unstable and generate Alfv\'en cyclotron waves. The Alfv\'en cyclotron waves scatter pick-up protons to a spherical shell distribution function that thickens over that time owing to the expansion-driven cooling. The Alf\'ven cyclotron waves heat solar wind protons in the perpendicular direction (with respect to the ambient magnetic field) through the cyclotron resonance. At later times, the Alfv\'en cyclotron waves become parametrically unstable and the generated ion acoustic waves heat protons in the parallel direction through the Landau resonance. The resulting heating of the solar wind protons is efficient on the expansion time scale., Comment: ApJ, 7 pages, 7 figures

Proton and alpha particle collisional transport coefficients (isotropization, relative deceleration frequencies and heating rates) at 1 AU are quantified using the WIND/SWE data. In agreement with previous studies the ion-ion Coulomb collisions are generally important for slow solar wind streams and tend to reduce the temperature anisotropies, the differential streaming and the differences between proton and alpha particle temperatures. In slow solar wind streams the Coulomb collisions between protons and alpha particles are important for the overall proton energetics as well as for the relative deceleration between the two species. It is also shown that ion temperature anisotropies and differential streaming need to be generally taken into account for evaluation of the collisional transport coefficients., Comment: ApJ, 5 pages, 4 figures

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

Supercritical collisionless perpendicular shocks have an average macrostructure determined primarily by the dynamics of ions specularly reflected at the magnetic ramp. Within the overall macrostructure, instabilities, both linear and nonlinear, generate fluctuations and microstructure. To identify the sources of such microstructure, high-resolution two- and three-dimensional simulations have been carried out using the hybrid method, wherein the ions are treated as particles and the electron response is modelled as a massless fluid. We confirm the results of earlier 2-D simulations showing both field-parallel aligned propagating fluctuations and fluctuations carried by the reflected-gyrating ions. In addition, it is shown that, for 2-D simulations of the shock coplanarity plane, the presence of short-wavelength fluctuations in all magnetic components is associated with the ion Weibel instability driven at the upstream edge of the foot by the reflected-gyrating ions. In 3-D simulations we show for the first time that the dominant microstructure is due to a coupling between field-parallel propagating fluctuations in the ramp and the motion of the reflected ions. This results in a pattern of fluctuations counter-propagating across the surface of the shock at an angle inclined to the magnetic field direction, due to a combination of field-parallel motion at the Alfv\'en speed of the ramp, and motion in the sense of gyration of the reflected ions., Comment: Submitted to Journal of Plasma Physics

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

The measure of the third-order structure function, Y, is employed in the solar wind to compute the cascade rate of turbulence. In the absence of a mean field B0=0, Y is expected to be isotropic (radial) and independent of the direction of increments, so its measure yields directly the cascade rate. For turbulence with mean field, as in the solar wind, Y is expected to become more two dimensional (2D), that is, to have larger perpendicular components, loosing the above simple symmetry. To get the cascade rate one should compute the flux of Y, which is not feasible with single-spacecraft data, thus measurements rely upon assumptions about the unknown symmetry. We use direct numerical simulations (DNS) of magneto-hydrodynamic (MHD) turbulence to characterize the anisotropy of Y. We find that for strong guide field B0=5 the degree of two-dimensionalization depends on the relative importance of shear and pseudo polarizations (the two components of an Alfv\'en mode in incompressible MHD). The anisotropy also shows up in the inertial range. The more Y is 2D, the more the inertial range extent differs along parallel and perpendicular directions. We finally test the two methods employed in observations and find that the so-obtained cascade rate may depend on the angle between B0 and the direction of increments. Both methods yield a vanishing cascade rate along the parallel direction, contrary to observations, suggesting a weaker anisotropy of solar wind turbulence compared to our DNS. This could be due to a weaker mean field and/or to solar wind expansion., Comment: Some text editing and typos corrected, 13 pages, 6 figures, to be published in ApJ

Kinetic instabilities in weakly collisional, high beta plasmas are investigated using two-dimensional hybrid expanding box simulations with Coulomb collisions modeled through the Langevin equation (corresponding to the Fokker-Planck one). The expansion drives a parallel or perpendicular temperature anisotropy (depending on the orientation of the ambient magnetic field). For the chosen parameters the Coulomb collisions are important with respect to the driver but are not strong enough to keep the system stable with respect to instabilities driven by the proton temperature anisotropy. In the case of the parallel temperature anisotropy the dominant oblique fire hose instability efficiently reduces the anisotropy in a quasilinear manner. In the case of the perpendicular temperature anisotropy the dominant mirror instability generates coherent compressive structures which scatter protons and reduce the temperature anisotropy. For both the cases the instabilities generate temporarily enough wave energy so that the corresponding (anomalous) transport coefficients dominate over the collisional ones and their properties are similar to those in collisionless plasmas., Comment: in press in J. Plasma Phys

The proton temperature anisotropy in the solar wind exhibits apparent bounds which are compatible with the theoretical constraints imposed by temperature-anisotropy driven kinetic instabilities. Recent statistical analyses based on conditional averaging indicate that near these theoretical constraints the solar wind protons have typically enhanced temperatures and a weaker collisionality. Here we carefully analyze the solar wind data and show that these results are a consequence of superposition of multiple correlations in the solar wind, namely, they mostly result from the correlation between the proton temperature and the solar wind velocity and from the superimposed anti-correlation between the proton temperature anisotropy and the proton parallel beta in the fast solar wind. Colder and more collisional data are distributed around temperature isotropy whereas hotter and less collisional data have a wider range of the temperature anisotropy anti-correlated with the proton parallel beta with signatures of constraints owing to the temperature-anisotropy driven instabilities. However, most of the hot and weakly collisional data, including the hottest and least collisional ones, lies far from the marginal stability regions. Consequently, we conclude that there is no clear relation between the enhanced temperatures and instability constraints and that the conditional averaging used for these analyses must be used carefully and need to be well tested., Comment: Accepted for publication in ApJL

Driven by the support and interest of the international space plasma community to examine simultaneous physical plasma scales and their interactions, the Cross-Scale Mission concept was submitted and accepted as an ESA Cosmic Vision M-class candidate mission. This report presents an overview of the assessment study phase of the 7 ESA spacecraft Cross-Scale mission. Where appropriate, discussion of the benefit of international collaboration with the SCOPE mission, as well as other interested parties, is included., Comment: Assessment Study Report as part of the ESA Cosmic Vision M-Class mission selection process

UNKNOWN

UNKNOWN