The Maxwell-Boltzmann (MB) distribution, see  , 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 ; propagation of longitudinal electron waves in a collisionless, homogeneous, isotropic plasma, whose velocity distribution function is a truncated MB ; fast ion production in laser plasma ; the release of a dust particle from a plasma-facing wall ; an explanation of an anomaly in the Dark Matter (DAMA) experiment ; a distorted MB distribution of epithermal ions observed associated with the collapse of energetic ions ; and deviations to MB distribution that could have observable effects which can be measured trough the vapor spectroscopy at an interface . 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
where a is a parameter and v denotes the velocity, see  . Conversion to the physics is done by introducing the variable a, which is defined as
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
where the index p stands for physics. The distribution function (DF), , is
The average value or mean, , is
the variance, , is
The rth moment about the origin for the MB distribution is, , is
is the gamma function, see . The root-mean-square speed , , can be obtained from this formula by inserting
see Equations (7-10-16) in . This equation allows us to derive the temperature once the root-mean-square speed is measured
The coefficient of variation (CV) is
which is constant. The first three rth moments about the mean for the MB distribution, , are
The mode is at
An approximate expression for the median can be obtained by a Taylor series of the DF around the mode. The approximation formula is
which has a percent error, , of in respect to the numerical value. The entropy is
where is the Euler-Mascheroni constant, which is defined as
see  for more details. The coefficient of skewness is
and the coefficient of kurtosis is
According to , 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
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
and is the error function, which is defined as
see . 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
The average value , is
The rth moment about the origin for the double truncated MB distribution is, ,
where is the Whittaker M function, see . The root-mean-square speed, , can be obtained from this formula by inserting , and is
The variance is defined as
and has the following explicit form
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
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 , 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  , 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
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
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 .
The average value of escape is defined as
In this integral, the following changes are made to the variables
where is the mode as represented by Equation (21). The flux is now
For more details see    . On adopting the parameters of Table 1 the Jeans escape flux for hydrogen is
The Jeans escape flux for Earth at varies between ; see  or Figure 1 in . and , see . 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
This integral can be solved by introducing the change of variable as given by Equation (48)
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
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.
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.
 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.
 Kishimoto, Y., Mima, K., Watanabe, T. and Nishikawa, K. (1983) Analysis of Fast-Ion velocity Distributions in Laser Plasmas with a Truncated Maxwellian Velocity Distribution of Hot Electrons. The Physics of Fluids, 26, 2308.
 Tomita, Y., Smirnov, R., Nakamura, H., Zhu, S., Takizuka, T. and Tskhakaya, D. (2007) Effect of Truncation of Electron Velocity Distribution on Release of Dust Particle from Plasma-Facing Wall. Journal of Nuclear Materials, 363, 264-269.
 Ida, K., Kobayashi, T., Yoshinuma, M., Akiyama, T., Tokuzawa, T., Tsuchiya, H., Itoh, K. and LHD Experiment Group (2017) Observation of Distorted Maxwell-Boltzmann Distribution of Epithermal Ions in LHD. Physics of Plasmas, 24, Article ID: 122502.
 Todorov, P., De Aquino Carvalho, J.C., Maurin, I., Laliotis, A. and Bloch, D. (2019) Search for Deviations from the Ideal Maxwell-Boltzmann Distribution for a Gas at an Interface. Proceedings of the SPIE, 11047, Article ID: 110470.
 Vidal-Madjar, A., Blamont, J.E. and Phissamay, B. (1974) Evolution with Solar Activity of the Atomic Hydrogen Density at 100 Kilometers of Altitude. Journal of Geophysical Research, 79, 233-241.
 Bertaux, J. (1974) Lhydrogène atomique dans lexosphère terrestre: Mesures dintensité et de largeur de raie de lémission lyman alpha à bord du satellite ogo 5 et interprétation. These détat, Université Paris.