[lammps-users] Why AIREBO for graphite not stable?

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:

What do papers using AIREBO for graphite say
is the correct cutoff?