| 查看: 788 | 回复: 1 | |||
ra2ghgzh木虫 (正式写手)
学痴
|
[求助]
Gromacs中最大distance violation超过1nm，而Dis. Rest.为零该怎么办？
|
|
请教各位专家。 我想在使用Gromacs时限制原子间距在0.1-0.3nm之间，所以在top文件中写入： [ distance_restraints ] ; ai aj type index type’ low up1 up2 fac 5 28 1 0 1 0.1 0.3 2.0 10000.0 6 29 1 1 1 0.1 0.3 2.0 10000.0 然后在mdp文件中写入： disre = Simple disre-weighting = Conservative disre-mixed = no disre-fc = 1000 disre-tau = 0 nstdisreout = 100 运行500ps后，使用gmx disre得到了drmax.xvg。结果显示，很多时候间距都超过了限制。再检查log文件，却发现Dis. Rest.一直为零，也就是说虽然违反了距离限制，但是在没有在能量上反映出来。但是总体来说，加上了distance_restraints后还是有一定效果的，原子间距的确比没加时小了很多。我想知道是不是哪里的设置出了问题，应该如何修改。 谢谢各位。 |
回复此楼» 猜你喜欢
材料工程专硕求调剂
已经有5人回复
-
296求调剂
已经有10人回复
-
332求调剂
已经有5人回复
-
318求调剂
已经有7人回复
-
求调剂 有机化学考研356分
已经有10人回复
-
279求调剂
已经有5人回复
-
283求调剂（080500）
已经有12人回复
-
考研调剂
已经有4人回复
-
300求调剂，材料科学英一数二
已经有14人回复
-
化学308分调剂
已经有6人回复
2楼2018-11-05 07:24:37
10