ファンデルワールスの状態方程式

ファンデルワールスの状態方程式とは実在気体の、圧力と体積の関係を近似するためにファン・デル・ワールスが提案した式。

概要

理想気体の状態方程式によれば気体の体積はそれが受ける圧力とそれ自体の温度(とその気体分子のモル数)にのみ依存することになる。よって、理想気体の法則は次のような性質を持つ気体に対してのみ適用できる。

1. 気体分子は体積を有しない
2. 分子間には相互作用はない

1.の問題を体積排除効果、2.の問題を分子間の引力効果と呼ぶことにする。明らかに、体積排除効果も分子間の引力効果も満たす気体は存在しない。だが、理想気体の状態方程式は高温あるいは10atm以下の低圧ではかなり有効である。気体分子間の間隔が大きいほど、相互作用も気体分子そのものの体積も小さくなるからだ。実在気体では圧力の増加とともに理想気体の法則からはずれ、圧力と体積の積は一定ではなくなる。その傾向は気体の種類によっても異なり、同一気体については低温・高圧であるほどそのずれが大きくなる。実在気体の理想気体からのずれを圧縮率因子(Zで表す場合が多い)で表すことがある。理想気体の状態方程式のずれを、圧縮因子をプロットすることで測定できる。

$z = {PV \over nRT}$

ファン・デル・ワールスは以上の点を考慮して、1873年に下記のような実在気体の状態方程式を提出した。

$\left( P+\frac{a}{V^2}\right) \left( V-b\right) = RT$

これを1モルの気体に対するファン・デル・ワールスの状態式という。 ここに、a と b はファンデルワールス定数といい、気体の種類によって定まる定数である（a は上記の2に、b は上記の1に由来する定数である）。 　

ファン・デル・ワールスの状態方程式はジュール・トムソン効果や気液相転移について正しい振る舞いを再現できる上、解析的扱いが易しいため頻繁に用いられるが、（理論的背景はあっても）あくまで1つのモデル（近似式）であり厳密に気体の振る舞いを再現できている訳ではない。また、ビリアル展開のように系統的に近似の精度を上げていく事はできないのが欠点ではある。

熱力学的性質

気体がファン･デル･ワールスの状態方程式に従う場合、$\left( \frac{\partial U}{\partial V} \right)_{T} = \frac{a}{V^2}$なので理想気体と異なり内部エネルギーは体積にも依存する。定式化すると以下のようになる。

内部エネルギー：$U = C_V T - \frac{a}{V} + U_0$

（C_Vは定積モル比熱、U₀は定数）

プロット

分子間の引力効果について、気体の1分子が持つ相互作用の有効範囲である体積を V₀ , V₀ の物質量を N₀ とすると, N₀ 個の分子から、2つの分子間の相互作用の組み合わせは、

${}_{N_0}{\rm C}_{2} = {N_0! \over (N_0 - 2)!2!}$

${} = {N_0(N_0 - 1) \over 2}$

${} = {N_0^2 \over 2}$

${} = {1 \over 2}n^2{V_0^2 \over V^2}$

${} = {V_0^2 \over 2} \times \left({n \over V}\right)^2$

個々の分子が容器に及ぼす圧力は、壁と分子の衝突の頻度および分子によって壁に伝えられる運動量に依存する。どちらも分子間力によって減少する。この式から、圧力の減少分は、V₀ と密度 ${n \over V}$ に依存することが分かる。

${V_0^2 \over 2} = a$

と定義すると、a は分子の種類によって定まる比例定数である。a は b と共にファンデルワールス定数と呼ばれる。

臨界定数との関係

ファンデルワールスの式による等温線

ファンデルワールス定数 a および b は物質の臨界点、すなわち臨界温度 T_c 、臨界圧力 P_c および臨界体積 V_c と以下の関係にある。これはファンデルワールスの状態方程式を V で偏微分し、等温線の極大値および極小値が互いに接近して消失する点を求めることにより得られる。

$T_c = \frac{8 a}{27 bR}$

$P_c = \frac{a}{27 b^2}$

V_c = 3b

ただし、ファンデルワールスの状態方程式は近似式であるから、実測によるファンデルワールス定数は臨界定数から計算した値とは厳密には一致しない。

主な気体のファンデルワールス定数と臨界定数は以下の通りである。

ファンデルワールス定数および臨界定数

	a / Pa m6 mol−2	b / m3 mol−1	Tc / K	Pc / Pa	Vc / m3 mol−1
空気	1.35 × 10−3	36.6 × 10−6	132.5	3.766 × 106	88.1 × 10−6
ヘリウム He	3.45 × 10−3	23.8 × 10−6	5.201	0.227 × 106	57.5 × 10−6
水素 H2	24.8 × 10−3	26.7 × 10−6	33.2	1.316 × 106	63.8 × 10−6
窒素 N2	141 × 10−3	39.2 × 10−6	126.20	3.400 × 106	89.2 × 10−6
酸素 O2	138 × 10−3	31.9 × 10−6	154.58	5.043 × 106	73.4 × 10−6
二酸化炭素 CO2	365 × 10−3	42.8 × 10−6	304.21	7.383 × 106	94.4 × 10−6
水蒸気 H2O	553 × 10−3	33.0 × 10−6	647.30	22.12 × 106	57.1 × 10−6

ビリアル方程式

ファンデルワールスの状態方程式をビリアル展開するとビリアル方程式（ビリアルほうていしき、virial equation）が得られる。ビリアル方程式はビリアルの式、ビリヤルの式、ビリアル状態式あるいはオネスの状態方程式、オンネスの状態方程式とも呼ばれる。

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

圧縮率因子 Z で表すと以下のような一般式になる。ここで B_V、C_V、D_V …はビリアル係数と呼ばれる。V_m はモル体積である。

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

ファンデルワールスの状態方程式の第2ビリアル係数 B_V は b − a / RT 、第3ビリアル係数 C_V は b²、第4ビリアル係数 D_V は b³となる。 第2ビリアル係数を実験的に求めれば、温度に依存する部分と定数部分とからaとbの定数を決定する事ができる。

関連項目

• ディーテリチの状態方程式