Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Extended Cyclotron Resonant Heating of the Turbulent Solar Wind

Trevor A. Bowen tbowen@berkeley.edu Space Sciences Laboratory, University of California, Berkeley, CA 94720-7450, USA Ivan Vasko William B. Hanson Center for Space Sciences, University of Texas at Dallas, Richardson, TX, USA Stuart D. Bale Space Sciences Laboratory, University of California, Berkeley, CA 94720-7450, USA Physics Department, University of California, Berkeley, CA 94720-7300, USA Benjamin D. G. Chandran Department of Physics & Astronomy, University of New Hampshire, Durham, NH 03824, USA Alexandros Chasapis Laboratory for Atmospheric and Space Physics, University of Colorado, Boulder, 80303, CO, USA Thierry Dudok de Wit LPC2E, CNRS and University of Orléans, Orléans, France ISSI, Bern, Switzerland Alfred Mallet Space Sciences Laboratory, University of California, Berkeley, CA 94720-7450, USA Michael McManus Space Sciences Laboratory, University of California, Berkeley, CA 94720-7450, USA Romain Meyrand Department of Physics, University of Otago, 730 Cumberland St., Dunedin 9016, New Zealand Marc Pulupa Space Sciences Laboratory, University of California, Berkeley, CA 94720-7450, USA Jonathan Squire Department of Physics, University of Otago, 730 Cumberland St., Dunedin 9016, New Zealand
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-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 15R55Rsimilar-toabsent15subscript𝑅direct-product55subscript𝑅direct-product\sim 15R_{\odot}-55R_{\odot}∼ 15 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT - 55 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. 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 Kármán 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.

1 Introduction

Weakly collisional plasmas, which are common in astrophysical environments, are fundamentally governed by kinetic processes (Marsch, 2006). Our understanding of kinetic processes responsible for turbulent dissipation, heating, and energy transfer in collisionless environments is relatively incomplete, and necessary to explain phenomena such as solar wind acceleration and coronal heating (Parker, 1958; Richardson et al., 1995; Hellinger et al., 2013; Fox et al., 2016), and analogous astrophysical processes.

Recent work on the near-Sun solar wind has highlighted the significant presence of circularly polarized ion-scale waves (Bale et al., 2019; Bowen et al., 2020a), and their association with non-thermal features in particle distributions (Verniero et al., 2020; Klein et al., 2021; Verniero et al., 2022). These waves are characterized by their quasi-parallel propagation along the magnetic field (Jian et al., 2014; Boardsen et al., 2015; Bowen et al., 2020a; Liu et al., 2023). Electric field measurements suggest that these wave predominantly propagate outward from the sun (Bowen et al., 2020b), similar to the observed propagation direction of larger-scale Alfvénic turbulent fluctuations (Roberts et al., 1987; Tu & Marsch, 1995; Bavassano et al., 1998; McManus et al., 2020).

While ion-scale waves are often associated with processes related to kinetic plasma distributions, e.g. instabilities and resonant damping, (Gary, 1993; Isenberg & Lee, 1996; Hollweg & Markovskii, 2002; Marsch, 2006; Klein et al., 2018, 2021), they may also serve as a mechanism to transfer turbulent energy via cyclotron resonance to particle thermal motion (Hollweg & Johnson, 1988; Tu & Marsch, 1997; Cranmer, 2000; Hollweg & Isenberg, 2002; Cranmer, 2014). Signatures of quasilinear cyclotron resonance in the observed proton distribution functions have been suggested in various spacecraft observations (Marsch & Tu, 2001a; He et al., 2015; Verniero et al., 2022; Bowen et al., 2022). Furthermore, steepening of turbulent spectra at ion cyclotron resonant scales has been interpreted as a signature of cyclotron resonant damping (Denskat et al., 1983; Woodham et al., 2018; Lotz et al., 2023) and is correlated with the presence of circularly polarized signatures (Goldstein et al., 1994; Leamon et al., 1998; He et al., 2011; Lion et al., 2016; Zhao et al., 2021). These helical signatures are both correlated with proton temperature anisotropy as well as the turbulent amplitudes (Telloni et al., 2019), suggesting they may be associated with turbulent dissipation.

The quasilinear damping of the in situ population of ion cyclotron waves (ICW) has been measured with a heating rate accounting for 10-20% of the turbulent energy flux (Bowen et al., 2022). Further studies have shown the direct transfer of energy from waves to protons using wave-particle correlation methods (Vech et al., 2020; Luo et al., 2022). The presence of these waves has been observed to correlate both to turbulent features such as the large-scale cross helicity and sub-ion-scale intermittency (Bowen et al., 2023).

In this Letter, we study an extended stream of fast solar wind observed by PSP from 15similar-toabsent15\sim 15∼ 15 to 55R55subscript𝑅direct-product55R_{\odot}55 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT with persistent signatures of left-hand polarized waves. We show that left-handed polarization directly corresponds to turbulent steepening in the ion-kinetic scale transition range (Sahraoui et al., 2009; Kiyani et al., 2009; Bowen et al., 2020c; Duan et al., 2020, 2021), which has historically been interpreted as a signature of dissipation (Denskat et al., 1983; Goldstein et al., 1994; Leamon et al., 1998; Smith et al., 2006, 2012; Lion et al., 2016; Bowen et al., 2020c, 2023). Using the cold plasma dispersion to determine the internal energy and Poynting flux of the measured wave spectrum suggests that the energy to generate waves is stored within the turbulent fluctuations, indicating that the waves are a pathway to ion-scale turbulent dissipation. By applying drifting bi-Maxwellian fits to the distribution function, we estimate the empirical quasilinear heating rate of the local ion-scale waves (Kennel & Engelmann, 1966; Isenberg & Lee, 1996; Bowen et al., 2022). We similarly estimate local turbulent dissipation rates via turbulent amplitudes and the Von Kármán decay law to determine the decay of the largest, outer-scale turbulent fluctuations (Hossain et al., 1995; Wan et al., 2012; Bandyopadhyay et al., 2020; Wu et al., 2022). We demonstrate strong correlations between the quasilinear heating and turbulent dissipation rates. These strong correlations are evident in both global scaling as well as local fluctuations in the turbulent cascade and quasilinear heating rates. These results suggest that cyclotron heating plays an important role in extended solar wind heating. These results may have strong implications for the nature of ion-scale heating in the corona and other weakly-collisonal astrophysical plasmas.

Refer to caption
Figure 1: a) Measured spectral densities from Nov 16-20, 2021 computed for 6737 intervals of 128 s. Every 40th spectra is plotted. Colors represent radial distance of PSP from the sun, ranging from 1455similar-toabsent1455\sim 14-55∼ 14 - 55 Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT, with light colors representing small distances becoming darker to indicate increasing distance. b) Local spectral index α𝛼\alphaitalic_α of measured spectral densities computed in moving window. c) Three power-law fit (3PL) to observed spectra for each interval. d) Local spectral index α𝛼\alphaitalic_α of 3PL spectral densities computed in moving window.

2 Methods & Results

Refer to caption
Figure 2: a) Frequencies corresponding to ion-kinetic scales computed with the Taylor Hypothesis for 6737 intervals as function of solar radius, ion-inertial scale fdisubscript𝑓subscript𝑑𝑖f_{d_{i}}italic_f start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT (red), ion gyroscale fρisubscript𝑓subscript𝜌𝑖f_{\rho_{i}}italic_f start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT (orange), cyclotron resonant scale fΩsubscript𝑓Ωf_{\Omega}italic_f start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT (blue). Additionally a correction to the Taylor Hypothesis approximation to the cyclotron scale fΩsubscriptsuperscript𝑓Ωf^{\prime}_{\Omega}italic_f start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT (green). Break points from spectral fits are also shown fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT (black) and fTKsubscript𝑓𝑇𝐾f_{TK}italic_f start_POSTSUBSCRIPT italic_T italic_K end_POSTSUBSCRIPT (teal). The gyrofrequency Ω/2πΩ2𝜋\Omega/2\piroman_Ω / 2 italic_π is shown in purple. b) Histograms show the distribution of frequency ratios fdi/fITsubscript𝑓subscript𝑑𝑖subscript𝑓𝐼𝑇f_{{d}_{i}}/f_{IT}italic_f start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT / italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT (red), fΩ/fITsubscript𝑓Ωsubscript𝑓𝐼𝑇f_{\Omega}/f_{IT}italic_f start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT / italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT (blue), fΩ/fITsubscriptsuperscript𝑓Ωsubscript𝑓𝐼𝑇f^{\prime}_{\Omega}/f_{IT}italic_f start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT / italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT (green) and fρi/fTKsubscript𝑓subscript𝜌𝑖subscript𝑓𝑇𝐾f_{{\rho}_{i}}/f_{TK}italic_f start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT / italic_f start_POSTSUBSCRIPT italic_T italic_K end_POSTSUBSCRIPT (orange).
Refer to caption
Figure 3: a) Fractional left-handed power B~LHσ2/B~total2subscriptsuperscript~𝐵2𝐿𝐻𝜎subscriptsuperscript~𝐵2𝑡𝑜𝑡𝑎𝑙\tilde{B}^{2}_{LH\sigma}/\tilde{B}^{2}_{total}over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT / over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t italic_o italic_t italic_a italic_l end_POSTSUBSCRIPT computed from wavelet-transform for each of the 6737 intervals as a function of kdisubscript𝑘parallel-tosubscript𝑑𝑖k_{\parallel}d_{i}italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT. b) Spectrogram of left-handed circularly polarized power as function of radial distance; fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT and fTKsubscript𝑓𝑇𝐾f_{TK}italic_f start_POSTSUBSCRIPT italic_T italic_K end_POSTSUBSCRIPT for each measurement are plotted in white, as well as the cyclotron scale fΩsubscript𝑓Ωf_{\Omega}italic_f start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT in green; teal lines show frequencies corresponding to spacecraft reaction wheels (Bowen et al., 2020a). c) internal energy UBsubscript𝑈𝐵U_{B}italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT and d) Poynting flux of waves in a heliocentric inertial frame SHCIsubscript𝑆HCIS_{\rm HCI}italic_S start_POSTSUBSCRIPT roman_HCI end_POSTSUBSCRIPT computed from the cold-plasma dispersion relation.

We focus on a stream of fast solar wind observed in the inner heliosphere by Parker Solar Probe (PSP) from 11/16/2021 to 11/20/2021. During this interval PSP was flying inwards towards the sun in near co-rotation with the solar surface, sampling a relatively singular source region over the range of 14-55 RR\odotitalic_R ⊙ (Badman et al., 2023; Davis et al., 2023). We use measurements from the PSP Solar wind Electron Alpha and Proton experiment’s Solar Probe Analyzer (Livi et al., 2022, SPAN). The PSP FIELDS experiment (Bale et al., 2016) provides measurements of the magnetic field from DC to sub-ion kinetic scales using merged search coil (Jannet et al., 2021) and fluxgate magnetometer (SCaM) measurements (Bowen et al., 2020d).

We break the 4 day stream into a set of intervals of 128 s with 50% overlap to study the evolution of the turbulent spectra alongside ion-scale waves. SCaM data are sensitive well into the sub-ion kinetic scales (Dudok de Wit et al., 2022), but are only available for two axes; as a result, the SCaM data are only used to measure the shape of the spectra. The fluxgate magnetometer data provides three component magnetic field measurements that enable comprehensive study of the properties of observed waves. The fluxgate magnetometer and SCaM data are re-sampled to a uniform 146.4845 Sa/sec rate. Data were discarded when discontinuities or changes in instrumental modes led to artifacts. In total 6737 intervals of 128 s were analyzed.

Turbulent Spectra

We compute spectral densities E(f)𝐸𝑓E(f)italic_E ( italic_f ), with units nT2/Hz, from the two-component SCaM data through ensemble averaging over 4096 point FFTs in each 128s interval. Each spectral density is interpolated onto 320 logarithmically spaced frequencies. Figure 1(a) shows spectral densities computed from the strem intervals, though for clarity only every 40th spectra is plotted. The colors corresponding to radial distance from the sun, which is between 14.1-55.0 Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT; lighter yellow colors are closer to the Sun, whereas darker colors correspond to further heliocentric distances.

An estimate of the local spectral index for each frequency is obtained by performing a linear least-square fit of log10E(f)subscriptlog10𝐸𝑓\text{log}_{10}E(f)log start_POSTSUBSCRIPT 10 end_POSTSUBSCRIPT italic_E ( italic_f ) in a moving window consisting of 60 neighboring logarithmically spaced frequencies. The slope of a linear least-square fit in log space provides a measure of the local spectral index α(f)𝛼𝑓\alpha(f)italic_α ( italic_f ). Figure 1(b) shows the locally measured α𝛼\alphaitalic_α as a function of frequency for spectra in Figure 1(a). The observed spectra are consistent with three-power-law spectra often reported in the solar wind (Sahraoui et al., 2009; Alexandrova et al., 2008; Bowen et al., 2020c): each spectra has a characteristic low-frequency inertial range scaling αIsubscript𝛼𝐼\alpha_{I}italic_α start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT, a transition range with steepened index αTsubscript𝛼𝑇\alpha_{T}italic_α start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT, followed by flattening to an index αKsubscript𝛼𝐾\alpha_{K}italic_α start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT at higher frequencies. The shape of the spectra is consistent but shifts to higher frequencies closer to the Sun. The index of the steep transition range spectra varies, and can have values of up to -10. The extremely steep transition range indices in this stream are predominantly associated with parallel spectra as shown in Duan et al. (2021). As we demonstrate, this steepening is also associated with significant ICW populations (Bowen et al., 2020a, 2023).

To determine the slope of the transition range αTsubscript𝛼𝑇\alpha_{T}italic_α start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT we implement a piecewise, three power-law (3PL) fit, following methods developed in Bowen et al. (2020c). Each of the three ranges is modeled as a power-law spectra with indices αIsubscript𝛼𝐼\alpha_{I}italic_α start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT, αTsubscript𝛼𝑇\alpha_{T}italic_α start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT, and αKsubscript𝛼𝐾\alpha_{K}italic_α start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT separated by break-point frequencies fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT and fTKsubscript𝑓𝑇𝐾f_{TK}italic_f start_POSTSUBSCRIPT italic_T italic_K end_POSTSUBSCRIPT, referring to the respective inertial-transition and transition-kinetic breaks. Fig. 1(c) shows the 3PL fits, corresponding to the spectra in Figure 1(a). Additionally, we apply the procedure to measure the local spectral index α(f)𝛼𝑓\alpha(f)italic_α ( italic_f ) to the 3PL-fits, with the results shown in Figure 1(d). There is good agreement between the observations and the 3PL fits.

Break Scales & Cyclotron Resonance

Refer to caption
Figure 4: a) Transition range spectral index αTsubscript𝛼𝑇\alpha_{T}italic_α start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT against maximum circular polarization max[B~LHσ2/B~total2]maxdelimited-[]subscriptsuperscript~𝐵2𝐿𝐻𝜎subscriptsuperscript~𝐵2𝑡𝑜𝑡𝑎𝑙\text{max}[\tilde{B}^{2}_{LH\sigma}/\tilde{B}^{2}_{total}]max [ over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT / over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t italic_o italic_t italic_a italic_l end_POSTSUBSCRIPT ]. b) Total internal energy of cyclotron waves UBsubscript𝑈𝐵U_{B}italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT plotted against αTsubscript𝛼𝑇\alpha_{T}italic_α start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT . c) Total internal energy of cyclotron waves UBsubscript𝑈𝐵U_{B}italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT, which is bandwidth limited and filtered by left hand polarization at ion-kinetic scales, plotted against total rms amplitude of magnetic field dBrms2𝑑subscriptsuperscript𝐵2𝑟𝑚𝑠dB^{2}_{rms}italic_d italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_r italic_m italic_s end_POSTSUBSCRIPT computed at 128 seconds. Data are colored by Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT with larger heliocentric distances in darker shades. Solid lines show the linear regression for 1020R1020subscript𝑅direct-product10-20R_{\odot}10 - 20 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT (yellow), 2030R2030subscript𝑅direct-product20-30R_{\odot}20 - 30 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT (orange), 3040R3040subscript𝑅direct-product30-40R_{\odot}30 - 40 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT (red), and 4050R4050subscript𝑅direct-product40-50R_{\odot}40 - 50 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT (purple), darker colors correspond to larger distances.

We compare the measured break scales fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT and fTKsubscript𝑓𝑇𝐾f_{TK}italic_f start_POSTSUBSCRIPT italic_T italic_K end_POSTSUBSCRIPT to various physical ion-kinetic scales: the ion gyroradius, ρi=vthi/Ωisubscript𝜌𝑖subscript𝑣𝑡𝑖subscriptΩ𝑖\rho_{i}=v_{thi}/\Omega_{i}italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_v start_POSTSUBSCRIPT italic_t italic_h italic_i end_POSTSUBSCRIPT / roman_Ω start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, where the ion thermal speed is vthi=2kBTi/misubscript𝑣𝑡𝑖2subscript𝑘𝐵subscript𝑇𝑖subscript𝑚𝑖v_{thi}=\sqrt{2k_{B}T_{i}/m_{i}}italic_v start_POSTSUBSCRIPT italic_t italic_h italic_i end_POSTSUBSCRIPT = square-root start_ARG 2 italic_k start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG and the ion gyrofrequency is Ωi=eB0/misubscriptΩ𝑖𝑒subscript𝐵0subscript𝑚𝑖\Omega_{i}=eB_{0}/m_{i}roman_Ω start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_e italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT; and the ion inertial scale di=ρi/βisubscript𝑑𝑖subscript𝜌𝑖subscript𝛽𝑖d_{i}=\rho_{i}/\sqrt{\beta_{i}}italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / square-root start_ARG italic_β start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG, where βi=2n0μ0Ti/B02subscript𝛽𝑖2subscript𝑛0subscript𝜇0subscript𝑇𝑖superscriptsubscript𝐵02\beta_{i}=2n_{0}\mu_{0}T_{i}/B_{0}^{2}italic_β start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = 2 italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. The fundamental charge is given as e𝑒eitalic_e, μ0subscript𝜇0\mu_{0}italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the permeability of free space, Tisubscript𝑇𝑖T_{i}italic_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is the proton temperature, B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the average background magnetic field, n0subscript𝑛0n_{0}italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the background number density. The wavenumbers corresponding to the ion-inertial and ion gyroscale are kρi=ρi1subscript𝑘subscript𝜌𝑖superscriptsubscript𝜌𝑖1k_{\rho_{i}}={\rho_{i}}^{-1}italic_k start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT = italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT and kdi=di1subscript𝑘subscript𝑑𝑖superscriptsubscript𝑑𝑖1k_{d_{i}}=d_{i}^{-1}italic_k start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT = italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. We approximate the wavenumber corresponding to cyclotron resonant interactions as kΩ=Ωi/(vA+vthi)subscript𝑘ΩsubscriptΩ𝑖subscript𝑣𝐴subscript𝑣𝑡𝑖k_{\Omega}=\Omega_{i}/(v_{A}+v_{thi})italic_k start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT = roman_Ω start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / ( italic_v start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT + italic_v start_POSTSUBSCRIPT italic_t italic_h italic_i end_POSTSUBSCRIPT ), where vAsubscript𝑣𝐴v_{A}italic_v start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT is the Alfvén speed vA=B0/mu0min0subscript𝑣𝐴subscript𝐵0𝑚subscript𝑢0subscript𝑚𝑖subscript𝑛0v_{A}=B_{0}/\sqrt{mu_{0}m_{i}n_{0}}italic_v start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / square-root start_ARG italic_m italic_u start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG (Leamon et al., 1998; Wicks et al., 2016; Woodham et al., 2018; Bowen et al., 2020a). This approximation for kΩsubscript𝑘Ωk_{\Omega}italic_k start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT is derived via the cyclotron resonance condition between outward going ICW and protons flowing inwards in the plasma frame, ω+kvth=Ωi𝜔subscript𝑘parallel-tosubscript𝑣𝑡subscriptΩ𝑖\omega+{k}_{\parallel}v_{th}=\Omega_{i}italic_ω + italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_t italic_h end_POSTSUBSCRIPT = roman_Ω start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT (Leamon et al., 1998), and a low frequency limit ω=kvA𝜔subscript𝑘parallel-tosubscript𝑣𝐴\omega=k_{\parallel}v_{A}italic_ω = italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT to the ion-cyclotron dispersion for parallel propagating waves. The cold-plasma dispersion relation is

(ωICWΩi)2=[kdi2(k2di2+4kdi)]2.superscriptsubscript𝜔𝐼𝐶𝑊subscriptΩ𝑖2superscriptdelimited-[]subscript𝑘parallel-tosubscript𝑑𝑖2superscriptsubscript𝑘parallel-to2superscriptsubscript𝑑𝑖24subscript𝑘parallel-tosubscript𝑑𝑖2\left(\frac{\omega_{ICW}}{\Omega_{i}}\right)^{2}=\left[\frac{k_{\parallel}d_{i% }}{2}\left(\sqrt{k_{\parallel}^{2}d_{i}^{2}+4}-k_{\parallel}d_{i}\right)\right% ]^{2}.( divide start_ARG italic_ω start_POSTSUBSCRIPT italic_I italic_C italic_W end_POSTSUBSCRIPT end_ARG start_ARG roman_Ω start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = [ divide start_ARG italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ( square-root start_ARG italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 end_ARG - italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (1)

For the ion-scales kdisubscript𝑘subscript𝑑𝑖k_{d_{i}}italic_k start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT, kρisubscript𝑘subscript𝜌𝑖k_{\rho_{i}}italic_k start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT, and kΩsubscript𝑘Ωk_{\Omega}italic_k start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT, each wavenumber is converted to an effective spacecraft frequency f𝑓fitalic_f first using the Taylor hypothesis, 2πfkvsw2𝜋𝑓𝑘subscript𝑣𝑠𝑤2\pi f\approx{k}{v}_{sw}2 italic_π italic_f ≈ italic_k italic_v start_POSTSUBSCRIPT italic_s italic_w end_POSTSUBSCRIPT to associate ion-kinetic scales with spacecraft frequencies fdisubscript𝑓subscript𝑑𝑖f_{{d}_{i}}italic_f start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT,fρisubscript𝑓subscript𝜌𝑖f_{{\rho}_{i}}italic_f start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT, and fΩsubscript𝑓Ωf_{\Omega}italic_f start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT. For the cyclotron-resonant scale, kΩsubscript𝑘Ωk_{\Omega}italic_k start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT, we also consider a correction to the Taylor hypothesis made by incorporating the Doppler shift equation

2πf=ω(k)+𝐤𝐯sw,2𝜋𝑓𝜔𝑘𝐤subscript𝐯𝑠𝑤2\pi f=\omega(k)+\mathbf{k}\cdot{\mathbf{v}}_{sw},2 italic_π italic_f = italic_ω ( italic_k ) + bold_k ⋅ bold_v start_POSTSUBSCRIPT italic_s italic_w end_POSTSUBSCRIPT , (2)

to improve on our estimate of the frequency corresponding to cyclotron resonant interactions as

2πfΩ=ωICW(kΩ)+kΩvsw,2𝜋subscriptsuperscript𝑓Ωsubscript𝜔𝐼𝐶𝑊subscript𝑘Ωsubscript𝑘Ωsubscript𝑣𝑠𝑤2\pi f^{\prime}_{\Omega}=\omega_{ICW}(k_{\Omega})+{k_{\Omega}}{{v}}_{sw},2 italic_π italic_f start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT = italic_ω start_POSTSUBSCRIPT italic_I italic_C italic_W end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT ) + italic_k start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_s italic_w end_POSTSUBSCRIPT , (3)

which includes the contribution from ωICWsubscript𝜔𝐼𝐶𝑊\omega_{ICW}italic_ω start_POSTSUBSCRIPT italic_I italic_C italic_W end_POSTSUBSCRIPT evaluated at kΩsubscript𝑘Ωk_{\Omega}italic_k start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT. For calculating fΩsubscriptsuperscript𝑓Ωf^{\prime}_{\Omega}italic_f start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT, we assume that the angle between the mean magnetic field and flow direction, θBVsubscript𝜃𝐵𝑉\theta_{BV}italic_θ start_POSTSUBSCRIPT italic_B italic_V end_POSTSUBSCRIPT, is small such that cosθBV1𝑐𝑜𝑠subscript𝜃𝐵𝑉1cos\theta_{BV}\approx 1italic_c italic_o italic_s italic_θ start_POSTSUBSCRIPT italic_B italic_V end_POSTSUBSCRIPT ≈ 1 and thus 𝐤𝐯swkΩvsw𝐤subscript𝐯𝑠𝑤subscript𝑘Ωsubscript𝑣𝑠𝑤\mathbf{k}\cdot{\mathbf{v}}_{sw}\approx{k_{\Omega}}{{v}}_{sw}bold_k ⋅ bold_v start_POSTSUBSCRIPT italic_s italic_w end_POSTSUBSCRIPT ≈ italic_k start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_s italic_w end_POSTSUBSCRIPT for parallel propagating ICWs.

Fig. 2(a) shows the computed kinetic-scale frequencies alongside the measured breaks from the 3PL fitting algorithm against Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. Fig. 2(b) shows distributions ratios of fdi/fITsubscript𝑓subscript𝑑𝑖subscript𝑓𝐼𝑇f_{{d}_{i}}/f_{IT}italic_f start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT / italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT fΩ/fITsubscript𝑓Ωsubscript𝑓𝐼𝑇f_{\Omega}/f_{IT}italic_f start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT / italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT, fΩ/fITsubscriptsuperscript𝑓Ωsubscript𝑓𝐼𝑇f^{\prime}_{\Omega}/f_{IT}italic_f start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT / italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT and fρi/fTKsubscript𝑓subscript𝜌𝑖subscript𝑓𝑇𝐾f_{{\rho}_{i}}/f_{TK}italic_f start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT / italic_f start_POSTSUBSCRIPT italic_T italic_K end_POSTSUBSCRIPT. These distributions provide a measure of agreement between break frequencies and the measured ion kinetic scales. Recent work has suggested that the cyclotron scale fΩsubscript𝑓Ωf_{\Omega}italic_f start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT corresponds well with fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT (Woodham et al., 2018; Vech et al., 2018; Duan et al., 2020; Lotz et al., 2023), here we find that fΩsubscript𝑓Ωf_{\Omega}italic_f start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT is statistically lower than fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT; however, the corrected fΩsubscriptsuperscript𝑓Ωf^{\prime}_{\Omega}italic_f start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT frequency well approximates the fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT break. The fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT transition range break additionally approximates fdisubscript𝑓subscript𝑑𝑖f_{{d}_{i}}italic_f start_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT, such that we cannot distinguish the inertial scale from the cyclotron resonant scale. We find that fρisubscript𝑓subscript𝜌𝑖f_{{\rho}_{i}}italic_f start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT is within the transition range, which is consistent with previous results (Bowen et al., 2020c).

Importantly, fΩsubscript𝑓Ωf_{\Omega}italic_f start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT and fΩsubscriptsuperscript𝑓Ωf^{\prime}_{\Omega}italic_f start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT are estimates to a single cyclotron resonant scale corresponding to resonance between particles at the thermal speed vthsubscript𝑣𝑡v_{th}italic_v start_POSTSUBSCRIPT italic_t italic_h end_POSTSUBSCRIPT, with low-frequency ICWs, ω=kvA𝜔𝑘subscript𝑣𝐴\omega=kv_{A}italic_ω = italic_k italic_v start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT. However, observations from the solar wind suggest that ICWs occur over a range of frequencies corresponding to kdi1𝑘subscript𝑑𝑖1kd_{i}\approx 1italic_k italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ≈ 1 (Bowen et al., 2020a). Empirical determination of circular polarization as a function of spacecraft frequency along with the cold plasma dispersion for parallel propagating ICWs Eq. (1) provides information regarding the wavenumbers at which cyclotron resonant waves occur (Bowen et al., 2020b).

A 64-scale Morlet wavelet transform is applied to each interval to compute a power spectral density B~2(f,t)superscript~𝐵2𝑓𝑡\tilde{B}^{2}(f,t)over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_f , italic_t ) with units nT2/Hz (Farge, 1992; Dudok de Wit et al., 2013; Bowen et al., 2020a). We extract circular polarization of the field via the reduced magnetic-helicity (Howes & Quataert, 2010),

σB(f,t)=2Im(B~1B~2)/(B~12+B~22).subscript𝜎𝐵𝑓𝑡2Imsubscript~𝐵perpendicular-toabsent1superscriptsubscript~𝐵perpendicular-toabsent2superscriptsubscript~𝐵perpendicular-toabsent12superscriptsubscript~𝐵perpendicular-toabsent22\displaystyle\sigma_{B}(f,t)=-2\text{Im}(\tilde{B}_{\perp 1}\tilde{B}_{\perp 2% }^{*})/(\tilde{B}_{\perp 1}^{2}+\tilde{B}_{\perp 2}^{2}).italic_σ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_f , italic_t ) = - 2 Im ( over~ start_ARG italic_B end_ARG start_POSTSUBSCRIPT ⟂ 1 end_POSTSUBSCRIPT over~ start_ARG italic_B end_ARG start_POSTSUBSCRIPT ⟂ 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ) / ( over~ start_ARG italic_B end_ARG start_POSTSUBSCRIPT ⟂ 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + over~ start_ARG italic_B end_ARG start_POSTSUBSCRIPT ⟂ 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) . (4)

of the wavelet coefficients perpendicular to the mean field with left/right handed polarization represented by positive/negative σBsubscript𝜎𝐵\sigma_{B}italic_σ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT. We calculate the left-handed polarized power spectra B~LHσ2subscriptsuperscript~𝐵2𝐿𝐻𝜎\tilde{B}^{2}_{LH\sigma}over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT by filtering out power with σB<0.9subscript𝜎𝐵0.9\sigma_{B}<0.9italic_σ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT < 0.9. The normalized fractional left handed spectra are computed as B~LHσ2/B~total2subscriptsuperscript~𝐵2𝐿𝐻𝜎subscriptsuperscript~𝐵2𝑡𝑜𝑡𝑎𝑙\tilde{B}^{2}_{LH\sigma}/\tilde{B}^{2}_{total}over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT / over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t italic_o italic_t italic_a italic_l end_POSTSUBSCRIPT. Very little right handed polarization is present in this stream.

Based on previous measurements of the phase speed of circularly polarized waves (Bowen et al., 2020b), which show a strong statistical preference for outward-propagation, we can assume that the waves are outward-propagating ICWs, travelling parallel to the mean field. The wavelet transform gives circular polarization as a function of the spacecraft frame frequency, f𝑓fitalic_f. The use of the Doppler shift Eq. (2) in the parallel-propagating limit gives

2πf=ωICW±kvswcosθBV,2𝜋𝑓plus-or-minussubscript𝜔𝐼𝐶𝑊subscript𝑘parallel-tosubscript𝑣𝑠𝑤𝑐𝑜𝑠subscript𝜃𝐵𝑉2\pi f=\omega_{ICW}\pm k_{\parallel}v_{sw}cos\theta_{BV},2 italic_π italic_f = italic_ω start_POSTSUBSCRIPT italic_I italic_C italic_W end_POSTSUBSCRIPT ± italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_s italic_w end_POSTSUBSCRIPT italic_c italic_o italic_s italic_θ start_POSTSUBSCRIPT italic_B italic_V end_POSTSUBSCRIPT , (5)

where we no longer assume either the Taylor Hypothesis or that cosθBV1𝑐𝑜𝑠subscript𝜃𝐵𝑉1cos\theta_{BV}\approx 1italic_c italic_o italic_s italic_θ start_POSTSUBSCRIPT italic_B italic_V end_POSTSUBSCRIPT ≈ 1. The combination of Eq. (5) with Eq. (1) determines the corresponding parallel wave-number ksubscript𝑘parallel-tok_{\parallel}italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT for each spacecraft frequency f𝑓fitalic_f (Bowen et al., 2020b). Figure 3(a) shows the normalized B~LHσ2/B~total2subscriptsuperscript~𝐵2𝐿𝐻𝜎subscriptsuperscript~𝐵2𝑡𝑜𝑡𝑎𝑙\tilde{B}^{2}_{LH\sigma}/\tilde{B}^{2}_{total}over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT / over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t italic_o italic_t italic_a italic_l end_POSTSUBSCRIPT as function of kdisubscript𝑘parallel-tosubscript𝑑𝑖k_{\parallel}d_{i}italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT computed from combining Eq. (1) with Eq. (5). At all radial distances the ICW population appears at kdi=1subscript𝑘parallel-tosubscript𝑑𝑖1k_{\parallel}d_{i}=1italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = 1, where ICWs are inherently dispersive, with a finite bandwidth in wave number ranging from approximately kdi0.34similar-tosubscript𝑘parallel-tosubscript𝑑𝑖0.34k_{\parallel}d_{i}\sim 0.3-4italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∼ 0.3 - 4. The finite bandwidth indicates that, at each time, a range of particle velocities simultaneously satisfy the ICW resonance condition. It is important to note that the spacecraft reaction wheels, highlighted in teal lines in Figure 3(b-d), contribute circularly polarized power at higher kdisubscript𝑘parallel-tosubscript𝑑𝑖k_{\parallel}d_{i}italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT (Bowen et al., 2020a).

Fig. 3(b) shows B~LHσ2/B~total2subscriptsuperscript~𝐵2𝐿𝐻𝜎subscriptsuperscript~𝐵2𝑡𝑜𝑡𝑎𝑙\tilde{B}^{2}_{LH\sigma}/\tilde{B}^{2}_{total}over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT / over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t italic_o italic_t italic_a italic_l end_POSTSUBSCRIPT as a function of spacecraft frequency at each interval studied as a function of solar radial distance. We plot both the fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT and fTKsubscript𝑓𝑇𝐾f_{TK}italic_f start_POSTSUBSCRIPT italic_T italic_K end_POSTSUBSCRIPT breaks. The fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT break agrees very well with the circularly polarized regime, while the fTKsubscript𝑓𝑇𝐾f_{TK}italic_f start_POSTSUBSCRIPT italic_T italic_K end_POSTSUBSCRIPT break bounds the circular polarization at higher frequencies, indicating that the transition range is entirely circularly polarized. We additionally plot the cyclotron scale fΩsubscriptsuperscript𝑓Ωf^{\prime}_{\Omega}italic_f start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT, which approximates the fITsubscript𝑓𝐼𝑇f_{IT}italic_f start_POSTSUBSCRIPT italic_I italic_T end_POSTSUBSCRIPT break.

Cyclotron Waves and Turbulent Dissipation

Using the cold-plasma dispersion and the Doppler shift Eqs. (1) & (2) and our observations of B~LHσ2(kdi)subscriptsuperscript~𝐵2𝐿𝐻𝜎subscript𝑘parallel-tosubscript𝑑𝑖\tilde{B}^{2}_{LH\sigma}(k_{\parallel}d_{i})over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) we measure both the internal energy density of the waves UBsubscript𝑈𝐵U_{B}italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT and the magnitude of the Poynting flux computed in an Heliocentric inertial frame (HCI), SHCIsubscript𝑆HCIS_{\rm HCI}italic_S start_POSTSUBSCRIPT roman_HCI end_POSTSUBSCRIPT assuming the ICW propagation is entirely parallel the radial solar wind flow (Karpman, 1974; Shklyar & Matsumoto, 2009):

UB=B~LHσ2μ0vphvg=B~LHσ22μ0(2Ωiω)(Ωiω),subscriptsuperscript𝑈𝐵subscriptsuperscript~𝐵2𝐿𝐻𝜎subscript𝜇0subscript𝑣𝑝subscript𝑣𝑔subscriptsuperscript~𝐵2𝐿𝐻𝜎2subscript𝜇02subscriptΩ𝑖𝜔subscriptΩ𝑖𝜔\displaystyle U^{\prime}_{B}=\frac{\tilde{B}^{2}_{LH\sigma}}{\mu_{0}}\frac{v_{% ph}}{v_{g}}=\frac{\tilde{B}^{2}_{LH\sigma}}{2\mu_{0}}\frac{(2\Omega_{i}-\omega% )}{(\Omega_{i}-\omega)},italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT = divide start_ARG over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG divide start_ARG italic_v start_POSTSUBSCRIPT italic_p italic_h end_POSTSUBSCRIPT end_ARG start_ARG italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT end_ARG = divide start_ARG over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG divide start_ARG ( 2 roman_Ω start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_ω ) end_ARG start_ARG ( roman_Ω start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_ω ) end_ARG , (6)
SHCI=(vsw+vg)UB,superscriptsubscript𝑆HCIsubscript𝑣𝑠𝑤subscript𝑣𝑔superscriptsubscript𝑈𝐵\displaystyle S_{\rm HCI}^{\prime}=(v_{sw}+v_{g})U_{B}^{\prime},italic_S start_POSTSUBSCRIPT roman_HCI end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = ( italic_v start_POSTSUBSCRIPT italic_s italic_w end_POSTSUBSCRIPT + italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ) italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , (7)

where vphsubscript𝑣𝑝v_{ph}italic_v start_POSTSUBSCRIPT italic_p italic_h end_POSTSUBSCRIPT and vgsubscript𝑣𝑔v_{g}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT are the plasma frame phase and group velocities respectively. As the units of B~LHσ2subscriptsuperscript~𝐵2𝐿𝐻𝜎{\tilde{B}^{2}_{LH\sigma}}over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT are nT2/Hz both UBsubscriptsuperscript𝑈𝐵U^{\prime}_{B}italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT and SHCIsuperscriptsubscript𝑆HCIS_{\rm HCI}^{\prime}italic_S start_POSTSUBSCRIPT roman_HCI end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT are defined as spectral densities. The internal energy UBsubscript𝑈𝐵U_{B}italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT and Poynting flux magnitude SHCIsubscript𝑆HCIS_{\rm HCI}italic_S start_POSTSUBSCRIPT roman_HCI end_POSTSUBSCRIPT corresponding to each wavelet coefficient are obtained by multiplying UBsubscriptsuperscript𝑈𝐵U^{\prime}_{B}italic_U start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT and SHCIsuperscriptsubscript𝑆HCIS_{\rm HCI}^{\prime}italic_S start_POSTSUBSCRIPT roman_HCI end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT by the bandwidth of each wavelet Δf=dfdkΔkΔ𝑓𝑑𝑓𝑑subscript𝑘parallel-toΔsubscript𝑘parallel-to\Delta f=\frac{df}{dk_{\parallel}}\Delta k_{\parallel}roman_Δ italic_f = divide start_ARG italic_d italic_f end_ARG start_ARG italic_d italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT end_ARG roman_Δ italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT, where the derivative df/dk𝑑𝑓𝑑subscript𝑘parallel-to{df}/{dk_{\parallel}}italic_d italic_f / italic_d italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT and ΔkΔsubscript𝑘parallel-to\Delta k_{\parallel}roman_Δ italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT are obtained from Eq. 5. Figure 3 (c&d) show UBsubscript𝑈𝐵U_{B}italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT and SHCIsubscript𝑆HCIS_{\rm HCI}italic_S start_POSTSUBSCRIPT roman_HCI end_POSTSUBSCRIPT at each frequency over the interval.

The total Poynting flux and internal energy can be integrated at each time to give total values Utotsubscript𝑈𝑡𝑜𝑡U_{tot}italic_U start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT and Stotsubscript𝑆𝑡𝑜𝑡S_{tot}italic_S start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT for each interval. Due to the sporadic noise that occurs at low frequencies, as well as high frequency contributions to polarization from the spacecraft reaction wheels, we zero out contributions to Stotsubscript𝑆𝑡𝑜𝑡S_{tot}italic_S start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT and Utotsubscript𝑈𝑡𝑜𝑡U_{tot}italic_U start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT where kdi<0.3𝑘subscript𝑑𝑖0.3kd_{i}<0.3italic_k italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT < 0.3 and kdi>5𝑘subscript𝑑𝑖5kd_{i}>5italic_k italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT > 5.

We use Utotsubscript𝑈𝑡𝑜𝑡U_{tot}italic_U start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT and BLHσmax=max[B~LHσ2/B~total2]subscriptsuperscript𝐵max𝐿𝐻𝜎maxdelimited-[]subscriptsuperscript~𝐵2𝐿𝐻𝜎subscriptsuperscript~𝐵2𝑡𝑜𝑡𝑎𝑙B^{\text{max}}_{LH\sigma}=\text{max}[\tilde{B}^{2}_{LH\sigma}/\tilde{B}^{2}_{% total}]italic_B start_POSTSUPERSCRIPT max end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT = max [ over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT / over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t italic_o italic_t italic_a italic_l end_POSTSUBSCRIPT ] as proxies for the ICWs in each interval. Fig. 4 shows how these quantities relate to parameters associated with turbulence, αTsubscript𝛼𝑇\alpha_{T}italic_α start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT and as well as the total RMS turbulent amplitude δB2delimited-⟨⟩𝛿superscript𝐵2\langle\delta B^{2}\rangle⟨ italic_δ italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ⟩ computed in each 128 s interval. We consider only intervals with 150<θBV<170superscript150subscript𝜃𝐵𝑉superscript170150^{\circ}<\theta_{BV}<170^{\circ}150 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT < italic_θ start_POSTSUBSCRIPT italic_B italic_V end_POSTSUBSCRIPT < 170 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT in order to control for effects associated with anisotropy, which affect observational signatures of both the waves and the turbulence (Chen et al., 2010; Bowen et al., 2020a); of the total 6736 intervals, 3172 occur with 150<θBV<170superscript150subscript𝜃𝐵𝑉superscript170150^{\circ}<\theta_{BV}<170^{\circ}150 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT < italic_θ start_POSTSUBSCRIPT italic_B italic_V end_POSTSUBSCRIPT < 170 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT. In each panel of Fig. 4 data is plotted with the color scale in Fig. 1 with lighter colors corresponding to intervals closer to the sun. We compute lines of best fit for data within 1020R1020subscript𝑅direct-product10-20R_{\odot}10 - 20 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT, 2030R2030subscript𝑅direct-product20-30R_{\odot}20 - 30 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT, 3040R3040subscript𝑅direct-product30-40R_{\odot}30 - 40 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT, and 4050R4050subscript𝑅direct-product40-50R_{\odot}40 - 50 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT, which are shown in each panel of Fig. 4. Fig. 4(a) shows that the transition range slope αTsubscript𝛼𝑇\alpha_{T}italic_α start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT is strongly correlated with the level of circular polarization, BLHσmaxsubscriptsuperscript𝐵max𝐿𝐻𝜎B^{\text{max}}_{LH\sigma}italic_B start_POSTSUPERSCRIPT max end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT (Bowen et al., 2023), with similar trends observed at all solar radii. Fig. 4(b) shows that the internal energy of the waves UBsubscript𝑈𝐵U_{B}italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT is anti-correlated with αTsubscript𝛼𝑇\alpha_{T}italic_α start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT, such that steeper slopes contain greater amounts of ICW energy. This correlation is present at all Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT, though at larger distances the maximum measured Utotsubscript𝑈𝑡𝑜𝑡U_{tot}italic_U start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT are similar to the lowest Utotsubscript𝑈𝑡𝑜𝑡U_{tot}italic_U start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT close to the sun, such that constraining the data by RR\odotitalic_R ⊙ is necessary to see the correlations. Fig. 4(c) shows that the internal energy contained within the circularly polarized ion-scale transition range is globally a function of the amplitude of the turbulent fluctuations dBrms2𝑑superscriptsubscript𝐵𝑟𝑚𝑠2dB_{rms}^{2}italic_d italic_B start_POSTSUBSCRIPT italic_r italic_m italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT computed at the 128 s scale, consistent with previous results in Shankarappa et al. (2023).

ICW Heating Rates

Our analysis of the transition range reveals signatures of circularly polarized ion-cyclotron resonant waves; we apply well established principles of energy conservation and quasilinear heating to connect these waves to turbulent dissipation.

The Poynting theorem states that the change in internal energy of the ICW population is

UBt=𝐒QICW+ϵICW,subscript𝑈𝐵𝑡𝐒subscript𝑄ICWsubscriptitalic-ϵICW\frac{\partial U_{B}}{\partial t}=-\nabla\cdot\mathbf{S}-Q_{\rm ICW}+\epsilon_% {\rm ICW},divide start_ARG ∂ italic_U start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_t end_ARG = - ∇ ⋅ bold_S - italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT + italic_ϵ start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT , (8)

where QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT is the dissipation rate of the waves and ϵICWsubscriptitalic-ϵICW\epsilon_{\rm ICW}italic_ϵ start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT is a driving term. Assuming a steady state, we compute the three terms on the right side of Eq. (8). We approximate 𝐒𝐒\nabla\cdot\mathbf{S}∇ ⋅ bold_S in the heliocentric inertial frame as the radial derivative

𝐒1r2rr2Stot.𝐒1superscript𝑟2𝑟superscript𝑟2subscript𝑆𝑡𝑜𝑡\nabla\cdot\mathbf{S}\approx\frac{1}{r^{2}}\frac{\partial}{\partial r}r^{2}S_{% tot}.∇ ⋅ bold_S ≈ divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG ∂ end_ARG start_ARG ∂ italic_r end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_S start_POSTSUBSCRIPT italic_t italic_o italic_t end_POSTSUBSCRIPT . (9)

We estimate QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT using the quasilinear heating rate (Kennel & Engelmann, 1966) with the empirically observed spectrum of cyclotron waves

I(k)=B~LHσ2B02dfdk.𝐼subscript𝑘parallel-tosubscriptsuperscript~𝐵2𝐿𝐻𝜎superscriptsubscript𝐵02𝑑𝑓𝑑subscript𝑘parallel-to\displaystyle I(k_{\parallel})=\frac{\tilde{B}^{2}_{LH\sigma}}{B_{0}^{2}}\frac% {df}{dk_{\parallel}}.italic_I ( italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) = divide start_ARG over~ start_ARG italic_B end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_H italic_σ end_POSTSUBSCRIPT end_ARG start_ARG italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG italic_d italic_f end_ARG start_ARG italic_d italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT end_ARG . (10)

For each SPANi measurement, we perform a drifting bi-Maxwellian fit in the proton core frame with the form

gp(v,v)=ncπ3/2wc,2wc,exp[v2wc,2v2wc,2]+nbπ3/2wb,2wb,exp[v2wb,2(vvD)2wb,2],\displaystyle\begin{split}g_{p}(v_{\perp},v_{\parallel})=&\frac{n_{c}}{\pi^{3/% 2}w_{c,\perp}^{2}w_{c,\parallel}}\text{exp}\left[-\frac{v_{\perp}^{2}}{w_{c,% \perp}^{2}}-\frac{v_{\parallel}^{2}}{w_{c,\parallel}^{2}}\right]\\ &+\frac{n_{b}}{\pi^{3/2}w_{b,\perp}^{2}w_{b,\parallel}}\text{exp}\left[-\frac{% v_{\perp}^{2}}{w_{b,\perp}^{2}}-\frac{(v_{\parallel}-v_{D})^{2}}{w_{b,% \parallel}^{2}}\right],\end{split}start_ROW start_CELL italic_g start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) = end_CELL start_CELL divide start_ARG italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_π start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT italic_w start_POSTSUBSCRIPT italic_c , ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_w start_POSTSUBSCRIPT italic_c , ∥ end_POSTSUBSCRIPT end_ARG exp [ - divide start_ARG italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUBSCRIPT italic_c , ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUBSCRIPT italic_c , ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL + divide start_ARG italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT end_ARG start_ARG italic_π start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT italic_w start_POSTSUBSCRIPT italic_b , ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_w start_POSTSUBSCRIPT italic_b , ∥ end_POSTSUBSCRIPT end_ARG exp [ - divide start_ARG italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUBSCRIPT italic_b , ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT - italic_v start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUBSCRIPT italic_b , ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] , end_CELL end_ROW (11)

which has anisotropic thermal speeds perpendicular and parallel the mean field (perpendicular-to\perp and )\parallel)∥ ) for the beam and core (subscripts, b and c): wc,,wc,,wb,,wb,w_{c,\perp},w_{c,\parallel},w_{b,\perp},w_{b,\parallel}italic_w start_POSTSUBSCRIPT italic_c , ⟂ end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT italic_c , ∥ end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT italic_b , ⟂ end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT italic_b , ∥ end_POSTSUBSCRIPT and a relative drift, vDsubscript𝑣𝐷v_{D}italic_v start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT, parallel to the mean magnetic field (Marsch, 2006; Klein et al., 2021). We compute the average fit parameters over the 128s interval such that an average gyrotropic g¯(v,v)¯𝑔subscript𝑣perpendicular-tosubscript𝑣parallel-to\bar{g}(v_{\perp},v_{\parallel})over¯ start_ARG italic_g end_ARG ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) is computed. Using g¯(v,v)¯𝑔subscript𝑣perpendicular-tosubscript𝑣parallel-to\bar{g}(v_{\perp},v_{\parallel})over¯ start_ARG italic_g end_ARG ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ), and the observed spectrum of cyclotron waves I(k)𝐼subscript𝑘parallel-toI(k_{\parallel})italic_I ( italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ), we use the cold plasma dispersion, Eq. 1, with wave numbers defined from Eq. 5, to determine the volumetric heating rate in the local interval:

QICW=mpv22g¯td3𝐯subscript𝑄ICWsubscript𝑚𝑝superscript𝑣22¯𝑔𝑡superscript𝑑3𝐯\displaystyle Q_{\rm ICW}=\int\frac{m_{p}v^{2}}{2}\frac{\partial\bar{g}}{% \partial t}d^{3}\mathbf{v}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT = ∫ divide start_ARG italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG divide start_ARG ∂ over¯ start_ARG italic_g end_ARG end_ARG start_ARG ∂ italic_t end_ARG italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT bold_v
=πe24mpd3𝐯{v20dk1v\displaystyle=\frac{\pi e^{2}}{4m_{p}}\int d^{3}\mathbf{v}\Bigg{\{}v^{2}\int_{% 0}^{\infty}dk_{\parallel}\frac{1}{v_{\perp}}= divide start_ARG italic_π italic_e start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT end_ARG ∫ italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT bold_v { italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT divide start_ARG 1 end_ARG start_ARG italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT end_ARG G^k[vδ(ωkkvΩp)\displaystyle\hat{G}_{k}\Bigg{[}v_{\perp}\delta(\omega_{k}-k_{\parallel}v_{% \parallel}-\Omega_{p})over^ start_ARG italic_G end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT [ italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT italic_δ ( italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT - italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT - roman_Ω start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT )
×ωk2k2c2I(k)G^kg¯(v,v)]},\displaystyle\times\frac{\omega_{k}^{2}}{k_{\parallel}^{2}c^{2}}I(k_{\parallel% })\hat{G}_{k}\bar{g}(v_{\perp},v_{\parallel})\Bigg{]}\Bigg{\}},× divide start_ARG italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_I ( italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) over^ start_ARG italic_G end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT over¯ start_ARG italic_g end_ARG ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) ] } , (12)

with

G^k=(1kvωk)v+kvωkvsubscript^𝐺𝑘1subscript𝑘parallel-tosubscript𝑣parallel-tosubscript𝜔𝑘subscript𝑣perpendicular-tosubscript𝑘parallel-tosubscript𝑣perpendicular-tosubscript𝜔𝑘subscript𝑣parallel-to\displaystyle\hat{G}_{k}=(1-\frac{k_{\parallel}v_{\parallel}}{\omega_{k}})% \frac{\partial}{\partial{v_{\perp}}}+\frac{k_{\parallel}v_{\perp}}{\omega_{k}}% \frac{\partial}{\partial{v_{\parallel}}}over^ start_ARG italic_G end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = ( 1 - divide start_ARG italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT end_ARG start_ARG italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG ) divide start_ARG ∂ end_ARG start_ARG ∂ italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT end_ARG + divide start_ARG italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT end_ARG start_ARG italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG divide start_ARG ∂ end_ARG start_ARG ∂ italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT end_ARG (13)

(Kennel & Engelmann, 1966). A gyrotropic differential volumetric heating rate Q(v,v)𝑄subscript𝑣perpendicular-tosubscript𝑣parallel-toQ(v_{\perp},v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) can be defined from QICW=Q(v,v)𝑑v𝑑vsubscript𝑄ICW𝑄subscript𝑣perpendicular-tosubscript𝑣parallel-todifferential-dsubscript𝑣parallel-todifferential-dsubscript𝑣perpendicular-toQ_{\rm ICW}=\int Q(v_{\perp},v_{\parallel})dv_{\parallel}dv_{\perp}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT = ∫ italic_Q ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) italic_d italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT. We outline how to calculate QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT and Q(v,v)𝑄subscript𝑣perpendicular-tosubscript𝑣parallel-toQ(v_{\perp},v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) in Appendix A through using the δ𝛿\deltaitalic_δ-function to evaluate the integral over wave number and the subsequent numerical integration of Q(v,v)𝑄subscript𝑣perpendicular-tosubscript𝑣parallel-toQ(v_{\perp},v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) over velocity coordinates.

Figure 5(a-b) shows Q(v)𝑄subscript𝑣parallel-toQ(v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ), the volumetric heating rate as function of resonant parallel velocity, which is computed from integrating Q(v,v)𝑄subscript𝑣perpendicular-tosubscript𝑣parallel-toQ(v_{\perp},v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) over perpendicular velocities, as a function of solar radius. Figure 5(a) shows positive Q(v)Δv𝑄subscript𝑣parallel-toΔsubscript𝑣parallel-toQ(v_{\parallel})\Delta v_{\parallel}italic_Q ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) roman_Δ italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT, at resonant velocities less than the thermal speed (shown in yellow), corresponding to wave-absorption and consequent plasma heating. The ΔvΔsubscript𝑣parallel-to\Delta v_{\parallel}roman_Δ italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT term, which is the resolution of the numerical integration 1 km/s, is included to normalize the differential Q(v)𝑄subscript𝑣parallel-toQ(v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) in terms of a volumetric heating rate W/m3. Figure 5(b) shows wave emission, which, consistent with Bowen et al. (2022), is significantly less energetically relevant than the regions with positive Q(v)Δv𝑄subscript𝑣parallel-toΔsubscript𝑣parallel-toQ(v_{\parallel})\Delta v_{\parallel}italic_Q ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) roman_Δ italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT wave-absorption, which indicates net resonant-heating of the plasma.

Importantly, the use of a biMaxwellian model to approximate the distribution may not capture non-thermal features that resonate with the observed populations of waves (Dum et al., 1980; Viñas & Gurgiolo, 2009). In Appendix B, we implement two well understood non-parametric techniques to model the observed distributions, Hermite polynomials and radial basis functions, to verify the independence of our results from the biMaxwellian model. While there are slight differences in the level of heating predicted by each model, the qualitative interpretation of ICW resonance heating is found in each case and the average values closely follow that of the drifting biMaxwellian.

Connecting ICW Heating to the Turbulent Cascade

The importance of resonant heating can be understood by comparing the net heating rates QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT to the turbulent cascade rate ϵitalic-ϵ\epsilonitalic_ϵ. Following (Bandyopadhyay et al., 2020) we compute the energy cascade rate of the turbulence assuming a von Kármán decay law.

ϵ±=α(δz±)2δzL±,superscriptitalic-ϵplus-or-minus𝛼superscript𝛿superscript𝑧plus-or-minus2𝛿superscript𝑧minus-or-plussuperscript𝐿plus-or-minus\epsilon^{\prime\pm}=\alpha\frac{(\delta z^{\pm})^{2}\delta z^{\mp}}{L^{\pm}},italic_ϵ start_POSTSUPERSCRIPT ′ ± end_POSTSUPERSCRIPT = italic_α divide start_ARG ( italic_δ italic_z start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_δ italic_z start_POSTSUPERSCRIPT ∓ end_POSTSUPERSCRIPT end_ARG start_ARG italic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT end_ARG , (14)

where δz±𝛿superscript𝑧plus-or-minus\delta z^{\pm}italic_δ italic_z start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT are root mean square fluctuation amplitudes of the Elsasser variables δ𝐳±=δ𝐯±δ𝐁/μ0ρ.𝛿superscript𝐳plus-or-minusplus-or-minus𝛿𝐯𝛿𝐁subscript𝜇0𝜌\delta\mathbf{z}^{\pm}=\delta\mathbf{v}\pm\delta\mathbf{B}/\sqrt{\mu_{0}\rho}.italic_δ bold_z start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT = italic_δ bold_v ± italic_δ bold_B / square-root start_ARG italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ρ end_ARG . The correlation length L±superscript𝐿plus-or-minusL^{\pm}italic_L start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT is determined from the solar-wind advected spatial-scale corresponding to the time-lag for which the correlation reduces below a factor e1superscript𝑒1e^{-1}italic_e start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. The coefficient α=0.03𝛼0.03\alpha=0.03italic_α = 0.03 is determined from numerical means (Hossain et al., 1995; Wan et al., 2012; Usmanov et al., 2014; Bandyopadhyay et al., 2018) and chosen to agree with Bandyopadhyay et al. (2020) and Wu et al. (2022). However, this value of α𝛼\alphaitalic_α is derived assuming low cross helicity, while the near-Sun solar wind is characterized by imbalance with σc1subscript𝜎𝑐1\sigma_{c}\approx 1italic_σ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ≈ 1 (McManus et al., 2020). Note the units for ϵsuperscriptitalic-ϵ\epsilon^{\prime}italic_ϵ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT in equation 14 are of J/s/kg; we define ϵitalic-ϵ\epsilonitalic_ϵ with units of W/m3, which are comparable to the volumetric heating rates in Equation 2, through multiplying ϵsuperscriptitalic-ϵ\epsilon^{\prime}italic_ϵ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT by the background mass density ϵ=ϵmin0.italic-ϵsuperscriptitalic-ϵsubscript𝑚𝑖subscript𝑛0\epsilon=\epsilon^{\prime}m_{i}n_{0}.italic_ϵ = italic_ϵ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT .

Measurement of ϵ+superscriptitalic-ϵ\epsilon^{\prime+}italic_ϵ start_POSTSUPERSCRIPT ′ + end_POSTSUPERSCRIPT via Eq. (14), which we assume approximates the total cascade rate ϵsuperscriptitalic-ϵ\epsilon^{\prime}italic_ϵ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, requires the Elsasser amplitude δz±𝛿superscript𝑧plus-or-minus\delta z^{\pm}italic_δ italic_z start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT at the outer, energy containing scales, of the turbulence. Davis et al. (2023) have studied the evolution of the outer scales in this same stream and found that the outer-scale spectral break evolves from 102superscript10210^{-2}10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT Hz close to 14 Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT to 103superscript10310^{-3}10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT Hz at 50Rsimilar-toabsent50subscript𝑅direct-product\sim 50R_{\odot}∼ 50 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. Accordingly, we compute the cascade rate ϵ+superscriptitalic-ϵ\epsilon^{+}italic_ϵ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT in successive 30 minute intervals (corresponding to 5×104absent5superscript104\approx 5\times 10^{-4}≈ 5 × 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT Hz) with 50%percent5050\%50 % overlap. In each 30 minute interval we then compute the average heating rate QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ and 𝐒delimited-⟨⟩𝐒\langle\nabla\cdot\mathbf{S}\rangle⟨ ∇ ⋅ bold_S ⟩.

Figure 5(c) shows ϵitalic-ϵ\epsilonitalic_ϵ, QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ and 𝐒delimited-⟨⟩𝐒\langle\nabla\cdot\mathbf{S}\rangle⟨ ∇ ⋅ bold_S ⟩ as a function of radial distance. A surprisingly good agreement is found between ϵitalic-ϵ\epsilonitalic_ϵ and QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩, while the divergence of the Poynting flux term is significantly smaller. The relatively small 𝐒delimited-⟨⟩𝐒\langle\nabla\cdot\mathbf{S}\rangle⟨ ∇ ⋅ bold_S ⟩ indicates that the wave dynamics are dominated by the source, ϵICWsubscriptitalic-ϵICW\epsilon_{\rm ICW}italic_ϵ start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT, and dissipation, QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT, and that the waves do not propagate significantly. This result is consistent with previous studies arguing for in situ driving of these waves (Bowen et al., 2020a; Vech et al., 2020; Liu et al., 2023).

We find a Spearman, non-parametric, ranked correlation between ϵitalic-ϵ\epsilonitalic_ϵ and QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ of 0.79; we further compute the Pearson-correlation of the logarithmic quantities as 0.88. Both of these values indicate strong correlations. Furthermore, we fit each of ϵitalic-ϵ\epsilonitalic_ϵ and QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ to a power-law in Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT noted as ϵsuperscriptitalic-ϵ\epsilon^{\prime}italic_ϵ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT and QICWsuperscriptsubscript𝑄ICWQ_{\rm ICW}^{\prime}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT. We compute the local variations from the global power-law trends using ϵϵitalic-ϵsuperscriptitalic-ϵ\epsilon-\epsilon^{\prime}italic_ϵ - italic_ϵ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT and QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ -QICWsuperscriptsubscript𝑄ICW\ Q_{\rm ICW}^{\prime}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT as a measure of how closely the local quasilinear heating rate follows the turbulent energy cascade rate. Figure 5(d) shows the local fluctuations QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ -QICWsuperscriptsubscript𝑄ICWQ_{\rm ICW}^{\prime}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT against ϵϵitalic-ϵsuperscriptitalic-ϵ\epsilon-\epsilon^{\prime}italic_ϵ - italic_ϵ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT. We measure a Spearman ranked-correlation of 0.56.

Additionally, we interpolate ϵitalic-ϵ\epsilonitalic_ϵ, measured at a 30 minute cadence, onto the 128 second cadence of the QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT measurements. Figure 5(e) shows a 2D-histogram of QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT against the interpolated ϵitalic-ϵ\epsilonitalic_ϵ; the distribution is column-normalized to the measured cascade rate. Good correlations are obtained, with a Spearman ranked-correlation of 0.58. These results suggest that characteristically, the quasilinear heating rate is proportional to the energy cascade rate with a constant of proportionality near unity. Figure 5(f) shows a 2D-histogram of QICW/ϵsubscript𝑄ICWitalic-ϵQ_{\rm ICW}/\epsilonitalic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT / italic_ϵ, where the interpolated ϵitalic-ϵ\epsilonitalic_ϵ is again used, against solar radius Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT; the distribution is column-normalized to Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT. Implicitly, this suggests that much of the turbulent cascade flux enters the ICW population such that QICW=χϵ+subscript𝑄ICW𝜒superscriptitalic-ϵQ_{\rm ICW}=\chi\epsilon^{+}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT = italic_χ italic_ϵ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT with χ0.11𝜒0.11\chi\approx 0.1-1italic_χ ≈ 0.1 - 1 over the vast majority of the observed stream. There is some radial trend, with 0.1<χ<0.50.1𝜒0.50.1<\chi<0.50.1 < italic_χ < 0.5 at 1520Rabsent1520subscript𝑅direct-product\approx 15-20R_{\odot}≈ 15 - 20 italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT and χ1𝜒1\chi\approx 1italic_χ ≈ 1 at larger distances. Values with χ>1𝜒1\chi>1italic_χ > 1 suggest dissipation in excess of the turbulent cascade rate, we believe this is predominantly due to uncertainty in the measurements of QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT and ϵ.italic-ϵ\epsilon.italic_ϵ .

Refer to caption
Figure 5: a) Volumetric heating rates, Q(v)Δv𝑄subscript𝑣parallel-toΔsubscript𝑣parallel-toQ(v_{\parallel})\Delta v_{\parallel}italic_Q ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) roman_Δ italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT, with Δv=1km/sΔsubscript𝑣parallel-to1km/s\Delta v_{\parallel}=1\textrm{km/s}roman_Δ italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT = 1 km/s, which have been integrated over perpendicular velocity, as a function of resonant parallel velocity and solar radial distance. b) Same as panel a, but for negative values of Q(v)Δv𝑄subscript𝑣parallel-toΔsubscript𝑣parallel-toQ(v_{\parallel})\Delta v_{\parallel}italic_Q ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) roman_Δ italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT. Yellow dots in panels a& b show the measured parallel thermal speed. c) Measured turbulent cascade rate ϵitalic-ϵ\epsilonitalic_ϵ from von Kármán decay law (black), average cyclotron heating rate from quasilinear theory QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ (red), and 𝐒delimited-⟨⟩𝐒\langle\nabla\cdot\mathbf{S}\rangle⟨ ∇ ⋅ bold_S ⟩ for total cyclotron waves (blue). Power-law fits, ϵitalic-ϵ\epsilonitalic_ϵ’ and QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT’, are shown in dashed lines. d) QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ -QICWsuperscriptsubscript𝑄ICWQ_{\rm ICW}^{\prime}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT plotted against ϵϵitalic-ϵsuperscriptitalic-ϵ\epsilon-\epsilon^{\prime}italic_ϵ - italic_ϵ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, where + and \triangle corresponds to data with QICWdelimited-⟨⟩subscript𝑄ICW\langle Q_{\rm ICW}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ -QICWsuperscriptsubscript𝑄ICWQ_{\rm ICW}^{\prime}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT and ϵϵitalic-ϵsuperscriptitalic-ϵ\epsilon-\epsilon^{\prime}italic_ϵ - italic_ϵ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT respectively greater and less than zero. e) Two-dimensional histogram of QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT against ϵitalic-ϵ\epsilonitalic_ϵ interpolated onto 128 second cadence. f) Two-dimensional histogram of QICW/ϵsubscript𝑄ICWitalic-ϵQ_{\rm ICW}/\epsilonitalic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT / italic_ϵ against Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT; dashed black lines show QICW/ϵsubscript𝑄ICWitalic-ϵQ_{\rm ICW}/\epsilonitalic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT / italic_ϵ at levels of 0.1, 0.5, and 1; the average value of QICW/ϵsubscript𝑄ICWitalic-ϵQ_{\rm ICW}/\epsilonitalic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT / italic_ϵ is shown in black; QICW/ϵ>1subscript𝑄ICWitalic-ϵ1Q_{\rm ICW}/\epsilon>1italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT / italic_ϵ > 1 indicating quasilinear heating at levels greater than the turbulent cascade may arise due to uncertainty in the estimates of both QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT and ϵitalic-ϵ\epsilonitalic_ϵ. The distributions in panels e&f are column normalized to the maximum value in each data column.

3 Conclusions

The nature of the steepening of turbulent energy spectrum at ion-kinetic scales has long been debated and is an important signature in understanding collisionless turbulent dissipation (Denskat et al., 1983; Smith et al., 1990; Goldstein et al., 1994; Leamon et al., 1998). Studies of the location of the ion-kinetic scale break frequency suggest that the most likely candidate scale corresponds to resonance between outward going ICWs and the thermal ion population (Woodham et al., 2018; Vech et al., 2018; Duan et al., 2020; Lotz et al., 2023). Previous observations from 1 AU have suggested that ICW play a role in dissipating solar wind turbulence (Leamon et al., 1998; Lion et al., 2016; Woodham et al., 2019; Telloni et al., 2019; Zhao et al., 2021; Luo et al., 2022). In this Letter, we concretely demonstrate that ion-scale spectral steepening is associated with circular polarization. The break between the inerital and transition ranges follows the regime of ICWs, and the spectral break between the transition range and subion-scale turbulence bounds the circularly polarized waves exactly at higher frequencies. We suggest that the transition range corresponds to cyclotron waves that are associated with turbulent dissipation and ion-scale heating (Bowen et al., 2023). Our results show that internal energy and momentum of the waves is correlated with the amplitude of the turbulent fluctuations, Figure 3, suggesting that the energy contained in, and transported by, the waves originates from the turbulent fluctuations. These results are consistent with earlier observations from 1 AU suggesting that the level of steepening in the turbulent cascade may relate to the level of Alfvénicity and the turbulent cascade rate (Smith et al., 2006; Bruno et al., 2014).

In this Letter, we compute two measures of energy transfer: first, the quasilinear heating rate of kinetic scale ICWs, QICWdelimited-⟨⟩subscript𝑄ICW\langle{Q_{\rm ICW}}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩, and second, the von Kármán turbulent decay rate ϵ+superscriptitalic-ϵ\epsilon^{+}italic_ϵ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT. The measurements of QICWdelimited-⟨⟩subscript𝑄ICW\langle{Q_{\rm ICW}}\rangle⟨ italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT ⟩ are obtained through integrating kinetic phase space densities measured by SPANi, while ϵ+superscriptitalic-ϵ\epsilon^{+}italic_ϵ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT is determined from outer scale turbulent fluctuations. It is striking that these quantitative estimates of energy transfer, ϵ+superscriptitalic-ϵ\epsilon^{+}italic_ϵ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and QICWsubscript𝑄ICW{Q_{\rm ICW}}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT, obtained via entirely different methods at vastly different spatial scales show significant correlations. Such strong correlations suggest that parallel-cyclotron resonance dissipates significant amounts of turbulent energy into the solar-wind ion populations. The correlation between energy flux and quasilinear heating rate is observed in both the global scaling of the quantities as well as in local variations: i.e., fluctuations in the turbulent cascade rate are typically accompanied by correlated variations in the quasilinear heating rate. These novel measurements, which support a radially extended cyclotron heating mechanism, are important in understanding how turbulent dissipation results in heating the expanding solar wind. These results further indicate significant progress on PSP’s objective “to trace the flow of energy that heats and accelerates the solar corona and solar wind” (Fox et al., 2016). Furthermore, the extended ICW heating, which is observed from 15 to 55 Rsubscript𝑅direct-productR_{\odot}italic_R start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT, may have significant implications for the role of ICW heating in the corona, (Hollweg & Johnson, 1988; Cranmer, 2000).

Recent simulations by Squire et al. (2022) suggest that the dissipation of strongly imbalanced Alfvénic turbulence may result in polarization signatures similar to those measured in the solar wind (Podesta & Gary, 2011; He et al., 2011; Huang et al., 2020). The Squire et al. (2022) simulations, and the underlying idea of a helicity barrier whereby imbalanced turbulence is prevented from cascading to small scales due to a conserved generalized helicity (Meyrand et al., 2021), provides a novel framework to understand the connection between turbulence and ion-cyclotron waves. These theoretical ideas have found support in observations showing that ion-scale ICWs preferentially occur when large-scale fluctuations are highly Alfvénic (Bowen et al., 2023) and that sub-ion scale intermittency depends largely on the level of ICWs present, which is in turn correlated to the cross-helicity, suggesting that ICWs play a role in the dissipation of imbalanced turbulence.

The idea that parallel propagating ICWs dissipate a predominantly perpendicular cascade challenges our current understandings of solar wind turbulence and dissipation. It has long been suggested that the turbulence should drive fluctuations into smaller perpendicular scales (Shebalin et al., 1983), and observations regularly indicate significant anisotropy perpendicular the mean field (Horbury et al., 2008; Chen et al., 2010; Duan et al., 2021). The hybrid simulations of the helicity barrier by Squire et al. (2022) suggest that oblique ICWs are the primary heating mechanism and that parallel ICWs are emitted as a secondary process (Chandran et al., 2010) through the Alfvén/ion cyclotron instability (Gary, 1993).

While the secondary emission of parallel ICWs from oblique-ICW heating provides a mechanism to generate parallel ICWs from turbulent heating, in such a mechanism, the plasma is heated by oblique ICWs but cooled by the parallel ICWs (since they are emitted as an instability). In contrast, our measurements indicate that the parallel ICWs are robustly heating the plasma, suggesting they should be directly driven by the turbulence, as opposed to generated via oblique ICW heating. If our present measurements, consistent with previous observations (Bowen et al., 2022), are correct, this presents a conundrum, given the difficulty of sourcing quasi-parallel waves from perpendicular turbulent structures. Our use of non-parametric representations of g(𝐯)𝑔𝐯g(\mathbf{v})italic_g ( bold_v ) (Appendix B) further complicates this issue, as our measurement of ICW heating is not simply due to the incorrect parameterization of the SPANi observations with a drifting-bi-Maxwellian fit.

A possible source of systematic error is in the assumption of a cold-plasma. We assume a cold plasma dispersion (Stix, 1992), that is likely not entirely accurate for the solar wind. Implementation of warm-plasma dispersion solvers with arbitrary distribution functions (Verscharen et al., 2017; Walters et al., 2023) may enable greater understanding of parallel wave-generation. Furthermore, inclusion of beam populations (Verniero et al., 2020; Ofman et al., 2022) and α𝛼\alphaitalic_α-particles may affect instabilities (McManus et al., 2023) and may help further understanding of these dynamics.

Outside of uncertainties in the distribution function and wave-dispersion relations, a possible explanation for the in situ production of ICWs that result in solar wind heating could include the adiabatic evolution associated with expansion. Adiabatic evolution associated with expansion (Chew et al., 1956) should drive distribution functions that are unstable and emitting ICWs towards stability. The extent that expansion may enable ICW heating rather than cooling within the helicity-barrier framework, remains largely unstudied.

Furthermore, estimates of the cascade rate are notoriously difficult to estimate via single spacecraft measurements (Bandyopadhyay et al., 2018, 2020). Care must be taken in understanding these cascade rates, especially in strongly imbalanced states, in which a stationary energy flux may not exist, or a negative cascade rate may occur (Smith et al., 2009; Meyrand et al., 2021). This is especially important in the regime of imbalanced turbulence frequently observed by PSP (Meyrand et al., 2021; Squire et al., 2022). However, recently Wu et al. (2022) found good agreement between the Von Kármán decay rates with measured perpendicular proton heating rates, suggesting that the turbulent dissipation can be studied via these methods. While care has to be taken in further studies of cascade rates, the general correspondence between the measured heating rates and the von Kármán decay rates obtained in this present work and by Wu et al. (2022) is promising.

In any case, these robust observations of extended cyclotron resonant process in the inner-heliosphere suggest that kinetic scale waves are strongly coupled to the turbulent cascade and provide an important pathway to the dissipation of turbulence in the solar wind and corona.

4 Acknowledgements

T.A.B. acknowledges NASA Grant No. 80NSSC24K0272. The work of I.V. was supported by NASA grant No. 80NSSC22K1634. B.D.G.C. acknowledges the support of NASA grant 80NSSC24K0171.

References

  • Alexandrova et al. (2008) Alexandrova, O., Carbone, V., Veltri, P., & Sorriso-Valvo, L. 2008, The Astrophysical Journal, 674, 1153, doi: 10.1086/524056
  • Badman et al. (2023) Badman, S. T., Riley, P., Jones, S. I., et al. 2023, arXiv e-prints, arXiv:2303.04852, doi: 10.48550/arXiv.2303.04852
  • Bale et al. (2019) Bale, S., Badman, S., Bonnell, J., et al. 2019, Nature, 1
  • Bale et al. (2016) Bale, S. D., Goetz, K., Harvey, P. R., et al. 2016, Space Science Rev., 204, 49, doi: 10.1007/s11214-016-0244-5
  • Bandyopadhyay et al. (2018) Bandyopadhyay, R., Oughton, S., Wan, M., et al. 2018, Physical Review X, 8, 041052, doi: 10.1103/PhysRevX.8.041052
  • Bandyopadhyay et al. (2020) Bandyopadhyay, R., Goldstein, M. L., Maruca, B. A., et al. 2020, ApJs, 246, 48, doi: 10.3847/1538-4365/ab5dae
  • Bavassano et al. (1998) Bavassano, B., Pietropaolo, E., & Bruno, R. 1998, JGR, 103, 6521, doi: 10.1029/97JA03029
  • Boardsen et al. (2015) Boardsen, S. A., Jian, L. K., Raines, J. L., et al. 2015, Journal of Geophysical Research (Space Physics), 120, 10,207, doi: 10.1002/2015JA021506
  • Bowen et al. (2023) Bowen, T. A., Chandran, B. D. G., Klein, K. G., et al. 2023, in 2023 XXXVth General Assembly and Scientific Symposium of the International Union of Radio Science (URSI GASS, 335, doi: 10.23919/URSIGASS57860.2023.10265538
  • Bowen et al. (2020a) Bowen, T. A., Mallet, A., Huang, J., et al. 2020a, ApJS, 246, 66, doi: 10.3847/1538-4365/ab6c65
  • Bowen et al. (2020b) Bowen, T. A., Bale, S. D., Bonnell, J. W., et al. 2020b, ApJ, 899, 74, doi: 10.3847/1538-4357/ab9f37
  • Bowen et al. (2020c) Bowen, T. A., Mallet, A., Bale, S. D., et al. 2020c, PRL, 125, 025102, doi: 10.1103/PhysRevLett.125.025102
  • Bowen et al. (2020d) Bowen, T. A., Bale, S. D., Bonnell, J. W., et al. 2020d, Journal of Geophysical Research (Space Physics), 125, e27813, doi: 10.1029/2020JA027813
  • Bowen et al. (2022) Bowen, T. A., Chandran, B. D. G., Squire, J., et al. 2022, PRL, 129, 165101, doi: 10.1103/PhysRevLett.129.165101
  • Broomhead & Lowe (1988) Broomhead, D., & Lowe, D. 1988, Complex Systems, 2, 321
  • Bruno et al. (2014) Bruno, R., Trenchi, L., & Telloni, D. 2014, ApJL, 793, L15, doi: 10.1088/2041-8205/793/1/L15
  • Chandran et al. (2010) Chandran, B. D. G., Pongkitiwanichakul, P., Isenberg, P. A., et al. 2010, ApJ, 722, 710, doi: 10.1088/0004-637X/722/1/710
  • Chen et al. (2010) Chen, C. H. K., Horbury, T. S., Schekochihin, A. A., et al. 2010, Physical Review Letters, 104, 255002, doi: 10.1103/PhysRevLett.104.255002
  • Chew et al. (1956) Chew, G. F., Goldberger, M. L., & Low, F. E. 1956, Proceedings of the Royal Society of London Series A, 236, 112, doi: 10.1098/rspa.1956.0116
  • Cranmer (2000) Cranmer, S. R. 2000, The Astrophysical Journal, 532, 1197, doi: 10.1086/308620
  • Cranmer (2014) —. 2014, ApJS, 213, 16, doi: 10.1088/0067-0049/213/1/16
  • Davis et al. (2023) Davis, N., Chandran, B. D. G., Bowen, T. A., et al. 2023, arXiv e-prints, arXiv:2303.01663, doi: 10.48550/arXiv.2303.01663
  • Denskat et al. (1983) Denskat, K. U., Beinroth, H. J., & Neubauer, F. M. 1983, Journal of Geophysics Zeitschrift Geophysik, 54, 60
  • Duan et al. (2021) Duan, D., He, J., Bowen, T. A., et al. 2021, ApJL, 915, L8, doi: 10.3847/2041-8213/ac07ac
  • Duan et al. (2020) Duan, D., Bowen, T. A., Chen, C. H. K., et al. 2020, ApJS, 246, 55, doi: 10.3847/1538-4365/ab672d
  • Dudok de Wit et al. (2013) Dudok de Wit, T., Alexandrova, O., Furno, I., Sorriso-Valvo, L., & Zimbardo, G. 2013, Space Sci. Rev., 178, 665, doi: 10.1007/s11214-013-9974-9
  • Dudok de Wit et al. (2022) Dudok de Wit, T., Krasnoselskikh, V. V., Agapitov, O., et al. 2022, Journal of Geophysical Research (Space Physics), 127, e30018, doi: 10.1029/2021JA030018
  • Dum et al. (1980) Dum, C. T., Marsch, E., & Pilipp, W. 1980, Journal of Plasma Physics, 23, 91, doi: 10.1017/S0022377800022170
  • Farge (1992) Farge, M. 1992, Annual Review of Fluid Mechanics, 24, 395, doi: 10.1146/annurev.fl.24.010192.002143
  • Fox et al. (2016) Fox, N. J., Velli, M. C., Bale, S. D., et al. 2016, Space Science Reviews, 204, 7, doi: 10.1007/s11214-015-0211-6
  • Gary (1993) Gary, S. P. 1993, Theory of Space Plasma Microinstabilities
  • Goldstein et al. (1994) Goldstein, M. L., Roberts, D. A., & Fitch, C. A. 1994, Journal of Geophysical Research, 99, 11519, doi: 10.1029/94JA00789
  • He et al. (2011) He, J., Marsch, E., Tu, C., Yao, S., & Tian, H. 2011, The Astrophysical Journal, 731, 85, doi: 10.1088/0004-637X/731/2/85
  • He et al. (2015) He, J., Wang, L., Tu, C., Marsch, E., & Zong, Q. 2015, ApJ, 800, L31, doi: 10.1088/2041-8205/800/2/L31
  • Hellinger et al. (2013) Hellinger, P., TráVníček, P. M., Štverák, Š., Matteini, L., & Velli, M. 2013, Journal of Geophysical Research (Space Physics), 118, 1351, doi: 10.1002/jgra.50107
  • Heuer & Marsch (2007) Heuer, M., & Marsch, E. 2007, Journal of Geophysical Research (Space Physics), 112, A03102, doi: 10.1029/2006JA011979
  • Hollweg & Isenberg (2002) Hollweg, J. V., & Isenberg, P. A. 2002, Journal of Geophysical Research (Space Physics), 107, 1147, doi: 10.1029/2001JA000270
  • Hollweg & Johnson (1988) Hollweg, J. V., & Johnson, W. 1988, JGR, 93, 9547, doi: 10.1029/JA093iA09p09547
  • Hollweg & Markovskii (2002) Hollweg, J. V., & Markovskii, S. A. 2002, Journal of Geophysical Research (Space Physics), 107, 1080, doi: 10.1029/2001JA000205
  • Horbury et al. (2008) Horbury, T. S., Forman, M., & Oughton, S. 2008, PRL, 101, 175005, doi: 10.1103/PhysRevLett.101.175005
  • Hossain et al. (1995) Hossain, M., Gray, P. C., Pontius, Duane H., J., Matthaeus, W. H., & Oughton, S. 1995, Physics of Fluids, 7, 2886, doi: 10.1063/1.868665
  • Howes & Quataert (2010) Howes, G. G., & Quataert, E. 2010, The Astrophysical Journal Letters, 709, L49, doi: 10.1088/2041-8205/709/1/L49
  • Huang et al. (2020) Huang, S. Y., Zhang, J., Sahraoui, F., et al. 2020, ApJL, 897, L3, doi: 10.3847/2041-8213/ab9abb
  • Isenberg & Lee (1996) Isenberg, P. A., & Lee, M. A. 1996, JGR, 101, 11055, doi: 10.1029/96JA00293
  • Jannet et al. (2021) Jannet, G., Dudok de Wit, T., Krasnoselskikh, V., et al. 2021, Journal of Geophysical Research (Space Physics), 126, e28543, doi: 10.1029/2020JA028543
  • Jian et al. (2014) Jian, L. K., Wei, H. Y., Russell, C. T., et al. 2014, ApJ, 786, 123, doi: 10.1088/0004-637X/786/2/123
  • Karpman (1974) Karpman, V. I. 1974, SSR, 16, 361, doi: 10.1007/BF00171564
  • Kennel & Engelmann (1966) Kennel, C. F., & Engelmann, F. 1966, Physics of Fluids, 9, 2377, doi: 10.1063/1.1761629
  • Kiyani et al. (2009) Kiyani, K. H., Chapman, S. C., Khotyaintsev, Y. V., Dunlop, M. W., & Sahraoui, F. 2009, Physical Review Letters, 103, 075006, doi: 10.1103/PhysRevLett.103.075006
  • Klein et al. (2018) Klein, K. G., Alterman, B. L., Stevens, M. L., Vech, D., & Kasper, J. C. 2018, PRL, 120, 205102, doi: 10.1103/PhysRevLett.120.205102
  • Klein et al. (2021) Klein, K. G., Verniero, J. L., Alterman, B., et al. 2021, ApJ, 909, 7, doi: 10.3847/1538-4357/abd7a0
  • Leamon et al. (1998) Leamon, R. J., Smith, C. W., Ness, N. F., Matthaeus, W. H., & Wong, H. K. 1998, Journal of Geophysical Research, 103, 4775, doi: 10.1029/97JA03394
  • Lion et al. (2016) Lion, S., Alexandrova, O., & Zaslavsky, A. 2016, The Astrophysical Journal, 824, 47, doi: 10.3847/0004-637X/824/1/47
  • Liu et al. (2023) Liu, W., Zhao, J., Wang, T., et al. 2023, ApJ, 951, 69, doi: 10.3847/1538-4357/acd53b
  • Livi et al. (2022) Livi, R., Larson, D. E., Kasper, J. C., et al. 2022, ApJ, 938, 138, doi: 10.3847/1538-4357/ac93f5
  • Lotz et al. (2023) Lotz, S., Nel, A. E., Wicks, R. T., et al. 2023, ApJ, 942, 93, doi: 10.3847/1538-4357/aca903
  • Luo et al. (2022) Luo, Q., Zhu, X., He, J., et al. 2022, ApJ, 928, 36, doi: 10.3847/1538-4357/ac52a9
  • Marsch (2006) Marsch, E. 2006, Living Reviews in Solar Physics, 3, 1, doi: 10.12942/lrsp-2006-1
  • Marsch & Goldstein (1983) Marsch, E., & Goldstein, H. 1983, J. Geophys. Res., 88, 9933, doi: 10.1029/JA088iA12p09933
  • Marsch et al. (1982) Marsch, E., Schwenn, R., Rosenbauer, H., et al. 1982, J. Geophys. Res., 87, 52, doi: 10.1029/JA087iA01p00052
  • Marsch & Tu (2001a) Marsch, E., & Tu, C. Y. 2001a, JGR, 106, 8357, doi: 10.1029/2000JA000414
  • Marsch & Tu (2001b) —. 2001b, JGR, 106, 227, doi: 10.1029/2000JA000042
  • McManus et al. (2020) McManus, M. D., Bowen, T. A., Mallet, A., et al. 2020, ApJS, 246, 67, doi: 10.3847/1538-4365/ab6dce
  • McManus et al. (2023) McManus, M. D., Klein, K. G., Larson, D., et al. 2023, arXiv e-prints, arXiv:2310.14136, doi: 10.48550/arXiv.2310.14136
  • Meyrand et al. (2021) Meyrand, R., Squire, J., Schekochihin, A. A., & Dorland, W. 2021, Journal of Plasma Physics, 87, 535870301, doi: 10.1017/S0022377821000489
  • Ofman et al. (2022) Ofman, L., Boardsen, S. A., Jian, L. K., Verniero, J. L., & Larson, D. 2022, ApJ, 926, 185, doi: 10.3847/1538-4357/ac402c
  • Parker (1958) Parker, E. N. 1958, ApJ, 128, 664, doi: 10.1086/146579
  • Podesta & Gary (2011) Podesta, J. J., & Gary, S. P. 2011, The Astrophysical Journal, 734, 15, doi: 10.1088/0004-637X/734/1/15
  • Richardson et al. (1995) Richardson, J. D., Paularena, K. I., Lazarus, A. J., & Belcher, J. W. 1995, Geophys. Research Letters, 22, 325, doi: 10.1029/94GL03273
  • Roberts et al. (1987) Roberts, D. A., Klein, L. W., Goldstein, M. L., & Matthaeus, W. H. 1987, JGR, 92, 11021, doi: 10.1029/JA092iA10p11021
  • Sahraoui et al. (2009) Sahraoui, F., Goldstein, M. L., Robert, P., & Khotyaintsev, Y. V. 2009, Physical Review Letters, 102, 231102, doi: 10.1103/PhysRevLett.102.231102
  • Shankarappa et al. (2023) Shankarappa, N., Klein, K. G., & Martinović, M. M. 2023, ApJ, 946, 85, doi: 10.3847/1538-4357/acb542
  • Shebalin et al. (1983) Shebalin, J. V., Matthaeus, W. H., & Montgomery, D. 1983, Journal of Plasma Physics, 29, 525, doi: 10.1017/S0022377800000933
  • Shklyar & Matsumoto (2009) Shklyar, D., & Matsumoto, H. 2009, Surveys in Geophysics, 30, 55, doi: 10.1007/s10712-009-9061-7
  • Smith et al. (2006) Smith, C. W., Hamilton, K., Vasquez, B. J., & Leamon, R. J. 2006, The Astrophysical Journal Letters, 645, L85, doi: 10.1086/506151
  • Smith et al. (2009) Smith, C. W., Stawarz, J. E., Vasquez, B. J., Forman, M. A., & MacBride, B. T. 2009, PRL, 103, 201101, doi: 10.1103/PhysRevLett.103.201101
  • Smith et al. (2012) Smith, C. W., Vasquez, B. J., & Hollweg, J. V. 2012, ApJ, 745, 8, doi: 10.1088/0004-637X/745/1/8
  • Smith et al. (1990) Smith, W. C., Matthaeus, H. W., & Ness, F. N. 1990, in International Cosmic Ray Conference, Vol. 5, International Cosmic Ray Conference, 280
  • Squire et al. (2022) Squire, J., Meyrand, R., Kunz, M. W., et al. 2022, Nature Astronomy, 6, 715, doi: 10.1038/s41550-022-01624-z
  • Stix (1992) Stix, T. H. 1992, Waves in plasmas
  • Telloni et al. (2019) Telloni, D., Carbone, F., Bruno, R., et al. 2019, ApJL, 885, L5, doi: 10.3847/2041-8213/ab4c44
  • Tu & Marsch (1995) Tu, C. Y., & Marsch, E. 1995, Space Sci Rev, 73, 1, doi: 10.1007/BF00748891
  • Tu & Marsch (1997) —. 1997, Sol Phys, 171, 363, doi: 10.1023/A:1004968327196
  • Tu & Marsch (2001) —. 2001, JGR, 106, 8233, doi: 10.1029/2000JA000024
  • Usmanov et al. (2014) Usmanov, A. V., Goldstein, M. L., & Matthaeus, W. H. 2014, ApJ, 788, 43, doi: 10.1088/0004-637X/788/1/43
  • Vech et al. (2018) Vech, D., Mallet, A., Klein, K. G., & Kasper, J. C. 2018, The Astrophysical Journal Letters, 855, L27, doi: 10.3847/2041-8213/aab351
  • Vech et al. (2020) Vech, D., Kasper, J. C., Klein, K. G., et al. 2020, ApJS, 246, 52, doi: 10.3847/1538-4365/ab60a2
  • Verniero et al. (2020) Verniero, J. L., Larson, D. E., Livi, R., et al. 2020, ApJS, 248, 5, doi: 10.3847/1538-4365/ab86af
  • Verniero et al. (2022) Verniero, J. L., Chandran, B. D. G., Larson, D. E., et al. 2022, ApJ, 924, 112, doi: 10.3847/1538-4357/ac36d5
  • Verscharen et al. (2017) Verscharen, D., Chen, C. H. K., & Wicks, R. T. 2017, ApJ, 840, 106, doi: 10.3847/1538-4357/aa6a56
  • Viñas & Gurgiolo (2009) Viñas, A. F., & Gurgiolo, C. 2009, Journal of Geophysical Research (Space Physics), 114, A01105, doi: 10.1029/2008JA013633
  • Walters et al. (2023) Walters, J., Klein, K. G., Lichko, E., et al. 2023, ApJ, 955, 97, doi: 10.3847/1538-4357/acf1fa
  • Wan et al. (2012) Wan, M., Matthaeus, W. H., Karimabadi, H., et al. 2012, PRL, 109, 195001, doi: 10.1103/PhysRevLett.109.195001
  • Wicks et al. (2016) Wicks, R. T., Alexander, R. L., Stevens, M., et al. 2016, ApJ, 819, 6, doi: 10.3847/0004-637X/819/1/6
  • Woodham et al. (2018) Woodham, L. D., Wicks, R. T., Verscharen, D., & Owen, C. J. 2018, The Astrophysical Journal, 856, 49, doi: 10.3847/1538-4357/aab03d
  • Woodham et al. (2019) Woodham, L. D., Wicks, R. T., Verscharen, D., et al. 2019, ApJL, 884, L53, doi: 10.3847/2041-8213/ab4adc
  • Wu et al. (2022) Wu, H., Tu, C., He, J., Wang, X., & Yang, L. 2022, ApJ, 926, 116, doi: 10.3847/1538-4357/ac4413
  • Wüest et al. (2007) Wüest, M., Evans, D. S., & von Steiger, R. 2007, Calibration of Particle Instruments in Space Physics
  • Zhao et al. (2021) Zhao, G. Q., Lin, Y., Wang, X. Y., et al. 2021, ApJ, 906, 123, doi: 10.3847/1538-4357/abca3b

Appendix A Numerical Determination of Quasilinear Heating Rate

Determination of QICWsubscript𝑄ICWQ_{\rm{ICW}}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT and Q(v,v)𝑄subscript𝑣perpendicular-tosubscript𝑣parallel-toQ(v_{\perp},v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) is obtained through integrating Eq 2 by parts to obtain

QICW=π2Ωp2mp𝑑v0𝑑v[0𝑑kv2δ(ωkkvΩp)ωk2k2I(k)G^kg¯(v,v)]subscript𝑄ICWsuperscript𝜋2superscriptsubscriptΩ𝑝2subscript𝑚𝑝superscriptsubscriptdifferential-dsubscript𝑣parallel-tosuperscriptsubscript0differential-dsubscript𝑣perpendicular-todelimited-[]superscriptsubscript0differential-dsubscript𝑘parallel-tosuperscriptsubscript𝑣perpendicular-to2𝛿subscript𝜔𝑘subscript𝑘parallel-tosubscript𝑣parallel-tosubscriptΩ𝑝superscriptsubscript𝜔𝑘2superscriptsubscript𝑘parallel-to2𝐼subscript𝑘parallel-tosubscript^𝐺𝑘¯𝑔subscript𝑣perpendicular-tosubscript𝑣parallel-to\displaystyle Q_{\rm ICW}=-\pi^{2}\Omega_{p}^{2}m_{p}\int_{-\infty}^{\infty}dv% _{\parallel}\int_{0}^{\infty}dv_{\perp}\Bigg{[}\int_{0}^{\infty}dk_{\parallel}% v_{\perp}^{2}\delta(\omega_{k}-k_{\parallel}v_{\parallel}-\Omega_{p})\frac{% \omega_{k}^{2}}{k_{\parallel}^{2}}I(k_{\parallel})\hat{G}_{k}\bar{g}(v_{\perp}% ,v_{\parallel})\Bigg{]}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT = - italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Ω start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT [ ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_δ ( italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT - italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT - roman_Ω start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) divide start_ARG italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_I ( italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) over^ start_ARG italic_G end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT over¯ start_ARG italic_g end_ARG ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) ] (A1)

The δlimit-from𝛿\delta-italic_δ -function in Eq 2& A1 corresponds to the ICW resonance condition, and defines a single parallel wave number that resonates with particles at each parallel velocity. Change of variables yields a δ𝛿\deltaitalic_δ-function for resonant wave number kres(v)subscript𝑘𝑟𝑒𝑠subscript𝑣parallel-tok_{res}(v_{\parallel})italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) as

δ(ωkkvΩp)=δ(kkres)|vg(kres)v|,𝛿subscript𝜔𝑘subscript𝑘parallel-tosubscript𝑣parallel-tosubscriptΩ𝑝𝛿subscript𝑘parallel-tosubscript𝑘𝑟𝑒𝑠subscript𝑣𝑔subscript𝑘𝑟𝑒𝑠subscript𝑣parallel-to\delta(\omega_{k}-k_{\parallel}v_{\parallel}-\Omega_{p})=\frac{\delta(k_{% \parallel}-k_{res})}{|v_{g}(k_{res})-v_{\parallel}|},italic_δ ( italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT - italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT - roman_Ω start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) = divide start_ARG italic_δ ( italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT - italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ) end_ARG start_ARG | italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ) - italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT | end_ARG , (A2)

where vg(kres)=ωk/k|kressubscript𝑣𝑔subscript𝑘𝑟𝑒𝑠evaluated-atsubscript𝜔𝑘𝑘subscript𝑘𝑟𝑒𝑠v_{g}(k_{res})=\partial\omega_{k}/\partial k|_{k_{res}}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ) = ∂ italic_ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT / ∂ italic_k | start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT is the group speed evaluated at each kressubscript𝑘𝑟𝑒𝑠k_{res}italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT. For positive ksubscript𝑘parallel-tok_{\parallel}italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT, taken as outward going ICWs, there is no resonance unless v<0subscript𝑣parallel-to0v_{\parallel}<0italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT < 0, limiting our concern only to negative thermal speeds. We use the δlimit-from𝛿\delta-italic_δ -function to substitute kres(v)subscript𝑘𝑟𝑒𝑠subscript𝑣parallel-tok_{res}(v_{\parallel})italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) for ksubscript𝑘parallel-tok_{\parallel}italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT, thereby evaluating the integral over wavenumber, Eq A1 becomes

QICW=π2Ωp2mp0𝑑v0𝑑v[v2|vg(kres)v|ωkres2kres2I(kres)G^kg¯(v,v)],subscript𝑄ICWsuperscript𝜋2superscriptsubscriptΩ𝑝2subscript𝑚𝑝superscriptsubscript0differential-dsubscript𝑣parallel-tosuperscriptsubscript0differential-dsubscript𝑣perpendicular-todelimited-[]superscriptsubscript𝑣perpendicular-to2subscript𝑣𝑔subscript𝑘𝑟𝑒𝑠subscript𝑣parallel-tosuperscriptsubscript𝜔subscript𝑘𝑟𝑒𝑠2superscriptsubscript𝑘𝑟𝑒𝑠2𝐼subscript𝑘𝑟𝑒𝑠subscript^𝐺𝑘¯𝑔subscript𝑣perpendicular-tosubscript𝑣parallel-to\displaystyle Q_{\rm ICW}=-\pi^{2}\Omega_{p}^{2}m_{p}\int_{-\infty}^{0}dv_{% \parallel}\int_{0}^{\infty}dv_{\perp}\Bigg{[}\frac{v_{\perp}^{2}}{|v_{g}(k_{% res})-v_{\parallel}|}\frac{\omega_{k_{res}}^{2}}{k_{res}^{2}}I(k_{res})\hat{G}% _{k}\bar{g}(v_{\perp},v_{\parallel})\Bigg{]},italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT = - italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Ω start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT italic_d italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_d italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT [ divide start_ARG italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG | italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ) - italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT | end_ARG divide start_ARG italic_ω start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_I ( italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ) over^ start_ARG italic_G end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT over¯ start_ARG italic_g end_ARG ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) ] , (A3)

where the spectrum I(k)𝐼subscript𝑘parallel-toI(k_{\parallel})italic_I ( italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) is interpolated from ksubscript𝑘parallel-tok_{\parallel}italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT to kressubscript𝑘𝑟𝑒𝑠k_{res}italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT, constructing I(kres)𝐼subscript𝑘𝑟𝑒𝑠I(k_{res})italic_I ( italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ). Eq A3 is entirely a function of vsubscript𝑣parallel-tov_{\parallel}italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT and vsubscript𝑣perpendicular-tov_{\perp}italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT and can be calculated via a numerical integration over velocity space. We further identify the differential heating rate

Q(v,v)=π2Ωp2mpv2|vg(kres)v|ωkres2kres2I(kres)G^kg¯(v,v).𝑄subscript𝑣perpendicular-tosubscript𝑣parallel-tosuperscript𝜋2superscriptsubscriptΩ𝑝2subscript𝑚𝑝superscriptsubscript𝑣perpendicular-to2subscript𝑣𝑔subscript𝑘𝑟𝑒𝑠subscript𝑣parallel-tosuperscriptsubscript𝜔subscript𝑘𝑟𝑒𝑠2superscriptsubscript𝑘𝑟𝑒𝑠2𝐼subscript𝑘𝑟𝑒𝑠subscript^𝐺𝑘¯𝑔subscript𝑣perpendicular-tosubscript𝑣parallel-to\displaystyle Q(v_{\perp},v_{\parallel})=-\pi^{2}\Omega_{p}^{2}m_{p}\frac{v_{% \perp}^{2}}{|v_{g}(k_{res})-v_{\parallel}|}\frac{\omega_{k_{res}}^{2}}{k_{res}% ^{2}}I(k_{res})\hat{G}_{k}\bar{g}(v_{\perp},v_{\parallel}).italic_Q ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) = - italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_Ω start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_m start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT divide start_ARG italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG | italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ( italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ) - italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT | end_ARG divide start_ARG italic_ω start_POSTSUBSCRIPT italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_I ( italic_k start_POSTSUBSCRIPT italic_r italic_e italic_s end_POSTSUBSCRIPT ) over^ start_ARG italic_G end_ARG start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT over¯ start_ARG italic_g end_ARG ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) . (A4)

To perform the numerical integral in Eq A3, we evaluate Q(v,v)𝑄subscript𝑣perpendicular-tosubscript𝑣parallel-toQ(v_{\perp},v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) on a 1 km/s ×\times× 1 km/s grid, which was found to be sufficient for convergence of Eq A3. The range of integration is set to -1000 km/s <v<absentsubscript𝑣parallel-toabsent<v_{\parallel}<< italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT < 0 km/s and 0 km/s <v<absentsubscript𝑣perpendicular-toabsent<v_{\perp}<< italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT < 1000 km/s. We have removed contributions to the spectra with kdi<0.3subscript𝑘parallel-tosubscript𝑑𝑖0.3k_{\parallel}d_{i}<0.3italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT < 0.3 and kdi>5subscript𝑘parallel-tosubscript𝑑𝑖5k_{\parallel}d_{i}>5italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT > 5, which for βp=0.1subscript𝛽𝑝0.1\beta_{p}=0.1italic_β start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = 0.1 corresponds to an approximate range of thermal speeds between 0.02 and 7.5 vthsubscriptsubscript𝑣𝑡parallel-to{v_{th}}_{\parallel}italic_v start_POSTSUBSCRIPT italic_t italic_h end_POSTSUBSCRIPT start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT. Figure 5(a&b) shows that the parallel thermal speeds hovers near 50 km/s over this entire interval. Thus, the 1000km/s limits are sufficiently large with respect to the thermal and Alfvén speeds for this stream that there are no significant contributions in portion of phase space out of the bounds of integration. The heating rate as a function of parallel resonant velocity Q(v)𝑄subscript𝑣parallel-toQ(v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) is computed by integrating Q(v,v)𝑄subscript𝑣perpendicular-tosubscript𝑣parallel-toQ(v_{\perp},v_{\parallel})italic_Q ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) over perpendicular velocities.

These methods have previously been used to obtain heating rates QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT of order 10%similar-toabsentpercent10\sim 10\%∼ 10 % of the turbulent cascade rate (Bowen et al., 2022) though this previously studied interval contained significantly less amounts of left-hand polarized waves than this stream. Our results shown in Fig 5 suggest that a great amount of the turbulent energy likely enters the plasma via cyclotron resonant heating.

Appendix B Application of Non-Parametric Representations

While in situ observations of the proton velocity distribution have long been recognized to have a mostly Gaussian shape (Marsch & Goldstein, 1983) and drifting biMaxwellian approximation has regularly been used to approximate ion velocity distributions (Marsch et al., 1982; Marsch & Tu, 2001b; Tu & Marsch, 2001; Heuer & Marsch, 2007; Klein et al., 2021; Verniero et al., 2022; McManus et al., 2023); however, it is clear that deviations from nonthermal structure in the velocity distribution can affect wave-particle resonant processes (Dum et al., 1980; Bowen et al., 2022; Walters et al., 2023). Furthermore, it is well known that the equilibrium kinetic contours for dispersive ICW resonance do not coincide with Maxwellian curvature (Isenberg & Lee, 1996), such that relaxation via ICW resonant processing cannot be entirely captured via bi-Maxwellian approximations (Heuer & Marsch, 2007). In this sense it is important to test the heating rates computed from the drifting biMaxwellian model in Equation 11 against heating rates computed from non-parametric models of the velocity distribution.

B.1 Hermite Polynomial Interpolation

Bowen et al. (2022) previously used linear least-square fits of orthogonal Hermite-polynomials to estimate non-Maxwellian features in the velocity distribution; following this previously developed method, we perform a linear least square fit to the average proton velocity distribution in each 128s interval using the Hermite polynomials Hnsubscript𝐻𝑛H_{n}italic_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT, and Hermite functions ϕmsubscriptitalic-ϕ𝑚\phi_{m}italic_ϕ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT

g¯H(v,v)=m,nAmnϕm(v/wth)ϕn(v/wth)\displaystyle\bar{g}^{H}(v_{\perp},v_{\parallel})=\sum_{m,n}A_{mn}\phi^{m}(v_{% \perp}/w_{\perp{th}})\phi^{n}(v_{\parallel}/w_{\parallel{th}})over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_m italic_n end_POSTSUBSCRIPT italic_ϕ start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT / italic_w start_POSTSUBSCRIPT ⟂ italic_t italic_h end_POSTSUBSCRIPT ) italic_ϕ start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT / italic_w start_POSTSUBSCRIPT ∥ italic_t italic_h end_POSTSUBSCRIPT ) (B1)
Hn(v)=(1)nev2dndxnev2subscript𝐻𝑛𝑣superscript1𝑛superscript𝑒superscript𝑣2superscript𝑑𝑛𝑑superscript𝑥𝑛superscript𝑒superscript𝑣2\displaystyle H_{n}(v)=(-1)^{n}e^{v^{2}}\frac{d^{n}}{dx^{n}}e^{-v^{2}}italic_H start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_v ) = ( - 1 ) start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT divide start_ARG italic_d start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT end_ARG start_ARG italic_d italic_x start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT end_ARG italic_e start_POSTSUPERSCRIPT - italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT (B2)
ϕm=Hm(v)2mπ1/2m!ev2.superscriptitalic-ϕ𝑚superscript𝐻𝑚𝑣superscript2𝑚superscript𝜋12𝑚superscript𝑒superscript𝑣2\displaystyle\phi^{m}=\frac{H^{m}(v)}{\sqrt{2^{m}\pi^{1/2}m!}}e^{-v^{2}}.italic_ϕ start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT = divide start_ARG italic_H start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ( italic_v ) end_ARG start_ARG square-root start_ARG 2 start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT italic_π start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT italic_m ! end_ARG end_ARG italic_e start_POSTSUPERSCRIPT - italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT . (B3)

We fit an orthogonal set of Hermite functions to each velocity distribution through constructing the matrix 𝐀𝐀\mathbf{A}bold_A, which consists of the transform coefficients Amnsubscript𝐴𝑚𝑛A_{mn}italic_A start_POSTSUBSCRIPT italic_m italic_n end_POSTSUBSCRIPT corresponding to each two-dimensional Hermite polynomial combination. Through linear least square fitting of ΦΦ\Phiroman_Φ, with elements ϕmϕnsuperscriptitalic-ϕ𝑚superscriptitalic-ϕ𝑛\phi^{m}\phi^{n}italic_ϕ start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT italic_ϕ start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT, to the observed velocity distribution, written as matrix 𝐠𝐠\mathbf{g}bold_g such that 𝐠=𝐀𝚽𝐠𝐀𝚽\mathbf{g}={\mathbf{A}}{\mathbf{\Phi}}bold_g = bold_A bold_Φ and least square fit inversion gives:

𝐀=𝐠𝚽T(𝚽𝚽T)1,𝐀𝐠superscript𝚽𝑇superscript𝚽superscript𝚽𝑇1\mathbf{A}=\mathbf{g}\mathbf{\Phi}^{T}({\mathbf{\Phi}}{\mathbf{\Phi}}^{T})^{-1},bold_A = bold_g bold_Φ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( bold_Φ bold_Φ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT , (B4)

where a singular value decomposition determines the pseudo-inverse of 𝚽𝚽T𝚽superscript𝚽𝑇{\mathbf{\Phi}}{\mathbf{\Phi}}^{T}bold_Φ bold_Φ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT. The matrix 𝐀𝐀\mathbf{A}bold_A then gives the best-fit coefficients for the Hermite polynomial composition.

We furthermore adopt errors on 𝐠𝐠\mathbf{g}bold_g associated with Poisson counting, such that the uncertainty at each velocity coordinate k𝑘kitalic_k is σgk=gksubscript𝜎subscript𝑔𝑘subscript𝑔𝑘\sigma_{g_{k}}=\sqrt{g_{k}}italic_σ start_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT = square-root start_ARG italic_g start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG (Wüest et al., 2007). Poisson noise between energy bins is uncorrelated such that we construct the diagonal-weight matrix as 𝐖𝐖\mathbf{W}bold_W with entries σ2superscript𝜎2\sigma^{2}italic_σ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT can be included in weighted linear least square fitting as 𝐠𝐖=𝐀𝚽𝐖𝐠𝐖𝐀𝚽𝐖\mathbf{gW}={\mathbf{A}}{\mathbf{\Phi}}\mathbf{W}bold_gW = bold_A bold_Φ bold_W such that

𝐀=𝐠𝐖𝚽T(𝚽𝐖𝚽T)1𝐀𝐠𝐖superscript𝚽𝑇superscript𝚽𝐖superscript𝚽𝑇1\mathbf{A=gW\Phi}^{T}\mathbf{(\Phi W}\mathbf{\Phi}^{T})^{-1}bold_A = bold_gW bold_Φ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( bold_Φ bold_W bold_Φ start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT (B5)

The coefficients 𝐀𝐀\mathbf{A}bold_A are then used to construct g¯H(v,v)superscript¯𝑔𝐻subscript𝑣perpendicular-tosubscript𝑣parallel-to\bar{g}^{H}(v_{\perp},v_{\parallel})over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) using Equation B1.

Using g¯H(v,v)superscript¯𝑔𝐻subscript𝑣perpendicular-tosubscript𝑣parallel-to\bar{g}^{H}(v_{\perp},v_{\parallel})over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ), which is an analytic, differentiable funcion, it is possible to separately estimate QICWsubscript𝑄ICWQ_{\rm ICW}italic_Q start_POSTSUBSCRIPT roman_ICW end_POSTSUBSCRIPT without the assumption of a drifting bi-Maxwellian structure to the velocity distribution.

B.2 Radial Basis Function Interpolation

We further include another approximation to the velocity distributions through interpolation via radial basis functions (RBFs) (Broomhead & Lowe, 1988; Bowen et al., 2023) to model the distribution function via a summed set of interpolating functions given by

gpRBF(v,v)=iNRBF1Riψi(ζi),superscriptsubscript𝑔𝑝𝑅𝐵𝐹subscript𝑣perpendicular-tosubscript𝑣parallel-tosuperscriptsubscript𝑖subscript𝑁𝑅𝐵𝐹1subscript𝑅𝑖subscript𝜓𝑖subscript𝜁𝑖\displaystyle{g_{p}}^{RBF}(v_{\perp},v_{\parallel})={\sum_{i}^{N_{RBF}-1}}R_{i% }\psi_{i}(\zeta_{i}),italic_g start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R italic_B italic_F end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT - 1 end_POSTSUPERSCRIPT italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_ψ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_ζ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , (B6)

with ζi=𝐯𝐯icsubscript𝜁𝑖𝐯subscriptsuperscript𝐯𝑐𝑖\zeta_{i}=\mathbf{v}-\mathbf{v}^{c}_{i}italic_ζ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = bold_v - bold_v start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, where 𝐯icsubscriptsuperscript𝐯𝑐𝑖\mathbf{v}^{c}_{i}bold_v start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is the center of each interpolating RBF, denoted by ψisubscript𝜓𝑖\psi_{i}italic_ψ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT. The RBF method requires choice of a basis function, ψ𝜓\psiitalic_ψ, which we choose to be an isotropic 2D bi-Maxwellian

ψi(v,v)=1π3wRBFi3exp[(vvic)2wRBFi2(vvic)2wRBFi2].\displaystyle\psi_{i}(v_{\perp},v_{\parallel})=\frac{1}{\sqrt{\pi^{3}}w_{RBFi}% ^{3}}\text{exp}\left[-\frac{(v_{\perp}-v^{c}_{\perp i})^{2}}{w_{RBFi}^{2}}-% \frac{(v_{\parallel}-v^{c}_{\parallel i})^{2}}{w_{RBFi}^{2}}\right].italic_ψ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_π start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG italic_w start_POSTSUBSCRIPT italic_R italic_B italic_F italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG exp [ - divide start_ARG ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT - italic_v start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ⟂ italic_i end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUBSCRIPT italic_R italic_B italic_F italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - divide start_ARG ( italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT - italic_v start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT start_POSTSUBSCRIPT ∥ italic_i end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_w start_POSTSUBSCRIPT italic_R italic_B italic_F italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] . (B7)

The number of interpolating functions, NRBFsubscript𝑁𝑅𝐵𝐹N_{RBF}italic_N start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT, must be specified along with the thermal speed wRBFisubscript𝑤𝑅𝐵𝐹𝑖w_{{RBFi}}italic_w start_POSTSUBSCRIPT italic_R italic_B italic_F italic_i end_POSTSUBSCRIPT and central location, 𝐯icsubscriptsuperscript𝐯𝑐𝑖\mathbf{v}^{c}_{i}bold_v start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, of each of the ψisubscript𝜓𝑖\psi_{i}italic_ψ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT implemented in the interpolation. In principle, wRBFisubscript𝑤𝑅𝐵𝐹𝑖w_{RBFi}italic_w start_POSTSUBSCRIPT italic_R italic_B italic_F italic_i end_POSTSUBSCRIPT does not need to be isotropic and can vary for each ψisubscript𝜓𝑖\psi_{i}italic_ψ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT; however, we determine that setting a Maxwellian-RBF at each measured SPAN energy bin with nonzero phase space density gives suitable results. Furthermore, we uniformly set wRBF=40subscript𝑤𝑅𝐵𝐹40w_{RBF}=40italic_w start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT = 40 km/s, which we find is suitable for interpolating the proton distribution.

The interpolated velocity distribution via radial basis functions is then given by

g¯RBF(v,v)=iNRBF1Riψi(ζ),superscript¯𝑔𝑅𝐵𝐹subscript𝑣perpendicular-tosubscript𝑣parallel-tosuperscriptsubscript𝑖subscript𝑁𝑅𝐵𝐹1subscript𝑅𝑖subscript𝜓𝑖𝜁\displaystyle\bar{g}^{RBF}(v_{\perp},v_{\parallel})={\sum_{i}^{N_{RBF}-1}}R_{i% }\psi_{i}(\zeta),over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_R italic_B italic_F end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT - 1 end_POSTSUPERSCRIPT italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_ψ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_ζ ) , (B8)

with ζ=vvc𝜁𝑣subscript𝑣𝑐\zeta=\vec{v}-\vec{v}_{c}italic_ζ = over→ start_ARG italic_v end_ARG - over→ start_ARG italic_v end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT. Determination of the weights Risubscript𝑅𝑖R_{i}italic_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is again performed through SVD estimation for the pseudoinverse giving a least square fit of g¯RBFsuperscript¯𝑔𝑅𝐵𝐹\bar{g}^{RBF}over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_R italic_B italic_F end_POSTSUPERSCRIPT to the observed velocity distribution. We again compute the RBF on each 128s average velocity distribution and include weighted errors in the least square fit corresponding to Poisson noise analogous to Eq B5.

B.3 Comparing Model Distributions

To verify the heating of the plasma via ICW resonant interactions, we recompute Equation 2 using the Hermite polynomial and radial basis function approximations to the proton distribution, g¯H(v,v)superscript¯𝑔𝐻subscript𝑣perpendicular-tosubscript𝑣parallel-to\bar{g}^{H}(v_{\perp},v_{\parallel})over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) and g¯RBF(v,v)superscript¯𝑔𝑅𝐵𝐹subscript𝑣perpendicular-tosubscript𝑣parallel-to\bar{g}^{RBF}(v_{\perp},v_{\parallel})over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_R italic_B italic_F end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ), we refer to the heating rates computed from Equation 2 via the drifting biMaxwellian, Hermite polynomials, and RBF interpolation as QBMsubscript𝑄𝐵𝑀Q_{BM}italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT, QHMTEsubscript𝑄𝐻𝑀𝑇𝐸Q_{HMTE}italic_Q start_POSTSUBSCRIPT italic_H italic_M italic_T italic_E end_POSTSUBSCRIPT, and QRBFsubscript𝑄𝑅𝐵𝐹Q_{RBF}italic_Q start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT respectively, we also compute the average heating rate of the three terms QAVGsubscript𝑄𝐴𝑉𝐺Q_{AVG}italic_Q start_POSTSUBSCRIPT italic_A italic_V italic_G end_POSTSUBSCRIPT. Note that the ICW spectrum I(k)𝐼subscript𝑘parallel-toI(k_{\parallel})italic_I ( italic_k start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) in unaffected by the choice of model as di=VA/ωcsubscript𝑑𝑖subscript𝑉𝐴subscript𝜔𝑐d_{i}=V_{A}/\omega_{c}italic_d start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT / italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is a function only of the magnetic field strength and total plasma density, which are independent of the kinetic phase-space distribution.

Refer to caption
Figure 6: a) 2D histogram of QHMTEsubscript𝑄𝐻𝑀𝑇𝐸Q_{HMTE}italic_Q start_POSTSUBSCRIPT italic_H italic_M italic_T italic_E end_POSTSUBSCRIPT against QBMsubscript𝑄𝐵𝑀Q_{BM}italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT for each of the 6737 studied intervals; data are column normalized and the black line shows QHMTE=QBMsubscript𝑄𝐻𝑀𝑇𝐸subscript𝑄𝐵𝑀Q_{HMTE}=Q_{BM}italic_Q start_POSTSUBSCRIPT italic_H italic_M italic_T italic_E end_POSTSUBSCRIPT = italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT. b) Same as panel a, but with QRBFsubscript𝑄𝑅𝐵𝐹Q_{RBF}italic_Q start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT against QBMsubscript𝑄𝐵𝑀Q_{BM}italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT. c) Same as panel a, but with the average of all three heating rates QAVGsubscript𝑄𝐴𝑉𝐺Q_{AVG}italic_Q start_POSTSUBSCRIPT italic_A italic_V italic_G end_POSTSUBSCRIPT against QBMsubscript𝑄𝐵𝑀Q_{BM}italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT. d) Histogram showing distributions of Pearson R2superscript𝑅2R^{2}italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT correlation between the 6737 modeled distributions and observed SPANi distributions; drifting biMaxwellian (BM) in black; Hermite (HMTE) in teal; radial basis function (RBF) in orange. e) Histogram showing distributions of χRBF2/χBM2superscriptsubscript𝜒𝑅𝐵𝐹2superscriptsubscript𝜒𝐵𝑀2\chi_{RBF}^{2}/\chi_{BM}^{2}italic_χ start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_χ start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (orange) and χHMTE2/χBM2superscriptsubscript𝜒𝐻𝑀𝑇𝐸2superscriptsubscript𝜒𝐵𝑀2\chi_{HMTE}^{2}/\chi_{BM}^{2}italic_χ start_POSTSUBSCRIPT italic_H italic_M italic_T italic_E end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_χ start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (teal) on a logarithmic scale.

The top row of Figure 6, panels (a-c), respectively show QHMTEsubscript𝑄𝐻𝑀𝑇𝐸Q_{HMTE}italic_Q start_POSTSUBSCRIPT italic_H italic_M italic_T italic_E end_POSTSUBSCRIPT, QRBFsubscript𝑄𝑅𝐵𝐹Q_{RBF}italic_Q start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT, and QAVGsubscript𝑄𝐴𝑉𝐺Q_{AVG}italic_Q start_POSTSUBSCRIPT italic_A italic_V italic_G end_POSTSUBSCRIPT against QBMsubscript𝑄𝐵𝑀Q_{BM}italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT. While all computed heating rates have the same qualitative behavior, i.e. positive valued indicating the absorption of ICWs and plasma heating via cyclotron resonance, QHMTEsubscript𝑄𝐻𝑀𝑇𝐸Q_{HMTE}italic_Q start_POSTSUBSCRIPT italic_H italic_M italic_T italic_E end_POSTSUBSCRIPT tends to systematically underestimate the heating rates predicted by QBMsubscript𝑄𝐵𝑀Q_{BM}italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT, whereas QRBFsubscript𝑄𝑅𝐵𝐹Q_{RBF}italic_Q start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT tends to overestimate QBMsubscript𝑄𝐵𝑀Q_{BM}italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT. The average of all heating rates QAVGsubscript𝑄𝐴𝑉𝐺Q_{AVG}italic_Q start_POSTSUBSCRIPT italic_A italic_V italic_G end_POSTSUBSCRIPT closely follows QBMsubscript𝑄𝐵𝑀Q_{BM}italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT, thus QBMsubscript𝑄𝐵𝑀Q_{BM}italic_Q start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT was chosen as the value reported in the body of the manuscript.

The bottom panels shows statistics on the quality of the various models. To compare the quality of the relative models we compute two parameters. Figure 6(d) shows the distribution of measured Pearson correlation coefficients between the logarithm of the observed distribution Log10gp(v,v)subscriptLog10subscript𝑔𝑝subscript𝑣perpendicular-tosubscript𝑣parallel-to\textrm{Log}_{10}g_{p}(v_{\perp},v_{\parallel})Log start_POSTSUBSCRIPT 10 end_POSTSUBSCRIPT italic_g start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ) and the three models Log10g¯BM(v,v)subscriptLog10superscript¯𝑔𝐵𝑀subscript𝑣perpendicular-tosubscript𝑣parallel-to\textrm{Log}_{10}\bar{g}^{BM}(v_{\perp},v_{\parallel})Log start_POSTSUBSCRIPT 10 end_POSTSUBSCRIPT over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_B italic_M end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ), Log10g¯H(v,v)subscriptLog10superscript¯𝑔𝐻subscript𝑣perpendicular-tosubscript𝑣parallel-to\textrm{Log}_{10}\bar{g}^{H}(v_{\perp},v_{\parallel})Log start_POSTSUBSCRIPT 10 end_POSTSUBSCRIPT over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ), and Log10g¯RBF(v,v)subscriptLog10superscript¯𝑔𝑅𝐵𝐹subscript𝑣perpendicular-tosubscript𝑣parallel-to\textrm{Log}_{10}\bar{g}^{RBF}(v_{\perp},v_{\parallel})Log start_POSTSUBSCRIPT 10 end_POSTSUBSCRIPT over¯ start_ARG italic_g end_ARG start_POSTSUPERSCRIPT italic_R italic_B italic_F end_POSTSUPERSCRIPT ( italic_v start_POSTSUBSCRIPT ⟂ end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT ∥ end_POSTSUBSCRIPT ). High levels of correlation indicate that the model correctly matches the observed distribution. The use of the logarithm of the distribution in computing the Pearson R2superscript𝑅2R^{2}italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT correlation coefficient allows for weighting of the entirety of phase space and includes the regions where ICW resonant interactions are most important; without the logarithmic weighting, the correlation coefficient is mostly dominated by the models performance at the peak of the distribution. We find R2superscript𝑅2R^{2}italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT is largest for the RBF interpolation and smallest for the parametric drifting biMaxwellian fit. Interestingly, the level of heating computed between these models remains relatively constant even though the R2superscript𝑅2R^{2}italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT values vary, this suggests that variations between the three models are mostly occurring in regions that are not resonant with the outward going ICWs. We speculate these differences are likely due to the models ability to approximate the ion-beam population.

We also compute the χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT of each of the models. Given the nature of the nonlinear fitting there is lack of clear number of degrees of freedom. Furthermore errors in the model are not likely to obey Gaussian statistics. For these reasons hypothesis testing via the χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT statistic is not possible. To determine the quality of the fits we compare the ratio of χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT to the value compared for the biMaxwellian, χBM2superscriptsubscript𝜒𝐵𝑀2\chi_{BM}^{2}italic_χ start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Figure 6(e) shows the distribution of the logarithmic ratio of χRBF2/χBM2superscriptsubscript𝜒𝑅𝐵𝐹2superscriptsubscript𝜒𝐵𝑀2\chi_{RBF}^{2}/\chi_{BM}^{2}italic_χ start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_χ start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and χHMTE2/χBM2superscriptsubscript𝜒𝐻𝑀𝑇𝐸2superscriptsubscript𝜒𝐵𝑀2\chi_{HMTE}^{2}/\chi_{BM}^{2}italic_χ start_POSTSUBSCRIPT italic_H italic_M italic_T italic_E end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_χ start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. The RBF model tends to perform much netter than the drifting biMaxwellian with a uniform χRBF2/χBM2<1superscriptsubscript𝜒𝑅𝐵𝐹2superscriptsubscript𝜒𝐵𝑀21\chi_{RBF}^{2}/\chi_{BM}^{2}<1italic_χ start_POSTSUBSCRIPT italic_R italic_B italic_F end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_χ start_POSTSUBSCRIPT italic_B italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT < 1. In terms of χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, the Hermite polynomial approximation tends to have a similar performance to the drifting biMaxwellian fit in this stream, though the Pearson correlation is significantly better than the biMaxwellian fit. Again, the ability of each of models to reproduce similar heating rates, but with significantly different χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT suggests that the main difference in the models may be in reproducing the ion-beam. Our future work will investigate further the role of the ion-beam in these processes.

In any case, the use of these various parametric and non-parametric models for the observed proton distribution function all recover heating rates that suggest that ICW are absorbed into the core of the proton distribution, resulting in plasma heating via ICW resonance. These results suggest that the model used may not be inherently important in understanding ICW heating of the proton core.