Dear LAMMPS users and developers,
Recently I was trying to use LAMMPS to study graphite using AIREBO potential. I constructed the MD domain with two graphene layers inside. However, when the simulation runs for a while, the two layers tend to slide w.r.t. each other instead of staying about their equilibrium. There is only one tunable parameter in pair_style airebo: LJ cutoff, so I felt I might have used a wrong value for it, could you please help me with this? Thanks!
The input to lammps: