The equivalent in special relativity (SR) of the Maxwell-Boltzmann (MB) distribution, see  , is the so called Maxwell-Jüttner distribution (MJ), see  . The MJ distribution has been recently revisited. We select some approaches among others: a model for the anisotropic MJ distribution , an astrophysical application of the MJ distribution to the energy distribution in radio jets , a new family of MJ distributions characterized by the parameter  and an application to counter-streaming beams of charged particles . The above approaches do not cover the determination of the statistical quantities of the MJ distribution. In this paper, the statistical parameters of the relativistic MB distribution are derived in Section 2 and those of the MJ distribution are derived in Section 3. Section 4 derives the spectral synchrotron emissivity in the framework of the two relativistic distributions here analyzed.
2. The Relativistic MB Distribution
The usual MB distribution, , for an ideal gas is
where m is the mass of the gas molecules, k is the Boltzmann constant and is the usual thermodynamic temperature. In SR, the total energy of a particle is
where m is the rest mass, c is the light velocity, is the Lorentz factor , and v is the velocity. The relativistic kinetic energy, , is
where the rest energy has been subtracted from the total energy, see formula (23.1) in . A relativistic MB distribution can be obtained from Equation (1) replacing the classical kinetic energy with the relativistic kinetic energy
where the relativistic temperature; T, is expressed in units; up to now the treatment is the same of  at pag. 665. The above relativistic PDF nolistsep
· has the velocity of the light as maximum velocity,
· becomes the usual MB distribution in the limit of low velocities,
· is not invariant for relativistic transformations.
2.1. Variable Lorentz Factor
We now change the variable of integration
The differential of the velocity, ,
and therefore the relativistic MB distribution in the variable is
where G is the Mejier G-function   ; Figure 1 reports the above PDF for three different temperatures.
The average value or mean, , is
the second moment about the origin is
the variance, is
The mode is the real solution of the following cubic equation in
which has the real solution
Figure 1. The PDF of the relativistic MB as function of for different values of T.
At the moment of writing a closed form for the distribution function (DF) which is
does not exists and we therefore present a numerical integration, see Figure 2.
The asymptotic behavior of the PDF, , is
The integration of the above approximate PDF gives an approximate DF which has a maximum percentage error of 7% in the interval when . The random numbers belonging to the relativistic MB can be generated through a numerical computation of the inverse function following the algorithm outlined in Sec. 4.9.1 of . The above PDF has only one parameter which can be derived approximating the average value with a Pade approximant
The above approximation in the interval has a percent error less than 1%. The inverse function allows to derive T as
Figure 2. The numerical DF of the relativistic MB as function of for different values of T.
Here is the sample mean defined as
formula which is useful to derive the variance of the sample
where are the n-data, see . An example of random generation of points is reported in Figure 3 where we imposed and we found from the generated random sample.
2.2. Variable Velocity
We now return to the variable velocity, the PDF is
where v is expressed in units. The mode is a solution of a sextic equation, see , in v
which has the following real solution
Figure 3. The DF for 3000 random points generated according to the relativistic MB (blue steps diagram) and the theoretical DF (red dashed line) see 3.
The position of the mode for the PDF in v is different from that one in , see Figure 4.
At the moment of writing the other statistical parameters cannot be presented in a closed form.
3. The Maxwell Jüttner Distribution
The PDF for the Maxwell Jüttner (MJ) distribution is
where , m is the mass of the gas molecules, k is the Boltzmann constant, is the usual thermodynamic temperature and is the Bessel function of second kind, see    . Figure 5 reports the above PDF for three different values of and Figure 6 displays the PDF as a 2-D contour.
The average value is
and the variance is
Figure 4. The position of the mode in the two PDFs: gamma variable (red full line) and v variable (dashed green line).
Figure 5. The PDF of the MJ as function of for different values of .
Figure 6. A 2-D contour for the MJ PDF as function of and .
The mode can be found by solving the following cubic equation
The real solution is
The asymptotic expansion of order 10 for the PDF is
The DF is evaluated with the following integral
which cannot be expressed in terms of special functions.
We now present some approximations for the distribution function A first approximation is given by a series expansion when, ad example,
which has a percent error less < 0.6% in interval when . A second approximation is given by an asymptotic expansion of order 50 for the PDF followed by the integration, see Figure 7. The parameter can be derived from the experimental sample once the average value is modeled by a Pade approximant and the inverse function is derived
An analogous formula allows to derive from the variance of the sample
Figure 7. The numerical MJ DF as function of for different values of .
An example of random generation of points is reported in Figure 8 where we imposed and we found from formula (30) and from formula (31).
We now change the variable of integration in , the PDF of the MJ is
where , see Figure 9. We have only one analytical result, the mode, which is found solving the following equation in
As an example when the mode is at and Figure 10 reports the mode as function of .
The DF of the MJ is given by the following integral
Figure 8. The PDF for 70,000 random points generated according to the MJ (blue steps) and the theoretical PDF (red dashed line), see formula (22).
Figure 9. The 3D surface of the MJ PDF as function of and .
Figure 10. The mode of the MJ distribution for different values of .
Figure 11. The mean of the MJ distribution for different values of .
Figure 12. The variance of the MJ distribution for different values of .
Figure 13. The numerical DF of the MJ distribution (red circles) and the Riemann approximation (green line) as function of .
with in [0, 1] which does not have an analytical expression. An approximation is given by the Riemann sums, see , when
see Figure 13. The above DF has a maximum percentage error of ≈10% at .
4. The Astrophysical Applications
This section reviews the synchrotron emissivity for a single relativistic electron, derives the spectral synchrotron emissivity for the two relativistic distributions here analyzed and models the observed synchrotron emission in some astrophysical sources.
4.1. Synchrotron Emissivity
The synchrotron emissivity of a single electron is
where, according to eqn.(8.58) in , e is the electron charge, B is the magnetic field, is the pitch angle, is the permittivity of free space, c is the light velocity, is the electron mass, x = is the ratio of the angular frequency ( ) to the critical angular frequency ( ) and
where is the modified Bessel function of second kind with order 5/3  . The modified Bessel function is also known as Basset function, modified Bessel function of the third kind or Macdonald function see pag. 527 in . The above function has the following analytical expression
where is a regularized hypergeometric function    . Figure 14 displays as function of x.
4.2. The Synchrotron Relativistic MB Distribution
We start from the PDF for the relativistic MB distribution as represented by Equation (7) and we perform the following first change of variable
Figure 14. F(x) as function of x with logarithmic axes.
where E is the relativistic energy. The resulting PDF in relativistic energy is
A second change of variable is
We know that where B is the magnetic field expressed in gauss and therefore the above PDF in frequency becomes
4.3. The Synchrotron Maxwell Jüttner Distribution
We start from the PDF for the Maxwell Jüttner distribution as given by Equation (22) and we perform two changes in variable as in the previous section. The resulting PDF in relativistic energy is
The second PDF in is
The astrophysical PDF in frequency for the Maxwell Jüttner distribution is
The mismatch between measured flux in Jy and theoretical flux, , can be obtained introducing a multiplicative constant C
4.4. The Spectrum of the Radio-Sources
As a first example we analyze the spectrum of an extended region around M87, see as example Figure 1 in  where the flux in Jy as function of the frequency is reported in the range . Figure 15 reports the measured and theoretical flux in the range for the quiet core of M87.
A second example is given by the radio sources with ultra steep spectra (USS) which are characterized by a spectral index, , lower than −1.30 when the radio flux, S, is proportional to , see . As a practical example we select the cluster Abell 1914 where the measured total flux densities at 150 MHz and 1.4 GHz are and which means . We now evaluate the theoretical spectral index of synchrotron emission for the relativistic MB distribution between 150 MHz and 1.4 GHz when B is fixed and T variable, see Figure 16 and Figure 17 when T and B are both variables. The two
Figure 15. Measured flux in Jy (green stars) of the quiet core of M87 and the theoretical flux (red line) for the Maxwell Jüttner distribution as given by Equation (48). The numerical parameters are , and .
Figure 16. The spectral index of the relativistic MB as function of T when .
Figure 17. The spectral index of the relativistic MB as function of T and B in gauss.
Figures above show that the theoretical spectral index is always smaller than −2 which can be considered as an asymptotic limit for high values of relativistic temperature. As an example when the spectral index is −2.17 when .
The relativistic MB distribution has been derived in  without any particular statistics: here we derived, when the main variable is the Lorentz factor , the constant of normalization, the average value, the second moment about the origin, the variance, the mode, the asymptotic behavior, an approximate expression for the average value as function of the temperature and an inverted expression for the temperature as function of average value.
We derived the following statistical parameters of the MJ distribution when is the main variable: average value, variance, mode, asymptotic expansion, two approximate expressions for the distribution function, a first evaluation of from the average value and a second evaluation of from the variance.
Following the usual argument which suggests a power law behavior for the spectral distribution of the synchrotron emission in presence of a power law distribution for the energy of the electrons, we derived the spectral distribution for the relativistic MB and MJ distributions which are now function of the selected generalized temperature and the magnetic field. Two astrophysical applications are given: the spectral distribution of emission in the core of M87 in the framework of the synchrotron emissivity and an explanation for the steep spectra sources in the framework of the synchrotron emissivity for the relativistic MJ distribution.
 Maxwell, J.C. (1860) V. Illustrations of the Dynamical Theory of Gases.—Part I. On the Motions and Collisions of Perfectly Elastic Spheres. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 19, 19-32.
 Aragón-Muñoz, L. and Chacón-Acosta, G. (2018) Modified Relativistic Jüttner-Like Distribution Functions with η-Parameter. Journal of Physics: Conference Series, 1030, Article ID: 012004.
 Sadegzadeh, S. and Mousavi, A. (2018) Maxwell-Jüttner Distributed Counterstreaming Magnetoplasmas—Parallel Propagation. Physics of Plasmas, 25, Article ID: 112107.
 Oldham, K.B., Myland, J. and Spanier, J. (2010) An Atlas of Functions: With Equator, the Atlas Function Calculator. Springer Science & Business Media, New York.
 Prieto, M.A., Fernández-Ontiveros, J.A., Markoff, S., Espada, D. and González-Martn, O. (2016) The Central Parsecs of M87: Jet Emission and an Elusive Accretion Disc. MNRAS, 457, 3801.
 De Breuck, C., van Breugel, W., Röttgering, H.J.A. and Miley, G. (2000) A Sample of 669 Ultra Steep Spectrum Radio Sources to Find High Redshift Radio Galaxies. Astronomy and Astrophysics Supplement Series, 143, 303.