Non-Thermal Solar Wind Electron Velocity Distribution Function

The quiet-time solar wind electrons feature non-thermal characteristics when viewed from the perspective of their velocity distribution functions. They typically have an appearance of being composed of a denser thermal “core” population plus a tenuous energetic “halo” population. At first, such a feature was empirically fitted with the kappa velocity space distribution function, but ever since the ground-breaking work by Tsallis, the space physics community has embraced the potential implication of the kappa distribution as reflecting the non-extensive nature of the space plasma. From the viewpoint of microscopic plasma theory, the formation of the non-thermal electron velocity distribution function can be interpreted in terms of the plasma being in a state of turbulent quasi-equilibrium. Such a finding brings forth the possible existence of a profound inter-relationship between the non-extensive statistical state and the turbulent quasi-equilibrium state. The present paper further develops the idea of solar wind electrons being in the turbulent equilibrium, but, unlike the previous model, which involves the electrostatic turbulence near the plasma oscillation frequency (i.e., Langmuir turbulence), the present paper considers the impact of transverse electromagnetic turbulence, particularly, the turbulence in the whistler-mode frequency range. It is found that the coupling of spontaneously emitted thermal fluctuations and the background turbulence leads to the formation of a non-thermal electron velocity distribution function of the type observed in the solar wind during quiet times. This demonstrates that the whistler-range turbulence represents an alternative mechanism for producing the kappa-like non-thermal distribution, especially close to the Sun and in the near-Earth space environment.


Introduction
In situ measurements of charged particles in the near-Earth space environment by artificial satellite became possible during the decade of the 1960s.It was realized then that the velocity space distributions of charged particles that make up the space plasma deviate from the expected Maxwell-Boltzmann-Gauss statistics; instead, the observed distributions typically feature a suprathermal (or non-thermal) component with inverse power-law "tail" characteristics for the suprathermal velocity regime, f ∝ v −γ for v ≫ α, where v represents the particle speed, f is the charged particle velocity distribution function, γ is the inverse power-law index, and α denotes the thermal speed [1][2][3].Recent inner heliospheric missions, the Parker Solar Probe and Solar Orbiter, further confirm that such a non-thermal feature persists even for heliospheric environments much closer to the Sun [4][5][6].The physical origin of such a feature was not understood then.Instead, Olbert and Vasyliunas [7][8][9][10] introduced an empirical model known as the kappa distribution, to fit the observation.Here, α = (2k B T/m) 1/2 is the Maxwellian thermal speed, meaning that α is the thermal speed had f (v) been given by the Maxwell-Boltzmann distribution.k B = 1.3806503 × 10 −23 m 2 kg s −2 K −1 is the Boltzmann constant, which can replaced by unity if we adopt the unit of eV for thermal energy.That is, if the temperature T is expressed in eV instead of Kelvins (K), then we may take k B = 1.Hereafter, we shall adopt such a convention.The mass of the charged particles is denoted by m.The free parameter κ determines the degree to which the observed distribution deviates from the Maxwellian-Boltzmann (MB or thermal) distribution in that if κ → ∞, then the model reduces to the thermal distribution, f MB (v) ∝ exp −v 2 /α 2 , while for v ≫ α, the kappa model depicts an inverse power-law velocity distribution, κ+1) .It is to be noted that, regardless of the value of the κ index, the kappa distribution approximates the MB distribution for v ≤ α-to be more precise, for a low v, the kappa distribution approaches the MB distribution with a sightly lower thermal speed, by a factor of κ/(κ + 1).That is, the kappa model naturally encompasses the quasi-Maxwellian feature in the "core" part of the velocity distribution characterized by v ≤ α and the inverse power-law tail portion of the distribution for the suprathermal regime, v ≫ α.
A sample non-thermal charged particle velocity distribution function in space is shown in Figure 1.Specifically, Figure 1 plots the typical electron velocity distribution function measured in the near-Earth space environment during quiet-time conditions.Figure 1 is a reproduction of Figure 4 of Ref. [11], and it shows two typical electron velocity distribution functions (eVDFs) in the solar wind at 1 au (astronomical unit) measured by the Wind/3DP electrostatic analyzers EESA-L and EESA-H.The left panels (a) and (c) show an eVDF in the slow solar wind (at 1995-06-19/00:06:38), and the right panels (b) and (d) show an eVDF in the fast solar wind (at 1995-06-19/23:13:59).The top panels (a) and (b) show cuts through the eVDF in one of the two directions perpendicular to the local magnetic field B: the diamonds are data points from EESA-L and the asterisks are data points from EESA-H.The dotted lines represent the one-count level for EESA-L and EESA-H.The blue dashed line in Figure 1a,b represents the sum of Maxwellian and kappa distributions (indicated in blue).The red line represents the fit to the measured perpendicular eVDF cut; the resulting fit parameters are indicated in red.The bottom panels (c) and (d) show cuts through the eVDF in the direction parallel to B. The perpendicular fit is shown in red, and the perpendicular fit parameters are used to initialize the parallel eVDF fit.The blue dashed line in Figure 1a,b represents the sum of Maxwellian and kappa distributions calculated using independent measurements of the core and halo densities and temperatures obtained from the fit of the spectrum of quasi-thermal fluctuations around the electron plasma frequency measured by the Wind/Waves electric field antennas.This "quasi-thermal noise" (QTN) technique is immune to spacecraft potential and therefore offers an independent and highly accurate measure of the core electron density and temperature, which are used as a reference to initiate the nonlinear least squares fitting of the measured VDF, resulting in the red curve fit, whose fit parameters are indicated in red as well.For more details, see Ref. [11] and Figure 4 therewith, including the accompanying description.
It is well known that the MB distribution corresponds to the maximum entropic (or the most probable) state as defined through the textbook Boltzmann-Gibbs (BG) definition for the entropy [12][13][14], namely, where dx is the spatial integration normalized to the total volume, dx → V −1 dx, and f (v) is the velocity distribution function.The Boltzmann-Gibbs (BG) entropy, which is additive and extensive, applies to an ideal gas or systems dictated by short-range interactions.The suitability of BG entropy for systems interacting through long-range forces, such as the plasma or gravitational systems, has been questioned since the inception of the BG entropy in the first place [15][16][17].The additive property relates to the BG entropy of a total system being equal to the entropies of subsystems.The extensivity means that the entropy is proportional to the total number of particles.The non-additive/non-extensive entropy, which presumably may be applicable to systems governed by long-range forces, violates these properties [18].The mathematical form of non-extensive entropy, which became well-known thanks to the work by Tsallis [19], was apparently independently discovered several times over, as entry 107 in Ref. [18], p. 347, describes.Specifically, it is mentioned there that several authors have independently rediscovered the form of entropy The list includes J. Havrda and F. Charvat, Kybernetika   For more details, see Ref. [11].
These earlier works notwithstanding, it is Tsallis's model [19] that is most well known, and it has triggered an explosive growth of recent interest in the topic of non-extensive thermostatics, in the space plasma context as well as in other applications [10,20,21].The celebrated Tsallis entropy in continuum form is defined by and the velocity distribution that corresponds to the maximum entropic (or the most probable) state is given by Upon identifying κ = 1/(1 − q) or alternatively κ = q/(1 − q), one finds that the solution reduces to either f Strictly speaking, neither is exactly identical to the kappa distribution since f κ is defined with a mixed κ and κ + 1-see Equation (1).Nonetheless, this convergence of Tsallis's nonextensive entropic principle and the kappa model has led to the space physics community embracing the notion that the space plasma may be in a state of non-extensive statistical quasi-equilibrium [10,20,[22][23][24].
From the microscopic plasma physics, it is known that the electron kappa distribution can be regarded as an end product of the weak electrostatic Langmuir turbulence [25,26].The initial findings involved a numerical study of a gentle weak electron beam-plasma (or bumpon-tail) instability and subsequent saturation of the Langmuir turbulence.It was found that the quasi-steady state of the Langmuir turbulence is characterized by the formation of a nonthermal, kappa-like velocity distribution function.Subsequently, more rigorous theoretical analysis revealed that the kappa distribution belongs to a family of unique solutions that characterize a steady-state electrostatic plasma turbulence [27,28].This finding implies that a profound inter-relationship may exist between the non-extensive statistical state and the turbulent quasi-equilibrium, but the precise mathematical formulation to establish such a connection does not yet exist at present.
The findings in Refs.[27,28] directly relate to the solar wind electrons [11,29], which can be interpreted as velocity distribution functions made of multiple subcomponents.The primary component is the quasi-Maxwellian core population (∼90-95% of the density, with ∼10 eV).The hotter and tenuous halo electron population can be distinguished from the core population by its distinct velocity profile, which can be modeled by an inverse power law.Other distinct populations can also be identified.For high-speed solar wind streams, a highly field-aligned strahl component can be separately classified from the halo electrons by their narrow pitch-angle distribution.The halo/strahl density is about ∼5-10% of the total density with a ∼50 eV energy range.Also, the highly energetic superhalo electrons (with a typical energy in the ∼2 keV range but extending up to 100 keV), which are observed in nearly all solar wind conditions, including the inner heliosphere [4,5] with a nearly invariant velocity power-law index, are a distinct component [29][30][31].The core, halo/strahl, and superhalo electron populations are sometimes associated with their respective slight temperature anisotropy, although the superhalo is almost completely isotropic, and relative drifts between the core and halo can also be detected.In the present discussion, however, we idealize the situation by considering that the velocity distribution is isotropic and without any net drifts.
In the present paper, we will first briefly overview the previous weak turbulence theory of electron kappa distribution [27], but, thereafter, we will discuss a new development, which involves the whistler-mode fluctuations and turbulence.For the near-Earth space plasma environment as well as for the inner heliosphere close to the Sun, the effects of wave-particle resonant interaction that involves the whistler-mode waves, instability, and fluctuations on the electrons are important [6,[32][33][34][35].As such, we consider the consequence of the electrons undergoing wave-particle resonant interactions with the background turbulence in the whistler-mode frequency range in the present paper.As will be shown, the impact of such interactions is none other than the formation of a non-thermal velocity distribution function for the electrons, which is not necessarily the kappa distribution but rather a more general one that must be generated by a numerical indefinite velocity integration.However, in the theoretical formalism of the present paper, it turns out that thermal fluctuations play an important role.A finite-temperature plasma constantly spontaneously emits and reabsorbs electromagnetic fluctuations-the fluctuation-dissipation theorem.A correct self-consistent theory of steady-state plasma particle velocity distribution based upon the steady-state Fokker-Planck particle kinetic equation thus requires the computation of thermal fluctuations.We thus begin the discourse by considering the thermal fluctuations emitted by the core electrons and the modification of the fluctuation spectrum by the presence of background turbulence.

Thermal Fluctuations Emitted by Maxwellian Core Electrons in the Background of Solar Wind Turbulence
In this section, we discuss the quasi-steady-state spectrum of the electrostatic and electromagnetic fluctuations in the background of solar wind turbulence.We assume that the thermal fluctuations are spontaneously emitted and reabsorbed predominantly by the Maxwellian core electrons.The background large-amplitude turbulence is assumed to be of the transverse electromagnetic type, with its characteristic frequency that encompasses the whistler-mode frequency range.The combined fluctuations and turbulence spectra determine the quasi-steady-state velocity distribution function for the solar wind halo electrons.As already discussed, the solar wind electrons are observed to be made of several distinct components, but the simplest description pertains to the two-component model, in which these electrons comprise dense Maxwellian core electrons and a tenuous but energetic halo electron population.It turns out that the halo electrons immersed in the field of thermal fluctuations alone will be organized in velocity space into a Maxwellian distribution.Thus, in this case, there will be no distinction between the core and halo so that both species will form one continuous thermal population.However, if there exist turbulent wave spectra for the whistler mode, then as the electrons interact with these combined spontaneously generated fluctuations and turbulence, they will organize into a non-thermal velocity distribution function, which manifests a clear demarcation between the core population and a tail component.The spontaneous emission is important because these background fluctuations provide the basis upon which non-thermal distribution can be built.
The electrostatic component of the spontaneous emission [36,37] is the well-known quasi-thermal noise [38], but the solar wind core electrons should also emit electromagnetic emissions as well, although a clear identification of such a transverse quasi-thermal noise is difficult because it will be partially occulted by the background turbulence.However, with improved future detection techniques, identifying the transverse quasi-thermal noise may become possible.Although we expect the electric and magnetic fields associated with the whistler-mode fluctuations to partially overlap with the frequency range of the background solar wind turbulence, the spectrum should extend to slightly higher frequencies so that with sufficiently sensitive instruments, the identification could be possible.Even with today's technology, if one analyzes the data with sufficient accuracy for the high-frequency end of the spectrum, one should be able to discern the characteristic signature associated with the whistler-mode thermal spectrum.Regardless, from a theoretical perspective, consideration of the emission of electromagnetic fluctuations in the whistler mode is important.In the presence of the combined background spectrum of Langmuir and whistlermode fluctuations as well as the whistler wave turbulence, it will be shown that the solar wind electrons naturally form a non-thermal velocity distribution function of the type observed in space, but it is not necessarily the kappa model in the analytic sense.Rather, the model distribution will be obtained by a numerical indefinite velocity integration.The first step in the present discussion is to consider the spectrum of electrostatic and electromagnetic fluctuations emitted by the thermal core electrons.In Ref. [39], the formulae for these fluctuation spectra are derived.For electromagnetic fluctuations propagating in a parallel direction with respect to the ambient magnetic field vector, the transverse electric and magnetic field spectra are designated as ⟨δE 2 ⊥ ⟩ k,ω and ⟨δB 2 ⊥ ⟩ k,ω , while for electrostatic fluctuations characterized by propagation parallel to the ambient magnetic field, the electric field spectrum is denoted by ⟨δE 2 ∥ ⟩ k,ω .These are given by [39][40][41].
where e is the unit electric charge; ω pe = (4πn/m e ) 1/2 e is the plasma frequency, with n and m e being the ambient density and electron mass, respectively; Ω e = eB/(m e c) is the electron cyclotron frequency, with B and c being the ambient magnetic field intensity and the speed of light, respectively; and k 2 m = Ω 2 e /α 2 e is the maximum perpendicular wave length, which results from the integration over the perpendicular wave number, with α e = (2T/m e ) 1/2 being the electron thermal speed [40].Here, f = f (v ⊥ , v ∥ ) represents the electron velocity distribution function (normalized to unity, dv f = 1), with v ⊥ and v ∥ denoting the velocity component perpendicular and parallel to the ambient magnetic field.The angular frequency and the parallel wave number are defined by ω and k, respectively.
For the Maxwellian thermal velocity distribution function, these are given as shown below: where is the plasma dispersion function with the prime indicating the derivative with respect to the argument.
Figure 2 plots the electrostatic and electromagnetic spectra, ⟨δE 2 ∥ ⟩ k,ω and ⟨δE 2 ⊥ ⟩ k,ω , respectively, computed from the theoretical formulae (4), versus ck/ω pe (horizontal axis) and ω/Ω e (vertical axis).The color scale is relative in that the maximum value for each panel is represented by red and the minimum intensity is plotted as a blue backdrop.The left-hand top and bottom panels correspond to the electrostatic and electromagnetic fluctuation spectra, respectively.The input parameters are ω pe /Ω e = 5 and β e = 1, where β e = 8πnT e /B 2 is the electron beta (ratio of electron thermal energy to the magnetic field energy).In order to verify that the theoretical formalism ( 4) is indeed reliable, we have also carried out a one-dimensional particle-in-cell (PIC) simulation.We have used a simulation box of L x = 512 c/Ω pe and n x = 4096 grid points, with 2000 particles per grid per species.The time step used was ∆t = 0.01 /ω pe , and the simulation ran until t = 2621.44/ω pe .The ratio of plasma frequency to electron gyro frequency was ω pe /Ω e = 5.This ratio is somewhat lower than the actual value typical of the solar wind at 1 au, which is close to ω pe /Ω e ∼ O(10)-O(10 2 ), but for the sake of illustration we have chosen a relatively low value of ω pe /Ω e .Otherwise, the spectral peak at ω ∼ ω pe associated with the electrostatic thermal fluctuations (upper panels) and the spectral characteristics associated with the transverse-mode fluctuations around the electron cyclotron frequency and below would have been separated by a wide gap, which would have made visual inspection quite challenging.Also, if the separation between the two frequencies is too high, it becomes very challenging for the simulations, too, because we need to resolve both time scales.Other parameters were electron and proton betas, which were taken to be β e = β p = 0.1.These choices are not atypical of the solar wind conditions at 1 au.The simulated electrostatic and electromagnetic fluctuation spectra are plotted in the top and bottom right-hand panels, respectively.As the readers may appreciate, the theoretical plots compare very well with the simulated spectra, which indicates that the theoretical method is a reliable tool for describing the spontaneously emitted thermal spectra in magnetized plasmas accurately.The electrostatic fluctuation spectrum is enhanced along the Langmuir wave dispersion curve but broadens in frequency somewhat for shorter wavelengths.In the simulated spectrum, the enhanced fluctuation along the Langmuir wave dispersion curve is broader than that of the theoretical spectrum, but, otherwise, the overall agreement is excellent.For the electromagnetic spectrum, it is seen that the fluctuation spectrum is enhanced along the whistler-mode dispersion curve, but the triangular (or conical) emission pattern that converges to the electron cyclotron frequency, ω = Ω e at the k → 0 limit, is also prominent in both the theoretical emission spectrum and the simulated spectrum.Such a feature is associated with the virtual (or higher-order) modes, that is, heavily damped solutions of the linear dispersion relation [42][43][44].Both the theoretical and simulated spectra accurately reproduce the emission characteristics associated with such modes.Note, however, that the simulation does not completely demonstrate the intensification of the higher-order mode as k approaches a zero value.This is owing to the limited resolution in the simulation spectrum.As will be discussed, this limitation further affects the k-integrated wave spectra for the electric and magnetic fields.
Shown in Figure 3 are wave number-integrated ( dk • • • ) spectra.The left-hand panel shows the k-integrated magnetic and transverse electric field fluctuation spectra dk⟨δB 2 ⟩ k,ω (red) and dk⟨δE 2 ⊥ ⟩ k,ω (blue) that were computed from theory, plotted against ω/Ω e .The right-hand panel displays the same spectra constructed from the PIC simulation result and integrated over the wave numbers.Both the theoretical and simulated spectra exhibit the behavior of increasing intensities, for both magnetic and electric spectra, over an increasing frequency, up to ω ∼ 0.5Ω e or so.However, some differences in the behavior are also evident.For instance, in the theoretical integrated spectra, both the electric and magnetic field intensify around the electron cyclotron frequency, ω ∼ Ω e , but the simulated spectra do not exhibit such a behavior.Clearly, the peak at ω/Ω e = 1 in the theoretical spectrum is associated with the contribution from the higher-order mode.In the simulated spectrum, the higher-order mode for the low k regime is not as clearly enhanced, which explains the absence of such a peak.This is due to the limited resolution in the simulation.Such an increasing behavior as a function of frequency for the fluctuation spectra in the lowfrequency regime is a characteristic of the plasma, and it is the baseline spectral behavior associated with the thermal motion of plasma particles.It is interesting to note that in many PIC simulations of low-frequency turbulence, such an increasing intensity can be seen at the high end of the simulation spectrum.In a typical kinetic simulation of the low-frequency turbulence, the MHD-like regime corresponding to ω 2 ≪ Ω 2 p ≪ Ω 2 e , where Ω i = eB/(m p c) is the proton cyclotron frequency, is characterized by a Kolmogorov type of inverse power-law spectrum, k −5/3 [45][46][47], but as k increases, in some cases, the intensity actually rises again [47,48].In the literature, such a behavior is not clearly explained nor understood.However, it is entirely possible that the simulation system is automatically generating the background thermal spectrum.As confirmed by Figures 2 and 3, the theoretical description of thermal fluctuations is consistent with the simulation result.Thus, we now focus on the analytical approach.Furthermore, henceforth, we are interested in the fluctuations associated with the eigen modes.For the electrostatic fluctuation, we are concerned with the spectral wave intensity along the Langmuir mode dispersion relation, ω = ω L (k), where e /(4ω 2 pe ) .Likewise, for the electromagnetic fluctuations, we pay attention to the whistler-mode dispersion relation, ω = ω W (k), where ω W = Ω e c 2 k 2 /(ω 2 pe + c 2 k 2 ).Then, by expanding the denominators by ϵ ∥ ≈ (ω where For more details regarding the derivation of this result, see [39][40][41].
In the solar wind, there exists a permanent low-frequency turbulence of a solar origin.Such turbulence is commonly believed to be generated on the surface of the Sun through various mechanisms, including the solar surface convection and small reconnection near the lower corona, and convected to outer space [49].The solar wind turbulence for a low-MHD frequency regime is hydromagnetic in nature and is characterized by a Kolmogorov-like inertial range spectrum but with a spectral break in the kinetic regime.That is, for the frequency range above the nominal proton cyclotron frequency and below the electron cyclotron frequency, Ω p < ω < |Ω e |, the turbulence exhibits a spectral break.Such a frequency range can be characterized as the whistler turbulence range.For an even higher frequency ω > |Ω e |, another spectral break is present.We may model such a multi-scale spectral behavior by adopting an analytical model first suggested by von Kármán [50] and generalizing to reflect the multiple spectral breaks, where l ≫ c/ω pi .Here, we explicitly extracted out the factor k 2 m since this is related to the integration over k ⊥ [40].The solar wind turbulence spectrum appears to behave as ω −5/3 in the frequency range corresponding to the MHD regime.If we make use of the Taylor hypothesis [51], then ω can be trivially replaced by k, but in the kinetic regime, beyond the ion skin depth, c/ω pi or shorter, and much more so for the electron skin depth, c/ω pe or shorter, the Taylor hypothesis may not be valid.Moreover, since we are interested in the parallel wave vector and the turbulence intensity integrated over k ⊥ , the inverse power-law index α may not be the same as that of the Kolmogorov value, namely α = 5/6.Nevertheless, we may model the MHD regime by the Kolmogorov type of spectrum.In any event, the model spectrum (8) describes a finite and maximum turbulence level at k = 0, and for 0 < k < ω pi /c it describes the k −2α behavior.For the wave number regime corresponding to ω pi /c < k < ω pe /c, the model depicts a k −2β behavior.For k > ω pe /c, the spectrum behaves as k −2γ .We illustrate this by choosing α = 5/6, β = 1.2, and γ = 2, which are admittedly arbitrary.We also choose the MHD scale factor l = 10 2 (m p /m e ), which is again arbitrary.
In Figure 4, we demonstrate the influence of the whistler-mode fluctuation spectrum on the background turbulence spectrum by considering the superposition of the model turbulence spectrum and the whistler-mode fluctuation spectrum, I turb (k) + I W (k), where the whistler-mode fluctuation spectrum I W (k) is defined in Equation ( 6) and the model turbulence spectrum I turb (k) is given by Equation (7). Figure 4 plots the spectral factor that defines the combined spectrum, namely, where α = 5/6, β = 1.2, γ = 2, and l = 10 2 (m p /m e ), as already noted above.The first term on the right-hand, q 2 /(1 + q 2 ) 3 , denotes the spontaneously emitted whistler-mode thermal fluctuation spectrum.The second term on the right-hand side is the model spectrum with multiple spectral breaks.In Figure 4, the dashed magenta-colored curve represents the spontaneously emitted whistler-mode fluctuation spectrum, q 2 /(1 + q 2 ) 3 .The dashed black curves represent the background turbulence spectrum without the influence of the fluctuation, for two cases of R = 10 6 and R = 10 7 .The total spectral factor for the two cases is plotted with thick blue (R = 10 6 ) and red (R = 10 7 ) curves.It is evident that the model turbulence spectrum (8) depicts a flat spectrum for a q → 0 regime; a Kolmogorov-type of spectrum in the "MHD" regime, (Ml) −1/2 < q < M −1/2 ; a slightly steeper spectrum of k −2.4 in the "kinetic proton" regime, M −1/2 < q < 1; and a yet steeper spectrum, k −4 , in the whistler turbulence regime, q > 1.It is in this wave number regime where the presence of the thermal fluctuation spectrum should be discernible.Specifically, in the case of a relatively low turbulence level, as indicated by R = 10 6 , we expect that the actual solar wind turbulence should reveal the presence of the fluctuation.However, for higher turbulence levels (as denoted by R = 10 7 ), the intensity of fluctuation will be partially hidden so that a clear identification might not be so straightforward.Model spectral factor, S(q), versus q in logarithmic horizontal and vertical scales.The magenta dashes represent the spontaneous whistler-mode fluctuation, while the black dashes denote the background turbulence for two cases of R = 10 6 and R = 10 7 .The combined spectra are plotted with thick blue (R = 10 6 ) and red (R = 10 7 ) curves.For a low turbulence level (R = 10 6 ), the presence of fluctuation should be more evident, but if the turbulence level is high (R = 10 7 ), then the fluctuation will be partially hidden.
As an example of actual solar wind turbulence spectra measured in the near-Earth environment, we reproduce a figure taken from Ref. [52].The result is Figure 5, which is constructed from the measurements made by Cluster spacecraft.The location of the spacecraft is at 1 au during a quiet-time condition on 30 January 2003.The detailed discussion of the instrumentation and data analysis method can be found in Ref. [52], but the main focus of the present paper is bring the readers' attention to the spectral flattening behavior for the high-frequency end of the turbulence spectra, especially for the electric field.According to the theory-see Equation ( 6)-and the model spectrum shown in Figure 4, the spontaneously emitted thermal fluctuations should affect the high end of the solar wind turbulence spectra, especially if the turbulence level is sufficiently low.Admittedly, just what exactly it means by "sufficiently low" is not entirely clear, and further study is called for.Nevertheless, the identification of the spontaneous quasi-thermal whistler-mode fluctuations based on observation could be an intriguing and innovative research topic.In any case, Figure 5 displays the Kolmogorov-like spectrum in the lowfrequency band while also showing a spectral break at frequency f b , which represents the "break" frequency for the transition of one spectral slope to another.This frequency could be associated with the kinetic proton effects.The whistler-mode thermal noise, however, is supposed to be associated with the electron kinetic effects, which are believed to be related to a much higher frequency.However, before one could reach such a frequency, the instrument noise floor would contaminate the data, so it is very challenging to delineate the noise effects versus the baseline thermal noise.Specifically, a key element that should be accounted for before one can definitely extract the theoretical signature, i.e., the whistler-mode thermal fluctuation, from the data is for the model turbulence spectrum of the type shown in Figure 4 to be translated into the spacecraft frame frequency using the appropriate solar wind speed and electron inertial length, as well as to properly scale the normalized amplitude S(q) to physical units.This includes translating the "R" parameter into actual units.The flattening of the E-field spectrum shown in Figure 5 could be entirely due to the instrumental artifact.In spite of this, however, the thermal noise associated with the whistler-mode fluctuations could partly contribute to the observed flattening of the spectrum, if not for this particular event, then at least for some other events.As will be shown in the next section, the combination of the quasi-thermal whistler noise and the background turbulence can account for the observed non-thermal electron velocity distribution function.We thus proceed with the discussion of the theory for the formation of the electron velocity distribution function under the influence of background whistler-mode turbulence and the quasi-thermal noise spectrum, which contains both electrostatic Langmuir-type and whistlermode-type electromagnetic fluctuations.

Formation of Kappa Electron Distribution by Langmuir Turbulence
In this section, we briefly overview the previous theory of kappa electron distribution by Langmuir turbulence advanced by Yoon [27].The full discourse of this theory is quite complex and requires a detailed exposition of kinetic weak plasma turbulence theory [37,[53][54][55][56][57][58][59][60][61][62][63][64][65], but, in its essence, it boils down to the modification of the spontaneously emitted Langmuir fluctuations to reflect the influence of the steady-state weak Langmuir turbulence spectrum.It was shown by considering the balance of the nonlinear wave kinetic equation for Langmuir turbulence that, in the steady state, the electrostatic fluctuation spectrum should be modified to include the effects of turbulence in the following form: The modification factor k 2 L /k 2 leads to the kappa electron distribution function when this spectrum is inserted into the diffusion coefficient of the steady-state electron distribution function computed from the kinetic theory.
Reference [40] derives the Fokker-Planck electron kinetic equation with waves and fluctuations that have a wave vector lying in the parallel direction defined with the ambient magnetic field vector.We summarize the equation for the electron velocity distribution function f , where the right-hand side of the kinetic equation is expressed in a velocity-space spherical coordinate system, in which v = v 2 ⊥ + v 2 ∥ is the magnitude and µ = v ∥ /v is the cosine of the pitch angle.Under the assumption of primarily electrostatic interaction, the velocity space friction and diffusion coefficients are given by We assume steady state, ∂/∂t → 0, and isotropy, ∂ f /∂µ = 0.Then, we average over µ.Then, we obtain the steady-state solution for the electron velocity distribution function, We should note that this type of steady-state solution of the Fokker-Planck equation is found in the literature [66][67][68], so the basic concept is not new.Making use of the property , where I L (k) is given by Equation ( 9), then we have From this, we obtain the desired electron kappa velocity distribution function, if we identify In this version of the theory, the formation of a non-thermal (kappa) electron distribution is attributed to the Langmuir turbulence in the asymptotical steady state.According to this theory, no clear separation of the core and halo electrons is made, but, instead, both populations are treated as a single kappa distribution function with the low end of the velocity spectrum mimicking the Maxwellian thermal core, while the suprathermal highvelocity regime represents the inverse power-law tail population.The brief overview of this section is not new, and a full discourse can be found in Refs.[27,37,[53][54][55][56][57][58][59][60][61][62][63][64][65].In the remaining part of the present manuscript, we put forth a new model for which the role of whistler turbulence is emphasized.

Formation of Non-Thermal Electron Distribution by Combined Background Turbulence and Thermal Fluctuations
Section 2 discussed the thermal fluctuations spontaneously emitted by Maxwellian core electrons.We also discussed the effects of pre-existing solar wind turbulence and how the combined model may relate to the existing literature on low-frequency turbulence simulations.We also discussed how the effects of baseline quasi-thermal spontaneous emission fluctuations may impact the observations, although we noted that the unambiguous identification of the predicted spectral features associated with thermal fluctuations in the observation may depend on the level of turbulence.In this section, we proceed to discuss the combined impact of the quasi-thermal whistler-mode fluctuations and the background turbulence on the electron velocity distribution function.
In Section 3, the steady-state electron distribution function subject only to the Langmuir turbulence was discussed, and it was shown that the result is the kappa electron velocity distribution.For the kappa model, however, no distinction is made between the core and halo populations.Moreover, the spontaneously emitted transverse fluctuations in the whistler-mode frequency are ignored.Further, the presence of background solar wind turbulence is not taken into consideration either.As such, the kappa distribution and Langmuir turbulence problem may pertain to the outer heliosphere where the local ambient magnetic field strength is sufficiently low so that the whistler-mode frequency range effects can be ignored and the underlying plasma may be treated as essentially unmagnetized.
For the near-Earth space environment, however, the whistler-mode dynamics may be an integral part of wave-particle interaction with the electrons [32,33].Thus, this section discusses the formation of the non-thermal electron velocity distribution function in the presence of spontaneous thermal fluctuations in both the longitudinal Langmuir and transverse whistler modes and also under the influence of background turbulence [52,[69][70][71][72][73][74].As will be shown, under such a physical environment, the self-consistent steady-state solution for the electron velocity distribution function will be characterized by a distinct core and halo populations, which is consistent with observations.
The notion of the pre-existing whistler-mode turbulence affecting the solar wind electron dynamics, resulting in a non-thermal velocity distribution function, has been discussed in the literature.For instance, Refs.[75][76][77][78] carried out extensive and detailed numerical simulation based on the quasilinear velocity diffusion theory where the diffusion coefficient is computed from the model whistler wave turbulence spectrum.It was shown in these references that the resonant wave-particle interaction between the solar wind electrons and the pre-existing turbulence in the whistler-mode frequency range leads to the gradual formation of a non-thermal energetic tail.The present paper is similar in conceptual background in that we are also seeking to find a non-thermal electron velocity distribution function that is a result of resonant wave-particle interaction with the background whistler wave turbulence.However, the main distinction between the present approach and those of previous works is twofold.Firstly, unlike the previous works, which relate to the dynamical evolution of the velocity distribution function, we are concerned with the steady-state solution.This aspect leads to the second distinction.That is, in order to obtain the steady-state solution, the effects of spontaneous thermal fluctuation are essential.The effects of spontaneous thermal fluctuations and the related velocity friction effects are not considered in the dynamical theories of solar wind electron distribution function in the above references.For dynamical problems, the velocity friction, which is intimately related to the spontaneous thermal fluctuations, is indeed relatively unimportant, but for the theory of an asymptotic steady state, the balance of velocity friction coefficient A and velocity diffusion coefficient D is crucially important-see Equation (12).Reference [67], however, considered a steady-state solution of a magnetospheric electron velocity distribution in resonant wave-particle interaction with the background whistler-mode waves.Their solution is very similar in conceptual background and mathematical methodology to the present work, except that, in their approach, the velocity friction coefficient is replaced by the collisional drag coefficient.
We again start from the Fokker-Planck electron kinetic equation with waves and fluctuations where a wave vector lying in the parallel direction is defined with the ambient magnetic field vector, that is, Equation (10) or the steady-state solution (12), except that now the electrons are immersed in the bath of thermal fluctuations of both Langmuir and whistler types and also the background turbulence.As a result, the velocity friction and the diffusion coefficients now contain contributions from both longitudinal and transverse modes, For a steady state, the formal solution ( 12) is still applicable, with the coefficient A = and approximating the resonance delta function by δ(ω W − kvµ − Ω e ) ≈ δ(kvµ + Ω e ), we may proceed with the computation of generalized A and D coefficients.For the present purpose, we adopt the whistler-mode spectrum by considering the effects of thermal fluctuations and the background turbulence spectrum, that is, Equation ( 8), but in a simplified form.In particular, we are interested in the frequency range that is sufficiently higher than both the MHD scale and proton kinetic scale but is sufficiently below the electron cyclotron frequency.Thus, in such a low-frequency limit relative to the electron cyclotron frequency, ω 2 ≪ Ω 2 e , the whistler-mode fluctuations can be approximated by m T e /(4π)](ωΩ e /ω 2 pe ) and m T e /(4π)], which, upon making use of the low-frequency version of the dispersion relation, ω = Ω e (ck/ω pe ) 2 , can be written as This result, together with the Langmuir fluctuation spectrum, I L (k) ∼ k 2 m T e /(4π), will be inserted into the expressions for A and D.
We may also simplify the model of the turbulence given by Equation ( 8) in that we only focus on the portion of the background turbulence spectrum corresponding to the whistlermode range inverse power law, which we simplify by ∝ k −2ν .If we thus superpose this simplified background turbulence spectrum to the approximate form of the spontaneous emission spectrum corresponding to the low-frequency whistler-mode thermal emission, then we may adopt a simplified form of the combined spectrum, Here, k 2 W is an appropriate parameter for correct dimensionality, which can be adjusted.This parameter effectively dictates the level of background turbulence as well.In applying the above model, we reiterate that the model spectrum ( 18) is meant for the whistler-mode frequency range satisfying Ω 2 i < ω 2 < Ω 2 e .As such, we confine the width of wave numbers, c∆k/ω pe , roughly corresponding to the above frequency limitation.In an earlier attempt to incorporate the solar wind turbulence effects into the model whistler-mode spectrum, Ref. [41] adopted a model where the thermal fluctuation spectrum was modified to reflect the inverse power-law feature, namely, m T e /(4π)](Ω 2 e /ω 2 pe )(c 2 k 2 /ω 2 pe ) 1−β , where β is a control parameter that can be chosen as 0 in the case of purely spontaneous emission and as β = 1 + ν if we wish to model the overall spectral profile to behave as an inverse power law, I W (k) ∝ k −2ν .However, we now realize that the more proper way to model the combined spontaneously emitted quasi-thermal whistler-mode spectrum and the background pre-existing turbulence is the linear superposition (8), which we simplify as shown in Equation (18).Thus, in the present section, we take the total whistler-mode spectral intensity to possess the proportionality dictated by the functional relationship, For the Langmuir mode spectrum, however, we only consider the thermal fluctuation, which is distinct from the previous section.Recall that in Section 3 we included the steady-state Langmuir turbulence factor, (k L /k) 2 , in the Langmuir turbulence spectrum (9), which led to the electron kappa distribution.In the present section, we are concerned with an alternative theory of a non-thermal, generalized kappa distribution, which is based upon the notion of background whistler-mode turbulence.In short, the transverse and longitudinal electric field spectral intensities adopted in the present discussion are defined by Inserting this into the generalized coefficients ( 16 The integral 1 0 dµ (1 − µ 2 )µ −3 is formally divergent.To regularize the divergence, we introduce the lower limit, The other µ integral is evaluated in a straightforward manner: Making use of all this, we have Here, we have made use of (c/α e ) 2 (Ω e /ω pe ) 2 = B 2 0 /(8πn 0 T e ) = 1/β e .This is a threeparameter model distribution, with ν, a, and κ W being the adjustable parameters.If we consider the limit of a → ∞, then we have In this limit, the contribution from electrostatic Langmuir-mode fluctuation, that is, the term x ′ 2 in both the numerator and denominator within the integrand, is ignored.This limiting form can be termed the W-only distribution.However, if we take the limit of a → 0, then we simply have the Maxwell-Boltzmann (MB) distribution.In this limit, the contribution from the whistlermode related terms are ignored, and, thus, this limit can be termed the L-only distribution.Another interesting limit is when κ W → ∞.In this limit, the contribution from the background turbulence disappears, and the resulting distribution is that of the MD distribution again.The parameters a and κ W , in turn, are determined by µ min , ω pe /Ω e , β e , ν, and ck W /Ω e .The parameter ω 2 pe /Ω 2 e can be determined from the solar wind data.Also, β e is known from the data.The fitting parameters ck W /Ω e and ν relate to the spectral profile of the solar wind turbulence in the whistler-mode frequency range.Thus, these parameters can also be determined from observational properties.The truly free parameter is Λ, which is determined from the choice of µ min .Let us consider the resonance condition, kvµ + Ω e = 0 or µ = −Ω e /(kv).We are interested in the minimum value for µ.In the formal µ integral, this is taken to be µ min = 0, but this means either k → ∞ or v → ∞, neither of which are physical.For whistler turbulence and fluctuations in the low-frequency limit, we choose the maximum k by ck max /ω pe ∼ 1 or so.For the velocity v, we generally determine the maximum value to be sufficiently higher than the thermal speed, v max ≫ α e .From this, we may see that µ min ∼ −Ω e /(k max v max ), which, while small, can have a substantial range of freedom.If, for instance, we choose µ min ∼ 10 −6 or so, then we obtain Λ ∼ 10 10 .However, if we choose µ max to be approximately 10 −4 or so, then we have Λ ∼ 10 6 and so on and so forth.With this information, let us consider the ratio a/κ W .If we choose ν = 1, which implies k −2 spectral behavior associated with the whistler frequency range turbulence, then we have a Suppose we take (ω pe /Ω e ) 2 β e = 10 2 and ck W /Ω e = 0.2.Then, by choosing Λ = 2.25 × 10 6 or so, we arrive at a/κ W ∼ 3. If, however, (ω pe /Ω e ) 2 β e = 10 4 , then the choice of Λ = 2.25 × 10 10 leads to the similar value of a/κ W ∼ 3.In the solar wind, the ratio ω pe /Ω e can be quite high, ranging from O(10) to O(10 2 ).The electron beta value in the solar wind can range from β ∼ O(10 −2 ) to O(1) or so, hence the above two choices of parameters, (ω pe /Ω e ) 2 β e = 10 2 and (ω pe /Ω e ) 2 β e = 10 4 .The choice of ck W /Ω e = 0.2 relates to the turbulence property in that this number represents the maximum effective range of whistlermode turbulence in the wave number space.Since the low-frequency whistler mode is characterized by ck W /Ω e < 1, such a choice is eminently reasonable.The above estimation of the crucial dimensionless parameter a/κ W , of course, is a rough exercise, and more precise attempts should be made by surveying the 1 au data.However, in view of the uncertainty associated with the lower limit of the cosine of the pitch angle, µ min , we defer the more accurate attempts for future.With these considerations, we construct the asymptotic electron velocity distribution function (21) by performing a numerical indefinite integration over the dimensionless velocity x = v/α e .The result is displayed in Figure 6, where we display on the left-hand panel the case for a/κ W = 3 with a = 30 and κ W = 10.For all the examples, we restrict ourselves to ν = 1.As visual guides, we plot the so-called L-only and W-only limiting case distributions.We also plot the inverse power-law velocity slop v −6.5 .In the solar wind, such an asymptotic high-velocity tail distribution is often observed [29].On the right-hand panel, we show the velocity distribution by varying the parameter a, which ranges from a = 1 to 10 to 20 to 30.Other parameters are fixed: κ W = 10 and ν = 1. Figure 6 thus demonstrates that the combined effects of background turbulence and finite spontaneously emitted fluctuations are capable of producing the electron velocity distribution function that remarkably resembles the observed distribution in the solar wind.We should note, however, that there exists a certain degree of freedom in our choice of parameters a and κ W , in particular, their ratio, a/κ W , which turns out to be important for determining the shape of the velocity distribution function, as the right-hand panel of Figure 6 indicates.

Summary
The main purpose of the present paper was to put forth a first principle theory of the steady-state electron velocity distribution function with non-thermal characteristics, which resembles the quiet-time solar wind electron distribution detected in the near-Earth space environment.Unlike the previous model [27], which invoked the steady-state Langmuir turbulence and the accompanying kappa distribution, the present paper employed the combined influence of the background solar wind turbulence in the whistler frequency range, as well as the quasi-thermal electromagnetic and electrostatic fluctuations.The resulting model electron distribution function was given in terms of an exponential function of an indefinite velocity integral, which does not in general lend itself to further closed-form analytical manipulations-Equation ( 21)-but must, in general, be computed by numerical means.Under a reasonable set of assumptions and input conditions, we have found that the resulting numerical calculation leads to a velocity distribution function whose profile is reminiscent of the measured distribution in space.
The formal electron velocity distribution function in the steady-state was given by Equation (21), and this mathematical expression contains the effects of background fluctuations as well as the impact of the pre-existing whistler-mode turbulence.This contrasts to the formal solution (14), which reflects the influence of electrostatic thermal fluctuations and the enhanced Langmuir wave turbulence.It is possible to construct a more general formal distribution by combining the effects of both the electrostatic and electromagnetic thermal fluctuations on the electrons, as well as the influence of enhanced electrostatic (Langmuir) and electromagnetic (whistler) turbulence intensities.The result is Here, κ L = 2ω 2 pe /(k L α) 2 , as defined in Equation (15).With this form of the electron distribution function, we now summarily discuss the various limits.Suppose that we ignore the influence of background turbulence altogether.This amounts to taking κ W → ∞ and κ L → ∞, which leads to Thus, in the absence of turbulence, we obtain the MB distribution, which is as expected.
Ignoring the influence of electromagnetic whistler modes, both the thermal fluctuations and the enhanced pre-existing turbulence, is equivalent to taking the limit of a → 0, which leads to the generalized form of the L-only distribution-see Equation ( 23)-which also happens to be the same kappa distribution discussed in Equation ( 14) If, however, we are to ignore the influence of electrostatic modes, both the thermal fluctuations and enhanced Langmuir turbulence, then all we need to do is consider the limit of a ≫ 1, This form of limiting solution was referred to as the W-only distribution in Equation ( 22), but this stand-alone solution is not a realistic model since the electrostatic fluctuations cannot simply be ignored.Nevertheless, at least as a mathematical exercise, one can certainly entertain such a limit.In our final solution (21), we have considered the limit of κ L → ∞, while other parameters, a and κ W , are considered finite.The overall concept of charged particles maintaining a steady-state wave-particle interaction with steady-state turbulence and fluctuations is an example of a stationary state far from equilibrium [79].Such a state, in turn, may be considered as an example of the non-extensive statistical state [18,19].It is in this regard that the present paper is relevant to the Special Issue "Nonadditive Entropies and Nonextensive Statistical Mechanics".The fact that the space plasma, which is governed by a long-range electromagnetic force, frequently exhibits a kappa-like non-thermal distribution function is in a way not too surprising in that, thanks to Tsallis's pioneering work, we now have a rather insightful understanding that any system with a long-range interaction is likely to be governed by a non-extensive, non-additive statistical principle.The present paper, as with the related earlier work [27], provides the physical "mechanism" that leads to a concrete example of a kappa-like non-thermal phase space distribution function.

Figure 1 .
Figure 1.Reproduced from Figure 4 of Ref. [11]: Two typical electron velocity distribution functions (eVDFs) measured by EESA-L and EESA-H onboard Wind spacecraft at 1 au in the slow solar windpanels (a,c)-and in the fast solar wind-panels (b,d).The top panels (a,b) show cuts through the eVDF in one of the two directions perpendicular to the local magnetic field B. The bottom panels (c,d) show cuts through the eVDF in the direction parallel to B. Explanations for the different lines in the figure is given in the main text.For more details, see Ref. [11].

Figure 4 .
Figure 4.Model spectral factor, S(q), versus q in logarithmic horizontal and vertical scales.The magenta dashes represent the spontaneous whistler-mode fluctuation, while the black dashes denote the background turbulence for two cases of R = 10 6 and R = 10 7 .The combined spectra are plotted with thick blue (R = 10 6 ) and red (R = 10 7 ) curves.For a low turbulence level (R = 10 6 ), the presence of fluctuation should be more evident, but if the turbulence level is high (R = 10 7 ), then the fluctuation will be partially hidden.

1 − 1 1 − 1
dµA v and D = dµD vv now containing the influence of whistler-mode fluctuations as well as the background turbulence.Making use of Im [ω 2 ), we obtain the desired coefficients A and D, which are now given by

Figure 6 .
Figure 6.[Left] The three-parameter electron velocity distribution function f versus x = v/α e , for a = 30, κ W = 10, and ν = 1.[Right] The variation in the input parameter a, which ranges from a = 1 to 10 to 20 to 30.