Approximating Maximally Localized Wannier Functions with Position Scaling-Eigenfunction

Yuji Hamai,Katsunori Wakabayashi

DOI: https://doi.org/10.1103/PhysRevB.108.245413

2023-12-15

Abstract:Position scaling-eigenfunctions are generated by transforming compactly supported orthonormal scaling functions and utilized for faster alternatives to maximally localized Wannier functions (MLWFs). The position scaling-eigenfunctions are first applied to numerical procedures solving Schrödinger and Maxwell's equations, and the solutions well agree with preceding results. Subsequently, by projecting the position scaling-eigenfunctions onto the space spanned by the Bloch functions, approximated MLWFs are obtained. They show good agreements with preceding results using MLWFs. In addition, analytical explanations of the agreements and an estimate of the error associated with the approximation are provided.

Mesoscale and Nanoscale Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to approximate maximally localized Wannier functions (MLWFs) by using position - scaling - eigenfunctions. Specifically, the author proposes a method to generate position - scaling - eigenfunctions by transforming compactly supported orthonormal scaling functions and applies them to numerically solve the Schrödinger equation and the Maxwell equation. These position - scaling - eigenfunctions are then projected onto the space spanned by Bloch functions to obtain approximate MLWFs. The paper shows that these approximate MLWFs are in good agreement with the results in previous studies, and provides an analytical explanation for this agreement as well as error estimates related to the approximation.

### Main problems and solutions
1. **Generate position - scaling - eigenfunctions**:
- The author generates position - scaling - eigenfunctions by transforming compactly supported orthonormal scaling functions. These eigenfunctions have strong locality and are at least twice - differentiable, making them ideal basis functions.
- Formula representation:
\[
|x_p\rangle = \sum_n |n\rangle \langle n | x_p \rangle = \sum_n \xi_{n,p} |n\rangle
\]
where \(\xi_{n,p}=\langle n | x_p \rangle\).

2. **Calculate matrix elements of operators**:
- The paper elaborates on how to calculate the matrix elements of the momentum operator and the kinetic energy operator. These matrix elements are crucial for solving the Schrödinger equation and the Maxwell equation.
- The matrix element \(Q_r\) of the kinetic energy operator can be calculated by the following formula:
\[
Q_r=\langle x_r | \hat{p}^2/2 | x_0 \rangle = \sum_{k,l} \xi_{k - r} \xi_l P_{k,l}=\sum_{k,q} \xi_{k - r} \xi_{k - q} P_q
\]

3. **Approximate MLWFs**:
- By projecting the position - scaling - eigenfunctions onto the composite band space, the author obtains approximate MLWFs. These approximate MLWFs are very close to the actual MLWFs in shape and position.
- Approximation formula:
\[
|W_s^l\rangle \approx \frac{1}{\sqrt{\langle x_C^p | x_C^p \rangle}} |x_C^p\rangle
\]
where \(|x_C^p\rangle\) is the position - scaling - eigenfunction projected onto the composite band space.

4. **Error analysis**:
- The paper also provides error estimates between the approximate MLWFs and the actual MLWFs, which is very important for evaluating the effectiveness of the approximation method.
- Error formula:
\[
\langle x | x_C^p \rangle=\langle x_p | \tilde{x}_C^p \rangle+O\left(\frac{\sigma_0^2}{N^2}\right)
\]

### Conclusion
By using position - scaling - eigenfunctions, the paper provides an efficient and accurate method to approximate MLWFs. This method is not only numerically consistent with the results of the traditional MLWFs method, but also has obvious advantages in computational efficiency. In addition, the paper also provides detailed theoretical analysis and error estimates, providing a solid theoretical basis for the reliability and applicability of this method.