ビリアル展開

ビリアル展開（ビリアルてんかい、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_{P}P+C_{P}P^{2}+...$

で表される。ここで、P は圧力、V_m は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 _{n=1}^{\infty }{\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_{1}dr_{2}$

$\beta _{2}={\frac {1}{2V}}\int f_{12}f_{23}f_{31}dr_{1}dr_{2}dr_{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_{1}dr_{2}dr_{3}dr_{4}$

のようになる。このように、$\beta _{1}$は2分子間、$\beta _{2}$は3分子間、$\beta _{3}$は4分子間の相互作用を表していることがわかる。
これより、第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$と表される。

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

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

$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}$は分子間の排除体積効果に関係している。

関連項目

• 状態方程式