求助：GROMACS 计算SPE体系的 锂阳离子 MSD 呈周期性震荡而非单调上升？

masjikx289

各位老师好！

我正在复现: Phys. Chem. Chem. Phys.,2023, 25, 29894这篇文献 中的 LiEtTFSI-G99 体系（O/Li = 20，80 °C）的 Li⁺ 电导率计算。体系、力场、盒子尺寸、离子数均按原文设定，但在 NVT 恒温模拟后 对 MSD 采样时，Li⁺ 的 MSD 曲线呈明显周期性震荡，而原文图 为平滑单调上升。显然聚阴离子和Li+之间的耦合非常强，没有游离的Li+产生，这显然与原文中的Li+解离度为67.5%是矛盾的。

我的MSD曲线：

原文献中的MSD曲线：

已尝试的排查

1）原子间作用力参数设置

2）itp文件中的力场设置

3）更换不同的ns(correction factors)

请问还有哪些原因导致我的体系Li+与聚阴离子的耦合非常强，无法解离，以致于无法生成合理的MSD曲线？附录中有1）原文献及其SI，2）我的nvt gromacs模拟算例