熱力学ポテンシャル

熱力学
古典的カルノー熱機関(英語版)
状態(英語版)
過程(英語版)
サイクル
系の特性
注: 斜体は共役変数(英語版)を示す。
状態の関数
過程関数(英語版)
材料特性(英語版)
比熱容量  c = {\displaystyle c=}
T {\displaystyle T} S {\displaystyle \partial S}
N {\displaystyle N} T {\displaystyle \partial T}
圧縮率  β = {\displaystyle \beta =-}
1 {\displaystyle 1} V {\displaystyle \partial V}
V {\displaystyle V} p {\displaystyle \partial p}
熱膨張  α = {\displaystyle \alpha =}
1 {\displaystyle 1} V {\displaystyle \partial V}
V {\displaystyle V} T {\displaystyle \partial T}
方程式(英語版)
熱力学ポテンシャル
  • 歴史
  • 文化
歴史
  • 一般(英語版)
  • (英語版)
  • エントロピー(英語版)
  • 気体の法則
  • 永久機関(英語版)
哲学(英語版)
理論
重要文献
年表
  • 芸術
  • 教育
  • マクスウェルの熱力学的表面(英語版)
  • エネルギー拡散としてのエントロピー(英語版)
統計力学


熱力学 · 気体分子運動論
粒子統計
マクスウェル=ボルツマン

ボース=アインシュタイン
フェルミ=ディラック
パラ · エニオン · 組み紐(英語版)

アンサンブル
ミクロカノニカルアンサンブル

カノニカルアンサンブル
グランドカノニカルアンサンブル
等温定圧アンサンブル
等エンタルピー-定圧

熱力学
気体の法則(英語版) · カルノーサイクル

デュロン=プティの法則

模型
デバイ · アインシュタイン · イジング
熱力学ポテンシャル
内部エネルギー
エンタルピー
ヘルムホルツの自由エネルギー
ギブズの自由エネルギー
グランドポテンシャル
科学者
マクスウェル · ギブズ · ボルツマン · アインシュタイン · オンサーガー · ウィルソン · 久保亮五 · カダノフ · フィッシャー · 川崎恭治 · パリージ · エドワーズ · ローレンツ · 蔵本由紀 · ジャルジンスキー

熱力学ポテンシャル(ねつりきがくポテンシャル、英語: thermodynamic potential)とは、熱力学において、系の平衡状態における熱力学的性質の情報を全て持つ示量性状態量である。完全な熱力学関数とも呼ばれる[1]

ウィラード・ギブズ基本的な方程式 (fundamental equations)と呼んでいた[2]

概要

「熱力学的性質の情報を全て持つ」とは全ての状態量がこの関数から(偏微分等の組み合わせにより)与えられるという意味である。言い換えれば、完全な熱力学関数が与えられればそこから状態方程式や熱容量などの系の性質が決まる[3]。熱力学からは関数形に制約(凸性など)を与えるが、具体的な関数形は実験的に決められるか、統計力学から導出するなど、熱力学以外から与えられる[4]

熱力学ポテンシャルの一つである内部エネルギー U は、エントロピー S体積 V、各成分の物質量 N = {Ni}、あるいはその他の示量性状態量[注釈 1] X を変数に持つ関数 U(S, N, V, X) として表されたときに完全な熱力学関数となる。このことはエネルギー表示と呼ばれることがある[5]。このとき、各変数による偏微分は

( U S ) V , N , X = T ( S , V , N , X ) , ( U V ) S , N , X = p ( S , V , N , X ) , ( U N i ) S , V , X = μ i ( S , V , N , X ) , ( U X ) S , V , N = x ( S , V , N , X ) {\displaystyle {\begin{aligned}\left({\frac {\partial U}{\partial S}}\right)_{V,N,X}&=T(S,V,N,X),\\\left({\frac {\partial U}{\partial V}}\right)_{S,N,X}&=-p(S,V,N,X),\\\left({\frac {\partial U}{\partial N_{i}}}\right)_{S,V,X}&=\mu _{i}(S,V,N,X),\\\left({\frac {\partial U}{\partial X}}\right)_{S,V,N}&=x(S,V,N,X)\end{aligned}}}

となり、熱力学温度 T圧力 p、各成分の化学ポテンシャル μ = {μi}、及び X に対応する示強性状態量[注釈 2] x が得られる。これらはエネルギー表示の示強性状態量[5]と呼ばれる。これらの状態量から熱容量圧縮率なども計算される。内部エネルギーの全微分

d U = T d S p d V + i μ i d N i + x d X {\displaystyle dU=T\,dS-p\,dV+\sum _{i}\mu _{i}\,dN_{i}+x\,dX}

となる。

エネルギー表示の U の全微分の式を変形すれば

d S = 1 T d U + p T d V i μ i T d N i x T d X {\displaystyle dS={\frac {1}{T}}\,dU+{\frac {p}{T}}\,dV-\sum _{i}{\frac {\mu _{i}}{T}}\,dN_{i}-{\frac {x}{T}}\,dX}

となり、エントロピー S は内部エネルギー U、体積 V、物質量 N、及びその他の示量性状態量 X を変数にもつ関数 S(U, V, N, X) として表されたときに完全な熱力学関数となる。このことをエントロピー表示と呼ぶ場合がある[5]。この表示でも熱力学的性質の情報としてはエネルギー表示と等価であり、各変数による偏微分を求めることによりエントロピー表示の示強変数[5]、たとえば

( S U ) V , N , X = 1 T ( U , V , N , X ) {\displaystyle \left({\frac {\partial S}{\partial U}}\right)_{V,N,X}={\frac {1}{T(U,V,N,X)}}}

などが得られる。しかし統計力学との関係からこの表示は重要である。ボルツマンの原理により分配関数と直接的に関係づけられるのはエントロピー表示の熱力学ポテンシャルである。

ルジャンドル変換

ボルンの熱力学的正方形[6](en:Thermodynamic square)は、ルジャンドル変換で繋がっている完全な熱力学関数とその変数を覚えやすくまとめたものである。

完全な熱力学関数には自然な独立変数[4]の組があり、同じ状態量であっても、変数が異なればそれは完全な熱力学関数とはならない。例えば内部エネルギー U は、エントロピー S に替えて温度 T を変数に持つときには完全な熱力学関数とはならない。系の平衡状態を指定する状態変数の組が (T, V, N, X) である場合は、ルジャンドル変換

d F = d ( U T S ) = S d T p d V + i μ i d N i + x d X {\displaystyle dF=d(U-TS)=-S\,dT-p\,dV+\sum _{i}\mu _{i}\,dN_{i}+x\,dX}

によってヘルムホルツエネルギー F(T, V, N, X) が完全な熱力学関数となる。

エントロピーに対してもルジャンドル変換を考えることができて

d Ψ = d ( S β U ) = U d β p T d V + i μ i T d N i + x T d X {\displaystyle d\Psi =d(S-\beta U)=-U\,d\beta -{\frac {p}{T}}\,dV+\sum _{i}{\frac {\mu _{i}}{T}}\,dN_{i}+{\frac {x}{T}}\,dX}

などの完全な熱力学関数を導入することができる[注釈 3]。なお、この関数はヘルムホルツエネルギーと Ψ = −F/T の関係にある[7]

熱力学ポテンシャルとその変数の例[7]
表示 熱力学ポテンシャル 記号と定義 自然な変数 全微分形
エネルギー表示 内部エネルギー U (S, V, N) dU = TdSpdV + μdN
エンタルピー H = U + pV (S, p, N) dH = TdS + Vdp + μdN
ヘルムホルツエネルギー F = UTS (T, V, N) dF = − SdTpdV + μdN
ギブズエネルギー G = F + pV (T, p, N) dG = − SdT + Vdp + μdN
グランドポテンシャル J = FμN (T, V, μ) dJ = − SdT - pdVNdμ
エントロピー表示 エントロピー S (U, V, N) dS = (1/T)dU + (p/T)dV − (μ/T)dN
マシュー関数(英語版) Ψ = SU/T = −F/T (β, V, N)[注釈 3] dΨ = −Udβ + (p/T)dV − (μ/T)dN
Planck関数 Φ = Ψ − (p/T)V = −G/T (β, p/T, N)[注釈 3] dΦ = −Hdβ − (V/T)dp − (μ/T)dN
Kramers関数 q = Ψ + αN = −J/T (β, V, α)[注釈 3] dq = − Udβ + (p/T)dV + Ndα

ギブズ-デュエムの関係

詳細は「ギブズ-デュエムの式」を参照

系のスケール変換を考えると、内部エネルギー U、エントロピー S、体積 V、物質量 N の示量性から、任意の λ > 0 に対し

U ( λ S , λ V , λ N ) = λ U ( S , V , N ) {\displaystyle U(\lambda S,\lambda V,\lambda N)=\lambda U(S,V,N)}

という1次同次性が成り立つ。このことから

U = U S S + U V V + i U N i N i = T S p V + i μ i N i {\displaystyle U={\frac {\partial U}{\partial S}}S+{\frac {\partial U}{\partial V}}V+\sum _{i}{\frac {\partial U}{\partial N_{i}}}N_{i}=TS-pV+\sum _{i}\mu _{i}N_{i}}

の関係が導かれる。他にも

H = U + p V = T S + i μ i N i , F = U T S = p V + i μ i N i , G = F + p V = i μ i N i , J = F i μ i N i = p V {\displaystyle {\begin{aligned}H&=U+pV=TS+\sum _{i}\mu _{i}N_{i},\\F&=U-TS=-pV+\sum _{i}\mu _{i}N_{i},\\G&=F+pV=\sum _{i}\mu _{i}N_{i},\\J&=F-\sum _{i}\mu _{i}N_{i}=-pV\end{aligned}}}

などの関係式が得られる。

また、この式を微分すると

S d T V d p + i N i d μ i = 0 {\displaystyle S\,dT-V\,dp+\sum _{i}N_{i}\,d\mu _{i}=0}

の関係式が得られる。この関係式をギブズ-デュエムの関係と言い、示強性状態量の組 (T, p, μ) を系の平衡状態を指定する状態変数として選ぶことは出来ないことを表している。

平衡状態の安定性

系が温度 Tex の外界と接しているとき、熱力学第二法則から、系に変化が起きるとき

δ Q T ex δ S {\displaystyle \delta 'Q\leq T_{\text{ex}}\delta S}

である。一方、エネルギー保存則から

δ Q = δ U + δ W = δ U + p ex δ V μ ex δ N {\displaystyle \delta 'Q=\delta U+\delta 'W=\delta U+p_{\text{ex}}\delta V-\mu _{\text{ex}}\delta N}

である。pex は外界の圧力、μex は外界の化学ポテンシャルである。これらをまとめると、

δ U T e x δ S + p e x δ V μ e x δ N 0 {\displaystyle \delta U-T_{\mathrm {ex} }\delta S+p_{\mathrm {ex} }\delta V-\mu _{\mathrm {ex} }\delta N\leq 0}

となる。系が平衡状態にあるとき、変化が起こらないので、

δ U T e x δ S + p e x δ V μ e x δ N > 0 {\displaystyle \delta U-T_{\mathrm {ex} }\delta S+p_{\mathrm {ex} }\delta V-\mu _{\mathrm {ex} }\delta N>0}

である。 これが成り立つ条件は、1次変分について

( U S T e x ) δ S + ( U V + p e x ) δ V + ( U N μ e x ) δ N = 0 {\displaystyle \left({\frac {\partial U}{\partial S}}-T_{\mathrm {ex} }\right)\delta S+\left({\frac {\partial U}{\partial V}}+p_{\mathrm {ex} }\right)\delta V+\left({\frac {\partial U}{\partial N}}-\mu _{\mathrm {ex} }\right)\delta N=0}

および、2次変分について

1 2 ( 2 U S 2 δ S 2 + 2 U V 2 δ V 2 + 2 U N 2 δ N 2 ) + 2 U S V δ S δ V + 2 U S N δ S δ N + 2 U V N δ V δ N > 0 {\displaystyle {\frac {1}{2}}\left({\frac {\partial ^{2}U}{\partial S^{2}}}\delta S^{2}+{\frac {\partial ^{2}U}{\partial V^{2}}}\delta V^{2}+{\frac {\partial ^{2}U}{\partial N^{2}}}\delta N^{2}\right)+{\frac {\partial ^{2}U}{\partial S\partial V}}\delta S\delta V+{\frac {\partial ^{2}U}{\partial S\partial N}}\delta S\delta N+{\frac {\partial ^{2}U}{\partial V\partial N}}\delta V\delta N>0}

である。

1次変分の条件から

  • 系が熱を交換するとき(δS ≠ 0)、 U S T = T ex {\displaystyle {\frac {\partial U}{\partial S}}\equiv T=T_{\text{ex}}}
  • 系の体積が変化するとき(δV ≠ 0)、 U V p = p ex {\displaystyle {\frac {\partial U}{\partial V}}\equiv -p=-p_{\text{ex}}}
  • 系が物質を交換するとき(δN ≠ 0)、 U N μ = μ ex {\displaystyle {\frac {\partial U}{\partial N}}\equiv \mu =\mu _{\text{ex}}}

などの平衡条件が得られる。

2次変分の条件からは

2 U S 2 = ( T S ) V = T C V 0 {\displaystyle {\frac {\partial ^{2}U}{\partial S^{2}}}=\left({\frac {\partial T}{\partial S}}\right)_{V}={\frac {T}{C_{V}}}\geq 0}

2 U V 2 = ( p V ) S = 1 V κ S 0 {\displaystyle {\frac {\partial ^{2}U}{\partial V^{2}}}=-\left({\frac {\partial p}{\partial V}}\right)_{S}={\frac {1}{V\kappa _{S}}}\geq 0}

2 U S 2 2 U V 2 ( 2 U S V ) 2 = T C p 1 V κ S = T C V 1 V κ T 0 {\displaystyle {\frac {\partial ^{2}U}{\partial S^{2}}}{\frac {\partial ^{2}U}{\partial V^{2}}}-\left({\frac {\partial ^{2}U}{\partial S\partial V}}\right)^{2}={\frac {T}{C_{p}}}{\frac {1}{V\kappa _{S}}}={\frac {T}{C_{V}}}{\frac {1}{V\kappa _{T}}}\geq 0}

などが得られる[8]。ここで、CV定積熱容量Cp定圧熱容量κS断熱圧縮率κT等温圧縮率である。

平衡状態の安定性から熱力学ポテンシャルは一般に凸関数となる。

統計力学との関係

完全な熱力学関数は分配関数と関係付けられる。巨視的な熱力学と微視的な統計力学を結びつける関係である。

S ( E , N , V ) = k ln W ( E , N , V ) {\displaystyle S(E,N,V)=k\ln W(E,N,V)}

F ( β , N , V ) = 1 β ln Z ( β , N , V ) {\displaystyle F(\beta ,N,V)=-{\frac {1}{\beta }}\ln Z(\beta ,N,V)}

J ( β , μ , V ) = 1 β ln Ξ ( β , μ , V ) {\displaystyle J(\beta ,\mu ,V)=-{\frac {1}{\beta }}\ln \Xi (\beta ,\mu ,V)}

ここで、β=1/kT逆温度である。kボルツマン定数である。

具体例

脚注

[脚注の使い方]

注釈

  1. ^ 例えば分極 P磁化 M など。
  2. ^ 例えば分極 P に対応する外部電場 E や磁化 M に対応する外部磁場 H など。
  3. ^ a b c d β=1/Tα=μ/T を表す。

出典

  1. ^ 田崎『熱力学』 p.120
  2. ^ 田崎『熱力学』 p.17
  3. ^ 佐々真一 著、兵頭俊夫 編『熱力学入門』共立出版、2000年、76頁。ISBN 4-320-03347-7。 
  4. ^ a b 久保『熱学・統計力学』 p.88
  5. ^ a b c d 清水『熱力学の基礎I』pp.94-98
  6. ^ 清水『熱力学の基礎II』p.2
  7. ^ a b 久保『熱学・統計力学』 p.90
  8. ^ 久保『熱学・統計力学』 p.97

参考文献

  • 久保亮五『大学演習 熱学・統計力学 修訂版』裳華房、1998年。ISBN 978-4-7853-8032-8。 
  • 田崎晴明『熱力学 現代的な視点から』培風館〈新物理学シリーズ〉、2000年。ISBN 4-563-02432-5。 
  • 清水明『熱力学の基礎I』(2版)東京大学出版会、2021年。ISBN 978-4-13-062622-4。 
  • 清水明『熱力学の基礎II』(2版)東京大学出版会、2021年。ISBN 978-4-13-062623-1。 
典拠管理データベース ウィキデータを編集
全般
  • FAST
国立図書館
  • ドイツ
  • イスラエル
  • アメリカ