GROMACS version: 2019.4
GROMACS modification: Yes with plumed
I have performed an MD run. The simulation has run smoothly without errors or warnings about inconsistent shifts (though it is not finished yet, as it stopped due to maxh option).
When I try to generate a trajectory with:
gmx_mpi trjconv -f prod_5/md.xtc -s prod_5/md.tpr -pbc whole -o prod5_prot.xtc
I get the following error:
Program: gmx trjconv, version 2019.4 Source file: src/gromacs/pbcutil/mshift.cpp (line 917) Fatal error: There are inconsistent shifts over periodic boundaries in a molecule type consisting of 543 atoms. The longest distance involved in such interactions is 4.286 nm which is above half the box length. This molecule type consists of muliple parts, e.g. monomers, that are connected by interactions that are not chemical bonds, e.g. restraints. Such systems can not be treated. The only solution is increasing the box size.
Does that mean there was a problem with the trajectory? Is there a way to solve it or I should start a new simulation from scratch. If so, what can be the cause of this error? I have run many similar simulations with other peptides without encountering this error.
Thanks for your time.