New Probability Distributions in Astrophysics: III. The Truncated Maxwell-Boltzmann Distribution ()
1. Introduction
The Maxwell-Boltzmann (MB) distribution, see [1] [2], is a powerful tool to explain the kinetic theory of gases. The range in velocity of this distribution spans the interval
, which produces several problems:
1) The maximum velocity of a gas cannot be greater than the velocity of light, c.
2) The kinetic theory is developed in a classical environment, which means that the involved velocities should be smaller than ≈1/10c.
These items point toward the hypothesis of an upper bound in velocity for the MB. We will now report some approaches, including an upper bound in velocity: the ion velocities parallel to the magnetic field in a low density surface of a ionized plasma [3]; propagation of longitudinal electron waves in a collisionless, homogeneous, isotropic plasma, whose velocity distribution function is a truncated MB [4]; fast ion production in laser plasma [5]; the release of a dust particle from a plasma-facing wall [6]; an explanation of an anomaly in the Dark Matter (DAMA) experiment [7]; a distorted MB distribution of epithermal ions observed associated with the collapse of energetic ions [8]; and deviations to MB distribution that could have observable effects which can be measured trough the vapor spectroscopy at an interface [9]. However, these approaches do not clearly cover the effect of introducing a lower and an upper boundary in the MB distribution, which is the subject that will be analyzed in this paper.
This paper is structured as follows. Section 2 reviews the basic statistics of the MB distribution and it derives a new approximate expression for the median. Section 3 introduces the double truncated MB and it derives the connected statistics. Section 4 derives the relationship for root-mean-square speed versus temperature in the double truncated MB. Finally, Section 5.2 derives a new formula for Jeans flux in the exosphere.
2. The Maxwell-Boltzmann Distribution
Let V be a random variable defined in
; the MB probability density function (PDF),
, is
(1)
where a is a parameter and v denotes the velocity, see [1] [2]. Conversion to the physics is done by introducing the variable a, which is defined as
(2)
where m is the mass of the gas molecules, k is the Boltzmann constant and T is the thermodynamic temperature. With this change of variable, the MB PDF is
(3)
where the index p stands for physics. The distribution function (DF),
, is
(4)
(5)
The average value or mean,
, is
(6)
(7)
the variance,
, is
(8)
(9)
The rth moment about the origin for the MB distribution is,
, is
(10)
(11)
where
(12)
is the gamma function, see [10]. The root-mean-square speed ,
, can be obtained from this formula by inserting
(13)
(14)
see Equations (7-10-16) in [11]. This equation allows us to derive the temperature once the root-mean-square speed is measured
(15)
The coefficient of variation (CV) is
(16)
which is constant. The first three rth moments about the mean for the MB distribution,
, are
(17)
(18)
(19)
The mode is at
(20)
(21)
An approximate expression for the median can be obtained by a Taylor series of the DF around the mode. The approximation formula is
(22)
(23)
which has a percent error,
, of
in respect to the numerical value. The entropy is
(24)
(25)
where
is the Euler-Mascheroni constant, which is defined as
(26)
see [10] for more details. The coefficient of skewness is
(27)
and the coefficient of kurtosis is
(28)
According to [12], a random number generation can be obtained via inverse transform sampling when the distribution function or cumulative distribution function,
, is known: 1) a pseudo number generator gives a random number R between zero and one; 2) the inverse function
is evaluated; and 3) the procedure is repeated for different values of R. In our case, the inverse function should be evaluated in a numerical way by solving for v the following nonlinear equation
(29)
(30)
where
and
are the two DF represented by Equations (4) and (5). As a practical example, by inserting in Equation (29)
and
, we obtain in a numerical way
.
3. The Double Truncated Maxwell-Boltzmann Distribution
Let V be a random variable that is defined in
; the double truncated version of the Maxwell-Boltzmann PDF,
, is
(31)
where
(32)
where
(33)
and
is the error function, which is defined as
(34)
see [10]. The physical meaning of a is still represented by Equation (2); however, due to the tendency to obtain complicated expressions, we will omit the double notation. The DF,
, is
(35)
The average value
, is
(36)
The rth moment about the origin for the double truncated MB distribution is,
,
(37)
where
(38)
where
is the Whittaker M function, see [10]. The root-mean-square speed,
, can be obtained from this formula by inserting
, and is
(39)
where
(40)
The variance
is defined as
(41)
and has the following explicit form
(42)
Although the coefficients of skewness and kurtosis for the truncated MB exist, they have a complicated expression.
4. A Laboratory Application
The temperature as a function of root-mean-square speed for the MB is given by Equation (15). In the truncated MB distribution, the temperature can be found by solving the following nonlinear equation
(43)
where
is not a theoretical variable but is the root-mean-square speed measured in the laboratory and
is given by Equation (39). The laboratory measures of
started with [13], where a
at 400˚C was found for a metallic vapor. In the truncated MB distribution, there are three parameters that can be measured in the laboratory from a kinematical point of view, as follows: the lowest velocity,
; the highest velocity,
; and the root-mean-square speed,
. Setting for simplicity
, we will now explore the effect of the variation of
on the root-mean-square speed; see Figure 1. The first example of the influence of the upper limit in velocity on the temperature is given by potassium gas [14] [15], in which molecular mass is 6.492429890 × 10−26 kg. In Figure 2, we evaluate in a numerical way the temperature when
and
is variable in the case of a measured value of
.
The second example is given by diatomic nitrogen, N2, in which molecular mass is 4.651737684 × 10−26 kg. In Figure 3, we evaluate the temperature when
and
is a variable in the case of a measured value of
.
5. The Jeans Escape
The standard formula for the escape of molecules from the exosphere is reviewed in the framework of the MB distribution. A new formula for the Jeans escape is derived in the framework of the truncated MB.
5.1. The Standard Case
In the exosphere, a molecule of mass m and velocity
is free to escape when
(44)
where G is the Newtonian gravitational constant, M is the mass of the Earth,
is the radius of the exosphere, R is the radius of the Earth and H is the altitude of the exosphere. The flux of the molecules that are living in the exosphere
is
(45)
where
is the number of molecules per unit volume and
is the average velocity of escape. In the presence of a given number of molecules per unit volume, the standard MB distribution in velocities in a unit volume,
, is
Figure 1. The theoretical root-mean-square speed as a function of the upper limit in velocity (continuous line) and standard value of the temperature (dotted line) when
and
.
Figure 2. Temperature as a function of the upper limit in velocity for Potassium (continuous line) and standard value of the temperature (dotted line) when
and
.
Figure 3. Temperature as a function of the upper limit in velocity for diatomic nitrogen, N2, (continuous line) and standard value of the temperature (dotted line) when
and
.
(46)
The average value of escape is defined as
(47)
In this integral, the following changes are made to the variables
(48)
Therefore,
(49)
with
(50)
where
is the mode as represented by Equation (21). The flux is now
(51)
For more details see [16] [17] [18] [19]. On adopting the parameters of Table 1 the Jeans escape flux for hydrogen is
(52)
and
(53)
The Jeans escape flux for Earth at
varies between
; see [20] or Figure 1 in [21]. and
, see [22]. Therefore, our choice of parameters is compatible with the suggested interval in flux.
5.2. The Truncated Case
The average value of escape for a truncated MB distribution,
, is
(54)
This integral can be solved by introducing the change of variable as given by Equation (48)
(55)
where
is the lower value of
and
is the upper value of
. The flux of the molecules that are living the exosphere in the truncated MB,
, is
(56)
The increasing flux of molecules is outlined when one parameter,
, is variable; see Figure 4. In other words, an increase in
produces an increase in the flux of the molecules. The dependence of the flux when two parameters are variable,
and
, is reported in Figure 5.
Table 1. Adopted physical parameters for the exosphere.
Figure 4. The flux of molecules as a function of
with parameters as in Table 1,
and
.
Figure 5. The flux of molecules as a function of
and
with parameters as in Table 1.
These Jeans escape fluxes for Earth are compatible with the observed values that were reported in Section 5.1.
6. Conclusion
This paper derived analytical formulae for the following quantities for a double truncated MB distribution: the PDF, the DF, the average value, the rth moment about the origin, the root-mean-square speed and the variance. The traditional correspondence between root-mean-square speed and temperature is replaced by the nonlinear Equation (43). The new formula (56) for the Jeans escape flux of molecules from an atmosphere is now a function of the lower and upper boundary in velocity.