[論文レビュー] Individual ion species chemical potentials in the Mean Spherical Approximation

Johan S. Høye, Dirk Gillespie|arXiv (Cornell University)|May 27, 2022

thermodynamics and calorimetric analyses参考文献 13被引用数 7

ひとこと要約

この論文は、個々のイオン種の化学ポテンシャルに欠落していた平均球近似（MSA）における項を導出しており、理論において長年の欠落を是正している。グローバルな電荷中性を仮定せずに基礎方程式から再導出した結果、平均的には相殺するとみなされていたイオン価数 $z_i$ に比例する項が、正確な個々の種の熱力学を記述するために不可欠であることが示された。主な結果として、$\mu_i$ に対する完全で物理的に整合性のあるMSAの定式化が得られ、高濃度領域におけるモンテカルロシミュレーションとの一致が著しく向上した。

ABSTRACT

The Mean Spherical Approximation (MSA) is a commonly-used thermodynamic theory for computing the energetics of ions in the primitive model (i.e., charged hard-sphere ions in a background dielectric). For the excess chemical potential, however, the early MSA formulations (which were widely adopted) only included the terms needed to compute the mean excess chemical potential (or the mean activity coefficient). Other terms for the chemical potential $\mu_i$ of individual species $i$ were not included because they sum to $0$ in the mean chemical potential. Here, we derive these terms to give a complete MSA formulation of the chemical potential. The result is a simple additive term for $\mu_i$ that we show is a qualitative improvement over the previous MSA version. In addition, our derivation shows that the MSA's assumption of global charge neutrality is not strictly necessary, so that the MSA is also valid for systems close to neutrality.

研究の動機と目的

• この論文は、個々のイオン種の化学ポテンシャルに対するMSA定式化において、長年にわたり無視されてきた価数依存項（$z_i$）の欠落を扱う。
• 標準MSAにおける熱力学的不整合を是正することを目的とし、$z_i$ に比例する項が平均化学ポテンシャルでは相殺されるため無視されるが、それらは正確な個々の種の熱力学を記述するために不可欠である。
• 電気二重層のような不均一系の正確なモデル化を可能にするために、$\mu_i$ に対する完全で物理的に整合性のあるMSA定式化を提供することを目的とする。
• 本研究は、グローバルな電荷中性とは限らないが、それに近い系に対してもMSAが有効であることを示しており、その適用範囲を拡大する。

提案手法

• 著者らは、ブロムの元来の定式化からMSAを再導出し、以前は $+z_i \cdot \text{const}$ 項にまとめていた項を明確に分離した。
• 中心イオンの周囲では局所的電荷中性が成立することを示し、グローバルな電荷中性が成り立たない場合でも、不均一系へのMSAの適用が正当化されることを示した。
• 導出には、オーレンシュタイン＝ツェルニケ方程式、フーリエ変換、およびバクサー因数分解を用いて、直接相関関数 $c_{ij}(r)$ と総相関関数 $h_{ij}(r)$ を解いた。
• 重要なステップとして、$J_{ij}(r)$ を含む積分の再評価があり、以前は電荷中性を過剰に早期に仮定していたが、著者らはこの段階での $r^3$ 項が漸近的に小さいものの物理的に意味を持つことを示した。
• 背景電荷をイオンコアに埋め込むことで有効MSA系を定義し、有効系内で電荷中性を保ちながら内部エネルギーと化学ポテンシャルを導出した。
• 最終的に、個々の種の電気的過剰化学ポテンシャル $\mu_i$ の式は、$\beta\mu_i = \beta u_{0i} + \beta\delta\mu_{0i} - 2z_i u^*$ として得られ、ここで $u^*$ は系の構造因子と密度から導かれる項である。

実験結果

リサーチクエスチョン

• RQ1 標準MSA定式化における個々のイオン化学ポテンシャルが高濃度電解質で失敗する理由は何か？標準定式化には何が欠けているのか？
• RQ2 電荷中性が完全に成立しない、例えば局所的電荷不均衡を示す系に対してもMSAを適用可能か？
• RQ3 以前に $\mu_i$ から省かれていた $z_i$-依存項の物理的起源とその重要性は何か？
• RQ4 グローバルな電荷中性を仮定せず、第一原理から一貫性のある完全なMSA定式化 $\mu_i$ をどのように導出できるか？
• RQ5 これらの欠落していた $z_i$-依存項を含めることによる化学ポテンシャル予測の精度向上の定量的影響は何か？

主な発見

• 標準MSA定式化に欠落していた項はイオン価数 $z_i$ に比例しており、それらを含めることで個々のイオン化学ポテンシャルの記述が質的に改善される。
• 導出結果から、各イオンの周囲の局所的電荷中性が保たれる限り、グローバルな電荷中性とは限らない系に対してもMSAが有効であることが示された。
• 新しいMSA定式化における $\mu_i$ は $\beta\mu_i = \beta u_{0i} + \beta\delta\mu_{0i} - 2z_i u^*$ で与えられ、ここで $u^*$ は系の相関関数から導かれる構造依存項である。
• $-2z_i u^*$ 項の含めることで、特に標準MSAが失敗する高濃度電解質領域において、モンテカルロシミュレーションとの一致が著しく向上した。
• $z_i$-依存項は平均化学ポテンシャルでは相殺されるため、長年にわたり無視されてきたが、不均一系における正確な個々の種の熱力学には不可欠である。
• 導出により、標準MSAの内部エネルギーおよび自由エネルギー式は、有効価数 $z_i^\text{eff}$ を用いる限り有効であり、補正項 $\Delta E$ を適切に取り入れれば成立することが確認された。