GROMACS version: 2020.3
GROMACS modification: No
Here post your question
Dear all,
When I run the script gmxpbsa0.sh of the GMXPBSA2.1, I already get the following error:
Fatal error:
Index[5885] 5887 is larger than the number of atoms in the
trajectory file (5886)
whatever I use the trajectory containing the protein-ligand complex and the water and that only containing the protein -ligand complex. Do you have any idea how this error is caused?
Hi,
in general such a error occurs, when your select a index group that, as dimension (number of atoms), does not correspond to the number of atoms in your trajectory.
Using an index file generated from a trajectory frame usually solves the problem. To generate an index file use, gmx make_ndx.
Thank you for you response. I use the md.gro which was generated simutaneously with the trajectory file MD.xtc in the MD simulation.
gmx make_ndx -f md.gro
This index file must have the same dimension with the trajectory, doesn’t it?
But the same error mentioned above comes out. So there might be other reason for it?
Hi,
Please note that not necessary the gmx mdrun output files, gro and and xtc , contain the same number of atoms.
What is saved in xtc file, depends on mdp option compressed-x-grps (group(s) to write to the compressed trajectory file, by default the whole system is written).
I have no experience in using the third-party tool GMXPBSA2.1, maybe someone else in the list has or you can directly contact the authors of the script.
even i’m getting same error while everything in my INPUT.dat file seems perfect. I ran the example mmpbsa and the run went smooth with no errors but when i introduce my trajectory files i get the error
See, i have set the c path, g path, a path and also the root directory accordingly and the example runs went fine. I have made a folder where INPUT.dat is present and its name is present in the INPUT.dat, yet i get the error saying the directory is not present. i don’t know where to ask about this since there is no forums made specially for gmxpbsa2.1. Thanks for the reply