ギブズの相律 シミュレーター 戻る
相平衡シミュレーター

ギブズの相律 シミュレーター — 相平衡の自由度

ギブズの相律 F = C − P + 2 − N により、成分数 C・相数 P・追加拘束数 N から相平衡の自由度 F、1 つの強度量固定下の F、2 つの強度量固定下の F、共存可能な最大相数 P_max を実時間に計算します。T-P 相図の概念図と F-P 線図で不変点・単変点・二変点を可視化します。

パラメータ設定
成分数 C
相数 P
追加拘束数 N
P_max 判定用拘束数

既定値は C=2、P=2、N=0、P_max 判定用拘束数=0。F = C − P + 2 − N で計算し、F < 0 は物理的に不可能(指定相数 P が大きすぎる)と判定します。1 つの強度量固定下では F = C − P + 1、2 つの強度量固定下では F = C − P で計算します。

計算結果
自由度 F
1 強度量固定下の F
2 強度量固定下の F
最大相数 P_max
T-P 相図の概念図(水)

横軸=温度 T/縦軸=圧力 P/青=固相(F=2)/緑=液相(F=2)/橙=気相(F=2)/白線=相境界(F=1)/黄点=三重点(F=0)/赤点=臨界点(相律外)

F-P 線図(C と N 固定)

横軸=相数 P/縦軸=自由度 F/青=F = C − P + 2 − N の直線/赤領域=F < 0(物理的に不可能)/黄マーカー=現在の (P, F)

理論・主要公式

ギブズの相律は、相平衡における示強変数の自由度を成分数と相数から決定します:

$$F = C - P + 2 - N$$

$F$ は自由度(独立に変化できる強度量の数)、$C$ は独立成分数、$P$ は共存する相の数、$+2$ は温度 $T$ と圧力 $P$ の 2 つの強度量、$N$ は追加で固定する拘束数です。$F=0$ は不変点、$F=1$ は単変点(線)、$F=2$ は二変点(領域)です。共存可能な最大相数は次式で与えられます:

$$P_{\max} = C + 2 - N$$

$N=0$ かつ純物質($C=1$)では $P_{\max}=3$ で、これが三重点を 1 点だけ持つことの理由です。一定圧力下($N=1$)の二成分系($C=2$)では:

$$F = C - P + 1$$

これは Y-X 相図や T-X 相図で液相線・固相線の上では F=1(共存線)、領域内では F=2 となる二相共存条件の根拠です。

ギブズの相律 シミュレーターとは

🙋
「ギブズの相律」って F = C − P + 2 で覚えたんですけど、この「自由度 F」って具体的に何を意味してるんですか?気温と気圧を測れるってことですか?
🎓
いい質問だ。F は「相平衡を崩さずに独立に変化できる示強変数(温度・圧力・組成など)の数」を意味する。例えば C=1(水だけ)で P=1(液相のみ)なら F = 1 − 1 + 2 = 2 で、T と P を独立に動かしても水のままだ。同じ C=1 で P=2(液と気が共存)なら F = 1 で、T を決めれば飽和蒸気圧が一意に決まる。本ツールの既定値 C=2、P=2、N=0 では F=2、つまり二変点領域。
🙋
じゃあ「三重点」って F=0 の不変点ってことですか?水の三重点が 0.01°C、611 Pa って習いましたけど、なんでそんなにピッタリ決まるんですか?
🎓
そのとおり、三重点は F=0 の不変点だ。C=1、P=3(固液気共存)を入れると F = 1 − 3 + 2 = 0、つまり自由度ゼロで T と P が一意に決まる。本ツールで C=1、P=3 にすると F=0 が出る。SI 単位系では水の三重点 273.16 K = 0.01°C が温度の定義(旧定義)に使われていたほど精密に再現できる点で、これは相律が保証する自然法則上の必然だ。
🙋
「拘束数 N」のスライダーって何ですか?「一定圧力下では F が 1 減る」みたいな話と関係ありますか?
🎓
まさにそれだ。N=1 にすると圧力を 1 atm に固定したことになり、F = C − P + 2 − 1 = C − P + 1 になる。例えば Fe-C 合金(C=2)で常圧下の二相共存(P=2)なら F = 2 − 2 + 1 = 1、つまり温度を決めれば組成が一意に決まる。これが金属の状態図(T-X 図)の理論的根拠だ。N=2 にすれば T と P 両方固定で F = C − P。本ツールでスライダーを動かすと、それぞれの統計カードに値が出るようにしてある。
🙋
最後に、F < 0 になったら「不可能」って表示されてますけど、これって相数が多すぎるってことですか?
🎓
そう、これも相律の重要な帰結だ。共存可能な最大相数は P_max = C + 2 − N で、純物質(C=1)かつ N=0 なら最大 3 相(三重点まで)、二成分系(C=2)なら最大 4 相、三成分系で最大 5 相となる。それを超える相を共存させようとすると F < 0 となり「物理的に不可能」になる。Gibbs(1875 年)の業績はこの単純な不等式で「世界中の相図がなぜそういう形をしているか」を一網打尽に説明した点にある。本ツールで C と P を極端な組み合わせにすると赤領域(F<0)が見える。

よくある質問

ギブズの相律は、定温・定圧条件下の相平衡で「独立に変化できる強度量(示強変数)の数」すなわち自由度 F を、成分数 C と相数 P から決定する基本則です。式は F = C − P + 2 − N で、+2 は温度 T と圧力 P の 2 つの強度量に対応し、N は追加で固定する拘束数(一定圧力下の実験なら 1 など)です。本ツールの既定値(C=2、P=2、N=0)では F = 2 − 2 + 2 − 0 = 2 となり、T と P を独立に動かせる二変点領域です。
F = 0 は不変点と呼ばれ、温度も圧力も値が一意に決まる「点」で表される平衡状態です。水の三重点(0.01°C、611.657 Pa)が代表例で、ここでは固相・液相・気相が同時に共存します。F = 1 は単変点で、温度を決めれば圧力が一意に決まる「線」上の平衡(蒸気圧曲線、融解曲線、昇華曲線など)です。F = 2 は二変点で、温度と圧力を独立に動かせる「領域」で、純粋な相が一つだけ存在する状態です。
拘束数 N は追加で固定する強度量の数を表します。例えば実験で圧力を 1 atm に固定すると N = 1、温度と圧力を両方固定すると N = 2 となり、自由度はそれぞれ 1、2 だけ減少します。1 atm 固定下の二成分系では F = C − P + 1 となり相律最大相数は C + 1。最大相数 P_max は F = 0 となる条件から P_max = C + 2 − N で、C = 2 かつ拘束なしなら 4 相、C = 1 の純物質では P_max = 3 相(三重点まで)と分かります。
水のような一成分系(C=1)では相律は F = 3 − P となり、P=1 は二変点(領域)、P=2 は単変点(線)、P=3 は不変点(点)です。三重点(0.01°C、611.657 Pa)は固液気の 3 相が同時に共存する F=0 の点で、相図上で必ず一意に決まります。臨界点(373.946°C、22.064 MPa)は液相と気相の区別が消失する特異点で、相律の枠組み外(連続的に超臨界流体に移行)として扱われます。本ツールで C=1、P=3 にすると F=0(三重点条件)が確認できます。

実世界での応用

金属合金の状態図設計(Fe-C 系、Al-Cu 系など):常圧下(N=1)の二成分系(C=2)では相律は F = 3 − P となります。Fe-C 合金鋳鉄・鋼の状態図ではオーステナイト+セメンタイト+液相の三相共存点(共晶 1148°C、4.3wt%C、共析 727°C、0.76wt%C)が F=0 の不変点として現れ、ここを境に組織が劇的に変化します。本ツールに C=2、N=1、P=3 を入れると F=0 と表示され、なぜ共晶・共析点が一点に決まるか理解できます。熱処理プロセス設計(焼入れ、焼戻し、時効硬化)の理論的基盤です。

蒸留・抽出プラントの設計(Aspen Plus、ASPEN HYSYS):化学プロセス工学では三成分系(C=3、例:メタノール-エタノール-水)の常圧蒸留塔設計に相律が直接効きます。F = 3 − P + 1 = 4 − P で、気液二相(P=2)なら F=2、すなわち温度と一成分の組成が独立、残り二成分は従属。共沸点(azeotrope)はこの自由度の縮退点で、相律の制約から「越えられない蒸留壁(distillation boundary)」が生じます。本ツールで C=3 にして相図の理解を深め、抽出蒸留や圧力変動蒸留が必要な系を見極められます。

地球科学・マグマの結晶分化:地殻のケイ酸塩鉱物(SiO₂、Al₂O₃、CaO、Na₂O、MgO 等)の多成分系の相平衡に相律が応用されます。N. L. Bowen の反応系列で有名な玄武岩マグマの分化過程は、温度低下に伴ってかんらん石→輝石→角閃石→雲母の順に晶出することを Mg-Fe-Si-O 四成分系の相律(F = 6 − P)で説明します。地球深部の超高圧下の鉱物相転移も同じ枠組みで議論され、本ツールで C=4、P=2〜5 にスイープすると相分離の自由度がどう推移するか観察できます。

食品科学・凍結乾燥(freeze drying):食品工業の凍結乾燥工程では水の三重点(0.01°C、611 Pa)以下に圧力を下げて氷を直接気化(昇華)させます。本ツールで C=1、P=2(固+気)にすると F=1 となり、温度を決めれば昇華圧が一意に決まることが分かります。コーヒー・医薬品・宇宙食の品質はこの単純な相律で支配されており、運転点を三重点より十分低圧側に保つことが「形を壊さず水を抜く」鍵となります。

よくある誤解と注意点

最も多い誤解が、「成分数 C はそこに存在する物質の種類の数」というものです。C は厳密には「化学反応式で互いに独立に決められる独立成分の数」で、化学反応や電気的中性条件などの拘束を引き算した値です。例えば水を H₂、O₂、H₂O の 3 種で表しても、化学平衡 2H₂ + O₂ ⇌ 2H₂O の拘束 1 つで C = 3 − 1 = 2 となります(さらに化学量論初期条件があれば C = 1)。本ツールの C は「独立成分数」として扱う前提でモデル化していますので、複雑な反応系では化学量論を整理してから入力してください。

次に多いのが、「相律は熱力学第二法則と独立」という混同です。実は相律は化学ポテンシャル μ_i の連続性条件(相間で μ_i が等しい)と独立変数の数え上げから直接導かれる純粋に「変数の数の問題」であり、第二法則の本質的内容ではありません。Gibbs(1875 年)の元論文 “On the Equilibrium of Heterogeneous Substances” の中で、Gibbs エネルギー最小化の必要条件として提示されました。本ツールの計算ロジックも単純な四則演算のみで、熱力学パラメータは一切使っていません。

最後に、「相律は超臨界流体や液晶にも適用できる」という拡大解釈の限界です。臨界点では液相と気相の区別が連続的に消失するため、Gibbs の相律は厳密には適用できません(一次相転移を前提とした枠組みのため)。同様に液晶のネマチック相・スメクチック相間の二次相転移、超伝導転移、合金の規則-不規則転移なども、相律の素朴な適用ではなく Landau の対称性破れの理論や臨界指数の議論が必要です。本ツールは一次相転移に限定した教育用ツールであり、超臨界流体や臨界現象の解析には別の枠組みを参照してください。