熱力学的な系 の内部エネルギー
U
{\displaystyle U}
の微小変化
d
U
{\displaystyle dU}
は、熱力学第一法則 より、
d
U
=
d
′
Q
+
d
′
W
{\displaystyle dU=d'Q+d'W}
である。ここで、
d
′
Q
{\displaystyle d'Q}
は外部から系に流れる熱量、
d
′
W
{\displaystyle d'W}
は外部が系にする仕事である。系のエントロピー を
S
{\displaystyle S}
、熱力学的温度 を
T
{\displaystyle T}
とすると、
d
′
Q
=
T
d
S
{\displaystyle d'Q=TdS}
であるから、上式は、
d
U
=
T
d
S
+
d
′
W
{\displaystyle dU=TdS+d'W}
である。系が外部と粒子の出入りがない場合では、
d
′
W
{\displaystyle d'W}
は力学的仕事
d
′
W
M
=
−
P
d
V
{\displaystyle d'W_{M}=-PdV}
(
P
{\displaystyle P}
は系の圧力 、
V
{\displaystyle V}
は系の体積)に等しく、
d
U
=
T
d
S
−
P
d
V
{\displaystyle dU=TdS-PdV}
となる[ 2] が、粒子の出入りがある場合では、
d
′
W
{\displaystyle d'W}
に化学的仕事[ 3]
d
′
W
C
,
i
=
μ
i
d
N
i
{\displaystyle d'W_{C,i}=\mu _{i}dN_{i}}
が加わる(ここでi は粒子の種類(成分)を識別する記号である)。すなわち、
d
′
W
=
d
′
W
M
+
∑
i
d
′
W
C
,
i
{\displaystyle d'W=d'W_{M}+\sum _{i}d'W_{C,i}}
であり、
d
U
{\displaystyle dU}
は
d
U
=
T
d
S
−
P
d
V
+
∑
i
μ
i
d
N
i
{\displaystyle dU=TdS-PdV+\sum _{i}\mu _{i}dN_{i}}
(*)
となる[ 3] 。ここで、
d
N
i
{\displaystyle dN_{i}}
は成分i の物質量
N
i
{\displaystyle N_{i}}
の微小変化を表しており、
μ
i
{\displaystyle \mu _{i}}
は
μ
i
(
S
,
V
,
N
)
=
(
∂
U
∂
N
i
)
S
,
V
,
N
j
≠
i
{\displaystyle \mu _{i}(S,V,{\boldsymbol {N}})=\left({\frac {\partial U}{\partial N_{i}}}\right)_{S,V,N_{j\neq i}}}
で定義される、化学ポテンシャルと呼ばれる量[ 4] である。(全微分 も参照)ここで、括弧に付く添え字はその変数を一定にして偏微分することを意味する。また、
N
j
≠
i
{\displaystyle N_{j\neq i}}
は
i
{\displaystyle i}
以外の全ての成分の物質量を表す。
化学ポテンシャルは様々な変数の組の関数として、また様々な熱力学ポテンシャル の偏微分の関数として表現される。例えば、F を系のヘルムホルツエネルギー とすると、成分
i
{\displaystyle i}
の化学ポテンシャルは
μ
i
(
T
,
V
,
N
)
=
(
∂
F
∂
N
i
)
T
,
V
,
N
j
≠
i
{\displaystyle \mu _{i}(T,V,{\boldsymbol {N}})=\left({\frac {\partial F}{\partial N_{i}}}\right)_{T,V,N_{j\neq i}}}
と表される[ 5] 。これは次のようにして示される。まず、
F
=
U
−
T
S
{\displaystyle F=U-TS}
なので、その微小変化は、
d
F
=
d
U
−
S
d
T
−
T
d
S
{\displaystyle dF=dU-SdT-TdS}
である。ここで(*)のdU を代入すると
d
F
(
T
,
V
,
N
)
=
−
S
d
T
−
P
d
V
+
∑
i
μ
i
d
N
i
{\displaystyle dF\left(T,V,{\boldsymbol {N}}\right)=-SdT-PdV+\sum _{i}\mu _{i}dN_{i}}
であるから、上式が成り立つことが言える。ここで、U からF への変換はルジャンドル変換 となっている。同様に、系のギブズエネルギー G=U-TS+PV とエンタルピー H=U+PV に対して、
μ
i
(
T
,
P
,
N
)
=
(
∂
G
∂
N
i
)
T
,
P
,
N
j
≠
i
{\displaystyle \mu _{i}(T,P,{\boldsymbol {N}})=\left({\frac {\partial G}{\partial N_{i}}}\right)_{T,P,N_{j\neq i}}}
μ
i
(
S
,
P
,
N
)
=
(
∂
H
∂
N
i
)
S
,
P
,
N
j
≠
i
{\displaystyle \mu _{i}(S,P,{\boldsymbol {N}})=\left({\frac {\partial H}{\partial N_{i}}}\right)_{S,P,N_{j\neq i}}}
も示される。また、(*)より、エントロピー S の微小変化は
d
S
=
1
T
d
U
+
P
T
d
V
−
1
T
∑
i
μ
i
d
N
i
{\displaystyle dS={\frac {1}{T}}dU+{\frac {P}{T}}dV-{\frac {1}{T}}\sum _{i}\mu _{i}dN_{i}}
なので、
μ
i
(
U
,
V
,
N
)
=
−
T
(
∂
S
∂
N
i
)
U
,
V
,
N
j
≠
i
{\displaystyle \mu _{i}(U,V,{\boldsymbol {N}})=-T\left({\frac {\partial S}{\partial N_{i}}}\right)_{U,V,N_{j\neq i}}}
も成り立つ。
温度 T と圧力 p 、および物質量 の組 N = (N 1 , N 2 ,...) により平衡状態が指定される場合での化学ポテンシャルは
μ
i
(
N
;
T
,
p
)
=
(
∂
G
∂
N
i
)
T
,
p
,
N
j
≠
i
{\displaystyle \mu _{i}({\boldsymbol {N}};T,p)=\left({\frac {\partial G}{\partial N_{i}}}\right)_{T,p,N_{j\neq i}}}
で与えられる[ 6] 。
このように温度と圧力と、1成分を除いた物質量を固定した示量変数の偏微分は部分モル量(partial molal quantity )と呼ばれ[ 8] [ 9] 、この意味で化学ポテンシャルは部分モルギブズエネルギーに等しい。
系のスケール変換を考えれば、ギブスエネルギー と物質量の示量性 、及び温度と圧力の示強性 から、スケール・パラメータ λ に対して
λ
G
(
N
;
T
,
p
)
=
G
(
λ
N
;
T
,
p
)
{\displaystyle \lambda G({\boldsymbol {N}};T,p)=G(\lambda {\boldsymbol {N}};T,p)}
が成り立つ。これを λ について微分すれば
G
(
N
;
T
,
p
)
=
∑
i
N
i
∂
G
∂
N
i
|
N
=
λ
N
=
∑
i
N
i
μ
i
(
λ
N
;
T
,
p
)
{\displaystyle G({\boldsymbol {N}};T,p)=\sum _{i}N_{i}\left.{\frac {\partial G}{\partial N_{i}}}\right|_{N=\lambda N}=\sum _{i}N_{i}\,\mu _{i}(\lambda {\boldsymbol {N}};T,p)}
であり、λ = 1 と置けば
G
(
N
;
T
,
p
)
=
∑
i
N
i
μ
i
(
N
;
T
,
p
)
{\displaystyle G({\boldsymbol {N}};T,p)=\sum _{i}N_{i}\,\mu _{i}({\boldsymbol {N}};T,p)}
の関係が得られる。各成分の化学ポテンシャルとその成分の物質量の積の総和がギブズエネルギーとなる。
特に単一成分系では
G
(
N
;
T
,
p
)
=
N
μ
(
T
,
p
)
{\displaystyle G(N;T,p)=N\mu (T,p)}
であり、ギブズエネルギーは物質量に比例 し、化学ポテンシャルは物質量に依らない。
つまり1成分系では温度と圧力が等しければ化学ポテンシャルは等しい。これは自由に熱を通し自由に動くことができる壁に穴を開けても、平衡状態は変化しない(壁の両側でマクロな物質量は変化しない)ことを意味する[ 10] 。
理想気体 のモル体積は V m = RT /p であり、これを積分すると
μ
(
T
,
p
)
=
μ
∘
(
T
)
+
R
T
ln
p
p
∘
{\displaystyle \mu (T,p)=\mu ^{\circ }(T)+RT\ln {\frac {p}{p^{\circ }}}}
となる。ここで p ° は標準状態圧力 、μ ° は標準化学ポテンシャルである。
実在気体の場合はモル体積をビリアル展開 で
V
m
(
T
,
p
)
=
R
T
p
+
B
(
T
)
+
C
(
T
)
p
+
O
(
p
2
)
{\displaystyle V_{\text{m}}(T,p)={\frac {RT}{p}}+B(T)+C(T)\,p+O(p^{2})}
と表わしたものを積分すれば
μ
(
T
,
p
)
=
μ
∘
(
T
)
+
R
T
ln
p
p
∘
+
B
(
T
)
p
+
C
(
T
)
p
2
2
+
O
(
p
3
)
{\displaystyle \mu (T,p)=\mu ^{\circ }(T)+RT\ln {\frac {p}{p^{\circ }}}+B(T)\,p+C(T)\,{\frac {p^{2}}{2}}+O(p^{3})}
となる。標準化学ポテンシャルは
μ
∘
(
T
)
=
lim
p
→
0
{
μ
(
T
,
p
)
−
R
T
ln
p
p
∘
}
{\displaystyle \mu ^{\circ }(T)=\lim _{p\to 0}\left\{\mu (T,p)-RT\ln {\frac {p}{p^{\circ }}}\right\}}
で定義される[ 11] 。
また、フガシティー f を用いることで実在気体の化学ポテンシャルを
μ
(
T
,
p
)
=
μ
∘
(
T
)
+
R
T
ln
f
(
T
,
p
)
p
∘
{\displaystyle \mu (T,p)=\mu ^{\circ }(T)+RT\ln {\frac {f(T,p)}{p^{\circ }}}}
と表わすこともできる。ビリアル展開の形と比較すれば、フガシティーは
f
(
T
,
p
)
=
p
e
B
p
/
R
T
+
O
(
p
2
)
{\displaystyle f(T,p)=p\,e^{Bp/RT+O(p^{2})}}
である[ 11] 。
混合物の組成をモル分率 xi の組で表したとき、理想混合系の化学ポテンシャルは
μ
i
id
(
x
;
T
,
p
)
=
μ
i
∗
(
T
,
p
)
+
R
T
ln
x
i
{\displaystyle \mu _{i}^{\text{id}}({\boldsymbol {x}};T,p)=\mu _{i}^{*}(T,p)+RT\ln x_{i}}
で表される[ 12] 。ここで * は純物質 における量を表している。
特に理想混合気体では
μ
i
id
(
x
;
T
,
p
)
=
μ
i
∘
(
T
)
+
R
T
ln
x
i
p
p
∘
{\displaystyle \mu _{i}^{\text{id}}({\boldsymbol {x}};T,p)=\mu _{i}^{\circ }(T)+RT\ln {\frac {x_{i}\,p}{p^{\circ }}}}
となり、純粋系での圧力をその組成の分圧 xi p で置き換えた形となる。
実在気体の混合系では、分圧をフガシティーへ置き換えて表される。
実在の固溶体 や実在溶液 では、モル分率を活量 へ置き換えて
μ
i
(
x
;
T
,
p
)
=
μ
i
∗
(
T
,
p
)
+
R
T
ln
a
i
(
x
;
T
,
p
)
{\displaystyle \mu _{i}({\boldsymbol {x}};T,p)=\mu _{i}^{*}(T,p)+RT\ln a_{i}({\boldsymbol {x}};T,p)}
で表わされる[ 13] 。
理想溶液 において、溶質 i の濃度が質量モル濃度 bi で表されるときの化学ポテンシャルは
μ
i
(
b
;
T
,
p
)
=
μ
b
,
i
∘
(
T
,
p
)
+
R
T
ln
b
i
b
∘
{\displaystyle \mu _{i}({\boldsymbol {b}};T,p)=\mu _{b,i}^{\circ }(T,p)+RT\ln {\frac {b_{i}}{b^{\circ }}}}
で表される。ここで b ° は標準質量モル濃度であり、通常 b ° = 1 mol /kg に選ばれる。
溶質の濃度がモル濃度 ci で表されるときの化学ポテンシャルは
μ
i
(
c
;
T
,
p
)
=
μ
c
,
i
∘
(
T
,
p
)
+
R
T
ln
c
i
c
∘
{\displaystyle \mu _{i}({\boldsymbol {c}};T,p)=\mu _{c,i}^{\circ }(T,p)+RT\ln {\frac {c_{i}}{c^{\circ }}}}
で表される[ 14] 。ここで c ° は標準モル濃度であり、通常 c ° = 1 mol/L に選ばれる。
実在溶液の場合は活量 を用いることで、それぞれに
μ
i
(
b
;
T
,
p
)
=
μ
b
,
i
∘
(
T
,
p
)
+
R
T
ln
a
b
,
i
(
b
;
T
,
p
)
{\displaystyle \mu _{i}({\boldsymbol {b}};T,p)=\mu _{b,i}^{\circ }(T,p)+RT\ln a_{b,i}({\boldsymbol {b}};T,p)}
μ
i
(
c
;
T
,
p
)
=
μ
c
,
i
∘
(
T
,
p
)
+
R
T
ln
a
c
,
i
(
c
;
T
,
p
)
{\displaystyle \mu _{i}({\boldsymbol {c}};T,p)=\mu _{c,i}^{\circ }(T,p)+RT\ln a_{c,i}({\boldsymbol {c}};T,p)}
と表わすことができる[ 15] [ 14] 。
無限希釈の極限 b → 0 あるいは c → 0 で理想溶液に漸近するので、標準化学ポテンシャルは
μ
b
,
i
∘
(
T
,
p
)
=
lim
b
→
0
{
μ
i
(
b
;
T
,
p
)
−
R
T
ln
b
i
b
∘
}
{\displaystyle \mu _{b,i}^{\circ }(T,p)=\lim _{b\to 0}\left\{\mu _{i}({\boldsymbol {b}};T,p)-RT\ln {\frac {b_{i}}{b^{\circ }}}\right\}}
μ
c
,
i
∘
(
T
,
p
)
=
lim
c
→
0
{
μ
i
(
c
;
T
,
p
)
−
R
T
ln
c
i
c
∘
}
{\displaystyle \mu _{c,i}^{\circ }(T,p)=\lim _{c\to 0}\left\{\mu _{i}({\boldsymbol {c}};T,p)-RT\ln {\frac {c_{i}}{c^{\circ }}}\right\}}
で定義される[ 15] 。
溶質の濃度が質量モル濃度 b で表されるときの溶媒のモル分率は
x
slv
=
1
1
+
M
slv
∑
slt
b
{\displaystyle x_{\text{slv}}={\frac {1}{1+M_{\text{slv}}\sum _{\text{slt}}b}}}
なので、理想溶液における溶媒の化学ポテンシャルは
μ
slv
id
(
b
;
T
,
p
)
=
μ
slv
∗
(
T
,
p
)
−
R
T
ln
(
1
+
M
slv
∑
slt
b
i
)
{\displaystyle \mu _{\text{slv}}^{\text{id}}({\boldsymbol {b}};T,p)=\mu _{\text{slv}}^{*}(T,p)-RT\ln \left(1+M_{\text{slv}}\sum _{\text{slt}}b_{i}\right)}
である。実在溶液においては活量で置き換えて
μ
slv
(
b
;
T
,
p
)
=
μ
slv
∗
(
T
,
p
)
+
R
T
ln
a
slv
(
b
;
T
,
p
)
=
μ
slv
∗
(
T
,
p
)
−
ϕ
R
T
M
slv
∑
slt
b
i
{\displaystyle {\begin{aligned}\mu _{\text{slv}}({\boldsymbol {b}};T,p)&=\mu _{\text{slv}}^{*}(T,p)+RT\ln a_{\text{slv}}({\boldsymbol {b}};T,p)\\&=\mu _{\text{slv}}^{*}(T,p)-\phi RTM_{\text{slv}}\sum _{\text{slt}}b_{i}\end{aligned}}}
となる。
ここで
ϕ
(
b
;
T
,
p
)
=
−
ln
a
slv
(
b
;
T
,
p
)
M
slv
∑
slt
b
{\displaystyle \phi ({\boldsymbol {b}};T,p)=-{\frac {\ln a_{\text{slv}}({\boldsymbol {b}};T,p)}{M_{\text{slv}}\sum _{\text{slt}}b}}}
は浸透係数である。
このときギブズエネルギーは
G
=
N
slt
μ
slv
+
∑
slt
N
i
μ
i
=
N
slt
μ
slv
∗
+
∑
slt
N
i
{
μ
i
−
ϕ
R
T
}
=
N
slt
μ
slv
∗
+
∑
slt
N
i
μ
i
app
{\displaystyle {\begin{aligned}G&=N_{\text{slt}}\,\mu _{\text{slv}}+\sum _{\text{slt}}N_{i}\,\mu _{i}\\&=N_{\text{slt}}\,\mu _{\text{slv}}^{*}+\sum _{\text{slt}}N_{i}\left\{\mu _{i}-\phi RT\right\}\\&=N_{\text{slt}}\,\mu _{\text{slv}}^{*}+\sum _{\text{slt}}N_{i}\,\mu _{i}^{\text{app}}\end{aligned}}}
となる。ここで app は見かけの量 (apparent molar property ) を表している。
化学量論数 νi で表される化学反応 において、反応進行度 を ξ とすれば、物質量は
N
i
(
ξ
)
=
N
0
,
i
+
ν
i
ξ
{\displaystyle N_{i}(\xi )=N_{0,i}+\nu _{i}\,\xi }
と表わされる。等温等圧条件下ではギブズエネルギー が減少する方向に変化が進行し、平衡状態においてギブズエネルギーが極小となる。従って
Δ
r
G
≡
d
d
ξ
G
(
N
(
ξ
)
;
T
,
p
)
=
0
{\displaystyle \Delta _{\text{r}}G\equiv {\frac {d}{d\xi }}G({\boldsymbol {N}}(\xi );T,p)=0}
を満たす ξ において化学平衡となる。
反応のギブズエネルギーは化学ポテンシャルを用いて
Δ
r
G
=
∑
i
d
N
i
d
ξ
(
∂
G
∂
N
i
)
T
,
p
,
N
j
≠
i
=
∑
i
ν
i
μ
i
(
N
(
ξ
)
;
T
,
p
)
{\displaystyle \Delta _{\text{r}}G=\sum _{i}{\frac {dN_{i}}{d\xi }}\,\left({\frac {\partial G}{\partial N_{i}}}\right)_{T,p,N_{j\neq i}}=\sum _{i}\nu _{i}\,\mu _{i}({\boldsymbol {N}}(\xi );T,p)}
と書くことができて、理想混合気体においては
Δ
r
G
=
∑
i
ν
i
μ
i
∘
(
T
)
+
R
T
∑
i
ln
(
p
x
i
(
ξ
)
p
∘
)
ν
i
{\displaystyle \Delta _{\text{r}}G=\sum _{i}\nu _{i}\,\mu _{i}^{\circ }(T)+RT\sum _{i}\ln \left({\frac {px_{i}(\xi )}{p^{\circ }}}\right)^{\nu _{i}}}
となる。
標準平衡定数 を
R
T
ln
K
∘
(
T
;
ν
)
=
−
∑
i
ν
i
μ
i
∘
(
T
)
{\displaystyle RT\ln K^{\circ }(T;{\boldsymbol {\nu }})=-\sum _{i}\nu _{i}\,\mu _{i}^{\circ }(T)}
で定義すれば、平衡の条件は
K
∘
(
T
;
ν
)
=
∏
i
(
p
x
i
(
ξ
)
p
∘
)
ν
i
{\displaystyle K^{\circ }(T;{\boldsymbol {\nu }})=\prod _{i}\left({\frac {px_{i}(\xi )}{p^{\circ }}}\right)^{\nu _{i}}}
となる。