熱力学ポテンシャル （ねつりきがくポテンシャル、英語 : 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}}}}
が導かれる。
![{\displaystyle S[U,N,V]=N\sigma ^{*}+NcR\ln {\frac {U-N\mu ^{*}}{NcRT^{*}}}+NR\ln {\frac {V}{NRT^{*}/p^{\circ }}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ec71c30926c6bb94b05542c61f326cfabbe821f3)
![{\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 }}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1de77e2ec4854d195023c28d0535cb7c9d4c0c14)
