Abstract:Maximally localized Wannier functions use the gauge freedom of Bloch wavefunctions to define the optimally smooth subspace with matrix elements that depend smoothly on crystal momentum. The associated Wannier functions are real-space localized, a feature often used to Fourier interpolate periodic observables in reciprocal space on ultradense momentum grids. However, Fourier interpolation cannot handle non-periodic quantities in reciprocal space, such as the oscillator strength matrix elements, which are crucial for the evaluation of optical properties. We show that a direct multidimensional interpolation in the optimally smooth subspace yields comparable accuracy with respect to Fourier interpolation at a similar or lower computational cost. This approach can also interpolate and extrapolate non-periodic observables, enabling the calculation of optical properties on ultradense momentum grids. Finally, we underline that direct interpolation in the optimally smooth subspace can be employed for periodic and non-periodic tensors of any order without any information on the position of the Wannier centers in real space.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the interpolation and extrapolation of non - periodic matrix elements in reciprocal space. Specifically, the traditional Wannier interpolation method can only handle observables with Brillouin - zone periodicity, such as electronic structure, Berry phase, Hall conductivity, and electron - phonon interaction. However, for non - periodic observables, such as oscillator strength matrix elements that determine optical absorption and photoluminescence properties, there is currently a lack of effective interpolation or extrapolation techniques. The solution to these problems is crucial for improving the efficiency and accuracy in material calculations.
The paper proposes a method for direct multi - dimensional interpolation within the optimal smoothing subspace. This method can handle not only periodic observables but also non - periodic matrix elements. Through this method, optical properties can be calculated on an ultra - dense momentum grid while maintaining a computational cost comparable to or even lower than that of Fourier interpolation. Moreover, this method can be applied to periodic and non - periodic tensors of any order without the need to know the position information of Wannier centers in real space.
The core of this method lies in using the maximally - localized Wannier functions (MLWF) to define an optimal smoothing subspace, in which the dependence of Bloch wave functions on crystal momentum is smooth. By directly performing multi - dimensional interpolation within this subspace, the step of converting to real space required in traditional methods can be avoided, thereby effectively handling non - periodic matrix elements.
The paper verifies the effectiveness of this method through its application in MgB₂ and monolayer electron - doped MoS₂. The results show that for periodic observables, multi - dimensional direct interpolation is slightly better than the standard interpolation method, and the computational overhead is comparable or lower; for non - periodic matrix elements, this method can obtain accurate interpolation results at a reasonable cost, which is difficult to achieve using Fourier interpolation.