熱力学ポテンシャル (ねつりきがくポテンシャル、英語 : thermodynamic potential )または完全な熱力学関数 [ 1] [ 注釈 1] とは、熱力学 において、ある示量性 状態量を別の状態量の組を変数とする関数で表したもの[ 2] で、系の平衡状態 における熱力学的性質の情報を全て持つもの(後述)である。
例えば、T 、V 、N を変数とする(関数としての)ヘルムホルツエネルギーF(T,V,N) やT 、p 、N を変数とするギブズエネルギーG(T,p,N) は熱力学ポテンシャルである[ 3] 。(ここで、T は温度、V は体積、N は物質量、p は圧力)
「熱力学的性質の情報を全て持つ」とは全ての状態量がこの関数から(偏微分 等の組み合わせにより)与えられるという意味である。言い換えれば、完全な熱力学関数が与えられればそこから状態方程式 や熱容量などの系の性質が決まる[ 4] 。熱力学からは関数形に制約(凸性 など)を与えるが、具体的な関数形は実験的に決められるか、統計力学 から導出するなど、熱力学以外から与えられる[ 5] 。
熱力学ポテンシャルの一つである内部エネルギー U は、エントロピー S 、体積 V 、各成分の物質量 N = {N i } 、あるいはその他の示量性状態量[ 注釈 2] X を変数に持つ関数 U (S , N , V , X ) として表されたときに完全な熱力学関数となる。このことはエネルギー表示 と呼ばれることがある[ 6] 。このとき、各変数による偏微分は
(
∂
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 に対応する示強性状態量[ 注釈 3] x が得られる。これらはエネルギー表示の示強性状態量 [ 6] と呼ばれる。これらの状態量から熱容量 や圧縮率 なども計算される。内部エネルギーの全微分 は
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 ) として表されたときに完全な熱力学関数となる。このことをエントロピー表示 と呼ぶ場合がある[ 6] 。この表示でも熱力学的性質の情報としてはエネルギー表示と等価であり、各変数による偏微分を求めることによりエントロピー表示の示強変数 [ 6] 、たとえば
(
∂
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)}}}
などが得られる。しかし統計力学との関係からこの表示は重要である。ボルツマンの原理 により分配関数と直接的に関係づけられるのはエントロピー表示の熱力学ポテンシャルである。
ボルンの熱力学的正方形[ 7] (en:Thermodynamic square )は、ルジャンドル変換で繋がっている完全な熱力学関数とその変数を覚えやすくまとめたものである。
完全な熱力学関数には自然な独立変数 [ 5] の組があり、同じ状態量であっても、変数が異なればそれは完全な熱力学関数とはならない。例えば内部エネルギー 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}
などの完全な熱力学関数を導入することができる[ 注釈 4] 。なお、この関数はヘルムホルツエネルギーと Ψ = −F /T の関係にある[ 8] 。
熱力学ポテンシャルとその変数の例[ 8]
表示
熱力学ポテンシャル
記号と定義
自然な変数
全微分形
エネルギー表示
内部エネルギー
U
(S , V , N )
dU = TdS − pdV μdN
エンタルピー
H = U pV
(S , p , N )
dH = TdS Vdp μdN
ヘルムホルツエネルギー
F = U − TS
(T , V , N )
dF = − SdT − pdV μdN
ギブズエネルギー
G = F pV
(T , p , N )
dG = − SdT Vdp μdN
グランドポテンシャル
J = F − μN
(T , V , μ )
dJ = − SdT - pdV − Ndμ
エントロピー表示
エントロピー
S
(U , V , N )
dS = (1/T )dU (p /T )dV − (μ /T )dN
マシュー関数 (英語版 )
Ψ = S − U /T = −F /T
(β , V , N ) [ 注釈 4]
d Ψ = −Udβ (p /T )dV − (μ /T )dN
Planck関数
Φ = Ψ − (p /T )V = −G /T
(β , p /T , N ) [ 注釈 4]
d Φ = −Hdβ − (V /T )dp − (μ /T )dN
Kramers関数
q = Ψ αN = −J /T
(β , V , α ) [ 注釈 4]
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 , μ ) を系の平衡状態を指定する状態変数として選ぶことは出来ないことを表している。
系が温度 T ex の外界と接しているとき、熱力学第二法則 から、系に変化が起きるとき
δ
′
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}
である。p ex は外界の圧力、μ 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}
などが得られる[ 9] 。ここで、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 はボルツマン定数 である。
この節の
加筆 が望まれています。
主に: 理想気体や光子気体などの例 (2020年12月 )
理想気体 のエントロピーは
S
[
U
,
N
,
V
]
=
N
σ
∗
N
c
R
ln
U
−
N
μ
∗
N
c
R
T
∗
N
R
ln
V
N
R
T
∗
/
p
∘
{\displaystyle S[U,N,V]=N\sigma ^{*} NcR\ln {\frac {U-N\mu ^{*}}{NcRT^{*}}} NR\ln {\frac {V}{NRT^{*}/p^{\circ }}}}
で与えられる。ここで σ * , T * , μ * , p ° はそれぞれエントロピー、温度、化学ポテンシャル、圧力の基準を与える適当な定数であり、積分定数の任意性を含んでいる。
また R はモル気体定数 であり、c は理想気体を特徴付ける唯一のパラメータで、気体粒子の自由度 の半分に相当する。
エネルギーと物質量、体積による偏微分が
1
T
=
∂
S
∂
U
=
N
c
R
U
−
N
μ
∗
{\displaystyle {\frac {1}{T}}={\frac {\partial S}{\partial U}}={\frac {NcR}{U-N\mu ^{*}}}}
μ
T
=
−
∂
S
∂
N
=
μ
∗
T
−
σ
∗
(
c
1
)
R
−
c
R
ln
U
−
N
μ
∗
N
c
R
T
∗
−
R
ln
V
N
R
T
∗
/
p
∘
{\displaystyle {\frac {\mu }{T}}=-{\frac {\partial S}{\partial N}}={\frac {\mu ^{*}}{T}}-\sigma ^{*} (c 1)R-cR\ln {\frac {U-N\mu ^{*}}{NcRT^{*}}}-R\ln {\frac {V}{NRT^{*}/p^{\circ }}}}
p
T
=
∂
S
∂
V
=
N
R
V
{\displaystyle {\frac {p}{T}}={\frac {\partial S}{\partial V}}={\frac {NR}{V}}}
となるので、理想気体の関係式として
U
=
N
μ
∗
N
c
R
T
{\displaystyle U=N\mu ^{*} NcRT}
p
V
=
N
R
T
{\displaystyle pV=NRT}
μ
=
μ
∗
−
T
σ
∗
(
c
1
)
R
T
(
1
−
ln
T
T
∗
)
R
T
ln
p
p
∘
{\displaystyle \mu =\mu ^{*}-T\sigma ^{*} (c 1)RT\left(1-\ln {\frac {T}{T^{*}}}\right) RT\ln {\frac {p}{p^{\circ }}}}
が導かれる。
ファン・デル・ワールス気体 のエントロピーは
S
[
U
,
N
,
V
]
=
N
σ
∗
N
c
R
ln
U
−
N
μ
∗
N
2
a
/
V
N
c
R
T
∗
N
R
ln
V
−
N
b
N
R
T
∗
/
p
∘
{\displaystyle S[U,N,V]=N\sigma ^{*} NcR\ln {\frac {U-N\mu ^{*} N^{2}a/V}{NcRT^{*}}} NR\ln {\frac {V-Nb}{NRT^{*}/p^{\circ }}}}
で与えられる。二つのファン・デル・ワールス定数 a, b を導入して理想気体からのずれが表現されている。
偏微分が
1
T
=
∂
S
∂
U
=
N
c
R
U
−
N
μ
∗
N
2
a
/
V
{\displaystyle {\frac {1}{T}}={\frac {\partial S}{\partial U}}={\frac {NcR}{U-N\mu ^{*} N^{2}a/V}}}
p
T
=
∂
S
∂
V
=
N
R
V
−
N
b
−
1
T
a
(
V
/
N
)
2
{\displaystyle {\frac {p}{T}}={\frac {\partial S}{\partial V}}={\frac {NR}{V-Nb}}-{\frac {1}{T}}{\frac {a}{(V/N)^{2}}}}
となるので
U
=
N
μ
∗
N
c
R
T
−
N
a
V
/
N
{\displaystyle U=N\mu ^{*} NcRT-{\frac {Na}{V/N}}}
p
=
R
T
(
V
/
N
)
−
b
−
a
(
V
/
N
)
2
{\displaystyle p={\frac {RT}{(V/N)-b}}-{\frac {a}{(V/N)^{2}}}}
が導かれる。
^ ウィラード・ギブズ は基本的な方程式 (fundamental equations)と呼んでいた
^ 例えば分極 P や磁化 M など。
^ 例えば分極 P に対応する外部電場 E や磁化 M に対応する外部磁場 H など。
^ a b c d β =1/T 、α =μ /T を表す。
^ 田崎『熱力学 現代的な視点から』 p.120
^ 田崎『熱力学 現代的な視点から』p.32、ただし、ここでは完全かどうか、または表す状態量が示量性かどうかに言及していない、一般的な熱力学関数の説明である
^ 田崎『熱力学 現代的な視点から』Fについてはp.120、Gについてはp.157
^ 佐々真一 著、兵頭俊夫 編『熱力学入門』共立出版、2000年、76頁。ISBN 4-320-03347-7 。
^ a b 久保『熱学・統計力学』 p.88
^ a b c d 清水『熱力学の基礎I』pp.94-98
^ 清水『熱力学の基礎II』p.2
^ a b 久保『熱学・統計力学』 p.90
^ 久保『熱学・統計力学』 p.97