QM/MM methods for crystalline defects. Part 1: Locality of the tight binding model

Huajie Chen,Christoph Ortner

DOI: https://doi.org/10.48550/arXiv.1505.05541

2015-06-18

Abstract:The tight binding model is a minimal electronic structure model for molecular modelling and simulation. We show that the total energy in this model can be decomposed into site energies, that is, into contributions from each atomic site whose influence on their environment decays exponentially. This result lays the foundation for a rigorous analysis of QM/MM coupling schemes.

Numerical Analysis

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to explore the locality problem of the Tight Binding Model in crystal defect simulation and provide a theoretical basis for the QM/MM (Quantum Mechanics/Molecular Mechanics) coupling scheme. Specifically:

1. **Decomposition of total energy**:
The author shows that the total energy of the Tight Binding Model can be decomposed into the energy contributions of each atomic site, and the influence of each site on its surrounding environment decays exponentially with distance. This can be expressed as:
\[
E(y)=\sum_{\ell} E_{\ell}(y)
\]
where \(E_{\ell}(y)\) represents the energy of the \(\ell\) -th atomic site.

2. **Strong locality condition**:
In order for the QM/MM coupling method to be effective, a strong locality condition needs to be satisfied. That is, the force decays rapidly with distance:
\[
\left|\frac{\partial f_{n}}{\partial y_{m}}\right|\to 0\quad\text{when}\quad r_{nm}\to 0
\]
Here \(r_{nm} = |y_{n}-y_{m}|\), and \(f_{n}\) represents the force acting on the atom at position \(y_{n}\).

3. **Stronger energy locality**:
The author proves an energy locality condition stronger than the above - mentioned condition, that is:
\[
\left|\frac{\partial E_{\ell}(y)}{\partial y_{n}}\right|\lesssim e^{-\gamma r_{\ell n}}, \quad\left|\frac{\partial^{2} E_{\ell}(y)}{\partial y_{n}\partial y_{m}}\right|\lesssim e^{-\gamma(r_{\ell n}+r_{\ell m})}
\]
where \(\gamma> 0\), and there are similar results for higher - order derivatives.

4. **Applications and significance**:
This locality result has broad application prospects, for example:
- It provides a solid theoretical basis for the concept of inter - atomic potential.
- The analytical difference between the Tight Binding Model and the inter - atomic potential model is small, so many analytical tools developed for the inter - atomic potential model can be applied.
- Defect models in infinite crystalline solids can be constructed and analyzed rigorously.
- The construction and analysis of the extended atom/continuum multiscale scheme can be carried out, especially the Cauchy - Born continuum limit analysis can be directly applied to the Tight Binding Model.

5. **Main motivation**:
The ultimate goal is to provide new formulas and analysis methods for the QM/MM coupling scheme of crystal defects. By using the above - mentioned energy decomposition and locality conditions, the computational cost can be reduced while ensuring accuracy.

In summary, through rigorous mathematical analysis, this paper proves the strong locality of the Tight Binding Model, thus laying a theoretical foundation for the application of the QM/MM coupling method in material modeling.