计算化学公社

标题: 求助：GROMACS 计算SPE体系的锂阳离子 MSD 呈周期性震荡而非单调上升？ [打印本页]

作者
Author: masjikx289 时间: 2025-9-1 16:21
标题: 求助：GROMACS 计算SPE体系的锂阳离子 MSD 呈周期性震荡而非单调上升？

各位老师好！

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

(, 下载次数 Times of downloads: 0)

我的MSD曲线：

(, 下载次数 Times of downloads: 0)

原文献中的MSD曲线：

(, 下载次数 Times of downloads: 0)

已尝试的排查

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

2）itp文件中的力场设置

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)