- #36
PeterDonis
Mentor
- 47,580
- 23,850
MikeW said:Not a group of molecules, a single SF6 molecule.
A single molecule doesn't have a velocity distribution function, which is what the M-B distribution is supposed to model in this connection.
MikeW said:Initial M-B Distribution assumed at bottom for upward vector
Distribution of what? There's only one molecule. It doesn't have a velocity distribution. See above.
If what you are trying to model is a probability distribution for a single molecule having a given velocity, why do you think a M-B distribution is the right one?
MikeW said:looses velocity due to gravity
As others have pointed out, you are assuming that the molecule has no collisions as it rises. But your claims about temperature only have meaning to begin with for a large number of molecules, not a single molecule, in which case the single molecule you are focusing on will have many collisions as it rises. Your analysis does not take that into account.
MikeW said:different distribution at top BEFORE touching top surface.
"Different distribution", yes. But that does not just mean "a M-B distribution with a different peak" (even if we assume your initial assumption of a M-B distribution at the bottom is correct, which it might not be, and that the molecule has no collisions as it rises, which it will--see above).
Write the M-B distribution for a system with a single degree of freedom as follows (where I have assumed the temperature ##T## is in energy units, i.e., Boltzmann's constant ##k## is 1):
$$
f(E) = C \frac{E}{T^{2/3}} \exp \left( - \frac{E}{T} \right)
$$
where ##C## is a combination of constant factors that I won't bother to write down explicitly. Your assumption is that at the bottom, there is some ##T = T_b## such that the distribution satisfies this equation, i.e., we have
$$
f_b(E) = C \frac{E}{T_b^{2/3}} \exp \left( - \frac{E}{T_b} \right)
$$
Now we shift ##E## by a constant amount, the potential energy difference ##\Delta U## from the bottom to the top. Your claim is that the resulting distribution will also satisfy the above equation at some lower temperature ##T_t##. But that can't possibly be the case, because of the factor ##E / T^{2/3}## in front of the exponential. That is, in order for your claim to be true, ##T_t## would have to satisfy two mutually inconsistent conditions: it would have to be equal to ##T_b - \Delta U## (to satisfy the argument of the exponential) but also equal to ## \left( T_b - \Delta U \right)^{2/3}## (to satisfy the factor in front of the exponential). That is impossible, so your claim must be false. A similar argument applies when going from top to bottom.