  A key starting assumption in many classical interatomic potential models for
materials is a site energy decomposition of the potential energy surface into
contributions that only depend on a small neighbourhood. Under a natural
stability condition, we construct such a spatial decomposition for
self-consistent tight binding models, extending recent results for linear tight
binding models to the non-linear setting.
