我是用sobtop产生各种分子的itp,然后进行组合得到混合团簇的top。各种分子的itp文件都是采用sobtop主页(
http://sobereva.com/soft/Sobtop/)例2的方法产生的。事先都用Gaussian/B3LYP/6-31G*优化到稳定结构,均指认了GAFF原子类型,LJ参数和成键相关参数都是GAFF力场的,在屏幕上没有提示有任何bond、angle、dihedral参数缺失,也添加了RESP2(0.5)原子电荷。任何一种分子做成多分子团簇跑em都很快收敛(5步之内),但是任何两种分子混合团簇跑em都是3步即崩溃。