I introduced dislocation in LAMMPS by removing 1/2 plane. but in FCC Aluminum, I want to create dislocation by dividing the box into two half crystal with N and N-1 plane and then minimizing it, according to the Osetsky method. I have no idea about how to create dislocation by this method, although it is mentioned in the paper that I attached here.
How do I model this with MD? thanks to help me.
The Osetsky method paper: https://doi.org/10.1088/0965-0393/11/4/302