SharkYYX2025 发表于 2026-5-18 14:13
我做相关模拟也这样。这个是VMD判断原子成键的时候，判据把C-P键当成不连接了。实际上只是没显示。
如果你 ...

好的谢谢，果然md步长dt减少之后成功了，牛

我做相关模拟也这样。这个是VMD判断原子成键的时候，判据把C-P键当成不连接了。实际上只是没显示。
如果你的模拟中，观察到分子没有从这里断开，那成键关系没问题。如果断开了就是有问题。可以具体贴一下报错信息，方便判断。（这个地方没显示成键，和模拟崩溃的因果不明确）
Sobtop  http://sobereva.com/soft/Sobtop/      FAQ 8写了模拟崩溃的可能原因，可以挨个试试。当前建议md步长dt用0.001试试，我实测这种含有碳杂原子的分子，用2fs很容易崩溃

一个分子的话，在CPK等显示方式下，可以手动添加成键，在菜单选Mouse  Add/remove bonds，点击C P原子中心；一堆分子的话，不好手动添加，可以按照下图选择动态键。C-P cutoff一般需要调到1.9 2.0这样


如果模拟中DCB变成了两个片段，那可以检查itp。sobtop一开始，用mol2来生成拓扑文件，具有连接关系，这样没问题。

这是我的计算文件，InO.itp太大了发不了