化学平衡 計算ツール｜Kc/Kp/Ksp・Le Chatelierの原理・ICE表対応

平衡定数の定義と単位系

基本式(質量作用の法則)

aA + bB ⇌ cC + dD → Kc = [C]^c × [D]^d ÷ ([A]^a × [B]^b)

可逆反応が平衡状態に達したとき、生成物の濃度積を反応物の濃度積で割った値は温度一定なら定数になります。これが質量作用の法則(1864年グルベルグ・ワーゲ)で、平衡定数Kc(cはconcentration=濃度)の定義式です。指数a,b,c,dは化学反応式の量論係数です。

Kc(濃度平衡定数)とKp(圧平衡定数)

Kp = Kc × (RT)^Δn (Δn=気体生成物mol数-気体反応物mol数)

気体反応ではKcの代わりに分圧を使ったKpも定義されます。両者はΔn(気体分子数の変化)を使って相互変換可能。工業化学では圧力測定が容易なためKpが多用されます。

平衡と速度の関係

平衡定数は正反応速度定数と逆反応速度定数の比 K = kf/kr。速度論と熱力学の橋渡しをする関係式で、化学工学の反応器設計の基礎となります。触媒はkfとkrを同じ倍率で速めるため、平衡は変えず到達時間のみ短縮します。

K定数の種類

化学平衡には目的に応じた複数の平衡定数があります。

本ツールはKc(濃度平衡定数)の簡易計算です。Kp・Ksp等は各専用ツールを参照してください。

代表的な平衡定数値(25℃)

ICE表(平衡濃度計算の手順)

ICE表(Initial-Change-Equilibrium)は平衡濃度を系統的に計算する定番手法。

この表からK=[C]/([A][B])=x/(1-x)² を解いてxを求めます。二次方程式を解く形になるため、K>1(平衡が右)か K<1(平衡が左)かで近似式の使い分けができます。本ツールは1:1:1の反応 A+B⇌C を想定し、xを近似解で表示しています。

Le Chatelierの原理

「平衡状態にある系に外的変化を加えると、その変化を打ち消す方向に平衡が移動する」というのがLe Chatelier(ルシャトリエ)の原理(1884年)。工業合成の設計指針として今も現役です。

温度依存(Van't Hoff式)

ln(K₂/K₁) = -ΔH°/R × (1/T₂ - 1/T₁)

平衡定数の温度依存性はVan't Hoff式で表現されます。ΔH°は反応エンタルピー(J/mol)、Rは気体定数(8.314 J/(mol·K))、T₁,T₂は絶対温度(K)。発熱反応(ΔH<0)は低温で有利、吸熱反応(ΔH>0)は高温で有利という原理から、工業合成の温度最適化に直結します。

Haber-Bosch法(発熱)は熱力学的には低温有利ですが、速度論的に低温では反応が進まないため、実プラントでは400-500℃と200-300atmの妥協点で運転しています。

工業合成の事例

Haber-Bosch法(アンモニア合成)

N₂+3H₂⇌2NH₃ (ΔH=-92.4 kJ/mol)。世界の窒素肥料の80%を担う基幹プロセス。1913年にBASFがドイツで工業化、Fritz HaberとCarl Boschがノーベル賞受賞。年間約1.5億トン生産され、化学肥料として世界の食糧生産の40%を支えるとされます。

接触法(硫酸製造)

2SO₂+O₂⇌2SO₃ をV₂O₅触媒下400-500℃で進行させ、SO₃を吸収塔でH₂SO₄化。世界最大の化学製品(年産2.6億トン)で、肥料・化学工業・石油精製に不可欠。

エステル化(可逆反応)

RCOOH+R'OH⇌RCOOR'+H₂O。平衡定数Kが数十程度で、水を除去することでLe Chatelierの原理により反応を右に進めます。酢酸エチル・エステル系溶剤・可塑剤の大量生産に応用。

メタノール合成

CO+2H₂⇌CH₃OH (ΔH=-90 kJ/mol)。Cu/ZnO/Al₂O₃触媒下250℃・50-100atmで生産。バイオ燃料・化学原料として重要度が増しています。

計算例(ICE表と平衡濃度)

例1: A+B⇌C, Kc=10, [A]₀=[B]₀=1.0M

Kc = x / (1-x)² = 10 → 10(1-x)² = x → 10x² -21x +10 = 0 → x = (21±√(441-400))/20 = 0.729または1.371(不適)。平衡濃度 [A]=[B]=0.271M, [C]=0.729M。転化率72.9%。

例2: 2SO₂+O₂⇌2SO₃, Kp=100 atm⁻¹, 初期SO₂=2atm, O₂=1atm

ICE表で SO₂→2-2x, O₂→1-x, SO₃→2x。Kp = (2x)²/((2-2x)²(1-x)) = 100。数値解でx≈0.87、平衡ではSO₃=1.74 atm(87%転化)。硫酸製造の接触法の原理。

例3: 水の自己解離 H₂O⇌H⁺+OH⁻, Kw=1×10⁻¹⁴

純水では [H⁺]=[OH⁻]=√Kw = 1×10⁻⁷ M。よって純水のpH=7。温度上昇でKwは増加(60℃で約1×10⁻¹³)、中性pHは7より低くなります。

例4: AgClの溶解度(Ksp=1.8×10⁻¹⁰)

AgCl(s)⇌Ag⁺+Cl⁻, Ksp=[Ag⁺][Cl⁻]。純水中では [Ag⁺]=[Cl⁻]=s → s²=1.8×10⁻¹⁰ → s=1.34×10⁻⁵ M。0.1M NaCl中では [Cl⁻]≈0.1、[Ag⁺]=1.8×10⁻⁹ M(共通イオン効果で1万倍減)。

例5: 酢酸の解離度(Ka=1.75×10⁻⁵, C=0.1M)

CH₃COOH⇌CH₃COO⁻+H⁺, ICE表でC(1-α)=0.1(1-α), Cα=0.1α×2。Ka=Cα²/(1-α)≈Cα²(αが小さい時)。α=√(Ka/C)=√(1.75×10⁻⁴)=0.013(1.3%)。pH=-log(0.1×0.013)=2.88。

シーン別活用例

よくある間違い・注意点

• KcとKpを混同 — 気相反応ではKpが実用的、液相ではKc。両者は Kp=Kc(RT)^Δn の関係。Δn=0の反応(2HI⇌H₂+I₂等)ではKc=Kpだが、Δn≠0で数桁ズレます。
• 純液体・純固体を平衡式に含める — 純物質の活量は1と定義され、平衡式に登場しません。CaCO₃(s)⇌CaO(s)+CO₂(g)ではKp=P(CO₂)のみで、固体を分母に入れるのは誤りです。
• 温度依存を無視 — Kは温度で大きく変わります。25℃のKと400℃のKを混同すると数桁誤差。Van't Hoff式で温度補正が必須です。
• 触媒で平衡が動くと誤解 — 触媒は正反応・逆反応を同じ倍率で速めるため、平衡定数は変わりません。反応時間を短縮するだけです。
• 不活性ガス添加で平衡移動と誤解 — 定容条件で不活性ガス(Ar等)を加えても各成分の分圧は変わらず、平衡は動きません。定圧条件では別の理屈です。
• ICE表の変化量の符号ミス — 反応物は減少(-x)、生成物は増加(+x)、量論係数を掛けます。2A→Bなら反応物は-2x、生成物は+x。
• Kが1に近い場合の近似式使用 — K>>1(1000以上)なら生成物側完全反応、K<<1(0.001以下)なら反応物側停止と近似できるが、中間領域(0.01-100)では二次方程式を厳密に解く必要があります。

計算ノウハウ(実務上のコツ)

1. Kの大きさで近似戦略を変える

K>>10³なら完全反応、K<<10⁻³なら反応なしと近似。中間領域(10⁻³〜10³)は二次方程式で厳密解を求めます。目分量で近似できるかを最初に判定するとエレガントに解けます。

2. ICE表を必ず書く

複雑な多成分系でも、ICE表を丁寧に書けば量論係数の誤り・符号ミスを防げます。反応物は-x、生成物は+xで量論係数を必ず掛けます(例: 2A→Bなら-2xと+x)。

3. 単位に注意

Kcは(mol/L)^Δn、Kpは(atm)^Δn。Δn=気体生成物mol-気体反応物molで符号が変わります。Δn=0の場合Kc=Kp(単位なし)、それ以外は換算式で相互変換が必要です。

4. 温度依存の見積もり

Van't Hoff式 ln(K₂/K₁) = -ΔH°/R×(1/T₂-1/T₁) で温度依存を計算。ΔH°>0(吸熱)は高温有利、ΔH°<0(発熱)は低温有利。反応エンタルピー ΔH°は化学便覧・NIST WebBookに表がありますので実務で参照。

5. 共通イオン効果の活用

難溶性塩の沈殿促進には共通イオン(例: Cl⁻でAgCl沈殿促進)を使用。0.1M NaCl添加でAgClの溶解度は1万分の1に低下します。定量分析・沈殿分離・分別再結晶で頻用される手法です。

6. 活量係数の考慮

イオン強度I > 0.01 Mでは活量係数γが有意にズレるため、厳密計算では活量a=γ·[濃度]を使います。Debye-Hückel式で近似計算、電位差滴定や電気化学では必須の補正です。

関連する規格・法令

• JIS K 0050 ― 化学分析方法通則。平衡計算・pH測定の基本規定。
• JIS K 0102 ― 工場排水試験方法。溶解度・平衡計算に基づく分析法。
• IUPAC Recommendations ― 平衡定数・活量・熱力学量の国際定義。
• NIST Standard Reference Data ― 米国NISTが管理する化学熱力学データベース(平衡定数含む)。
• 化学工学プラント設計基準(JACT) ― 化学工学会が発行する反応器設計指針。
• 環境省 化審法・化管法 ― 化学物質管理法。水/大気/土壌中の平衡挙動に基づく分配係数計算。

化学平衡の科学史

質量作用の法則(1864年)

ノルウェーのカト・グルベルグ(1836-1902)とペーター・ワーゲ(1833-1900)が「反応速度は反応物濃度の積に比例する」ことを実験的に発見、後の平衡定数概念の礎に。当時は化学界で見過ごされましたが、後年再評価され、化学熱力学の基礎法則として確立。

ファントホッフの化学熱力学(1884年)

オランダのヤコブス・ファントホッフが「Etudes de Dynamique Chimique」で化学平衡の熱力学的定式化を提示、Van't Hoff式ln(K₂/K₁)=-ΔH°/R×(1/T₂-1/T₁)を導出。1901年第1回ノーベル化学賞受賞。

Le Chatelierの原理(1884年)

フランスのアンリ・ル・シャトリエ(1850-1936)が「平衡系に外部変化を加えるとそれを打ち消す方向に平衡が移動する」原理を発表。工業合成の設計指針として今も現役です。

Haber-Bosch法(1913年)

ドイツのフリッツ・ハーバー(1868-1934)とカール・ボッシュ(1874-1940)がアンモニア工業合成を確立、Le Chatelierの原理の実用化例として最大の成果。ハーバーは1918年、ボッシュは1931年ノーベル化学賞受賞。世界の窒素肥料生産の80%を担い、20世紀の食糧増産を支えました。

ギブスの化学熱力学(1876年)

米国のジョサイア・ウィラード・ギブス(1839-1903)が化学平衡を自由エネルギー ΔG=ΔH-TΔS で定式化、化学ポテンシャル・相律・化学平衡定数の熱力学的な統一理論を構築。現代の化学工学プロセスシミュレーションの理論基盤です。

現代の応用

Aspen Plus・HYSYS・COMSOL等の化学工学プロセスシミュレーターは、内部で平衡定数計算・熱力学モデル(UNIQUAC・NRTL・SRK等)を実装。石油精製・アンモニア合成・エチレン製造など巨大化学プラントの設計・最適化に日常的に使用されています。

参考文献・公的資料

より正確な平衡定数データや詳細な理論解説は、以下の公的機関・学会の資料を参照してください。

• NIST Chemistry WebBook ― 米国国立標準技術研究所の化学物性データベース。平衡定数・熱力学量を無料検索可能。
• 産業技術総合研究所 計量標準総合センター(NMIJ) ― 化学熱力学量の国家標準。認証標準物質(CRM)。
• IUPAC(国際純正・応用化学連合) ― 平衡定数・活量の国際定義。SC-Database(酸解離定数)を公開。
• 日本化学会 ― 化学便覧・化学工学便覧の平衡データ標準的引用元。
• 化学工学会(SCEJ) ― 化学工学プロセス設計・反応工学の学会。
• 日本分析化学会 ― 分析化学における平衡計算の実務。
• 日本産業標準調査会(JISC) ― JIS K 0050・K 0102等の公式索引。
• 文部科学省 高等学校学習指導要領(化学) ― 化学平衡・Le Chatelier原理の教育カリキュラム。
• 環境省 ― 環境化学・化学物質管理における平衡計算。
• 国立環境研究所(NIES) ― 大気・海洋・土壌の化学平衡研究。

よくある質問（FAQ）

関連計算ツール