Showing 1–50 of 90 results for author: Chandran, B D G

Scale-Dependent Dynamic Alignment in MHD Turbulence: Insights into Intermittency, Compressibility, and Imbalance Effects

Authors: Nikos Sioulas, Marco Velli, Alfred Mallet, Trevor A. Bowen, B. D. G. Chandran, Chen Shi, S. S. Cerri, Ioannis Liodis, Tamar Ervin, Davin E. Larson

Abstract: Scale-Dependent Dynamic Alignment (SDDA) in Elsässer field fluctuations is theorized to suppress nonlinearities and modulate the energy spectrum. Limited empirical evidence exists for SDDA within the solar wind turbulence's inertial range. We analyzed data from the WIND mission to assess the effects of compressibility, intermittency, and imbalance on SDDA. SDDA consistently appears at energy-conta… ▽ More Scale-Dependent Dynamic Alignment (SDDA) in Elsässer field fluctuations is theorized to suppress nonlinearities and modulate the energy spectrum. Limited empirical evidence exists for SDDA within the solar wind turbulence's inertial range. We analyzed data from the WIND mission to assess the effects of compressibility, intermittency, and imbalance on SDDA. SDDA consistently appears at energy-containing scales, with a trend toward misalignment at inertial scales. Compressible fluctuations show no increased alignment; however, their impact on SDDA's overall behavior is minimal. The alignment angles inversely correlate with field gradient intensity, likely due to "anomalous" or "counterpropagating" wave packet interactions. This suggests that SDDA originates from mutual shearing of Elsässer fields during imbalanced ($δ\boldsymbol{z}^{\pm} \gg δ\boldsymbol{z}^{\mp}$) interactions. Rigorous thresholding on field gradient intensity reveals SDDA signatures across much of the inertial range. The scaling of Elsässer increments' alignment angle, $Θ^{z}$, steepens with increasing global Alfvénic imbalance, while the angle between magnetic and velocity field increments, $Θ^{ub}$, becomes shallower. $Θ^{ub}$ only correlates with global Elsässer imbalance, steepening as the imbalance increases. Furthermore, increasing alignment in $Θ^{ub}$ persists deep into the inertial range of balanced intervals but collapses at large scales for imbalanced ones. Simplified theoretical analysis and modeling of high-frequency, low-amplitude noise in the velocity field indicate significant impacts on alignment angle measurements even at very low frequencies, with effects growing as global imbalance increases. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Extended Cyclotron Resonant Heating of the Turbulent Solar Wind

Authors: Trevor A. Bowen, Ivan Y. Vasko, Stuart D. Bale, Benjamin D. G. Chandran, Alexandros Chasapis, Thierry Dudok de Wit, Alfred Mallet, Michael McManus, Romain Meyrand, Marc Pulupa, Jonathan Squire

Abstract: Circularly polarized, nearly parallel propagating waves are prevalent in the solar wind at ion-kinetic scales. At these scales, the spectrum of turbulent fluctuations in the solar wind steepens, often called the transition-range, before flattening at sub-ion scales. Circularly polarized waves have been proposed as a mechanism to couple electromagnetic fluctuations to ion gyromotion, enabling ion-s… ▽ More Circularly polarized, nearly parallel propagating waves are prevalent in the solar wind at ion-kinetic scales. At these scales, the spectrum of turbulent fluctuations in the solar wind steepens, often called the transition-range, before flattening at sub-ion scales. Circularly polarized waves have been proposed as a mechanism to couple electromagnetic fluctuations to ion gyromotion, enabling ion-scale dissipation that results in observed ion-scale steepening. Here, we study Parker Solar Probe observations of an extended stream of fast solar wind ranging from 15-55 solar radii. We demonstrate that, throughout the stream, transition-range steepening at ion-scales is associated with the presence of significant left handed ion-kinetic scale waves, which are thought to be ion-cyclotron waves. We implement quasilinear theory to compute the rate at which ions are heated via cyclotron resonance with the observed circularly polarized waves given the empirically measured proton velocity distribution functions. We apply the Von Karman decay law to estimate the turbulent decay of the large-scale fluctuations, which is equal to the turbulent energy cascade rate. We find that the ion-cyclotron heating rates are correlated with, and amount to a significant fraction of, the turbulent energy cascade rate, implying that cyclotron heating is an important dissipation mechanism in the solar wind. △ Less

Submitted 14 June, 2024; originally announced June 2024.

3-Minute Oscillations in the Upper Corona: Evidence from Parker Solar Probe

Authors: Zesen Huang, Marco Velli, Chen Shi, Yingjie Zhu, B. D. G. Chandran, Victor Réville, Trevor Bowen, Nikos Sioulas, Marc Pulupa, Jia Huang, Sheng Huang

Abstract: Recent observations of Parker Solar Probe (PSP) from around the Alfvén surface have shown that the trace magnetic power spectrum density (PSD) is often characterized by a shallow-inertial double power law, where in the low frequency energy injection range, the power spectrum is shallow (flatter than $1/f$), and in the inertial range the spectrum is steep, with a scaling index of [1.5, 1.67]. Conse… ▽ More Recent observations of Parker Solar Probe (PSP) from around the Alfvén surface have shown that the trace magnetic power spectrum density (PSD) is often characterized by a shallow-inertial double power law, where in the low frequency energy injection range, the power spectrum is shallow (flatter than $1/f$), and in the inertial range the spectrum is steep, with a scaling index of [1.5, 1.67]. Consequently, close to the sun, the majority of the fluctuation energy concentrates in a small frequency range around the low frequency power spectral break. In this work, we conduct a systematic survey of PSP observations for the first 17 encounters to statistically study the energy behaviors of the magnetic fluctuations. Our results show that the center frequency of fluctuation energy systematically drifts to around 3-minute for the most pristine solar wind (smallest solar wind advection time). Moreover, the center frequency rapidly drifts to lower frequency as solar wind advection time increases, as expected for active turbulence. The concentration of fluctuation energy around 3-minutes suggests that Alfvénic fluctuations in solar wind might mostly be coming from resonant p-mode oscillations in the photosphere, though other potential sources are discussed. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Higher-Order Analysis of Three-Dimensional Anisotropy in Imbalanced Alfvénic Turbulence

Authors: Nikos Sioulas, Themistocles Zikopoulos, Chen Shi, Marco Velli, Trevor Bowen, Alfred Mallet, Luca Sorriso-Valvo, Andrea Verdini, B. D. G. Chandran, Mihailo M. Martinović, S. S. Cerri, Nooshin Davis, Corina Dunn

Abstract: We analyze in-situ observations of imbalanced solar wind turbulence to evaluate MHD turbulence models grounded in "Critical Balance" (CB) and "Scale-Dependent Dynamic Alignment" (SDDA). At energy injection scales, both outgoing and ingoing modes exhibit a weak cascade; a simultaneous tightening of SDDA is noted. Outgoing modes persist in a weak cascade across the inertial range, while ingoing mode… ▽ More We analyze in-situ observations of imbalanced solar wind turbulence to evaluate MHD turbulence models grounded in "Critical Balance" (CB) and "Scale-Dependent Dynamic Alignment" (SDDA). At energy injection scales, both outgoing and ingoing modes exhibit a weak cascade; a simultaneous tightening of SDDA is noted. Outgoing modes persist in a weak cascade across the inertial range, while ingoing modes shift to a strong cascade at $λ\approx 3 \times 10^{4} d_i$, with associated spectral scalings deviating from expected behavior due to "anomalous coherence" effects. The inertial range comprises two distinct sub-inertial segments. Beyond $λ\gtrsim 100 d_i$, eddies adopt a field-aligned tube topology, with SDDA signatures mainly evident in high amplitude fluctuations. The scaling exponents $ζ_{n}$ of the $n$-th order conditional structure functions, orthogonal to both the local mean field and fluctuation direction, align with the analytical models of Chandran et al. 2015 and Mallet et al. 2017, indicating "multifractal" statistics and strong intermittency; however, scaling in parallel and displacement components is more concave than predicted, possibly influenced by expansion effects. Below $λ\approx 100 d_i$, eddies become increasingly anisotropic, evolving into thin current sheet-like structures. Concurrently, $ζ_{n}$ scales linearly with order, marking a shift towards "monofractal" statistics. At $λ\approx 8 d_i$, the increase in aspect ratio halts, and the eddies become quasi-isotropic. This change may signal tearing instability, leading to reconnection, or result from energy redirection into the ion-cyclotron wave spectrum, aligning with the "helicity barrier". Our analysis utilizes 5-point structure functions, proving more effective than the traditional 2-point method in capturing steep scaling behaviors at smaller scales. △ Less

Submitted 5 April, 2024; originally announced April 2024.

On the Heating of the Slow Solar-Wind by Imbalanced Alfvén-Wave Turbulence from 0.06 au to 1 au: Parker Solar Probe and Solar Orbiter observations

Authors: Sofiane Bourouaine, Jean C. Perez, Benjamin D. G. Chandran, Vamsee K. Jagarlamudi, Nour E. Raouafi, Jasper S. Halekas

Abstract: In this work we analyze plasma and magnetic field data provided by the Parker Solar Probe (\emph{PSP}) and Solar Orbiter (\emph{SO}) missions to investigate the radial evolution of the heating of Alfvénic slow wind (ASW) by imbalanced Alfvén-Wave (AW) turbulent fluctuations from 0.06 au to 1 au. in our analysis we focus on slow solar-wind intervals with highly imbalanced and incompressible turbule… ▽ More In this work we analyze plasma and magnetic field data provided by the Parker Solar Probe (\emph{PSP}) and Solar Orbiter (\emph{SO}) missions to investigate the radial evolution of the heating of Alfvénic slow wind (ASW) by imbalanced Alfvén-Wave (AW) turbulent fluctuations from 0.06 au to 1 au. in our analysis we focus on slow solar-wind intervals with highly imbalanced and incompressible turbulence (i.e., magnetic compressibility $C_B=δB/B\leq 0.25$, plasma compressibility $C_n=δn/n\leq 0.25$ and normalized cross-helicity $σ_c\geq 0.65$). First, we estimate the AW turbulent dissipation rate from the wave energy equation and find that the radial profile trend is similar to the proton heating rate. Second, we find that the scaling of the empirical AW turbulent dissipation rate $Q_W$ obtained from the wave energy equation matches the scaling from the phenomenological AW turbulent dissipation rate $Q_{\rm CH09}$ (with $Q_{\rm CH09}\simeq 1.55 Q_W$) derived by~\cite{chandran09} based on the model of reflection-driven turbulence. Our results suggest that, as in the fast solar wind, AW turbulence plays a major role in the ion heating that occurs in incompressible slow-wind streams. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: This paper has been accepted for publication in the Astrophysical Journal Letters

The Effects of Non-Equilibrium Velocity Distributions on Alfvén Ion-Cyclotron Waves in the Solar Wind

Authors: Jada Walters, Kristopher G. Klein, Emily Lichko, Michael L. Stevens, Daniel Verscharen, Benjamin D. G. Chandran

Abstract: In this work, we investigate how the complex structure found in solar wind proton velocity distribution functions (VDFs), rather than the commonly assumed two-component bi-Maxwellian structure, affects the onset and evolution of parallel-propagating microinstabilities. We use the Arbitrary Linear Plasma Solver (ALPS), a numerical dispersion solver, to find the real frequencies and growth/damping r… ▽ More In this work, we investigate how the complex structure found in solar wind proton velocity distribution functions (VDFs), rather than the commonly assumed two-component bi-Maxwellian structure, affects the onset and evolution of parallel-propagating microinstabilities. We use the Arbitrary Linear Plasma Solver (ALPS), a numerical dispersion solver, to find the real frequencies and growth/damping rates of the Alfvén modes calculated for proton VDFs extracted from Wind spacecraft observations of the solar wind. We compare this wave behavior to that obtained by applying the same procedure to core-and-beam bi-Maxwellian fits of the Wind proton VDFs. We find several significant differences in the plasma waves obtained for the extracted data and bi-Maxwellian fits, including a strong dependence of the growth/damping rate on the shape of the VDF. By application of the quasilinear diffusion operator to these VDFs, we pinpoint resonantly interacting regions in velocity space where differences in VDF structure significantly affect the wave growth and damping rates. This demonstration of the sensitive dependence of Alfvén mode behavior on VDF structure may explain why the Alfvén ion-cyclotron instability thresholds predicted by linear theory for bi-Maxwellian models of solar wind proton background VDFs do not entirely constrain spacecraft observations of solar wind proton VDFs, such as those made by the Wind spacecraft. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 11 pages, 4 figures, 1 table. To be published in The Astrophysical Journal

Reflection-driven turbulence in the super-Alfvénic solar wind

Authors: Romain Meyrand, Jonathan Squire, Alfred Mallet, Benjamin D. G. Chandran

Abstract: In magnetized, stratified astrophysical environments such as the Sun's corona and solar wind, Alfvénic fluctuations ''reflect'' from background gradients, enabling nonlinear interactions and thus dissipation of their energy into heat. This process, termed ''reflection-driven turbulence,'' is thought to play a crucial role in coronal heating and solar-wind acceleration, explaining a range of observ… ▽ More In magnetized, stratified astrophysical environments such as the Sun's corona and solar wind, Alfvénic fluctuations ''reflect'' from background gradients, enabling nonlinear interactions and thus dissipation of their energy into heat. This process, termed ''reflection-driven turbulence,'' is thought to play a crucial role in coronal heating and solar-wind acceleration, explaining a range of observational correlations and constraints. Building on previous works focused on the inner heliosphere, here we study the basic physics of reflection-driven turbulence using reduced magnetohydrodynamics in an expanding box -- the simplest model that can capture the local turbulent plasma dynamics in the super-Alfvénic solar wind. Although idealized, our high-resolution simulations and simple theory reveal a rich phenomenology that is consistent with a diverse range of observations. Outwards-propagating fluctuations, which initially have high imbalance, decay nonlinearly to heat the plasma, becoming more balanced and magnetically dominated. Despite the high imbalance, the turbulence is strong because Elsässer collisions are suppressed by reflection, leading to ''anomalous coherence'' between the two Elsässer fields. This coherence, together with linear effects, causes the turbulence to anomalously grow the ''anastrophy'' (squared magnetic potential) as it decays, forcing the energy to rush to larger scales and forming a ''$1/f$-range'' energy spectrum as it does so. At late times, the expansion overcomes the nonlinear and Alfvénic physics, forming isolated, magnetically dominated ''Alfvén vortex'' structures that minimize their nonlinear dissipation. These results can plausibly explain the observed radial and wind-speed dependence of turbulence imbalance, residual energy, plasma heating, and fluctuation spectra, as well as making testable predictions for future observations. △ Less

Submitted 20 August, 2023; originally announced August 2023.

Mediation of Collisionless Turbulent Dissipation Through Cyclotron Resonance

Authors: Trevor A. Bowen, Stuart D. Bale, Benjamin D. G. Chandran, Alexandros Chasapis, Christopher H. K. Chen, Thierry Dudok de Wit, Alfred Mallet, Romain Meyrand, Jonathan Squire

Abstract: The dissipation of magnetized turbulence is fundamental to understanding energy transfer and heating in astrophysical systems. Collisionless interactions, such as resonant wave-particle process, are known to play a role in shaping turbulent astrophysical environments. Here, we present evidence for the mediation of turbulent dissipation in the solar wind by ion-cyclotron waves. Our results show tha… ▽ More The dissipation of magnetized turbulence is fundamental to understanding energy transfer and heating in astrophysical systems. Collisionless interactions, such as resonant wave-particle process, are known to play a role in shaping turbulent astrophysical environments. Here, we present evidence for the mediation of turbulent dissipation in the solar wind by ion-cyclotron waves. Our results show that ion-cyclotron waves interact strongly with magnetized turbulence, indicating that they serve as a major pathway for the dissipation of large-scale electromagnetic fluctuations. We further show that the presence of cyclotron waves significantly weakens observed signatures of intermittency in sub-ion-kinetic turbulence, which are known to be another pathway for dissipation. These observations results suggest that in the absence of cyclotron resonant waves, non-Gaussian, coherent structures are able to form at sub-ion-kinetic scales, and are likely responsible for turbulent heating. We further find that the cross helicity, i.e. the level of Alfvénicity of the fluctuations, correlates strongly with the presence of ion-scale waves, demonstrating that dissipation of collisionless plasma turbulence is not a universal process, but that the pathways to heating and dissipation at small scales are controlled by the properties of the large-scale turbulent fluctuations. We argue that these observations support the existence of a helicity barrier, in which highly Alfvénic, imbalanced, turbulence is prevented from cascading to sub-ion scales thus resulting in significant ion-cyclotron resonant heating. Our results may serve as a significant step in constraining the nature of turbulent heating in a wide variety of astrophysical systems. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Quantifying the Energy Budget in the Solar Wind from 13.3-100 Solar Radii

Authors: J. S. Halekas, S. D. Bale, M. Berthomier, B. D. G. Chandran, J. F. Drake, J. C. Kasper, K. G. Klein, D. E. Larson, R. Livi, M. P. Pulupa, M. L. Stevens, J. L. Verniero, P. Whittlesey

Abstract: A variety of energy sources, ranging from dynamic processes like magnetic reconnection and waves to quasi-steady terms like the plasma pressure, may contribute to the acceleration of the solar wind. We utilize a combination of charged particle and magnetic field observations from the Parker Solar Probe (PSP) to attempt to quantify the steady-state contribution of the proton pressure, the electric… ▽ More A variety of energy sources, ranging from dynamic processes like magnetic reconnection and waves to quasi-steady terms like the plasma pressure, may contribute to the acceleration of the solar wind. We utilize a combination of charged particle and magnetic field observations from the Parker Solar Probe (PSP) to attempt to quantify the steady-state contribution of the proton pressure, the electric potential, and the wave energy to the solar wind proton acceleration observed by PSP between 13.3 and ~100 solar radii (RS). The proton pressure provides a natural kinematic driver of the outflow. The ambipolar electric potential acts to couple the electron pressure to the protons, providing another definite proton acceleration term. Fluctuations and waves, while inherently dynamic, can act as an additional effective steady-state pressure term. To analyze the contributions of these terms, we utilize radial binning of single-point PSP measurements, as well as repeated crossings of the same stream at different distances on individual PSP orbits (i.e. "fast radial scans"). In agreement with previous work, we find that the electric potential contains sufficient energy to fully explain the acceleration of the slower wind streams. On the other hand, we find that the wave pressure plays an increasingly important role in the faster wind streams. The combination of these terms can explain the continuing acceleration of both slow and fast wind streams beyond 13.3 RS. △ Less

Submitted 22 May, 2023; originally announced May 2023.

The Evolution of the 1/f Range Within a Single Fast-Solar-Wind Stream Between 17.4 and 45.7 Solar Radii

Authors: Nooshin Davis, B. D. G. Chandran, T. A. Bowen, S. T. Badman, T. Dudok de Wit, C. H. K. Chen, S. D. Bale, Zesen Huang, Nikos Sioulas, Marco Velli

Abstract: The power spectrum of magnetic-field fluctuations in the fast solar wind ($V_{\rm SW}> 500 \mbox{ km} \mbox{ s}^{-1}$) at magnetohydrodynamic (MHD) scales is characterized by two different power laws on either side of a break frequency $f_{\rm b}$. The low-frequency range at frequencies $f$ smaller than $f_{\rm b}$ is often viewed as the energy reservoir that feeds the turbulent cascade at… ▽ More The power spectrum of magnetic-field fluctuations in the fast solar wind ($V_{\rm SW}> 500 \mbox{ km} \mbox{ s}^{-1}$) at magnetohydrodynamic (MHD) scales is characterized by two different power laws on either side of a break frequency $f_{\rm b}$. The low-frequency range at frequencies $f$ smaller than $f_{\rm b}$ is often viewed as the energy reservoir that feeds the turbulent cascade at $f>f_{\rm b}$. At heliocentric distances $r$ exceeding $60$ solar radii ($R_{\rm s}$), the power spectrum often has a $1/f$ scaling at $f<f_{\rm b}$; i.e., the spectral index is close to $-1$. In this study, measurements from the encounter $10$ of ${Parker Solar Probe}$ (PSP) with the Sun are used to investigate the evolution of the magnetic-field power spectrum at $f< f_{\rm b}$ at $r<60 R_{\rm s}$ during a fast radial scan of a single fast-solar-wind stream. We find that the spectral index in the low-frequency part of the spectrum decreases from approximately $-0.61$ to $-0.94$ as $r$ increases from $17.4 $ to $45.7$ solar radii. Our results suggest that the $1/f $ spectrum that is often seen at large $r$ in the fast solar wind is not produced at the Sun, but instead develops dynamically as the wind expands outward from the corona into the interplanetary medium. △ Less

Submitted 2 March, 2023; originally announced March 2023.

New Observations of Solar Wind 1/f Turbulence Spectrum from Parker Solar Probe

Authors: Zesen Huang, Nikos Sioulas, Chen Shi, Marco Velli, Trevor Bowen, Nooshin Davis, B. D. G. Chandran, Ning Kang, Xiaofei Shi, Jia Huang, Stuart D. Bale, J. C. Kasper, Davin E. Larson, Roberto Livi, P. L. Whittlesey, Ali Rahmati, Kristoff Paulson, M. Stevens, A. W. Case, Thierry Dudok de Wit, David M. Malaspina, J. W. Bonnell, Keith Goetz, Peter R. Harvey, Robert J. MacDowall

Abstract: The trace magnetic power spectrum in the solar wind is known to be characterized by a double power law at scales much larger than the proton gyro-radius, with flatter spectral exponents close to -1 found at the lower frequencies below an inertial range with indices closer to $[-1.5,-1.6]$. The origin of the $1/f$ range is still under debate. In this study, we selected 109 magnetically incompressib… ▽ More The trace magnetic power spectrum in the solar wind is known to be characterized by a double power law at scales much larger than the proton gyro-radius, with flatter spectral exponents close to -1 found at the lower frequencies below an inertial range with indices closer to $[-1.5,-1.6]$. The origin of the $1/f$ range is still under debate. In this study, we selected 109 magnetically incompressible solar wind intervals ($δ|\boldsymbol B|/|\boldsymbol B| \ll 1$) from Parker Solar Probe encounters 1 to 13 which display such double power laws, with the aim of understanding the statistics and radial evolution of the low frequency power spectral exponents from Alfvén point up to 0.3 AU. New observations from closer to the sun show that in the low frequency range solar wind turbulence can display spectra much shallower than $1/f$, evolving asymptotically to $1/f$ as advection time increases, indicating a dynamic origin for the $1/f$ range formation. We discuss the implications of this result on the Matteini et al. (2018) conjecture for the $1/f$ origin as well as example spectra displaying a triple power law consistent with the model proposed by Chandran et al. (2018), supporting the dynamic role of parametric decay in the young solar wind. Our results provide new constraints on the origin of the $1/f$ spectrum and further show the possibility of the coexistence of multiple formation mechanisms. △ Less

Submitted 23 May, 2023; v1 submitted 1 March, 2023; originally announced March 2023.

Comments: Accepted by ApJL

On the evolution of the Anisotropic Scaling of Magnetohydrodynamic Turbulence in the Inner Heliosphere

Authors: Nikos Sioulas, Marco Velli, Zesen Huang, Chen Shi, Trevor A. Bowen, B. D. G. Chandran, Ioannis Liodis, Nooshin Davis, Stuart D. Bale, T. S. Horbury, Thierry Dudok de Wit, Davin Larson, Justin Kasper, Christopher J. Owen, Michael L. Stevens, Anthony Case, Marc Pulupa, David M. Malaspina, J. W. Bonnell, Keith Goetz, Peter R. Harvey, Robert J. MacDowall

Abstract: We analyze a merged Parker Solar Probe ($PSP$) and Solar Orbiter ($SO$) dataset covering heliocentric distances $13 \ R_{\odot} \lesssim R \lesssim 220$ $R_{\odot}$ to investigate the radial evolution of power and spectral-index anisotropy in the wavevector space of solar wind turbulence. Our results show that anisotropic signatures of turbulence display a distinct radial evolution when fast,… ▽ More We analyze a merged Parker Solar Probe ($PSP$) and Solar Orbiter ($SO$) dataset covering heliocentric distances $13 \ R_{\odot} \lesssim R \lesssim 220$ $R_{\odot}$ to investigate the radial evolution of power and spectral-index anisotropy in the wavevector space of solar wind turbulence. Our results show that anisotropic signatures of turbulence display a distinct radial evolution when fast, $V_{sw} \geq ~ 400 ~km ~s^{-1}$, and slow, $V_{sw} \leq ~ 400 ~km ~s^{-1}$, wind streams are considered. The anisotropic properties of slow wind in Earth orbit are consistent with a ``critically balanced'' cascade, but both spectral-index anisotropy and power anisotropy diminish with decreasing heliographic distance. Fast streams are observed to roughly retain their near-Sun anisotropic properties, with the observed spectral index and power anisotropies being more consistent with a ``dynamically aligned'' type of cascade, though the lack of extended fast-wind intervals makes it difficult to accurately measure the anisotropic scaling. A high-resolution analysis during the first perihelion of PSP confirms the presence of two sub-ranges within the inertial range, which may be associated with the transition from weak to strong turbulence. The transition occurs at $κd_{i} \approx 6 \times 10^{-2}$, and signifies a shift from -5/3 to -2 and -3/2 to -1.57 scaling in parallel and perpendicular spectra, respectively. Our results provide strong observational constraints for anisotropic theories of MHD turbulence in the solar wind. △ Less

Submitted 20 March, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: Accepted to APJ

Parker Solar Probe: Four Years of Discoveries at Solar Cycle Minimum

Authors: N. E. Raouafi, L. Matteini, J. Squire, S. T. Badman, M. Velli, K. G. Klein, C. H. K. Chen, W. H. Matthaeus, A. Szabo, M. Linton, R. C. Allen, J. R. Szalay, R. Bruno, R. B. Decker, M. Akhavan-Tafti, O. V. Agapitov, S. D. Bale, R. Bandyopadhyay, K. Battams, L. Berčič, S. Bourouaine, T. Bowen, C. Cattell, B. D. G. Chandran, R. Chhiber , et al. (32 additional authors not shown)

Abstract: Launched on 12 Aug. 2018, NASA's Parker Solar Probe had completed 13 of its scheduled 24 orbits around the Sun by Nov. 2022. The mission's primary science goal is to determine the structure and dynamics of the Sun's coronal magnetic field, understand how the solar corona and wind are heated and accelerated, and determine what processes accelerate energetic particles. Parker Solar Probe returned a… ▽ More Launched on 12 Aug. 2018, NASA's Parker Solar Probe had completed 13 of its scheduled 24 orbits around the Sun by Nov. 2022. The mission's primary science goal is to determine the structure and dynamics of the Sun's coronal magnetic field, understand how the solar corona and wind are heated and accelerated, and determine what processes accelerate energetic particles. Parker Solar Probe returned a treasure trove of science data that far exceeded quality, significance, and quantity expectations, leading to a significant number of discoveries reported in nearly 700 peer-reviewed publications. The first four years of the 7-year primary mission duration have been mostly during solar minimum conditions with few major solar events. Starting with orbit 8 (i.e., 28 Apr. 2021), Parker flew through the magnetically dominated corona, i.e., sub-Alfvénic solar wind, which is one of the mission's primary objectives. In this paper, we present an overview of the scientific advances made mainly during the first four years of the Parker Solar Probe mission, which go well beyond the three science objectives that are: (1) Trace the flow of energy that heats and accelerates the solar corona and solar wind; (2) Determine the structure and dynamics of the plasma and magnetic fields at the sources of the solar wind; and (3) Explore mechanisms that accelerate and transport energetic particles. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 157 pages, 65 figures

Journal ref: Space Science Reviews 2023

Electron-driven instabilities in the solar wind

Authors: Daniel Verscharen, Benjamin D. G. Chandran, Elisabetta Boella, Jasper Halekas, Maria Elena Innocenti, Vamsee K. Jagarlamudi, Alfredo Micera, Viviane Pierrard, Stepan Stverak, Ivan Y. Vasko, Marco Velli, Phyllis L. Whittlesey

Abstract: 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-insta… ▽ More 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. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Frontiers in Astronomy and Space Sciences. 47 pages, 11 figures

The In Situ Signature of Cyclotron Resonant Heating

Authors: Trevor A. Bowen, Benjmin D. G. Chandran, Jonathan Squire, Stuart D. Bale, Die Duan, Kristopher G. Klein, Davin Larson, Alfred Mallet, Michael D. McManus, Romain Meyrand, Jaye L. Verniero, Lloyd D. Woodham

Abstract: The dissipation of magnetized turbulence is an important paradigm for describing heating and energy transfer in astrophysical environments such as the solar corona and wind; however, the specific collisionless processes behind dissipation and heating remain relatively unconstrained by measurements. Remote sensing observations have suggested the presence of strong temperature anisotropy in the sola… ▽ More The dissipation of magnetized turbulence is an important paradigm for describing heating and energy transfer in astrophysical environments such as the solar corona and wind; however, the specific collisionless processes behind dissipation and heating remain relatively unconstrained by measurements. Remote sensing observations have suggested the presence of strong temperature anisotropy in the solar corona consistent with cyclotron resonant heating. In the solar wind, in situ magnetic field measurements reveal the presence of cyclotron waves, while measured ion velocity distribution functions have hinted at the active presence of cyclotron resonance. Here, we present Parker Solar Probe observations that connect the presence of ion-cyclotron waves directly to signatures of resonant damping in observed proton-velocity distributions. We show that the observed cyclotron wave population coincides with both flattening in the phase space distribution predicted by resonant quasilinear diffusion and steepening in the turbulent spectra at the ion-cyclotron resonant scale. In measured velocity distribution functions where cyclotron resonant flattening is weaker, the distributions are nearly uniformly subject to ion-cyclotron wave damping rather than emission, indicating that the distributions can damp the observed wave population. These results are consistent with active cyclotron heating in the solar wind. △ Less

Submitted 28 November, 2022; v1 submitted 9 November, 2021; originally announced November 2021.

Strong perpendicular velocity-space in proton beams observed by Parker Solar Probe

Authors: J. L. Verniero, B. D. G. Chandran, D. E. Larson, K. Paulson, B. L. Alterman, S. Badman, S. D. Bale, J. W. Bonnell, T. A. Bowen, T. Dudok de Wit, J. C. Kasper, K. G. Klein, E. Lichko, R. Livi, M. D. McManus, A. Rahmati, D. Verscharen, J. Walters, P. L. Whittlesey

Abstract: The SWEAP instrument suite on Parker Solar Probe (PSP) has detected numerous proton beams associated with coherent, circularly polarized, ion-scale waves observed by PSP's FIELDS instrument suite. Measurements during PSP Encounters 4-8 revealed pronounced complex shapes in the proton velocity distribution functions (VDFs), in which the tip of the beam undergoes strong perpendicular diffusion, resu… ▽ More The SWEAP instrument suite on Parker Solar Probe (PSP) has detected numerous proton beams associated with coherent, circularly polarized, ion-scale waves observed by PSP's FIELDS instrument suite. Measurements during PSP Encounters 4-8 revealed pronounced complex shapes in the proton velocity distribution functions (VDFs), in which the tip of the beam undergoes strong perpendicular diffusion, resulting in VDF level contours that resemble a `hammerhead.' We refer to these proton beams, with their attendant `hammerhead' features, as the ion strahl. We present an example of these observations occurring simultaneously with a 7-hour ion-scale wave storm and show results from a preliminary attempt at quantifying the occurrence of ion-strahl broadening through 3-component ion-VDF fitting. We also provide a possible explanation of the ion perpendicular scattering based on quasilinear theory and the resonant scattering of beam ions by parallel-propagating, right circularly polarized, fast-magnetosonic/whistler waves. △ Less

Submitted 17 October, 2021; originally announced October 2021.

Comments: 13 pages, 8 figures, Under Review in the Astrophysical Journal

A solar source of Alfvénic magnetic field switchbacks: {\em in situ} remnants of magnetic funnels on supergranulation scales

Authors: S. D. Bale, T. S. Horbury, M. Velli, M. I. Desai, J. S. Halekas, M. D. McManus, O. Panasenco, S. T. Badman, T. A. Bowen, B. D. G. Chandran, J. F. Drake, J. C. Kasper, R. Laker, A. Mallet, L Matteini, T. D. Phan, N. E. Raouafi, J. Squire, L. D. Woodham, T. Wooley

Abstract: One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfvénic magnetic field reversals termed 'switchbacks'. These $δB_R/B \sim \mathcal{O}(1$) fluctuations occur on a range of timescales and in {\em patches} separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate tha… ▽ More One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfvénic magnetic field reversals termed 'switchbacks'. These $δB_R/B \sim \mathcal{O}(1$) fluctuations occur on a range of timescales and in {\em patches} separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate that patches of switchbacks are localized within the extensions of plasma structures originating at the base of the corona. These structures are characterized by an increase in alpha particle abundance, Mach number, plasma $β$ and pressure, and by depletions in the magnetic field magnitude and electron temperature. These intervals are in pressure-balance, implying stationary spatial structure, and the field depressions are consistent with overexpanded flux tubes. The structures are asymmetric in Carrington longitude with a steeper leading edge and a small ($\sim$1$^\circ$) edge of hotter plasma and enhanced magnetic field fluctuations. Some structures contain suprathermal ions to $\sim$85 keV that we argue are the energetic tail of the solar wind alpha population. The structures are separated in longitude by angular scales associated with supergranulation. This suggests that these switchbacks originate near the leading edge of the diverging magnetic field funnels associated with the network magnetic field - the primary wind sources. We propose an origin of the magnetic field switchbacks, hot plasma and suprathermals, alpha particles in interchange reconnection events just above the solar transition region and our measurements represent the extended regions of a turbulent outflow exhaust. △ Less

Submitted 2 September, 2021; originally announced September 2021.

Comments: 18 pages, 9 figures, submitted to Astrophys. J

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

Authors: Alfred Mallet, Benjamin D. G. Chandran

Abstract: 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… ▽ More 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. △ Less

Submitted 1 June, 2021; v1 submitted 28 May, 2021; originally announced May 2021.

Comments: 5 pages

How Alfvén waves energize the solar wind: heat vs work

Authors: Jean C. Perez, Benjamin D. G. Chandran, Kristopher G. Klein, Mihailo M. Martinović

Abstract: A growing body of evidence suggests that the solar wind is powered to a large extent by an Alfvén-wave (AW) energy flux. AWs energize the solar wind via two mechanisms: heating and work. We use high-resolution direct numerical simulations of reflection-driven AW turbulence (RDAWT) in a fast-solar-wind stream emanating from a coronal hole to investigate both mechanisms. In particular, we compute th… ▽ More A growing body of evidence suggests that the solar wind is powered to a large extent by an Alfvén-wave (AW) energy flux. AWs energize the solar wind via two mechanisms: heating and work. We use high-resolution direct numerical simulations of reflection-driven AW turbulence (RDAWT) in a fast-solar-wind stream emanating from a coronal hole to investigate both mechanisms. In particular, we compute the fraction of the AW power at the coronal base ($P_{\rm AWb}$) that is transferred to solar-wind particles via heating between the coronal base and heliocentric distance $r$, which we denote $χ_{\rm H}(r)$, and the fraction that is transferred via work, which we denote $χ_{\rm W}(r)$. We find that $χ_{\rm W}(r_{\rm A})$ ranges from 0.15 to 0.3, where $r_{\rm A}$ is the Alfvén critical point. This value is small compared to~one because the Alfvén speed $v_{\rm A} $ exceeds the outflow velocity $U$ at $r<r_{\rm A}$, so the AWs race through the plasma without doing much work. At $r>r_{\rm A}$, where $v_{\rm A} < U$, the AWs are in an approximate sense "stuck to the plasma", which helps them do pressure work as the plasma expands. However, much of the AW power has dissipated by the time the AWs reach $r=r_{\rm A}$, so the total rate at which AWs do work on the plasma at $r>r_{\rm A}$ is a modest fraction of $P_{\rm AWb}$. We find that heating is more effective than work at $r<r_{\rm A}$, with $χ_{\rm H}(r_{\rm A})$ ranging from 0.5 to 0.7. The reason that $χ_{\rm H} \geq 0.5$ in our simulations is that an appreciable fraction of the local AW power dissipates within each Alfvén-speed scale height in RDAWT, and there are a few Alfvén-speed scale heights between the coronal base and $r_{\rm A}$. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 15 pages, 2 figures, accepted for publication in the Journal of Plasma Physics

Multiscale Solar Wind Turbulence Properties inside and near Switchbacks measured by Parker Solar Probe

Authors: Mihailo M. Martinović, Kristopher G. Klein, Jia Huang, Benjamin D. G. Chandran, Justin C. Kasper, Emily Lichko, Trevor Bowen, Christopher H. K. Chen, Lorenzo Matteini, Michael Stevens, Anthony W. Case, Stuart D. Bale

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

Submitted 27 February, 2021; originally announced March 2021.

Turbulent generation of magnetic switchbacks in the Alfvénic solar wind

Authors: Munehito Shoda, Benjamin D. G. Chandran, Steven R. Cranmer

Abstract: One of the most important early results from the Parker Solar Probe (PSP) is the ubiquitous presence of magnetic switchbacks, whose origin is under debate. Using a three-dimensional direct numerical simulation of the equations of compressible magnetohydrodynamics from the corona to 40 solar radii, we investigate whether magnetic switchbacks emerge from granulation-driven Alfvén waves and turbulenc… ▽ More One of the most important early results from the Parker Solar Probe (PSP) is the ubiquitous presence of magnetic switchbacks, whose origin is under debate. Using a three-dimensional direct numerical simulation of the equations of compressible magnetohydrodynamics from the corona to 40 solar radii, we investigate whether magnetic switchbacks emerge from granulation-driven Alfvén waves and turbulence in the solar wind. The simulated solar wind is an Alfvénic slow-solar-wind stream with a radial profile consistent with various observations, including observations from PSP. As a natural consequence of Alfvén-wave turbulence, the simulation reproduced magnetic switchbacks with many of the same properties as observed switchbacks, including Alfvénic v-b correlation, spherical polarization (low magnetic compressibility), and a volume filling fraction that increases with radial distance. The analysis of propagation speed and scale length shows that the magnetic switchbacks are large-amplitude (nonlinear) Alfvén waves with discontinuities in the magnetic field direction. We directly compare our simulation with observations using a virtual flyby of PSP in our simulation domain. We conclude that at least some of the switchbacks observed by PSP are a natural consequence of the growth in amplitude of spherically polarized Alfvén waves as they propagate away from the Sun. △ Less

Submitted 16 June, 2021; v1 submitted 23 January, 2021; originally announced January 2021.

Comments: accepted for publication in The Astrophysical Journal

An approximate analytic solution to the coupled problems of coronal heating and solar-wind acceleration

Authors: Benjamin D. G. Chandran

Abstract: Between the base of the solar corona and the Alfven critical point, the solar-wind density decreases by approximately five orders of magnitude, but the temperature varies by a factor of only a few. In this paper, I show that such quasi-isothermal evolution out to the Alfven critical point is a generic property of outflows powered by reflection-driven Alfven-wave (AW) turbulence, in which outward-p… ▽ More Between the base of the solar corona and the Alfven critical point, the solar-wind density decreases by approximately five orders of magnitude, but the temperature varies by a factor of only a few. In this paper, I show that such quasi-isothermal evolution out to the Alfven critical point is a generic property of outflows powered by reflection-driven Alfven-wave (AW) turbulence, in which outward-propagating AWs partially reflect, and counter-propagating AWs interact to produce a cascade of fluctuation energy to small scales, which leads to turbulent heating. Approximating the sub-Alfvenic region as isothermal, I first present a simplified calculation of the mass outflow rate, asymptotic wind speed, and coronal temperature that neglects conductive losses and the wave pressure force. I then develop a more detailed model of the transition region, corona, and solar wind that accounts for the heat flux from the coronal base into the transition region and momentum deposition by AWs. I solve analytically for this heat flux by balancing, within the transition region, conductive heating against internal-energy losses from radiation, pdV work, and advection. The density at the coronal base is determined by locally balancing turbulent heating and radiative cooling. I solve the equations of the model analytically in two different parameter regimes. Analytic and numerical solutions to the model equations agree with a number of observations. △ Less

Submitted 22 April, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: 42 pages, 7 figures, accepted for publication in the Journal of Plasma Physics. (Updated version includes the same minor editing as the proofs of the published paper.)

The Near-Sun Streamer Belt Solar Wind: Turbulence and Solar Wind Acceleration

Authors: C. H. K. Chen, B. D. G. Chandran, L. D. Woodham, S. I. Jones-Mecholsky, J. C. Perez, S. Bourouaine, T. A. Bowen, K. G. Klein, M. Moncuquet, J. C. Kasper, S. D. Bale

Abstract: The fourth orbit of Parker Solar Probe (PSP) reached heliocentric distances down to 27.9 Rs, allowing solar wind turbulence and acceleration mechanisms to be studied in situ closer to the Sun than previously possible. The turbulence properties were found to be significantly different in the inbound and outbound portions of PSP's fourth solar encounter, likely due to the proximity to the heliospher… ▽ More The fourth orbit of Parker Solar Probe (PSP) reached heliocentric distances down to 27.9 Rs, allowing solar wind turbulence and acceleration mechanisms to be studied in situ closer to the Sun than previously possible. The turbulence properties were found to be significantly different in the inbound and outbound portions of PSP's fourth solar encounter, likely due to the proximity to the heliospheric current sheet (HCS) in the outbound period. Near the HCS, in the streamer belt wind, the turbulence was found to have lower amplitudes, higher magnetic compressibility, a steeper magnetic field spectrum (with spectral index close to -5/3 rather than -3/2), a lower Alfvénicity, and a "1/f" break at much lower frequencies. These are also features of slow wind at 1 au, suggesting the near-Sun streamer belt wind to be the prototypical slow solar wind. The transition in properties occurs at a predicted angular distance of ~4° from the HCS, suggesting ~8° as the full-width of the streamer belt wind at these distances. While the majority of the Alfvénic turbulence energy fluxes measured by PSP are consistent with those required for reflection-driven turbulence models of solar wind acceleration, the fluxes in the streamer belt are significantly lower than the model predictions, suggesting that additional mechanisms are necessary to explain the acceleration of the streamer belt solar wind. △ Less

Submitted 1 January, 2021; originally announced January 2021.

Journal ref: A&A 650, L3 (2021)

Electron heat flux in the near-Sun environment

Authors: J. S. Halekas, P. L. Whittlesey, D. E. Larson, D. McGinnis, S. D. Bale, M. Berthomier, A. W. Case, B. D. G. Chandran, J. C. Kasper, K. G. Klein, K. E. Korreck, R. Livi, R. J. MacDowall, M. Maksimovic, D. M. Malaspina, L. Matteini, M. P. Pulupa, M. L. Stevens

Abstract: We survey the electron heat flux observed by the Parker Solar Probe (PSP) in the near-Sun environment at heliocentric distances of 0.125-0.25 AU. We utilized measurements from the Solar Wind Electrons Alphas and Protons and FIELDS experiments to compute the solar wind electron heat flux and its components and to place these in context. The PSP observations reveal a number of trends in the electron… ▽ More We survey the electron heat flux observed by the Parker Solar Probe (PSP) in the near-Sun environment at heliocentric distances of 0.125-0.25 AU. We utilized measurements from the Solar Wind Electrons Alphas and Protons and FIELDS experiments to compute the solar wind electron heat flux and its components and to place these in context. The PSP observations reveal a number of trends in the electron heat flux signatures near the Sun. The magnitude of the heat flux is anticorrelated with solar wind speed, likely as a result of the lower saturation heat flux in the higher-speed wind. When divided by the saturation heat flux, the resulting normalized net heat flux is anticorrelated with plasma beta on all PSP orbits, which is consistent with the operation of collisionless heat flux regulation mechanisms. The net heat flux also decreases in very high beta regions in the vicinity of the heliospheric current sheet, but in most cases of this type the omnidirectional suprathermal electron flux remains at a comparable level or even increases, seemingly inconsistent with disconnection from the Sun. The measured heat flux values appear inconsistent with regulation primarily by collisional mechanisms near the Sun. Instead, the observed heat flux dependence on plasma beta and the distribution of suprathermal electron parameters are both consistent with theoretical instability thresholds associated with oblique whistler and magnetosonic modes. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Journal ref: A&A 650, A15 (2021)

In-situ switchback formation in the expanding solar wind

Authors: Jonathan Squire, Benjamin D. G. Chandran, Romain Meyrand

Abstract: 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 wave… ▽ More 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. △ Less

Submitted 6 February, 2020; v1 submitted 23 January, 2020; originally announced January 2020.

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

Authors: Trevor A. Bowen, Alfred Mallet, Stuart D. Bale, J. W. Bonnell, Anthony W. Case, Benjamin D. G. Chandran, Alexandros Chasapis, Christopher H. K. Chen, Die Duan, Thierry Dudok de Wit, Keith Goetz, Jasper Halekas, Peter R. Harvey, J. C. Kasper, Kelly E. Korreck, Davin Larson, Roberto Livi, Robert J. MacDowall, David M. Malaspina, Marc Pulupa, Michael Stevens, Phyllis Whittlesey

Abstract: 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 demonstrat… ▽ More 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. △ Less

Submitted 14 January, 2020; originally announced January 2020.

Journal ref: Phys. Rev. Lett. 125, 025102 (2020)

Cross Helicity Reversals In Magnetic Switchbacks

Authors: Michael D. McManus, Trevor A. Bowen, Alfred Mallet, Christopher H. K. Chen, Benjamin D. G. Chandran, Stuart D. Bale, Davin E. Larson, Thierry Dudok de Wit, Justin C. Kasper, Michael Stevens, Phyllis Whittlesey, Roberto Livi, Kelly E. Korreck, Keith Goetz, Peter R. Harvey, Marc Pulupa, Robert J. MacDowall, David M. Malaspina, Anthony W. Case, John W. Bonnell

Abstract: 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 c… ▽ More 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. △ Less

Submitted 17 December, 2019; originally announced December 2019.

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

Authors: Mihailo M. Martinović, Kristopher G. Klein, Justin C. Kasper, Anthony W. Case, Kelly E. Korreck, Davin Larson, Roberto Livi, Michael Stevens, Phyllis Whittlesey, Benjamin D. G. Chandran, Ben L. Alterman, Jia Huang, Christopher H. K. Chen, Stuart D. Bale, Marc Pulupa, David M. Malaspina, John W. Bonnell, Peter R. Harvey, Keith Goetz, Thierry Dudok de Wit, Robert J. MacDowall

Abstract: 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 t… ▽ More 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. △ Less

Submitted 5 December, 2019; originally announced December 2019.

The Evolution and Role of Solar Wind Turbulence in the Inner Heliosphere

Authors: C. H. K. Chen, S. D. Bale, J. W. Bonnell, D. Borovikov, T. A. Bowen, D. Burgess, A. W. Case, B. D. G. Chandran, T. Dudok de Wit, K. Goetz, P. R. Harvey, J. C. Kasper, K. G. Klein, K. E. Korreck, D. Larson, R. Livi, R. J. MacDowall, D. M. Malaspina, A. Mallet, M. D. McManus, M. Moncuquet, M. Pulupa, M. Stevens, P. Whittlesey

Abstract: 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 leve… ▽ More 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énic 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én 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én 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. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Reflection-driven MHD turbulence in the solar atmosphere and solar wind

Authors: Benjamin D. G. Chandran, Jean C. Perez

Abstract: 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 impos… ▽ More 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én 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. △ Less

Submitted 26 August, 2019; v1 submitted 2 August, 2019; originally announced August 2019.

Comments: 40 pages, 7 figures, accepted for publication in the Journal of Plasma Physics (JPP). Includes proof corrections

Self-induced Scattering of Strahl Electrons in the Solar Wind

Authors: Daniel Verscharen, Benjamin D. G. Chandran, Seong-Yeop Jeong, Chadi S. Salem, Marc P. Pulupa, Stuart D. Bale

Abstract: 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 microinstab… ▽ More 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 $β_{\mathrm c}$ regimes, where $β_{\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-$β_{\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 $β_{\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én-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}. △ Less

Submitted 3 December, 2019; v1 submitted 6 June, 2019; originally announced June 2019.

Comments: 11 pages, 11 figures

Journal ref: Astrophys. J. 886, 2, 2019

[Plasma 2020 Decadal] Multipoint Measurements of the Solar Wind: A Proposed Advance for Studying Magnetized Turbulence

Authors: K. G. Klein, O. Alexandrova, J. Bookbinder, D. Caprioli, A. W. Case, B. D. G. Chandran, L. J. Chen, T. Horbury, L. Jian, J. C. Kasper, O. Le Contel, B. A. Maruca, W. Matthaeus, A. Retino, O. Roberts, A. Schekochihin, R. Skoug, C. Smith, J. Steinberg, H. Spence, B. Vasquez, J. M. TenBarge, D. Verscharen, P. Whittlesey

Abstract: 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, provide… ▽ More 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. △ Less

Submitted 13 August, 2019; v1 submitted 13 March, 2019; originally announced March 2019.

Comments: 6 pages, submitted to the Plasma 2020 Decadal Survey Committee; version 2 includes movie of proposed orbit

Hybrid-Kinetic Simulations of Ion Heating in Alfvénic Turbulence

Authors: Lev Arzamasskiy, Matthew W. Kunz, Benjamin D. G. Chandran, Eliot Quataert

Abstract: We present three-dimensional, hybrid-kinetic numerical simulations of driven Alfvé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 a… ▽ More We present three-dimensional, hybrid-kinetic numerical simulations of driven Alfvé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é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. △ Less

Submitted 17 May, 2019; v1 submitted 30 January, 2019; originally announced January 2019.

Comments: 14 pages, 10 figures, accepted to ApJ

Stochastic proton heating by kinetic-Alfvén-wave turbulence in moderately high-$β$ plasmas

Authors: Ian W. Hoppock, Benjamin D. G. Chandran, Kristopher G. Klein, Alfred Mallet, Daniel Verscharen

Abstract: 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… ▽ More 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én-wave (AW) and kinetic-Alfvén-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 $β_{\rm p} \sim 1-30$, where $β_{\rm p}$ is the ratio of the proton pressure to the magnetic pressure. (We do not consider the $β_{\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 $β_{\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 $β_{\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. △ Less

Submitted 21 November, 2018; originally announced November 2018.

Comments: 22 pages, 5 figures, accepted for publication in the Journal of Plasma Physics

NHDS: The New Hampshire Dispersion Relation Solver

Authors: Daniel Verscharen, Benjamin D. G. Chandran

Abstract: 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. 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. △ Less

Submitted 26 April, 2018; originally announced April 2018.

Comments: 3 pages, 1 figure

Journal ref: Research Notes of the AAS 2, 13, 2018

ALPS: The Arbitrary Linear Plasma Solver

Authors: Daniel Verscharen, Kristopher G. Klein, Benjamin D. G. Chandran, Michael L. Stevens, Chadi S. Salem, Stuart D. Bale

Abstract: 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 o… ▽ More 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, $κ$-distributed, and Jüttner-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. △ Less

Submitted 13 March, 2018; originally announced March 2018.

Comments: 26 pages, 13 figures, submitted

Parametric Instability, Inverse Cascade, and the $1/f$ Range of Solar-Wind Turbulence

Authors: Benjamin D. G. Chandran

Abstract: In this paper, weak turbulence theory is used to investigate the nonlinear evolution of the parametric instability in 3D low-$β$ 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én wave quanta, and several exact solutions… ▽ More In this paper, weak turbulence theory is used to investigate the nonlinear evolution of the parametric instability in 3D low-$β$ 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én 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én waves that initially satisfy $e^+ \gg e^-$, where $e^+$ and $e^-$ are the frequency ($\,f$) spectra of Alfvén 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. △ Less

Submitted 26 December, 2017; originally announced December 2017.

Comments: 25 pages, 4 figures, accepted for publication in the Journal of Plasma Physics

Heating of Accretion-Disk Coronae and Jets by General Relativistic MHD Turbulence

Authors: Benjamin D. G. Chandran, Francois Foucart, Alexander Tchekhovskoy

Abstract: Turbulence in an accretion disk launches Alfvén waves (AWs) that propagate away from the disk along magnetic field lines. Because the Alfvén 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 for… ▽ More Turbulence in an accretion disk launches Alfvén waves (AWs) that propagate away from the disk along magnetic field lines. Because the Alfvén 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 (ρ_{\rm b}/ρ_{\rm d})^{1/2} β_{\rm net,d}^{-1/2}$ times the disk's radiative flux, where $ρ_{\rm b}$ and $ρ_{\rm d}$ are the mass densities at the coronal base and disk midplane, respectively, and $β_{\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. △ Less

Submitted 17 April, 2018; v1 submitted 19 July, 2017; originally announced July 2017.

Comments: 21 pages, no figures, accepted for publication in JPP. Version 2 contains more detail and a new section on a sub-grid model for incorporating turbulence into GRMHD simulations

Disruption of Alfvénic turbulence by magnetic reconnection in a collisionless plasma

Authors: A. Mallet, A. A. Schekochihin, B. D. G. Chandran

Abstract: We calculate the disruption scale $λ_{\rm D}$ at which sheet-like structures in dynamically aligned Alfvénic turbulence are destroyed by the onset of magnetic reconnection in a low-$β$ collisionless plasma. The scaling of $λ_{\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 d… ▽ More We calculate the disruption scale $λ_{\rm D}$ at which sheet-like structures in dynamically aligned Alfvénic turbulence are destroyed by the onset of magnetic reconnection in a low-$β$ collisionless plasma. The scaling of $λ_{\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 $λ_{\rm D}\sim L_\perp^{1/9}(d_eρ_s)^{4/9}$, where $d_e$ is the electron inertial scale, $ρ_s$ is the ion sound scale, and $L_\perp$ is the outer scale of the turbulence. When $β_e$ and $ρ_s/L_\perp$ are sufficiently small, the scale $λ_{\rm D}$ is larger than $ρ_s$ and there is a break in the energy spectrum at $λ_{\rm D}$, rather than at $ρ_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 $ρ_s$. Both of these fractions depend strongly on $β_e$, with the disrupted structures becoming more important at lower $β_e$. Finally, we predict that the energy spectrum between $λ_{\rm D}$ and $ρ_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. △ Less

Submitted 12 October, 2017; v1 submitted 18 July, 2017; originally announced July 2017.

Comments: 21 pages, 1 figure, accepted for publication in JPP

Disruption of sheetlike structures in Alfvénic turbulence by magnetic reconnection

Authors: A. Mallet, A. A. Schekochihin, B. D. G. Chandran

Abstract: We propose a mechanism whereby the intense, sheet-like structures naturally formed by dynamically aligning Alfvénic turbulence are destroyed by magnetic reconnection at a scale $\hatλ_{\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 is… ▽ More We propose a mechanism whereby the intense, sheet-like structures naturally formed by dynamically aligning Alfvénic turbulence are destroyed by magnetic reconnection at a scale $\hatλ_{\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λ_{\rm D}$ the turbulent dynamics change: at scales larger than $\hatλ_{\rm D}$, the turbulence exhibits scale-dependent dynamic alignment and a spectral index approximately equal to $-3/2$, while at scales smaller than $\hatλ_{\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λ_{\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λ_{\rm D} \sim S_{L_\perp}^{-0.6}$. The spectral index below $\hatλ_{\rm D}$ is bounded between $-5/3$ and $-2.3$. The final dissipation scale is at $\hatλ_{η,\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. △ Less

Submitted 26 July, 2017; v1 submitted 22 December, 2016; originally announced December 2016.

Comments: 10 pages, 4 figures, published in MNRAS 468 p.4862, replaced with published version

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

Authors: Daniel Verscharen, Benjamin D. G. Chandran, Kristopher G. Klein, Eliot Quataert

Abstract: 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… ▽ More 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. △ Less

Submitted 2 November, 2016; v1 submitted 23 May, 2016; originally announced May 2016.

Comments: 11 pages, published in ApJ

Journal ref: Astrophys. J. 831, 128, 2016

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

Authors: Kristopher G. Klein, Benjamin D. G. Chandran

Abstract: 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 p… ▽ More 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. △ Less

Submitted 16 February, 2016; originally announced February 2016.

Comments: 6 pages, 6 figures, accepted for publication in The Astrophysical Journal

Measures of Three-Dimensional Anisotropy and Intermittency in Strong Alfvénic Turbulence

Authors: A. Mallet, A. A. Schekochihin, B. D. G. Chandran, C. H. K. Chen, T. S. Horbury, R. T. Wicks, C. C. Greenan

Abstract: We measure the local anisotropy of numerically simulated strong Alfvénic 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 elongat… ▽ More We measure the local anisotropy of numerically simulated strong Alfvénic 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énic 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. △ Less

Submitted 4 December, 2015; originally announced December 2015.

Comments: 9 pages, 5 figures. Submitted to MNRAS

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

Authors: Benjamin D. G. Chandran, Jean C. Perez, Daniel Verscharen, Kristopher G. Klein, Alfred Mallet

Abstract: 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 ba… ▽ More 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. △ Less

Submitted 3 September, 2015; originally announced September 2015.

Comments: 9 pages, 1 figure, in press at ApJ

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

Authors: Kristopher G. Klein, Jean C. Perez, Daniel Verscharen, Alfred Mallet, Benjamin D. G. Chandran

Abstract: 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… ▽ More 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. △ Less

Submitted 18 February, 2015; v1 submitted 11 December, 2014; originally announced December 2014.

Comments: 5 pages, 1 figure, accepted in ApJ Letter

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

Authors: Daniel Verscharen, Benjamin D. G. Chandran, Sofiane Bourouaine, Joseph V. Hollweg

Abstract: 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 thi… ▽ More 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. △ Less

Submitted 16 June, 2015; v1 submitted 17 November, 2014; originally announced November 2014.

Comments: 14 pages, 10 figures, submitted to Astrophys. J

Journal ref: Astrophys. J. 806, 157, 2015

Refined critical balance in strong Alfvenic turbulence

Authors: A. Mallet, A. A. Schekochihin, B. D. G. Chandran

Abstract: 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… ▽ More 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. △ Less

Submitted 24 August, 2015; v1 submitted 21 June, 2014; originally announced June 2014.

Comments: revtex, 5 pages, 5 figures; replaced to match published version

Journal ref: MNRAS 449, L77 (2015)

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

Authors: Peera Pongkitiwanichakul, Benjamin D. G. Chandran

Abstract: 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 effec… ▽ More 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. △ Less

Submitted 3 June, 2014; originally announced June 2014.

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

Authors: Joseph V. Hollweg, Daniel Verscharen, Benjamin D. G. Chandran

Abstract: 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… ▽ More 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$. △ Less

Submitted 17 April, 2014; originally announced April 2014.

Comments: 35 pages, 5 figures, accepted for publication in Astrophys. J

Intermittency and Alignment in Strong RMHD Turbulence

Authors: Benjamin D. G. Chandran, Alexander A. Schekochihin, Alfred Mallet

Abstract: 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 parame… ▽ More 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. △ Less

Submitted 25 March, 2014; originally announced March 2014.

Comments: 12 pages, 4 figures, submitted to ApJ