We calculate the average momentum distribution and the average deviation from this distribution for a one-dimensional ideal gas of N atoms. The calculation is performed using analytical and numerical methods, and the results obtained using the two methods are compared. We use these results to show that in the limit of large N, almost all the possible microstates of an ideal gas have momentum distributions that are very close to the Maxwell distribution. We discuss the significance of this fact for understanding why an ideal gas approaches thermal equilibrium.

1.
The program used to perform the numerical experiment will be provided upon request.
2.
Computer simulations of the evolution of an ideal gas toward thermal equilibrium are discussed in
J.
Novak
and
A. B.
Bortz
, “
The evolution of the two-dimensional Maxwell–Boltzmann distribution
,”
Am. J. Phys.
38
(
12
),
1402
1406
(
1970
) and
A. D.
Boozer
, “
Time asymmetry in a dynamical model of the one-dimensional ideal gas
,”
Am. J. Phys.
76
(
11
),
1026
1030
(
2008
).
3.
The first property is established in
K.
Huang
,
Statistical Mechanics
, 2nd ed. (
Wiley
,
New York
,
1987
), Sec. 4.3.
4.
A close analog to the second property is sometimes established by employing the canonical rather than microcanonical ensemble. See, for example,
R. C.
Tolman
,
The Principles of Statistical Mechanics
(
Oxford U. P.
,
Oxford
,
1938
), p.
506
.
5.
An overview of various proofs of this result is given in
P. A.
Mello
and
T. A.
Brody
, “
A different proof of the Maxwell–Boltzmann distribution
,”
Am. J. Phys.
40
(
9
),
1239
1245
(
1972
).
6.
By an “N-sphere” we mean a sphere in N-dimensional Euclidean space. According to this convention a circle is a 2-sphere.
7.
Here θ(x) is the step function, defined such that θ(x)=1 for x>1, θ(x)=1/2 for x=0, and θ(x)=0 for x<0.
8.
Most treatments prove the equivalent of Eq. (5) only in the limit N. As shown in Sec. IV, it holds for arbitrary N.
9.
One can show that the expected number of trials is n(N)=2NN/ΩN, where ΩN, the solid angle for the unit N-sphere, is given by Eq. (11). The function n(N) grows very rapidly with increasing N, so this algorithm is practical only for small values of N.
10.
For simplicity we compute P¯(p) and D¯(p) only for values of p that are integer multiples of δp. This procedure is equivalent to making a histogram of the atom momenta in which the size of the momentum bins is δp.
11.
Equation (11) is derived in
M. E.
Peskin
and
D. V.
Schroeder
,
An Introduction to Quantum Field Theory
(
Addison-Wesley
,
Reading, MA
,
1995
), p.
249
.
12.
Equation (15) is derived in
R.
López-Ruiz
and
X.
Calbert
, “
Derivation of the Maxwell distribution from the microcanonical ensemble
,”
Am. J. Phys.
75
(
8
),
752
753
(
2007
), and in
J. R.
Ray
and
H. W.
Graben
, “
Small systems have non-Maxwellian momentum distributions in the microcanonical ensemble
,”
Phys. Rev. A
44
(
10
),
6905
6908
(
1991
).
[PubMed]
13.
The analog to Eq. (15) for a two-dimensional ideal gas is derived in
S.
Velasco
,
J. A.
White
, and
J.
Guemez
, “
Single particle energy and velocity distributions for finite simple systems in the microcanonical ensemble
,”
Eur. J. Phys.
14
(
4
),
166
170
(
1993
).
14.
One can show that Eq. (5) is also valid for N=2. This result must be established as a special case because Eq. (17) is not valid when N=2.
AAPT members receive access to the American Journal of Physics and The Physics Teacher as a member benefit. To learn more about this member benefit and becoming an AAPT member, visit the Joining AAPT page.