Calculate the probability of a molecule of N2 being in the ground vibrational state at 298 K.
The probability that a system occupies a given state at an instant of time and a specific temperature is given by the Boltzmann distribution.
P
i
=
exp
(
−
E
i
k
B
T
)
∑
j
exp
(
−
E
j
k
B
T
)
{\displaystyle P_{i}={\frac {\exp \left({\frac {-E_{i}}{k_{B}T}}\right)}{\sum _{j}\exp \left({\frac {-E_{j}}{k_{B}T}}\right)}}}
P
i
=
exp
(
−
E
i
k
B
T
)
Q
{\displaystyle P_{i}={\frac {\exp \left({\frac {-E_{i}}{k_{B}T}}\right)}{Q}}}
where:
i is the energy of the specific state, i, of interestkB is Boltzmann's constant, which equals
1.3806
×
10
−
34
{\displaystyle 1.3806\times 10^{-34}}
-1
T is the temperature in Kelvin The denominator of this function is known as the partition function, Q, which corresponds to the total number of accessible states of the molecule.
The closed form of the molecular vibrational partition function is given by:
q
v
i
b
=
1
1
−
e
−
h
ν
/
k
B
T
{\displaystyle q_{vib}={\frac {1}{1-e^{-h\nu /k_{B}T}}}}
where:
ν
{\displaystyle \nu }
2 in s-1 h is Planck's constant, which is
6.62607
×
10
−
34
{\displaystyle 6.62607\times 10^{-34}}
This is equivalent to Q since only the vibrational energy states are of interest and there is only one molecule of N2 . The equation for determining the partition function Q, from molecular partition functions, q, is given by:
Q
=
q
N
N
!
{\displaystyle Q={\frac {q^{N}}{N!}}}
where:
N is the number of molecules The fundamental vibrational frequency of N2 in wavenumbers,
ν
~
{\displaystyle {\tilde {\nu }}}
-1 [ 1]
The fundamental vibrational frequency in s-1 is given by:
ν
=
ν
~
×
c
{\displaystyle \nu ={\tilde {\nu }}\times c}
where
c is the speed of light, which is
2.9979
×
10
10
{\displaystyle 2.9979\times 10^{10}}
For N2 ,
ν
{\displaystyle \nu }
-1 ) \times (2.9979 \times 10^{10| cm/s) = 7.0708 \times 10^{13}</math>
For N2 at 298 K,
q
v
=
(
1
1
−
e
−
(
6.62607
×
10
−
34
J
s
×
7.0708
×
10
13
s
−
1
)
/
(
1.3806
×
10
−
23
J
K
−
1
×
298
K
)
)
=
1.000011333
{\displaystyle q_{v}=\left({\frac {1}{1-e^{-(6.62607\times 10^{-34}Js\times 7.0708\times 10^{13}s^{-1})/(1.3806\times 10^{-23}JK^{-1}\times 298K)}}}\right)=1.000011333}
The vibrational energy levels follow that of a quantum mechanical harmonic oscillator. The energy levels are represented by:
E
n
=
h
ν
(
n
+
1
2
)
{\displaystyle E_{n}=h\nu (n+{\frac {1}{2}})}
where:
n is the quantum vibrational number, which equals 0, 1, 2,... For the ground state (n=0), the energy becomes:
E
0
=
1
2
h
ν
{\displaystyle E_{0}={\frac {1}{2}}h\nu }
Since the vibrational zero point energy is not zero, the energy levels are defined relative to the n=0 level. This is used in the molecular partition function above and therefore, the ground state is regarded as having zero energy.
For N2 the probability of being in the ground state at 298K is:
P
0
=
e
−
E
0
/
k
B
T
q
v
{\displaystyle P_{0}={\frac {e^{-E_{0}/k_{B}T}}{q_{v}}}}
P
0
=
e
(
0
J
)
/
(
1.3806
×
10
−
23
×
298
K
)
1.000011333
{\displaystyle P_{0}={\frac {e^{(0J)/(1.3806\times 10^{-23}\times 298K)}}{1.000011333}}}
P
0
=
0.999988667
{\displaystyle P_{0}=0.999988667}
This means that at room temperature, the probability of a molecule of N2 being in the ground vibrational state is 99.9988667%.
↑ Lide, D. R., (84th ed.). (2003-2004). Handbook of Chemistry and Physics. pg.9-85.
Derive an equation for the population of rotational state i of a linear diatomic. Make a bar graph of the distribution of rotational states of N2 at 298 K.
1. Equation for the population of rotational state i
For a diatomic molecule, it can be approximated as a rigid rotor.
Solving Schrödinger equation of rigid rotors gives the energy levels of the molecule at state J:
E
(
J
)
=
B
J
(
J
+
1
)
{\displaystyle E(J)=BJ(J+1)}
where
J
{\displaystyle J}
B
{\displaystyle B}
-1 .
B
=
h
8
π
2
c
I
{\displaystyle B={\frac {h}{8\pi ^{2}cI}}}
where
h
{\displaystyle h}
Planck constant ;
c
{\displaystyle c}
I
{\displaystyle I}
I
=
μ
r
2
{\displaystyle I=\mu r^{2}}
where
μ
{\displaystyle \mu }
r
{\displaystyle r}
By the Maxwell–Boltzmann_distribution , the population of rotational state i comparing to ground state is:
N
i
N
0
=
g
i
g
0
(
e
−
E
i
/
k
B
T
e
−
E
0
/
k
B
T
)
{\displaystyle {\frac {N_{i}}{N_{0}}}={\frac {g_{i}}{g_{0}}}\left({\frac {e^{-E_{i}/k_{B}T}}{e^{-E_{0}/k_{B}T}}}\right)}
where
g
{\displaystyle g}
E
{\displaystyle E}
k
B
{\displaystyle k_{B}}
Boltzmann constant and
T
{\displaystyle T}
Substitute
E
i
=
h
c
B
i
(
i
+
1
)
{\displaystyle E_{i}=hcBi(i+1)}
E
0
=
h
c
B
0
(
0
+
1
)
=
0
{\displaystyle E_{0}=hcB0(0+1)=0}
g
i
=
2
i
+
1
{\displaystyle g_{i}=2i+1}
N
i
N
0
=
(
2
i
+
1
)
(
e
−
h
c
B
i
(
i
+
1
)
/
k
B
T
e
0
/
k
B
T
)
=
(
2
i
+
1
)
e
−
h
c
B
i
(
i
+
1
)
/
k
B
T
{\displaystyle {\frac {N_{i}}{N_{0}}}=(2i+1)\left({\frac {e^{-hcBi(i+1)/k_{B}T}}{e^{0/k_{B}T}}}\right)=(2i+1)e^{-hcBi(i+1)/k_{B}T}}
2. Distribution of rotational states of N2
edit
For N2 , the population of state i is:
N
i
N
0
=
(
2
i
+
1
)
e
−
h
c
B
i
(
i
+
1
)
/
k
B
T
{\displaystyle {\frac {N_{i}}{N_{0}}}=(2i+1)e^{-hcBi(i+1)/k_{B}T}}
Combine the constant part, define constant
a
=
h
c
B
/
k
B
T
=
h
2
8
π
2
k
B
T
μ
r
2
{\displaystyle a=hcB/k_{B}T={\frac {h^{2}}{8\pi ^{2}k_{B}T\mu r^{2}}}}
The reduced mass μ of Nitrogen is 7.00D a=1.16×10-26 kg and the bond length r of N2 is 110 pm = 1.10×10-10 m.
Substitute T = 298 K,
a
=
h
2
8
π
2
k
B
T
μ
r
2
=
9.60
×
10
−
3
{\displaystyle a={\frac {h^{2}}{8\pi ^{2}k_{B}T\mu r^{2}}}=9.60\times 10^{-3}}
By sum over states, qrot =1/a = 104, so the probability of N2 occupying the ground vibrational state at 298 K is 1/104 = 9.60×10-3
The bar graph of the distribution of rotational states of N2 at 298 K is:
Distribution of rotational states of nitrogen gas at 298 K.
Estimate the number of translational states that are available to a molecule of N2 in 1 m3 container at 298 K.
The equation that is used to determine translational states of the molecule of N2 at 298 K is shown below.
q
t
r
a
n
s
(
V
,
T
)
=
(
2
π
m
k
B
T
h
2
)
3
2
{\displaystyle q_{trans}(V,T)=\left({\frac {2\pi mk_{B}T}{h^{2}}}\right)^{\frac {3}{2}}}
V
{\displaystyle V}
Where
q
t
r
a
n
s
{\displaystyle q_{trans}}
m
{\displaystyle m}
k
B
{\displaystyle k_{B}}
(
1.38064852
×
10
−
23
J
×
K
−
1
)
{\displaystyle (1.38064852\times 10^{-23}J\times K^{-1})}
T
{\displaystyle T}
(
K
)
{\displaystyle (K)}
h
2
{\displaystyle h^{2}}
(
6.626
×
10
−
34
J
×
s
)
{\displaystyle (6.626\times 10^{-34}J\times s)}
V
{\displaystyle V}
m
3
{\displaystyle m^{3}}
The steps for solving this problem is shown below:
m =
28.0102
{\displaystyle 28.0102}
m =
(
28.0102
a
m
u
)
(
1.6605
×
10
−
27
k
g
)
{\displaystyle (28.0102amu)(1.6605\times 10^{-27}kg)}
m =
4.6512
×
10
−
26
k
g
{\displaystyle 4.6512\times 10^{-26}kg}
q
t
r
a
n
s
(
V
,
T
)
=
(
2
π
m
k
B
T
h
2
)
3
/
2
V
{\displaystyle q_{trans}(V,T)=\left({\frac {2\pi mk_{B}T}{h^{2}}}\right)^{3/2}V}
q
t
r
a
n
s
(
V
,
T
)
=
(
2
×
3.14159
×
(
4.6512
×
10
−
26
k
g
)
×
(
1.38064852
×
10
−
23
J
/
K
)
×
298
K
(
6.626
×
10
−
34
J
×
s
)
2
)
3
/
2
(
1.000
m
3
)
{\displaystyle q_{trans}(V,T)=\left({\frac {2\times 3.14159\times (4.6512\times 10^{-26}kg)\times (1.38064852\times 10^{-23}J/K)\times 298K}{{(6.626\times 10^{-34}J\times s})^{2}}}\right)^{3/2}(1.000m^{3})}
q
t
r
a
n
s
(
V
,
T
)
=
1.4322
×
10
32
{\displaystyle q_{trans}(V,T)=1.4322\times 10^{32}}
Therefore, there should be
1.4322
×
10
32
{\displaystyle 1.4322\times 10^{32}}
2 at 298 K.
Calculate the DeBroglie wavelength, rotational temperature, and vibrational temperature of N2 and Cl2 at 298 K.
The DeBroglie Wavelength
edit
Λ
=
(
2
π
m
k
B
T
h
2
)
−
1
/
2
{\displaystyle \Lambda =\left({\frac {2{\pi }mk_{B}T}{h^{2}}}\right)^{-1/2}}
where:
m is the mass of the molecule.Boltzmann's constant kB = 1.3806488×10-23 J K-1
Planck's constant h = 6.62606957×10-34 J s For N2 at 298K,
m
=
2
×
14.0067
u
{\displaystyle m=2\times 14.0067u}
m
=
28.0134
u
×
1.66054
×
10
−
27
k
g
u
{\displaystyle m=28.0134u\times {\frac {1.66054\times 10^{-27}kg}{u}}}
m
=
4.65173712
×
10
−
26
k
g
{\displaystyle m=4.65173712\times 10^{-26}kg}
Λ
=
(
2
π
(
4.65173712
×
10
−
26
k
g
)
(
1.3806488
×
10
−
23
J
K
−
1
)
(
298
K
)
(
6.62606957
×
10
−
34
J
s
)
2
)
−
1
/
2
{\displaystyle \Lambda =\left({\frac {2{\pi }(4.65173712\times 10^{-26}kg)(1.3806488\times 10^{-23}JK^{-1})(298K)}{(6.62606957\times 10^{-34}Js)^{2}}}\right)^{-1/2}}
Λ
=
(
1.20252611
×
10
−
45
k
g
J
4.39047979
×
10
−
67
J
2
s
2
)
−
1
/
2
{\displaystyle \Lambda =\left({\frac {1.20252611\times 10^{-45}kgJ}{4.39047979\times 10^{-67}J^{2}s^{2}}}\right)^{-1/2}}
Λ
=
(
2.738940097
×
10
−
21
k
g
J
−
1
s
−
2
)
−
1
/
2
{\displaystyle \Lambda =(2.738940097\times 10^{-21}kgJ^{-1}s^{-2})^{-1/2}}
Λ
=
1.9107714
×
10
−
11
m
{\displaystyle \Lambda =1.9107714\times 10^{-11}m}
For Cl2 at 298K,
m
=
2
×
35.453
u
{\displaystyle m=2\times 35.453u}
m
=
70.906
u
×
1.66054
×
10
−
27
k
g
u
{\displaystyle m=70.906u\times {\frac {1.66054\times 10^{-27}kg}{u}}}
m
=
1.17742249
×
10
−
25
k
g
{\displaystyle m=1.17742249\times 10^{-25}kg}
Λ
=
(
2
π
(
1.17742249
×
10
−
25
k
g
)
(
1.3806488
×
10
−
23
J
K
−
1
)
(
298
K
)
(
6.62606957
×
10
−
34
J
s
)
2
)
−
1
/
2
{\displaystyle \Lambda =\left({\frac {2{\pi }(1.17742249\times 10^{-25}kg)(1.3806488\times 10^{-23}JK^{-1})(298K)}{(6.62606957\times 10^{-34}Js)^{2}}}\right)^{-1/2}}
Λ
=
(
3.04376893
×
10
−
45
k
g
J
4.39047979
×
10
−
67
J
2
s
2
)
−
1
/
2
{\displaystyle \Lambda =\left({\frac {3.04376893\times 10^{-45}kgJ}{4.39047979\times 10^{-67}J^{2}s^{2}}}\right)^{-1/2}}
Λ
=
(
6.932656726
×
10
21
k
g
J
−
1
s
−
2
)
−
1
/
2
{\displaystyle \Lambda =(6.932656726\times 10^{21}kgJ^{-1}s^{-2})^{-1/2}}
Λ
=
1.20101976
×
10
−
11
m
{\displaystyle \Lambda =1.20101976\times 10^{-11}m}
Unit Conversion for the DeBroglie Wavelength:
(
k
g
J
−
1
s
−
2
)
−
1
/
2
=
(
k
g
k
g
m
2
s
−
2
s
2
)
−
1
/
2
{\displaystyle (kgJ^{-1}s^{-2})^{-1/2}=\left({\frac {kg}{kgm^{2}s^{-2}s^{2}}}\right)^{-1/2}}
(
k
g
J
−
1
s
−
2
)
−
1
/
2
=
(
m
−
2
)
−
1
/
2
{\displaystyle (kgJ^{-1}s^{-2})^{-1/2}=(m^{-2})^{-1/2}}
(
k
g
J
−
1
s
−
2
)
−
1
/
2
=
m
{\displaystyle (kgJ^{-1}s^{-2})^{-1/2}=m}
The Rotational Temperature
edit
Θ
r
=
ℏ
2
2
k
B
μ
r
e
2
{\displaystyle \Theta _{r}={\frac {\hbar ^{2}}{2k_{B}\mu r_{e}^{2}}}}
where:
Planck's constant
ℏ
=
h
2
π
=
6.62606957
×
10
−
34
J
s
2
π
=
1.054571726
×
10
−
34
J
s
{\displaystyle \hbar ={\frac {h}{2\pi }}={\frac {6.62606957\times 10^{-34}Js}{2\pi }}=1.054571726\times 10^{-34}Js}
μ is the reduced mass.re is the bond length between two atoms in a moleculeFor N2 ,
re = 1.09769Å = 1.09769×10-10 m [ 1]
μ
=
m
N
m
N
m
N
+
m
N
=
m
N
2
2
m
N
=
m
N
2
=
14.00307400529
u
2
=
7.001537003
u
{\displaystyle \mu ={\frac {m_{N}m_{N}}{m_{N}+m_{N}}}={\frac {m_{N}^{2}}{2m_{N}}}={\frac {m_{N}}{2}}={\frac {14.00307400529u}{2}}=7.001537003u}
μ
=
7.001537003
u
×
1.66054
×
10
−
27
k
g
u
{\displaystyle \mu =7.001537003u\times {\frac {1.66054\times 10^{-27}kg}{u}}}
μ
=
1.16263323
×
10
−
26
k
g
{\displaystyle \mu =1.16263323\times 10^{-26}kg}
Θ
r
=
(
1.054571726
×
10
−
34
J
s
)
2
2
(
1.3806488
×
10
−
23
J
K
−
1
)
(
1.16263323
×
10
−
26
k
g
)
(
1.09769
×
10
−
10
m
)
2
{\displaystyle \Theta _{r}={\frac {(1.054571726\times 10^{-34}Js)^{2}}{2(1.3806488\times 10^{-23}JK^{-1})(1.16263323\times 10^{-26}kg)(1.09769\times 10^{-10}m)^{2}}}}
Θ
r
=
1.11212153
×
10
−
68
J
2
s
2
3.86825738
×
10
−
69
J
K
−
1
k
g
m
2
{\displaystyle \Theta _{r}={\frac {1.11212153\times 10^{-68}J^{2}s^{2}}{3.86825738\times 10^{-69}JK^{-1}kgm^{2}}}}
Θ
r
=
2.874993624
J
s
2
K
k
g
−
1
m
−
2
{\displaystyle \Theta _{r}=2.874993624Js^{2}Kkg^{-1}m^{-2}}
Θ
r
=
2.874993624
K
{\displaystyle \Theta _{r}=2.874993624K}
For Cl2 ,
re = 1.988Å = 1.988×10-10 m [ 2]
μ
=
m
C
l
m
C
l
m
C
l
+
m
C
l
=
m
C
l
2
2
m
C
l
=
m
C
l
2
=
34.968852714
u
2
=
17.48442636
u
{\displaystyle \mu ={\frac {m_{Cl}m_{Cl}}{m_{Cl}+m_{Cl}}}={\frac {m_{Cl}^{2}}{2m_{Cl}}}={\frac {m_{Cl}}{2}}={\frac {34.968852714u}{2}}=17.48442636u}
μ
=
17.4842636
u
×
1.66054
×
10
−
27
k
g
u
{\displaystyle \mu =17.4842636u\times {\frac {1.66054\times 10^{-27}kg}{u}}}
μ
=
2.90335893
×
10
−
26
k
g
{\displaystyle \mu =2.90335893\times 10^{-26}kg}
Θ
r
=
(
1.054571726
×
10
−
34
J
s
)
2
2
(
1.3806488
×
10
−
23
J
K
−
1
)
(
2.90335893
×
10
−
26
k
g
)
(
1.988
×
10
−
10
m
)
2
{\displaystyle \Theta _{r}={\frac {(1.054571726\times 10^{-34}Js)^{2}}{2(1.3806488\times 10^{-23}JK^{-1})(2.90335893\times 10^{-26}kg)(1.988\times 10^{-10}m)^{2}}}}
Θ
r
=
1.11212153
×
10
−
68
J
2
s
2
3.16844888
×
10
−
68
J
K
−
1
k
g
m
2
{\displaystyle \Theta _{r}={\frac {1.11212153\times 10^{-68}J^{2}s^{2}}{3.16844888\times 10^{-68}JK^{-1}kgm^{2}}}}
Θ
r
=
0.350998729
J
s
2
K
k
g
−
1
m
−
2
{\displaystyle \Theta _{r}=0.350998729Js^{2}Kkg^{-1}m^{-2}}
Θ
r
=
0.350998729
K
{\displaystyle \Theta _{r}=0.350998729K}
Unit Conversion for the Rotational Temperature:
J
s
2
K
k
g
−
1
m
−
2
=
J
(
k
g
−
1
m
−
2
s
2
)
K
{\displaystyle Js^{2}Kkg^{-1}m^{-2}=J(kg^{-1}m^{-2}s^{2})K}
J
s
2
K
k
g
−
1
m
−
2
=
J
J
−
1
K
{\displaystyle Js^{2}Kkg^{-1}m^{-2}=JJ^{-1}K}
J
s
2
K
k
g
−
1
m
−
2
=
K
{\displaystyle Js^{2}Kkg^{-1}m^{-2}=K}
The Vibrational Temperature
edit
Θ
v
=
h
ν
k
B
{\displaystyle \Theta _{v}={\frac {h\nu }{k_{B}}}}
Θ
v
=
h
ν
~
c
k
B
{\displaystyle \Theta _{v}={\frac {h{\tilde {\nu }}c}{k_{B}}}}
where:
ν
~
{\displaystyle {\tilde {\nu }}}
The speed of light c = 2.99792458×1010 cm s-1 For N2 ,
ν
~
=
2358.57
c
m
−
1
{\displaystyle {\tilde {\nu }}=2358.57cm^{-1}}
[ 1]
Θ
v
=
(
6.62606957
×
10
−
34
J
s
)
(
2358.57
c
m
−
1
)
(
2.99792458
×
10
10
c
m
s
−
1
)
1.3806488
×
10
−
23
J
K
−
1
{\displaystyle \Theta _{v}={\frac {(6.62606957\times 10^{-34}Js)(2358.57cm^{-1})(2.99792458\times 10^{10}cms^{-1})}{1.3806488\times 10^{-23}JK^{-1}}}}
Θ
v
=
4.6851712
×
10
−
20
J
1.3806488
×
10
−
23
J
K
−
1
{\displaystyle \Theta _{v}={\frac {4.6851712\times 10^{-20}J}{1.3806488\times 10^{-23}JK^{-1}}}}
Θ
v
=
3393.456175
K
{\displaystyle \Theta _{v}=3393.456175K}
For Cl2 ,
ν
~
=
559.7
c
m
−
1
{\displaystyle {\tilde {\nu }}=559.7cm^{-1}}
[ 2]
Θ
v
=
(
6.62606957
×
10
−
34
J
s
)
(
559.7
c
m
−
1
)
(
2.99792458
×
10
10
c
m
s
−
1
)
1.3806488
×
10
−
23
J
K
−
1
{\displaystyle \Theta _{v}={\frac {(6.62606957\times 10^{-34}Js)(559.7cm^{-1})(2.99792458\times 10^{10}cms^{-1})}{1.3806488\times 10^{-23}JK^{-1}}}}
Θ
v
=
1.11181365
×
10
−
20
J
1.3806488
×
10
−
23
J
K
−
1
{\displaystyle \Theta _{v}={\frac {1.11181365\times 10^{-20}J}{1.3806488\times 10^{-23}JK^{-1}}}}
Θ
v
=
805.2834645
K
{\displaystyle \Theta _{v}=805.2834645K}
↑ a b Lide, D. R., (84th ed.). (2003-2004). Handbook of Chemistry and Physics. CRC Press. pg.9-85.
↑ a b Lide, D. R., (84th ed.). (2003-2004). Handbook of Chemistry and Physics. CRC Press. pg.9-83.
Molecule
N2
Cl2
Λ
{\displaystyle \Lambda }
1.91×10-11 m
1.20×10-11 m
Θ
r
{\displaystyle \Theta _{r}}
2.87 K
0.351 K
Θ
v
{\displaystyle \Theta _{v}}
3.39×103 K
805 K
At what temperature would the probability for N2 being the ground vibrational state be reduced to 50%?
50% of the N2 molecules are in the ground state, and 50% are in the excited state. So, the population can be expressed in the following equation:
P
0
=
0.5
{\displaystyle P_{0}=0.5}
The value for the
ν
~
{\displaystyle {\tilde {\nu }}}
2 molecule, taken from the NIST Webbook, is 2358.57 cm-1
This can be converted into the fundamental vibrational frequency by the relation,
ν
=
ν
~
×
c
{\displaystyle \nu ={\tilde {\nu }}\times c}
where:
ν
~
{\displaystyle {\tilde {\nu }}}
-1 .
c
{\displaystyle c}
10 cm/s By knowing this relation, the fundamental vibrational frequency can be calculated.
ν
=
ν
~
×
c
=
2358.57
c
m
−
1
×
2.998
×
10
10
c
m
/
s
=
7.07099
×
10
13
s
−
1
{\displaystyle \nu ={\tilde {\nu }}\times c=2358.57cm^{-1}\times 2.998\times 10^{10}cm/s=7.07099\times 10^{13}s^{-1}}
The population of a system can be expressed by the following equation:
P
0
=
exp
(
−
E
0
k
B
T
)
q
v
i
b
{\displaystyle P_{0}={\frac {\exp \left({\frac {-E_{0}}{k_{B}T}}\right)}{q_{vib}}}}
where:
P
0
{\displaystyle P_{0}}
E
0
{\displaystyle E_{0}}
k
B
{\displaystyle k_{B}}
-23 J/K
T
{\displaystyle T}
q
v
i
b
{\displaystyle q_{vib}}
q
v
i
b
{\displaystyle q_{vib}}
q
v
i
b
=
1
1
−
exp
(
−
h
ν
k
B
T
)
{\displaystyle q_{vib}={\frac {1}{1-\exp \left({\frac {-h\nu }{k_{B}T}}\right)}}}
where:
h is the Planck constant, 6.626*10-34 J*s
Knowing these relations, we can solve for the temperature of the system when the molecule is reduced to 50% ground state, which is when 50% of the molecules are in the ground state and 50% of the molecules are excited.
P
0
=
exp
−
E
0
k
B
T
q
v
i
b
=
50
%
=
0.5
{\displaystyle P_{0}={\frac {\exp {\frac {-E_{0}}{k_{B}T}}}{q_{vib}}}=50\%=0.5}
q
v
i
b
=
exp
(
−
E
0
k
B
T
)
P
0
{\displaystyle q_{vib}={\frac {\exp \left({\frac {-E_{0}}{k_{B}T}}\right)}{P_{0}}}}
q
v
i
b
=
exp
(
0
)
0.5
=
2
{\displaystyle q_{vib}={\frac {\exp(0)}{0.5}}=2}
Knowing the value of
q
v
i
b
{\displaystyle q_{vib}}
q
v
i
b
=
1
1
−
e
−
h
ν
k
B
T
{\displaystyle q_{vib}={\frac {1}{1-e^{\frac {-h\nu }{k_{B}T}}}}}
=
2
{\displaystyle =2}
2
=
1
1
−
exp
−
h
ν
k
B
T
{\displaystyle 2={\frac {1}{1-\exp ^{\frac {-h\nu }{k_{B}T}}}}}
2
(
1
−
exp
−
h
ν
k
B
T
)
=
1
{\displaystyle 2(1-\exp ^{\frac {-h\nu }{k_{B}T}})=1}
2
−
2
exp
−
h
ν
k
B
T
=
1
{\displaystyle 2-2\exp ^{\frac {-h\nu }{k_{B}T}}=1}
−
2
exp
(
−
h
ν
k
B
T
)
=
−
1
{\displaystyle -2\exp \left({\frac {-h\nu }{k_{B}T}}\right)=-1}
exp
(
−
h
ν
k
B
T
)
=
0.5
{\displaystyle \exp \left({\frac {-h\nu }{k_{B}T}}\right)=0.5}
ln
(
exp
(
−
h
ν
k
B
T
)
=
ln
(
0.5
)
{\displaystyle \ln(\exp \left({\frac {-h\nu }{k_{B}T}}\right)=\ln(0.5)}
−
h
ν
k
B
T
=
−
0.6931472
{\displaystyle {\frac {-h\nu }{k_{B}T}}=-0.6931472}
T
=
h
ν
0.6931472
k
B
{\displaystyle T={\frac {h\nu }{0.6931472k_{B}}}}
Therefore, by substituting the Planck constant, fundamental vibrational frequency and the Boltzmann constant into the equation, we can obtain the temperature for when half of the molecules are in the ground state.
T
=
6.626
×
10
−
34
J
⋅
s
×
7.07099
×
10
13
s
−
1
/
(
0.6931472
×
1.38064853
×
10
−
23
J
/
K
)
=
4895.792994
K
{\displaystyle T=6.626\times 10^{-34}J\cdot s\times 7.07099\times 10^{13}s^{-}1/(0.6931472\times 1.38064853\times 10^{-23}J/K)=4895.792994K}
The temperature at which 50% of the N2 molecules are in the ground vibrational state is 4895 K.
