Error when Using Command read_restart

Dear LAMMPS Users,

I am running a simulation using peri/lps model. In the input script, I used ‘restart 10000 test.restart’ command to generate restart files. Then I used the following script to restart a previous simulation:

read_restart test.restart.20000
neighbor 0.0010 bin
fix 1 all nve
fix 21 bottom setforce 0.0 0.0 0.0
fix 2 all indent 1e17 sphere 0.0000 v_y 0.0000 0.0015 units box
fix 3 all ave/time 1 1 100 f_2 f_2[1] f_2[2] f_2[3] file indent.restart
timestep 1.0e-7
run 10000

When I execute this script, I got an error message:

ERROR: Pair peri requires an atom map, see atom_modify (…/pair_peri_lps.cpp:431)

I copy & paste all the neighbor, fix, and timestep settings from the original script. Because the manual said those information is not stored in the restart file. I have no clue what causes this error. Any suggestions will be appreciated.

Best regards,
Yifei

Dear LAMMPS Users,

I am running a simulation using peri/lps model. In the input script, I used
'restart 10000 test.restart' command to generate restart files. Then I used
the following script to restart a previous simulation:

read_restart test.restart.20000
neighbor 0.0010 bin
fix 1 all nve
fix 21 bottom setforce 0.0 0.0 0.0
fix 2 all indent 1e17 sphere 0.0000 v_y 0.0000 0.0015 units box
fix 3 all ave/time 1 1 100 f_2 f_2[1] f_2[2] f_2[3] file indent.restart
timestep 1.0e-7
run 10000

When I execute this script, I got an error message:
ERROR: Pair peri requires an atom map, see atom_modify
(../pair_peri_lps.cpp:431)

I copy & paste all the neighbor, fix, and timestep settings from the
original script. Because the manual said those information is not stored in
the restart file. I have no clue what causes this error. Any suggestions
will be appreciated.

this is a very surprising statement, since LAMMPS already gave you a
suggestion right in the error message.
have a look at: http://lammps.sandia.gov/doc/atom_modify.html

axel.

Thank you very much for the suggestion. Now I add 'atom_modify map array’ to the top of the restart script. Then the previous error message is gone. But the simulation cannot continue. After the simulation stuck at step 20000 (when the restart file is generated) for a while, it gives an error message saying: Too many neighbor bins.

I am sure that I can run the simulation to at least 30000 steps without any problem. I wonder why I got ’too many neighbor bins’ error when I restart from 20000 step?

Best regards,
Yifei

Thank you very much for the suggestion. Now I add 'atom_modify map array’ to
the top of the restart script. Then the previous error message is gone. But
the simulation cannot continue. After the simulation stuck at step 20000
(when the restart file is generated) for a while, it gives an error message
saying: Too many neighbor bins.

please before posting questions about error messages, have a look at
this page (or the equivalent in the manual.pdf file corresponding to
the version of LAMMPS that you are using).

http://lammps.sandia.gov/doc/Section_errors.html#err_3

I am sure that I can run the simulation to at least 30000 steps without any
problem. I wonder why I got ’too many neighbor bins’ error when I restart
from 20000 step?

i don't have a crystal ball, so i don't know what is going on over
there. most likely your simulation is not running stable and generates
unphysical results.
i just made a test restarting a peri/lps run based on the input in
examples/peri/in.peri and this has no problem restarting.

axel.

Thanks. I will check the model carefully.

Have a nice day,
Yifei