技術(非IT系)

ギブズ-ヘルムホルツの式とは?導出方法も解説!(熱力学・状態関数・温度依存性・エントロピー・ギブズ自由エネルギーなど)

当サイトでは記事内に広告を含みます

ギブズ-ヘルムホルツの式は、熱力学における重要な関係式のひとつであり、ギブズ自由エネルギーの温度依存性を系統的に記述する式です。

化学反応の温度依存性・電気化学的起電力の温度変化・平衡定数の温度依存性(ファント・ホッフの式)など、幅広い熱力学的問題の解析に活用されます。

エンタルピー変化量ΔHを実験的に求めることでギブズ自由エネルギーの温度依存性を予測できるという実用的な側面も持ち合わせています。

本記事では、ギブズ-ヘルムホルツの式の定義から始まり、ヘルムホルツ自由エネルギーを用いた導出方法、化学熱力学への応用、電気化学・平衡定数との関係まで丁寧に解説していきます。

熱力学の理論的背景をしっかり把握したい方にとって、体系的な理解を深めるための内容となっているでしょう。

ギブズ-ヘルムホルツの式とは何か?定義と物理的意味

それではまず、ギブズ-ヘルムホルツの式の定義と物理的意味について解説していきます。

ギブズ-ヘルムホルツの式には複数の表現形があり、ギブズ自由エネルギーGとヘルムホルツ自由エネルギーFのそれぞれについて定式化されています。

【ギブズ-ヘルムホルツの式(ギブズ自由エネルギー版)】

(∂(G/T)/∂T)p = −H / T²

または等価な表現:

(∂(ΔG/T)/∂T)p = −ΔH / T²

G:ギブズ自由エネルギー(J)

H:エンタルピー(J)

T:絶対温度(K)

【ヘルムホルツ自由エネルギー版】

(∂(F/T)/∂T)V = −U / T²

F:ヘルムホルツ自由エネルギー(J)

U:内部エネルギー(J)

この式の物理的意味は、ギブズ自由エネルギーを温度Tで割った量(G/T)の温度微分が、エンタルピーHと温度二乗の比によって決まるということです。

実験的にエンタルピー変化ΔHを測定することで、異なる温度でのΔGを予測できるという実用的価値が非常に高い式といえるでしょう。

ギブズ自由エネルギーの全微分からの導出

ギブズ-ヘルムホルツの式は、ギブズ自由エネルギーの全微分式から導くことができます。

【ギブズ-ヘルムホルツの式の導出】

ギブズ自由エネルギーの定義:G = H − TS

→ エントロピー:S = −(∂G/∂T)p

G = H + T×(∂G/∂T)p ……(①)

G/T の温度微分を計算:

∂(G/T)/∂T = (1/T)×(∂G/∂T)p − G/T²

(∂G/∂T)p = −S を代入:

∂(G/T)/∂T = −S/T − G/T²

G = H − TS → −TS = G − H を利用:

∂(G/T)/∂T = (G − H)/T² − G/T² = −H/T²

∴ (∂(G/T)/∂T)p = −H/T²

この導出過程から、ギブズ-ヘルムホルツの式はエントロピーの温度微分関係とギブズ自由エネルギーの定義式G=H−TSから自然に導かれる恒等式であることがわかります。

追加の仮定なしに成立する厳密な熱力学関係式であるため、幅広い系に適用可能でしょう。

1/Tを変数とした積分形の活用

ギブズ-ヘルムホルツの式は1/Tを独立変数として書き直すことで、実用計算に使いやすい形になります。

1/T = x とおくと:

(∂(G/T)/∂(1/T))p = H

積分形:

ΔG(T₂)/T₂ = ΔG(T₁)/T₁ − ∫(T₁→T₂)(ΔH/T²)dT

ΔHが温度に依存しない(一定)と近似できる場合:

ΔG(T₂)/T₂ = ΔG(T₁)/T₁ + ΔH × (1/T₂ − 1/T₁)

ΔHが温度によらず一定という近似は狭い温度範囲では妥当であり、この積分形を使うことで基準温度T₁でのΔGから任意温度T₂でのΔGを計算できます。

標準生成エンタルピーΔH°と25℃でのΔG°が与えられれば、他の温度でのΔG°を比較的容易に推算できる点が化学熱力学の実務で重宝されているでしょう。

ヘルムホルツ自由エネルギー版の導出と応用

続いては、ヘルムホルツ自由エネルギー版のギブズ-ヘルムホルツの式の導出と応用について確認していきます。

等温・等体積系の解析に適したヘルムホルツ版は、統計力学・物性物理・高圧化学などの分野で特に有用です。

ヘルムホルツ版の導出

ヘルムホルツ自由エネルギーF = U − TSを使ってギブズ-ヘルムホルツの式を導出します。

F = U − TS

エントロピー:S = −(∂F/∂T)V

F/T の温度微分:

(∂(F/T)/∂T)V = (1/T)×(∂F/∂T)V − F/T²

= −S/T − F/T²

F = U − TS → −TS = F − U を利用:

= (F − U)/T² − F/T²

= −U/T²

∴ (∂(F/T)/∂T)V = −U/T²

ヘルムホルツ版では内部エネルギーUがギブズ版のエンタルピーHに対応しており、等体積条件(∂F系)での温度依存性を内部エネルギーUから計算できます。

統計力学の分配関数ZからF = −kBT ln Zを経由して内部エネルギーを計算し、さらにヘルムホルツ版ギブズ-ヘルムホルツ式で温度依存性を解析するという流れが理論物理での標準的な計算手順でしょう。

統計力学と分配関数への応用

統計力学では、分配関数Zを使って全ての熱力学量を計算します。

F = −kBT ln Z より:

F/T = −kB ln Z

(∂(F/T)/∂T)V = −kB(∂ ln Z/∂T)V = −U/T²

→ U = kBT² (∂ ln Z/∂T)V

これは内部エネルギーの統計力学的表現式です。

ヘルムホルツ版ギブズ-ヘルムホルツの式は、分配関数Zの温度微分から内部エネルギーUを導出する基礎となる恒等式として統計力学の理論体系を支えています。

固体の格子振動・電子気体・磁性体・量子気体など多くの物性計算でこの関係式が活用されているでしょう。

内部エネルギーとエントロピーの分離

ヘルムホルツ版ギブズ-ヘルムホルツの式を用いると、ヘルムホルツ自由エネルギーFをエネルギー的寄与(U)とエントロピー的寄与(TS)に分離して評価できます。

F = U − TS の分離評価:

U = F − T × (∂F/∂T)V = F + T² × (∂(F/T)/∂T)V

S = −(∂F/∂T)V

TS = −T × (∂F/∂T)V

エネルギー的寄与とエントロピー的寄与を定量的に分離することで、ゴム弾性・高分子溶液・タンパク質折りたたみなどのエントロピー駆動現象のメカニズムを詳細に解析できます。

これはヘルムホルツ版ギブズ-ヘルムホルツの式が材料科学・生物物理の研究で特に有用な理由のひとつでしょう。

化学熱力学でのギブズ-ヘルムホルツの式の応用

続いては、化学熱力学でのギブズ-ヘルムホルツの式の応用について確認していきます。

化学反応の温度依存性の解析・工業プロセスの最適温度の設計・燃料電池の起電力予測など、化学工学・反応工学の実務においてこの式は基盤的な役割を担っています。

化学反応のΔGの温度依存性計算

化学反応のギブズ自由エネルギー変化ΔG°(標準状態)の温度依存性を計算するうえでギブズ-ヘルムホルツの式は直接役立ちます。

【化学反応への適用例】

反応:N₂ + 3H₂ → 2NH₃(アンモニア合成)

298 K での値:ΔG°₂₉₈ = −32.9 kJ/mol、ΔH°₂₉₈ = −92.4 kJ/mol

500 K でのΔG°₅₀₀を計算(ΔHを一定と近似):

ΔG°₅₀₀/500 ≈ ΔG°₂₉₈/298 + ΔH°(1/500 − 1/298)

ΔG°₅₀₀/500 ≈ (−32900/298)+ (−92400)×(0.002 − 0.00336)

ΔG°₅₀₀/500 ≈ −110.4 + (−92400)×(−0.00136)

ΔG°₅₀₀/500 ≈ −110.4 + 125.7 ≈ 15.3 J/(mol·K)

ΔG°₅₀₀ ≈ 15.3 × 500 = 7650 J/mol ≈ +7.7 kJ/mol

この計算からアンモニア合成反応は298 Kでは自発的(ΔG°<0)ですが、温度を上げると自発性が失われる(ΔG°>0)傾向があることがわかります。

実際のハーバー・ボッシュ法では高温が必要な反応速度と低温が有利な熱力学平衡のトレードオフを考慮した最適温度(400〜500℃)が選択されており、これまさにギブズ-ヘルムホルツの式による解析の典型的応用でしょう。

ファント・ホッフの式との関係

ギブズ-ヘルムホルツの式と平衡定数Kpの温度依存性を記述するファント・ホッフの式は密接に関連しています。

ΔG° = −RT ln Kp の関係から:

ΔG°/T = −R ln Kp

ギブズ-ヘルムホルツの式:(∂(ΔG°/T)/∂T)p = −ΔH°/T²

→ −R × (∂ ln Kp/∂T)p = −ΔH°/T²

∴ (∂ ln Kp/∂T)p = ΔH°/(RT²)

これがファント・ホッフの式です。

ファント・ホッフの式はギブズ-ヘルムホルツの式から直接導出されることがわかり、両者が同じ熱力学的基盤に立脚しています。

ln Kp対1/Tのプロット(ファント・ホッフプロット)の傾きから反応のΔH°を実験的に求める手法は、化学・生化学・医薬品研究で標準的に用いられているでしょう。

電気化学での起電力の温度依存性

電池・燃料電池の起電力Eの温度依存性もギブズ-ヘルムホルツの式から導かれます。

電気化学の関係:ΔG = −nFE

(n:電子数、F:ファラデー定数)

ギブズ-ヘルムホルツの式に代入:

(∂(−nFE/T)/∂T)p = −ΔH/T²

起電力の温度係数:(∂E/∂T)p = ΔS/(nF)

ΔH = −nF〔E − T(∂E/∂T)p〕

起電力の温度依存性(∂E/∂T)を測定することで電池反応のΔSとΔHを同時に決定できるという実用的な電気熱量測定の基礎がここに含まれています。

燃料電池の設計では動作温度における起電力の予測と最適化にギブズ-ヘルムホルツの式の解析が活用されているでしょう。

ギブズ-ヘルムホルツの式を用いた実践的計算

続いては、ギブズ-ヘルムホルツの式を用いた実践的な計算手順について確認していきます。

実験データからΔHを求める方法・グラフ解析・コンピュータを使った数値計算など、実務で使える手法を紹介します。

ファント・ホッフプロットによるΔHの実験的決定

複数の温度での平衡定数Kpを測定してファント・ホッフプロット(ln Kp対1/T)を作成することで、反応エンタルピーΔH°を実験的に求めることができます。

【ファント・ホッフプロットの手順】

①複数の温度T₁、T₂、T₃…で平衡定数Kp₁、Kp₂、Kp₃…を測定

②縦軸にln Kp、横軸に1/Tをプロット

③直線の傾きm = −ΔH°/R から:

 ΔH° = −R × m(J/mol)

④y切片 b = ΔS°/R から:

 ΔS° = R × b(J/(mol·K))

⑤298 Kでの標準ギブズエネルギー:

 ΔG°₂₉₈ = ΔH° − 298 × ΔS°(J/mol)

ファント・ホッフプロットの直線性は反応のΔHが温度によらず一定(近似)であることを示しており、直線からのずれはΔHの温度依存性(熱容量差ΔCp)を示唆します。

タンパク質の熱変性・配位子結合・核酸二重鎖の融解など生化学系の研究でも広く用いられている信頼性の高い解析手法でしょう。

熱力学データベースを用いた計算例

NISTなどの熱力学データベースから標準生成エンタルピー・標準ギブズエネルギーを取得してギブズ-ヘルムホルツの式で温度依存性を計算する手順を確認しましょう。

反応 ΔH°₂₉₈(kJ/mol) ΔG°₂₉₈(kJ/mol) ΔS°₂₉₈(J/(mol·K))
H₂ + ½O₂ → H₂O(g) −241.8 −228.6 −44.3
C + O₂ → CO₂ −393.5 −394.4 +2.9
N₂ + 3H₂ → 2NH₃ −92.4 −32.9 −198.5

ΔS°の符号がマイナスの反応(生成物のエントロピーが小さい)は高温になるほどΔG°が増大して自発性が失われ、プラスの反応は高温で自発性が増すという傾向がギブズ-ヘルムホルツの式から直接読み取れます。

このような熱力学的考察が工業的な化学プロセスの温度条件設計の基盤となっているでしょう。

キルヒホッフの法則との組み合わせ

ΔHの温度依存性を考慮した正確な計算には、キルヒホッフの法則(熱容量差ΔCpによるΔHの温度補正)とギブズ-ヘルムホルツの式を組み合わせます。

キルヒホッフの法則:ΔH(T)= ΔH°₂₉₈ + ∫(298→T)ΔCp dT

これをギブズ-ヘルムホルツに代入:

ΔG(T₂)/T₂ = ΔG°₂₉₈/298 − ∫(298→T₂)(ΔH(T)/T²)dT

ΔCp = 一定の近似下:

ΔH(T)≈ ΔH°₂₉₈ + ΔCp(T − 298)

キルヒホッフ補正を加えることで、広い温度範囲にわたってより正確なΔG(T)の計算が可能となり、高温工業プロセス(製鉄・セメント・ガラス製造など)の熱力学設計精度が大幅に向上します。

現代の熱力学計算ソフトウェア(FactSage・HSC Chemistryなど)はこれらの補正を自動で処理してくれるため、実務では計算ツールを活用した効率的な解析が一般的でしょう。

まとめ

ギブズ-ヘルムホルツの式(∂(G/T)/∂T)p = −H/T² は、ギブズ自由エネルギーの温度依存性をエンタルピーと温度の関係で表した熱力学の重要な恒等式です。

ギブズ自由エネルギーの定義とエントロピーの関係式から厳密に導出でき、追加の仮定なしに成立する普遍的な関係式となっています。

ファント・ホッフの式・電池の起電力温度係数・化学反応の温度依存性計算など、化学熱力学・電気化学・生化学にわたる幅広い応用が確立されています。

ヘルムホルツ版(∂(F/T)/∂T)V = −U/T² は統計力学・物性物理の理論計算においてエネルギーとエントロピーの分離評価に役立つ重要な関係式でしょう。

実験データとの組み合わせではファント・ホッフプロットが標準的な解析手法として広く活用されており、反応の熱力学的特性を実験から導出する強力なツールとなっています。