求助：Amber进行MMGBSA计算时CalcError:failed with prmtop

氯酸钡

创建complex、receptor和ligand命令：

1. ante-MMPBSA.py -p com1_2QRZ.prmtop --radii mbondi2 -s ":WAT,Na+,MG" -c com_2QRZ.prmtop -r receptor_2QRZ.prmtop -n ":GTP" -l GTP_2QRZ.prmtop

复制代码

运行命令：

1. MMPBSA.py -O -i mmgbsa_2QRZ.in -o test_RESULT.dat -cp com_test.prmtop -rp receptortest.prmtop -lp GTPtest.prmtop -y test_2QRZ.nc

复制代码

创建轨迹文件命令：
parm 2qrzWT_solvent.prmtop
trajin md3.nc 24980 last 1
trajin md4.nc 0 10 1
autoimage
strip :WAT,Na+,MG
trajout test_2QRZ.nc

报错：
File "/lustre/home/acct-phazyx/phazyx/Amber20/Amber/amber20_src/bin/MMPBSA.py.MPI", line 100, in <module>
    app.run_mmpbsa()
  File "/lustre/home/acct-phazyx/phazyx/Amber20/Amber/amber20_src/lib/python2.7/site-packages/MMPBSA_mods/main.py", line 218, in run_mmpbsa
    self.calc_list.run(rank, self.stdout)
  File "/lustre/home/acct-phazyx/phazyx/Amber20/Amber/amber20_src/lib/python2.7/site-packages/MMPBSA_mods/calculation.py", line 82, in run
    calc.run(rank, stdout=stdout, stderr=stderr)
  File "/lustre/home/acct-phazyx/phazyx/Amber20/Amber/amber20_src/lib/python2.7/site-packages/MMPBSA_mods/calculation.py", line 375, in run
    Calculation.run(self, rank, stdout=self.output)
  File "/lustre/home/acct-phazyx/phazyx/Amber20/Amber/amber20_src/lib/python2.7/site-packages/MMPBSA_mods/calculation.py", line 157, in run
    self.prmtop))
CalcError: /lustre/home/acct-phazyx/phazyx/Amber20/Amber/amber20_src/bin/cpptraj failed with prmtop com.prmtop!
Error occured on rank 54.
Exiting. All files have been retained.

版本是amber20，python2.7
网上搜到的大部分答案都无法解决，求赐教
文件已附（轨迹取了30帧）

zfzhao

老师，您好，请问您计算MMGBSA的这个错误有解决吗，他是什么原因导致的，解决办法是？？

Huschein

我记得会有详细的log还是out输出吧，那个里面的会记载更详细的失败原因，光看这个看不出来

Sagaf

我遇到了同一个报错，igb为2的时候正常计算，换别的就报错，最后是生成复合物的拓扑文件时指定了set default PBRadii mbondi2，然后再用这个到处的拓扑文件生成受体和配体的对应文件，然后计算就正常进行了。