I’m a new user of LAMMPS. I believe I’ve used the correct input file to describe the 2 layer graphene structure. However, the result of energy minimization makes the graphene system shrink into 1 layer (as I see from the dump file). What is the best way to keep it 2 layers, or any command that I miss? I attached the files for you to take a look at. bilayer.lmp (32.7 KB) blg.in (7…