An ideal gas of non-interacting helium atoms is in a thermal equilibrium at the temperature T = 3 0 ∘ C . What's the probability that an atom in this gas has a speed near the speed of sound, i.e. has a speed in the range 3 4 0 ± 5 m/s ?
Details and assumptions
This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try
refreshing the page, (b) enabling javascript if it is disabled on your browser and,
finally, (c)
loading the
non-javascript version of this page
. We're sorry about the hassle.
Energies in J corresponding to 335 m/s and 345 m/s are:
E 1 = 2 1 m ( 3 3 5 2 ) = 3 . 7 2 9 2 4 × 1 0 − 2 2 E 2 = 2 1 m ( 3 4 5 2 ) = 3 . 9 5 5 2 0 0 7 5 × 1 0 − 2 2
Using k = 1 . 3 8 × 1 0 − 2 3 , T = 2 7 3 . 1 5 + 3 0 , m = 6 . 6 4 6 × 1 0 − 2 7 , required probability is:
∫ E 1 E 2 p ( E ) d E = 0 . 0 0 1 6 8 4 8 9
# from math import * p = lambda E: (2.0/((1.380649 pow(10,-23)) 303 sqrt(pi))) * sqrt(E/((1.380649 pow(10,-23)) 303)) * pow(e,(-E/((1.380649 pow(10,-23)) 303))) ll = 3.323 pow(10,-27) 335 335 ul = 3.323 pow(10,-27) 345 345 def integrate(f, l, u, dx): ans = 0.0 x=l while x<u: ans += f(x) x += dx return ans dx integrate(p,ll,ul,pow(10,-29))
I wasn't sure whether this would work. I tried finding the integral mathematically but couldn't get anywhere. So I thought I'd try it some other day and just wrote the above out of curiosity to see if it'd work. It's not very accurate, but it works for upto 5 decimal places.
Indefinite integral (in exact terms) is k T e ⎝ ⎛ e k T erf ( k T e ) − π 2 e − k T e ⎠ ⎞
But I guess we only need a 3sf answer?
The probability can be determined by evaluating the integral of the probability density function p ( E ) with lower and upper limits E 1 and E 2 , respectively.
E 1 = 2 1 m v 1 2 = 2 1 ( 6 . 6 4 6 × 1 0 − 2 7 k g ) ( 3 3 5 m / s ) 2 = 3 . 7 2 9 2 × 1 0 − 2 2 J . E 2 = 2 1 m v 2 2 = 2 1 ( 6 . 6 4 6 × 1 0 − 2 7 k g ) ( 3 4 5 m / s ) 2 = 3 . 9 5 5 2 × 1 0 − 2 2 J . Evaluating the integral, we get E 1 ∫ E 2 p ( E ) d E = 0 . 0 0 1 6 8 3 .
E = 2 1 m v 2
Minimum energy is when v = 3 3 5 m / s : E = 3 . 7 2 9 2 3 6 7 5 ∗ 1 0 − 2 2
Maximum energy is when v = 3 4 5 m / s : E = 3 . 9 5 5 2 0 0 7 5 ∗ 1 0 − 2 2
Integrate probability distribution function from minimum energy to maximum energy. I did it numerically.
The lowest possible energy that the atom can have to be in our acceptable range happens with the lowest possible speed and equals E min = 2 1 m v min 2 = 3 . 7 2 9 ⋅ 1 0 − 2 2 J
Similarly the maximum happens with the maximum speed E max = 2 1 m v max 2 = 3 . 9 5 5 ⋅ 1 0 − 2 2 J
Now that we have the energy bounds, we must simply integrate the PDF function in the range to find the probability that the energy (and by extension the speed) is in the desired range: p = ∫ E min E max k T π 2 k T E exp ( − k t E ) d E where k = 1 . 3 8 ⋅ 1 0 − 2 3 J ⋅ K − 1 and T = 2 7 3 + 3 0 K.
Evaluating gives us that p ≈ 0 . 0 0 0 1 6 8 3 .
According to kinetic gas theory and since we are dealing with an ideal gas, the total energy for an atom equals its kinetic energy. Therefore
E
=
(
m
v
2
)
/
2
, where
m
is the mass of a helium atom. The probability that an atom has a speed near the speed of sound, which means a speed between
v
1
=
3
3
5
m
/
s
and
v
2
=
3
4
5
m
/
s
is the probability that the atom's energy is between
E
1
=
(
m
v
1
2
)
/
2
=
3
.
7
2
9
×
1
0
−
2
2
J
and
E
2
=
(
m
v
2
2
)
/
2
=
3
.
9
5
5
×
1
0
−
2
2
J
. Given that the probability as a function of energy for an atom in an ideal gas is
p
(
E
)
(where
p
(
E
)
is the function given in details and assumptions) the desired probability is the integral of
p
(
E
)
from
E
1
to
E
2
. By doing so we find that the probability is
0
.
0
0
1
6
8
, keeping 5 decimal digits. And now just a few notes:
The temperature we substitute in the
p
(
E
)
formula needs to be in S.I. units so
T
=
3
0
+
2
7
3
.
1
5
=
3
0
3
.
1
5
K
.
That
k
in
p
(
E
)
formula is Boltzmann's constant and it is equal to
1
.
3
8
0
6
4
8
8
×
1
0
−
2
3
m
2
k
g
(
s
)
−
2
K
−
1
.
The energy of an atom with speed 335 is a = 2 1 6 . 6 4 6 ⋅ 1 0 − 2 7 ⋅ 3 3 5 2 The energy of an atom with speed 345 is b = 2 1 6 . 6 4 6 ⋅ 1 0 − 2 7 ⋅ 3 4 5 2 If we integrate p ( E ) from a to b we get our answer, 0 . 0 0 1 6 8 3 .
To find the probability that that the particle has an energy in a certain range, we simply integrate the probability density function over that range.The probability density function is a function of energy, so we must find the kinetic energies associated with the given velocities v 1 = 3 3 5 m / s , v 2 = 3 4 5 m / s .
$$K 1 = \frac{1}{2}m v 1^2=3.746 \times 10^{-22} \ \mathrm{J}$$
$$K 2 = \frac{1}{2}m v 2^2=3.973 \times 10^{-22} \ \mathrm{J}$$
So then the probability P is
$$P=\int {K 1}^{K_2} {\frac{2}{kT\sqrt{\pi}}\sqrt{\frac{E}{kT}}\exp{\left(-\frac{E}{kT}\right)}} \ dE$$
which (taking k = 1 . 3 8 × 1 0 − 2 3 J / K and T = 3 0 3 K ) can be evaluated numerically to be
$$ \fbox{0.001697}. $$
The probability for an atom to have a speed in a given range is P ( v 1 , v 2 ) = ∫ v 1 v 2 p ( v ) d v . As the speed can be bijectively mapped to the kinetic energy, we can also calculate the sought probability as P ( v 1 , v 2 ) = P ( E 1 , E 2 ) = ∫ E 1 E 2 p ( E ) d E .
On the other hand, energies are expressed through speeds as E = 2 1 m v 2 . Thus P ( E 1 , E 2 ) = k T π 2 k T E ∫ E 1 E 2 exp ( − E / k T ) d E , which, substituting the numbers, gives us P ( 3 3 5 , 3 4 5 ) = 0 . 1 6 8 % .
Problem Loading...
Note Loading...
Set Loading...
First, it is important to note two things:
Ideal gas has no internal interactions between particles. This means that all of the energy in the particle is kinetic energy.
Now, we will focus on energy range. Since the target speed is anywhere from 3 3 5 m/s to 3 4 5 m/s, we will have energy levels from
2 1 ( 6 . 6 4 6 × 1 0 − 2 7 ) ( 3 3 5 ) 2 ≈ 3 . 7 2 9 × 1 0 − 2 2 J
to
2 1 ( 6 . 6 4 6 × 1 0 − 2 7 ) ( 3 4 5 ) 2 ≈ 3 . 9 5 5 × 1 0 − 2 2 J
Next, we can integrate the probability density function, p ( E ) , from the lower energy level to the higher energy level.
∫ 3 . 7 2 9 × 1 0 − 2 2 3 . 9 5 5 × 1 0 − 2 2 ( k T π 2 × k T E × e − k T E ) d E
Solving this integral, we find the probability to be 0 . 0 0 1 6 8 3
Relevant Quantities
k is the Boltzmann constant or 1 . 3 8 × 1 0 − 2 3 m 2 k g s − 2 K − 1
T is the temperature in Kelvins . In this case, it is 3 0 3 . 1 5 ∘ K