面向电池析出反应的颗粒追踪方法

《The Journal of Physical Chemistry C》：Particle Tracking Methods for Battery Precipitation Reactions

　　

固-液界面处的析出与沉积反应在多种电池化学体系中发挥关键作用，涵盖锂离子电池、无负极电池、锌基电池及锂硫电池等体系。尽管已有文献建立了包含异相成核与生长（heterogeneous nucleation and growth, HNG）现象的模型，但现有研究尚未充分阐明用于追踪电极表面沉积物颗粒尺寸分布（particle size distribution, PSD）随时间演化的数值算法细节。研究人员评估了多种离散化与追踪PSD的方法，发现常用方法会导致PSD出现非物理性的展宽与平坦化现象。最终，研究人员提出了一种可在颗粒生长过程中保持PSD形态合理性的算法。

《The Journal of Physical Chemistry C》刊发的这项研究聚焦于电池异相成核与生长（HNG）模拟中颗粒尺寸分布（PSD）的数值追踪难题。固-液界面析出反应广泛存在于无负极锂金属、锂硫、锂空气等前沿电池体系中，其动力学行为直接决定电极形貌、阴极利用率及器件循环寿命，是预测电池性能与退化机制的核心要素。然而，现有文献虽普遍集成HNG动力学，却缺乏对PSD追踪算法的系统性描述与验证，导致模型假设与实际计算结果常出现偏差。

研究人员构建了简化的半球形颗粒生长模型，设定恒定径向生长速率（dr/dt = 0.25 μm·s^-1）与固定成核半径（r_nuc= 0.5 nm），确保理想PSD形状应保持不变，以此作为算法精度的基准。基于此框架，团队评估了三种数值算法：简单有限差分法、边界书签法及创新的“移动料斗”模型。

关键技术方法包括：1. 采用离散化粒径区间（bin）划分PSD，建立各区间颗粒数密度（N_i）的常微分方程组；2. 设计恒定与变形成核率两种工况，结合非维度时间（t?）分析PSD演化；3. 引入最大/最小半径约束与粒子均匀性假设检验，量化算法误差；4. 将算法植入锂空气电池模型进行实证验证，对比放电容量与计算效率。

研究结果如下：

3.1 算法1：简单有限差分法

该方法假定每个粒径区间内颗粒均匀分布。研究发现，由于区间内部同质性假设，颗粒在未达到区间上限前即被允许“扩散”至下一区间，导致PSD随时间异常展宽和平坦化。模拟结果显示颗粒最大半径远超理论值（1 μm），且尾部颗粒滞留时间过长，违背了恒定生长速率的物理前提。

3.2 算法2：“书签”法

该方法引入前导与后沿“书签”追踪极值粒径，成功限制了PSD的边界扩展。然而，算法仍继承区间均匀性假设，在变形成核率工况下，PSD在边界内依然发生不期望的平坦化。此外，书签跨越区间边界时引发残差突变，导致计算步长缩减，运行时间显著增加（达0.961秒）。

3.3 算法3：“移动料斗”模型

该模型颠覆传统思路，不再移动颗粒，而是移动“成核窗口”。成核仅发生在动态移动的特定区间，颗粒一旦沉积即固定于所在区间，仅通过增大该区间的代表性粒径（r_i）来反映生长。此方法彻底消除了区间内的伪扩散，无论成核率如何变化，PSD形状均保持恒定且与理论预测一致，计算耗时仅比算法1增加76%（0.338秒）。

3.4 对电池模拟的影响

在锂空气电池模型实证中，算法1与2的平均误差高达66.49%，而算法3精准复现了参考文献的PSD与放电容量（2.619 mAh）。研究表明，PSD的形态误差在循环模拟中会被放大，影响表面覆盖率计算及后续电化学性能预测。

讨论与结论部分指出，传统的有限差分法因强制颗粒在区间内均匀分布，必然引入非物理的PSD展宽。提出的“移动料斗”模型通过将生长转化为区间特征半径的更新，仅需新增一个关于最大半径（r₀）的微分方程，即以极低的计算成本解决了这一长期存在的数值稳定性问题。该算法具有普适性，可扩展至粒径依赖的生长速率场景。研究人员强调，准确的PSD追踪是确保电池模型可靠性的关键，建议在相关研究中明确记录所采用的数值算法。相关代码已在GitHub开源，以支持开放科学实践。