Thanks Aidan!
I was applying strain as “change_box all”. My initial thought was the repositioning of the atoms in the graphene sheet will come back to initial positions during the time integration between strain steps due to its finite nature. So the compressions may be due to the time integration is not long enough or the infinite dimension may have affected the repositioning.
However after I change the command to “change_box polymer”, I don’t see any compressions in the graphene sheet.
If I have a finite (in both directions) graphene sheet, is it still a good idea to use "change_box polymer” instead?
Best,
- A