Kinetic corrections from analytic nonMaxwellian distribution functions in magnetized plasmas.
Abstract
In magnetized plasma physics, almost all developed analytic theories assume a Maxwellian distribution function (MDF) and in some cases small deviations are described using the perturbation theory. The deviations with respect to the Maxwellian equilibrium, called kinetic effects, are required to be taken into account specially for fusion reactor plasmas. Generally, because the perturbation theory is not consistent with observed steadystate nonMaxwellians, these kinetic effects are numerically evaluated by very CPUexpensive codes, avoiding the analytic complexity of velocity phase space integrals. We develop here a new method based on analytic nonMaxwellian distribution functions constructed from nonorthogonal basis sets in order to (i) use as few parameters as possible, (ii) increase the efficiency to model numerical and experimental nonMaxwellians, (iii) help to understand unsolved problems such as diagnostics discrepancies from the physical interpretation of the parameters, and (iv) obtain analytic corrections due to kinetic effects given by a small number of terms and removing the numerical error of the evaluation of velocity phase space integrals. This work does not attempt to derive new physical effects even if it could be possible to discover one from the better understandings of some unsolved problems, but here we focus on the analytic prediction of kinetic corrections from analytic nonMaxwellians. As applications, examples of analytic kinetic corrections are shown for the secondary electron emission, the Langmuir probe characteristic curve, and the entropy. This is done by using three analytic representations of the distribution function: the Kappa (KDF), the bimodal or a new interpreted nonMaxwellian (INMDF) distribution function. The existence of INMDFs is proved by new understandings of the experimental discrepancy of the measured electron temperature between two diagnostics in JET. As main results, it is shown that (i) the empirical formula for the secondary electron emission is not consistent with a MDF due to the presence of superthermal particles, (ii) the superthermal particles can replace a diffusion parameter in the Langmuir probe current formula and (iii) the entropy can explicitly decrease in presence of sources only for the introduced INMDF without violating the second law of thermodynamics. Moreover, the first order entropy of an infinite number of superthermal tails stays the same as the entropy of a MDF. The latter demystifies the Maxwell’s demon by statistically describing nonisolated systems.
The observation of nonequilibrium distribution functions is possible in a large number of areas other than laboratory plasma physics such as in astrophysics plasmas, hydrodynamic fluids and molecular dynamics, atomic physics, condensed matter, chemical reactions and in statistics (see Refs. Scudder:1992 (); Kuhlen:2010 (); Leubner:2008 (); Leuthausser:1983 (); Riley:1975 ()). We focus here on laboratory plasma physics with charged particles in a magnetic field relevant to tokamaks and spherical tokamaks for the purpose of energy production by fusion energy. Some phenomena occurring in these devices break the symmetry of the Maxwellian distribution function (MDF) such as radiofrequency wave heating, neutral beam injection, ion orbit loss or simply boundary and external conditions (plasmasurface interaction, external magnetic field configuration, etc) (see Refs. Karney:1979 (); Jaeger:2006 (); Hirshman:1981 (); Dawson:1971 (); Jassby:1977 (); Shaing:1992 (); Battaglia:2013 (); Dahlin:2015 ()). Another phenomenon that is highly relevant to future fusion reactors is the selfheating by alpha particles produced by fusion reactions Matsuura:2009 (); Hay:2015 (). All these phenomena need to be considered in plasma physics in order to better reproduce the dynamics of current tokamaks and predict the confinement of future fusion reactors because to date, there is no rigorous and efficient theory describing nonMaxwellian distributions (NMDF) at finite collisionality. The ability to describe NMDFs is one of the most important unsolved problem in plasma physics. The number of studies about NMDF has been increasing in the last decades, based on some measurements of distribution function as well as on kinetic numerical simulations. Indeed, the Kappa distribution function (KDF) is usually observed in astrophysics Scudder:1992 (); Leubner:2008 (); Livadiotis_2014_Entropy () where different power laws in velocity phase space appear, but it is not usually observed in laboratory plasma physics. Moreover, bimodal distribution functions (i.e., the sum of two MDFs) are often used in order to better describe the presence of superthermal particles Popov:2009 (); Jaworski:2013 (), but we argue here that the sum of Maxwellians is valid only in two collisional regimes (i.e., when there is no collision or at infinite collisionality). The interpreted nonMaxwellian distribution function (INMDF) introduced here could be the first selfconsistent solution for the description of laboratory nonMaxwellians at finite collisionality. The initial purpose of the INMDF introduced here is to more effectively model already known effects. However, better understandings coming from the physical interpretation of the INMDF could help in the future to intuitively develop theories describing new physical phenomena. These three ways of describing nonMaxwellians can be generalized as needed by the readers in order to keep manageable velocity phase space integrals. Indeed, the analytic kinetic corrections obtained here are possible with these three nonMaxwellians and our developed method can be used with the creation of new distribution functions describing the phenomenon under investigation. In order to clarify to the readers the novelty of this work, we emphasize that in comparison to the literature where a large number of orthogonal basis functions are used for the numerical representation of NMDFs Peysson_2008_PoP (); Omotani_2013_PPCF (); Hirvijoki_2015_JPP (), we show the advantage to use nonorthogonal basis functions for the efficient analytic description of NMDFs. This allows at the same time the full analytic computation of the velocity phase space integrals with a small number of simple terms, as detailed in the appendices. Indeed, all required velocity phase space integrals which contain the introduced INMDFs are reduced to some special functions (e.g., the hypergeometric or the incomplete Gamma functions). However, because these functions are evaluated at specific values, we obtain for the first time all results as function of a small number of terms including simple polynomials and the exponential and Error functions. This method stays valid for a class of NMDFs obtained from our nonorthogonal basis sets. Our analytic description represents a real asset in plasma physics, leading to the resolution of the efficient unification between kinetic and fluid theories Izacard_2016_FischSymposium (); Izacard_2016_Paper2 (). In summary, this work is unique by its capability to analytically link a specific shape of NMDFs to a small number of fluid quantities observable in experiments. In order to represent the generality and universality of our work, we focus here on a selection of three relatively different analytic theories which are usually associated with MDFs in the literature:

As an example, there is a recurrent discrepancy between the transport simulated by fluid codes against its measurement in radiated detached divertor plasmas. The divertor plasma is a crucial component for the fusion energy because it is the main material component in direct contact with the plasma and can highly impact the cost efficiency due to the increase of required maintenance. The difference of transport and plasma profiles between experimental measurements and edge simulations seems to be explained by the presence of NMDFs Chankin_2007_NF_47a (); Chankin_2007_NF_47b (); Groth_2013_NF_53 (). The main issue of these edge simulations is the description of the radiation which is very sensitive to kinetic effects. One of the radiations that directly impact the edge (via the floating potential, the drag force on dust and impurity, or the neutral collision) is the secondary electron emission Bruining:1954 (); Sternglass_1957_PR (); Dionne:1975 (); MeyerVernet_1982_AA (); Mayer_1988_PRB (); Bacharis_2010_PoP_17 (); Raitses:2011 (); Raitses:2015 (); Shotorban_2015_PRE (). Recent works focused on the secondary electron emission by using a MDF and obtaining an empirical formula Bacharis_2010_PoP_17 (). We obtain here the analytic formulas for MDFs and NMDFs.

Another example is the discrepancy of the Langmuir probe interpretation Arslanbekov (); Sudit_1994_JAP (); Stangeby_1995_PPCF (); Batishchev_1997_PoP (); Demidov_2002_RSI (); Tskhahaya_2011_JNM (); Sheehan_2011_PoP (); Godyak_2011_JPD (); Jaworski_2012_FED (); Popov_2012_PSST (); Popov_2014_CPP (); Godyak_2015_JAP () with respect to the Thomson scattering measurements Luna_2003_RSI (); Beausang_2011_RSI () of the electron temperature which has been associated with the presence of superthermal particles Popov:2009 (); Jaworski:2013 () (i.e., bulk and superthermal populations, both at different thermodynamic equilibrium). The Langmuir probe is one of the most used diagnostic for low temperature plasmas. Advanced concepts allow the measurement of the electric potential, the electron density and temperature or the radial electric field. All of these measurements are possible by interpreting the signal and by using the assumption that the incident electrons are described by a MDF. Then, the presence of NMDF obviously results in discrepancies of the interpretation of physical quantities from the signal. We show here analytic results for MDFs and NMDFs in order to prepare future modified interpretations of experimental data.

Finally, the last example which is analytically based on MDFs is our understanding of the entropy at the thermodynamic equilibrium (i.e., the MDF). The entropy, viewed as the degree of disorder, is well known for a MDF. However, the experimental observations of NMDF steadystates has to be addressed in order to have a more general point of view of the entropy, similarly than in Refs. Kranys_1970_ARMA (); Livadiotis_2014_Entropy (). The analytic result shown here for an INMDF motivates the generalization of thermodynamics for nonisolated systems (extending Refs. Bizarro_2011_PRE (); Bizarro_2015_PoP ()) by using our INMDFs. Moreover, a significant new result shown here is the explicit decrease of the local entropy for some INMDFs.
Our work is drawn as follow. The first section describes the usual analytic representations (i.e., the KDF and the sum of MDF) that are well documented in the literature, and contains the foundations of new distribution functions (i.e., the INMDFs) relevant to some experimentally observed NMDFs. Then, this section focuses on the physical motivations of the introduced INMDF and its possible generalization. In Sec. II the kinetic corrections of the secondary electron emission, the Langmuir probe characteristic curve and the entropy are analytically predicted and shown with as relevant plasma physics parameters as possible. Finally, the conclusion is detailed in Sec. III.
I Choice of an analytic representation of nonMaxwellian distribution functions
This section deals with the velocity phase space representation of nonMaxwellian distribution functions (NMDFs) using specific analytic forms instead of the usual numerical discretization. In fact, the numerical approximation is not efficient for the description of some NMDFs because of the requirement to use a high number of terms. For example, it is not efficient to use a numerical approximation of a Maxwellian distribution function (MDF) in comparison to its analytic form which uses only 3 variables (density, fluid velocity and temperature). Indeed, this work argues the advantage of using a meshfree representation of the distribution function in the velocity phase space by using fluid parameters (called here the hidden variables) instead of the discretization in the velocity phase space. We adopt the notation of for the one dimension velocity phase space coordinate, for the three dimension position and for the time. The generalization to three dimensional velocity phase space is ongoing for isotropic distribution functions but some complexities appear for anisotropic ones Izacard_2016_Paper2 (). Further work will include anisotropic NMDFs. For simplicity of notations the division of the temperature by the mass is omitted in the distribution function.
In this section we show notations for existing distributions (the KDF and the sum of two MDFs) and we introduce a new class of NMDFs (i.e., the INMDF) which can be viewed as a nonorthogonal generalization of the sum of Maxwellians and Hermite polynomials.
i.1 The Kappa distributions
The KDFs are commonly used in astrophysics Scudder:1992 (); Leubner:2008 () to describe different tail power laws than the Maxwellian. The definition of the KDF is
(1) 
where
(2)  
(3) 
and is the Euler Gamma function. The fluid quantities, which are evolving in time and are functions of the position, are , (or ) and . For a shifted distribution function the fluid velocity appears by the change of variable in Eqs. (1), (4) and (5). The KDF recovers the MDF for . Many properties of the KDF have been published in the literature and are not reproduced here. Integrals of KDFs over the velocity phase space introduce functions detailed in App. A
(4)  
(5) 
By using these integrals, it is trivial to obtain all fluid moments of any shifted KDF centered at the fluid velocity as function of a finite sum of the function . The KDF is successfully used to describe different power laws in velocity phase space. However, because an over population of fast particles localized around a specific velocity is often observed, the sum of two Maxwellian can be a better approximation for these cases.
i.2 The sum of Maxwellians
The approximation of localized superthermal tails is commonly used Popov:2009 (); Jaworski:2013 (); Jaworski_2015_Private (). It is usually called the bimodal distribution function where two different temperatures dominate as function of the velocity coordinate . The sum of two shifted MDFs reads
(6)  
with
(7)  
(8) 
where the density, fluid velocity and temperature are respectively for the bulk and the superthermal (i.e., fast particles) populations , , and , , . The velocity integrals are directly given by the function detailed in App. B and defined by
(9) 
Another useful function is detailed in App. C and App. D and defined by
(10) 
The generalization with a sum of more than two MDFs is possible with the Radial Gaussian Basis Function (RGBF) which has successfully been used for artificial neural networks. Many properties of the bimodal distribution function and the RGBF have been published in the literature and are not reproduced here.
This representation is accurate for example when there is an external production of fast particles and when there is no interaction between the bulk plasma (of collisionality ) and these fast particles (of collisionality ). This means that the RGBF is valid only at two collisional regimes: the limit of no collision between thermalized and fast particles (i.e., ) or the limit of infinite collisionality (i.e., ). In other words, the RGBF is not consistent with interactions between the bulk plasma and the fast particles at finite collisionality because each of them become nonMaxwellian. The novel representation introduced below is a possible solution of this inconsistency.
i.3 The interpreted nonMaxwellian
The definition of the INMDF was introduced for the first time in one dimension by Refs. Izacard_2013_GA_Talk (); Izacard_2013_Emails_Candy (); Izacard_2016_FischSymposium ()
(11) 
with
(12)  
(13) 
and with given by Eq. (7) and
(14) 
where the fluid hidden variables are function of . The first part is a MDF and due to the second part , the INMDF cannot be described with a finite number of terms using any of the existing analytic representations. The second part is related to a physical interpretation and could be valid for all collisionality regimes because represents an enhancement of the energy of a population of particles. This new formula represents a real asset for the analytic modeling of NMDFs. More details of the difference between the INMDF and all existing analytic NMDFs are given in Sec. I.4. The moments of the nonMaxwellian are
(15)  
(16)  
(17)  
(18)  
(19)  
(20)  
(21)  
and generally given as function of . The usual moments for of the nonMaxwellian are also function of the hidden variables. We do not write these moments since it is not used here, but if needed, readers can easily obtain them from the given moments.
Nevertheless, we remark as expected for a MDF that when the odd moments are equal to 0 and the even moments are .
The coefficient has the dimension of a momentum (a density times a velocity) then has the dimension of a velocity, of a velocity, and of a temperature (square of a velocity, including the division of the temperature by the mass which is omitted by simplicity of notations).
We remark that, by using the method of hidden variables, there is no reason to resolve the inversion of the Eqs. (15)(21) in order to extract the hidden variables as function of the fluid moments. To date, the inversion of these formulas to extract the hidden variables as function of the moments is possible only for a MDF and it seems too constraining to use this criteria for the description of nonMaxwellian steadystate distribution functions observed experimentally.
Fig. 1.(a) represents an INMDF (solid blue curve) with a correction (dashed red curve) centered at the mean velocity with respect to the MDF (dashed black curve). This example corresponds to an asymmetric distribution function observed at least in presence of ion orbit losses Battaglia:2013 () and will be investigated in future work. Fig. 1.(b) represents the same correction but centered at a higher velocity coordinate (i.e., higher value of ).
The tail at high energy (centered at ) is unstable when the distribution function increases and is not monotonous. In another words, the distribution function is unstable when its first derivative is positive around . Indeed, the stability criteria of the INMDF is obtained when the derivative of the Maxwellian part is bigger than the derivative of the additional part (i.e. ). We found the stability criteria
(22) 
then, for the highest value at gives for the stability criteria
(23) 
For , similar stability criterion can be obtained for a range of velocity but are not detailed here. Moreover, the constraint to consider a strictly positive distribution function is given by the relation for all , particularly when, for , the additional nonMaxwellian part is minimum (for ). The minimum value of the additional nonMaxwellian part is obtained at
(24)  
Then the constraint of a strictly positive distribution function is given by , so
(25) 
We remark that for , similar relations can be obtained using the minimum value of obtained for (i.e., at ).
Before investigating corrections of some existing theories due to nonthermal population of particles, we have to understand how the previously given formulas could accurately describe deviations from MDFs. For that, the physical interpretation of the INMDF is given here.
The fluid quantity is called here the kinetic flux, the central flow and the width of the heat spread. Similar to the graphical interpretation of , and for the MDF, the graphical interpretation of the nonMaxwellian part , and is shown in Fig. 2.
Fig. 2 represents the displacement of a population of particles () from the darkgreen area to the lightgreen area. The hidden variable is called the kinetic flux since it happens in the velocity phase space and has a dimension of a particle flux (i.e., a density times a velocity). The central flow is the velocity where the INMDF equals the MDF. Finally, the width of the heat spread characterizes the width in velocity phase space of the superthermal population modified by an external source of energy (e.g., the current drive). In the following figures, the coefficients are used to parameterize the superthermal tail such that where represents approximately the percentage of superthermal particles over the total number of particles, represents the position of the central flow by the relation , and represents the ratio of the width of the heat spread over the temperature (i.e., is the ratio of the widths).
i.4 Generalization with other INMDFs
It is possible to use other formulas for the description of nonMaxwellians such as shown in
Fig. 3 where different NMDFs are obtained from
(26) 
with , where is the floor function and , , , , are fluid coefficients for all integer . The set of coefficients are the hidden variables. We notice that has a dimension of a density if is even and of a particle flux if is odd, and have the dimension of velocities and and of temperatures. This general form can describe asymmetric distribution functions which deviate from MDFs and the effective density can be different to (e.g., when we consider odd values of ). Moreover, because the hidden variables and as well as and can be different, in opposite to the Fig. 1, we can consider asymmetric corrections as shown by red curves in Fig. 4.
Further investigations using this generalized INMDFs given by Eq. (26) are possible. In summary, this formulation (i.e., the generalized INMDFs) can be related to a generalization of the Hermite polynomial representation Grad:1949a (); Grad:1949b (); Grad:1963 () (due to the presence of term) but each term is associated with a different MDF. The representation given by Eq. (26) forms a nonorthogonal basis in opposition to usual work found in the literature where orthogonal basis are used in order to project the distribution function on the basis set, but here this nonorthogonal basis is relevant to a significant reduction of the number of required hidden variables. Moreover, the completeness property is intrinsically inherited from the usual Hermite polynomial since Eq. (26) can be reduced to an Hermite polynomial with specific relation between all hidden variables (similar argument with a possible reduction to the RGBF, see below). The main advantage of this representation against the Hermite polynomial is that the required number of terms can be drastically reduced for the description of nonMaxwellians with high flows.
Some researchers argue that one needs to keep at least hundreds or thousands of terms with the usual Hermite polynomial representation (so at least hundreds or thousands independent fluid moments, see Ref. Izacard_2016_Paper2 ()) for an accurate description of nonMaxwellian plasmas. With our generalization given by Eq. (26), it is conceivable to keep much less hidden variables (related to much less independent fluid moments). Having less hidden variables is an argument of the advantage of the INMDF. Another argument is the physical interpretation (see the end of Sec. I.3) of the hidden variables. We remark that our general representation can also reproduce the sum of shifted MDFs (i.e., the RGBF) by imposing other constraints such as and .
The moments of the generalized distribution function given by Eq. (26) become
(27) 
where , is the combination and the function is introduced above. The introduced INMDF given by Eq. (11) or Eq. (26) is a new minimal formulation for a local nonMaxwellian deviation. This distribution function cannot be exactly obtained with a finite number of terms by any of the existing basis functions. In comparison to the sum of MDFs (i.e., the RGBF Hirvijoki:2015 ()), valid only at two collisionality limits, the INMDF is the first distribution function which can be selfconsistent with MDFs steady states at finite collisionality. More properties of the INMDF will be published later.
i.5 Physical reality of the INMDFs
The readers may be concerned about the physical reality of the newly introduced INMDFs with respect to the well observed KDF in astrophysics or the natural expansion with a sum of MDFs. In order to argue about the physical picture introduced by the displacement of population of particles in the velocity phase space of the INMDFs, we describe below a qualitative fitting of an indirect experimental measurement of the distribution function. In fact, to date there is no universal direct measurement of the velocity phase space variations of the distribution function but it is common to indirectly compute the distribution function from different diagnostics. Moreover, even if we choose to perform the qualitative fitting (without the necessity to use the experimental data) of only one published article on the electron temperature measurement discrepancy, it is at least possible (not shown here, see Ref. Izacard_2016_Paper2 ()) to fit NMDFs observed by FokkerPlanck numerical codes in presence of lower hybrid current drive or by particle in cell codes in presence of ion orbit losses.
Fig. 5 proves the capability of INMDFs to describe experimentally observed physical processes. With some density and temperature constraints in presence of NBI and ICRF in JET Luna_2003_RSI (); Beausang_2011_RSI () and TFTR Taylor_1996_PPCF (); Fidone_1996_PoP (); White_2012_NF (), it has been observed a discrepancy of the electron temperature between the interpretation from the electron cyclotron emission (ECE) and Thomson scattering (TS) due to kinetic effects Beausang_2011_RSI (). De La Luna, et al. and Beausang, et al. numerically found a model NMDF (dashed red curve) constructed from the spectrum of the ECE. In fact, the measurement of the ECE spectrum was not consistent with a MDF. De La Luna, et al. originally found the way to better recover this ECE spectrum by numerically modifying the distribution function and obtaining the called model NMDF. It turned out that this model NMDF resolves the discrepancy of the interpreted electron temperature by TS. The physical processes at the origin of this model NMDF are not understood. Our three over plots (blue, green and cyan dashed curves) in Fig. 5 can help to understand the origin of the detected NMDF. The blue curve corresponds to a MDF with a specific set of parameters. The green curve is obtained by adding with the parameters . Then, the cyan dashed curve is obtained by adding with . This means that we detect both, an enhancement (i.e., a heating ) and a reduction (i.e., a cooling ) of the energy of two different populations of particles. We notice that because we use interpreted parameters, we do not need to access the data (i.e., the density, fluid velocity and temperature of the background MDF). Our additional parameters () are the only way to modify the shape of the distribution function and the use of the terms () make our figure independent of the data. We prove here that the TSECE discrepancy is due to approximately of heated particles and of cooled particles, and both are attracted by a resonantlike process around . We highlight the fact that this result may not be the detection of a new physical effect but, at least, it helps to describe and understand the unsolved problem of the TSECE discrepancy. Future work could try to link our understandings with known heating and current drive resonant theories or may lead to the discovery of a new physical effect. FokkerPlanck codes may help to fully understand the details of the energy transfer. In summary, we do not attempt to describe here what is resonant with what in this complex unsolved problem where an electron NMDF, associated here with a kind of current drive resonance, is experimentally observed in presence of ion heating and current drive (NBI and ICRH).
However, even if we do not know yet the exact origin of these detected physical processes (heating and cooling), the analytic formula of the INMDF can be directly used to understand the effects of nonMaxwellian bulks on other diagnostics and on the transport and turbulence. Indeed, even without theoretically or numerically describing the physical processes at the origin of the heating and cooling down we observed here, we can still have a good description of the deviations of the NMDF with respect to a MDF and analytically predict and validate the kinetic effects from this observed NMDF on different diagnostics. Our capability to calibrate different experimental diagnostics is indispensable since we prove here that less than of superthermal particles can lead to a discrepancy of the electron temperature interpretation between TS and ECE when a MDF is assumed. It is specially indispensable for the prediction of the effects of superthermal tails on different diagnostics in ITER where a of selfheating by fusion reactions is planned in addition to some external heating. Future work will focus on our detected INMDF and will help to initiate standard experimental procedures to measure NMDFs and to resolve inconsistencies between diagnostics more or less sensitives to kinetic effects.
In the next section, the advantage of using one of these NMDFs is shown for the analytic computation of kinetic corrections on the secondary electron emission, the Langmuir probe characteristic curve and the entropy.
Ii Kinetic corrections due to NonMaxwellians
This section applies the NMDFs shown in the previous section on some existing theories relevant to current and future plasma devices. Analytic corrections of these theories are given as function of the fluid hidden variables. These corrections can be implemented in existing numerical simulations in order to better describe some kinetic effects when they need to be taken into account.
ii.1 Corrections for the secondary electron emission
In tokamaks scrapeofflayer (SOL), incident primary particles (i.e., charged or neutral particles) release secondary electron from a solid surface. The secondary electron emission (SEE) reduces the sheath potential and can modify many phenomena in the plasmamaterial interface Bruining:1954 (). The computation of the SEE is obtained by formulating the SEE yield for incoming particles of energy related to the material under investigation (Sternglass formula Bruining:1954 (); Sternglass_1957_PR (); Dionne:1975 ()) multiplied by the proportion of particles at this energy (i.e., the distribution function ) and has been under investigation in the literature MeyerVernet_1982_AA (); Mayer_1988_PRB (); Bacharis_2010_PoP_17 (); Raitses:2011 (); Shotorban_2015_PRE (). By integrating over the possible energy of incoming particles, the secondary electron emission with respect to the incoming particles distribution function becomes
(28) 
with , the constant coefficient , and is the kinetic energy. Some values of and have been obtained experimentally and are dependent of the material under consideration (e.g., see Refs. Bacharis_2010_PoP_17 (); Dionne:1975 ()). The value of corresponds to the energy associated to the value of , the maximum SEE measured and used to obtain a normalized . The empirical formula for a MDF obtained in Ref. Bacharis_2010_PoP_17 () is
(29) 
with , the coefficient are given adhoc for each species, and the dependence with respect to the fluid velocity disappears by using the assumption of a Maxwellian plasma with no mean flow. However, we can reformulate Eq. (28) to
(30) 
with , and .
Here we found the following exact analytic formula from Eq. (30) for a MDF given by Eq. (12)
(31) 
with and where the general form and some functions are detailed in App. C where
(32) 
Following the same analytic method, it has been possible to obtain the corrections from a KDF. However, we do not focus on this result here because it involves the hypergeometric PFQ function instead of due to of the presence of in the term and it may not be possible to obtain a simple analytic form in comparison to the following results. This may mean that the KDF is less natural than the sum of two MDFs or the INMDF (see results below) due to the complexity of the solution of the SEE . We found for the sum of two Maxwellians the following corrections
(33) 
and for the INMDF given by Eq. (11) the following corrections
(34) 
where the secondary electron emissions , and are functions of all fluid hidden variables and the quantities , , and are previously defined.
Fig. 6.(a) shows different NMDFs computed from the bimodal NMDF given by Eq. (6) and from the INMDF given by Eq. (11) for a temperature . The corrections of the secondary electron emission as function of the temperature are shown in Fig. 6.(b).
An important result is that the black curve obtained for a MDF with our exact formula does not match exactly the red curve of the empirical formula given in Ref. Bacharis_2010_PoP_17 (). The reason is due to an undetected presence of superthermal particles in Ref. Bacharis_2010_PoP_17 () because our analytic formula with the INMDF (cyan dashed curve) matches better the empirical formula (red curve) better in opposition to our analytic formula with the MDF (black curve). Moreover, the empirical formula is better matched for with (blue dashed curve) and for with (green dashed curve). We remark that only for the cyan curve in Fig. 6.(b) we used instead of , but it is possible to find other parameters which match the red empirical curve (e.g., keeping and multiplying and by 0.95 because is linear with respect to these two quantities). Finally, it seems that the superthermal particles of the empirical formula are clearly better described by the INMDF rather than by a sum of two MDFs (yellow solid curve is one of the best fit). Nevertheless, the goal is to predict the modifications of the secondary electron emission in presence of different superthermal particles. The effects of a tail on the secondary electron emission are significant specially for a low temperature (i.e., ) of the background plasma since a factor around can be observed (i.e., not shown here). Then, it is indispensable to take into account superthermal particles even if this population represents less than because, at least, the dynamics of dusts, neutrals and impurities are highly impacted by the secondary electron emission. It looks like our analytic prediction (with ) is able to observe two peaks in the SEE like it has been experimentally observed (see Ref. Raitses:2015 ()). Future work will investigate this two peaks observation, particularly the addition of other terms in Eq. (11), in order to avoid negative distribution function for all .
ii.2 Corrections for the Langmuir probes interpretations
A commonly used diagnostic of edge (cold) plasma is the Langmuir probe. By applying an electric potential scan on the Langmuir probe, we can measure the current of the probe as function of the applied potential. This give us the well known characteristic curve. The theory of the Langmuir probe has successfully been developed in the past century Stangeby_2000 (); Hutchinson_2002 () and properties of the plasma can be extracted from the characteristic curves such as the electron temperature, or the floating and plasma potentials by interpreting the measurement for a background Maxwellian plasma. However, a recurrent discrepancy of the interpretation of the plasma properties with other diagnostics (e.g., with Thomson scattering measurements Luna_2003_RSI (); Beausang_2011_RSI (); Jaworski_2015_Private () in attached plasma) has been observed. This discrepancy has been related to NMDFs as detailed in Refs. Popov:2009 (); Jaworski_2012_FED (); Jaworski:2013 () where bimodal distribution functions (distribution function of populations of particles at different temperatures) have been used. Other results are published in the literature describing effects of nonMaxwellians on the Langmuir probe measurements Arslanbekov (); Sudit_1994_JAP (); Stangeby_1995_PPCF (); Batishchev_1997_PoP (); Demidov_2002_RSI (); Tskhahaya_2011_JNM (); Sheehan_2011_PoP (); Godyak_2011_JPD (); Jaworski_2012_FED (); Popov_2012_PSST (); Popov_2014_CPP (); Godyak_2015_JAP (). The bimodal approximation is the first efficient way to describe NMDFs, but it assumes the superposition of two populations of particles at the Maxwellian equilibrium and the selfconsistency is very restricted since in this bimodal description the plasma is enough collisional to assure a MDF for each species but do not allow collisions between these two populations. In order to resolve this inconsistency, we propose to use the INMDF of the plasma (only one population of particles with effective temperature) which contains superthermal particles. In our case, all particles (viewed as one species) collide together and external sources or sinks of energy create nonMaxwellian steady states.
The general formula of the electron current Arslanbekov (); Demidov_2002_RSI (); Jaworski_2012_FED (); Popov_2014_CPP () can be written as a function of any distribution function of the kinetic energy of particles
(35) 
where is a diffusion parameter, a geometric parameter of the probe, is the difference between the applied potential at the probe and the plasma potential, the surface of the probe, and is the mass of electron. The classical regime is obtained by assuming a diffusionless limit (), and with for a spherical probes. This formula computes the flux of electron which have higher relative energy than the potential because electrons of lower energy cannot contribute to the current due to the Coulomb barrier generated by other electrons at the probe interface. With the change of coordinate and , the electron current becomes
(36) 
If we assume a background MDF given by Eq. (12), the electron current becomes
(37) 
with the definition of the function
(38) 
detailed in App. D. Then, from the KDF given by Eq. (1) the electron current becomes
(39) 
from the sum of two MDFs given by Eq. (6) the electron current becomes
(40)  
and from the INMDF given by Eq. (11) the electron current becomes
(41)  
As shown in Fig. 7, different NMDFs (see Fig. 7.(a)) with the same Maxwellian background (black curve) can significantly modify (see Fig. 7.(b)) the electric current measured by a Langmuir probe. This figure is obtained for parameters relevant to the SOL plasmas close to the divertor plates (i.e., , ). However, we remark that the amplitudes are not consistent with the experimental measurements by the Langmuir probes since the formula given by Eq. (35) (see Refs. Arslanbekov (); Demidov_2002_RSI (); Jaworski_2012_FED (); Popov_2014_CPP ()) is not dimensionally correct. Usually, results are given in arbitrary units in the literature Jaworski_2015_Private () in order to avoid this dimensional inconsistency. This dimensional issue can be investigated later since it is not the focus here.
Moreover, Eq. (36) (directly linked to the Druyvesteyn formula Arslanbekov (); Demidov_2002_RSI (); Popov_2012_PSST (); Popov_2014_CPP ()) is the diffusionless limit of Eq. (35). In comparison to other studies Jaworski_2015_Private (), the INMDF used in the diffusionless limit Eq. (36) can reproduce similar curves than the use of a MDF in the diffusional Eq. (35). This observation is not surprising because the particles diffusion describes an enhancement of the displacement of particles as well as the population of superthermal particles. This is a very important observation because instead of using adhoc diffusion parameter in order to reproduce experimental observations, the existence of NMDFs could replace those.
Finally, all common interpretations of the plasma parameters are modified by the NMDF. For example, the floating potential is defined by where is the current of the Langmuir probe generated by ions and commonly assumed to be constant due to ion mass and temperature involved. From the result given by Eq. (41) and its version for the ion current, it is possible to better evaluate the floating potential from the distribution function. Of course, all other quantities such as the plasma potential , the electron temperature and the radial electric field are interpreted differently due to NMDFs. Additional investigations will be reported later.
ii.3 Entropy decrease
In the thermodynamics theory, the entropy is interpreted as the degree of disorder of a system. A generalization of the entropy is required for the description of non thermodynamics equilibrium. By sharing similar perspectives than Refs. Kranys_1970_ARMA (); Livadiotis_2014_Entropy (), corrections of the entropy can be analytically computed from a NMDF. In an isolated plasma, the entropy can only increase due to collisions in order to reach its maximum value when the distribution function is Maxwellian. The definition of the statistical entropy developed by Boltzmann is
(42) 
where is the distribution function of the particles of density and is the Boltzmann constant. The maximal value of the entropy is obtained with a MDF (see App. F)
(43)  
(44) 
such that for isolated equilibrium plasmas. For other distribution functions than the MDF, the entropy can only be smaller to this value and increase (). However, in order to make a link with observed NMDFs for nonisolated plasmas, it makes sense to interpret the entropy by a level of sharing energy or information. The maximal entropy obtained with a MDF of an isolated plasma is reached when the interaction between all particles of the plasmas have shared the total energy of the system after a thermalization characteristic time . This level of sharing can decrease when an external source or sink of energy of a nonisolated plasma appears locally in the velocity phase space at a finite collisionality. The presence of a steadystate superthermal population of particles due to external sources introduces corrections of the entropy because at finite collisionality some particles can be sensitive to the external source of energy in time range , shorter than the one associated to the sharing of energy. The analytic computation of the corrections becomes possible using the INMDF written in the form with the MDF and the corrections previously defined (respectively by Eqs. (12) and (13)). The modified entropy due to the presence of a superthermal population reads
(45) 
where . Using the Taylor expansion of the logarithm of around (i.e., ) if we assume relatively small deviations from the Maxwellian (i.e., ), then, following the details of App. F, the entropy becomes with
(46) 
At the first order in the expansion of (i.e., the sum is of higher order, see App. F) and with the use of the following relations
(47)  
(48)  
(49)  
(50) 
when , and with the functions given in App. B, the correction of the entropy is
(51)  
(52) 
In comparison to Refs. Bizarro_2011_PRE (); Bizarro_2015_PoP () where the entropy can decrease in presence of frictions, here we assume an INMDF steadystate without friction. We found that the first order correction of the entropy can either be positive or negative as function of the hidden variables describing the INMDF. We remark that our INMDF is not related to an isolated system, so the inclusion of the collision and the source of energy (e.g., radiofrequency waves, neutral beam, runaway electron) would recover known thermodynamic results of isolated systems (e.g., similar to Ref. Bizarro_2015_PoP ()). This means that when an external source or sink of energy is turned on, the corrected entropy of the plasma () computed from an INMDF can increase or decrease with respect to the entropy obtained from a MDF. Both cases are possible without violating the second law of the thermodynamics since here the plasma is nonisolated and the total entropy increase due to the external source is not taken into account. In fact, instead of describing a complete theory of the source and the plasma, we rather focus on the statistical description of the plasma in presence of sources. Then, for a nonisolated system, the maximum value of the entropy is a higher value than as function of the plasma parameters. Moreover, since both variations of are possible, the intermediate case is the solution where the superthermal particles do not modify globally the entropy of the Maxwellian at the first order. There are solutions of . The first, for finite values of the density and temperature when , recovers the MDF. The second solution is obtained for and when
(53) 
The physical interpretation of this specific value of the width of the heat spread still needs to be understood and is still under investigation. However, with this specific INMDF profile, the statistical entropy of the plasma, i.e., without including the entropy increase due to the external sources, is locally constant. This means that the increase of the energy (which increases the entropy) is compensated exactly with the departure of the INMDF with respect to the MDF (which decreases the entropy). The time derivative of the entropy is found by using the Boltzmann equation and omitting the collision operator. It reads where . For the INMDF, with , and . Using the same expansion of the logarithm, we found that the time evolution of the entropy is modified at the first order by the spatial derivatives of
(54)  
(55) 
These relations can be used in order to compute the localized statistical entropy evolution of the plasma as a fluid quantity. As a numerical error criteria, for a MDF and without sources, we have to verify that where .
The increase or decrease of the entropy in time is given here as function of the spatial gradients of the the hidden variables. A very important interpretation of this result is following:
In many experiments (even for other areas than plasmas) nonMaxwellian steadystate distribution functions are observed. This means that in order for us to approximate this steadystate with the proposed INMDF or by creating new analytic NMDFs with as few hidden variables as possible, we have to verify that the local entropy of a nonisolated plasma (including sources) is constant in time because the entropy increase due to the sources is compensated by the entropy reduction of the NMDF steadystate with respect to the maximum entropy of a heated MDF. Then, there cannot be evolving local entropy at NMDF steadystate. From our result of written as function of the hidden variables, there are an infinite number of solutions which locally conserve the entropy at least at the first order in the expansion of small nonMaxwellian deviations. For example, we obtain one class of these solutions by trivially extracting the density as function of the other fluid hidden variables since the equation contains only one occurrence of . We can write this equation as with and which are both functions of . This class of solutions of locally constant entropy is obtained when , if (i.e., if and ). We have an infinite number of solutions because we have parameters to resolve the equation