Classical atomistic simulations of surfaces and heterogeneous interfaces with the charge-optimized many body (COMB) potentials | Litlas