Dear LAMMPS users,
I am having issues restarting rRESPA calculations.
LAMMPS version: LAMMPS (29 Sep 2021 - Update 1)
While trying to figure out the appropriate parameters for respa in NVE runs (no barostats, thermostats, fixed regions etc.), I encountered the following problem. When trying to restart from a restart file written during a respa run, LAMMPS reads fewer atoms than are present in the system and throws the error message:
ERROR: Did not assign all restart atoms correctly (…/read_restart.cpp:471)
Observations and additional tests:
The error happens even if one uses a restart file generated after just 1 integration step or using a small time step (1 fs). Interestingly, the number of atoms read from the restart in these two cases is also different. The error does not occur, if 1-level respa is used.
The classical explanation I found is that atom coordinates may be outside the bounds of the non-periodic box, actually, my box is periodic. Using the “remap” flag fixes the issue, though, suggesting atom coordinates may be indeed outside the bounds of the box.
I am wondering, though, why the error happens only in multilevel respa, even if the time step is half as large as in the one-level case.
After using remap, all energies in the first energy output of the new run a are identical to the last step from the previous run, as it should be.
However, the pressure is not.
Did anyone experience the same kind of behavior and what are your thoughts on that?
Is it intended that one has to use “remap” if one restarts respa runs?
Why would the pressure be different, or is it expected after a remap?
If the level of detail is not sufficient, I will provide all the inputs and descriptions to reproduce the behavior