A hybrid neural network (NN) and Lennard-Jones (LJ) model driver for multilayer two-dimensional materials. The NN term models short-range intralayer and orbital overlap interactions and the theoretically-motivated LJ term models long-range dispersion.