An analytical benchmark for MD codes: testing and correcting LAMMPS

In the attached paper, recently appeared on Computer Physics Communications, we have proposed an analytical benchmark and a simple consistent Mathematica program for graphene and carbon nanotubes, that may serve to test any molecular dynamics code implemented with REBO potentials. By exploiting the benchmark, we checked results produced by LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) when adopting the second generation Brenner potential, we made evident that this code in its current implementation produces results which are offset from those of the benchmark by a significant amount, and provided evidence of the reason.