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