你好，请问是哪一篇文献啊，能否分享下呢，十分感谢

sobereva 发表于 2024-4-25 06:28
sobtop主页http://sobereva.com/soft/Sobtop/#FAQ里明确说了

收到！感谢老师的回答，我再仔细看看去。

phosphorescence 发表于 2024-4-24 16:07
收到！感谢老师回复，还有两个问题想要请教老师。
问题一：想要使用packmol构建包含500个分子的模型，因 ...

sobtop主页http://sobereva.com/soft/Sobtop/#FAQ里明确说了

我强烈建议（这也是GROMACS的良好使用习惯），对体系中每种通过bond项（用于描述有显著共价特征的化学键，也包括配位键）连接到一起的片段，单独用Sobtop产生拓扑文件，并计算其原子电荷。比如离子液体，由于阴阳离子间是非共价作用，因此就应当对阳离子和阴离子分别照常产生拓扑文件，然后把二者的itp都include到主top文件里，并且把二者中的[ atomtypes ]也都合并、去重后挪到主top的[ defaults ]字段的后面去。

再顺带一提，有些零基础的GROMACS初学者，居然试图用Sobtop产生整个模拟体系（含一大堆分子）的拓扑文件，然后直接用GROMACS跑，这种做法是极端野蛮粗暴不可理喻的。这么做可能在Sobtop里耗时很高，甚至完全产生不了而崩溃，而且所有分子的拓扑信息都搅合在一起作为一个[ moleculetype ]出现时会特别难以管理和修改。

诸如[bonds]这种字段，都是对于它上面出现的[moleculetype]而言的，显然[moleculetype]里[atoms]中序号最多是多少，[bonds]里能出现的序号最多就是多少。

1 是
2 退火是10-13 ns期间
3 都属于。10ns的时候未必能达到平衡，所以没法说后10ns就已经能算作是产生相。但如果10ns的时候确实已经平衡了，那么后10ns可以当做产生相的前期。一般产生相的模拟时间尺度远长于平衡相模拟，所以后10ns不可能真正算作一般意义的完整的产生相。