The magnetization $|\Omega_{\mathrm e}|/\omega_{\mathrm{e}}$ is an important parameter in plasma astrophysics, where $\Omega_{\mathrm e}$ and $\omega_{\mathrm{e}}$ are the electron gyro-frequency and electron plasma frequency, respectively. It only depends on the mass ratio $m_{\mathrm i}/m_{\mathrm e}$ and the light-to-Alfv\'en speed ratio $c/v_{\mathrm{Ai}}$, where $m_{\mathrm i}$ ($m_{\mathrm e}$) is the ion (electron) mass, $c$ is the speed of light, and $v_{\mathrm{Ai}}$ is the ion Alfv\'en speed. Nonlinear numerical plasma models such as particle-in-cell simulations must often assume unrealistic values for $m_{\mathrm i}/m_{\mathrm e}$ and for $c/v_{\mathrm{Ai}}$. Because linear theory yields exact results for parametric scalings of wave properties at small amplitudes, we use linear theory to investigate the dispersion relations of Alfv\'en/ion-cyclotron and fast-magnetosonic/whistler waves as prime examples for collective plasma behaviour depending on $m_{\mathrm i}/m_{\mathrm e}$ and $c/v_{\mathrm{Ai}}$. We analyse their dependence on $m_{\mathrm i}/m_{\mathrm e}$ and $c/v_{\mathrm{Ai}}$ in quasi-parallel and quasi-perpendicular directions of propagation with respect to the background magnetic field for a plasma with $\beta_j\sim1$, where $\beta_j$ is the ratio of the thermal to magnetic pressure for species $j$. Although their dispersion relations are largely independent of $c/v_{\mathrm{Ai}}$ for $c/v_{\mathrm{Ai}}\gtrsim 10$, the mass ratio $m_{\mathrm i}/m_{\mathrm e}$ has a strong effect at scales smaller than the ion inertial length. Moreover, we study the impact of relativistic electron effects on the dispersion relations. Based on our results, we recommend aiming for a more realistic value of $m_{\mathrm i}/m_{\mathrm e}$ than for a more realistic value of $c/v_{\mathrm{Ai}}$ in non-relativistic plasma simulations if such a choice is necessary, although $\dots$

The solar wind is a magnetized plasma and as such exhibits collective plasma behavior associated with its characteristic spatial and temporal scales. The characteristic length scales include the size of the heliosphere, the collisional mean free paths of all species, their inertial lengths, their gyration radii, and their Debye lengths. The characteristic timescales include the expansion time, the collision times, and the periods associated with gyration, waves, and oscillations. We review the past and present research into the multi-scale nature of the solar wind based on in-situ spacecraft measurements and plasma theory. We emphasize that couplings of processes across scales are important for the global dynamics and thermodynamics of the solar wind. We describe methods to measure in-situ properties of particles and fields. We then discuss the role of expansion effects, non-equilibrium distribution functions, collisions, waves, turbulence, and kinetic microinstabilities for the multi-scale plasma evolution.

We investigate the scattering of strahl electrons by microinstabilities as a mechanism for creating the electron halo in the solar wind. We develop a mathematical framework for the description of electron-driven microinstabilities and discuss the associated physical mechanisms. We find that an instability of the oblique fast-magnetosonic/whistler (FM/W) mode is the best candidate for a microinstability that scatters strahl electrons into the halo. We derive approximate analytic expressions for the FM/W instability threshold in two different $\beta_{\mathrm c}$ regimes, where $\beta_{\mathrm c}$ is the ratio of the core electrons' thermal pressure to the magnetic pressure, and confirm the accuracy of these thresholds through comparison with numerical solutions to the hot-plasma dispersion relation. We find that the strahl-driven oblique FM/W instability creates copious FM/W waves under low-$\beta_{\mathrm c}$ conditions when $U_{0\mathrm s}\gtrsim 3w_{\mathrm c}$, where $U_{0\mathrm s}$ is the strahl speed and $w_{\mathrm c}$ is the thermal speed of the core electrons. These waves have a frequency of about half the local electron gyrofrequency. We also derive an analytic expression for the oblique FM/W instability for $\beta_{\mathrm c}\sim 1$. The comparison of our theoretical results with data from the \emph{Wind} spacecraft confirms the relevance of the oblique FM/W instability for the solar wind. The whistler heat-flux, ion-acoustic heat-flux, kinetic-Alfv\'en-wave heat-flux, and electrostatic electron-beam instabilities cannot fulfill the requirements for self-induced scattering of strahl electrons into the halo. We make predictions for the electron strahl close to the Sun, which will be tested by measurements from \emph{Parker Solar Probe} and \emph{Solar Orbiter}.

A multi-institutional, multi-national science team will soon submit a NASA proposal to build a constellation of spacecraft to fly into the near-Earth solar wind in a swarm spanning a multitude of scales in order to obtain critically needed measurements that will reveal the underlying dynamics of magnetized turbulence. This white paper, submitted to the Plasma 2020 Decadal Survey Committee, provides a brief overview of turbulent systems that constitute an area of compelling plasma physics research, including why this mission is needed, and how this mission will achieve the goal of revealing how energy is transferred across scales and boundaries in plasmas throughout the universe.

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

Using 1D and 2D particle-in-cell (PIC) simulations of a plasma with a growing magnetic field $\vec{B}$, we show that ions can be stochastically accelerated by the ion-cyclotron (IC) instability. As $\vec{B}$ grows, an ion pressure anisotropy $p_{\perp,i} > p_{||,i}$ arises, due to the adiabatic invariance of the ion magnetic moment ($p_{||,i}$ and $p_{\perp,i}$ are the ion pressures parallel and perpendicular to $\vec{B}$). When initially $\beta_i = 0.5$ ($\beta_i \equiv 8\pi p_i/|\vec{B}|^2$, where $p_i$ is the ion isotropic pressure), the pressure anisotropy is limited mainly by inelastic pitch-angle scattering provided by the IC instability, which in turn produces a non-thermal tail in the ion energy spectrum. After $\vec{B}$ is amplified by a factor $\sim 2.7$, this tail can be approximated as a power-law of index $\sim 3.4$ plus two non-thermal bumps, and accounts for $2-3\%$ of the ions and $\sim 18\%$ of their kinetic energy. On the contrary, when initially $\beta_i =2$, the ion scattering is dominated by the mirror instability and the acceleration is suppressed. This implies that efficient ion acceleration requires that initially $\beta_i \lesssim 1$. Although we focus on cases where $\vec{B}$ is amplified by plasma shear, we check that the acceleration occurs similarly if $\vec{B}$ grows due to plasma compression. Our results are valid in a sub-relativistic regime where the ion thermal energy is $\sim 10\%$ of the ion rest mass energy. This acceleration process can thus be relevant in the inner region of low-luminosity accretion flows around black holes.

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

This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.

Kinetic Alfv\'{e}n waves (KAWs) are the short-wavelength extension of the MHD Alfv\'{e}n-wave branch in the case of highly-oblique propagation with respect to the background magnetic field. Observations of space plasma show that small-scale turbulence is mainly KAW-like. We apply two theoretical approaches, collisional two-fluid theory and collisionless kinetic theory, to obtain predictions for the KAW polarizations depending on $\beta_\mathrm{p}$ (the ratio of the proton thermal pressure to the magnetic pressure) at the ion gyroscale in terms of fluctuations in density, bulk velocity, and pressure. We perform a wavelet analysis of MMS magnetosheath measurements and compare the observations with both theories. We find that the two-fluid theory predicts the observations better than kinetic theory, suggesting that the small-scale KAW-like fluctuations exhibit a fluid-like behavior in the magnetosheath although the plasma is weakly collisional. We also present predictions for the KAW polarizations in the inner heliosphere that are testable with Parker Solar Probe and Solar Orbiter.

Stochastic heating refers to an increase in the average magnetic moment of charged particles interacting with electromagnetic fluctuations whose frequencies are smaller than the particles' cyclotron frequencies. This type of heating arises when the amplitude of the gyroscale fluctuations exceeds a certain threshold, causing particle orbits in the plane perpendicular to the magnetic field to become stochastic rather than nearly periodic. We consider the stochastic heating of protons by Alfv\'en-wave (AW) and kinetic-Alfv\'en-wave (KAW) turbulence, which may make an important contribution to the heating of the solar wind. Using phenomenological arguments, we derive the stochastic-proton-heating rate in plasmas in which $\beta_{\rm p} \sim 1-30$, where $\beta_{\rm p}$ is the ratio of the proton pressure to the magnetic pressure. (We do not consider the $\beta_{\rm p} \gtrsim 30$ regime, in which KAWs at the proton gyroscale become non-propagating.) We test our formula for the stochastic-heating rate by numerically tracking test-particle protons interacting with a spectrum of randomly phased AWs and KAWs. Previous studies have demonstrated that at $\beta_{\rm p} \lesssim 1$, particles are energized primarily by time variations in the electrostatic potential and thermal-proton gyro-orbits are stochasticized primarily by gyroscale fluctuations in the electrostatic potential. In contrast, at $\beta_{\rm p} \gtrsim 1$, particles are energized primarily by the solenoidal component of the electric field and thermal-proton gyro-orbits are stochasticized primarily by gyroscale fluctuations in the magnetic field.

NHDS is the New Hampshire Dispersion Relation Solver. This article describes the numerics of the solver and its capabilities. The code is available for download on https://github.com/danielver02/NHDS.

The Arbitrary Linear Plasma Solver (ALPS) is a parallelised numerical code that solves the dispersion relation in a hot (even relativistic) magnetised plasma with an arbitrary number of particle species with arbitrary gyrotropic equilibrium distribution functions for any direction of wave propagation with respect to the background field. ALPS reads the background momentum distributions as tables of values on a $(p_{\perp},p_{\parallel})$ grid, where $p_{\perp}$ and $p_{\parallel }$ are the momentum coordinates in the directions perpendicular and parallel to the background magnetic field, respectively. We present the mathematical and numerical approach used by ALPS and introduce our algorithms for the handling of poles and the analytic continuation for the Landau contour integral. We then show test calculations of dispersion relations for a selection of stable and unstable configurations in Maxwellian, bi-Maxwellian, $\kappa$-distributed, and J\"uttner-distributed plasmas. These tests demonstrate that ALPS derives reliable plasma dispersion relations. ALPS will make it possible to determine the properties of waves and instabilities in the non-equilibrium plasmas that are frequently found in space, laboratory experiments, and numerical simulations.

We use particle-in-cell (PIC) simulations of a collisionless, electron-ion plasma with a decreasing background magnetic field, $B$, to study the effect of velocity-space instabilities on the viscous heating and thermal conduction of the plasma. If $B$ decreases, the adiabatic invariance of the magnetic moment gives rise to pressure anisotropies with $p_{||,j} > p_{\perp,j}$ ($p_{||,j}$ and $p_{\perp,j}$ represent the pressure of species $j$ ($=i$ or $e$) parallel and perpendicular to the magnetic field). Linear theory indicates that, for sufficiently large anisotropies, different velocity-space instabilities can be triggered. These instabilities, which grow on scales comparable to the electron and ion Larmor radii, in principle have the ability to pitch-angle scatter the particles, limiting the growth of the anisotropies. Our PIC simulations focus on the nonlinear, saturated regime of the instabilities. This is done through the permanent decrease of the magnetic field by an imposed shear in the plasma. Our results show that, in the regime $2 \lesssim \beta_j \lesssim 20$ ($\beta_j \equiv 8\pi p_j/B^2$), the saturated ion and electron pressure anisotropies are controlled by the combined effect of the oblique ion firehose (OIF) and the fast magnetosonic/whistler (FM/W) instabilities. These instabilities grow preferentially on the ion Larmor radius scale, and make the ion and electron pressure anisotropies nearly equal: $\Delta p_e/p_{||,e} \approx \Delta p_i/p_{||,i}$ (where $\Delta p_j=p_{\perp,j} - p_{||,j}$). We also quantify the thermal conduction of the plasma by directly calculating the mean free path of electrons along the mean magnetic field, which we find strongly depends on whether $B$ decreases or increases. Our results can be applied in studies of low collisionality plasmas such as the solar wind, the intracluster medium, and some accretion disks around black holes.

Observations in the solar wind suggest that the compressive component of inertial-range solar-wind turbulence is dominated by slow modes. The low collisionality of the solar wind allows for non-thermal features to survive, which suggests the requirement of a kinetic plasma description. The least-damped kinetic slow mode is associated with the ion-acoustic (IA) wave and a non-propagating (NP) mode. We derive analytical expressions for the IA-wave dispersion relation in an anisotropic plasma in the framework of gyrokinetics and then compare them to fully-kinetic numerical calculations, results from two-fluid theory, and MHD. This comparison shows major discrepancies in the predicted wave phase speeds from MHD and kinetic theory at moderate to high $\beta$. MHD and kinetic theory also dictate that all plasma normal modes exhibit a unique signature in terms of their polarization. We quantify the relative amplitude of fluctuations in the three lowest particle velocity moments associated with IA and NP modes in the gyrokinetic limit and compare these predictions with MHD results and in-situ observations of the solar-wind turbulence. The agreement between the observations of the wave polarization and our MHD predictions is better than the kinetic predictions, suggesting that the plasma behaves more like a fluid in the solar wind than expected.

Compressive fluctuations are a minor yet significant component of astrophysical plasma turbulence. In the solar wind, long-wavelength compressive slow-mode fluctuations lead to changes in $\beta_{\parallel \mathrm p}\equiv 8\pi n_{\mathrm p}k_{\mathrm B}T_{\parallel \mathrm p}/B^2$ and in $R_{\mathrm p}\equiv T_{\perp \mathrm p}/T_{\parallel \mathrm p}$, where $T_{\perp \mathrm p}$ and $T_{\parallel \mathrm p}$ are the perpendicular and parallel temperatures of the protons, $B$ is the magnetic field strength, and $n_{\mathrm p}$ is the proton density. If the amplitude of the compressive fluctuations is large enough, $R_{\mathrm p}$ crosses one or more instability thresholds for anisotropy-driven microinstabilities. The enhanced field fluctuations from these microinstabilities scatter the protons so as to reduce the anisotropy of the pressure tensor. We propose that this scattering drives the average value of $R_{\mathrm p}$ away from the marginal stability boundary until the fluctuating value of $R_{\mathrm p}$ stops crossing the boundary. We model this "fluctuating-anisotropy effect" using linear Vlasov--Maxwell theory to describe the large-scale compressive fluctuations. We argue that this effect can explain why, in the nearly collisionless solar wind, the average value of $R_{\mathrm p}$ is close to unity.

In low-collisionality plasmas, velocity-space instabilities are a key mechanism providing an effective collisionality for the plasma. We use particle-in-cell (PIC) simulations to study the interplay between electron and ion-scale velocity-space instabilities and their effect on electron pressure anisotropy, viscous heating, and thermal conduction. The adiabatic invariance of the magnetic moment in low-collisionality plasmas leads to pressure anisotropy, $p_{\perp,j} > p_{||,j}$, if the magnetic field $\vec{B}$ is amplified ($p_{\perp,j}$ and $p_{||,j}$ denote the pressure of species $j$ [electron, ion] perpendicular and parallel to $\vec{B}$). If the resulting anisotropy is large enough, it can in turn trigger small-scale plasma instabilities. Our PIC simulations explore the nonlinear regime of the mirror, ion-cyclotron, and electron whistler instabilities, through continuous amplification of the magnetic field $|\vec{B}|$ by an imposed shear in the plasma. In the regime $1 \lesssim \beta_j \lesssim 20$ ($\beta_j \equiv 8\pi p_j/|\vec{B}|^2$), the saturated electron pressure anisotropy, $\Delta p_e/p_{||,e}$, is determined mainly by the (electron-lengthscale) whistler marginal stability condition, with a modest factor of $\sim 1.5-2$ decrease due to the trapping of electrons by the mirrors. We explicitly calculate the mean free path of the electrons and ions along the mean magnetic field and provide a simple physical prescription for the mean free path and thermal conductivity in low-collisionality $\beta_j \gtrsim 1$ plasmas. Our results imply that velocity-space instabilities likely decrease the thermal conductivity of plasma in the outer parts of massive, hot, galaxy clusters. We also discuss the implications of our results for electron heating and thermal conduction in low-collisionality accretion flows onto black holes, including Sgr A* in the Galactic Center.

We consider a fast magnetosonic multifluid shock as a representation of the solar-wind termination shock. We assume the action of the transition happens in a three-step process: In the first step, the upstream supersonic solar-wind plasma is subject to a strong electric field that flashes up on a small distance scale $\Delta z\simeq U_1/ \Omega _{\mathrm e}$ (first part of the transition layer), where $\Omega_{\mathrm e}$ is the electron gyro-frequency and $U_1$ is the upstream speed. This electric field both decelerates the supersonic ion flow and accelerates the electrons up to high velocities. In this part of the transition region, the electric forces connected with the deceleration of the ion flow strongly dominate over the Lorentz forces. We, therefore, call this part the demagnetization region. In the second phase, Lorentz forces due to convected magnetic fields compete with the electric field, and the highly anisotropic and energetic electron distribution function is converted into a shell distribution with energetic shell electrons storing about 3/4 of the upstream ion kinetic energy. In the third phase, the plasma particles thermalize due to the relaxation of free energy by plasma instabilities. The first part of the transition region opens up a new thermodynamic degree of freedom never before taken into account for the electrons, since the electrons are usually considered to be enslaved to follow the behavior of the protons in all velocity moments like density, bulk velocity, and temperature. We show that electrons may be the downstream plasma fluid that dominates the downstream plasma pressure.

The interaction between Alfven-wave turbulence and the background solar wind affects the cross helicity in two ways. Non-WKB reflection converts outward-propagating Alfven waves into inward-propagating Alfven waves and vice versa, and the turbulence transfers momentum to the background flow. When both effects are accounted for, the total cross helicity is conserved. In the special case that the background density and flow speed are independent of time, the equations of cross-helicity conservation and total-energy conservation can be combined to recover a well-known equation derived by Heinemann and Olbert that has been interpreted as a non-WKB generalization of wave-action conservation. This latter equation (in contrast to cross-helicity and energy conservation) does not hold when the background varies in time.

In the majority of the literature on plasma shock waves, electrons play the role of "ghost particles," since their contribution to mass and momentum flows is negligible, and they have been treated as only taking care of the electric plasma neutrality. In some more recent papers, however, electrons play a new important role in the shock dynamics and thermodynamics, especially at the solar-wind termination shock. They react on the shock electric field in a very specific way, leading to suprathermal nonequilibrium distributions of the downstream electrons, which can be represented by a kappa distribution function. In this paper, we discuss why this anticipated hot electron population has not been seen by the plasma detectors of the Voyager spacecraft downstream of the solar-wind termination shock. We show that hot nonequilibrium electrons induce a strong negative electric charge-up of any spacecraft cruising through this downstream plasma environment. This charge reduces electron fluxes at the spacecraft detectors to nondetectable intensities. Furthermore, we show that the Debye length $\lambda _{\mathrm D}^{\kappa}$ grows to values of about $\lambda _{\mathrm D}^{\kappa}/\lambda _{\mathrm D}\simeq 10^{6}$ compared to the classical value $\lambda _{\mathrm D}$ in this hot-electron environment. This unusual condition allows for the propagation of a certain type of electrostatic plasma waves that, at very large wavelengths, allow us to determine the effective temperature of the suprathermal electrons directly by means of the phase velocity of these waves. At moderate wavelengths, the electron-acoustic dispersion relation leads to nonpropagating oscillations with the ion-plasma frequency $\omega _{\mathrm p}$ , instead of the traditional electron plasma frequency.

Protons and alpha particles in the fast solar wind are only weakly collisional and exhibit a number of non-equilibrium features, including relative drifts between particle species. Two non-collisional mechanisms have been proposed for limiting differential flow between alpha particles and protons: plasma instabilities and the rotational force. Both mechanisms decelerate the alpha particles. In this paper, we derive an analytic expression for the rate $Q_{\mathrm{flow}}$ at which energy is released by alpha-particle deceleration, accounting for azimuthal flow and conservation of total momentum. We show that instabilities control the deceleration of alpha particles at $r< r_{\mathrm{crit}}$, and the rotational force controls the deceleration of alpha particles at $r> r_{\mathrm{crit}}$, where $r_{\mathrm{crit}} \simeq 2.5 \,\mathrm{AU}$ in the fast solar wind in the ecliptic plane. We find that $Q_{\mathrm{flow}}$ is positive at $r<r_{\mathrm{crit}}$ and $Q_{\mathrm{flow}} = 0$ at $r\geq r_{\mathrm{crit}}$, consistent with the previous finding that the rotational force does not lead to a release of energy. We compare the value of~$Q_{\mathrm{flow}}$ at $r< r_{\mathrm{crit}}$ with empirical heating rates for protons and alpha particles, denoted $Q_{\mathrm{p}}$ and $Q_{\alpha}$, deduced from in-situ measurements of fast-wind streams from the \emph{Helios} and \emph{Ulysses} spacecraft. We find that $Q_{\mathrm{flow}}$ exceeds $Q_{\alpha}$ at $r < 1\,\mathrm{AU}$, and that $Q_{\mathrm{flow}}/Q_{\rm p}$ decreases with increasing distance from the Sun from a value of about one at $r=0.29 - 0.42\,\mathrm{AU}$ to about 1/4 at 1 AU. We conclude that the continuous energy input from alpha-particle deceleration at $r< r_{\mathrm{crit}}$ makes an important contribution to the heating of the fast solar wind.

The collisionless dissipation of anisotropic Alfv\'enic turbulence is a promising candidate to solve the solar wind heating problem. Extensive studies examined the kinetic properties of Alfv\'en waves in simple Maxwellian or bi-Maxwellian plasmas. However, the observed electron velocity distribution functions in the solar wind are more complex. In this study, we analyze the properties of kinetic Alfv\'en waves in a plasma with two drifting electron populations. We numerically solve the linearized Maxwell-Vlasov equations and find that the damping rate and the proton-electron energy partition for kinetic Alfv\'en waves are significantly modified in such plasmas, compared to plasmas without electron drifts. We suggest that electron drift is an important factor to take into account when considering the dissipation of Alfv\'enic turbulence in the solar wind or other $\beta \sim 1$ astrophysical plasmas.

The Solar Probe Plus (SPP) spacecraft will explore the near-Sun environment, reaching heliocentric distances less than $10 R_{\odot}$. Near Earth, spacecraft measurements of fluctuating velocities and magnetic fields taken in the time domain are translated into information about the spatial structure of the solar wind via Taylor's "frozen turbulence" hypothesis. Near the perihelion of SPP, however, the solar-wind speed is comparable to the Alfv\'en speed, and Taylor's hypothesis in its usual form does not apply. In this paper, we show that, under certain assumptions, a modified version of Taylor's hypothesis can be recovered in the near-Sun region. We consider only the transverse, non-compressive component of the fluctuations at length scales exceeding the proton gyroradius, and we describe these fluctuations using an approximate theoretical framework developed by Heinemann and Olbert. We show that fluctuations propagating away from the Sun in the plasma frame obey a relation analogous to Taylor's hypothesis when $V_{\rm sc,\perp} \gg z^-$ and $z^+ \gg z^-$, where $V_{\rm sc,\perp}$ is the component of the spacecraft velocity perpendicular to the mean magnetic field and $\bm{z}^+$ ($\bm{z}^-$) is the Elsasser variable corresponding to transverse, non-compressive fluctuations propagating away from (towards) the Sun in the plasma frame. Observations and simulations suggest that, in the near-Sun solar wind, the above inequalities are satisfied and $\bm{z}^+$ fluctuations account for most of the fluctuation energy. The modified form of Taylor's hypothesis that we derive may thus make it possible to characterize the spatial structure of the energetically dominant component of the turbulence encountered by SPP.

We use particle-in-cell (PIC) simulations to study the nonlinear evolution of ion velocity space instabilities in an idealized problem in which a background velocity shear continuously amplifies the magnetic field. We simulate the astrophysically relevant regime where the shear timescale is long compared to the ion cyclotron period, and the plasma beta is ~ 1-100. The background field amplification in our calculation is meant to mimic processes such as turbulent fluctuations or MHD-scale instabilities. The field amplification continuously drives a pressure anisotropy with the perpendicular pressure larger than the parallel pressure, and the plasma becomes unstable to the mirror and ion cyclotron instabilities. In all cases, the nonlinear state is dominated by the mirror instability, not the ion cyclotron instability, and the plasma pressure anisotropy saturates near the threshold for the linear mirror instability. The magnetic field fluctuations initially undergo exponential growth but saturate in a secular phase in which the fluctuations grow on the same timescale as the background magnetic field (with \delta B ~ 0.3 <B> in the secular phase). At early times, the ion magnetic moment is well-conserved but once the fluctuation amplitudes exceed \delta B ~ 0.1 <B>, the magnetic moment is no longer conserved but instead changes on a timescale comparable to that of the mean magnetic field. We discuss the implications of our results for low-collisionality astrophysical plasmas, including the near-Earth solar wind and low-luminosity accretion disks around black holes.

Ambient plasma electrons undergo strong heating in regions associated with compressive traveling interplanetary solar-wind bulk-velocity jumps due to their specific interactions with the jump-inherent electric fields. After thermalization of this energy gain per shock passage through the operation of the Buneman instability, strong electron heating occurs that substantially influences the radial electron temperature profile. We describe the reduction of the jump amplitude due to energy expended by the traveling jump structure. We consider three effects; namely energy loss due to heating of electrons, energy loss due to work done against the pick-up-ion pressure gradient, and an energy gain due to nonlinear jump steepening. Taking these effects into account, we show that the decrease in jump amplitude with solar distance is more pronounced when the initial jump amplitude is higher in the inner solar system. Independent of the initial jump amplitude, it eventually decreases with increasing distance to a value of the order of $\Delta U/U\simeq 0.1$ at the position of the heliospheric termination shock, where $\Delta U$ is the jump amplitude, and $U$ is the average solar-wind bulk velocity.The electron temperature, on the other hand, is strongly correlated with the initial jump amplitude, leading to electron temperatures between 6000 K and 20 000 K at distances beyond 50 AU. We compare our results with in-situ measurements of the electron-core temperature from the Ulysses spacecraft in the plane of the ecliptic for $1.5\, \mathrm{AU}\leq r \leq 5\,\mathrm{AU}$, where $r$ is the distance from the Sun. We find a very good agreement between our results and these observations, which corroborates our extrapolated predictions beyond $r=5\,\mathrm{AU}$.

The MHD slow mode wave has application to coronal seismology, MHD turbulence, and the solar wind where it can be produced by parametric instabilities. We consider analytically how a drifting ion species (e.g. He$^{++}$) affects the linear slow mode wave in a mainly electron-proton plasma, with potential consequences for the aforementioned applications. Our main conclusions are: 1. For wavevectors highly oblique to the magnetic field, we find solutions that are characterized by very small perturbations of total pressure. Thus, our results may help to distinguish the MHD slow mode from kinetic Alfv\'en waves and non-propagating pressure-balanced structures, which can also have very small total pressure perturbations. 2. For small ion concentrations, there are solutions that are similar to the usual slow mode in an electron-proton plasma, and solutions that are dominated by the drifting ions, but for small drifts the wave modes cannot be simply characterized. 3. Even with zero ion drift, the standard dispersion relation for the highly oblique slow mode cannot be used with the Alfv\'en speed computed using the summed proton and ion densities, and with the sound speed computed from the summed pressures and densities of all species. 4. The ions can drive a non-resonant instability under certain circumstances. For low plasma beta, the threshold drift can be less than that required to destabilize electromagnetic modes, but damping from the Landau resonance can eliminate this instability altogether, unless $T_{\mathrm e}/T_{\mathrm p}\gg1$.

Previous studies have shown that the observed temperature anisotropies of protons and alpha particles in the solar wind are constrained by theoretical thresholds for pressure-anisotropy-driven instabilities such as the Alfv\'en/ion-cyclotron (A/IC) and fast-magnetosonic/whistler (FM/W) instabilities. In this letter, we use a long period of in-situ measurements provided by the {\em Wind} spacecraft's Faraday cups to investigate the combined constraint on the alpha-proton differential flow velocity and the alpha-particle temperature anisotropy due to A/IC and FM/W instabilities. We show that the majority of the data are constrained to lie within the region of parameter space in which A/IC and FM/W waves are either stable or have extremely low growth rates. In the minority of observed cases in which the growth rate of the A/IC (FM/W) instability is comparatively large, we find relatively higher values of $T_{\perp\alpha}/T_{\perp p}$ ($T_{\parallel\alpha}/T_{\parallel p}$) when alpha-proton differential flow velocity is small, where $T_{\perp\alpha}$ and $T_{\perp p}$ ($T_{\parallel\alpha}$ and $T_{\parallel p}$) are the perpendicular (parallel) temperatures of alpha particles and protons. We conjecture that this observed feature might arise from preferential alpha-particle heating which can drive the alpha particles beyond the instability thresholds.

We present a method for studying the evolution of plasma turbulence by tracking dispersion relations in the energy spectrum in the wavenumber-frequency domain. We apply hybrid plasma simulations in a simplified two-dimensional geometry to demonstrate our method and its applicability to plasma turbulence in the ion kinetic regime. We identify four dispersion relations: ion-Bernstein waves, oblique whistler waves, oblique Alfv\'en/ion-cyclotron waves, and a zero-frequency mode. The energy partition and frequency broadening are evaluated for these modes. The method allows us to determine the evolution of decaying plasma turbulence in our restricted geometry and shows that it cascades along the dispersion relations during the early phase with an increasing broadening around the dispersion relations.

We investigate the conditions under which parallel-propagating Alfv\'en/ion-cyclotron (A/IC) waves and fast-magnetosonic/whistler (FM/W) waves are driven unstable by the differential flow and temperature anisotropy of alpha particles in the solar wind. We focus on the limit in which $w_{\parallel \alpha} \gtrsim 0.25 v_{\mathrm A}$, where $w_{\parallel \alpha} $ is the parallel alpha-particle thermal speed and $v_{\mathrm A}$ is the Alfv\'en speed. We derive analytic expressions for the instability thresholds of these waves, which show, e.g., how the minimum unstable alpha-particle beam speed depends upon $w_{\parallel \alpha}/v_{\mathrm A}$, the degree of alpha-particle temperature anisotropy, and the alpha-to-proton temperature ratio. We validate our analytical results using numerical solutions to the full hot-plasma dispersion relation. Consistent with previous work, we find that temperature anisotropy allows A/IC waves and FM/W waves to become unstable at significantly lower values of the alpha-particle beam speed $U_\alpha$ than in the isotropic-temperature case. Likewise, differential flow lowers the minimum temperature anisotropy needed to excite A/IC or FM/W waves relative to the case in which $U_\alpha =0$. We discuss the relevance of our results to alpha particles in the solar wind near 1 AU.

We extend previous theories of stochastic ion heating to account for the motion of ions along the magnetic field. We derive an analytic expression for the ion-to-proton perpendicular temperature ratio in the solar wind for any ion species, assuming that stochastic heating is the dominant ion heating mechanism. This expression describes how this temperature ratio depends upon the average velocity of the ions along the magnetic field direction and the ratio of the parallel proton pressure to the magnetic pressure. We compare our model with previously published measurements of alpha particles and protons from the WIND spacecraft. We find that stochastic heating offers a promising explanation for these measurements when the fractional cross helicity and Alfven ratio at the proton-gyroradius scale have values that are broadly consistent with solar-wind measurements.

We investigate the conditions under which parallel-propagating Alfv\'en/ion-cyclotron waves are driven unstable by an isotropic ($T_{\perp \alpha} = T_{\parallel\alpha}$) population of alpha particles drifting parallel to the magnetic field at an average speed $U_{\alpha}$ with respect to the protons. We derive an approximate analytic condition for the minimum value of $U_{\alpha}$ needed to excite this instability and refine this result using numerical solutions to the hot-plasma dispersion relation. When the alpha-particle number density is $\simeq 5%$ of the proton number density and the two species have similar thermal speeds, the instability requires that $\beta_{\rm p} \gtrsim 1$, where $\beta_{\rm p}$ is the ratio of the proton pressure to the magnetic pressure. For $1\lesssim \beta_{\mathrm p}\lesssim 12$, the minimum $U_{\alpha}$ needed to excite this instability ranges from $0.7v_{\mathrm A}$ to $0.9v_{\mathrm A}$, where $v_{\mathrm A}$ is the Alfv\'en speed. This threshold is smaller than the threshold of $\simeq 1.2v_{\mathrm A}$ for the parallel magnetosonic instability, which was previously thought to have the lowest threshold of the alpha-particle beam instabilities at $\beta_{\mathrm p}\gtrsim 0.5$. We discuss the role of the parallel Alfv\'enic drift instability for the evolution of the alpha-particle drift speed in the solar wind. We also analyze measurements from the \emph{Wind} spacecraft's Faraday cups and show that the $U_{\alpha}$ values measured in solar-wind streams with $T_{\perp \alpha}\approx T_{\parallel\alpha}$ are approximately bounded from above by the threshold of the parallel Alfv\'enic instability.

An ion beam can destabilize Alfv\'en/ion-cyclotron waves and magnetosonic/whistler waves if the beam speed is sufficiently large. Numerical solutions of the hot-plasma dispersion relation have previously shown that the minimum beam speed required to excite such instabilities is significantly smaller for oblique modes with $\vec k \times \vec B_0\neq 0$ than for parallel-propagating modes with $\vec k \times \vec B_0 = 0$, where $\vec k$ is the wavevector and $\vec B_0$ is the background magnetic field. In this paper, we explain this difference within the framework of quasilinear theory, focusing on low-$\beta$ plasmas. We begin by deriving, in the cold-plasma approximation, the dispersion relation and polarization properties of both oblique and parallel-propagating waves in the presence of an ion beam. We then show how the instability thresholds of the different wave branches can be deduced from the wave--particle resonance condition, the conservation of particle energy in the wave frame, the sign (positive or negative) of the wave energy, and the wave polarization. We also provide a graphical description of the different conditions under which Landau resonance and cyclotron resonance destabilize Alfv\'en/ion-cyclotron waves in the presence of an ion beam. We draw upon our results to discuss the types of instabilities that may limit the differential flow of alpha particles in the solar wind.

Certain types of plasma waves are known to become parametrically unstable under specific plasma conditions, in which the pump wave will decay into several daughter waves with different wavenumbers and frequencies. In the past, the related plasma instabilities have been treated analytically for various parameter regimes and by use of various numerical methods, yet the oblique propagation with respect to the background magnetic field has rarely been dealt with in two dimensions, mainly because of the high computational demand. Here we present a hybrid-simulation study of the parametric decay of a moderately oblique Alfv\'en wave having elliptical polarization. It is found that such a compressive wave can decay into waves with higher and lower wavenumbers than the pump.

The nature of solar wind turbulence in the dissipation range at scales much smaller than the large MHD scales remains under debate. Here a two-dimensional model based on the hybrid code abbreviated as A.I.K.E.F. is presented, which treats massive ions as particles obeying the kinetic Vlasov equation and massless electrons as a neutralizing fluid. Up to a certain wavenumber in the MHD regime, the numerical system is initialized by assuming a superposition of isotropic Alfv\'en waves with amplitudes that follow the empirically confirmed spectral law of Kolmogorov. Then turbulence develops and energy cascades into the dispersive spectral range, where also dissipative effects occur. Under typical solar wind conditions, weak turbulence develops as a superposition of normal modes in the kinetic regime. Spectral analysis in the direction parallel to the background magnetic field reveals a cascade of left-handed Alfv\'en/ion-cyclotron waves up to wave vectors where their resonant absorption sets in, as well as a continuing cascade of right-handed fast-mode and whistler waves. Perpendicular to the background field, a broad turbulent spectrum is found to be built up of fluctuations having a strong compressive component. Ion-Bernstein waves seem to be possible normal modes in this propagation direction for lower driving amplitudes. Also signatures of short-scale pressure-balanced structures (very oblique slow-mode waves) are found.

In this paper, we study the weakly-compressive high-frequency plasma waves which are superposed on a large-amplitude Alfv\'en wave in a multi-fluid plasma consisting of protons, electrons, and alpha particles. For these waves, the plasma environment is inhomogenous due to the presence of the low-frequency Alfv\'en wave with a large amplitude, a situation that may apply to space plasmas such as the solar corona and solar wind. The dispersion relation of the plasma waves is determined from a linear stability analysis using a new eigenvalue method that is employed to solve the set of differential wave equations which describe the propagation of plasma waves along the direction of the constant component of the Alfv\'en wave magnetic field. This approach also allows one to consider weak compressive effects. In the presence of the background Alfv\'en wave, the dispersion branches obtained differ significantly from the situation of a uniform plasma. Due to compressibility, acoustic waves are excited and couplings between various modes occur, and even an instability of the compressive mode. In a kinetic treatment, these plasma waves would be natural candidates for Landau-resonant wave-particle interactions, and may thus via their damping lead to particle heating.

Large-amplitude Alfv\'en waves are ubiquitous in space plasmas and a main component of magnetohydrodynamic (MHD) turbulence in the heliosphere. As pump waves they are prone to parametric instability by which they can generate cyclotron and acoustic daughter waves. Here we revisit a related process within the framework of the multi-fluid equations for a plasma consisting of many species. The nonlinear coupling of the Alfv\'en wave to acoustic waves is studied, and a set of compressive and coupled wave equations for the transverse magnetic field and longitudinal electric field is derived for waves propagating along the mean-field direction. It turns out that slightly compressive Alfv\'en waves exert, through induced gyro-radius and kinetic-energy modulations, an electromotive force on the particles in association with a longitudinal electric field, which has a potential that is given by the gradient of the transverse kinetic energy of the particles gyrating about the mean field. This in turn drives electric fluctuations (sound and ion-acoustic waves) along the mean magnetic field, which can nonlinearly react back on the transverse magnetic field. Mutually coupled Alfv\'en-cyclotron-acoustic waves are thus excited, a nonlinear process that can drive a cascade of wave energy in the plasma and may generate compressive microturbulence. These driven electric fluctuations might have consequences for the dissipation of MHD turbulence and, thus, for the heating and acceleration of particles in the solar wind.

The fast solar wind is a collisionless plasma permeated by plasma waves on many different scales. A plasma wave represents the natural interplay between the periodic changes of the electromagnetic field and the associated coherent motions of the plasma particles. In this paper, a model velocity distribution function is derived for a plasma in a single, coherent, large-amplitude wave. This model allows one to study the kinetic effects of wave motions on particle distributions. They are by in-situ spacecraft measured by counting, over a certain sampling time, the particles coming from various directions and having different energies. We compare our results with the measurements by the Helios spacecraft, and thus find that by assuming high wave activity we are able to explain key observed features of the measured distributions within the framework of our model. We also address the recent discussions on nonresonant wave--particle interactions and apparent heating. The applied time-averaging procedure leads to an apparent ion temperature anisotropy which is connected but not identical to the intrinsic temperature of the underlying distribution function.

Though pick-up ions (PUIs) are a well known phenomenon in the inner heliosphere their phase-space distribution nevertheless is a theoretically unsettled problem. Especially the question of how pick-up ions form their suprathermal tails, extending to far above their injection energies, still now is unsatistactorily answered. Though Fermi-2 velocity diffusion theories have revealed that such tails are populated, they nevertheless show that resulting population densities are much less than seen in observations showing power-laws with a velocity index of "-5". We first investigate here, whether or not observationally suggested power-laws can be the result of a quasi-equilibrium state between suprathermal ions and magnetohydrodynamic turbulences in energy exchange with eachother. We demonstrate that such an equilibrium cannot be established, since it would require too high pick-up ion pressures enforcing a shock-free deceleration of the solar wind. We furthermore show that Fermi-2 type energy diffusion in the outer heliosphere is too inefficient to determine the shape of the distribution function there. As we can show, however, power-laws beyond the injection threshold can be established, if the injection takes place at higher energies of the order of 100 keV. As we demonstrate here, such an injection is connected with modulated anomalous cosmic ray (ACR) particles at the lower end of their spectrum when they again start being convected outwards with the solar wind. Therefore, we refer to these particles as ACR-PUIs. In our quantitative calculation of the pick-up ion spectrum resulting under such conditions we in fact find again power-laws, however with a velocity power index of "-4" and fairly distance-independent spectral intensities.