[论文解读] 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）中单个离子种类的缺失化学势项，纠正了该理论中长期存在的遗漏。通过在不假设全局电中性的前提下重新推导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 中关于中心离子周围局部电中性的假设依然成立，从而验证了其在非均匀系统中的适用性。
• 推导过程使用了 Ornstein-Zernike 方程、傅里叶变换以及 Baxter 因子分解法，以求解直接关联函数 $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 此前被忽略的 $z_i$-依赖项的物理起源及其重要性是什么？
• RQ4 能否不依赖全局电中性假设，从第一性原理出发推导出一个一致且完整的 $\mu_i$ MSA 公式？
• RQ5 包含这些缺失的 $z_i$-依赖项对化学势预测精度的定量影响是什么？

主要发现

• 标准 MSA 公式中缺失的项与离子价数 $z_i$ 成比例，其引入可显著改善对单个离子化学势的定性描述。
• 推导表明，只要每个离子周围的局部电中性得以保持，MSA 在接近但不严格满足全局电中性的系统中依然有效。
• 新的 $\mu_i$ MSA 公式为 $\beta\mu_i = \beta u_{0i} + \beta\delta\mu_{0i} - 2z_i u^*$，其中 $u^*$ 是由系统关联函数导出的结构依赖项。
• 引入 $-2z_i u^*$ 项后，与蒙特卡罗模拟的吻合度显著提高，尤其在标准 MSA 失效的高电解质浓度区域。
• 此前忽略 $z_i$-依赖项的原因是其在平均化学势中相互抵消，但在非均匀系统中，这些项对准确描述单个物种热力学至关重要。
• 推导结果证实，当使用有效价数 $z_i^\text{eff}$ 且正确计入校正项 $\Delta E$ 时，标准 MSA 的内能和吉布斯自由能表达式依然有效。