QMMM interface related: ignore QM-MM bonds while forming bond lists, update charge matrix so that the total charge constraint only applies to QM atoms.
Currently, updated matrix works fine for non-periodic simulations but energy fluctuates for periodic simulations