GROMACS REST(replica exchange with solute tempering) 计算报错

紫藤竹

我现在在做一个蛋白质的REST(replica exchange withsolute tempering)计算重新定义type B   grompp时报错

Cannot automatically perturb a torsion withmultiple terms to different

form.

Please specify perturbed parametersmanually for this torsion in your

topology!

使用charmm27力场在ffbonded.itp中dihedraltypes有的相同原子类型的二面角会有不同参数

HS   S     CT2 CT1 9     0.00 1.00416   1

HS   S     CT2 CT1 9     0.00 0.6276     2

HS   S     CT2 CT1 9     0.00 1.12968   3

protein.itp中对于有多种二面角参数是怎么处理

36   43    45    47    9   0.00         6.6944     1

36   43    45    47    9   180.00     10.46      2

我直接在itp的dihedral添加多种参数这种错误消失了，不知道这种方法对不对

或者是这种

36   43    45    47    9   0.00         6.6944     1    180.00 10.46      2

同样报错