锂电池固态电解质锂离子扩散 MD 模拟案例 | 分子动力学模拟实战

LLZO固态电解质中锂离子扩散的分子动力学模拟

这个项目在启动之前，已经有一个明确的结论需要验证。

实验团队测得Li₆.₅La₃Zr₁.₅Al₀.₅O₁₂(Al掺杂LLZO)在室温下的离子电导率约为8×10⁻⁴ S/cm，他们需要计算端给出对应的锂离子扩散系数，反推激活能，并分析Al掺杂对锂离子扩散路径的影响。这是一个典型的实验-计算协同研究场景。

麻烦从力场选择开始。

LLZO是一个多组元氧化物体系(Li-La-Zr-Al-O)，可用的文献力场并不统一。项目先后测试了三个来源的力场参数，其中一个来自2019年的一篇JPCB，另一个来自Materials Today Physics上针对石榴石型氧化物的专用参数集。

初始测试的标准是用力场跑出的晶格常数和弹性常数，与实验值的偏差控制在3%以内。前两个力场分别在La-O键长(偏差6.1%)和Li-O-Li弯曲势(偏差11.3%)上超标，最终选定的第三套参数在所有结构指标上通过验证。

这一步用了约一周时间。很多项目在这里没有做严格验证就直接开跑，结果出来的扩散系数往往偏低一个数量级，和实验对不上。

力场确定之后，温度窗口的设计是第二个关键决策。

室温下LLZO的锂离子扩散很慢，直接在300 K下跑MD，需要几十纳秒才能积累足够的统计样本，计算资源不允许。项目采用了”高温外推”策略：在600 K、700 K、800 K、900 K四个温度点各跑5 ns的NVT模拟，用Arrhenius关系外推到室温扩散系数。

“`

# GROMACS mdp关键参数

integrator    = md

dt            = 0.001       ; 1 fs

nsteps        = 5000000     ; 5 ns

tcoupl        = V-rescale

ref_t         = 700

tau_t         = 0.1

constraints   = none        ; 离子体系不加约束

“`

MSD(均方位移)计算时有一个细节：Li离子的MSD曲线在前0.5 ns有一段亚扩散行为(slope < 1)，这是固态体系中粒子在局域势阱中振荡的正常表现。扩散系数的拟合必须跳过这段，从MSD曲线的线性区间(0.5~4 ns)取斜率，否则D值会被系统性低估约20-30%。

四个温度点的扩散系数：

| 温度 (K) | D_Li (×10⁻¹⁰ cm²/s) |

|———-|———————-|

| 600      | 1.83                 |

| 700      | 4.76                 |

| 800      | 10.21                |

| 900      | 19.87                |

Arrhenius拟合得到激活能 Ea = 0.28 eV，外推到300 K的扩散系数为 2.1×10⁻¹² cm²/s，换算离子电导率约为 7.3×10⁻⁴ S/cm，与实验值(8×10⁻⁴ S/cm)偏差在10%以内。

Al掺杂的作用在轨迹分析中得到了直观体现：Al³⁺取代Zr⁴⁺位置后，局部形成锂空位，Li离子的跳跃路径由原来的单一通道扩展为多路径网络，这是掺杂提升电导率的微观机制。

这个数字层面的吻合，不是一开始就注定的。力场验证、温度窗口设计、MSD曲线的线性区间选取，每一步都可以把最终结果推向另一个方向。

如需了解类似固态电池体系的分子动力学模拟服务，欢迎访问科研学术网(https://www.keyanxueshu.com)查看更多案例。