Wannier interpolation of reciprocal-space periodic and nonperiodic matrix elements in the optimally smooth subspace

Maximally localized Wannier functions use the gauge freedom of Bloch wave functions 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 nonperiodic 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 nonperiodic 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 nonperiodic tensors of any order without any information on the position of the Wannier centers in real space.