ビリアル展開

出典: フリー百科事典『ウィキペディア(Wikipedia)』
移動: 案内検索

ビリアル展開(ビリアルてんかい、virial expansion)とは、実在気体の圧力(主に圧縮因子の形で)や浸透圧を、温度圧力に依存する様子を解析的に表すためにモル体積の逆数の冪級数に展開することである。ヘイケ・カメルリング・オネス1901年に提出した。ビリアル (virial) とはを意味するラテン語である。

実在気体の状態方程式について[編集]

理想気体の場合は圧縮因子が1であるが、実在気体はそうではないのでそのずれを補正する項として圧力(P)や体積の逆数(\frac{1}{V})の冪級数で表した(つまりビリアル展開した)のがビリアル方程式である。

Z = \frac{PV_m}{RT} = 1 + \frac{B_V}{V_m} + \frac{C_V}{V_m^2} + ...

または P の冪級数では

Z = \frac{PV_m}{RT} = 1 + B_PP + C_PP^2 + ...

で表される。ここで、P圧力Vm は1モルあたりの体積(モル体積)、R気体定数T温度である。B, C, ... は温度など分子間の相互作用に依存し、実験的に求められる各温度での気体ごとの定数で、ビリアル係数という。それぞれ第2ビリアル係数 (second virial coefficient)、第3ビリアル係数、……と呼ばれる。理想気体の場合、または実在気体でも圧力0の極限では Z は1になり、圧力が上がるごとに、高次の Pの寄与が大きくなる。 それぞれのビリアル係数は温度の関数である。 第2項\frac{B}{V_m}は2分子間相互作用に、第3項\frac{C}{V_m^2}は3分子間の相互作用に由来している。 ビリアル方程式は、マイアー夫妻のクラスター展開の理論(1940年)によると、

\frac{Pv}{kT} = 1 - \sum^\infty_{n = 1} \frac{n}{n+1}\beta_{n}\rhoと表せる。(\rho = \frac{1}{v}

ここで、v = \frac{V}{N_A} 、k = \frac{R}{N_A} なので、\frac{PV_m}{RT} = \frac{Pv}{kT} である。 \beta_{n}既約クラスター積分と呼ばれるもので以下のように定式化される。

\beta_{n} = \frac{1}{n!}\frac{1}{V}\int \Sigma^{(s)} \Pi f_{ij}dr_1...dr_n,  f_{ij} = exp\left( \frac{- U_{ij}}{kT} \right) - 1, (Uは分子間ポテンシャル。添字は分子の番号を表す。)
なお、この分子間ポテンシャルUは、実験値にあうようなものがいくつか提案されている。

上に挙げた既約クラスター積分は、具体的に計算すると

\beta_1 = \frac{1}{V} \int f_{12}dr_1dr_2
\beta_2 = \frac{1}{2V} \int f_{12}f_{23}f_{31}dr_1dr_2dr_3
\beta_3 = \frac{1}{3!V} \int [3f_{12}f_{23}f_{34}f_{41} + 6f_{12}f_{23}f_{34}f_{41}f_{13} + f_{12}f_{23}f_{34}f_{41}f_{13}f_{24}]dr_1dr_2dr_3dr_4

のようになる。このように、\beta_1は2分子間、\beta_2は3分子間、\beta_3は4分子間の相互作用を表していることがわかる[1]
これより、第2ビリアル係数は、

B_V = -\frac{N_A}{2}\beta_1 = - \frac{N_A}{2} \int f_{12} dr_{12} = 2 \pi N_A \int_0^{\infty} r^2 \left[1 - exp\left(- \frac{U_{12}}{kT}\right) \right]drと表される。

この式は、ファン・デル・ワールスの状態方程式に現れる物質係数abをミクロに導くときに重要となる。


気体がファンデルワールスの状態方程式に従うとするならば、圧縮因子 Z は以下のようになる[2]

Z = \frac{PV_m}{RT} = \frac{1}{1 - b / V_m} - \frac{a}{RTV_m}

またx < |1|のとき以下の二項式展開が成立する。

(1 - x)^{-1} = 1 + x + x^2  + x^3 + ...

上の二項式展開を用いて状態方程式の \frac{1}{1 - b / V_m} の項を級数に展開し、圧縮因子 Z を用いた式で表すと以下のようになる。

Z = \frac{PV_m}{RT} = 1 +  \left(b - \frac{a}{RT} \right) \frac{1}{V_m} + \frac{b^2}{V_m^2} + ...

この式を使うと、実験で求めた第2ビリアル係数の定数部分からbが、温度に反比例する部分からaが求められる。

浸透圧について[編集]

浸透圧\Piも、実在気体の状態方程式と同様にビリアル展開する事が出来る。

\frac{\Pi}{cRT} = \frac{1}{M} + B_2 c + B_3 c^2 + ...

cは質量濃度、Mは分子量である。
このとき、第2ビリアル定数B_2は分子間の排除体積効果に関係している。

参考文献[編集]

  1. ^ 物理学辞典 三訂版、培風館 ISBN 4-563-02094-X
  2. ^ Gordon M. Barrow著、大門寛、堂免一成訳  『バーロー物理化学』 東京化学同人、1999年

関連項目[編集]