[AMBER] problem for MMPBSA.py in AmberTools12

From: zhenquan hu <zhenquanhu.yahoo.com.cn>
Date: Fri, 3 Aug 2012 21:49:15 +0800 (CST)

Dear Amber Users,

I found a problem on running mmpbsa job in Amber12/AmberTools12


MMPBSA.py -O -i mmpbsa.in -o f.dat -sp com_ions.top -cp com.top -rp rec.top -lp lig.top -y md.trj
Reading command-line arguments and input files...
Loading and checking parameter files for compatibility...
mmpbsa_py_energy found! Using /opt/amber12/bin/mmpbsa_py_energy
cpptraj found! Using /opt/amber12/bin/cpptraj
Preparing trajectories for simulation...
4 frames were processed by cpptraj for use in calculation.

Beginning PB calculations with /opt/amber12/bin/mmpbsa_py_energy
  calculating complex contribution...
CalcError: /opt/amber12/bin/mmpbsa_py_energy failed with prmtop com.top!
Exiting. All files have been retained.


I checked the program in $AMBERHOME/AmberTools/test by
./test_at_serial.sh and make test.mmpbsa
and all passed

Can anyone give me some suggestions?


Below are my parameters for the MM-PB/SA job
#mmpbsa.in   =========================
PB test
&general
   startframe=1, endframe=4
/
&pb
   istrng=0.150,
/
#===================================

Topology file generated:
# leap.in      ===========================
source leaprc.ff12SB
source leaprc.gaff
loadamberparams frcmod.ionsjc_tip3p
loadamberparams lig.frcmod
lig = loadMol2 lig.mol2
rec = loadpdb rec.pdb
com = combine {rec lig}
set default PBRadii mbondi2
saveamberparm com com.top com.crd
saveamberparm rec rec.top rec.crd
saveamberparm lig lig.top lig.crd
solvateoct com TIP3PBOX 10
addions com Na+ 1 Cl- 1
addions com Na+ 1 Cl- 1
...
addions com Na+ 1 Cl- 1
saveamberparm com com_ions.top com_ions.crd
quit
#=====================================

MD simulation with pmemd.MPI
#  md.in       =============================
Production
 &cntrl
   imin = 0, irest = 1, ntx = 5,
   ntb = 2, ntp = 1, taup = 0.2
   es_cutoff = 8.0, vdw_cutoff = 9.0,
   ntc = 2, ntf = 2, tempi = 300.0, temp0 = 300.0
   ntt = 3, gamma_ln = 2.0, nstlim = 500000, dt = 0.002 ioutfm = 1,
   ntpr = 100, ntwx = 1000, ntwr = 1000
   ntr=0   ig = -1
#=====================================
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Aug 03 2012 - 07:00:04 PDT
Custom Search