计算化学公社

标题: MSD计算聚合物电解质体系锂离子扩散系数重复性差，平行实验计算的值差别很大 [打印本页]

作者
Author: shiyuanfan 时间: yesterday 09:38
标题: MSD计算聚合物电解质体系锂离子扩散系数重复性差，平行实验计算的值差别很大
本帖最后由 shiyuanfan 于 2026-1-24 15:41 编辑

各位老师，我的体系是PEO/LiTFSI（16：1），用的是GAFF2力场，sobtop生成的拓扑文件，锂离子采用的MERZ参数。平衡相模拟参数是
annealing_time = 0 1000 3000 4000 5000 6000 8000 9000 10000 11000 13000 14000 15000 20000;
psannealing_temp = 0 500 500 400 400 500 500 400 400 500 500 400 303.15 303.15
产生相进行了10ns的303.15K模拟，取2-7ns进行拟合
在通过MSD计算得到扩散系数，进行平行重复，发现每次重复的结果偏差都很大，3次得到的扩散系数分别为0.0002626 (+/- 0.0000) (1e-5 cm^2/s)，0.0004117 (+/- 0.0000) (1e-5 cm^2/s)，0.0003591 (+/- 0.0000) (1e-5 cm^2/s）差别很大
我想法请问这种平行实验的误差是正常的吗？通过优化什么模拟参数能够降低这种误差呢？通过多次模拟取平均值是否可行？
感谢各位老师！
以下是我的模拟文件

作者
Author: sobereva 时间: yesterday 15:00
确保平衡相到最后确实已经平衡了，并且把产生相增加到50ns
步长建议用2 fs，节约一倍耗时

没有结构文件，没法在结构方面进行评论

作者
Author: shiyuanfan 时间: yesterday 15:35
本帖最后由 shiyuanfan 于 2026-1-24 15:44 编辑

sobereva 发表于 2026-1-24 15:00
确保平衡相到最后确实已经平衡了，并且把产生相增加到50ns
步长建议用2 fs，节约一倍耗时

谢谢站长，我试一下。结构文件没注意上传失败了，我又添加了一个压缩包，里面有eq.gro和三次prod的结构文件。

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