多成分系の熱平衡

Prerequisite

ここまで熱力学系に含まれる物質の種類は1種類としてきた.
たとえば酸素と窒素など,気体が複数種類ある場合の平衡状態はどのように記述できるだろうか.
この節では複数種の物質は化学反応しないとする.

まず理想気体を仮定して気体分子どうしの相互作用を無視する.
気体 1  と気体 2  がそれぞれ N_1,\,N_2  個含まれる系を場合を考える.
系の粒子数N=N_1+N_2
系の温度と体積を T,\,V  とすると内部エネルギーは加法性により,

で与えられる.
ここで c_i>1  は定数.
圧力は状態方程式により

で与えられる.

一方で気体 1  の状態方程式より

が成立する.
x_1=N_1/N  という因子は系に含まれる気体 1  の割合である.
P_1  を気体 1  分圧 (partial pressure) という.

理想気体のエントロピーも加法性により

で与えられる.
ここで S_{i0}  は定数.
U_i=c_iNk_{\mathrm{B}}T  を代入して整理すると

となる.
ここで S'_{i0}=S_{i0}+c_ik_{\mathrm{B}}\ln(c_ik_{\mathrm{B}})  とおいた.

混合過程

2種類の異なる気体の系と1種類のみの系の違いを明確にしておこう.
部分系 1  2  を考えて,それぞれに気体 1  2  N_1,\,N_2  ずつ入っている.
系の間には粒子は通さないが可動かつ透熱の壁があって,温度 T  ,圧力 P  の平衡状態にあるとする.
この系全体のエントロピーは加法性により

で与えられる.
ただし U=U_1+U_2  V=V_1+V_2

間の壁を取り除いて2つの気体を混合する過程を考える.
この過程では仕事が 0  であり断熱過程なので Q=0
よって \delta U=0  で十分時間が経つと温度と圧力は T,\,P  のままの新たな平衡状態に達する.
この系のエントロピーは

となる.
遷移前後のちがいは体積に添字があるかないかである.
混合過程におけるエントロピーの変化分は

となる.
温度と圧力は一定だから状態方程式より

が成り立っていることから

混合エントロピー

と書き換えられる.
これは2種類の気体の混合によるエントロピーの増分であり混合エントロピー (entropy of mixing) とよばれる.

実際部分系の両方に同じ気体を入れた場合,遷移後のエントロピーは同じ函数 S(U,V,N)  によって計算できるから,

でありエントロピー函数の示量性によってこの2つは等しいことがわかる.
よって \delta S = 0  で混合によるエントロピーの増加はない.
エントロピー増大の法則から \delta S>0  の断熱過程は不可逆である.
同種の気体の混合では再び壁を挿入すれば元の状態に戻せるが,異種の場合はそれは(断熱的には)できない.

一般の気体(液体含む)に対して p  種類の物質を含む多成分系のエントロピーを

とおく.
一般には気体分子どうしが相互作用するので単なる和で書くことはできない:

微分形式では

とかける.ここで

は気体 i  の化学ポテンシャルである.

多成分系の平衡状態

部分系 1  2  にそれぞれ p  種類の気体が入っている系を考える.
部分系どうしは壁によっていくつかの物理量が隔絶されている.
この壁は断熱/透熱や可動/不動などを自由に選べた.
ここで新たな壁として特定の種類の気体だけ透過するような壁を用意しよう.
いま透熱・可動かつ気体 i  だけを透過する壁とすると,エントロピー最大の原理から

が導かれる.
壁をさらに工夫して p  種類のうち q  種類の気体 i_1,\cdots,i_q  のみを透過するようにした場合は,温度・圧力の平衡に加えて

が成り立つ.
つまり平衡状態ではやりとりのある気体の化学ポテンシャルが等しくなる.

Gibbsの自由エネルギーが G(T,P,N)=NG(T,P,1)=N\mu(T,P)  で与えられるように,1種類の気体では温度と圧力が定まると,化学ポテンシャルも1つに定まり,Gibbsの自由エネルギーも決まる.
2種類の気体のある2成分系ではまずGibbsの自由エネルギーが

とかける.ここで

定義から明らかなように各化学ポテンシャルは示強変数なので

となるので化学ポテンシャルは温度と圧力に加えて気体 1  2  の粒子数の比 N_2/N_1  によって決まる.
p  種類ある多成分系では,

となるので温度と圧力,さらに p-1  個の気体 1  との粒子数の比によって決まる.
平衡状態を記述する熱力学変数はGibbsの自由エネルギーを選んだ場合,温度と圧力と粒子数で p+2  個ある.
しかし実質的には 1  つ余分で p+1  個で決定できる.

熱力学的自由度という.
より一般論は相転移の章で議論する.

Problems

\textsc{Problem1.}

p 種類の理想気体が混合された系のHelmholtzの自由エネルギー F とGibbsの自由エネルギー G を求めよ.
また各気体の化学ポテンシャルを計算せよ.

\textsc{Solution.}

まずこの系の内部エネルギーとエントロピーは T,\,V,\,N_1,\cdots,N_p  の函数として,

ここで c_i,\,S_{i0}  は気体ごとに決まる定数.
よってLegendre変換 F=U-TS  によって,

ここで F_{i0}=c_i-S_{i0}/k_{\mathrm{B}}  とおいた.

さらにLegendre変換 G=F+PV  と状態方程式 PV=\sum_jN_jk_{\mathrm{B}}T  によって,

ここで G_{i0}=F_{i0} +1-\ln k_{\mathrm{B}}  とおいた.
分圧 P_i = N_iP/\sum_jN_j = x_iP  とおくと,

ともかける.

最後に化学ポテンシャルは, \sum_jN_j\ln P_j  N_i  微分が 0  であることに注意して

と計算される.

コメントを残す