这是本文档旧的修订版!
版本:2017.1
硫是用于可再充电能量存储设备中较为有前途的阴极材料。电池放电期间,锂离子从富锂阳极迁移到硫阴极,并与硫反应(阴极锂化):
$$x \textrm{Li} + \textrm{S} \rightarrow \textrm{Li}_x \textrm{S}$$
开路电压(OCV)是该锂化过程中电化学驱动力。
在本教程中,我们将利用分子动力学(MD)模拟来计算电池放电期间的 OCV 曲线。我们采用了 ATK-ForceField 计算器和在参考文献 [IOB+15] 中提到并由伊斯兰教同事开发的 ReaxFF 力场。模拟退火法可用于模拟一系列成分为 x 的无定形 LixS 化合物,并将所得 OCV 曲线与参考文献 [IOB+15] 中的图 3 进行比较。
参考文献 [IOB+15] 中的图 3 是由采用混合了大规模典型的蒙特卡洛/分子动力学(GC-MC / MD)方法获得的数据制作而成,该方法目前在 QuantumATK 中尚无法实现。我们没有使用 GC-MC 将 Li 原子置于低能位,而是在这里使用了一种模拟退火的方法。
忽略焓和熵对能量的贡献,OCV 根据总能量差计算得来:
$$\textrm{OCV} = -\frac{E_{\textrm{Li}_x \textrm{S}} - nE_\textrm{Li} - mE_\textrm{S}}{n} ,$$
在这里,$n$($m$) 是 LixS 中锂(硫)原子的数量(如$x=n/m$),$E_\textrm{Li}$ 和 $E_\textrm{S}$ 分别为每个纯锂和硫晶体的总能量。
我们首先考虑 Li0.4S,并描述为计算纯 S 部分锂化到 Li0.4S 的开路电压的所需步骤:
然后利用 ATK 的 Python 脚本自动模拟更多范围的 LixS 化合物。然后计算得到完整的 OCV 曲线,并研究不同成分的径向分布函数:
本教程提供了所有可供下载的 QuantumATK Python 脚本。您也可以在此处下载压缩文件:↓ scripts.zip。
这是一个高级教程,因此不会详细地介绍如何操作 ATK-ForceField 计算器和 QuantumATK 分析工具的相关基础知识。如有需要,请参阅教程 Molecular Dynamics: Basics。
在电池工作条件下,阴极材料不会如期待的保持完美的晶相;它最可能是无定形的。QuantumATK 提供了多种用于生成非晶结构的工具,例如 Amorhpous Prebuilder 和 PackMol 工具。本教程中我们使用后者。