分子动力学模拟蛋白质配体复合物时mdp步长求助

别说离开我

老师们好，在进行分子动力学模拟蛋白质配体复合物时，我的NVT热浴步长设置为2fs，但后续的NPT压浴和md的mdp文件设置为了1fs，这样合理吗，如果在论文中精确地表示说明出来，会不会不影响，出的各种RMSD,RMSF的图都很完美，还是说平衡时间长一点可以更稳定

别说离开我

student0618 发表于 2026-5-10 19:59
平衡时间长度看体系，不能一概而论。柔性高、体系大的，自然需时更长。

一般来说NVT非必要，直接NPT限制 ...

好的感谢老师解答，还是对这些掌握不透彻，我仔细看一下这篇博文

student0618

平衡时间长度看体系，不能一概而论。柔性高、体系大的，自然需时更长。

一般来说NVT非必要，直接NPT限制性动力学 + NPT 平衡相及产生相就可以了，见谈谈做分子动力学模拟什么时候适合用NPT和NVT系综
http://bbs.keinsci.com/forum.php ... 5&fromuid=64740
(出处: 计算化学公社)


步长一般是基于highest frequency vibration决定的，因此用了shake/lincs来constraint氢就可以增加到2fs，也因此如今全原子普遍用2fs步长。

用1fs通常是因为初始结构不理想先用较短步长平衡。生产相较少用，但也不是不行，也有少数特殊情况会用。就是很多时这样用其实有点浪费，一样的时间只模拟一半的模拟时长。