Increasing simulation box makes bond error

Hello everyone,

After equilibration, I was trying to merge the following simulation box with another box with bigger dimension. But when I increase the dimension of this box by manualy changing the xyz values in the data file, the bond shows strange picture. I know this is because some atoms are on the periodic wall.

How do I avoid this strange results?

normal
error

Regards
Nur

Please see this recent discussion for some strategies that you may be able to adapt for your case.

Another possible procedure:

  • use write_dump to write out unwrapped coordinates
  • use set command to set all image flags to 0
  • use change_box to change your box dimensions
  • use read_dump with box no to read back the coordinates, but retain the enlarged box
  • use reset_atoms image to reset image flags for the dump data while keeping molecules intact
  • use write_data to write out the updated system as data file
1 Like