GROMACS分子动力学模拟：一个离子液体体系中锂离子传输的机理研究

离子液体作为新型电解质的潜力已经不需要再论证，但锂离子在离子液体中的传输行为比传统碳酸酯电解质复杂得多——高黏度、强配位、非经典输运特征，这些因素让实验手段很难给出微观层面的解释。GROMACS分子动力学模拟在这里提供了一个直接观察离子运动的窗口。

项目背景：高浓度锂盐的传输悖论

团队处理过一个1-乙基-3-甲基咪唑鎓双(三氟甲基磺酰)亚胺盐（EMIM-TFSI）+ 双(三氟甲基磺酰)亚胺锂（LiTFSI）体系的锂离子输运研究项目。客户正在开发基于离子液体的固态电池电解质，需要知道锂盐浓度对锂离子电导率的影响规律。

实验数据显示，在LiTFSI摩尔分数从0.05增加到0.20时，锂离子电导率先升后降——低浓度时增加锂盐会提高载流子浓度，但高浓度时黏度急剧增大反而降低了离子迁移率。这个”先升后降”的转折点在哪里，实验需要大量样品反复测试。分子动力学模拟可以用一组参数扫描直接给出不同浓度下的离子扩散系数，辅助实验快速定位最优浓度。

力场参数与体系构建

EMIM阳离子和TFSI阴离子的力场参数基于OPLS-AA力场框架，Lennard-Jones参数和部分电荷取自文献中经过量子化学计算验证的参数集。Li⁺的力场使用Joung-Cheatham参数——这是专门针对水溶液和有机溶剂中Li⁺优化过的参数，对第一溶剂化壳层的描述精度优于通用的OPLS参数。

团队构建了四个不同LiTFSI摩尔分数的模拟体系：x=0.05、0.10、0.15、0.20，每个体系约包含200个EMIM-TFSI离子对和相应数量的LiTFSI。使用GROMACS 2023版本，在NPT系综下（313 K，1 bar）各运行100 ns的平衡模拟。

均方位移分析揭示非经典输运

锂离子的均方位移（MSD）分析给出了令人意外的结果。在低浓度（x=0.05）下，锂离子的MSD曲线呈线性增长趋势，扩散系数为2.1×10⁻⁷ cm²/s，符合经典Fick扩散行为。但随着浓度增加到x=0.15和x=0.20，MSD曲线出现了明显的亚扩散行为——在长时间尺度上MSD的增长偏离线性，斜率逐渐减小。

这意味着高浓度下锂离子的运动不再是自由的随机行走，而是受到了周围离子网络的束缚。锂离子的配位分析证实了这一点：在x=0.05时，每个Li⁺平均配位4.2个TFSI阴离子中的O原子，配位交换频率约0.8 ns⁻¹；在x=0.20时，配位数增加到5.8个，配位交换频率下降到0.15 ns⁻¹——锂离子被更紧密地”锁”在了配位网络中。

配位环境的可视化分析

为了更直观地展示锂离子的配位环境，团队做了配位数的径向分布函数（RDF）分析。Li⁺-O(TFSI)的RDF第一峰位于2.05 Å，与文献中X射线衍射数据一致。随着浓度升高，第一峰的高度增加、宽度减小，说明配位更加紧密但更加刚性。

连续配位寿命（continuous residence time）的分析给出了另一个关键信息。在x=0.05时，Li⁺-TFSI配位的平均连续寿命为0.8 ns；在x=0.20时增长到4.5 ns。配位寿命越长，锂离子挣脱一个配位环境、跳跃到下一个”宿位”的时间间隔就越大，宏观上表现为扩散系数的下降。

四个浓度的锂离子扩散系数汇总：x=0.05时2.1×10⁻⁷、x=0.10时1.4×10⁻⁷、x=0.15时0.6×10⁻⁷、x=0.20时0.2×10⁻⁷ cm²/s。这个趋势与实验电导率的变化方向完全一致——模拟预测的最优浓度区间在x=0.08-0.12，客户据此安排了针对性的实验验证，最终确认x=0.10为最优配比。

GROMACS在这个项目中的计算效率很高——每个100 ns模拟在单节点（32核）上只需约12小时，四个浓度总共不到3天。作为一款开源的分子动力学软件，GROMACS在离子液体、电解质、聚合物溶液等软物质体系中都有着成熟的应用案例，参数可复现性也很好。

如需了解更多的科研计算案例与项目经验，可以参考 https://www.keyanxueshu.com/ 上的技术文章。