統計熱力学入門 第2章 §2.6 ★★ 極小原理の第二の応用。反応が「途中で止まる」理由は混合エントロピーにある。
化学ポテンシャルの役割
多成分系では、成分 i i i の粒子数 N i N_i N i ごとに化学ポテンシャルμ i = ( ∂ G / ∂ N i ) T , p , N j ≠ i \mu_i=(\partial G/\partial N_i)_{T,p,N_{j\ne i}} μ i = ( ∂ G / ∂ N i ) T , p , N j = i が定義される。μ i \mu_i μ i は「T , p T,p T , p を保ったまま成分 i i i の粒子を1個系に足すときの G G G の増分」であり、粒子の移動先を選ぶ判定量になる。粒子は μ i \mu_i μ i の高い場所から低い場所へ移ると G G G が下がるから、粒子のやり取りが許されるあらゆる平衡(相平衡・化学平衡・膜平衡)は μ \mu μ の釣り合いで書ける 。前節の相共存条件(式(2.9) )はその一例だった。本節は化学反応へ適用する。
反応の平衡条件
化学反応を量論係数 ν i \nu_i ν i で ∑ i ν i A i = 0 \sum_i\nu_i\mathrm A_i=0 ∑ i ν i A i = 0 と書く(生成物側を正にとる。例:N 2 O 4 ⇌ 2 N O 2 \mathrm{N_2O_4}\rightleftharpoons2\,\mathrm{NO_2} N 2 O 4 ⇌ 2 N O 2 ならν N 2 O 4 = − 1 \nu_{\mathrm{N_2O_4}}=-1 ν N 2 O 4 = − 1 、ν N O 2 = + 2 \nu_{\mathrm{NO_2}}=+2 ν N O 2 = + 2 )。反応が微小量 d ξ \mathrm{d}\xi d ξ 進むと各成分は d N i = ν i d ξ \mathrm{d} N_i=\nu_i\,\mathrm{d}\xi d N i = ν i d ξ ずつ変わる(ξ \xi ξ は反応進行度 )。T , p T,p T , p 固定での G G G の変化は
d G = ∑ i μ i d N i = ( ∑ i ν i μ i ) d ξ \mathrm{d} G = \sum_i\mu_i\,\mathrm{d} N_i = \Bigl(\sum_i\nu_i\,\mu_i\Bigr)\mathrm{d}\xi d G = i ∑ μ i d N i = ( i ∑ ν i μ i ) d ξ
だから、G G G 最小の条件として
title=★★★ 化学平衡の条件
∑ i ν i μ i = 0 (2.13) \boxed{\;\sum_i\nu_i\,\mu_i = 0\;}\tag{2.13} i ∑ ν i μ i = 0 ( 2.13 )
反応式の「左辺の μ \mu μ の合計=右辺の μ \mu μ の合計」(量論係数の重みつき)。相平衡の μ 1 = μ 2 \mu_1=\mu_2 μ 1 = μ 2 と同じ形の、μ \mu μ の釣り合いである。
理想気体混合と質量作用の法則
条件(2.13) を実用式にするには μ i \mu_i μ i の具体形が要る。理想気体の混合では、成分 i i i の化学ポテンシャルが分圧 p i p_i p i を用いて
μ i ( T , p i ) = μ i ∘ ( T ) + R T ln p i p ∘ (2.14) \mu_i(T,p_i) = \mu_i^\circ(T) + RT\ln\frac{p_i}{p^\circ}\tag{2.14} μ i ( T , p i ) = μ i ∘ ( T ) + R T ln p ∘ p i ( 2.14 )
と書ける(1モルあたり。p ∘ p^\circ p ∘ は基準圧力、μ i ∘ \mu_i^\circ μ i ∘ は基準圧力での純物質の値)。対数項の由来は混合エントロピーである。全圧 p p p の混合気体の中で成分 i i i は分圧 p i < p p_i<p p i < p しか感じておらず、その分だけ「薄まって」エントロピーの得をしている。式(2.14) を条件(2.13) に代入して整理すると
∑ i ν i μ i ∘ + R T ∑ i ν i ln p i p ∘ = 0 ⟺ ∏ i ( p i p ∘ ) ν i = exp ( − Δ r G ∘ R T ) , \sum_i\nu_i\,\mu_i^\circ + RT\sum_i\nu_i\ln\frac{p_i}{p^\circ}=0
\;\Longleftrightarrow\;
\prod_i\Bigl(\frac{p_i}{p^\circ}\Bigr)^{\nu_i}
= \exp\Bigl(-\frac{\Delta_{\mathrm r}G^\circ}{RT}\Bigr), i ∑ ν i μ i ∘ + R T i ∑ ν i ln p ∘ p i = 0 ⟺ i ∏ ( p ∘ p i ) ν i = exp ( − R T Δ r G ∘ ) ,
ここで Δ r G ∘ ≡ ∑ i ν i μ i ∘ \Delta_{\mathrm r}G^\circ\equiv\sum_i\nu_i\mu_i^\circ Δ r G ∘ ≡ ∑ i ν i μ i ∘ は標準反応 Gibbs エネルギーである。
title=★★ 質量作用の法則
K p ( T ) ≡ ∏ i ( p i p ∘ ) ν i = exp ( − Δ r G ∘ ( T ) R T ) (2.15) \boxed{\;K_p(T) \equiv \prod_i\Bigl(\frac{p_i}{p^\circ}\Bigr)^{\nu_i}
= \exp\Bigl(-\frac{\Delta_{\mathrm r}G^\circ(T)}{RT}\Bigr)\;}\tag{2.15} K p ( T ) ≡ i ∏ ( p ∘ p i ) ν i = exp ( − R T Δ r G ∘ ( T ) ) ( 2.15 )
平衡での分圧の組合せは、温度だけで決まる定数 K p K_p K p に拘束される。
ここにも第2.4 節の釣り合いがそのまま再演されている。Δ r G ∘ \Delta_{\mathrm r}G^\circ Δ r G ∘ はエネルギー項の役(反応がどちら向きに「落ちたい」か)を、式(2.14) の対数項は混合エントロピー項の役を演じる。Δ r G ∘ > 0 \Delta_{\mathrm r}G^\circ>0 Δ r G ∘ > 0 (生成物側が損)でも反応は完全にゼロにはならない。結晶中の空孔の例とまったく同じ理由で、生成物がごく僅かのとき混合エントロピーの微係数が発散し、少量だけ作る方が必ず G G G の得になるからである。化学反応が「完全に右」でも「完全に左」でもない途中で止まるのは、混合エントロピーの仕業 であり、その止まる位置を温度で重みづけして指定するのが式(2.15) である。
例題(label=ex:no2)
四酸化二窒素の解離平衡オリジナル 無色の気体 N 2 O 4 \mathrm{N_2O_4} N 2 O 4 は赤褐色の N O 2 \mathrm{NO_2} N O 2 へ解離する:N 2 O 4 ⇌ 2 N O 2 \mathrm{N_2O_4}\rightleftharpoons2\,\mathrm{NO_2} N 2 O 4 ⇌ 2 N O 2 。温度 298.15 K 298.15\,\mathrm K 298.15 K でこの反応の標準反応 Gibbs エネルギーはΔ r G ∘ = + 4.73 k J / m o l \Delta_{\mathrm r}G^\circ=+4.73\,\mathrm{kJ/mol} Δ r G ∘ = + 4.73 kJ/mol である。はじめ純粋な N 2 O 4 \mathrm{N_2O_4} N 2 O 4 を 1 mol 用意し、全圧 p = p ∘ = 1 b a r p=p^\circ=1\,\mathrm{bar} p = p ∘ = 1 bar に保つ。解離した割合(解離度)を α \alpha α とする。
(1) 平衡定数 K p K_p K p を α \alpha α と p / p ∘ p/p^\circ p / p ∘ で表せ。
(2) 298.15 K 298.15\,\mathrm K 298.15 K 、p = p ∘ p=p^\circ p = p ∘ での α \alpha α を求めよ。
(3) 全圧を上げると α \alpha α はどう変わるか。式(2.15) から示し、Le Chatelier の原理と整合することを述べよ。
解答 (1) 物質量は N 2 O 4 \mathrm{N_2O_4} N 2 O 4 が 1 − α 1-\alpha 1 − α 、N O 2 \mathrm{NO_2} N O 2 が 2 α 2\alpha 2 α 、合計 1 + α 1+\alpha 1 + α mol。分圧はモル分率と全圧の積だから
p N 2 O 4 = 1 − α 1 + α p , p N O 2 = 2 α 1 + α p . p_{\mathrm{N_2O_4}} = \frac{1-\alpha}{1+\alpha}\,p,\qquad
p_{\mathrm{NO_2}} = \frac{2\alpha}{1+\alpha}\,p. p N 2 O 4 = 1 + α 1 − α p , p N O 2 = 1 + α 2 α p .
量論係数(N O 2 \mathrm{NO_2} N O 2 が + 2 +2 + 2 、N 2 O 4 \mathrm{N_2O_4} N 2 O 4 が − 1 -1 − 1 )を式(2.15) に入れて
K p = ( p N O 2 / p ∘ ) 2 p N 2 O 4 / p ∘ = 4 α 2 ( 1 + α ) ( 1 − α ) p p ∘ = 4 α 2 1 − α 2 p p ∘ . K_p = \frac{(p_{\mathrm{NO_2}}/p^\circ)^2}{p_{\mathrm{N_2O_4}}/p^\circ}
= \frac{4\alpha^2}{(1+\alpha)(1-\alpha)}\,\frac{p}{p^\circ}
= \frac{4\alpha^2}{1-\alpha^2}\,\frac{p}{p^\circ}. K p = p N 2 O 4 / p ∘ ( p N O 2 / p ∘ ) 2 = ( 1 + α ) ( 1 − α ) 4 α 2 p ∘ p = 1 − α 2 4 α 2 p ∘ p .
(2) まず平衡定数を求める:
K p = exp ( − 4730 8.314 × 298.15 ) = e − 1.908 = 0.148. K_p = \exp\Bigl(-\frac{4730}{8.314\times298.15}\Bigr) = \mathrm{e}^{-1.908} = 0.148. K p = exp ( − 8.314 × 298.15 4730 ) = e − 1.908 = 0.148.
p = p ∘ p=p^\circ p = p ∘ のとき 4 α 2 / ( 1 − α 2 ) = K p 4\alpha^2/(1-\alpha^2)=K_p 4 α 2 / ( 1 − α 2 ) = K p を α \alpha α について解くと
α = K p 4 + K p = 0.148 4.148 = 0.189. \alpha = \sqrt{\frac{K_p}{4+K_p}} = \sqrt{\frac{0.148}{4.148}} = 0.189. α = 4 + K p K p = 4.148 0.148 = 0.189.
約 19% が解離した状態で平衡になる。Δ r G ∘ > 0 \Delta_{\mathrm r}G^\circ>0 Δ r G ∘ > 0 なのに解離がゼロでない点、かつ K p < 1 K_p<1 K p < 1 なので解離しきりもしない点に、本文の「混合エントロピーが途中で止める」機構が現れている。G ( α ) G(\alpha) G ( α ) を直接数値最小化しても同じ α = 0.189 \alpha=0.189 α = 0.189 が得られることを検算 [9] で確認した。
(3) K p K_p K p は温度だけの定数だから、p p p を増やすと積 4 α 2 1 − α 2 p p ∘ \dfrac{4\alpha^2}{1-\alpha^2}\dfrac{p}{p^\circ} 1 − α 2 4 α 2 p ∘ p を一定に保つために 4 α 2 1 − α 2 \dfrac{4\alpha^2}{1-\alpha^2} 1 − α 2 4 α 2 が減る、すなわち α \alpha α が減らねばならない。解離は分子数を 1 個から 2 個に増やす反応だから、「圧縮すると分子数を減らす向き(解離を戻す向き)へずれる」という Le Chatelier の原理の定量版になっている。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com