1. Introduction
Many natural plasmas (e.g. interstellar medium, galaxy clusters, black hole accretion disks, solar wind) are in a weakly collisional state, where the particle collision frequency $\nu _{\mathrm {coll}}$ is smaller than other characteristic frequencies (e.g. proton gyrofrequency $\varOmega _{\mathrm {p}}$, inverse magnetic-field correlation time $1/\tau _c$ etc. (Quataert Reference Quataert2003; Marsch Reference Marsch2006; Schekochihin & Cowley Reference Schekochihin and Cowley2006; Schekochihin et al. Reference Schekochihin, Cowley, Dorland, Hammett, Howes, Quataert and Tatsuno2009). Thus, the regimes of the plasma motions, characterised by $\omega$, can span from collisional (fluid) $\omega \ll \nu _{\mathrm {coll}}$ to collisionless $\omega \gg \nu _{\mathrm {coll}}$, where $\omega$ is the temporal frequency of a fluctuating plasma property (Schekochihin et al. Reference Schekochihin, Cowley, Kulsrud, Hammett and Sharma2005, Reference Schekochihin, Cowley, Dorland, Hammett, Howes, Quataert and Tatsuno2009). Knowledge of the transition scale $\omega \sim \nu _{\mathrm {coll}}$ is vital to understand the behaviour of astrophysical plasmas.
The escaping solar corona, known as the solar wind, expands into interplanetary space as a super-Alfvénic and turbulent plasma (Parker Reference Parker1958; Bruno & Carbone Reference Bruno and Carbone2013; Verscharen, Klein & Maruca Reference Verscharen, Klein and Maruca2019). In situ measurements of particle distribution functions and electromagnetic fields enable fundamental plasma physics observations (Chen Reference Chen2016). The Spitzer–Härm proton–proton collision frequency $\nu _{\mathrm {p},\mathrm {p}}^{\mathrm {SH}}$ decreases with radial distance from the Sun and, by a few solar radii, is much smaller than other characteristic frequencies. In principle, the dynamics should be described by collisionless plasma equations. For reference, at 1 AU, typical frequencies are $\nu _{\mathrm {p},\mathrm {p}}^{\mathrm {SH}} \approx 4 \times 10^{-7} \ \mathrm {s}^{-1}$, $\varOmega _{\mathrm {p}} \approx 10^{-1} \ \mathrm {s}^{-1}$, $1/\tau _c \approx 10^{-6} \ \mathrm {s}^{-1}$ (see Appendix A for the calculation of the collision frequency; Spitzer Reference Spitzer1962; Huba Reference Huba1983; Matthaeus et al. Reference Matthaeus, Dasso, Weygand, Kivelson and Osman2010; Verscharen et al. Reference Verscharen, Klein and Maruca2019).
Despite the weak collisionality of the solar wind, many aspects appear to be described by fluid equations: magnetohydrodynamic (MHD) turbulence theory predicts the shape of power spectra (e.g. magnetic field, proton density; Coleman Reference Coleman1968; Matthaeus & Goldstein Reference Matthaeus and Goldstein1982; Tu & Marsch Reference Tu and Marsch1995; Goldreich & Sridhar Reference Goldreich and Sridhar1997; Bruno & Carbone Reference Bruno and Carbone2013), spatial transport (Zank, Matthaeus & Smith Reference Zank, Matthaeus and Smith1996; Matthaeus et al. Reference Matthaeus, Zank, Smith and Oughton1999) and the proton heating rate by the energy cascade (MacBride, Smith & Forman Reference MacBride, Smith and Forman2008; Stawarz et al. Reference Stawarz, Smith, Vasquez, Forman and MacBride2009; Coburn et al. Reference Coburn, Smith, Vasquez, Stawarz and Forman2012). This success is certainly due, in part, to the dominance of Alfvénic fluctuations, which have identical properties in the collisionless and fluid limits on scales above the proton gyroradius, implying much of the turbulent energy cascade is basically insensitive to the plasma's collisionality (Kulsrud, Sagdeev & Rosenbluth Reference Kulsrud, Sagdeev and Rosenbluth1980; Schekochihin et al. Reference Schekochihin, Cowley, Dorland, Hammett, Howes, Quataert and Tatsuno2009). However, observations show a strong correlation between the density and thermal pressure (i.e. compressive fluctuations), indicating a polytropic equation of state (Marsch et al. Reference Marsch, Muehlhaeuser, Rosenbauer and Schwenn1983; Totten, Freeman & Arya Reference Totten, Freeman and Arya1995; Verscharen, Chen & Wicks Reference Verscharen, Chen and Wicks2017; Nicolaou et al. Reference Nicolaou, Livadiotis, Wicks, Verscharen and Maruca2020). Moreover, compressive fluctuations generally display the MHD slow-mode polarisation (anticorrelated magnetic and thermal pressure; Verscharen et al. Reference Verscharen, Chen and Wicks2017), rather than being severely damped, as expected in a collisionless plasma (Barnes Reference Barnes1966). These modes are routinely detected at a range of scales (Tu & Marsch Reference Tu and Marsch1995; Kellogg & Horbury Reference Kellogg and Horbury2005; Yao et al. Reference Yao, He, Marsch, Tu, Pedersen, Rème and Trotignon2011; Howes et al. Reference Howes, Bale, Klein, Chen, Salem and TenBarge2012; Klein et al. Reference Klein, Howes, TenBarge, Bale, Chen and Salem2012; Yao et al. Reference Yao, He, Tu, Wang and Marsch2013a,Reference Yao, He, Tu, Wang and Marschb) following a power law predicted from the MHD equations (Montgomery, Brown & Matthaeus Reference Montgomery, Brown and Matthaeus1987; Marsch & Tu Reference Marsch and Tu1990; Lithwick & Goldreich Reference Lithwick and Goldreich2001; Schekochihin et al. Reference Schekochihin, Cowley, Dorland, Hammett, Howes, Quataert and Tatsuno2009).
While the Spitzer–Härm collision frequency appears incompatible with the fluid-like behaviour of the solar wind, weakly collisional plasmas are also subject to collisionless relaxation processes that prevent extreme departure from equilibrium (Nishida Reference Nishida1969; Griffel & Davis Reference Griffel and Davis1969; Hamasaki & Krall Reference Hamasaki and Krall1973; Gary, Yin & Winske Reference Gary, Yin and Winske2000; Yoon Reference Yoon2017). Solar wind observations present substantial evidence of temperature anisotropy instabilities constraining the particle distribution functions (Kasper, Lazarus & Gary Reference Kasper, Lazarus and Gary2002; Tu & Marsch Reference Tu and Marsch2002; Marsch Reference Marsch2006; Hellinger et al. Reference Hellinger, Trávníček, Kasper and Lazarus2006; Bale et al. Reference Bale, Kasper, Howes, Quataert, Salem and Sundkvist2009; Chen et al. Reference Chen, Matteini, Schekochihin, Stevens, Salem, Maruca, Kunz and Bale2016; Yoon Reference Yoon2017). These processes can play a similar role to collisions viz., they are effective collision processes.
This Letter presents a measurement of the effective mean free path of the solar wind by comparing observations of compressive wave-mode polarisation with numerical solutions of varying effective collision frequency. It is shown that the transition from fluid to collisionless dynamics in the solar wind occurs at scales several orders of magnitude below the classical Spitzer–Härm mean free path, explaining the fluid-like behaviour of the weakly collisional solar wind.
2. Theory and numerical solutions
The kinetic MHD equations with the Bhatnagar–Gross–Krook (BGK) collision operator (Bhatnagar, Gross & Krook Reference Bhatnagar, Gross and Krook1954; Gross & Krook Reference Gross and Krook1956) produce dispersion relations and plasma fluctuations (e.g. magnetic field and pressure) that span between the collisionless and collisional limits (Kulsrud et al. Reference Kulsrud, Sagdeev and Rosenbluth1980; Snyder, Hammett & Dorland Reference Snyder, Hammett and Dorland1997; Sharma, Hammett & Quataert Reference Sharma, Hammett and Quataert2003; Chandran et al. Reference Chandran, Dennis, Quataert and Bale2011). They describe a non-relativistic, magnetised plasma of arbitrary collision frequency (Kulsrud et al. Reference Kulsrud, Sagdeev and Rosenbluth1980; Sharma et al. Reference Sharma, Hammett and Quataert2003). See Appendix B for details on the equations and BGK collision operator. The specific equations, which we refer to as kinetic MHD-BGK (KMHD-BGK), model both the proton and electron responses with the kinetic equation so that Landau and Barnes dampings are retained. The BGK operator is used here to model relaxation processes, not particle collisions, so we use the language of an effective proton collision frequency $\nu _{\mathrm {eff}}$ or mean free path $\lambda _{\mathrm {mfp}}^{\mathrm {eff}} = v_{\mathrm {th}}^{\mathrm {p}} / \nu _{\mathrm {eff}}$, where the proton thermal speed is $v_{\mathrm {th}}^{\mathrm {p}}$.
Assuming plasma motions are slow compared with the gyrofrequency $\varOmega _p$, the second moment of the kinetic equation and the ideal induction equation leads to
where $\mathrm {d}/\mathrm {d}t$ is the convective derivative and the quantities are the proton density $n_{\mathrm {p}}$, magnetic-field strength $B$, parallel (perpendicular) proton pressure $p_{\parallel }^{\mathrm {p}} \ (p_{\perp }^{\mathrm {p}})$, field parallel flux of parallel (perpendicular) proton heat $q_{\parallel }^{\mathrm {p}} \ (q_{\perp }^{\mathrm {p}})$ and the unit magnetic-field vector $\hat {\boldsymbol {b}} = \boldsymbol {B} / B$ (Chew, Goldberger & Low Reference Chew, Goldberger and Low1956; Hunana et al. Reference Hunana, Tenerani, Zank, Khomenko, Goldstein, Webb, Cally, Collados, Velli and Adhikari2019). The Alfvén speed is $v_A= B/\sqrt {4 {\rm \pi}n_{\mathrm {p}} m_{\mathrm {p}}}$, the proton gyroradius is $\rho _{\mathrm {p}} = v_{\mathrm {th}}^{\mathrm {p}} / \varOmega _{\mathrm {p}}$ and the ion-acoustic speed is $c_{\mathrm {s}} = \sqrt {(3 k_{\mathrm {B}} T_{\parallel }^{\mathrm {p}} + k_{\mathrm {B}} T_{\parallel }^{\mathrm {e}})/m_{\mathrm {p}}}$, where the parallel proton (electron) temperature is $T_{\parallel }^{\mathrm {p}} (T_{\parallel }^{\mathrm {e}})$.
Equations (2.1) are often discussed when the right-hand sides are zero and are then referred to as the double adiabatic equations or Chew–Goldberger–Low (CGL) invariants (Chew et al. Reference Chew, Goldberger and Low1956). The focus here is on how the CGL invariants are broken, for example, by the heat flux terms in the collisionless limit, and by the effective collisional terms ($\propto \nu _{\mathrm {eff}}$). Therefore, the relative non-conservation of the CGL invariants provides a sensitive test of the equation of state.
Cross-correlations and amplitude ratios, which can be measured, are constructed from the left-hand sides of (2.1),
where $\delta \chi = \chi - \langle \chi \rangle$ is the fluctuation about the average $\langle \chi \rangle$. They describe the relative non-conservation of the CGL invariants. The method compares predictions for (2.2) derived from the slow-mode eigenmodes of the linearised KMHD-BGK system (e.g. $\delta p_{\perp }^{\mathrm {p}}, \, \delta B$ etc.) with solar wind measurements, since the slow mode is the dominant compressive mode at large scales (Howes et al. Reference Howes, Bale, Klein, Chen, Salem and TenBarge2012). A description of the equations and the linear system of equations appears in Appendix B.
The model's free parameters are the propagation angle $\theta _{\hat {\boldsymbol {b}},\hat {\boldsymbol {k}}}$ and proton effective mean free path $\lambda _{\mathrm {mfp}}^{\mathrm {eff}}$; they are determined by fitting to solar wind observations. The wavenumber $k$ and the proton beta $\beta = 8 {\rm \pi}p^{\mathrm {p}} / B^{2}$, where $p^{\mathrm {p}} = 2p^{\mathrm {p}}_{\perp }/3 + p^{\mathrm {p}}_{\parallel }/3$, are set to measured values. The species temperature ratio is set to a typical value for the solar wind $T_{\mathrm {p}} / T_{\mathrm {e}} = 1$ and the effective mean free path species ratio is set to $\lambda _{\mathrm {mfp}}^{\mathrm {eff}} / \lambda _{\mathrm {mfp,electrons}}^{\mathrm {eff}} = 1$ (see Appendix B for details). At small $\beta$, these two ratios $T_{\mathrm {p}} / T_{\mathrm {e}}, \, \lambda _{\mathrm {mfp}}^{\mathrm {eff}} / \lambda _{\mathrm {mfp,electrons}}^{\mathrm {eff}}$ have an insignificant influence on the correlations, (2.2), and nearly no influence at large $\beta$.
Figure 1 demonstrates the ability of the KMHD-BGK equations to resolve the dynamics of the compressive slow mode from collisional (lighter blue) to collisionless (black). Numerical predictions for (2.2a), (2.2c) (bottom panels of figure 1) show distinct differences at $\beta > 1$ for different $k_{\parallel } \lambda _{\mathrm {mfp}}^{\mathrm {eff}}$, which can be compared with observations. The MHD ($k_{\parallel } \lambda _{\mathrm {mfp}}^{\mathrm {eff}} \ll 1$) and collisionless ($k_{\parallel } \lambda _{\mathrm {mfp}}^{\mathrm {eff}} \gg 1$) limits are illustrated in magenta, for $C_{\perp }$ (bottom right panel) these two limits produce similar trends, therefore it is necessary to make comparisons at multiple $k_{\parallel }$ to measure $\lambda _{\mathrm {mfp}}^{\mathrm {eff}}$.
3. Measurements
The dataset consists of Wind spacecraft measurements of the pristine solar wind during years 2005–2010. The electrostatic analyser, 3DP, records onboard moments of the proton density, velocity and pressure tensor, and the magnetometer MFI records the magnetic field, at a nominal ${\sim }3s$ cadence (Lepping et al. Reference Lepping, Acũna, Burlaga, Farrell, Slavin, Schatten, Mariani, Ness, Neubauer and Whang1995; Lin et al. Reference Lin, Anderson, Ashford, Carlson, Curtis, Ergun, Larson, McFadden, McCarthy and Parks1995).
The dataset is restricted to time intervals satisfying three criteria: (i) 95 % of the data are available (the remaining are then linearly interpolated); (ii) the median density must be greater than 1 particle per cm$^{3}$; and (iii) the average norm of the non-gyrotropic tensor ($\mathit{\boldsymbol{\Pi}} ^{\mathrm {p}} =\boldsymbol{ \mathsf{ p}}^{\mathrm {p}} - \hat {\boldsymbol {b}} \hat {\boldsymbol {b}} \, p^{\mathrm {p}}_{\parallel } - (1 - \hat {\boldsymbol {b}} \hat {\boldsymbol {b}} ) p_{\perp }^{\mathrm {p}}$), must be less than 30 % of the average norm of the pressure tensor $\boldsymbol{ \mathsf{ p}}^{\mathrm {p}}$. The final point here is satisfied $\sim$94 % of the time at the ${\sim }3s$ time interval and more so at longer time intervals.
To probe a set of wavenumbers, the four quantities in (2.2) are measured, along with the average radial solar wind velocity $\langle V_{\mathrm {SW}} \rangle$ and the average proton beta for a set of time intervals $\tau = [30\ \mathrm {s}, 1\ \mathrm {min.}, 2\ \mathrm {mins.}, \ldots, 128\ \mathrm {mins.}]$. The time scales are converted to wavenumber $k_{\mathrm {SW}} = 1/ \tau \langle V_{\mathrm {SW}} \rangle$ via Taylor's frozen-in-flow (TFF) assumption (Taylor Reference Taylor1938). Outliers in the distribution of $k_{\mathrm {SW}}$ are removed and then three bins of equal probability density are obtained where the median of each bin is $k_{\mathrm {SW}} = [0.288, 1.41, 6.34] \times 10^{-5} \ \mathrm {km}^{-1}$, which lies within the inertial range of the magnetic-field power spectrum at 1 AU (Kiyani, Osman & Chapman Reference Kiyani, Osman and Chapman2015). The choice of bins provides enough separation in wavenumber to resolve differences in the measured quantities (e.g. figure 1) and sufficient sampling; the wavenumber bins contain $[2.98, 16.6, 70.0] \times 10^{5}$ samples.
For bin $k_{\mathrm {SW}} = 0.288 \times 10^{-5} \ \mathrm {km}^{-1}$ the $\beta$-conditioned probability functions of (2.2), mapped to a common colour bar, are displayed in figure 2. The $\beta$-trend lines in magenta (see caption) capture statistically significant differences between $\beta \lessgtr 1$. From the correlations $C_{\parallel }, C_{\perp }$ it is clear that the CGL invariants are rarely conserved $(C_{\parallel }, C_{\perp } = 1)$, but display similar trends to the theoretical expectations seen in figure 1. The amplitude ratios $A_{\parallel }, A_{\perp }$ demonstrate a relative decrease in fluctuation amplitude of the pressure components at $\beta > 1$.
4. Comparison of measurements and numerical solutions
The theoretical predictions for (2.2) from the numerical model in § 2 are compared with the observations in figure 2 to determine the most probable effective mean free path $\lambda _{\mathrm {mfp}}^{\mathrm {eff}}$ and propagation angles $\theta _{\hat {\boldsymbol {b}},\hat {\boldsymbol {k}}}$. The numerical predictions can be fitted to the observations by altering the parameters (e.g. the effective mean free path), but a degeneracy in parametrisation must be dealt with. The numerical solutions primarily depend on $k_{\parallel } \lambda _{\mathrm {mfp}}^{\mathrm {eff}} = k \cos (\theta _{\hat {\boldsymbol {b}},\hat {\boldsymbol {k}}}) \lambda _{\mathrm {mfp}}^{\mathrm {eff}}$ (Sharma et al. Reference Sharma, Hammett and Quataert2003), which implies that $\lambda _{\mathrm {mfp}}^{\mathrm {eff}}$ and $\theta _{\hat {\boldsymbol {b}},\hat {\boldsymbol {k}}}$ are degenerate. To address this, a scale-dependent anisotropy model $(k_{\parallel } \sim k_{\perp }^{\alpha })$ is introduced, which relates $k$ and $\theta _{\hat {\boldsymbol {b}},\hat {\boldsymbol {k}}}$,
where $k_{\mathrm {iso}}$ is the isotropic wavenumber ($k = k_{\mathrm {iso}}$ when $k_{\perp } = k_{\parallel }$) and $\alpha$ is the anisotropy exponent, generalised from turbulence models (Goldreich & Sridhar Reference Goldreich and Sridhar1995). The wavenumber model provides $\theta _{\hat {\boldsymbol {b}},\hat {\boldsymbol {k}}}$, given $k$, parameterised by $\alpha, k_{\mathrm {iso}}$, so that $\lambda _{\mathrm {mfp}}^{\mathrm {eff}}$ can be determined at the given wavenumber. Then comparing solutions parameterised by $\lambda _{\mathrm {mfp}}^{\mathrm {eff}}$, $\alpha$, $k_{\mathrm {iso}}$ at multiple wavenumbers $k = k_{\mathrm {SW}}$ clears the degeneracy and allows all three parameters to be measured. Additionally, the model permits quantification of the observed increase of obliqueness with wavenumber of the compressive fluctuations (Chen et al. Reference Chen, Mallet, Schekochihin, Horbury, Wicks and Bale2012; Chen Reference Chen2016).
Finally, the predictions of (2.2) from the numerical solutions are normalised to the measured $\beta$-conditioned mean value (dashed magenta lines in figure 2) of $C_{\parallel }, C_{\perp }, A_{\parallel }, A_{\perp }$ at $\beta \simeq 10^{-1}$. This is to account for the fact that linear wave properties are only approximately observed in strong turbulence (Chen Reference Chen2016; Grošelj et al. Reference Grošelj, Chen, Mallet, Samtaney, Schneider and Jenko2019).
The ranges $\alpha = [0.05,1.0]$, $k_{\mathrm {iso}} = [5 \times 10^{-9}, 5 \times 10^{-7}]\ \mathrm {km}^{-1}$ and $\lambda _{\mathrm {mfp}}^{\mathrm {eff}} = [3.5 \times 10^{4}, 2.1 \times 10^{6}] \ \mathrm {km}$ are chosen for computing numerical solutions. The ranges of $\alpha, k_{\mathrm {iso}}$ are consistent with previous observations (Chen et al. Reference Chen, Mallet, Schekochihin, Horbury, Wicks and Bale2012; Chen Reference Chen2016). The range of $\lambda _{\mathrm {mfp}}^{\mathrm {eff}}$ returns numerical solutions of (2.2) that compare qualitatively well with the observations (seen in figure 2). The Spitzer–Härm mean free path returns the (collisionless) ion-acoustic dispersion relation which is inconsistent with the measurements.
To make a quantitative comparison, we compute the ‘goodness of fit’
where $\hat {y}_i$ ($\bar {y}_i$) is the local numerical solution (local measured mean), summed over $i$, denoting the $i$th $\beta$-bin. Here, $R(k_{\mathrm {SW}}; \alpha, k_{\mathrm {iso}}, \lambda _{\mathrm {mfp}}^{\mathrm {eff}})$ is calculated for each wavenumber $k_{\mathrm {SW}}$, where the mean $\bar {y}_i$ is respective to the wavenumber bin. The $R$-values are inverted for unnormalised weights ($w = R^{-1}$), divided by the maximum weight ($w_{\mathrm {max}}$), then summed over wavenumber $\mathcal {W}(\alpha, k_{\mathrm {iso}}, \lambda _{\mathrm {mfp}}^{\mathrm {eff}}) = \sum _k w(k; \alpha, k_{\mathrm {iso}}, \lambda _{\mathrm {mfp}}^{\mathrm {eff}})/ w_{\mathrm {max}}$ to break the aforementioned degeneracy. From this volume, weighted geometric means $\mu _x$, covariances $\sigma ^{2}_{x,y}$, and two sigma confidence intervals $\mathrm {CI}_{x}$ are calculated (see Appendix C for the statistics; Norris Reference Norris1940; Kendall & Stuart Reference Kendall and Stuart1977). This is the method employed to measure the quantities $\lambda _{\mathrm {mfp}}^{\mathrm {eff}}, \alpha, k_{\mathrm {iso}}$, providing the main results of the Letter.
To visualise the weighted parameter space for $C_{\perp }$, figure 3 illustrates the weight volume $\mathcal {W}(\alpha, k_{\mathrm {iso}}, \lambda _{\mathrm {mfp}}^{\mathrm {eff}})$ numerically integrated over each parameter axis $\chi$,
where $\chi _n, \chi _0$ are limits of the range. The weighted means in figure 3 lie in the maximum regions of $\mathcal {W}_{\chi }$, within the confidence intervals, indicating the weighted geometric statistics are a good representation of the observations.
To check the scale dependence, figure 4 displays the observed $\beta$-conditioned means of (2.2) and the numerical solutions corresponding to the maximum $\mathcal {W}$. The numerical solutions and observations trend similarly with wavenumber, indicating the scale dependence of the effective collisionality has been well modelled. The parameters of the maximum (recorded in figure 4) do not correspond exactly to the weighted geometric means of $\mathcal {W}$ (seen in figure 3), reflecting the statistical nature of the measured quantities.
The method of calculating statistics for $\lambda _{\mathrm {mfp}}^{\mathrm {eff}}$, $\alpha$, $k_{\mathrm {iso}}$ displayed in figure 3 for $C_{\perp }$ produces similar statistics for $C_{\parallel }, A_{\parallel }, A_{\perp }$ (see Appendix C). Therefore, in table 1, combined statistics are reported. The measured effective mean-free-path and mean proton thermal speed (measured with this dataset) give an effective collision frequency of $\nu _{\mathrm {eff}} = v_{\mathrm {th}}^{\mathrm {p}} / \lambda _{\mathrm {mfp}}^{\mathrm {eff}}$ = 1.11 $\times 10^{-4}$ s$^{-1}$.
The transition frequency, where $\nu _{\mathrm {eff}} \simeq \omega$, can be estimated with $\nu _{\mathrm {eff}}$ and the ion-acoustic dispersion relation $\omega _{\mathrm {IA}} = k_{\parallel } \, c_{\mathrm {s}}$, giving the parallel transition wavenumber $k^{\mathrm {trans}}_{\parallel } = v_{\mathrm {th}}^{\mathrm {p}}/c_{\mathrm {s}} \, \lambda _{\mathrm {mfp}}^{\mathrm {eff}}$. Using the wavenumber model, (4.1), the transition wavenumber is
(more details are provided in Appendix D). Inserting the combined statistics from table 1, using a typical value of $v_{\mathrm {th}}^{\mathrm {p}} / c_{\mathrm {s}} = \sqrt {1/2}$ for the solar wind, and using the TFF assumption, the transition wavenumber in spacecraft-frame frequency at 1 AU is $\langle V_{\mathrm {SW}} \rangle \, k^{\mathrm {trans}} = f^{\mathrm {trans}} =0.19$ Hz and $\mathrm {CI}_{f^{\mathrm {trans}}} = [0.046, 0.33 ]$ Hz. The uncertainties are propagated from $V_{\mathrm {SW}}$ and the four estimates of $k^{\mathrm {trans}}$ from $C_{\parallel }, \, A_{\parallel }, \, C_{\perp }, \, A_{\perp }$.
5. Discussion
The relative non-conservation of the CGL invariants has been measured and the behaviour has been modelled with the slow-mode branch of the KMHD-BGK equations to measure the effective mean free path of the solar wind protons and the scale dependence of the slow-mode wavenumber anisotropy (table 1 reports the statistics of these measurements). The primary result of this Letter is the measured effective proton mean free path, which is ${\sim }10^{3}$ times smaller than the Spitzer–Härm mean free path ($\lambda _{\mathrm {mfp}}^{\mathrm {SH}} = 1.14 \times 10^{8} \ \mathrm {km}$, see Appendix A). Therefore, the fluid-like range in the solar wind extends to much smaller scales than the prediction based on particle collisions. In addition, the scale-dependent anisotropy of the compressive fluctuations ($\alpha \simeq 0.4$) is consistent with previous measurements (Chen et al. Reference Chen, Mallet, Schekochihin, Horbury, Wicks and Bale2012; Chen Reference Chen2016), being more anisotropic than the Alfvénic fluctuations, which follow the critical balance value $\alpha _{\mathrm {CB}} = 2/3$ (Goldreich & Sridhar Reference Goldreich and Sridhar1997).
The measured transition frequency, the scale between fluid behaviour ($f \ll f^{\mathrm {trans}}$) and collisionless behaviour ($f \gg f^{\mathrm {trans}}$), of $f^{\mathrm {trans}} = 0.19$ Hz is at the well-known break in power law $(k_{\perp } \rho _{\mathrm {p}} \sim 1)$ of the magnetic-field power spectrum at 1 AU (Leamon et al. Reference Leamon, Smith, Ness, Matthaeus and Wong1998; Kiyani et al. Reference Kiyani, Osman and Chapman2015; Verscharen et al. Reference Verscharen, Klein and Maruca2019). These measurements therefore justify the use of fluid MHD theory at larger scales $(k_{\perp } \rho _{\mathrm {p}} < 1)$ (Tu & Marsch Reference Tu and Marsch1995; Goldreich & Sridhar Reference Goldreich and Sridhar1997; Bruno & Carbone Reference Bruno and Carbone2013; Chen Reference Chen2016). If the result $k_{\perp } \rho _{\mathrm {p}} \simeq k_{\parallel } \lambda ^{\mathrm {eff}}_{\mathrm {mfp}}$ turns out to be a general property of weakly collisional plasma, this provides a simple parameterisation for the effective collisionality of astrophysical plasmas.
The method employed in this Letter relies on linear theory to describe cross-correlations and amplitude ratios of various plasma properties. While there is evidence that such linear quantities are approximately preserved in solar wind turbulence (e.g. Chen Reference Chen2016; Verscharen et al. Reference Verscharen, Chen and Wicks2017; Grošelj et al. Reference Grošelj, Chen, Mallet, Samtaney, Schneider and Jenko2019), and the numerical predictions in this Letter match the observations well, it would be interesting to study in the future the degree to which nonlinear effects may also contribute to these correlations. The BGK operator has been used in this work as a simple way to model effective collision processes, however, it is possible that the use of alternative collision operators may alter the numerical value of the effective mean free path obtained, although we would expect the order of magnitude result to hold. It should also be noted that the transition scale obtained here is consistent with previous observations of fluid-like behaviour in the solar wind above the ion gyroscale (Verscharen et al. Reference Verscharen, Chen and Wicks2017).
An important neglected effect is heating, which is needed to describe the solar wind temperature profile (Verma, Roberts & Goldstein Reference Verma, Roberts and Goldstein1995; Vasquez et al. Reference Vasquez, Smith, Hamilton, MacBride and Leamon2007). Heating would also break the CGL invariants through additional terms appearing on the right-hand side of (2.1) (Chandran et al. Reference Chandran, Dennis, Quataert and Bale2011; Hellinger et al. Reference Hellinger, Trávníček, Štverák, Matteini and Velli2013). The importance of this simplification can be addressed by comparing our measured effective collision time with the heating time inferred in previous papers (Vasquez et al. Reference Vasquez, Smith, Hamilton, MacBride and Leamon2007; Hellinger et al. Reference Hellinger, Trávníček, Štverák, Matteini and Velli2013). The heating time can be estimated from $Q_{\mathrm {heat}} / k_{\mathrm {B}} T_{\mathrm {p}} = \nu ^{\mathrm {heat}} \approx 10^{-6}\ {\rm s}^{-1}$, which is approximately $10^{2}$ times smaller than the effective collision frequency ($\nu _{\mathrm {eff}}$) measured in this Letter. This suggests isotropisation dominates over heating, which is not unreasonable given the extreme anisotropic evolution the CGL invariants would dictate due to expansion ($T_{\perp }^{\mathrm {p}} \propto R^{-2}, \; T_{\parallel }^{\mathrm {p}} \propto \mathrm {const.}$, where $R$ is the radial distance from the Sun; Matteini et al. Reference Matteini, Hellinger, Landi, Trávníček and Velli2012).
There are many possible mechanisms that could lead to the observed effective collision frequency. It is well known that large departures from the Maxwellian velocity distribution function invalidate the Spitzer–Härm approach (Marsch Reference Marsch2006), in particular, large gradients in velocity space lead to fast collisional thermalisation (Pezzi, Valentini & Veltri Reference Pezzi, Valentini and Veltri2016). However, given that the departure from Spitzer–Härm is a factor ${\approx }10^{3}$, and the fact that our analysis is based on the properties of the low-order moments of the velocity distribution function (which would not be strongly influenced by collisional effects; Pezzi et al. Reference Pezzi, Perrone, Servidio, Valentini, Sorriso-Valvo and Veltri2019), collisional processes are very unlikely to lead to the observed behaviour, implicating collisionless physics.
Previous works have suggested a range of possible effective collision processes that can arise in collisionless plasmas. These include wave–particle interactions (Kellogg Reference Kellogg2000; Graham et al. Reference Graham, Khotyaintsev, André, Vaivads, Divin, Drake, Norgren, Contel, Lindqvist and Rager2022), instabilities (Gary et al. Reference Gary, Yin and Winske2000; Yoon Reference Yoon2017) and the plasma wave echo (Schekochihin et al. Reference Schekochihin, Parker, Highcock, Dellar, Dorland and Hammett2016; Meyrand et al. Reference Meyrand, Kanekar, Dorland and Schekochihin2019). They have long been studied theoretically and numerically (Coroniti & Eviatar Reference Coroniti and Eviatar1977; Schekochihin & Cowley Reference Schekochihin and Cowley2006; Kunz, Schekochihin & Stone Reference Kunz, Schekochihin and Stone2014; Helander, Strumik & Schekochihin Reference Helander, Strumik and Schekochihin2016; Kunz, Stone & Quataert Reference Kunz, Stone and Quataert2016; Rincon et al. Reference Rincon, Califano, Schekochihin and Valentini2016; Squire, Quataert & Kunz Reference Squire, Quataert and Kunz2017), but it is an open question as to the relevant role of the various mechanisms and how they are activated (Verscharen et al. Reference Verscharen, Chandran, Klein and Quataert2016; Squire et al. Reference Squire, Kunz, Quataert and Schekochihin2017; Kunz et al. Reference Kunz, Squire, Schekochihin and Quataert2020). Therefore, further studies are necessary to assess exactly what key physics of a weakly collisional plasma leads to the measured effective collisionality, since most astrophysical plasmas, being multi-scale and turbulent, will support effective collision mechanisms (Zhuravleva et al. Reference Zhuravleva, Churazov, Schekochihin, Allen, Vikhlinin and Werner2019). The measurements presented here provide constraints to be satisfied by theories of effective collision processes.
Acknowledgements
We would like to thank S.D. Bale for pointing us to the dataset and L.B. Wilson III for helping us understand the dataset.
Editor Thierry Passot thanks the referees for their advice in evaluating this article.
Funding
J.T.C. was supported by a QMUL Principal Studentship. C.H.K.C. was supported by UKRI Future Leaders Fellowship MR/W007657/1, STFC Ernest Rutherford Fellowship ST/N003748/2, and STFC Consolidated Grant ST/T00018X/1. Support for J.S. was provided by Rutherford Discovery Fellowship RDF-U001804, which is managed through the Royal Society Te Apārangi.
Declaration of interests
The authors report no conflict of interest.
Appendix A. Collision length and time scales
Following the unit convention of Huba (Reference Huba1983), the Spitzer–Härm proton–proton collision frequency (Spitzer Reference Spitzer1962) for a proton–electron plasma with $T_{\mathrm {p}} \leqslant T_{\mathrm {e}}$, where $T_{\mathrm {p}} \, (T_{\mathrm {e}})$ is the proton (electron) temperature, is written
where $n_{\mathrm {p}} \ (\mathrm {cm}^{-3})$ is the proton number density, $k_{\mathrm {B}} \, T_{\mathrm {p}}$ is in eV and the Coulomb logarithm is $\lambda$. The Coulomb logarithm for proton–proton collisions
The dataset described in § 3 provides the following averages:
where the proton thermal speed is $v_{\mathrm {th}}^{\mathrm {p}}$. With these measurements collision scales can be calculated
where $\lambda _{\mathrm {mfp}}^{\mathrm {SH}}$ is the Spitzer–Härm proton–proton mean free path.
Appendix B. Linear collisional-kinetic MHD
The kinetic magnetohydrodynamic (KMHD) equations are found in Kulsrud et al. (Reference Kulsrud, Sagdeev and Rosenbluth1980) where a BGK collision operator can be added (Snyder et al. Reference Snyder, Hammett and Dorland1997) to study wave modes of arbitrary collision frequency (Sharma et al. Reference Sharma, Hammett and Quataert2003). The BGK collision operator is defined as
where $F_s$ is the equilibrium distribution function, assumed to be a Maxwellian, parameterised by the density ($n_s$) and temperature ($T_s$) to conserve particle number and energy. The BGK collision operator approximates any process that restores the distribution function ($f_s$) to the equilibrium ($F_s$) at a rate of $\nu ^{s}$.
The so-called drift kinetic equation is derived by transforming into the bulk velocity frame and assuming the distribution function is gyrotropic so the background electric can be ignored (Kulsrud et al. Reference Kulsrud, Sagdeev and Rosenbluth1980). The equations are linearised on a static background and Fourier transformed to produce
where $\omega$ is the complex frequency, $\tilde {u}_{\perp } \ (\tilde {u}_{\parallel })$ is the perpendicular (parallel) bulk velocity, $k_{\perp } \ (k_{\parallel })$ is the perpendicular (parallel) wavenumber, $\beta$ is the proton beta, $\tilde {b}_{\perp } \ (\tilde {b}_{\parallel })$ is the perpendicular (parallel) magnetic field, $\tilde {p}^{ \mathrm {s}}_{\perp } (\tilde {p}^{\mathrm {s}}_{\parallel })$ is the perpendicular (parallel) species ‘s’ pressure and $\tilde {n}$ is the density (quasi-neutrality). The tilde denotes the Fourier amplitude. The equations are closed by taking density and pressure moments of the linear Fourier analysed drift kinetic equation, which produces
with the definitions
and $\mathcal {R}^{s} = 1 + \zeta ^{s} \mathcal {Z}^{s}$. The plasma dispersion function is $\mathcal {Z}^{s}$ (Fried & Conte Reference Fried and Conte2015). Throughout the text $\beta$ is the proton beta and $\beta _{\mathrm {e}}$ is the electron beta. The kinetic equation for the electrons must be solved since the electron pressure appears in the momentum equation. The parameters relevant to the electrons are $T_{\mathrm {p}} / T_{\mathrm {e}}$ which is set to 1 and the effective mean free path species ratio $\lambda _{\mathrm {mfp}}^{\mathrm {eff}} / \lambda _{\mathrm {mfp,electrons}}$ which is set to 1, where the mean free path of the electrons is $\lambda _{\mathrm {mfp,electrons}}$. These equations form a linear system of equations that can be solved numerically and provide all the numerical solutions of the Letter.
Appendix C. Weighted geometric statistics
This appendix follows Norris (Reference Norris1940); Kendall & Stuart (Reference Kendall and Stuart1977). If the observations $x_i$ have unnormalised weights $w_i$, the definitions follow of the weighted geometric mean
weighted geometric covariance matrix
and the two sigma standard deviation
where $\mathrm {CI}^{x}$ is the weighted geometric confidence interval. The statistics detailed here are used to calculate the main results of the Letter, which are reported in table 1.
Figure 5 shows that the statistics of $\alpha, k_{\mathrm {iso}}, \lambda _{\mathrm {mfp}}^{\mathrm {eff}}$ for $C_{\perp }, C_{\parallel }, A_{\perp }, A_{\parallel }$ are similar so that combined statistics are reported in table 1. Figure 5(b,d,f) shows the normalised weighted geometric covariance between the model parameters
These terms are small.
Appendix D. Wavenumber model
The model introduced here $k_{\parallel } \sim k_{\perp }^{\alpha }$ is generalised from the critical balance model of Alfvénic turbulence (Goldreich & Sridhar Reference Goldreich and Sridhar1995). In this Letter it is used to model the compressive wave propagation angle $\theta _{\hat {{\boldsymbol {b}}}, \hat {{\boldsymbol {k}}}}$. To ensure the isotropic wavenumber ($k_{\mathrm {iso}} = k$ when $k_{\perp } = k_{\parallel }$) is defined correctly
which leads to
which appears as (4.1). Here, $k$ depends on $\theta _{\hat {{\boldsymbol {b}}}, \hat {{\boldsymbol {k}}}}$ parametrised by $\alpha \in [0,1), \ k_{\mathrm {iso}}$. Equation (D2) can be inverted on $\theta _{\hat {{\boldsymbol {b}}}, \hat {{\boldsymbol {k}}}} \in [0,90^{\circ })$ for $k$.
The wavenumber model is also used in the derivation of (4.4) of the Letter. Just above (4.4) the relation
is argued to define $k_{\parallel }^{\mathrm {trans}}$, which can be compared with measurements with the full wavenumber $k^{\mathrm {trans}}$. Using the model, (D2), to write
solving for $\theta _{\hat {{\boldsymbol {b}}}, \hat {{\boldsymbol {k}}}}^{\mathrm {trans}}$
then $k^{\mathrm {trans}}$ can be written
Using trigonometric identities the transition wavenumber is,
This appears as (4.4) of the letter.