50 ns后RMSD看似稳定但结构被“拉丝”，mdp参数该怎么排查？

模小助

这里有个 mdp 参数细节，想看看大家通常怎么设。

我做了一个50ns的模拟，看rmsd应该是稳定了，但是抽了几帧结构发现跑拉丝了。Run parametersintegrator = md ;

主要想确认：50 ns后RMSD看似稳定但结构被“拉丝”，mdp参数该怎么排查？

admin

“RMSD稳定但结构被拉丝”通常首先怀疑 PBC 可视化问题，而不是 RMSD 本身。建议先对轨迹做 -pbc mol/-center 处理，再看结构是否仍然异常。

如果处理 PBC 后依旧有长键或原子被真正拉开，再检查约束、时间步长、LINCS warning、温压耦合和初始接触。RMSD 只反映选定原子经过拟合后的整体偏差，不能证明每个局部结构都合理，所以最好同时看键长异常、能量、Rg 和可视化轨迹。