Source author record

Benjamin D. G. Chandran

Benjamin D. G. Chandran appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2022arXiv

Electron-driven instabilities in the solar wind

The electrons are an essential particle species in the solar wind. They often exhibit non-equilibrium features in their velocity distribution function. These include temperature anisotropies, tails (kurtosis), and reflectional asymmetries (skewness), which contribute a significant heat flux to the solar wind. If these non-equilibrium features are sufficiently strong, they drive kinetic micro-instabilities. We develop a semi-graphical framework based on the equations of quasi-linear theory to describe electron-driven instabilities in the solar wind. We apply our framework to resonant instabilities driven by temperature anisotropies. These include the electron whistler anisotropy instability and the propagating electron firehose instability. We then describe resonant instabilities driven by reflectional asymmetries in the electron distribution function. These include the electron/ion-acoustic, kinetic Alfvén heat-flux, Langmuir, electron-beam, electron/ion-cyclotron, electron/electron-acoustic, whistler heat-flux, oblique fast-magnetosonic/whistler, lower-hybrid fan, and electron-deficit whistler instability. We briefly comment on non-resonant instabilities driven by electron temperature anisotropies such as the mirror-mode and the non-propagating firehose instability. We conclude our review with a list of open research topics in the field of electron-driven instabilities in the solar wind.

preprint2021arXiv

Exact nonlinear solutions for three-dimensional Alfvén-wave packets in relativistic magnetohydrodynamics

We show that large-amplitude, non-planar, Alfvén wave (AW) packets are exact nonlinear solutions of the relativistic MHD equations when the total magnetic-field strength in the local fluid rest frame ($b$) is a constant. We derive analytic expressions relating the components of the fluctuating velocity and magnetic field. We also show that these constant-$b$ AWs propagate without distortion at the relativistic Alfvén velocity and never steepen into shocks. These findings and the observed abundance of large-amplitude, constant-$b$ AWs in the solar wind suggest that such waves may be present in relativistic outflows around compact astrophysical objects.

preprint2020arXiv

In-situ switchback formation in the expanding solar wind

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énic 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.

preprint2020arXiv

Inner-Heliosphere Signatures of Ion-Scale Dissipation and Nonlinear Interaction

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én waves.

preprint2019arXiv

Cross Helicity Reversals In Magnetic Switchbacks

We consider 2D joint distributions of normalised residual energy $σ_r(s,t)$ and cross helicity $σ_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énic, 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 $σ_c$ as PSP gets closer to the Sun in later orbits.

preprint2019arXiv

The Enhancement of Proton Stochastic Heating in the near-Sun Solar Wind

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.

preprint2016arXiv

Collisionless Isotropization of the Solar-Wind Protons by Compressive Fluctuations and Plasma Instabilities

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 $β_{\parallel \mathrm p}\equiv 8π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.

preprint2016arXiv

Evolution of The Proton Velocity Distribution due to Stochastic Heating in the Near-Sun Solar Wind

We investigate how the proton distribution function evolves when the protons undergo stochastic heating by strong, low-frequency, Alfvén-wave turbulence under the assumption that $β$ 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.

preprint2015arXiv

A Modified Version of Taylor's Hypothesis for Solar Probe Plus Observations

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én 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.

preprint2015arXiv

Deceleration of Alpha Particles in the Solar Wind by Instabilities and the Rotational Force: Implications for Heating, Azimuthal Flow, and the Parker Spiral Magnetic Field

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_α$, 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_α$ 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.

preprint2015arXiv

On the Conservation of Cross Helicity and Wave Action in Solar-Wind Models with Non-WKB Alfven Wave Reflection

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.

preprint2014arXiv

Acceleration of Relativistic Electrons by MHD Turbulence: Implications for Non-thermal Emission from Black Hole Accretion Disks

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.

preprint2014arXiv

Intermittency and Alignment in Strong RMHD Turbulence

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 $λ$ are sheet-like with a volume filling factor proportional to $λ$, 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.

preprint2014arXiv

Magnetohydrodynamic Slow Mode with Drifting He$^{++}$: Implications for Coronal Seismology and the Solar Wind

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

preprint2014arXiv

Stochastic Acceleration of Electrons by Fast Magnetosonic Waves in Solar Flares: the Effects of Anisotropy in Velocity andWavenumber Space

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.

preprint2013arXiv

A Parallel-Propagating Alfvénic Ion-Beam Instability in the High-Beta Solar Wind

We investigate the conditions under which parallel-propagating Alfvén/ion-cyclotron waves are driven unstable by an isotropic ($T_{\perp α} = T_{\parallelα}$) population of alpha particles drifting parallel to the magnetic field at an average speed $U_α$ with respect to the protons. We derive an approximate analytic condition for the minimum value of $U_α$ 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 $β_{\rm p} \gtrsim 1$, where $β_{\rm p}$ is the ratio of the proton pressure to the magnetic pressure. For $1\lesssim β_{\mathrm p}\lesssim 12$, the minimum $U_α$ needed to excite this instability ranges from $0.7v_{\mathrm A}$ to $0.9v_{\mathrm A}$, where $v_{\mathrm A}$ is the Alfvén 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 $β_{\mathrm p}\gtrsim 0.5$. We discuss the role of the parallel Alfvénic 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_α$ values measured in solar-wind streams with $T_{\perp α}\approx T_{\parallelα}$ are approximately bounded from above by the threshold of the parallel Alfvénic instability.

preprint2013arXiv

Direct Numerical Simulations of Reflection-Driven, Reduced MHD Turbulence from the Sun to the Alfven Critical Point

We present direct numerical simulations of inhomogeneous reduced magnetohydrodynamic (RMHD) turbulence between the Sun and the Alfvén 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én waves ($z^+$ fluctuations) launched from the Sun, which undergo partial non-WKB reflection to produce sunward-propagating Alfvén waves ($z^-$ fluctuations). We present ten simulations with different values of the correlation time $τ_{\rm c\,\sun}^+$ and perpendicular correlation length $L_{\perp \sun}$ of outward-propagating Alfvén 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én 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^{-α^\pm}$, where $k_\perp$ is the wavenumber in the plane perpendicular to $\vec{B}_0$. In our simulation with the smallest value of $τ_{\rm c\,\sun}^+$ ($\sim 2 \mbox{min}$) and largest value of $L_{\perp \sun}$ ($2\times 10^4 \mbox{km}$), we find that $α^+$ 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 $τ_{\rm c\,\sun}^+$ exhibit alignment between the contours of constant $ϕ^+$, $ϕ^-$, $Ω_0^+$, and $Ω_0^-$, where $ϕ^\pm$ are the Elsässer potentials and $Ω_0^\pm$ are the outer-scale parallel Elsässer vorticities.

preprint2013arXiv

Instabilities Driven by the Drift and Temperature Anisotropy of Alpha Particles in the Solar Wind

We investigate the conditions under which parallel-propagating Alfvén/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 α} \gtrsim 0.25 v_{\mathrm A}$, where $w_{\parallel α} $ is the parallel alpha-particle thermal speed and $v_{\mathrm A}$ is the Alfvén 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 α}/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_α$ 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_α=0$. We discuss the relevance of our results to alpha particles in the solar wind near 1 AU.

preprint2013arXiv

Limits on Alpha Particle Temperature Anisotropy and Differential Flow from Kinetic Instabilities: Solar Wind Observations

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én/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α}/T_{\perp p}$ ($T_{\parallelα}/T_{\parallel p}$) when alpha-proton differential flow velocity is small, where $T_{\perpα}$ and $T_{\perp p}$ ($T_{\parallelα}$ 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.

preprint2013arXiv

Observational Test of Stochastic Heating in Low-$β$ Fast Solar Wind Streams

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-$β$ 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, $δB_{\rm p}$, near the proton gyroradius scale $ρ_{\rm p}$. We then evaluate the stochastic heating rate $Q_{\perp \rm stoch}$ using the measured value of $δ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-$β$ fast-wind streams.

preprint2013arXiv

Perpendicular Ion Heating by Reduced Magnetohydrodynamic Turbulence

Recent theoretical studies argue that the rate of stochastic ion heating in low-frequency Alfvén-wave turbulence is given by $Q_\perp = c_1 [(δu)^3 /ρ] \exp(-c_2/ε)$, where $δu$ is the rms turbulent velocity at the scale of the ion gyroradius $ρ$, $ε= δ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én 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én speed. We discuss the importance of these results for perpendicular ion heating in the solar wind.

preprint2013arXiv

The Dispersion Relations and Instability Thresholds of Oblique Plasma Modes in the Presence of an Ion Beam

An ion beam can destabilize Alfvén/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-$β$ 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én/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.

preprint2013arXiv

The Efficiency of Second-Order Fermi Acceleration by Weakly Compressible MHD Turbulence

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.

preprint2012arXiv

Resonance Broadening and Heating of Charged Particles in Magnetohydrodynamic Turbulence

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.

preprint2011arXiv

Incorporating Kinetic Physics into a Two-Fluid Solar-Wind Model with Temperature Anisotropy and Low-Frequency Alfven-Wave Turbulence

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.

preprint2010arXiv

Alfven Wave Turbulence and Perpendicular Ion Temperatures in Coronal Holes

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.

preprint2010arXiv

Perpendicular Ion Heating by Low-Frequency Alfven-Wave Turbulence in the Solar Wind

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.

preprint2009arXiv

Alfven Wave Reflection and Turbulent Heating in the Solar Wind from 1 Solar Radius to 1 AU: an Analytical Treatment

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.

preprint2009arXiv

Constraining Low-Frequency Alfvenic Turbulence in the Solar Wind Using Density Fluctuation Measurements

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.

preprint2009arXiv

Resonant Interactions Between Protons and Oblique Alfvén/Ion-Cyclotron Waves

Resonant interactions between ions and Alfvén/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 $θ$ 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.