化学ポテンシャル (かがくポテンシャル、英語 : chemical potential )は、熱力学 で用いられる示強性 状態量 の一つで、浸透圧 や相平衡 、化学反応 のようなマクロな物質量の移動が伴う現象で重要となる物理量 である。 推奨される量記号は、μ (ミュー)である。
化学ポテンシャルはアメリカ の化学者 ウィラード・ギブズ により導入された概念である。
化学ポテンシャルは、物質の多寡により系が潜在的に持つエネルギー の大きさの尺度となる量である。 例えば、半透膜 で隔てられた二つの系の間に濃度 差が有った場合、浸透圧 が生じ仕事 を為す事が出来る。 また、物質が増減する化学反応 では熱の出入り(発熱反応 、吸熱反応 )を伴う。 このように、物質が存在することにより系は潜在的にエネルギーを持つ。 その系 に含まれるある成分の単位物質量 あたりのギブスエネルギー がその成分の化学ポテンシャルに相当する。
温度 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}}}
で定義される[ 1] 。ここで括弧に付く添え字はその変数を固定して偏微分することを意味しており、N j ≠i は成分 i とは異なる残りの成分を固定して偏微分することを意味している。 このように温度と圧力と、1成分を除いた物質量を固定した示量変数の偏微分は部分モル量(partial molal quantity )と呼ばれ[ 3] [ 4] 、この意味で化学ポテンシャルは部分モルギブズエネルギーに等しい。
この定義においてギブズエネルギー G はT, p, N を変数とする関数であり、このときギブズエネルギーは熱力学ポテンシャル (完全な熱力学関数)である。変数が異なる場合にギブズエネルギーの偏微分は化学ポテンシャルとならず、別の熱力学ポテンシャルの偏微分として表される。いずれの定義であっても化学ポテンシャルとしては同じである。
温度 T と体積 V 、および物質量の組 N を自然な変数とする熱力学ポテンシャルはヘルムホルツエネルギー F であり、化学ポテンシャルが
μ i ( V , N ; T ) = ( ∂ F ∂ N i ) T , V , N j ≠ i {\displaystyle \mu _{i}(V,{\boldsymbol {N}};T)=\left({\frac {\partial F}{\partial N_{i}}}\right)_{T,V,N_{j\neq i}}}
で与えられる[ 5] [ 6] 。 また、内部エネルギー U は S, V, N を、エンタルピー H は S, p, N を、それぞれに自然な変数とする熱力学ポテンシャルであり、化学ポテンシャルが
μ 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}}}
μ i ( S , N ; p ) = ( ∂ H ∂ N i ) S , p , N j ≠ i {\displaystyle \mu _{i}(S,{\boldsymbol {N}};p)=\left({\frac {\partial H}{\partial N_{i}}}\right)_{S,p,N_{j\neq i}}}
で与えられる[ 6] 。
エントロピー S は内部エネルギー U,V,N を自然な変数とする熱力学ポテンシャルである。化学ポテンシャルは
μ 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}}}
で与えられる。
系の物質量が変化するとき、内部エネルギーU やギブズエネルギーG の全微分において、化学ポテンシャルが現れる。内部エネルギーの全微分は
d U = T d S − p d V + ∑ i μ i d N i {\displaystyle dU=T\,dS-p\,dV+\sum _{i}\mu _{i}\,dN_{i}}
であり[ 6] [ 7] 、ルジャンドル変換により内部エネルギーからギブズエネルギーへ変換すれば、全微分は
d G = − S d T + V d p + ∑ i μ i d N i {\displaystyle dG=-S\,dT+V\,dp+\sum _{i}\mu _{i}\,dN_{i}}
となる[ 1] [ 6] [ 8] 。(化学ポテンシャルは物質量 に共役 な示強変数である[ 9] 。)
化学ポテンシャルの物理的な意味は、同じ示強性である圧力との対応を考えるとわかりやすい。 [独自研究? ] たとえば圧力 (示強性)は、熱力学的な系の体積 (示量性)を少し変えたときに外界が感じる『手ごたえ』 [要追加記述 ] である。この関係性を化学ポテンシャルに当てはめてみると、化学ポテンシャル(示強性)とは、熱力学的な系の物質量(示量性)を少し変えたときの『手ごたえ』と考えることができる[ 5] 。よって [要追加記述 ] 平衡状態に向かうときは、化学ポテンシャルが等しくなるように物質量は移動する[ 10] 。
また電磁気学 において電荷 q とその移動を司る静電ポテンシャル φ との積がポテンシャルエネルギー qφ である。この関係性を化学ポテンシャルに当てはめてみると[ 11] 、マクロな物質量N の移動を司るポテンシャルが化学ポテンシャルμ であり[ 12] 、それらの積であるギブズエネルギーNμ はポテンシャルエネルギーのような量だと考えることもできる。 [誰? ] ただし実際にはミクロな粒子間にある複雑な相互作用などの結果としてマクロな化学ポテンシャルは決まると考えられ、力学におけるポテンシャル [要追加記述 ] と熱力学における化学ポテンシャルはかなり異なり同一視することはできない[ 12] 。
系のスケール変換を考えれば、ギブスエネルギー と物質量の示量性 、及び温度と圧力の示強性 から、スケール・パラメータ λ に対して
λ 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\}}
で定義される[ 13] 。 また、フガシティー 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})}}
である[ 13] 。
混合物の組成をモル分率 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}}
で表される[ 14] 。ここで * は純物質 における量を表している。
特に理想混合気体では
μ 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)}
で表わされる[ 15] 。
理想溶液 において、溶質 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 }}}}
で表される[ 16] 。ここで 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)}
と表わすことができる[ 17] [ 16] 。 無限希釈の極限 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\}}
で定義される[ 17] 。
溶質の濃度が質量モル濃度 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}}
^ a b バーロー『物理化学(上)』 pp.233-235, §8.10 ^ バーロー『物理化学(上)』 p.163 ^ Kirkwood & Oppenheim p.9 ^ a b 田崎『熱力学』 p.120 ^ a b c d Kirkwood & Oppenheim p.52-53, §6-2. ^ 清水『熱力学の基礎』 p.110 ^ 田崎『熱力学』 p.178 ^ 清水『熱力学の基礎』 p.107 ^ a b 佐々真一 『熱力学入門』共立出版 、2000年。ISBN 978-4320033474 。 ^ 田崎『熱力学』 ^ a b 清水『熱力学の基礎』 ^ 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. 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 。 北孝文、鈴木久男(監修)『演習しよう熱・統計力学』数理工学社、2018年。ISBN 978-4-86481-053-1 。