[AMBER] error running mmpbsa

From: giulia palermo <giulia.palermo83.gmail.com>
Date: Fri, 23 Jun 2017 13:32:54 -0700

Dear all,

I am doing some MMPBSA calculations using the dollowing input and command
line:

--- mmpbsa.inp
&general
   endframe=2201, keep_files=2, interval=6,
/
&gb
  igb=2, saltcon=0.100,
/
&pb
  istrng=0.100,
/

--- command line
MMPBSA.py -O -i mmpbsa.in -o FINAL_RESULTS_MMPBSA.dat -cp
./0_make_tutto/tutto.top -rp ./0_make_azzurro/azzurro.top -lp
./0_make_rosa/rosa.top -y tutto.mdcrd

However, the calculation crashes and I am having the following error:

Loading and checking parameter files for compatibility...
mmpbsa_py_energy found! Using
/net/software/pkg/amber/amber14-test/bin/mmpbsa_py_energy
cpptraj found! Using /net/software/pkg/amber/amber14-test/bin/cpptraj
Preparing trajectories for simulation...
Error: Third line contains invalid atom index.

*Error: -17.043Error: Could not set up tutto.mdcrd for reading.Error:
Could not set up input trajectory 'tutto.mdcrd'.*
  File "/net/software/pkg/amber/amber14-test/bin/MMPBSA.py", line 103, in
<module>
    app.file_setup()
  File
"/net/software/pkg/amber/amber14-test/lib/python2.7/site-packages/MMPBSA_mods/main.py",
line 156, in file_setup
    self.mpi_size, str(external_progs['cpptraj']), self.pre)
  File
"/net/software/pkg/amber/amber14-test/lib/python2.7/site-packages/MMPBSA_mods/make_trajs.py",
line 66, in make_trajectories
    traj = Trajectory(FILES.complex_prmtop, FILES.mdcrd, cpptraj)
  File
"/net/software/pkg/amber/amber14-test/lib/python2.7/site-packages/MMPBSA_mods/make_trajs.py",
line 467, in __init__
    self.Query()
  File
"/net/software/pkg/amber/amber14-test/lib/python2.7/site-packages/MMPBSA_mods/make_trajs.py",
line 602, in Query
    raise TrajError('%s failed when querying %s' % (self.exe, traj))
TrajError: /net/software/pkg/amber/amber14-test/bin/cpptraj failed when
querying tutto.mdcrd
Exiting. All files have been retained.


I have successfully preformed MMPBSA calculations in the past and I can not
understand what is gouing wrong.
It seems that it can not read the trajectory of the complex. But I have
checked it and it seems OK.
Do you have any idea?

thank you

Giulia
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Jun 23 2017 - 14:00:02 PDT
Custom Search