New Probability Distributions in Astrophysics: IV. The Relativistic Maxwell-Boltzmann Distribution ()
1. Introduction
The equivalent in special relativity (SR) of the Maxwell-Boltzmann (MB) distribution, see [1] [2], is the so called Maxwell-Jüttner distribution (MJ), see [3] [4]. The MJ distribution has been recently revisited. We select some approaches among others: a model for the anisotropic MJ distribution [5], an astrophysical application of the MJ distribution to the energy distribution in radio jets [6], a new family of MJ distributions characterized by the parameter
[7] and an application to counter-streaming beams of charged particles [8]. 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
(1)
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
(2)
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
(3)
where the rest energy has been subtracted from the total energy, see formula (23.1) in [9]. A relativistic MB distribution can be obtained from Equation (1) replacing the classical kinetic energy
with the relativistic kinetic energy
(4)
where the relativistic temperature; T, is expressed in
units; up to now the treatment is the same of [10] 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
(5)
The differential of the velocity,
,
(6)
and therefore the relativistic MB distribution in the variable
is
(7)
where G is the Mejier G-function [11] [12] [13]; Figure 1 reports the above PDF for three different temperatures.
The average value or mean,
, is
(8)
the second moment about the origin is
(9)
the variance,
is
(10)
The mode is the real solution of the following cubic equation in
(11)
which has the real solution
Figure 1. The PDF of the relativistic MB as function of
for different values of T.
(12)
At the moment of writing a closed form for the distribution function (DF) which is
(13)
does not exists and we therefore present a numerical integration, see Figure 2.
The asymptotic behavior of the PDF,
, is
(14)
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 [14]. The above PDF has only one parameter which can be derived approximating the average value with a Pade approximant
(15)
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.
(16)
Here
is the sample mean defined as
(17)
formula which is useful to derive the variance of the sample
(18)
where
are the n-data, see [15]. 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
(19)
where v is expressed in
units. The mode is a solution of a sextic equation, see [16], in v
(20)
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.
(21)
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
(22)
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 [3] [4] [5] [6]. 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
(23)
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
.
(24)
The mode can be found by solving the following cubic equation
(25)
The real solution is
(26)
The asymptotic expansion of order 10 for the PDF is
(27)
The DF is evaluated with the following integral
(28)
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,
(29)
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
(30)
An analogous formula allows to derive
from the variance
of the sample
(31)
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).
Variable β
We now change the variable of integration
in
, the PDF of the MJ is
(32)
where
, see Figure 9. We have only one analytical result, the mode, which is found solving the following equation in
(33)
As an example when
the mode is at
and Figure 10 reports the mode as function of
.
The mean and the variance of the MJ distribution does not have an analytical expression and they are reported in a numerical way, see Figure 11 and Figure 12.
The DF of the MJ is given by the following integral
(34)
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 [17], when
(35)
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
(36)
where, according to eqn.(8.58) in [18], 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
(37)
where
is the modified Bessel function of second kind with order 5/3 [13] [19]. The modified Bessel function is also known as Basset function, modified Bessel function of the third kind or Macdonald function see pag. 527 in [20]. The above function has the following analytical expression
(38)
where
is a regularized hypergeometric function [19] [21] [22] [13]. 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
(39)
Figure 14. F(x) as function of x with logarithmic axes.
where E is the relativistic energy. The resulting PDF in relativistic energy is
(40)
A second change of variable is
(41)
produces
(42)
where
(43)
We know that
where B is the magnetic field expressed in gauss and therefore the above PDF in frequency becomes
(44)
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
(45)
The second PDF in
is
(46)
The astrophysical PDF in frequency for the Maxwell Jüttner distribution is
(47)
The mismatch between measured flux in Jy and theoretical flux,
, can be obtained introducing a multiplicative constant C
(48)
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 [23] 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 [24]. 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
.
5. Conclusions
The relativistic MB distribution has been derived in [10] 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.