Hi dear LAMMPS users,
I wan to model CSL (Coincidence Site lattice) grain boundaries, here for example sigma 5 (310) with θ=36.87 tilt angle for Aluminum. I create two crystals and rotate each for θ/2 and I delete the overlapped atoms, and after visualising with VMD all the time I have two atoms in the base of kite-shaped structure that normally are supposed to be one to have my desired structure but instead I have two atoms very close to each other connected with a bond.
Can anybody help please?
here is my script: