Numerical Analysis of the Multiple Scattering Theory for Electronic Structure Calculations

Xiaoxu Li,Huajie Chen,Xingyu Gao

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

2021-11-14

Abstract:The multiple scattering theory (MST) is one of the most widely used methods in electronic structure calculations. It features a perfect separation between the atomic configurations and site potentials, and hence provides an efficient way to simulate defected and disordered systems. This work studies the MST methods from a numerical point of view and shows the convergence with respect to the truncation of the angular momentum summations, which is a fundamental approximation parameter for all MST methods. We provide both rigorous analysis and numerical experiments to illustrate the efficiency of the MST methods within the angular momentum representations.

Numerical Analysis,Computational Physics

What problem does this paper attempt to address?

This paper attempts to solve the problems of convergence and error estimation regarding angular momentum truncation in the multiple - scattering theory (MST) for electronic structure calculations. Specifically, the author studies the MST method from the perspective of numerical analysis and shows the influence of angular momentum summation truncation on the convergence of the MST method. The following are the core problems and objectives of the paper:

### Problems the paper attempts to solve

1. **Numerical analysis of the MST method**:
- The MST method is a widely used method for electronic structure calculations, especially suitable for simulating defective and disordered systems. However, there has been no strict mathematical analysis for the numerical analysis of the MST method, especially regarding the convergence and error estimation of angular momentum truncation.

2. **Convergence of angular momentum truncation**:
- In the MST method, the accuracy of the solution depends on the truncation of angular momentum summation. Therefore, it is crucial to study the convergence of angular momentum truncation on the solution, which determines the reliability and efficiency of MST simulation results.

3. **Error estimation**:
- The author aims to derive a priori error estimates of the MST solution (i.e., the solution of the KKR eigen - equation) with respect to angular momentum truncation. Specifically, they hope to prove the super - algebraic convergence rate of eigenpairs under certain assumptions.

### Main contributions

- **Strict mathematical analysis**: The author provides a strict mathematical analysis, filling the gap in the numerical analysis of the MST method.
- **Super - algebraic convergence rate**: It is proved that under certain assumptions, the eigenvalues and eigenfunctions of the MST solution have a super - algebraic convergence rate with respect to angular momentum truncation.
- **Numerical experiment verification**: The results of theoretical analysis are verified by numerical experiments, further supporting their conclusions.

### Formula summary

The key formulas involved in the paper include:

- Super - algebraic convergence rate of the KKR eigen - equation:
\[
|E_L - E|+\|u_L - u\| \leq C_s L^{1 - s}\quad\forall s > 1
\]
where \(E\) and \(u\) are energy (eigenvalue) and wave function (eigenfunction), \(L\) is the angular momentum truncation, and the constant \(C_s\) depends only on the system and \(s\).

- Definition of the t - matrix:
\[
t_\ell(z):=\frac{b_\ell(z)}{a_\ell(z)}\quad\text{for}\quad\ell\in\mathbb{N}
\]

- Weak continuity condition:
\[
\int_{\Gamma_n}[v]_{\Gamma_n}Y_{\ell m}(\hat{r})d\hat{r} = 0\quad\forall 0\leq\ell\leq L,|m|\leq\ell
\]

Through these formulas and analysis, the author provides a solid theoretical basis for the numerical reliability of the MST method and guidance for parameter selection in practical applications.