and Maxwell Velocity Distribution

A widely applied concept in the physics of gases in general and plasmas in particular is that of a Local Thermodynamic Equilibrium (LTE) and the associated Maxwell velocity distribution. Yet in some branches of physics, like astrophysics for instance, the time scale of elastic collisions of particles, as given by the corresponding collision frequency, is often longer than the time scale of inelastic collisions and/or dynamic changes, and then these concepts can then strictly speaking not be applied anymore (i.e. collisional relaxation does not fully occur). In these non-LTE cases, the energy and velocity distribution has to be calculated self-consistently by means of a detailed equilibrium model. This is usually well recognized in practice with regard to the population of discrete atomic energy levels (which prove in this case not to follow the usual Boltzmann distribution), but the velocity distribution of atoms/ions and molecules is often even in these cases assumed to be a Maxwell distribution, (with all components of the gas having essentially the same temperature). Even where 'non-Maxwellian' distributions are considered, this is often only in the sense that a high energy tail is added to the Maxwell distribution, with the concept of a 'temperature' otherwise still being applied to the situation. For dilute gases, and especially plasmas (partially ionized gases), inelastic collisions are however on the contrary dominant for the whole distribution, as is evident from ionospheric electron spectra shown below

These distributions have been computed with a detailed equilibrium model, and they obviously show no sign of a Maxwell/Boltzmann energy distribution throughout the spectrum (their form being mainly determined by radiative recombination and inelastic collisions; see the page Electron Spectrum in the Ionosphere for more information). So in this case (and many others) collisional relaxation towards a Maxwellian velocity distribution simply doesn't take place for some constituents (i.e. the notion of a 'temperature' can strictly speaking not be applied) , and neither can an equipartition of energies between the constituents be assumed to hold here (the latter is especially true if the constituents have very different masses (e.g. electrons and ions), as then the energy exchange for an elastic collision is very small). One could be tempted to conclude from this that if inelastic collisions and dynamic processes are absent (i.e.if only elastic collisions affect the particle energies), the distribution function must necessarily be Maxwellian. This is indeed the conclusion at which Maxwell himself arrived in his 'Dynamical Theory of Gases'. However, in his derivation he makes an assumption which is equivalent to implying a normal distribution (i.e a Gaussian- (Maxwell-) distribution) in the first place (namely by assuming that the initial velocities of the two molecules are perpendicular to each other). Subsequent derivations of the Maxwell velocity distribution using Boltzmann's H-theorem are similarly incorrect in claiming that the Maxwell distribution is not only a sufficient but also a necessary condition for a collisional equilibrium. For instance, Saha in his book 'A Treatise on Heat' (pp.827-833) claims that dH/dt=0 would be synonymous with df/dt=0 as H is merely an integral of f over velocity. It should be clear however that generally ∫F(x)dx=0 does not imply F(x)=0. The conclusion is only correct if one restricts F(x) correspondingly. This is exactly what Saha did in his derivation as he implicitly assumed an isotropic distribution function by assuming the indices of the molecules as inter-changeable. This suggests already that the required condition for the equilibrium distribution being Maxwellian is that it is isotropic, and it is shown in the following that this is indeed the case: assuming a perfectly elastic head-on- (i.e. a 1D-) collision of two particles with masses m_{1} and m_{2}, the energy- and momentum conservation laws result in the well-known relationship for the velocity v_{1}' of particle 1 after the collision:
_{1} and K_{2} are the kinetic energies of the particles before the collision.
Since energy conservation demands that

Eq.(5) is the result for a 1D- (head-on) collision. This can easily be generalized to a 3D situation, by noting that the average energy transfer for a 3-D collision is a factor 8/5π/4 smaller than that for a head-on collision (as is evident from a comparison of Eqs.(3) and (11) on my page Coulomb Collision Cross Section and Coulomb Logarithm). Multiplying this value to the factor in front of the square bracket gives thus the corresponding 3D result_{1}=m_{2} or <v_{1}>^{.}<v_{2}>=0 (or both). In this case Eq.(6) simplifies to
_{1}>=0 and <v_{2}>=0), these will then approach a Maxwell distribution (with an identical temperature), i.e. collisional relaxation will take place. Even if the initial velocities are not isotropic, it is obvious that they become so after some time if one has a closed system. However, if the system is not closed, and one has significant gain- and loss processes determining the particle distribution function, the latter will not be Maxwellian even in equilibrium (this was effectively already addressed above by mentioning the effect of inelastic collisions, but it could also simply be gains and losses due to gas entering or leaving a given volume). This circumstance is, according to Eq.(7), particularly relevant if the masses m_{1} and m_{2} are very different, as then the energy change is very small and it takes a long time for the energy difference to become smaller. In this case, the equilibrium distribution will in general not be given by a Maxwellian. By noting that the condition for a general equilibrium is <K_{2}'>-<K_{1}'>= <K_{2}>-<K_{1}> , this requires according to Eq.(6) merely
_{1}=m_{2}, the two constituents will have different average kinetic energies in equilibrium if <v_{1}> and <v_{2}> are different from zero. And even if the latter are zero, the kinetic energies may stay different for constituents with much different masses, as then, according to Eq.(7), the energies will change very little in a collision, and collisional equilibrium (through elastic collisions) will not be established. The equilibrium will in this case be determined by other processes (e.g. inelastic collisions and recombination as for the electron spectra shown in the figure above).

This shows that the Maxwell distribution is a much stronger assumption than that of an equilibrium in general. The former would in terms of the distribution function f require that*Note: the consideration on this page has very much benefited from an extended email discussion with Dr. P. Radhakrishnamurty*

These distributions have been computed with a detailed equilibrium model, and they obviously show no sign of a Maxwell/Boltzmann energy distribution throughout the spectrum (their form being mainly determined by radiative recombination and inelastic collisions; see the page Electron Spectrum in the Ionosphere for more information). So in this case (and many others) collisional relaxation towards a Maxwellian velocity distribution simply doesn't take place for some constituents (i.e. the notion of a 'temperature' can strictly speaking not be applied) , and neither can an equipartition of energies between the constituents be assumed to hold here (the latter is especially true if the constituents have very different masses (e.g. electrons and ions), as then the energy exchange for an elastic collision is very small). One could be tempted to conclude from this that if inelastic collisions and dynamic processes are absent (i.e.if only elastic collisions affect the particle energies), the distribution function must necessarily be Maxwellian. This is indeed the conclusion at which Maxwell himself arrived in his 'Dynamical Theory of Gases'. However, in his derivation he makes an assumption which is equivalent to implying a normal distribution (i.e a Gaussian- (Maxwell-) distribution) in the first place (namely by assuming that the initial velocities of the two molecules are perpendicular to each other). Subsequent derivations of the Maxwell velocity distribution using Boltzmann's H-theorem are similarly incorrect in claiming that the Maxwell distribution is not only a sufficient but also a necessary condition for a collisional equilibrium. For instance, Saha in his book 'A Treatise on Heat' (pp.827-833) claims that dH/dt=0 would be synonymous with df/dt=0 as H is merely an integral of f over velocity. It should be clear however that generally ∫F(x)dx=0 does not imply F(x)=0. The conclusion is only correct if one restricts F(x) correspondingly. This is exactly what Saha did in his derivation as he implicitly assumed an isotropic distribution function by assuming the indices of the molecules as inter-changeable. This suggests already that the required condition for the equilibrium distribution being Maxwellian is that it is isotropic, and it is shown in the following that this is indeed the case: assuming a perfectly elastic head-on- (i.e. a 1D-) collision of two particles with masses m

(1) v_{1}' = [2m_{2}v_{2} + (m_{1}-m_{2})v_{1}]/ (m_{1}+m_{2}) ,

(2) K_{1}' = m_{1}/2^{.}v_{1}'^{2} = K_{1} + 4^{.}m_{1}m_{2}/(m_{1}+m_{2})^{2} ^{.}[K_{2}-K_{1} +v_{1}^{.}v_{2}^{.}(m_{1}-m_{2})/2] ,

(3) K_{1}'+K_{2}' = K_{1}+K_{2} = K = const. ,

(4) K_{2}'-K_{1}'= K_{2}-K_{1} - 8^{.}m_{1}m_{2}/(m_{1}+m_{2})^{2} ^{.}[K_{2}-K_{1} +v_{1}^{.}v_{2}^{.}(m_{1}-m_{2})/2] ,

(5) <K_{2}'>-<K_{1}'>= <K_{2}>-<K_{1}> -8^{.}m_{1}m_{2}/(m_{1}+m_{2})^{2} ^{.}[<K_{2}>-<K_{1}> +<v_{1}>^{.}<v_{2}>^{.}(m_{1}-m_{2})/2] ,

Eq.(5) is the result for a 1D- (head-on) collision. This can easily be generalized to a 3D situation, by noting that the average energy transfer for a 3-D collision is a factor 8/5π/4 smaller than that for a head-on collision (as is evident from a comparison of Eqs.(3) and (11) on my page Coulomb Collision Cross Section and Coulomb Logarithm). Multiplying this value to the factor in front of the square bracket gives thus the corresponding 3D result

(6) <K_{2}'>-<K_{1}'>= <K_{2}>-<K_{1}> -16/5π ^{.}m_{1}m_{2}/(m_{1}+m_{2})^{2} ^{.}[<K_{2}>-<K_{1}> +<v_{1}>^{.}<v_{2}>^{.}(m_{1}-m_{2})/2] .

(7) <K_{2}'>-<K_{1}'>= (<K_{2}>-<K_{1}>) ^{.}(1 -16/5π ^{.}m_{1}m_{2}/(m_{1}+m_{2})^{2}) ,

(8) <K_{2}>-<K_{1}> +<v_{1}>^{.}<v_{2}>^{.}(m_{1}-m_{2})/2 = 0 ,

This shows that the Maxwell distribution is a much stronger assumption than that of an equilibrium in general. The former would in terms of the distribution function f require that

(9) f(K_{1}')^{.}f(K_{2}')=f(K_{1})^{.}f(K_{2})

(10) e^{-k1'}^{.}e^{-k2'} = e^{-(k1'+k2')} = e^{-(k1+k2)} = e^{-k1}^{.}e^{-k2} .

Thomas Smid (M.Sc. Physics, Ph.D. Astronomy)

email: