I am planning to model shock simulation for an energetic material, TKX-50. I want to do it by shock applying shock-absorbing boundary condition. According to the literature, to conduct this method, a rigid piston has to be created at one side of the simulation box and particles have to be moved toward the piston from another side. How can I prepare a rigid piston first in LAMMPS?