300 Kで熱平衡にある粒子について、基底準位に1状態、励起準位に同じエネルギーの3状態がある模型を考えます。励起準位にいる粒子は全体の何%で、一つの励起状態には何%いるでしょうか。 各状態のボルツマン重みを計算してから、四つの状態を合計した分母で割ります。
一つずつ比べると、高いエネルギーの状態は占有されにくくなります。それでも三状態をまとめた人数はどうなるか、状態ごとの割合と準位の合計を分けて答えます。

1.「同じエネルギーの状態が3つ」を数える
熱平衡にある、互いの相互作用を無視できる粒子の内部状態を、二つのエネルギー準位だけで表す仮想モデルを考えます。基底準位を0、励起準位を1とし、基底に属する状態は1つ、励起に属する状態は3つ。縮退度gは、一つのエネルギーに属する異なる状態の数で、g₀ = 1、g₁ = 3です。g₁ = 3は「粒子が3個いる」という意味ではありません。
このモデルではT₀=300 Kを基準に、モルあたりのエネルギー差をΔE_m = RT₀と与えました。T₀と温度Tの単位はK、気体定数RはJ mol⁻¹ K⁻¹、ΔE_mはJ mol⁻¹。エネルギー差は固定してあり、温度を変えるたびにRTへ設定し直すことはしません。まずT = T₀で計算します。

2.一つの状態の重みから、準位の人数の比へ
熱平衡のボルツマン分布では、基底のエネルギーを0として、一つの励起状態の重みはexp(−ΔE_m/RT)です。指数の中は同じJ mol⁻¹どうしを割るので無次元。T = T₀ならΔE_m/RT₀ = 1で、各励起状態の重みはe⁻¹ ≈ 0.367879となります。原典が粒子1個あたりの差ΔεとkBTで書く式も、ΔE_m = NAΔε、R = NAkBとすれば同じ比です。
同じエネルギーの三状態なので、それぞれの重みは等しくなります。縮退度3を掛けるのは、粒子を3個と決めたからではなく、同じ重みの項が三つあるからです。
N₀は基底準位にいる粒子の合計、N₁は励起準位の3状態にいる粒子の合計です。同じ温度、同じ全体の分布で比べるため、全体を規格化する共通の分母は比で消えます。
T = T₀では N₁/N₀ = 3e⁻¹ = 3/e
≈ 1.10364
一つの励起状態/基底状態 = e⁻¹ ≈ 0.367879
ペンタ高いエネルギーは占有されにくいのに、励起準位の合計が大きくなってよいの?
イコール一つずつなら励起状態は基底状態の約0.368倍です。でも同じエネルギーの励起状態が3つあるので、合計は3倍して約1.104倍になります。比較の単位を「状態」と「準位の合計」で分けましょう。
ペンタでは、高い方へ反転分布が起きているということ?
イコールここでは正の温度の熱平衡で、状態一つあたりの占有は高い方が少ないです。縮退した3状態の合計が大きいだけです。準位全体の人数をそのまま比べて、状態あたりの占有が逆転したと判断はできません。
3.比1.104を、全体に対する割合へ直す
N₁/N₀の分母は基底の人数です。全体に対する励起の割合P₁を知りたいなら、分母をN₀+N₁へ替えます。この模型の四つの内部状態の重みを全部足し、q = 1+3/eと書きます。求める割合の分母は全粒子数なので、各状態・各準位の重みをこの共通のqで割ります。P₀、P₁は単位のない割合で、和は1です。
P₀ = N₀/(N₀+N₁)= 1/q ≈ 0.47537
P₁ = N₁/(N₀+N₁)=(3/e)/q ≈ 0.52463
一つの励起状態の割合 = P₁/3 ≈ 0.17488
検算:P₀+P₁ = 1


ペンタN₁/N₀が1.104なら、励起の割合は110.4%?
イコールその比の分母は全体ではなく基底です。基底の重み1に対して励起の重みが約1.104です。全体は約2.104なので、励起の割合は1.104を2.104で割った約52.46%です。
ペンタ52.46%を3で割った17.49%も、励起の中での割合なの?
イコール今の17.49%の分母も全体です。励起の中だけを分母にすれば、同じ重みの3状態はそれぞれ3分の1です。全体のうち17.49%と、励起の中の3分の1は別の割合です。
元の重みに戻して確認します。丸める前の式を使えば、P₁/P₀ = [(3/e)/q]/[1/q] = 3/eで、励起準位と基底準位の人数比に一致します。一つの励起状態では(P₁/3)/P₀ = 1/e。三状態の合計が基底より多くても、各励起状態は基底より少ないままです。
さらにP₀+3(P₁/3) = (1+3/e)/q = 1。四つの状態の割合を足すと全体100%へ戻ります。図2の上段の三本を足したものが下段の励起の棒であることも、この式で確かめられます。
4.温度を下げるとき、固定するのはエネルギー差
このモデルではg₀、g₁とΔE_mを固定して温度だけ変えます。ΔE_m = RT₀を代入すると、一般の温度TではN₁/N₀ = 3 exp(−T₀/T)。Tが低くなると指数の負の大きさが増し、励起準位の合計も小さくなります。温度を変えた後も熱平衡に達したときの割合を計算します。
合計N₁がN₀を上回る条件は3 exp(−T₀/T)>1。両辺の自然対数を取るとln 3−T₀/T>0なので、T > T₀/ln 3 ≈ 0.91024 T₀、本例では約273.07 Kより高い温度です。境界で準位の合計は等しくなりますが、そのときも各励起状態の占有は基底状態の3分の1です。
5.自分で確かめる:T₀の半分へ冷やすと?
同じg₀ = 1、g₁ = 3と、固定したΔE_m = RT₀で、温度をT = T₀/2 = 150 Kにします。①N₁/N₀、②励起準位の全体に対する割合P₁、③一つの励起状態の全体に対する割合を求めてください。指数を−1のままにした答えは、何を変えてしまったことになるでしょうか。
解答 ΔE_m/RT = RT₀/{R(T₀/2)}= 2。①3e⁻² ≈ 0.40601。②P₁ = 3e⁻²/(1+3e⁻²)≈ 0.28877(28.88%)。③P₁/3 ≈ 0.09626(9.63%)。基底P₀ ≈ 0.71123と励起P₁の和は1です。指数を−1のままにすると、温度を半分にしたのにエネルギー差まで半分のRTへ替えた計算になり、指定の固定差を保っていません。
参考資料
MIT 5.62 Lecture 11:状態ごとのボルツマン重み、縮退度、準位の人数比(PDF物理ページ4〜6/講義ページ3〜5)


コメント