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 nonlinear setting.

We consider atomistic geometry relaxation in the context of linear tight binding models for point defects. A limiting model as Fermi-temperature is sent to zero is formulated, and an exponential rate of convergence for the nuclei configuration is established. We also formulate the thermodynamic limit model at zero Fermi-temperature, extending the results of [H. Chen, J. Lu and C. Ortner, Thermodynamic limit of crystal defects with finite temperature tight binding, Arch. Ration. Mech. Anal. 230 (2018) 701–733]. We discuss the non-trivial relationship between taking zero temperature and thermodynamic limits in the finite Fermi-temperature models.

The tight binding model is a minimalistic electronic structure model for predicting properties of materials and molecules. For insulators at zero Fermi-temperature we show that the potential energy surface of this model can be decomposed into exponentially localised site energy contributions, thus providing qualitatively sharp estimates on the interatomic interaction range which justifies a range of multi-scale models. For insulators at finite Fermitemperature we obtain locality estimates that are uniform in the zero-temperature limit. A particular feature of all our results is that they depend only weakly on the point spectrum. Numerical tests confirm our analytical results.This work extends and strengthens (Chen, Ortner 2016) and (Chen, Lu, Ortner 2018) for finite temperature models.

We show that the local density of states (LDOS) of a wide class of tight-binding models has a weak body-order expansion. Specifically, we prove that the resulting body-order expansion for analytic observables such as the electron density or the energy has an exponential rate of convergence both at finite Fermi-temperature as well as for insulators at zero Fermi-temperature. We discuss potential consequences of this observation for modelling the potential energy landscape, as well as for solving the electronic structure problem.