Recent near-sun solar-wind observations from Parker Solar Probe have found a highly dynamic magnetic environment, permeated by abrupt radial-field reversals, or "switchbacks." We show that many features of the observed turbulence are reproduced by a spectrum of Alfv\'enic fluctuations advected by a radially expanding flow. Starting from simple superpositions of low-amplitude outward-propagating waves, our expanding-box compressible MHD simulations naturally develop switchbacks because (i) the normalized amplitude of waves grows due to expansion and (ii) fluctuations evolve towards spherical polarization (i.e., nearly constant field strength). These results suggest that switchbacks form in-situ in the expanding solar wind and are not indicative of impulsive processes in the chromosphere or corona.

We perform a statistical study of the turbulent power spectrum at inertial and kinetic scales observed during the first perihelion encounter of Parker Solar Probe. We find that often there is an extremely steep scaling range of the power spectrum just above the ion-kinetic scales, similar to prior observations at 1 AU, with a power-law index of around $-4$. Based on our measurements, we demonstrate that either a significant ($>50\%$) fraction of the total turbulent energy flux is dissipated in this range of scales, or the characteristic nonlinear interaction time of the turbulence decreases dramatically from the expectation based solely on the dispersive nature of nonlinearly interacting kinetic Alfv\'en waves.

We consider 2D joint distributions of normalised residual energy $\sigma_r(s,t)$ and cross helicity $\sigma_c(s,t)$ during one day of Parker Solar Probe's (PSP's) first encounter as a function of wavelet scale $s$. The broad features of the distributions are similar to previous observations made by HELIOS in slow solar wind, namely well correlated and fairly Alfv\'enic, except for a population with negative cross helicity which is seen at shorter wavelet scales. We show that this population is due to the presence of magnetic switchbacks, brief periods where the magnetic field polarity reverses. Such switchbacks have been observed before, both in HELIOS data and in Ulysses data in the polar solar wind. Their abundance and short timescales as seen by PSP in its first encounter is a new observation, and their precise origin is still unknown. By analysing these MHD invariants as a function of wavelet scale we show that MHD waves do indeed follow the local mean magnetic field through switchbacks, with net Elsasser flux propagating inward during the field reversal, and that they therefore must be local kinks in the magnetic field and not due to small regions of opposite polarity on the surface of the Sun. Such observations are important to keep in mind as computing cross helicity without taking into account the effect of switchbacks may result in spurious underestimation of $\sigma_c$ as PSP gets closer to the Sun in later orbits.

Stochastic heating is a non-linear heating mechanism driven by the violation of magnetic moment invariance due to large-amplitude turbulent fluctuations producing diffusion of ions towards higher kinetic energies in the direction perpendicular to the magnetic field. It is frequently invoked as a mechanism responsible for the heating of ions in the solar wind. Here, we quantify for the first time the proton stochastic heating rate $Q_\perp$ at radial distances from the Sun as close as $0.16$ au, using measurements from the first two Parker Solar Probe encounters. Our results for both the amplitude and radial trend of the heating rate, $Q_\perp \propto r^{-2.5}$, agree with previous results based on the Helios data set at heliocentric distances from 0.3 to 0.9 au. Also in agreement with previous results, $Q_\perp$ is significantly larger in the fast solar wind than in the slow solar wind. We identify the tendency in fast solar wind for cuts of the core proton velocity distribution transverse to the magnetic field to exhibit a flat-top shape. The observed distribution agrees with previous theoretical predictions for fast solar wind where stochastic heating is the dominant heating mechanism.

The first two orbits of the Parker Solar Probe (PSP) spacecraft have enabled the first in situ measurements of the solar wind down to a heliocentric distance of 0.17 au (or 36 Rs). Here, we present an analysis of this data to study solar wind turbulence at 0.17 au and its evolution out to 1 au. While many features remain similar, key differences at 0.17 au include: increased turbulence energy levels by more than an order of magnitude, a magnetic field spectral index of -3/2 matching that of the velocity and both Elsasser fields, a lower magnetic compressibility consistent with a smaller slow-mode kinetic energy fraction, and a much smaller outer scale that has had time for substantial nonlinear processing. There is also an overall increase in the dominance of outward-propagating Alfv\'enic fluctuations compared to inward-propagating ones, and the radial variation of the inward component is consistent with its generation by reflection from the large-scale gradient in Alfv\'en speed. The energy flux in this turbulence at 0.17 au was found to be ~10% of that in the bulk solar wind kinetic energy, becoming ~40% when extrapolated to the Alfv\'en point, and both the fraction and rate of increase of this flux towards the Sun is consistent with turbulence-driven models in which the solar wind is powered by this flux.

Understanding the physical processes in the solar wind and corona which actively contribute to heating, acceleration, and dissipation is a primary objective of NASA's Parker Solar Probe (PSP) mission. Observations of coherent electromagnetic waves at ion scales suggests that linear cyclotron resonance and non-linear processes are dynamically relevant in the inner heliosphere. A wavelet-based statistical study of coherent waves in the first perihelion encounter of PSP demonstrates the presence of transverse electromagnetic waves at ion resonant scales which are observed in 30-50\% of radial field intervals. Average wave amplitudes of approximately 4 nT are measured, while the mean duration of wave events is of order 20 seconds; however long duration wave events can exist without interruption on hour-long timescales. Though ion scale waves are preferentially observed during intervals with a radial mean magnetic field, we show that measurement constraints, associated with single spacecraft sampling of quasi-parallel waves superposed with anisotropic turbulence, render the measured quasi-parallel ion-wave spectrum unobservable when the mean magnetic field is oblique to the solar wind flow; these results imply that the occurrence of coherent ion-scale waves is not limited to a radial field configuration. The lack of strong radial scaling of characteristic wave amplitudes and duration suggests that the waves are generated {\em{in-situ}} through plasma instabilities. Additionally, observations of proton distribution functions indicate that temperature anisotropy may drive the observed ion-scale waves.

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

We present 3D numerical simulations and an analytic model of reflection-driven MHD turbulence in the solar wind. Our simulations describe transverse, non-compressive MHD fluctuations within a narrow magnetic flux tube that extends from the photosphere out to a heliocentric distance $r$ of 21 solar radii $(R_s)$. We launch outward-propagating "$z^+$ fluctuations" into the simulation domain by imposing a randomly evolving photospheric velocity field. As these fluctuations propagate away from the Sun, they undergo partial reflection, producing inward-propagating "$z^-$ fluctuations." Counter-propagating fluctuations subsequently interact, causing fluctuation energy to cascade to small scales and dissipate. Our analytic model incorporates alignment, allows for strongly or weakly turbulent nonlinear interactions, and divides the $z^+$ fluctuations into two populations with different characteristic radial correlation lengths. The inertial-range power spectra in our simulations evolve toward a $k_\perp^{-3/2}$ scaling at $r>10 R_s$, where $k_\perp$ is the wave-vector component perpendicular to the background magnetic field. In two of our simulations, the $z^+$ power spectra are much flatter between the coronal base and $r \simeq 4 R_s$. We argue that these spectral scalings are caused by: (1) high-pass filtering in the upper chromosphere; (2) the anomalous coherence of inertial-range $z^-$ fluctuations in a reference frame propagating outwards with the $z^+$ fluctuations; and (3) the change in the sign of the radial derivative of the Alfv\'en speed at $r=r_m \simeq 1.7 R_s$, which disrupts this anomalous coherence between $r=r_m$ and $r\simeq 2r_m$. At $r>1.3 R_s$, the turbulent heating rate in our simulations is comparable to the heating rate in a previously developed solar-wind model that agreed with a number of observational constraints.

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.

We present three-dimensional, hybrid-kinetic numerical simulations of driven Alfv\'{e}n-wave turbulence of relevance to the collisionless near-Earth solar wind. Special attention is paid to the spectral transition that occurs near the ion-Larmor scale and to the origins of preferential perpendicular ion heating and of non-thermal wings in the parallel distribution function. Several novel diagnostics are used to show that the ion heating rate increases as the kinetic-Alfv\'{e}n-wave fluctuations, which comprise the majority of the sub-ion-Larmor turbulent cascade, attain near-ion-cyclotron frequencies. We find that ${\approx}75$-$80\%$ of the cascade energy goes into heating the ions, broadly consistent with the near-Earth solar wind. This heating is accompanied by clear velocity-space signatures in the particle energization rates and the distribution functions, including a flattened core in the perpendicular-velocity distribution and non-Maxwellian wings in the parallel-velocity distribution. The latter are attributed to transit-time damping and the pitch-angle scattering of perpendicularly heated particles into the parallel direction. Accompanying these features is a steepening of the spectral index of sub-ion-Larmor magnetic-field fluctuations beyond the canonical $-2.8$, as field energy is transferred to thermal energy. These predictions may be tested by measurements in the near-Earth solar wind.

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.

Turbulence in an accretion disk launches Alfv\'en waves (AWs) that propagate away from the disk along magnetic field lines. Because the Alfv\'en speed varies with distance from the disk, the AWs undergo partial non-WKB reflection, and counter-propagating AWs subsequently interact, causing AW energy to cascade to small scales and dissipate. To investigate this process, we introduce an Elsasser-like formulation of general relativistic MHD (GRMHD) and develop the theory of general relativistic reduced MHD in an inhomogeneous medium. We then derive a set of equations for the mean-square AW amplitude $M_+$ and turbulent heating rate $Q$ under the assumption that, in the plasma rest frame, AWs propagating away from the disk are much more energetic than AWs propagating toward the disk. For the case in which the background flow is axisymmetric and time-independent, we solve these equations analytically to determine $M_+$ and $Q$ as functions of position. We find that, for an idealized thin disk threaded by a large-scale poloidal magnetic field, the AW energy flux is $\sim (\rho_{\rm b}/\rho_{\rm d})^{1/2} \beta_{\rm net,d}^{-1/2}$ times the disk's radiative flux, where $\rho_{\rm b}$ and $\rho_{\rm d}$ are the mass densities at the coronal base and disk midplane, respectively, and $\beta_{\rm net,d}$ is the ratio (evaluated at the disk midplane) of plasma-plus-radiation pressure to the pressure of the average vertical magnetic field. This energy flux could have a significant impact on disk coronae and outflows. To lay the groundwork for future global simulations of turbulent disk coronae and jets, we derive a set of averaged GRMHD equations that account for reflection-driven AW turbulence using a sub-grid model.

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.

In this paper, weak turbulence theory is used to investigate the nonlinear evolution of the parametric instability in 3D low-$\beta$ plasmas at wavelengths much greater than the ion inertial length under the assumption that slow magnetosonic waves are strongly damped. It is shown analytically that the parametric instability leads to an inverse cascade of Alfv\'en wave quanta, and several exact solutions to the wave kinetic equations are presented. The main results of the paper concern the parametric decay of Alfv\'en waves that initially satisfy $e^+ \gg e^-$, where $e^+$ and $e^-$ are the frequency ($\,f$) spectra of Alfv\'en waves propagating in opposite directions along the magnetic field lines. If $e^+$ initially has a peak frequency $f_0$ (at which $f e^+$ is maximized) and an "infrared" scaling $f^p$ at smaller $f$ with $-1 < p < 1$, then $e^+$ acquires an $f^{-1}$ scaling throughout a range of frequencies that spreads out in both directions from $f_0$. At the same time, $e^-$ acquires an $f^{-2}$ scaling within this same frequency range. If the plasma parameters and infrared $e^+$ spectrum are chosen to match conditions in the fast solar wind at a heliocentric distance of 0.3 astronomical units (AU), then the nonlinear evolution of the parametric instability leads to an $e^+$ spectrum that matches fast-wind measurements from the Helios spacecraft at 0.3 AU, including the observed $f^{-1}$ scaling at $f \gtrsim 3 \times 10^{-4} \mbox{ Hz}$. The results of this paper suggest that the $f^{-1}$ spectrum seen by Helios in the fast solar wind at $f \gtrsim 3\times 10^{-4} \mbox{ Hz}$ is produced in situ by parametric decay and that the $f^{-1}$ range of $e^+$ extends over an increasingly narrow range of frequencies as $r$ decreases below 0.3 AU. This prediction will be tested by measurements from the Parker Solar Probe.

We calculate the disruption scale $\lambda_{\rm D}$ at which sheet-like structures in dynamically aligned Alfv\'enic turbulence are destroyed by the onset of magnetic reconnection in a low-$\beta$ collisionless plasma. The scaling of $\lambda_{\rm D}$ depends on the order of the statistics being considered, with more intense structures being disrupted at larger scales. The disruption scale for the structures that dominate the energy spectrum is $\lambda_{\rm D}\sim L_\perp^{1/9}(d_e\rho_s)^{4/9}$, where $d_e$ is the electron inertial scale, $\rho_s$ is the ion sound scale, and $L_\perp$ is the outer scale of the turbulence. When $\beta_e$ and $\rho_s/L_\perp$ are sufficiently small, the scale $\lambda_{\rm D}$ is larger than $\rho_s$ and there is a break in the energy spectrum at $\lambda_{\rm D}$, rather than at $\rho_s$. We propose that the fluctuations produced by the disruption are circularised flux ropes, which may have already been observed in the solar wind. We predict the relationship between the amplitude and radius of these structures and quantify the importance of the disruption process to the cascade in terms of the filling fraction of undisrupted structures and the fractional reduction of the energy contained in them at the ion sound scale $\rho_s$. Both of these fractions depend strongly on $\beta_e$, with the disrupted structures becoming more important at lower $\beta_e$. Finally, we predict that the energy spectrum between $\lambda_{\rm D}$ and $\rho_s$ is steeper than $k_\perp^{-3}$, when this range exists. Such a steep "transition range" is sometimes observed in short intervals of solar-wind turbulence. The onset of collisionless magnetic reconnection may therefore significantly affect the nature of plasma turbulence around the ion gyroscale.

We propose a mechanism whereby the intense, sheet-like structures naturally formed by dynamically aligning Alfv\'enic turbulence are destroyed by magnetic reconnection at a scale $\hat{\lambda}_{\rm D}$, larger than the dissipation scale predicted by models of intermittent, dynamically aligning turbulence. The reconnection process proceeds in several stages: first, a linear tearing mode with $N$ magnetic islands grows and saturates, and then the $X$-points between these islands collapse into secondary current sheets, which then reconnect until the original structure is destroyed. This effectively imposes an upper limit on the anisotropy of the structures within the perpendicular plane, which means that at scale $\hat{\lambda}_{\rm D}$ the turbulent dynamics change: at scales larger than $\hat{\lambda}_{\rm D}$, the turbulence exhibits scale-dependent dynamic alignment and a spectral index approximately equal to $-3/2$, while at scales smaller than $\hat{\lambda}_{\rm D}$, the turbulent structures undergo a succession of disruptions due to reconnection, limiting dynamic alignment, steepening the effective spectral index and changing the final dissipation scale. The scaling of $\hat{\lambda}_{\rm D}$ with the Lundquist (magnetic Reynolds) number $S_{L_\perp}$ depends on the order of the statistics being considered, and on the specific model of intermittency; the transition between the two regimes in the energy spectrum is predicted at approximately $\hat{\lambda}_{\rm D} \sim S_{L_\perp}^{-0.6}$. The spectral index below $\hat{\lambda}_{\rm D}$ is bounded between $-5/3$ and $-2.3$. The final dissipation scale is at $\hat{\lambda}_{\eta,\infty}\sim S_{L_\perp}^{-3/4}$, the same as the Kolmogorov scale arising in theories of turbulence that do not involve scale-dependent dynamic alignment.

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.

We investigate how the proton distribution function evolves when the protons undergo stochastic heating by strong, low-frequency, Alfv\'en-wave turbulence under the assumption that $\beta$ is small. We apply our analysis to protons undergoing stochastic heating in the supersonic fast solar wind and obtain proton distributions at heliocentric distances ranging from 4 to 30 solar radii. We find that the proton distribution develops non-Gaussian structure with a flat core and steep tail. For $r >5 \ R_{\rm S}$, the proton distribution is well approximated by a modified Moyal distribution. Comparisons with future measurements from \emph{Solar Probe Plus} could be used to test whether stochastic heating is occurring in the solar-wind acceleration region.

We measure the local anisotropy of numerically simulated strong Alfv\'enic turbulence with respect to two local, physically relevant directions: along the local mean magnetic field and along the local direction of one of the fluctuating Elsasser fields. We find significant scaling anisotropy with respect to both these directions: the fluctuations are "ribbon-like" --- statistically, they are elongated along both the mean magnetic field and the fluctuating field. The latter form of anisotropy is due to scale-dependent alignment of the fluctuating fields. The intermittent scalings of the $n$th-order conditional structure functions in the direction perpendicular to both the local mean field and the fluctuations agree well with the theory of Chandran et al. 2015, while the parallel scalings are consistent with those implied by the critical-balance conjecture. We quantify the relationship between the perpendicular scalings and those in the fluctuation and parallel directions, and find that the scaling exponent of the perpendicular anisotropy (i.e., of the aspect ratio of the Alfv\'enic structures in the plane perpendicular to the mean magnetic field) depends on the amplitude of the fluctuations. This is shown to be equivalent to the anticorrelation of fluctuation amplitude and alignment at each scale. The dependence of the anisotropy on amplitude is shown to be more significant for the anisotropy between the perpendicular and fluctuation-direction scales than it is between the perpendicular and parallel scales.

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.

We present numerical evidence that in strong Alfvenic turbulence, the critical balance principle---equality of the nonlinear decorrelation and linear propagation times---is scale invariant, in the sense that the probability distribution of the ratio of these times is independent of scale. This result only holds if the local alignment of the Elsasser fields is taken into account in calculating the nonlinear time. At any given scale, the degree of alignment is found to increase with fluctuation amplitude, supporting the idea that the cause of alignment is mutual dynamical shearing of Elsasser fields. The scale-invariance of critical balance (while all other quantities of interest are strongly intermittent, i.e., have scale-dependent distributions) suggests that it is the most robust of the scaling principles used to describe Alfvenic turbulence. The quality afforded by situ fluctuation measurements in the solar wind allows for direct verification of this fundamental principle.

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 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 analytic estimates and numerical simulations of test particles interacting with magnetohydrodynamic (MHD) turbulence to show that subsonic MHD turbulence produces efficient second-order Fermi acceleration of relativistic particles. This acceleration is not well-described by standard quasi-linear theory but is a consequence of resonance broadening of wave-particle interactions in MHD turbulence. We provide momentum diffusion coefficients that can be used for astrophysical and heliospheric applications and discuss the implications of our results for accretion flows onto black holes. In particular, we show that particle acceleration by subsonic turbulence in radiatively inefficient accretion flows can produce a non-thermal tail in the electron distribution function that is likely important for modeling and interpreting the emission from low luminosity systems such as Sgr A* and M87.

We develop a model for stochastic acceleration of electrons in solar flares. As in several previous models, the electrons are accelerated by turbulent fast magnetosonic waves ("fast waves") via transit-time-damping (TTD) interactions. (In TTD interactions, fast waves act like moving magnetic mirrors that push the electrons parallel or anti-parallel to the magnetic field). We also include the effects of Coulomb collisions and the waves' parallel electric fields. Unlike previous models, our model is two-dimensional in both momentum space and wavenumber space and takes into account the anisotropy of the wave power spectrum $F_k$ and electron distribution function $f_{\rm e}$. We use weak turbulence theory and quasilinear theory to obtain a set of equations that describes the coupled evolution of $F_k$ and $f_{\rm e}$. We solve these equations numerically and find that the electron distribution function develops a power-law-like non-thermal tail within a restricted range of energies $E\in (E_{\rm nt}, E_{\rm max})$. We obtain approximate analytic expressions for $E_{\rm nt}$ and $E_{\rm max}$, which describe how these minimum and maximum energies depend upon parameters such as the electron number density and the rate at which fast-wave energy is injected into the acceleration region at large scales. We contrast our results with previous studies that assume that $F_k$ and $f_{\rm e}$ are isotropic, and we compare one of our numerical calculations with the time-dependent hard-x-ray spectrum observed during the June 27, 1980 flare. In our numerical calculations, the electron energy spectra are softer (steeper) than in models with isotropic $F_k$ and $f_{\rm e}$ and closer to the values inferred from observations of solar flares.

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

We develop an analytic model of intermittent, three-dimensional, strong, reduced magnetohydrodynamic turbulence with zero cross helicity. We take the fluctuation amplitudes to have a log-Poisson distribution and incorporate into the model a new phenomenology of scale-dependent dynamic alignment. The log-Poisson distribution in our model is characterized by two parameters. To calculate these parameters, we make use of three assumptions: that the energy cascade rate is independent of scale within the inertial range, that the most intense coherent structures at scale $\lambda$ are sheet-like with a volume filling factor proportional to $\lambda$, and that most of the cascade power arises from interactions between exceptionally intense fluctuations and much weaker fluctuations. We then compute the scalings of the power spectrum, the kurtosis, higher-order structure functions, and three different average alignment angles. These scalings appear to be consistent with numerous results from direct numerical simulations.

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 direct numerical simulations of inhomogeneous reduced magnetohydrodynamic (RMHD) turbulence between the Sun and the Alfv\'en critical point. These are the first such simulations that take into account the solar-wind outflow velocity and the radial inhomogeneity of the background solar wind without approximating the nonlinear terms in the governing equations. RMHD turbulence is driven by outward-propagating Alfv\'en waves ($z^+$ fluctuations) launched from the Sun, which undergo partial non-WKB reflection to produce sunward-propagating Alfv\'en waves ($z^-$ fluctuations). We present ten simulations with different values of the correlation time $\tau_{\rm c\,\sun}^+$ and perpendicular correlation length $L_{\perp \sun}$ of outward-propagating Alfv\'en waves (AWs) at the coronal base. We find that between 15% and 33% of the $z^+$ energy launched into the corona dissipates between the coronal base and Alfv\'en critical point. Between 33% and 40% of this input energy goes into work on the solar-wind outflow, and between 22% and 36% escapes as $z^+$ fluctuations through the simulation boundary at $r=r_{\rm A}$. The $z^\pm$ power spectra scale like $k_\perp^{-\alpha^\pm}$, where $k_\perp$ is the wavenumber in the plane perpendicular to $\vec{B}_0$. In our simulation with the smallest value of $\tau_{\rm c\,\sun}^+$ ($\sim 2 \mbox{min}$) and largest value of $L_{\perp \sun}$ ($2\times 10^4 \mbox{km}$), we find that $\alpha^+$ decreases approximately linearly with increasing $\ln(r)$, reaching a value of 1.3 at $r=11.1 R_{\sun}$. Our simulations with larger values of $\tau_{\rm c\,\sun}^+$ exhibit alignment between the contours of constant $\phi^+$, $\phi^-$, $\Omega_0^+$, and $\Omega_0^-$, where $\phi^\pm$ are the Els\"asser potentials and $\Omega_0^\pm$ are the outer-scale parallel Els\"asser vorticities.

Recent theoretical studies argue that the rate of stochastic ion heating in low-frequency Alfv\'en-wave turbulence is given by $Q_\perp = c_1 [(\delta u)^3 /\rho] \exp(-c_2/\epsilon)$, where $\delta u$ is the rms turbulent velocity at the scale of the ion gyroradius $\rho$, $\epsilon = \delta u/v_{\perp \rm i}$, $v_{\perp \rm i}$ is the perpendicular ion thermal speed, and $c_1$ and $c_2$ are dimensionless constants. We test this theoretical result by numerically simulating test particles interacting with strong reduced magnetohydrodynamic (RMHD) turbulence. The heating rates in our simulations are well fit by this formula. The best-fit values of $c_1$ are $\sim 1$. The best-fit values of $c_2$ decrease (i.e., stochastic heating becomes more effective) as the grid size and Reynolds number of the RMHD simulations increase. As an example, in a $1024^2 \times 256$ RMHD simulation with a dissipation wavenumber of order the inverse ion gyroradius, we find $c_2 = 0.21$. We show that stochastic heating is significantly stronger in strong RMHD turbulence than in a field of randomly phased Alfv\'en waves with the same power spectrum, because coherent structures in strong RMHD turbulence increase orbit stochasticity in the regions where ions are heated most strongly. We find that $c_1$ increases by a factor of $\sim 3$ while $c_2$ changes very little as the ion thermal speed increases from values $\ll v_{\rm A}$ to values $\sim v_{\rm A}$, where $v_{\rm A}$ is the Alfv\'en speed. We discuss the importance of these results for perpendicular ion heating in the solar wind.

We investigate the effects of pitch-angle scattering on the efficiency of particle heating and acceleration by MHD turbulence using phenomenological estimates and simulations of non-relativistic test particles interacting with strong, subsonic MHD turbulence. We include an imposed pitch-angle scattering rate, which is meant to approximate the effects of high frequency plasma waves and/or velocity space instabilities. We focus on plasma parameters similar to those found in the near-Earth solar wind, though most of our results are more broadly applicable. An important control parameter is the size of the particle mean free path lambda_{mfp} relative to the scale of the turbulent fluctuations L. For small scattering rates, particles interact quasi-resonantly with turbulent fluctuations in magnetic field strength. Scattering increases the long-term efficiency of this resonant heating by factors of a few-10, but the distribution function does not develop a significant non-thermal power-law tail. For higher scattering rates, the interaction between particles and turbulent fluctuations becomes non-resonant, governed by particles heating and cooling adiabatically as they encounter turbulent density fluctuations. Rapid pitch-angle scattering can produce a power-law tail in the proton distribution function but this requires fine-tuning of parameters. Moreover, in the near-Earth solar wind, a significant power-law tail cannot develop by this mechanism because the particle acceleration timescales are longer than the adiabatic cooling timescale set by the expansion of the solar wind. Our results thus imply that MHD-scale turbulent fluctuations are unlikely to be the origin of the v^{-5} tail in the proton distribution function observed in the solar wind.

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.

Spacecraft measurements show that protons undergo substantial perpendicular heating during their transit from the Sun to the outer heliosphere. In this paper, we use {\em Helios~2} measurements to investigate whether stochastic heating by low-frequency turbulence is capable of explaining this perpendicular heating. We analyze {\em Helios~2} magnetic-field measurements in low-$\beta$ fast-solar-wind streams between heliocentric distances $r=0.29$ AU and $r=0.64$ AU to determine the rms amplitude of the fluctuating magnetic field, $\delta B_{\rm p}$, near the proton gyroradius scale $\rho_{\rm p}$. We then evaluate the stochastic heating rate $Q_{\perp \rm stoch}$ using the measured value of $\delta B_{\rm p}$ and a previously published analytical formula for $Q_{\perp \rm stoch}$. Using {\em Helios} measurements we estimate the `empirical' perpendicular heating rate $Q_{\perp \rm emp} = (k_{\rm B}/m_{\rm p}) B V (d/dr) (T_{\perp \rm p}/B)$ that is needed to explain the $T_{\perp \rm p}$ profile. We find that $Q_{\perp \rm stoch} \sim Q_{\perp \rm emp}$, but only if a key dimensionless constant appearing in the formula for $Q_{\perp \rm stoch}$ lies within a certain range of values. This range is approximately the same throughout the radial interval that we analyze and is consistent with the results of numerical simulations of the stochastic heating of test particles in reduced magnetohydrodynamic turbulence. These results support the hypothesis that stochastic heating accounts for much of the perpendicular proton heating occurring in low-$\beta$ fast-wind streams.

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.

The heating, acceleration, and pitch-angle scattering of charged particles by MHD turbulence are important in a wide range of astrophysical environments, including the solar wind, accreting black holes, and galaxy clusters. We simulate the interaction of high-gyrofrequency test particles with fully dynamical simulations of subsonic MHD turbulence, focusing on the parameter regime with beta ~ 1, where beta is the ratio of gas to magnetic pressure. We use the simulation results to calibrate analytical expressions for test particle velocity-space diffusion coefficients and provide simple fits that can be used in other work. The test particle velocity diffusion in our simulations is due to a combination of two processes: interactions between particles and magnetic compressions in the turbulence (as in linear transit-time damping; TTD) and what we refer to as Fermi Type-B (FTB) interactions, in which charged particles moving on field lines may be thought of as beads spiralling around moving wires. We show that test particle heating rates are consistent with a TTD resonance which is broadened according to a decorrelation prescription that is Gaussian in time. TTD dominates the heating for v_s >> v_A (e.g. electrons), where v_s is the thermal speed of species s and v_A is the Alfven speed, while FTB dominates for v_s << v_A (e.g. minor ions). Proton heating rates for beta ~ 1 are comparable to the turbulent cascade rate. Finally, we show that velocity diffusion of collisionless, large gyrofrequency particles due to large-scale MHD turbulence does not produce a power-law distribution function.

We develop a 1D solar-wind model that includes separate energy equations for the electrons and protons, proton temperature anisotropy, collisional and collisionless heat flux, and an analytical treatment of low-frequency, reflection-driven, Alfven-wave turbulence. To partition the turbulent heating between electron heating, parallel proton heating, and perpendicular proton heating, we employ results from the theories of linear wave damping and nonlinear stochastic heating. We account for mirror and oblique firehose instabilities by increasing the proton pitch-angle scattering rate when the proton temperature anisotropy exceeds the threshold for either instability. We numerically integrate the equations of the model forward in time until a steady state is reached, focusing on two fast-solar-wind-like solutions. These solutions are consistent with a number of observations, supporting the idea that Alfven-wave turbulence plays an important role in the origin of the solar wind.

In February 2011 we proposed a NASA Explorer Mission of Opportunity program to develop and operate a large-aperture ultraviolet coronagraph spectrometer called the Coronal Physics Investigator (CPI) as an attached International Space Station (ISS) payload. The primary goal of this program is to identify and characterize the physical processes that heat and accelerate the primary and secondary components of the fast and slow solar wind. In addition, CPI can make key measurements needed to understand CMEs. UVCS/SOHO allowed us to identify what additional measurements need to be made to answer the fundamental questions about how solar wind streams are produced, and CPI's next-generation capabilities were designed specifically to make those measurements. Compared to previous instruments, CPI provides unprecedented sensitivity, a wavelength range extending from 25.7 to 126 nm, higher temporal resolution, and the capability to measure line profiles of He II, N V, Ne VII, Ne VIII, Si VIII, S IX, Ar VIII, Ca IX, and Fe X, never before seen in coronal holes above 1.3 solar radii. CPI will constrain the properties and effects of coronal MHD waves by (1) observing many ions over a large range of charge and mass, (2) providing simultaneous measurements of proton and electron temperatures to probe turbulent dissipation mechanisms, and (3) measuring amplitudes of low-frequency compressive fluctuations. CPI is an internally occulted ultraviolet coronagraph that provides the required high sensitivity without the need for a deployable boom, and with all technically mature hardware including an ICCD detector. A highly experienced Explorer and ISS contractor, L-3 Com Integrated Optical Systems and Com Systems East, will provide the tracking and pointing system as well as the instrument, and the integration to the ISS.

Understanding the physical processes responsible for accelerating the solar wind requires detailed measurements of the collisionless plasma in the extended solar corona. Some key clues about these processes have come from instruments that combine the power of an ultraviolet (UV) spectrometer with an occulted telescope. This combination enables measurements of ion emission lines far from the bright solar disk, where most of the solar wind acceleration occurs. Although the UVCS instrument on SOHO made several key discoveries, many questions remain unanswered because its capabilities were limited. This white paper summarizes these past achievements and also describes what can be accomplished with next-generation instrumentation of this kind.

Low-frequency Alfven-wave turbulence causes ion trajectories to become chaotic, or "stochastic," when the turbulence amplitude is sufficiently large. Stochastic orbits enable ions to absorb energy from the turbulence, increasing the perpendicular ion temperature even when the fluctuation frequencies are too small for a cyclotron resonance to occur. In this paper, an analytic expression for the stochastic heating rate is used in conjunction with an observationally constrained turbulence model to obtain an analytic formula for the perpendicular ion temperature as a function of heliocentric distance r, ion mass, and ion charge in coronal holes for values of r between 2 and 15 solar radii (Rs). The resulting temperature profiles provide a good fit to observations of protons and OVI ions at 2Rs < r < 3Rs from the Ultraviolet Coronagraph Spectrometer (UVCS). Stochastic heating also offers a natural explanation for several detailed features of the UVCS observations, including the preferential and anisotropic heating of minor ions, the rapid radial increase in the OVI temperature between 1.6Rs and 1.9Rs, and the abrupt flattening of the OVI temperature profile as r increases above 1.9Rs.

We consider ion heating by turbulent Alfven waves (AWs) and kinetic Alfven waves (KAWs) with perpendicular wavelengths comparable to the ion gyroradius and frequencies smaller than the ion cyclotron frequency. When the turbulence amplitude exceeds a certain threshold, an ion's orbit becomes chaotic. The ion then interacts stochastically with the time-varying electrostatic potential, and the ion's energy undergoes a random walk. Using phenomenological arguments, we derive an analytic expression for the rates at which different ion species are heated, which we test by simulating test particles interacting with a spectrum of randomly phased AWs and KAWs. We find that the stochastic heating rate depends sensitively on the quantity epsilon = dv/vperp, where vperp is the component of the ion velocity perpendicular to the background magnetic field B0, and dv (dB) is the rms amplitude of the velocity (magnetic-field) fluctuations at the gyroradius scale. In the case of thermal protons, when epsilon << eps1, where eps1 is a constant, a proton's magnetic moment is nearly conserved and stochastic heating is extremely weak. However, when epsilon > eps1, the proton heating rate exceeds the cascade power that would be present in strong balanced KAW turbulence with the same value of dv, and magnetic-moment conservation is violated. For the random-phase waves in our test-particle simulations, eps1 is approximately 0.2. For protons in low-beta plasmas, epsilon is approximately dB/B0 divided by the square root of beta, and epsilon can exceed eps1 even when dB/B0 << eps1. At comparable temperatures, alpha particles and minor ions have larger values of epsilon than protons and are heated more efficiently as a result. We discuss the implications of our results for ion heating in coronal holes and the solar wind.

Resonant interactions between ions and Alfv\'en/ion-cyclotron (A/IC) waves may play an important role in the heating and acceleration of the fast solar wind. Although such interactions have been studied extensively for "parallel" waves, whose wave vectors ${\bf k}$ are aligned with the background magnetic field ${\bf B}_0$, much less is known about interactions between ions and oblique A/IC waves, for which the angle $\theta$ between ${\bf k}$ and ${\bf B}_0$ is nonzero. In this paper, we present new numerical results on resonant cyclotron interactions between protons and oblique A/IC waves in collisionless low-beta plasmas such as the solar corona. We find that if some mechanism generates oblique high-frequency A/IC waves, then these waves initially modify the proton distribution function in such a way that it becomes unstable to parallel waves. Parallel waves are then amplified to the point that they dominate the wave energy at the large parallel wave numbers at which the waves resonate with the particles. Pitch-angle scattering by these waves then causes the plasma to evolve towards a state in which the proton distribution is constant along a particular set of nested "scattering surfaces" in velocity space, whose shapes have been calculated previously. As the distribution function approaches this state, the imaginary part of the frequency of parallel A/IC waves drops continuously towards zero, but oblique waves continue to undergo cyclotron damping while simultaneously causing protons to diffuse across these kinetic shells to higher energies. We conclude that oblique A/IC waves can be more effective at heating protons than parallel A/IC waves, because for oblique waves the plasma does not relax towards a state in which proton damping of oblique A/IC waves ceases.

We study the propagation, reflection, and turbulent dissipation of Alfven waves in coronal holes and the solar wind. We start with the Heinemann-Olbert equations, which describe non-compressive magnetohydrodynamic fluctuations in an inhomogeneous medium with a background flow parallel to the background magnetic field. Following the approach of Dmitruk et al, we model the nonlinear terms in these equations using a simple phenomenology for the cascade and dissipation of wave energy, and assume that there is much more energy in waves propagating away from the Sun than waves propagating towards the Sun. We then solve the equations analytically for waves with periods of hours and longer to obtain expressions for the wave amplitudes and turbulent heating rate as a function of heliocentric distance. We also develop a second approximate model that includes waves with periods of roughly one minute to one hour, which undergo less reflection than the longer-period waves, and compare our models to observations. Our models generalize the phenomenological model of Dmitruk et al by accounting for the solar wind velocity, so that the turbulent heating rate can be evaluated from the coronal base out past the Alfven critical point - that is, throughout the region in which most of the heating and acceleration occurs. The simple analytical expressions that we obtain can be used to incorporate Alfven-wave reflection and turbulent heating into fluid models of the solar wind.

The intracluster medium (ICM) is stably stratified in the hydrodynamic sense with the entropy $s$ increasing outwards. However, thermal conduction along magnetic field lines fundamentally changes the stability of the ICM, leading to the "heat-flux buoyancy instability" when $dT/dr>0$ and the "magnetothermal instability" when $dT/dr<0$. The ICM is thus buoyantly unstable regardless of the signs of $dT/dr$ and $ds/dr$. On the other hand, these temperature-gradient-driven instabilities saturate by reorienting the magnetic field (perpendicular to $\hat{\bf r}$ when $dT/dr>0$ and parallel to $\hat{\bf r}$ when $dT/dr<0$), without generating sustained convection. We show that after an anisotropically conducting plasma reaches this nonlinearly stable magnetic configuration, it experiences a buoyant restoring force that resists further distortions of the magnetic field. This restoring force is analogous to the buoyant restoring force experienced by a stably stratified adiabatic plasma. We argue that in order for a driving mechanism (e.g, galaxy motions or cosmic-ray buoyancy) to overcome this restoring force and generate turbulence in the ICM, the strength of the driving must exceed a threshold, corresponding to turbulent velocities $\gtrsim 10 -100 {km/s}$. For weaker driving, the ICM remains in its nonlinearly stable magnetic configuration, and turbulent mixing is effectively absent. We discuss the implications of these findings for the turbulent diffusion of metals and heat in the ICM.

One proposed mechanism for heating the solar wind, from close to the sun to beyond 10 AU, invokes low-frequency, oblique, Alfven-wave turbulence. Because small-scale oblique Alfven waves (kinetic Alfven waves) are compressive, the measured density fluctuations in the solar wind place an upper limit on the amplitude of kinetic Alfven waves and hence an upper limit on the rate at which the solar wind can be heated by low-frequency, Alfvenic turbulence. We evaluate this upper limit for both coronal holes at 5 solar radii and in the near-Earth solar wind. At both radii, the upper limit we find is consistent with models in which the solar wind is heated by low-frequency Alfvenic turbulence. At 1 AU, the upper limit on the turbulent heating rate derived from the measured density fluctuations is within a factor of 2 of the measured solar wind heating rate. Thus if low-frequency Alfvenic turbulence contributes to heating the near-Earth solar wind, kinetic Alfven waves must be one of the dominant sources of solar wind density fluctuations at frequencies of order 1 Hz. We also present a simple argument for why density fluctuation measurements do appear to rule out models in which the solar wind is heated by non-turbulent high-frequency waves ``sweeping'' through the ion-cyclotron resonance, but are compatible with heating by low-frequency Alfvenic turbulence.

Different results for the cascade power in strong, incompressible MHD turbulence with nonzero cross helicity appear in the literature. In this paper, we discuss the conditions under which these different results are valid. We define z+ to be the rms amplitude of Alfven waves propagating parallel to the background magnetic field, and z- to be the rms amplitude of Alfven waves propagating anti-parallel to the background magnetic field. Nonzero cross helicity implies that z+ and z- differ, and we take z- to be less than z+. We find that the mechanism that generates the z- fluctuations strongly affects the cascade power, because it controls the coherence time for interactions between oppositely directed wave packets at the outer scale. In particular, for fixed values of z+ and z-, the cascade power is in many cases larger when the z- fluctuations are generated by the reflection of z+ fluctuations than when the z- fluctuations are generated by forcing that is only weakly correlated with the z+ fluctuations.

(Abridged abstract) A theory of incompressible MHD turbulence recently developed by Boldyrev predicts the existence of a scale dependent angle of alignment between velocity and magnetic field fluctuations that is proportional to the lengthscale of the fluctuations to the power 1/4. In this study, plasma and magnetic field data from the Wind spacecraft are used to investigate the angle between velocity and magnetic field fluctuations in the solar wind as a function of the timescale of the fluctuations and to look for the power law scaling predicted by Boldyrev.

Using a linear stability analysis and two and three-dimensional nonlinear simulations, we study the physics of buoyancy instabilities in a combined thermal and relativistic (cosmic ray) plasma, motivated by the application to clusters of galaxies. We argue that cosmic ray diffusion is likely to be slow compared to the buoyancy time on large length scales, so that cosmic rays are effectively adiabatic. If the cosmic ray pressure $p_{cr}$ is $\gtrsim 25 %$ of the thermal pressure, and the cosmic ray entropy ($p_{\rm cr}/\rho^{4/3}$; $\rho$ is the thermal plasma density) decreases outwards, cosmic rays drive an adiabatic convective instability analogous to Schwarzschild convection in stars. Global simulations of galaxy cluster cores show that this instability saturates by reducing the cosmic ray entropy gradient and driving efficient convection and turbulent mixing. At larger radii in cluster cores, the thermal plasma is unstable to the heat flux-driven buoyancy instability (HBI), a convective instability generated by anisotropic thermal conduction and a background conductive heat flux. Cosmic-ray driven convection and the HBI may contribute to redistributing metals produced by Type 1a supernovae in clusters. Our calculations demonstrate that adiabatic simulations of galaxy clusters can artificially suppress the mixing of thermal and relativistic plasma; anisotropic thermal conduction allows more efficient mixing, which may contribute to cosmic rays being distributed throughout the cluster volume.

In this Letter, weak turbulence theory is used to investigate interactions among Alfven waves and fast and slow magnetosonic waves in collisionless low-beta plasmas. The wave kinetic equations are derived from the equations of magnetohydrodynamics, and extra terms are then added to model collisionless damping. These equations are used to provide a quantitative description of a variety of nonlinear processes, including "parallel" and "perpendicular" energy cascade, energy transfer between wave types, "phase mixing," and the generation of back-scattered Alfven waves.

We report the results of a local stability analysis for a magnetized, gravitationally stratified plasma containing cosmic rays. We account for cosmic-ray diffusion and thermal conduction parallel to the magnetic field and allow beta to take any value, where p is the plasma pressure and B is the magnetic field strength. We take the gravitational acceleration to be in the -z-direction and the equilibrium magnetic field to be in the y-direction, and we derive the dispersion relation for small-amplitude instabilities and waves in the large-|k_x| limit. We use the Routh-Hurwitz criterion to show analytically that the necessary and sufficient criterion for stability in this limit is n k_B dT/dz + dp_cr/dz + (1/8pi)dB^2/dz > 0, where T is the temperature, n is the number density of thermal particles, and p_cr is the cosmic-ray pressure. We present approximate analytical solutions for the normal modes in the low- and high-diffusivity limits, show that they are consistent with the derived stability criterion, and compare them to numerical results obtained from the full, unapproximated, dispersion relation. Our results extend earlier analyses of buoyancy instabilities in galaxy-cluster plasmas to the beta <= 1 regime. Our results also extend earlier analyses of the Parker instability to account for anisotropic thermal conduction, and show that the interstellar medium is more unstable to the Parker instability than was predicted by previous studies in which the thermal plasma was treated as adiabatic.

In clusters of galaxies, the specific entropy of intracluster plasma increases outwards. Nevertheless, a number of recent studies have shown that the intracluster medium is subject to buoyancy instabilities due to the effects of cosmic rays and anisotropic thermal conduction. In this paper, we present a new numerical algorithm for simulating such instabilities. This numerical method treats the cosmic rays as a fluid, accounts for the diffusion of heat and cosmic rays along magnetic field lines, and enforces the condition that the temperature and cosmic-ray pressure remain positive. We carry out several tests to ensure the accuracy of the code, including the detailed matching of analytic results for the eigenfunctions and growth rates of linear buoyancy instabilities. This numerical scheme will be useful for simulating convection driven by cosmic-ray buoyancy in galaxy cluster plasmas and may also be useful for other applications, including fusion plasmas, the interstellar medium, and supernovae remnants.

Central peaks in the iron abundance of intracluster plasma are a common feature of cooling-core galaxy clusters. These abundance peaks have a much broader profile than the stars of the central brightest cluster galaxy (BCG), which produce the excess iron, indicating that metal-enriched plasma is transported out of the BCG by some process such as turbulent diffusion. At the same time, cooling-core clusters are likely heated by central active galactic nuclei (AGNs) by means of jets, cosmic-ray bubbles, and convection. The recent AGN-driven convection model of Chandran&Rasera predicts the turbulent velocity profile in a steady-state cluster in which radiative cooling is balanced by heating from a combination of AGN-driven convection and thermal conduction. We use the velocity profiles as input into an advection/diffusion model for the transport of metals in the intracluster medium, taking the iron to be injected by the BCG. We compare the results of our model to XMM and Chandra observations of eight clusters. Assuming a constant turbulence level over a cluster's lifetime, the turbulent velocities in the model can explain the observed abundance profiles in only five of the eight clusters. However, we go on to develop an analytic fit of the turbulent velocity profile as a function of the AGN power. We then deduce for each cluster the average AGN power (during the past 10 Gyr) required to match the abundance profiles. The required average values are between 10^43 and 2.10^44 erg/s, while the present AGN powers span a much larger range from 6.10^41 (Virgo) to 2.10^44 erg/s (Hydra A). Our results suggest that AGN-driven convection can account for the observed abundance profiles if the AGN power varies over a cluster's lifetime between Virgo-like and Hydra-A-like values, with average values in the above-quoted range.

This paper proposes a new phenomenology for strong incompressible MHD turbulence with nonzero cross helicity. This phenomenology is then developed into a quantitative Fokker-Planck model that describes the time evolution of the anisotropic power spectra of the fluctuations propagating parallel and anti-parallel to the background magnetic field. It is found that in steady state the power spectra of the magnetic field and total energy are steeper than a Kolmogorov spectrum and become increasingly steep as C/E increases, where C is the cross helicity and E is the fluctuation energy. Increasing C with fixed E increases the time required for energy to cascade to smaller scales, reduces the cascade power, and increases the anisotropy of the small-scale fluctuations. The implications of these results for the solar wind and solar corona are discussed in some detail.

A number of studies have shown that the convective stability criterion for the intracluster medium (ICM) is very different from the Schwarzchild criterion due to the effects of anisotropic thermal conduction and cosmic rays. Building on these studies, we develop a model of the ICM in which a central active galactic nucleus (AGN) accretes hot intracluster plasma at the Bondi rate and produces cosmic rays that cause the ICM to become convectively unstable. The resulting convection heats the intracluster plasma and regulates its temperature and density profiles. By adjusting a single parameter in the model (the size of the cosmic-ray acceleration region), we are able to achieve a good match to the observed density and temperature profiles in a sample of eight clusters. Our results suggest that convection is an important process in cluster cores. An interesting feature of our solutions is that the cooling rate is more sharply peaked about the cluster center than is the convective heating rate. As a result, in several of the clusters in our sample, a compact cooling flow arises in the central region with a size R that is typically a few kpc. The cooling flow matches onto a Bondi flow at smaller radii. The mass accretion rate in the Bondi flow is equal to, and controlled by, the rate at which mass flows in through the cooling flow. Our solutions suggest that the AGN regulates the mass accretion rate in these clusters by controlling R: if the AGN power rises above the equilibrium level, R decreases, the mass accretion rate drops, and the AGN power drops back down to the equilibrium level.

This Letter presents a calculation of the power spectra of weakly turbulent Alfven waves and fast magnetosonic waves ("fast waves") in low-beta plasmas. It is shown that three-wave interactions transfer energy to high-frequency fast waves and, to a lesser extent, high-frequency Alfven waves. MHD turbulence is thus a promising mechanism for producing the high-frequency waves needed to explain the anisotropic heating of minor ions in the solar corona.

This paper describes how active galactic nuclei can heat galaxy-cluster plasmas by driving convection in the intracluster medium. A model is proposed in which a central supermassive black hole accretes intracluster plasma at the Bondi rate and powers a radio source. The central radio source produces cosmic rays, which mix into the thermal plasma. The cosmic-ray luminosity is proportional to the mass accretion rate. The cosmic-ray pressure gradient drives convection, which causes plasma heating. It is assumed that plasma heating balances radiative cooling. The plasma heating rate is self-regulating because the Bondi accretion rate is a decreasing function of the specific entropy near the cluster center; if heating exceeds cooling, then the central entropy increases, the Bondi rate and cosmic-ray luminosity decrease, and the convective heating rate is reduced. This paper focuses on the role of intracluster magnetic fields, which affect convection by causing heat and cosmic rays to diffuse primarily along magnetic field lines. A new stability criterion is derived for convection in a thermal-plasma/cosmic-ray fluid, and equations for the average fluid properties in a convective cluster are obtained with the use of a nonlocal two-fluid mixing length theory. Numerical solutions of the model equations compare reasonably well with observations without requiring fine tuning of the model parameters.

Turbulent heating may play an important role in galaxy-cluster plasmas, but if turbulent heating is to balance radiative cooling in a quasi-steady state, some mechanism must set the turbulent velocity to the required value. This paper explores one possible regulating mechanism associated with an active galactic nucleus at cluster center. A steady-state model for the intracluster medium is presented in which radiative cooling is balanced by a combination of turbulent heating and thermal conduction. The turbulence is generated by convection driven by the buoyancy of cosmic rays produced by a central radio source. The cosmic-ray luminosity is powered by the accretion of intracluster plasma onto a central black hole. The model makes the rather extreme assumption that the cosmic rays and thermal plasma are completely mixed. Although the intracluster medium is convectively unstable near cluster center in the model solutions, the specific entropy of the thermal plasma still increases outwards because of the cosmic-ray modification to the stability criterion. The model provides a self-consistent calculation of the turbulent velocity as a function of position, but fails to reproduce the steep central density profiles observed in clusters. The principal difficulty is that in order for the fully mixed intracluster medium to become convectively unstable, the cosmic-ray pressure must become comparable to or greater than the thermal pressure within the convective region. The large cosmic-ray pressure gradient then provides much of the support against gravity, reducing the thermal pressure gradient near cluster center and decreasing the central plasma density gradient. A more realistic AGN-feedback model of intracluster turbulence in which relativistic and thermal plasmas are only partially mixed may have greater success.

We study relativistic particles undergoing surfing acceleration at perpendicular shocks. We assume that particles undergo diffusion in the component of momentum perpendicular to the shock plane due to moderate fluctuations in the shock electric and magnetic fields. We show that dN/dE, the number of surfing-accelerated particles per unit energy, attains a power-law form, dN/dE \propto E^{-b}. We calculate b analytically in the limit of weak momentum diffusion, and use Monte Carlo test-particle calculations to evaluate b in the weak, moderate, and strong momentum-diffusion limits.

This brief invited review provides an introduction to key ideas in the theory of incompressible magnetohydrodynamic (MHD) turbulence.

Energetic particles that undergo strong pitch-angle scattering and diffuse through a plasma containing strong compressible MHD turbulence undergo diffusion in momentum space with diffusion coefficient Dp. In this paper, the contribution of slow modes to Dp is calculated assuming the rms turbulent velocity is of order the Alfven speed. The energy spectrum of accelerated particles is derived assuming slow modes make the dominant contribution to Dp, taking into account Coulomb losses and particle escape from the acceleration region with an energy-independent escape time. The results are applied to solar flares.

Fast particles diffusing along magnetic field lines in a turbulent plasma can diffuse through and then return to the same eddy many times before the eddy is randomized in the turbulent flow. This leads to an enhancement of particle acceleration by large-scale compressible turbulence relative to previous estimates in which isotropic particle diffusion is assumed.

We investigate field-line separation in strong MHD turbulence using direct numerical simulations. We find that in the static-magnetic-field approximation the thermal conductivity in galaxy clusters is reduced by a factor of about 50 relative to the Spitzer thermal conductivity of a non-magnetized plasma. This value is too small for heat conduction to balance radiative cooling in clusters.

We investigate field-line separation in strong MHD turbulence analytically and with direct numerical simulations. We find that in the static-magnetic-field approximation the thermal conductivity in galaxy clusters is reduced by a factor of about 5-10 relative to the Spitzer thermal conductivity of a non-magnetized plasma. We also estimate how the thermal conductivity would be affected by efficient turbulent resistivity.

Observations of narrow radio-emitting filaments near the Galactic center have been interpreted in previous studies as evidence of a pervasive vertical (i.e. perpendicular to the Galactic plane) milliGauss magnetic field in the central 150 pc of the Galaxy. A simple cylindrically symmetric model for the equilibrium in this central region is proposed in which horizontal (i.e. parallel to the Galactic plane) magnetic fields embedded in an annular band of partially ionized molecular material of radius 150 pc are wrapped around vertical magnetic fields threading low-density hot plasma. The central vertical magnetic field, which has a pressure that significantly exceeds the thermal pressure of the medium, is confined by the weight of the molecular material. The stability of this equilibrium is studied indirectly by analyzing a uniformly rotating cylinder of infinite extent along the z axis in cylindrical coordinates (r,theta,z), with low-density plasma and an axial magnetic field at r< 150 pc, high-density plasma and an azimuthal field at r> 150 pc, and a gravitational acceleration g* proportional to r directed in the negative-r-hat direction. The density profile and gravity tend to destabilize the plasma, but the plasma tends to be stabilized by rotation and magnetic tension--since the interface between the high and low-density plasmas can not be perturbed without bending either the horizontal or vertical field. It is shown analytically that when beta= 8(pi)p/B^2 is small and the dense plasma is supported against gravity primarily by rotation, the necessary and sufficient condition for stability to k_z=0 modes is |g| < (2|Omega| a), where g = g* - Omega^2 r is the effective gravity, Omega is the uniform angular velocity, and "a" is the sound speed in the dense plasma.

We study the propagation of radio waves through a medium possessing density fluctuations that are elongated along the ambient magnetic field and described by an anisotropic Goldreich-Sridhar power spectrum. We derive general formulas for the wave phase structure function, visibility, angular broadening, diffraction-pattern length scales, and scintillation time scale for arbitrary distributions of turbulence along the line of sight, and specialize these formulas to idealized cases.

Scattering rates for a Goldreich-Sridhar (GS) spectrum of anisotropic, incompressible, magnetohydrodynamic turbulence are calculated in the quasilinear approximation. Because the small-scale fluctuations are constrained to have wave vectors nearly perpendicular to the background magnetic field, scattering is too weak to provide either the mean free paths commonly used in Galactic cosmic-ray propagation models or the mean free paths required for acceleration of cosmic rays at quasi-parallel shocks. Where strong pitch-angle scattering occurs, it is due to fluctuations not described by the GS spectrum, such as fluctuations generated by streaming cosmic rays.

Recent studies suggest that when magnetohydrodynamic (MHD) turbulence is excited by stirring a plasma at large scales, the cascade of energy from large to small scales is anisotropic, in the sense that small-scale fluctuations satisfy the inequality k_parallel << k_perp, where k_parallel and k_perp are, respectively, the components of a fluctuation's wave vector parallel and perpendicular to the background magnetic field. Such anisotropic fluctuations are very inefficient at scattering cosmic rays. Results based on the quasilinear approximation for scattering of cosmic rays by anisotropic MHD turbulence are presented and explained. The important role played by molecular-cloud magnetic mirrors in confining and isotropizing cosmic rays when scattering is weak is also discussed.