Dual MaxwellianKappa modeling of the solar wind electrons: new clues on the temperature of Kappa populations
^{1} Centre for Mathematical Plasma Astrophysics, Celestijnenlaan 200B, 3001 Leuven, Belgium
email: mlazar@tp4.rub.de
^{2} Institut für Theoretische Physik, Lehrstuhl IV: Weltraum und Astrophysik, RuhrUniversität Bochum, 44780 Bochum, Germany
^{3} Royal Belgian Institute for Space Aeronomy, 3 av. Circulaire, 1180 Brussels, Belgium
^{4} Université Catholique de Louvain, Georges Lemaître Centre for Earth and Climate Research (TECLIM), Place Louis Pasteur 3, 1348 LouvainLaNeuve, Belgium
^{5} Theoretical Physics Research Group, Physics Department, Faculty of Science, Mansoura University, 35516 Mansoura, Egypt
^{6} Research Department of Complex Plamas, RuhrUniversität Bochum, 44780 Bochum, Germany
Received: 5 December 2016
Accepted: 31 March 2017
Context. Recent studies on Kappa distribution functions invoked in space plasma applications have emphasized two alternative approaches that may assume the temperature parameter either dependent or independent of the powerindex κ. Each of them can obtain justification in different scenarios involving Kappadistributed plasmas, but direct evidence supporting either of these two alternatives with measurements from laboratory or natural plasmas is not available yet.
Aims. This paper aims to provide more facts on this intriguing issue from direct fitting measurements of suprathermal electron populations present in the solar wind, as well as from their destabilizing effects predicted by these two alternative approaches.
Methods. Two fitting models are contrasted, namely, the global Kappa and the dual MaxwellianKappa models, which are currently invoked in theory and observations. The destabilizing effects of suprathermal electrons are characterized on the basis of a kinetic approach that accounts for the microscopic details of the velocity distribution.
Results. In order to be relevant, the model is chosen to accurately reproduce the observed distributions and this is achieved by a dual MaxwellianKappa distribution function. A statistical survey indicates a κdependent temperature of the suprathermal (halo) electrons for any heliocentric distance. Only for this approach are the instabilities driven by the temperature anisotropy found to be systematically stimulated by the abundance of suprathermal populations, thus lowering the values of κindex.
Key words: solar wind / plasmas / instabilities
© ESO, 2017
1. Introduction
Kappa modeling is a widely exploited technique for space plasma diagnosis and kinetic analysis (Pierrard & Lazar 2010; Livadiotis & McComas 2013). The Kappa (κ)distribution function is nearly Maxwellian at low energies and decreases as a powerlaw at higher energies, enabling a satisfactory overall (global) description of the velocity distributions of plasma particles measured in the solar wind (Vasyliunas 1968; Christon et al. 1989; Collier et al. 1996; Maksimovic et al. 1997). The powerindex κ quantifies the presence of suprathermal (nonMaxwellian) populations, which enhance the highenergy tails of the distributions. Alternatively, the same Kappa powerlaws have also been used to reproduce only the highenergy tails (also known as the halo component), while the core of the distribution is well fitted by a standard Maxwellian (Maksimovic et al. 2005; Štverák et al. 2008; Pierrard et al. 2016).
Recent studies on global Kappa models have emphasized two distinct alternatives, namely, assuming the temperature (fitting) parameter, either dependent or independent of the κindex (Lazar et al. 2015a, 2016a). Modeling with a κindependent temperature is convenient in computations, and it is therefore very much invoked in theoretical predictions on dispersion and stability properties of Kappa distributed plasmas; we refer to the reviews by Hellberg et al. (2005) and Pierrard & Lazar (2010). The effects of a varying κ (including the Maxwellian limit κ → ∞) are also easier to parametrize if temperature is taken as constant. However, for studies intended to outline the effects of suprathermal populations, Lazar et al. (2015a, 2016a) have shown that only the approach assuming a κdependent temperature may provide a fairly rigorous comparison of the global Kappa describing the observed distribution, with the cooler Maxwellian κ → ∞ limit reproducing the core of the distribution. Moreover, only in this case are the instabilities driven by the kinetic anisotropies of plasma particles systematically stimulated by the suprathermals (i.e., lowering the values of κ), confirming the expectation that these populations should contribute with an additional free energy (Lazar et al. 2015a, 2016b). From a different perspective, including, for instance, the origin of Kappa distributions, there are processes responsible for the formation of suprathermal tails, for example, particle energization by the wave turbulence, which may imply a κindependent evolution of the velocity distributions (Yoon 2014; Lazar et al. 2016a). In this case, the highenergy tails are counterbalanced by a lowvelocity enhancement, so that the temperature remains indeed independent of parameter κ.
A global Kappa distribution implies a reduced number of parameters and is, therefore, easy to manipulate in observational and theoretical analyses. However, a global Kappa does not always provide a good fit to the observed distributions, as we see here below. For example, the analysis in Sect. 2 (also Figs. 1 and 2), suggests that a global Kappa is just a simplified (or idealized) approach that artificially constrains the core and suprathermal populations to be described by the same parameters, for example, density, temperature, and anisotropy (not fully justified yet, although some studies in this direction have been done by Leubner 2004; Lazar et al. 2015a, 2016a). It may be, therefore, that fitting measurements of the observed distributions with global Kappa models can not provide relevant evidence to support any of the two alternatives to temperature, either dependent or independent of the powerindex κ.
Fig. 1 Nonstreaming electron velocity distribution F [cm^{6} s^{3}] measured by Ulysses in the slow wind in 2002 on DOY 15 03:33:42: a) parallel and perpendicular cuts; b) Maxwellian fit to the core and Kappa fit to the halo suprathermal tails; c) combined (resultant) MaxwellianKappa model fitting the entire distribution; d) global Kappa model fitting the entire distribution. 

Open with DEXTER 
Fig. 2 Nonstreaming electron velocity distribution F [cm^{6} s^{3}] measured by Ulysses in the slow wind in 2002 on DOY 19 08:19:49: a) parallel and perpendicular cuts; b) Maxwellian fit to the core and Kappa fit to the halo suprathermal tails; c) combined (resultant) MaxwellianKappa model fitting the entire distribution; d) global Kappa model fitting the entire distribution. 

Open with DEXTER 
Fig. 3 Asymmetric electron distribution F [cm^{6} s^{3}] measured by Ulysses in 2002 DOY 18 18:50:36: a) Maxwellian fit to the core, Kappa fit the halo suprathermal tails, and a drifting Kappa fit to the streaming (strahl) component; b) combined model fitting the entire distribution. 

Open with DEXTER 
More complex models that combine distinct components are, in general, more accurate in reproducing the distributions measured in the solar wind, and may, therefore, unveil details of the observed distributions. For instance, the nonstreaming distributions detected in the slow winds and at sufficiently large distances from the Sun (Maksimovic et al. 2005) are expected to be better described by a dual model comprising a standard Maxwellian distribution function to reproduce the core at low energies, and a Kappa powerlaw to fit the suprathermal tails (halo component) of the distribution (see Figs. 1 and 2). In the fast winds, the observed distributions become asymmetric (see Fig. 3), and to describe them we need an additional component (called strahl) streaming along the magnetic field direction. Such complex models have already been invoked in observational analyses in order to quantify the main physical properties of these components, for example, densities, relative drifts, temperatures, and even temperature anisotropies (Maksimovic et al. 2005; NievesChinchilla & Viñas 2008; Štverák et al. 2008; Viñas et al. 2010; Pierrard et al. 2016). The halo component can be called suprathermal as long it is well reproduced by Kappa powerlaws with sufficiently low values of κ, that is, κ< 10. Otherwise, for higher values, such as κ> 10, the halo component appears to be thermalized and well described by a Maxwellian, but it remains a distinct component, being, in general, less dense but hotter than the core. In this case, the overall distribution can be described by a twoMaxwellian model, as a limiting κ → ∞ case of a dual MaxwellianKappa representation. Although the existence of a halo component with κ> 10 is not so evident from observations, a twoMaxwellian model has already been invoked in theoretical predictions on the dispersion and stability of space plasmas (Gary et al. 1975, 1994) as well as in observational analyses (Feldman et al. 1975; Pilipp et al. 1987) and in comparisons with global Kappa (Maksimovic et al. 1997). Note that this Maxwellian limit (κ → ∞) of a Kappa describing only the halo component is completely distinct from the Maxwellian fit to the core, and a correlation between these two is not justified in this case (as for a global Kappa modeling, where the Maxwellian limit may need to reproduce the core of the distribution to compare with the global Kappa and emphasize the effects of the suprathermals).
When the core and suprathermal components are anisotropic, for example, due to temperature anisotropies or relative drifts, a dual MaxwellianKappa modeling becomes complicated to deal with, even in numerical computations. Recently, important progress has been made for overcoming this challenge in studies of the stability and relaxation of a dual corehalo plasma with intrinsic temperature anisotropies (Lazar et al. 2014, 2015b; Shaaban et al. 2016). Both interpretations were adopted for the anisotropic halo modeled by a biKappa with temperatures being either dependent on or independent of the powerindex κ. The results from these studies were again favorable to a Kappa model with a κdependent temperature, since only this approach provided a systematic stimulation of the instabilities in the presence of suprathermals.
The present paper aims to provide additional insight into these properties of Kappa distributed plasmas as they are revealed by observational data. As discussed above, there are two fitting models widely invoked in theory and observations, namely, the global Kappa and the dual MaxwellianKappa models. In Sect. 2 we describe them mathematically and qualitatively analyze their accuracy in reproducing the solar wind observational data. Particularly relevant for our analysis is the Kappa model that better reproduces the observed distributions. This model is then considered for fitting to more than 100 000 velocity distributions measured in the ecliptic within a wide interval of heliocentric distance between 0.3 and 4.0 AU. Determined as fitting parameters, the temperature and the powerindex κ are subjected to a statistical analysis that enables us to extract the relationship between them. The κdependency indicated by the observations for the temperature of Kappa populations is used in Sect. 3 for a kinetic analysis of the electromagnetic instabilities driven by the temperature anisotropy. These instabilities can explain the enhanced fluctuations observed in the solar wind; we refer to the review by Zimbardo et al. (2010) and references therein. The results of our study are summarized and discussed in Sect. 4.
2. Refined models from the observations
2.1. Global Kappa versus dual MaxwellianKappa
Fluxes of plasma particles measured insitu in the solar wind are transformed into the frame of bulk flow, such that, in the absence of any differential streaming, the velocity distributions are almost symmetric in the Sunward and antiSunward directions, without indications of any additional strahl (or streaming) component (see examples in Figs. 1 and 2). One intriguing feature that makes them impossible to describe with standard Maxwellian models is the enhanced highenergy tails, known as suprathermal tails. Instead, the Kappa (Lorentzian) powerlaws are used to reproduce the suprathermal tails of the observed distributions, and this can be done in two distinct ways. Thus, a Kappa powerlaw (subscript κ) is used either as a global model (subscript GK) to fit the entire distribution (1)or to describe only the suprathermal halo (subscript h) component (2)In the second case, the core (subscript c) fitting is a Maxwellian (subscript M) (3)and the resulting model to describe the entire distribution f is a dual MaxwellianKappa (subscript MK) (4)Since kinetic modeling deals with homogeneous plasma systems in general, in the above n_{c} = N_{c}/N and n_{h} = N_{h}/N are the relative number densities of the core and halo components, respectively (where N is the total number of plasma particles of a given species). All the distribution functions are normalized to unity. In magnetized plasmas from space, the velocity distributions are in general gyrotropic (without a major anisotropy in the plane transverse to the magnetic field) allowing us to describe them in polar coordinates (v_{x},v_{y},v_{z}) = (v_{⊥}cosφ,v_{⊥}sinφ,v_{∥}) with ∥ and ⊥ denoting directions relative to the stationary magnetic field. Gyrotropic distribution functions reduce the analysis to only two variables in velocity space and enable us to quantify the principal anisotropies reported by the observations, for example, the biaxis temperature anisotropies T_{⊥} ≠ T_{∥} combined, or not, with the fieldaligned streams (also known as strahls).
Relevant for the nonstreaming distributions observed in the solar wind is the wellknown biKappa distribution function (Summers & Thorne 1991; Lazar et al. 2015a) (5)where v_{∥} and v_{⊥} denote particle velocity parallel and perpendicular with respect to a largescale magnetic field, m is particle mass, k_{B} is the Boltzmann constant, Γ is the Gamma function, the powerindex κ ∈ (3/2,∞ ], and T_{∥ , ⊥} and θ_{∥ , ⊥} are the corresponding temperatures and thermal velocities, which are related by (6)(7)This biKappa model is invoked in both fitting techniques, either as a global model to describe the entire distribution (Vasyliunas 1968; Christon et al. 1989; Collier et al. 1996; Maksimovic et al. 1997), or to fit partially only the suprathermal tails of the distribution (Maksimovic et al. 2005; NievesChinchilla & Viñas 2008; Štverák et al. 2008; Pierrard et al. 2016).
The lowenergy core is well reproduced by a biMaxwellian (e.g., in Pilipp et al. 1987) (8)where w_{∥ , ⊥} are thermal velocities defined by the temperatures as moments of second order (9)(10)The socalled global Kappa approach from Eq. (1) implies a reduced number of parameters and is, therefore, preferred in computations, while a more complex, dual MaxwellianKappa as defined in Eq. (4) is expected to better reproduce the observations, although direct comparisons with global Kappa fittings are not reported yet (at least to our knowledge). In Figs. 1 and 2 we compare these two fitting models as obtained for the electron velocity distributions measured by Ulysses in January 2002 (dots), two events on DOY 15 03:33:42 and DOY 19 08:19:49, respectively. These distributions are chosen to be almost isotropic, that is, with insignificant differences between parallel (red dots) and perpendicular (black dots) cuts (see top panels a), such that fittings to the perpendicular cut are representative for the entire distribution. The velocity distributions in Figs. 1 and 2 are constructed based on the Ulysses/SWOOPS electron pitchangle distributions published by the ESA Web service^{1}. These data are in the solar wind frame and have been corrected for the spacecraft potential (Scime et al. 1994).
Fig. 4 Scatter plots (blue dots) of the halo temperature T [K] vs. κindex, measured in the ecliptic (0.3 <R< 4 AU): total temperature (top), and the components, parallel (middle) and perpendicular temperature (bottom). The temperatures predicted by the Eqs. (15)–(17) as functions of κ parameter are overplotted with dashed lines. 

Open with DEXTER 
The global Kappa fit is displayed in the bottom panels (solid blue lines), and the dual MaxwellianKappa fit is explicitly shown in the middle panels: in the second panel from the top we have the core (solid red lines) and halo (blue dashed lines) fits and these two are then combined as a dual MaxwellianKappa fit in the third panel (solid green lines). Now, a direct comparison of these two fitting models becomes possible simply by analyzing the last two bottom panels in Figs. 1 and 2. It is obvious that in both events the observed distribution is better reproduced by a dual MaxwellianKappa approach than a global Kappa. We note that in Fig. 2 a global Kappa provides a better fit to the observed distribution than in Fig. 1, but in both these two cases the global Kappa fits are visibly less accurate than those obtained with a dual MaxwellianKappa approach. In this case the effects of suprathermals may simply be outlined by a direct comparison between the dual model f_{MK} and the Maxwellian fit f_{M} to the core, and there is no need to consider the Maxwellian limit κ → ∞ of the Kappa component.
The fact that a dual model is more accurate may simply suggest the existence of two distinct components, namely, a lowenergy core population and a suprathermal Kappadistributed halo. There is probably no other possibility to make a distinction between these components, which may have different origins in the solar corona (Pierrard et al. 1999) or in the solar wind (Maksimovic et al. 2005). On the other hand, these two components are expected to be strongly correlated, especially by the collective effects of plasma particles which should tend to diminish differences between them. However, this is not always confirmed or evident from observations, which show, for instance, an electron halo enhancing with the radial distance from the Sun at the expense of the more anisotropic strahl, which is diminished in almost the same measure (Maksimovic et al. 2005). In this case an explanation may be offered by the instabilities and amplified fluctuations selfgenerated by the strahl, which can isotropize the strahl, scattering the particles and thus feeding the more isotropic suprathermal halo. A dual corehalo model is relevant not only in the slow winds when the flowing (bulk) speed is sufficiently low (i.e., V_{SW}< 360 km s^{1}), and the relative density of the strahl is negligibly small n_{s}<n_{h}<n_{c} (Lazar et al. 2015b), but according to Maksimovic et al. (2005) it may also describe any slow or fast wind conditions at sufficiently large heliocentric distances, for example, beyond 1 AU where the strahl component is markedly reduced. A dual MaxwellianKappa remains a major component in the fast winds, when the velocity distribution exhibits an additional streaming component (called strahl). In Fig. 3 we display an example of velocity distribution of electrons measured during a fast solar wind by Ulysses in January 2002. In this case, the best fitting model contains three components: a Maxwellian core, a Kappa halo, and a driftingKappa strahl (or beam of particles).
2.2. The Kappa (halo) temperature from observations
If a κdependency exists for the temperature of Kappa population, then this must be revealed by the observational data. In order to be relevant, the Kappa model should accurately reproduce the measured distribution, and here above we have demonstrated that this condition is satisfied by a dual MaxwellianKappa model, with Kappa distribution function fitting only the suprathermal (halo) tails.
Štverák et al. (2008) have used the same dual MaxwellianKappa model from Eq. (4) to determine the principal moments of the electron distributions measured in the solar wind. The dual model was fitted to more than 100 000 velocity distributions measured by three different spacecraft missions (Helios 1, Cluster II, and Ulysses) in the ecliptic at different heliocentric distances in the interval 0.3–3.95 AUs. Details about the electron analyzers and the methods of correction and reconstruction of the 3D velocity distributions are given in Štverák et al. (2008) and some references therein. Here we make use of the same dataset, and investigate the temperature parameters determined for the Kappa (halo) populations. Figure 4 displays scatter plots (with blue dots) of the temperatures determined for the electron halo (subscript h) versus the corresponding values of the powerindex κ: total temperature (T_{h}) in the top panel, parallel (T_{h, ∥}) and perpendicular (T_{h, ⊥}) temperatures in the middle and bottom panels, respectively. One dimensional fits provide temperature components T_{h, ∥} and T_{h, ⊥}, enabling us to calculate T_{h} = T_{h, ∥}/ 3 + 2T_{h, ⊥}/ 3.
For each panel in Fig. 4 the spread of data points is very similar and we can analyze them generically. The data points concentrate along a curve which clearly shows that temperature is not constant but varies as a function of κ. Thus, the temperature of Kappa population decreases with increasing powerindex κ, with the allure of an asymptotic decrease towards the Maxwellian limit κ → ∞. In this limit of a very large κ → ∞, the biKappa reproducing the halo components reduces to a biMaxwellian (similar to that reproducing the core component). However, Lazar et al. (2016a) have shown that any Kappa distribution function admits two distinct Maxwellian limits:

(a)
A cooler Maxwellian limit corresponding to the approachwith a κdependent temperature (11)with thermal velocities θ_{∥ , ⊥} determined in this case by lower temperatures

(b)
The second limit corresponds to the approach with a constant temperature, not dependent on κ: (14)where and .
The observations indicate a κdependent temperature, when the components T_{∥ , ⊥} can be defined in terms of their Maxwellian limits as follows (15)(16)and a similar relation is obtained for the total temperature (17)Using these relations to fit to the observations (see the dashed curves in Fig. 4), we find for the Maxwellian limits K. We should point out that the accumulation of data aligns with these dashed curves very well, which confirms the validity of κdependency found in Eqs. (15)–(17). However, the powerindex κ describing the electron halo population in the solar wind does not take very large values, being limited to κ< 10. Therefore, for a rigorous interpretation of the solar wind electrons, a twoMaxwellian model does not seem to be justified, but can be invoked as a limiting case of a dual MaxwellianKappa in order to depict and investigate the effects of the suprathermal populations.
Fig. 5 Scatter plots (blue dots) of the halo (total) temperature T [K] vs. κindex, measured in different intervals of heliocentric distances R (selected from the data set in Fig. 4): a) 0.3 <R< 0.4 AU; b) <0.6 AU; c) ~1 AU (0.9 AU AU); d) R> 1.1 AU. Fitting curves are explained in the text. 

Open with DEXTER 
Both parameters T and κ analyzed here characterize the halo (Kappa) electron component of the solar wind plasma, and should therefore undergo a variation with the radial (heliocentric) distance from the Sun, being naturally conditioned by the solar wind evolution during its expansion in space. Thus, recent studies of the same set of data (Pierrard et al. 2016) have shown significant differences between the core and halo temperatures, and implicitly between their evolutions with heliocentric distance; the core temperature decreases (with a tendency of stabilizing after 2 AU) while the halo temperature increases with radial distance from the Sun (with an apparent saturation at approx. 3 AU). Normally, the temperature should indeed cool down with the distance from the Sun, as the core temperature is confirming. The increase of the halo temperature is somehow unusual, although it is in perfect agreement with the evolution of the other parameters reported by different observations (Maksimovic et al. 2005; Pierrard et al. 2016), namely, the increase of halo population (i.e., the increase of its relative number density n_{h}/n_{c}) at the same time with a decrease of κ with radial distance. The halo population seems to be enhanced at the expense of the strahl population, which is eventually scattered and isotropized by the selfgenerated instabilities (collisions are inefficient at large distances in the solar wind). At this point it therefore becomes interesting to check whether or not such processes may have an influence on the κdependency shown by the halo temperature in the observations.
Figure 5 is intended to unveil the existence of these influences if they exist, by presenting in detail a radial profile of the halo temperature as a function of the powerindex κ. We made a selection of the data from different heliocentric distances, and plot in four panels with the distance increasing from top panel (a) to bottom panel (d) as follows: (a) 0.3 <R< 0.4 AU; (b) <0.6 AU; (c) ~1 AU ( AU); (d) R> 1.1 AU. At small radial distances, 0.3 <R< 0.4 AU, the data accumulates and aligns very well to the curve defined by Eq. (17) for a κvariation of the halo temperature. With increasing distance, the data points spread above and below this fitting curve, most probably as an effect of some mechanisms at work in the solar wind expansion, such as cooling, magnetic focusing, and particle scattering by the fluctuations, and so on. However, up to 1 AU, the data aligns well, on average, on the κdependency in Eq. (17) with the same asymptotic (Maxwellian) limit K. We cannot claim the same thing for the data collected beyond 1 AU, which, on average, seems to align to the same κdependency law as the one given in Eq. (17), but with a lower Maxwellian (asymptotic) limit K (see the solid (red) line in panel (d)). Moreover, the spread of data above the dashed line suggests the existence of other hotter populations that may fit to another similar law but with a higher Maxwellian (asymptotic) limit K (see the dotted (blue) line in panel (d)). In order to explain these distinct populations indicated by different fitting curves (and leading to a spread of data along the average fitting) we would plausibly need to invoke the different origins of these populations. For instance, scattering of the beaming strahl by the selfgenerated instabilities may provide hotter Kappa populations (Maksimovic et al. 2005), as also discussed here above. We can conclude stating that the observations in the solar wind show clear evidences that support a κdependency of the temperature of Kappadistributed populations, as the one explicitly shown above in Eqs. (15)–(17).
3. Kappa electrons: destabilizing effects
In this section we extend our comparative study to the effects of suprathermal (halo) electrons on the temperature anisotropy instabilities as they are predicted by the same dual MaxwellianKappa model when the temperatures of Kappa populations are dependent or independent of κ. Also considered is a comparison with the Maxwellian limit κ → ∞ that enables us to emphasize the effects of suprathermal electrons which are expected to stimulate the instability.
According to the reports on the temperature anisotropy in the solar wind (Štverák et al. 2008; Pierrard et al. 2016) both the core and halo populations may exhibit anisotropies, and the halo is, in general, more anisotropic being hotter and less dense than the core. Therefore, in order to depict the effects of the suprathermal populations, here we minimize the influence of the Maxwellian core assuming isotropic with T_{c, ∥} = T_{c, ⊥} = T_{c} (i.e., A_{c} = T_{c, ⊥}/T_{c, ∥} = 1). In the solar wind, the electron halo may exhibit both deviations from isotropy, namely, an excess of temperature in parallel direction A_{e,h}< 1 that can drive the socalled electron firehose (EFH) instability, or an excess of temperature in perpendicular direction A_{e,h}> 1 that can be at the origin of the whistler instability (WI), also known as the electromagnetic electroncyclotron (EMEC) instability. For the sake of simplicity, here we investigate only the effects of suprathermal electrons on the parallel electromagnetic modes (assuming a fieldaligned wavevector k = k_{∥}). We have analyzed the unstable solutions for a significant number of cases, keeping constant the core parameters for average values measured in the solar wind, for example, A_{e,c} = 1, β_{e,c} = 1, the relative halocore density η = 0.05, and varying the halo parameters, A_{e,h}, β_{e,h}, and κ.
3.1. Firehose instability
We first discuss the electron firehose instability (EFHI), which is driven by a temperature anisotropy A_{e,h} = T_{h, ⊥}/T_{h, ∥}< 1 of the electron halo component. The contribution of protons is minimized by considering them isotropic (both the proton core and halo components), such that the dispersion relation for the EFHI can take the following form
(18)where , , κ = κ_{e,p}, α = (1−1.5 /κ)^{0.5}, μ = m_{p}/m_{e} is the protonelectron mass ratio, η is the halocore relative density, A_{j,h} and A_{j,c} are the halo and core temperature anisotropies for protons (subscript j = p) or electrons (subscript j = e), respectively, and are the plasma beta parameters for different plasma components, Z_{j,M} is the standard dispersion function for (bi)Maxwellian distributed plasmas, and Z_{j,κ} is the modified dispersion function for (bi)Kappa distributed plasmas (see Appendix A for the definitions of these functions).
Fig. 6 EFHI: growth rates (top) and real frequencies (bottom) derived for a dual model comprising a Maxwellian core and a biKappa halo (κ = 2) with κdependent temperatures (solid lines), κindependent temperatures (dashed lines), or their biMaxwellian limit (κ → ∞, dottedlines). The plasma parameters for electrons are given explicitly in each panel and for protons are A_{p,c} = A_{p,h} = 1, β_{p,c} = 1, β_{p,h} = 4. 

Open with DEXTER 
Fig. 7 The same as in Fig. 6, but for a lower plasma beta β_{e,h} = 1. The plasma parameters for electrons are given explicitly in each panel and for protons are A_{p,c} = A_{p,h} = 1, β_{p,c} = β_{p,h} = 1. 

Open with DEXTER 
Fig. 8 WI: growth rates (top) and real frequencies (bottom) derived for a dual model comprising a Maxwellian core and a biKappa halo (κ = 2) with κdependent temperatures (solid lines), κindependent temperatures (dashed lines), or their biMaxwellian limit (κ → ∞, dottedlines). The electron plasma parameters are given explicitly in each panel. 

Open with DEXTER 
Fig. 9 The same as in Fig. 8, but for a lower halo anisotropy A_{e,h} = 1.1, and a higher plasma beta β_{e,h} = 1. 

Open with DEXTER 
Since the plasma beta parameter is an important factor in triggering the EFH instability by the anisotropic electrons with A_{e}< 1, the unstable solutions displayed in Figs. 6 and 7 are derived for two distinct conditions, respectively; for a high halo plasma beta β_{e,h} = 4, and for a lower β_{e,h} = 1. For each of these two cases, the comparison invokes both Kappa approaches (κ = 2) with temperatures dependent or independent of κ, as well as their Maxwellian limit (κ → ∞). The growth rates (top panels) predicted by a Kappa model with κdependent temperatures are markedly enhanced in the presence of suprathermals. For instance, in Fig. 7 only this model predicts an instability, while the system becomes stable within the other two approaches. A systematic stimulation of the EFH instability in the presence of suprathermal electrons is confirmed only by a Kappa model with κdependent temperatures. The corresponding wavefrequencies (bottom panels) also show distinct and important variations with the Kappa models and the powerindex κ.
3.2. Electron cyclotron (whistler) instability
In the opposite case, an excess of perpendicular temperature A_{e,h} = T_{h, ⊥}/T_{h, ∥}> 1 may drive the whistler instability. Since the wave frequency of the whistler modes is high (ω ≫ Ω_{p}), the protons do not react and the dispersion relation that describes the whistler instability (WI) in a dual MaxwellianKappa plasma reads (Lazar et al. 2015b) (19)where K = kc/ω_{h,e}, W = ω/  Ω_{e} , and the other quantities are the same as explained above for the EFHI.
Figures 8 and 9 present two representative conditions for the WI, respectively, when the halo plasma beta is low β_{e,h} = 0.05 but the halo anisotropy is high A_{e,h} = 3, and for a higher β_{e,h} = 1 and a lower anisotropy A_{e,h} = 1.1. For even higher betas the instability can be triggered by very low anisotropies close to the isotropy condition A_{e,h} → 1. Top panels display the growthrates, and the corresponding wavefrequencies are plotted in the bottom panels. The unstable solutions are derived for both Kappa approaches (κ = 2), with temperatures dependent or independent of κ, and for their Maxwellian limit (κ → ∞) enabling us to emphasize the effects of suprathermal populations. Comparing the growth rates from different models, a systematic stimulation of the instability in the presence of suprathermals can be inferred only for the approach with κdependent temperatures. The other Kappa approach assuming κindependent temperature leads to an irregular variation of the growth rates with κ. Unlike the growthrates, the wave frequencies are not markedly influenced by the Kappa models and the variation of κ.
4. Conclusions
The temperature and powerindex κ are the main fitting parameters in a Kappa distribution model widely invoked to describe the velocity distributions of plasma particles and their dynamics. Recent studies have raised an important question on two alternative approaches, considering the temperature of Kappa populations either dependent or independent of κ. Various scenarios can provide justification for each of these two approaches, but observational evidences supporting any of them are not reported yet. Here we have shown for the first time that Kappa populations observed in the solar wind exhibit a direct interdependence between T and κ. Thus, the temperature decreases with increasing κ(20)converging asymptotically to a limit that can be considered the cooler Maxwellian limit , defined (by the components) in Eqs. (12), (13). This is exactly the variation predicted by the previous studies (Lazar et al. 2015a, 2016a; Pierrard et al. 2016). It is important to emphasize that observational data invoked in the present paper are obtained with a refined dual MaxwelliaKappa model, which is more accurate than a global Kappa in reproducing the observed distributions and implicitly the Kappa populations. We have also used these two approaches with temperatures dependent or independent of κ to describe and compare the effects of suprathermal electrons on two electromagnetic instabilities driven by the temperature anisotropy, namely, the firehose and whistler instabilities. Only the approach with κdependent temperatures may confirm the expectation predicting a systematic stimulation of the instabilities by increasing the presence of suprathermals, while a Kappa approach with constant temperatures leads to a questionable variation of the growth rates with κ.
http://ufa.esac.esa.int/ufa/, author Ruth Skoug, Los Alamos National Laboratory.
Acknowledgments
The authors acknowledge use of the Ulysses/SWOOPS electron pitchangle distributions published by the ESA Web service (http://ufa.esac.esa.int/ufa/, author Ruth Skoug, Los Alamos National Laboratory), and support from the Katholieke Universiteit Leuven, the RuhrUniversität Bochum, and the Deutsche Forschungsgemeinschaft (DFG). These results were obtained in the framework of the projects GOA/2015014 (KU Leuven), G.0A23.16N (FWOVlaanderen) and C 90347 (ESA Prodex). This research has been funded by the Interuniversity Attraction Poles Programme initiated by the Belgian Science Policy Office (IAP P7/08 CHARM). S.M. Shaaban would like to thank the Egyptian Ministry of Higher Education for supporting his research activities. Thanks are due to Š. Štverák for providing the observational data.
References
 Christon, S. P., Williams, D. J., Mitchell, D. G., Frank, L. A., & Huang, C. Y. 1989, J. Geophys. Res., 94, 13409 [NASA ADS] [CrossRef] [Google Scholar]
 Collier, M. R., Hamilton, D. C., Gloeckler, G., Bochsler, P., & Sheldon, R. B. 1996, Geophys. Res. Lett., 23, 1191 [NASA ADS] [CrossRef] [Google Scholar]
 Feldman, W. C., Asbridge, J. R., Bame, S. J., Montgomery, M. D., & Gary, S. P. 1975, J. Geophys. Res., 80, 4181 [NASA ADS] [CrossRef] [Google Scholar]
 Fried, B. D., & Conte, S. D. 1961, The plasma dispersion function: the Hilbert transform of the Gaussian (New York: Academic Press) [Google Scholar]
 Gary, S. P., Feldman, W. C., Forslund, D. W., & Montgomery, M. D. 1975, J. Geophys. Res., 80, 4197 [NASA ADS] [CrossRef] [Google Scholar]
 Gary, S. P., Moldwin, M. B., Thomsen, M. F., Winske, D., & McComas, D. J. 1994, J. Geophys. Res., 99, 23 [Google Scholar]
 Hellberg, M., Mace, R., & Cattaert, T. 2005, Space Sci. Rev., 121, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Lazar, M., Schlickeiser, R., & Shukla, P. K. 2008, Phys. Plasmas, 15, 042103 [NASA ADS] [CrossRef] [Google Scholar]
 Lazar, M., Poedts, S., & Schlickeiser, R. 2014, J. Geophys. Res. (Space Physics), 119, 9395 [NASA ADS] [CrossRef] [Google Scholar]
 Lazar, M., Poedts, S., & Fichtner, H. 2015a, A&A, 582, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lazar, M., Poedts, S., Schlickeiser, R., & Dumitrache, C. 2015b, MNRAS, 446, 3022 [NASA ADS] [CrossRef] [Google Scholar]
 Lazar, M., Fichtner, H., & Yoon, P. H. 2016a, A&A, 589, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lazar, M., Shaaban, S. M., Poedts, S., & Štverák, Š. 2016b, MNRAS [Google Scholar]
 Leubner, M. P. 2004, ApJ, 604, 469 [NASA ADS] [CrossRef] [Google Scholar]
 Livadiotis, G., & McComas, D. J. 2013, Space Sci. Rev., 175, 183 [NASA ADS] [CrossRef] [Google Scholar]
 Maksimovic, M., Pierrard, V., & Riley, P. 1997, Geophys. Res. Lett., 24, 1151 [NASA ADS] [CrossRef] [Google Scholar]
 Maksimovic, M., Zouganelis, I., Chaufray, J.Y., et al. 2005, J. Geophys. Res. (Space Physics), 110, A09104 [NASA ADS] [CrossRef] [Google Scholar]
 NievesChinchilla, T., & Viñas, A. F. 2008, J. Geophys. Res. (Space Physics), 113, A02105 [NASA ADS] [CrossRef] [Google Scholar]
 Pierrard, V., & Lazar, M. 2010, Sol. Phys., 267, 153 [NASA ADS] [CrossRef] [Google Scholar]
 Pierrard, V., Maksimovic, M., & Lemaire, J. 1999, J. Geophys. Res., 104, 17021 [NASA ADS] [CrossRef] [Google Scholar]
 Pierrard, V., Lazar, M., Poedts, S., et al. 2016, Sol. Phys., 291, 2165 [NASA ADS] [CrossRef] [Google Scholar]
 Pilipp, W. G., Muehlhaeuser, K.H., Miggenrieder, H., Montgomery, M. D., & Rosenbauer, H. 1987, J. Geophys. Res., 92, 1075 [NASA ADS] [CrossRef] [Google Scholar]
 Scime, E. E., Phillips, J. L., & Bame, S. J. 1994, J. Geophys. Res., 99, 14 [Google Scholar]
 Shaaban, S. M., Lazar, M., Poedts, S., & Elhanbaly, A. 2016, J. Geophys. Res.: Space Physics, 121, 6031 [NASA ADS] [CrossRef] [Google Scholar]
 Summers, D., & Thorne, R. M. 1991, Phys. Fluids B, 3, 1835 [NASA ADS] [CrossRef] [Google Scholar]
 Štverák, Š., Trávníček, P., Maksimovic, M., et al. 2008, J. Geophys. Res. (Space Physics), 113, A03103 [NASA ADS] [CrossRef] [Google Scholar]
 Vasyliunas, V. M. 1968, J. Geophys. Res., 73, 2839 [NASA ADS] [CrossRef] [Google Scholar]
 Viñas, A., Gurgiolo, C., NievesChinchilla, T., Gary, S. P., & Goldstein, M. L. 2010, Twelfth International Solar Wind Conference, 1216, 265 [NASA ADS] [Google Scholar]
 Yoon, P. H. 2014, J. Geophys. Res. (Space Physics), 119, 7074 [NASA ADS] [CrossRef] [Google Scholar]
 Zimbardo, G., Greco, A., SorrisoValvo, L., et al. 2010, Space Sci. Rev., 156, 89 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Plasma dispersion functions
For the biMaxwellian distributed plasma components (e.g., the core of different species, i.e., electrons with subscript j = e and protons with subscript j = p), the plasma dispersion function in Eqs. (18) and (19) takes the standard form (Fried & Conte 1961) (A.1)
of argument . For the biKappa distributed components, in the dispersion relations (18) and (19) we use the modified Kappa dispersion function (Lazar et al. 2008) (A.2)of argument In these equations ± denotes the circular polarizations, righthanded (RH) and lefthanded (LH), respectively, and Ω_{j} is the (nonrelativistic) gyrofrequency.
All Figures
Fig. 1 Nonstreaming electron velocity distribution F [cm^{6} s^{3}] measured by Ulysses in the slow wind in 2002 on DOY 15 03:33:42: a) parallel and perpendicular cuts; b) Maxwellian fit to the core and Kappa fit to the halo suprathermal tails; c) combined (resultant) MaxwellianKappa model fitting the entire distribution; d) global Kappa model fitting the entire distribution. 

Open with DEXTER  
In the text 
Fig. 2 Nonstreaming electron velocity distribution F [cm^{6} s^{3}] measured by Ulysses in the slow wind in 2002 on DOY 19 08:19:49: a) parallel and perpendicular cuts; b) Maxwellian fit to the core and Kappa fit to the halo suprathermal tails; c) combined (resultant) MaxwellianKappa model fitting the entire distribution; d) global Kappa model fitting the entire distribution. 

Open with DEXTER  
In the text 
Fig. 3 Asymmetric electron distribution F [cm^{6} s^{3}] measured by Ulysses in 2002 DOY 18 18:50:36: a) Maxwellian fit to the core, Kappa fit the halo suprathermal tails, and a drifting Kappa fit to the streaming (strahl) component; b) combined model fitting the entire distribution. 

Open with DEXTER  
In the text 
Fig. 4 Scatter plots (blue dots) of the halo temperature T [K] vs. κindex, measured in the ecliptic (0.3 <R< 4 AU): total temperature (top), and the components, parallel (middle) and perpendicular temperature (bottom). The temperatures predicted by the Eqs. (15)–(17) as functions of κ parameter are overplotted with dashed lines. 

Open with DEXTER  
In the text 
Fig. 5 Scatter plots (blue dots) of the halo (total) temperature T [K] vs. κindex, measured in different intervals of heliocentric distances R (selected from the data set in Fig. 4): a) 0.3 <R< 0.4 AU; b) <0.6 AU; c) ~1 AU (0.9 AU AU); d) R> 1.1 AU. Fitting curves are explained in the text. 

Open with DEXTER  
In the text 
Fig. 6 EFHI: growth rates (top) and real frequencies (bottom) derived for a dual model comprising a Maxwellian core and a biKappa halo (κ = 2) with κdependent temperatures (solid lines), κindependent temperatures (dashed lines), or their biMaxwellian limit (κ → ∞, dottedlines). The plasma parameters for electrons are given explicitly in each panel and for protons are A_{p,c} = A_{p,h} = 1, β_{p,c} = 1, β_{p,h} = 4. 

Open with DEXTER  
In the text 
Fig. 7 The same as in Fig. 6, but for a lower plasma beta β_{e,h} = 1. The plasma parameters for electrons are given explicitly in each panel and for protons are A_{p,c} = A_{p,h} = 1, β_{p,c} = β_{p,h} = 1. 

Open with DEXTER  
In the text 
Fig. 8 WI: growth rates (top) and real frequencies (bottom) derived for a dual model comprising a Maxwellian core and a biKappa halo (κ = 2) with κdependent temperatures (solid lines), κindependent temperatures (dashed lines), or their biMaxwellian limit (κ → ∞, dottedlines). The electron plasma parameters are given explicitly in each panel. 

Open with DEXTER  
In the text 
Fig. 9 The same as in Fig. 8, but for a lower halo anisotropy A_{e,h} = 1.1, and a higher plasma beta β_{e,h} = 1. 

Open with DEXTER  
In the text 