化学ポテンシャル (かがくポテンシャル、英語 : chemical potential )は、熱力学 で用いられる示強性 状態量 の一つで、浸透圧 や相平衡 、化学反応 のようなマクロな物質量の移動が伴う現象で重要となる物理量 である。 推奨される量記号は、μ (ミュー)である。
化学ポテンシャルの概念および記号と用語は、ウィラード・ギブズ の1876年の論文『不均一な物質系の平衡に就いて 』で導入された。
化学ポテンシャルは、物質の多寡により系が潜在的に持つエネルギー の大きさの尺度となる量である。 例えば、半透膜 で隔てられた二つの系の間に濃度 差が有った場合、浸透圧 が生じ仕事 を為す事が出来る。 また、物質が増減する化学反応 では熱の出入り(発熱反応 、吸熱反応 )を伴う。 このように、物質が存在することにより系は潜在的にエネルギーを持つ。 その系 に含まれるある成分の単位物質量 あたりのギブスエネルギー がその成分の化学ポテンシャルに相当する。
熱力学的な系 の内部エネルギー 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] 。
温度T 、圧力p を変数としたときの化学ポテンシャルの偏微分は
( ∂ μ ∂ T ) p = 1 N ( ∂ G ∂ T ) p , N = − S N {\displaystyle \left({\frac {\partial \mu }{\partial T}}\right)_{p}={\frac {1}{N}}\left({\frac {\partial G}{\partial T}}\right)_{p,N}=-{\frac {S}{N}}}
( ∂ μ ∂ p ) T = 1 N ( ∂ G ∂ p ) T , N = V N {\displaystyle \left({\frac {\partial \mu }{\partial p}}\right)_{T}={\frac {1}{N}}\left({\frac {\partial G}{\partial p}}\right)_{T,N}={\frac {V}{N}}}
となる。
「ギブズ・デュエムの式 」も参照
理想気体 のモル体積は 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}}}
となる。
モル数でなく、粒子 数としての化学ポテンシャルμ も考えることができる。固体電子論 における電子系(例:電子ガス )でも化学ポテンシャルを定義することができ、特に温度T = 0 Kにおける化学ポテンシャル μ のことを、フェルミエネルギー ε F と呼ぶ場合がある。
μ T = 0 = ϵ F {\displaystyle \mu _{T=0}=\epsilon _{F}}
^ 清水 明『熱力学の基礎』p.135 ^ a b 清水 明『熱力学の基礎』p.158 ただし、多成分の場合に拡張して書いてある。 ^ 清水 明『熱力学の基礎』p.107 ^ 田崎 p.174 ^ バーロー『物理化学(上)』 pp.233-235, §8.10 ^ バーロー『物理化学(上)』 p.163 ^ Kirkwood & Oppenheim p.9 ^ 佐々真一 『熱力学入門』共立出版 、2000年。ISBN 978-4320033474 。 ^ a b Kirkwood & Oppenheim p.89-92, §7-6. ^ バーロー『物理化学(下)』 p.612 ^ バーロー『物理化学(下)』 pp.619-621, §19.4 ^ a b バーロー『物理化学(下)』 pp.621-624, §19.5 ^ a b Kirkwood & Oppenheim p.160-164, §11-1. Gibbs, J. Willard (1876-05-01). “On the Equilibrium of Heterogeneous Substances. First Part.”. the Transactions of the Connecticut Academy of Arts and Sciences (Connecticut Academy of Arts and Sciences) 3 : 108–248. J.G. Kirkwood, I. Oppenheim (1961). Chemical Thermodynamics . McGraw-Hill G.M.バーロー『物理化学』 上巻、藤代亮一 訳(第1版)、東京化学同人、1968年。 G.M.バーロー『物理化学』 下巻、藤代亮一 訳(第1版)、東京化学同人、1968年。 田崎晴明『熱力学 - 現代的な視点から』培風館〈新物理学シリーズ〉、2000年。ISBN 4-563-02432-5 。 清水明『熱力学の基礎』東京大学出版会、2007年。ISBN 978-4-13-062609-5 。