When we run the molecule dynamic on a twin Fe crystals model with GULP module, the upper or lower periodic boundary may shift into the crystals model, even if we fixed the top and bottom atoms. Can Biovia explain why this may happen and teach us how to fix the upper and lower boundaries? By the way, the Materials Studio version is 2020.The tested version of ms2019 has no problem.