Why is the most probable energy different from the speed?

  • Thread starter Thread starter JohnnyGui
  • Start date Start date
  • Tags Tags
    Energy Speed
AI Thread Summary
The discussion centers on the distinction between the most probable energy and the most probable speed in the context of the Maxwell-Boltzmann distribution. It is noted that the most probable energy level does not equal the most probable speed squared multiplied by half the mass, leading to confusion about the physical interpretation of this difference. Participants explore the idea that energy levels may encompass a range of speeds, suggesting that an infinitesimal change in energy could cover a broader range of speed values than a similar change in speed. The mathematical relationships governing these distributions are emphasized, highlighting that averages and most probable values are not the same. Overall, the conversation seeks to clarify the physical implications of these statistical mechanics concepts.
JohnnyGui
Messages
796
Reaction score
51
It appears that the most probable energy level according to the Maxwell-Boltzmann distribution is not equal to the most probable speed squared multiplied by ##\frac{1}{2}m##. The most probable speed has a different value.
$$E_{max} = \frac{k_BT}{2}$$
$$v_{max} = \sqrt{\frac{2k_BT}{m}}$$

I am trying to understand this physically rather than mathematically, but I have a hard time comprehending this.

An energy level with the highest probability means it has the most particles out of all other energy levels, therefore that same largest number of particles has the speed corresponding to that energy level.

If I’m imagining energy levels as containers and I draw a specific number of particles in each container, there is no way I can distribute the particles in a way to have the most probable speed in an energy container different from the most probable energy. Energy is tied to speed.

The only way I can think of how this can be possible is if in the continuous approach a certain energy level covers a range of different speed values. If that’s the case then I’d deduce further that the reason is because an infinitesimally small ##dE## covers a larger range of speed values dan an infinitesimally small speed ##dv## does.

Is this reasoning correct?
 
Last edited:
Physics news on Phys.org
Where are you getting your formulas for ##E_{max}## and ##v_{max}## from?
 
You take averages, and thus [EDIT corrected in view of #5]
$$\langle |\vec{v}| \rangle=\left \langle \sqrt{\vec{v}^2} \right \rangle \neq \sqrt{\langle \vec{v}^2 \rangle}.$$
 
Last edited:
  • Like
Likes Tazerfish, etotheipi, Ibix and 1 other person
Note that vanhees' comment is not only applicable to this one problem. This happens all the time with distributions.
 
  • Like
Likes vanhees71
vanhees71 said:
You take averages, and thus
$$\langle |\vec{v}| \rangle=\left \langle \sqrt{\vec{v}}^2 \right \rangle \neq \sqrt{\langle \vec{v}^2 \rangle}.$$
Minor quibble - shouldn't the middle expression have the square under the square root?
 
  • Like
Likes vanhees71
of course.
 
  • Like
Likes Ibix
JohnnyGui said:
If I’m imagining energy levels as containers and I draw a specific number of particles in each container, there is no way I can distribute the particles in a way to have the most probable speed in an energy container different from the most probable energy. Energy is tied to speed.
Of course there is. Try it. Here's a simple Python program with a triangular distribution. <v> is different from sqrt(<v^2>). The output is:
Mean velocity = 6.3313
Mean energy = 45.1507
Velocity of mean energy = 6.7194

Python:
import numpy as np

left = 0.0
mode = 9.0
right = 10.0
nSamples = 1000000
velocities = []
energies = []
for i in range(nSamples):
    v =  np.random.triangular(left, mode, right)
    velocities.append(v)
    energies.append(v * v)
   
vbar = np.mean(velocities)
ebar = np.mean(energies)

print("Mean velocity = %.4f"%vbar)
print("Mean energy = %.4f"%ebar)
print("Velocity of mean energy = %.4f"%(np.sqrt(ebar)))
 
  • Like
Likes vanhees71
Major quibble "most probable" is not the same as "average".
 
  • Like
Likes Tazerfish and JohnnyGui
PeterDonis said:
Where are you getting your formulas for EmaxEmax and vmaxvmax from?

They can be found here and http://people.physics.tamu.edu/syeager/teaching/222/hw4solution.pdf. For the 2nd link, please search for the terms "most probable".

@vanhees71 @Ibix @Vanadium 50 @phyzguy I'm aware that the average speed is not the same as the root-mean-square speed. I don't see how this has a link with the most probable speed. Please see @PeroK 's post #8 above.

The formula for those averages also differ from the most probable speed as shown here.

PeroK said:
Major quibble "most probable" is not the same as "average".

This is exactly what I thought.
 
Last edited:
  • #10
The "most probable speed", i.e., the ##v##, where the probability distribution ##f(v)## has a maximum is in general again different from the other two considered "average speeds".
 
  • Like
Likes Tazerfish
  • #11
JohnnyGui said:
They can be found here and http://people.physics.tamu.edu/syeager/teaching/222/hw4solution.pdf. For the 2nd link, please search for the terms "most probable".

@vanhees71 @Ibix @Vanadium 50 I'm aware that the average speed is not the same as the root-mean-square speed. I don't see how this has a link with the most probable speed.

The formula for those averages also differ from the most probable speed as shown here.

Please see @PeroK 's post #8 above.
This is exactly what I thought.

Do all particles have the same mass?
 
  • Like
Likes JohnnyGui
  • #12
vanhees71 said:
The "most probable speed", i.e., the ##v##, where the probability distribution ##f(v)## has a maximum is in general again different from the other two considered "average speeds".

That's what I'm saying. My question in that regard is how it can be physically explained that the most probable energy does not correspond to the most probable speed that belongs to that energy.

@PeroK Yes
 
  • #13
A distribution is a distribution, i.e., it transforms under changes of the independent variable. Let ##f(v)## be the velocity distribution. For the Maxwell distribution, it's ##f(v) = N v^2 \exp[-m v^2/(2 kT)]##.

Then you have ##E=m v^2/2m##. From this you get
$$\tilde{f}(E) \mathrm{d} E =f(v) \mathrm{d} v.$$
Now ##v=\sqrt{2mE}## and thus ##\mathrm{d} v=\mathrm{d} E \sqrt{m/(2E)}##, leading to
$$\tilde{f}(E) \mathrm{d} E = f(v) \sqrt{m/(2 E)} \mathrm{d} E \; \Rightarrow \; \tilde{f}(E)=\sqrt{\frac{m}{2 E}} f(\sqrt{2 m E}).$$
 
  • #14
vanhees71 said:
A distribution is a distribution, i.e., it transforms under changes of the independent variable. Let ##f(v)## be the velocity distribution. For the Maxwell distribution, it's ##f(v) = N v^2 \exp[-m v^2/(2 kT)]##.

Then you have ##E=m v^2/2m##. From this you get
$$\tilde{f}(E) \mathrm{d} E =f(v) \mathrm{d} v.$$
Now ##v=\sqrt{2mE}## and thus ##\mathrm{d} v=\mathrm{d} E \sqrt{m/(2E)}##, leading to
$$\tilde{f}(E) \mathrm{d} E = f(v) \sqrt{m/(2 E)} \mathrm{d} E \; \Rightarrow \; \tilde{f}(E)=\sqrt{\frac{m}{2 E}} f(\sqrt{2 m E}).$$

I understand that the formulas are different mathematically, but my question is how this can be imagined physically, for example using my scenario to distribute the particles in such a way to fulfill that criteria.

If that's not possible then how about the following reasoning?
JohnnyGui said:
The only way I can think of how this can be possible is if in the continuous approach a certain energy level covers a range of different speed values. If that’s the case then I’d deduce further that the reason is because an infinitesimally small dEdEdE covers a larger range of speed values dan an infinitesimally small speed dvdvdv does.
 
  • #15
JohnnyGui said:
The only way I can think of how this can be possible is if in the continuous approach a certain energy level covers a range of different speed values. If that’s the case then I’d deduce further that the reason is because an infinitesimally small ##dE## covers a larger range of speed values dan an infinitesimally small speed ##dv## does.

Is this reasoning correct?

Here's a slightly artificial example:

Consider a large finite set of particles, clustered round ##v = 0## and ##v = v_0##, with more around ##v_0##. If you square these numbers, the cluster round ##0## stays close to zero, but the cluster around ##v_0## spreads out more. And the density around ##v^2 = 0## may be greater than the density around ##v^2 = v_0^2##.

To use a bit of simple maths:

##v = \Delta v \approx 0 \ \Rightarrow \ v^2 = (\Delta v)^2 \approx 0##

But:

##v = v_0 + \Delta v \ \Rightarrow \ v^2 = v_0^2 + 2v_0 \Delta v##

Which is spread much further away from ##v_0^2##.
 
  • Like
Likes vanhees71
  • #16
PeroK said:
Here's a slightly artificial example:

Consider a large finite set of particles, clustered round ##v = 0## and ##v = v_0##, with more around ##v_0##. If you square these numbers, the cluster round ##0## stays close to zero, but the cluster around ##v_0## spreads out more. And the density around ##v^2 = 0## may be greater than the density around ##v^2 = v_0^2##.

To use a bit of simple maths:

##v = \Delta v \approx 0 \ \Rightarrow \ v^2 = (\Delta v)^2 \approx 0##

But:

##v = v_0 + \Delta v \ \Rightarrow \ v^2 = v_0^2 + 2v_0 \Delta v##

Which is spread much further away from ##v_0^2##.

Ah, so in the context of my reasoning. A ##dE## starting from low energy values (thus low speeds) would cover more particles of the same speed and a ##dE## starting from higher energy (higher speeds) covers only a portion of the particles with the same speed because they are more spread out when squared?
 
  • Like
Likes vanhees71 and PeroK
  • #17
@PeroK I noticed an inner clash afterwards, because I also know that ##dE = mv \cdot dv##. Doesn't this mean that ##dE## get wider at higher speeds and thus covers the increasing spread of the speed anyway?
 
  • #18
JohnnyGui said:
@PeroK I noticed an inner clash afterwards, because I also know that ##dE = mv \cdot dv##. Doesn't this mean that ##dE## get wider at higher speeds and thus covers the increasing spread of the speed anyway?
Isn't that just another way of saying what's already been said?
 
  • #19
PeroK said:
Isn't that just another way of saying what's already been said?

I don't think so because I said in post #16 that ##dE## covers only a portion of higher speeds because of the spread, but if ##dE## gets wider at higher speeds, it would compensate for that spread and thus cover all of them anyway. I.e. the "density" within ##dE## stays the same.
 
  • #20
JohnnyGui said:
I said in post #16 that ##dE## covers only a portion of higher speeds because of the spread, but if ##dE## gets wider at higher speeds, it would compensate for that spread and thus cover all of them anyway. I.e. the "density" within ##dE## stays the same.

The fact that these look different to you in vague ordinary language does not mean you have actually proposed two different possible answers to your question. The mathematical fact that ##dE = m v dv## is the answer to your question. If you think that means "gets wider at higher speeds", and you think that's something different from "covers only a portion of higher speeds", then what do you think "covers only a portion of higher speeds" means mathematically that makes it different?
 
  • #21
PeterDonis said:
The fact that these look different to you in vague ordinary language does not mean you have actually proposed two different possible answers to your question.

How an answer "looks" is also influenced by how one understands it. I'm not solely basing the difference on the "looks" of my answers, but on my own understandings when writing them.
As I see them, there is a difference since one answer says "covers a portion" and the other says "covers all"

PeterDonis said:
The mathematical fact that dE=mvdvdE=mvdvdE = m v dv is the answer to your question. If you think that means "gets wider at higher speeds", and you think that's something different from "covers only a portion of higher speeds", then what do you think "covers only a portion of higher speeds" means mathematically that makes it different?

"covers only a portion" in the sense that the value ##dE## stays as a constant and does not get wider with increasing speed.
 
  • #22
JohnnyGui said:
one answer says "covers a portion" and the other says "covers all"

The only one who used the words "covers a portion" was you. What mathematical fact are you basing that on?

JohnnyGui said:
"covers only a portion" in the sense that the value ##dE## stays as a constant

So what mathematical fact justifies this description? Note that nothing @PeroK said in any of his posts implies this.
 
  • #23
PeterDonis said:
The only one who used the words "covers a portion" was you. What mathematical fact are you basing that on?

The fact that he, as well as another user, liked my post #16 containing that very description of ##dE## as a statement I asked about would make me assume that my understanding about ##dE## was correct. Hence me asking about it after.

PeterDonis said:
So what mathematical fact justifies this description? Note that nothing @PeroK said in any of his posts implies this.

Perhaps I based my description on his mathematical fact because I understood it in another way than what he meant?

When reading him talking about the density of speed squared values changes around different ##v##'s, I concluded that the reason the probability gets influenced by that is because a ##dE## would eventually not cover all the speed squared values when they're spread more out, thus ignoring a number of particles that had speeds very close to a certain (high) ##v##.

I discarded this thought when I noticed that ##dE## gets wider with increasing speeds, thus thinking not ignoring a number of particles that are close to a certain (high) ##v## and and taking them all into account for the probability. Which brought me back to the initial OP question.
 
  • #24
JohnnyGui said:
When reading him talking about the density of speed squared values changes around different ##v##'s, I concluded

...the wrong thing. What he meant was simply that ##dE = m v dv##.

This is a good example of why it's better to use math, not ordinary language, when talking about physics.
 
  • Like
Likes vanhees71
  • #25
JohnnyGui said:
he, as well as another user, liked my post #16

Because they thought you were just trying to say ##dE = m v dv## in different words. If they had known you were describing something different and inconsistent with ##dE = m v dv##, they would not have liked your post.

Again, an example of why it's better to just use math instead of vague ordinary language.
 
  • #26
PeterDonis said:
...the wrong thing. What he meant was simply that ##dE = m v dv##

You asked about what I based it on and I took the time to explain to you how I understood it so you could perhaps pinpoint how I should understand it with ##dE = mvdv##.
Merely substituting my whole explanation with "...the wrong thing" and that he meant something that I already mentioned noticing it in the meantime won't help me understand this any better.

PeterDonis said:
Because they thought you were just trying to say ##dE = m v dv## in different words. If they had known you were describing something different and inconsistent with ##dE = m v dv##, they would not have liked your post.

I would not have known that. I would've expected mentioning "##dE## covering a portion" would be inconsistent with that formula.
 
Last edited:
  • #27
JohnnyGui said:
how I should understand it with ##dE = mvdv##.

Do you see how ##dE = m v dv## explains why the most probable energy is not the same as the most probable velocity?
 
  • #28
PeterDonis said:
Do you see how ##dE = m v dv## explains why the most probable energy is not the same as the most probable velocity?
It should be even clearer if we replace the phrase "most probable" with "greatest probability density".
 
  • Like
Likes JohnnyGui and Ibix
  • #29
PeroK said:
It should be even clearer if we replace the phrase "most probable" with "greatest probability density".
And if that doesn't help you, @JohnnyGui, draw the distribution twice as a function of ##v##. In one case, slice it up into bands of width ##\Delta v##. In the other, slice it up into bands of width ##\Delta E## - note that you'll have to determine the velocities corresponding to ##n\Delta E##. Then use the trapezium rule to approximate the area under the curve in each bin in both cases. Is the highest value in the same place?

Then imagine repeating the exercise with smaller and smaller ##\Delta v## and ##\Delta E## - eventually in the limit as ##\Delta\rightarrow d##. Do you expect the answers to ever be the same?
 
  • Like
Likes JohnnyGui
  • #30
@Ibix & @PeroK I think using your info's helped me a bit by explaining it as follows:

For a certain ##E_0## and a corresponding ##v_0##, the following must be valid:$$f(E_0)\cdot dE = f(v_0) \cdot dv$$
For ##f(E)##, the maximum number of particles lies within the range ##\frac{k_BT}{2} \geq \frac{k_BT}{2} + dE##.
For ##f(v)##, the maximum number of particles lies within the range ##\sqrt{\frac{2k_BT}{m}} \geq \sqrt{\frac{2k_BT}{m}}+dv## which, in energy terms, is equal to a range of ##k_BT \geq k_BT + \frac{dE}{\sqrt{2mE}}##

The maximum number of particles differ because for ##f(E)## a range width of ##dE## has been used while for ##f(v)## a range width of ##\frac{dE}{\sqrt{2mE}}## has been used.

Does this make sense?
 
  • Like
Likes Philip Koeck
  • #31
The complete solution is in #13. I still don't know, what else you need.
 
  • #32
vanhees71 said:
The complete solution is in #13. I still don't know, what else you need.

Never mind, I figured it out, thanks. My calcuations in post #30 are wrong, I didn't pay attention when doing them. I was overlooking the fact that the relation between ##dE## and ##dv## is a function of ##v## while the probability itself is also a function of ##v##. If the relation between ##dE## and ##dv## were a constant, then then ##E_{max}## would correspond to ##v_{max}## .

Is there actually a mathematical way to explain why the maximum probability density in terms of energy ##E_{max}## is exactly half the energy that corresponds to ##v_{max}##?
 
  • Like
Likes Philip Koeck
  • #33
JohnnyGui said:
I figured it out, thanks.

It doesn't look like it, since you come right back with the same question:

JohnnyGui said:
Is there a mathematical way to explain why the maximum probability density in terms of energy ##E_{max}## is exactly half the energy that corresponds to ##v_{max}##?

The math has already been discussed quite sufficiently in this thread. I think the discussion has run its course.

Thread closed.
 
  • Like
Likes vanhees71
Back
Top