各位大佬、社长好，
我希望对单层碳纳米管进行模拟，首先用VMD生成了结构文件
然后，按照sobtop的生成碳纳米管的实例教程，使用sobtop生成了itp、rtp文件，然后将rtp文件copy到对应力场目录下的
rtp、atp以及ffnonbonded等文件中，使用的是amber14sb_parmbsc1.ff力场。
用这个方式模拟小分子是没有问题的，但是模拟碳纳米管第一步pdb2gmx生成结构文件时，就出现了几百行内容一致的报错
Duplicate line found in or between hackblock and rtp entries
经过和同学交流后，发现上述报错可能源于sobtop生成的rtp文件中[bond]字段后碳原子没有被编号，无法被程序识别。
比如rtp文件中的[bond]字段内容如下（后面二面角同样的）
[ bonds ]
; name_i name_j r0 (nm) k (kJ/mol/nm^2)
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
C C 0.139840 3.858485E+05 ; Prebuilt ca-ca
可是sobtop生成的小分子的rtp，都对这些有很好的编号，也能模拟下来。
不知道是不是我的哪个操作失误了，导致碳管出现这样的问题...
#后来我用sobtop时，不再用pdb作为输入文件，而改用包含键连关系的mol2文件作为输入文件，并且手动指定原子类型，
发现这些键连的编号关系都填进去。
虽然这个报错解决了，但是接下来的模拟依然出现了Fatal error
Fatal error:
Atom C in residue CNT 1 was not found in rtp entry CNT with 816 atoms
while sorting atoms.
有点懵。。
继续检查了.gro文件后发现，依然是编号问题。。
故将pdb文件导入gaussview，去除H，保存为mol2文件，并输入sobtop再生成.gro文件后
此时终于解决了问题。。。