Misinterpretation of the Maxwell velocity distribution?

In summary, the most probable speed for a particle in a Maxwell-Boltzmann distribution is not zero, as it is determined by the Boltzmann factor multiplied by the density of states. This density of states is proportional to the energy, resulting in a peak at a non-zero speed. However, the average velocity is zero, and the distribution does not predict that there are no stopped molecules in the population. The distribution for speed is derived from the phase-space distribution and takes into account the appropriate Jacobians. Overall, the distribution of speed in a monoatomic gas is expected to follow a half Gaussian rather than a bell curve.
  • #1
Kairos
177
14
I suddenly have a problem with the famous Maxwell velocity distribution. The maximum of this bell-shaped curve is commonly interpreted as the most probable velocity for a particle (see for instance https://en.wikipedia.org/wiki/Maxwell–Boltzmann_distribution: "The most probable speed, vp, is the speed most likely to be possessed by any molecule..")
But according to Boltzmann's distribution the most probable energy for a particle is zero (decreasing exponential: P(E=Ei) = Exp(-Ei/kT)/Σ), so we should expect that the most probable velocity for a randomly picked particle is zero as well? but the probability of velocity zero is zero in the Maxwell distribution!
thank you in advance if you can enlighten me
 
Physics news on Phys.org
  • #2
The Boltzmann factor is ## e^{-E/(k_B T)} ##, but to get the Maxwellian distribution, this gets multiplied by the density of states per unit energy interval ## n(E) ##. If you look over the derivation carefully, you will see how this factor ## n(E) ## arises in the computation of the number of states that is first counted in momentum space. The factor ## n(E) ## is proportional to ## E^{1/2} ##, making the distribution function equal to zero at ## E=0 ##.
 
  • Like
Likes Kairos
  • #3
Indeed and my question is about the physical meaning of this Maxwellian distribution.
Very concretely: the most probable speed of any particle is ## \sqrt{2k_BT/m} ## or zero ?
 
  • #4
The most probable speed is ## v_{mp}=\sqrt{\frac{2k_BT}{m}} ##. The average speed ## \bar{v}=\sqrt{\frac{8 k_B T}{\pi \, m}} ##, and the r.m.s. velocity is ## v_{rms}=\sqrt{\frac{3 k_B T}{m}} ##. It can be helpful to memorize each of these. ## \\ ## And no, the most probable speed is not zero, basically for the same reason that the most probable energy is not zero either. ## \\ ## Edit: It might interest you that the average velocity is zero. ## \\ ## It might also interest you, because the density of states is uniform w.r.t. volume on such a graph, if we make a 3-D graph with x-velocity, y-velocity, and z-velocity on the 3 axes, and then plotted (with density proportional to the shading) a probability density function on it=the Boltzmann factor, then you would indeed find the most probable velocity would, in fact, be at ## \vec{v}=(v_x,v_y, v_z)=(0,0,0) ##.
 
Last edited:
  • #5
Charles Link said:
The most probable speed is vmp=√2kBTm

this result is indeed the peak of the Maxwell-Boltzmann distribution but is it really the most probable speed of a molecule? I don't think the ordinate of this distribution is the probability of velocity. If that were the case, this distribution would also predict that there are no stopped molecules in the population, which I do not find obvious. Do you think that no molecule is motionless in a box of air in equilibrium?

Charles Link said:
the most probable energy is not zero
Boltzmann energy probability is ## P(E=Ej)= e^{-Ej/(k_B T)}/\sum_i e^{-Ei/(k_B T)} ##. It is maximal when ## Ej=0 ## isn't it?
 
  • #6
I recommend you study this some more. See also https://chem.libretexts.org/Textboo...hase_Kinetics/Maxwell-Boltzmann_Distributions ## \\ ## To answer your last question, the sum is over the states ## j ##. The result is the distribution function. In counting states, they are counted in momentum space. See also https://www.physicsforums.com/threads/boltzmann-vs-maxwell-distribution.918232/#post-5788233 In determining the density of states, see posts 2 and 4 of this "link". In counting the density of states in energy space, ## E=p^2/(2m) ##, the result is a continuous function rather than a discrete one, and this continuous function is proportional to ## E^{1/2} ##. The result is that the density of states that get counted in the sum over ## j ## has value zero at ## E=0 ##.
 
Last edited:
  • Like
Likes Kairos
  • #7
You must be careful to distinguish "speed" and "velocity". Anyway, velocity is not the best concept in kinetic theory, because as all of physics from a fundamental point of view everything happens in phase space, and there momentum rather than the velocity is the "natural" variable.

For an ideal gas the Boltzmann distribution in non-relativistic approximation indeed reads
$$f(\vec{p})=\exp[-(\vec{p}^2/(2 m)-\mu)/T],$$
where ##\mu## is the chemical potential and ##T## the temperature. This distribution gives the number of particles in a phase-space cell ##\mathrm{d}^3 r \mathrm{d}^3p## to be
$$\mathrm{d} N = \mathrm{d}^3 r \mathrm{d}^3p \frac{1}{(2 \pi \hbar)^3} f(\vec{p}).$$
Now if you want distributions for other quantities, which are functions of phase-space positions ##(\vec{r},\vec{p})## you have to consider the appropriate Jacobians for the phase-space volume.

E.g., if you want the distribution for speed you need
$$v^2=\vec{p}^2/m^2,$$
which leads to
$$\mathrm{d}^3 p = \mathrm{d} p \mathrm{d}^2 \Omega p^2.$$
Now
$$p^2=m^2 v^2 \; \Rightarrow \; \mathrm{d}p p=m^2 \mathrm{d} v v.$$
Then you get
$$\mathrm{d}^3 p =\mathrm{d} v \mathrm{d}^2 \Omega m^3 v^2,$$
which implies
$$\mathrm{d} N = \mathrm{d}^3 r \mathrm{d} v 4 \pi m^3 v^2 \frac{1}{(2 \pi \hbar)^3} \exp[-(m v^2/2-\mu)/T].$$
Note that the proper treatment of the Jacobian always yields the right dimension: ##\mathrm{d} N## must be dimensionless, and that's the case for both the phase-space distribution as well as the speed distribution derived from it (thanks to the Planck action quantum ##\hbar## in the formula!).
 
  • Like
Likes Demystifier and Charles Link
  • #8
Charles Link said:

I hadn't seen that thread thank you; my question is indeed close to Toby's:
The particle energy E is distributed over a Boltzmann distribution. In a purely translational system like a monoatomic gas, all the energy is kinetic (## E=mv^2/2 ##), so ## v ## is not expected to follow a bell curve but a half gaussian.

By (mis?)interpreting the Maxwell bell curve as the probabilities of values of ## v ##, the immobile molecules are absent from the system, whereas with the half Gaussian, the most frequent value of ## v ## for a molecule is zero: a radically different result, hence my question: do you think that motionless molecules are present or absent in a gas?
 
  • #9
Again you must be careful. The energy distribution is again different. Here you need to transform from ##\vec{p}## to ##E##. Since
$$E=p^2/(2m) \; \Rightarrow \; \mathrm{d} E=\mathrm{d} p p/m$$
you get
$$\mathrm{d}^3 p = \mathrm{d} p \mathrm{d}^2 \Omega p^2 =\mathrm{d} E \mathrm{d}^2 \Omega m \sqrt{2m E},$$
and integrating over the solid angle you get
$$\mathrm{d} N =\mathrm{d}^3 r \mathrm{d} E 4 \pi m \sqrt{2mE} \frac{1}{(2 \pi \hbar)^3} \exp[-(E-\mu)/T].$$
I don't understand what you mean by "half Gaussian". The Boltzmann distribution is a distribution for ##\vec{p}## and not ##p##. See my previous posting!

Distributions are named distributions not only because of the physical meaning of distributions as some quantity per unit space or unit phase space but also mathematically, because they transform under variable transformations with the Jacobian of the transformation. The Jacobian is thus very important. Otherwise you get the problems you are asking about!
 
  • #10
You can use the speed distribution result to compute the fraction using an integral, that have speeds that are very small, e.g. speed ##v<.01 \sqrt{\frac{k_B T}{m}} ##. You can compare that to the fraction that have ## \sqrt{\frac{2 k_BT}{m}}<v<\sqrt{\frac{2 k_B T}{m}}+.01 \sqrt{\frac{k_B T}{m}} ##. You will find the fraction with small speeds to be very small, perhaps surprisingly small. I think Heisenberg's uncertainty principle keeps us from saying that they are ever completely motionless. ## \\ ## And change the .01 to a .001 and you will see an even more dramatic effect. The fraction will be extremely small for the molecules that have that slow of a speed compared to the number that are at the same size interval, but at the most probable speed.
 
  • #11
vanhees71 said:
I don't understand what you mean by "half Gaussian"

##
P(v)= \sqrt{\frac{2m}{\pi k_{B}T}} \ \large{e}^{-\frac{mv^{2}}{2 k_{B}T}}
##

Charles Link said:
the fraction with small speeds to be very small, perhaps surprisingly small.

To be very clear: do you mean "the fraction of particles with small speeds is very small ? Thanks for that clear answer, actually in agreement with the widespread interpretation of Maxwell's distribution... but I have trouble convincing myself of it!
 
  • #12
Yes, the fraction of motionless particles is vanishingly small. If you have trouble visualising it, think of the following.

Start with a box (you can keep it 2D for simplicity) in which you some balls, in the absence of gravity and friction. The balls can move around, collide with each other and the walls. Now, start with all balls motionless, except for one to which you impart a great force. Can you see how, after a relatively short while, you basically will not have a single ball standing still?
 
  • Like
Likes Charles Link
  • #13
Kairos said:
##
P(v)= \sqrt{\frac{2m}{\pi k_{B}T}} \ \large{e}^{-\frac{mv^{2}}{2 k_{B}T}}
##
This is a Gaussian (normalized) distribution, i.e., it gives the probability to find a particle with one component of the velocity ##v##. Here ##v \in \mathbb{R}##. It's not the correct probability distribution for the speed, which is ##|\vec{v}| \geq 0##.

To be very clear: do you mean "the fraction of particles with small speeds is very small ? Thanks for that clear answer, actually in agreement with the widespread interpretation of Maxwell's distribution... but I have trouble convincing myself of it!
I mean what I wrote: You have to carefully take into account the Jacobian when changing the variable of which you want to describe the distribution (no matter whether it's a probability distribution or a particle-number distribution).
 
  • #14
DrClaude said:
Yes, the fraction of motionless particles is vanishingly small. If you have trouble visualising it, think of the following.

Start with a box (you can keep it 2D for simplicity) in which you some balls, in the absence of gravity and friction. The balls can move around, collide with each other and the walls. Now, start with all balls motionless, except for one to which you impart a great force. Can you see how, after a relatively short while, you basically will not have a single ball standing still?

I see but if all the balls move, none is energyless, which is not predicted by Boltzmann's energy distribution over the particles?
Consider another way to see the balls in the box: all are moving and colliding, but following a collision, a fast ball can be brutally slowered and even stop ?
 
  • #15
vanhees71 said:
it gives the probability to find a particle with one component of the velocity ##v##.
It is the velocity norm (or the component along the trajectory) such that ## Ec=mv^2/2 ##
 
  • #16
Kairos said:
I see but if all the balls move, none is energyless, which is not predicted by Boltzmann's energy distribution over the particles?
You get a similar distribution for energy, since there is a simple relation between speed and kinetic energy. I think you are forgetting to factor in the degeneracy factor (or density of states). You get a similar situation when consider the distribution of rotational states of molecules in a gas, see https://commons.wikimedia.org/wiki/...otational_states_of_nitrogen_gas_at_298_K.png
Kairos said:
Consider another way to see the balls in the box: all are moving and colliding, but following a collision, a fast ball can be brutally slowered and even stop ?
The probability of a collision where one of the balls would stop is extremely small.
 
  • Like
Likes Charles Link
  • #17
Kairos said:
I see but if all the balls move, none is energyless, which is not predicted by Boltzmann's energy distribution over the particles?
Consider another way to see the balls in the box: all are moving and colliding, but following a collision, a fast ball can be brutally slowered and even stop ?
You need to remember that there is motion in 3-dimensions for these particles. If the motion of the particles were all in one dimension, there would be a tremendously higher probability that the particle would come to a stop after a collision. If it was just one-dimensional motion, the speed distribution function would be the half-Gaussian that you previously mentioned. ## \\ ## You could observe the motion in just one dimension and have an artificial speed function for that one-dimensional motion, and I think that might be exactly what you are trying to do. ## \\ ## In computing the speed distribution function, all 3 dimensions are part of that distribution. In the process, it picks up a ## v^2 ## factor on the Boltzmann factor in the conversion of counting states: ## dp_x \, dp_y \, dp_z=4 \pi \, p^2 \, dp ## with a ## (4 \pi ) \, m^3 ## factored out to give ## v^2 \, dv ##.
 
Last edited:
  • #18
Kairos said:
To be very clear: do you mean "the fraction of particles with small speeds is very small ? Thanks for that clear answer, actually in agreement with the widespread interpretation of Maxwell's distribution... but I have trouble convincing myself of it!

In making the step from the Boltzmann “energy distribution” to the Maxwell “speed distribution”, the Boltzmann energy distribution function must be multiplied by an velocity dependent factor to account for the density of velocity states available to the particles. Have a look at http://hyperphysics.phy-astr.gsu.edu/hbase/Kinetic/maxspe.html#c4
 
  • Like
Likes Charles Link
  • #19
Lord Jestocost said:
In making the step from the Boltzmann “energy distribution” to the Maxwell “speed distribution”, the Boltzmann energy distribution function must be multiplied by an velocity dependent factor to account for the density of velocity states available to the particles. Have a look at http://hyperphysics.phy-astr.gsu.edu/hbase/Kinetic/maxspe.html#c4
Yes. And that ## v^2 ## factor results for the reason I mentioned in post 17.
 
  • Like
Likes Lord Jestocost
  • #20
thank you very much for all your kind answers; I will think about them deeply... but slowly
Regards
 
  • Like
Likes Charles Link
  • #21
Consider the statement in #4:
Charles Link said:
It might also interest you, because the density of states is uniform w.r.t. volume on such a graph, if we make a 3-D graph with x-velocity, y-velocity, and z-velocity on the 3 axes, and then plotted (with density proportional to the shading) a probability density function on it=the Boltzmann factor, then you would indeed find the most probable velocity would, in fact, be at ## \vec{v}=(v_x,v_y, v_z)=(0,0,0)## .
The most probable velocity - that is, the most probable individual combination of vx, vy and vz values - is (0,0,0). Any other particular combination is less likely, but for a given speed v > 0 there are many combinations that give the same v, e.g. for v = √3 you could have (1,1,1) or (1,-1,-1) or (√2,1,0) or (0,√3,0) or many more. The probability of v in the interval v → v + dv is proportional to 4πv2dv f(v). That is why v = 0 is so improbable.
Likewise in energy terms, the density of states increases with energy, so with appropriate coordinate transformation, the distribution function is proportional to E1/2e-E/kT, and it does not have a maximum at E = 0.
 
  • Like
Likes Kairos, vanhees71, Demystifier and 1 other person
  • #22
Lord Jestocost said:
In making the step from the Boltzmann “energy distribution” to the Maxwell “speed distribution”, the Boltzmann energy distribution function must be multiplied by an velocity dependent factor to account for the density of velocity states available to the particles.
I'm confused by what you are saying here. When I hear "distribution," I figure that this means a probability distribution function. In that case, even the energy distribution will contain a factor ##\sqrt{E}## coming from the density of states (as @mjc123 just mentioned while I was typing this). I think you are meaning the Boltzmann factor in stead of "Boltzmann energy distribution."
 
  • Like
Likes Lord Jestocost and Kairos
  • #23
mjc123 said:
The most probable velocity - that is, the most probable individual combination of vx, vy and vz values - is (0,0,0). Any other particular combination is less likely, but for a given speed v > 0 there are many combinations that give the same v, e.g. for v = √3 you could have (1,1,1) or (1,-1,-1) or (√2,1,0) or (0,√3,0) or many more. The probability of v in the interval v → v + dv is proportional to 4πv2dv f(v).

understood. Now imagine a graph whose ordinate would be "fractions of particles" (not density of velocity states) and the abscissa would be "velocity norms" (in which the different 3D combinations you enumerate collapse). According to you, this graph would be bell-shaped or monotonically decreasing?

mjc123 said:
That is why v = 0 is so improbable

v=0 is indeed reduced to a single combination of 3D components, but to me not improbable because the velocity component combinations are equiprobable only for given norms, owing to the principle of symmetry, while the norms are directly related to the energy distribution
 
  • #24
@Kairos For 3-D combinations, what gets compared is the volume of a sphere at the origin of radius ## dv ##, resulting in ## V_1=\frac{4}{3} \pi (dv)^3 ##,compared to the volume of a spherical shell of surface area ## 4 \pi \, v^2 ## and thickness ## dv ##, resulting in ## V_2=4 \pi \, v^2 \, dv ##. The spherical shell has tremendously more volume. The spherical shell also has a non-unity Boltzmann factor, but the overwhelming difference is the volumes of the two. ## \\ ## In the previous post, you ask, is the curve bell-shaped? The curve is ## f(v)=A v^2 e^{-m v^2/(2 k_B T)} ##, (## v>0 ##). It looks a little like a bell, but it is non-Gaussian.
 
Last edited:
  • #25
Hi, I rethought this question after your advices. If we momentarily drop the calculations with the shell volumes and return to physics: Boltzmann says that in an ideal monoatomic gas, there is a certain fraction of particles without energy, then after conversion into velocities, all the molecules start moving! (see #12) This calculation would therefore have the magical effect of giving energy to energyless molecules. For me, this paradox is unacceptable and is explained by a wrong interpretation of the Maxwell distribution ordinate: it does not represent the fractions of particles at the different velocities, but velocity densities (a more abstract notion!) where the runs at the same speed from right to left and from left to right are distinguished (as explained in #21). I can't give up the idea that the velocity occupied by the largest number of molecules could be 0 and not ## \sqrt{2k_BT/m} ##.
 
  • #26
Kairos said:
Boltzmann says that in an ideal monoatomic gas, there is a certain fraction of particles without energy, then after conversion into velocities, all the molecules start moving! (see #12)
That's what I was addressing in #22. Your statement is incorrect, as you are confusing the Boltzmann factor with the distribution of energy. Yes, the probability of being in a specific state of energy ##E_s## is ##\propto e^{-\beta E_s}##, but the probability of being in any state of energy ##E##, for an ideal monatomic gas, is ##\propto E^{1/2} e^{-\beta E}##, because you have to take into account how many states of that energy there are.

For an ideal gas particle, there is a single state of zero energy, so the probability of being in that sate is extremely small compared to other states of higher energy (up to a certain limit).
 
  • Like
Likes Kairos and Charles Link
  • #27
ok, I missed this important subtlety. many thanks!
 
  • #28
mjc123 said:
Consider the statement in #4:

The most probable velocity - that is, the most probable individual combination of vx, vy and vz values - is (0,0,0). Any other particular combination is less likely, but for a given speed v > 0 there are many combinations that give the same v, e.g. for v = √3 you could have (1,1,1) or (1,-1,-1) or (√2,1,0) or (0,√3,0) or many more. The probability of v in the interval v → v + dv is proportional to 4πv2dv f(v). That is why v = 0 is so improbable.
Likewise in energy terms, the density of states increases with energy, so with appropriate coordinate transformation, the distribution function is proportional to E1/2e-E/kT, and it does not have a maximum at E = 0.
That may be especially confusing to people whose native language is not English, because in some languages (such as my own Croatian language) the speed and velocity are represented by the same word, so such people may fail to distinguish velocity ##{\bf v}## from speed ##v=|{\bf v}|##. In short, the velocity ##{\bf v}=(0,0,0)## is the most probable velocity, but the speed ##v=0## is not the most probable speed.
 
  • #29
vanhees71 said:
You must be careful to distinguish "speed" and "velocity".
Does German language has separate words for speed and velocity? (The Croatian language doesn't.)
 
  • #30
Kairos said:
It is the velocity norm (or the component along the trajectory) such that ## Ec=mv^2/2 ##
It's better to write it in full vectorial form to make clear which quantity we talk about. It's also better to use momentum instead of velocity (although in Newtonian physics it's only a factor ##m## difference). The phase-space distribution is
$$\frac{\mathrm{d}^6 N}{\mathrm{d}^3 r \mathrm{d}^3 p}=\frac{1}{(2 \pi \hbar)^3} \exp \left [-\frac{1}{k_{\text{B}} T} \left (\frac{\vec{p}^2}{2m}-\mu \right ) \right].$$
Here I also introduced the Boltzmann constant. Before I used natural units, where temperatures are measured in energy units. Maybe it's better to use the SI to make everything clear concerning dimensions.

I've explained in detail, how to get the distributions wrt. the speed ##v=|\vec{v}|## or energy ##E=\vec{p}^2/(2m)## in my previous postings in this thread.
 
  • #31
Demystifier said:
Does German language has separate words for speed and velocity? (The Croatian language doesn't.)
No. In German unfortunately both are named with the same word "Geschwindigkeit". English has much more words than other European languages and thus is sometimes easier to express things more accurately with less words. Another example is that in English you can easily distinguish the sciences from the humanities, while in German both are called "Wissenschaft"...

Of course, the only adequate language to express physical laws after all is mathematics, and this is international.
 
  • Like
Likes Demystifier and Charles Link
  • #32
vanhees71 said:
Another example is that in English you can easily distinguish the sciences from the humanities, while in German both are called "Wissenschaft"...
The same with Croatian, both are "znanost".
vanhees71 said:
English has much more words than other European languages and thus is sometimes easier to express things more accurately with less words.
There is an interesting counterexample; English has the same word for singular and plural versions of "you". In Croatian those are different words. How about German?
 
  • #33
DrClaude said:
I think you are meaning the Boltzmann factor in stead of "Boltzmann energy distribution."

You are right!
 
  • #34
Demystifier said:
The same with Croatian, both are "znanost".

There is an interesting counterexample; English has the same word for singular and plural versions of "you". In Croatian those are different words. How about German?
In German we have "du" (singular) and "ihr" (plural).

I'm not sure, why in modern English the singular version "though" is not used anymore.
 
  • Like
Likes Demystifier
  • #35
DrClaude said:
distribution of energy. Yes, the probability of being in a specific state of energy EsEsE_s is ∝e−βEs∝e−βEs\propto e^{-\beta E_s}, but the probability of being in any state of energy EEE, for an ideal monatomic gas, is ∝E1/2e−βE∝E1/2e−βE\propto E^{1/2} e^{-\beta E}, because you have to take into account how many states of that energy there are.

I trouble you again for a detail: if I am not mistaken, the maximum of this energy distribution is ## k_B T/2 ##. Since the most probable velocity is ## \sqrt{2k_B T/m}##, ## E_{mp} \neq \frac{1}{2}mv_{mp}^{2} ## by a factor 2. Is that normal? (I guess a calculation with momenta would not give a different result)
 

Similar threads

Replies
3
Views
1K
  • Electromagnetism
Replies
4
Views
2K
Replies
4
Views
3K
Replies
8
Views
1K
  • Electromagnetism
Replies
2
Views
1K
  • Electromagnetism
Replies
11
Views
3K
Replies
1
Views
2K
  • Electromagnetism
Replies
6
Views
4K
Back
Top