No Arabic abstract
In a tight-binding lattice model with $n$ orbitals (single-particle states) per site, Wannier functions are $n$-component vector functions of position that fall off rapidly away from some location, and such that a set of them in some sense span all states in a given energy band or set of bands; compactly-supported Wannier functions are such functions that vanish outside a bounded region. They arise not only in band theory, but also in connection with tensor-network states for non-interacting fermion systems, and for flat-band Hamiltonians with strictly short-range hopping matrix elements. In earlier work, it was proved that for general complex band structures (vector bundles) or general complex Hamiltonians---that is, class A in the ten-fold classification of Hamiltonians and band structures---a set of compactly-supported Wannier functions can span the vector bundle only if the bundle is topologically trivial, in any dimension $d$ of space, even when use of an overcomplete set of such functions is permitted. This implied that, for a free-fermion tensor network state with a non-trivial bundle in class A, any strictly short-range parent Hamiltonian must be gapless. Here, this result is extended to all ten symmetry classes of band structures without additional crystallographic symmetries, with the result that in general the non-trivial bundles that can arise from compactly-supported Wannier-type functions are those that may possess, in each of $d$ directions, the non-trivial winding that can occur in the same symmetry class in one dimension, but nothing else. The results are obtained from a very natural usage of algebraic $K$-theory, based on a ring of polynomials in $e^{pm ik_x}$, $e^{pm ik_y}$, . . . , which occur as entries in the Fourier-transformed Wannier functions.
Orbital magnetic susceptibility involves rich physics such as interband effects despite of its conceptual simplicity. In order to appreciate the rich physics related to the orbital magnetic susceptibility, it is essential to derive a formula to decompose the susceptibility into the contributions from each band. Here, we propose a scheme to perform this decomposition using the modified Wannier functions. The derived formula nicely decomposes the susceptibility into intraband and interband contributions, and from the other aspect, into itinerant and local contributions. The validity of the formula is tested in a couple of simple models. Interestingly, it is revealed that the quality of the decomposition depends on the degree of localization of the used Wannier functions. The formula here complements another formula using Bloch functions, or the formula derived in the semiclassical theory, which deepens our understanding of the orbital magnetic susceptibility and may serve as a foundation of a better computational method. The relationship to the Berry curvature in the present scheme is also clarified.
We present an algorithm for the adaptive tetrahedral integration over the Brillouin zone of crystalline materials, and apply it to compute the optical conductivity, dc conductivity, and thermopower. For these quantities, whose contributions are often localized in small portions of the Brillouin zone, adaptive integration is especially relevant. Our implementation, the woptic package, is tied into the wien2wannier framework and allows including a many-body self energy, e.g. from dynamical mean-field theory (DMFT). Wannier functions and dipole matrix elements are computed with the DFT package Wien2k and Wannier90. For illustration, we show DFT results for fcc-Al and DMFT results for the correlated metal SrVO$_3$.
Within this paper we outline a method able to generate truly minimal basis sets which describe either a group of bands, a band, or even just the occupied part of a band accurately. These basis sets are the so-called NMTOs, Muffin Tin Orbitals of order N. For an isolated set of bands, symmetrical orthonormalization of the NMTOs yields a set of Wannier functions which are atom-centered and localized by construction. They are not necessarily maximally localized, but may be transformed into those Wannier functions. For bands which overlap others, Wannier-like functions can be generated. It is shown that NMTOs give a chemical understanding of an extended system. In particular, orbitals for the pi and sigma bands in an insulator, boron nitride, and a semi-metal, graphite, will be considered. In addition, we illustrate that it is possible to obtain Wannier-like functions for only the occupied states in a metallic system by generating NMTOs for cesium. Finally, we visualize the pressure-induced s to d transition.
The modern theory of polarization allows for the determination of the macroscopic end charge of a truncated one-dimensional insulator, modulo the charge quantum $e$, from a knowledge of bulk properties alone. A more subtle problem is the determination of the corner charge of a two-dimensional insulator, modulo $e$, from a knowledge of bulk and edge properties alone. While previous works have tended to focus on the quantization of corner charge in the presence of symmetries, here we focus on the case that the only bulk symmetry is inversion, so that the corner charge can take arbitrary values. We develop a Wannier-based formalism that allows the corner charge to be predicted, modulo $e$, only from calculations on ribbon geometries of two different orientations. We elucidate the dependence of the interior quadrupole and edge dipole contributions upon the gauge used to construct the Wannier functions, finding that while these are individually gauge-dependent, their sum is gauge-independent. From this we conclude that the edge polarization is not by itself a physical observable, and that any Wannier-based method for computing the corner charge requires the use of a common gauge throughout the calculation. We satisfy this constraint using two Wannier construction procedures, one based on projection and another based on a gauge-consistent nested Wannier construction. We validate our theory by demonstrating the correct prediction of corner charge for several tight-binding models. We comment on the relations between our approach and previous ones that have appeared in the literature.
Assume that samples of a filtered version of a function in a shift-invariant space are avalaible. This work deals with the existence of a sampling formula involving these samples and having reconstruction functions with compact support. Thus, low computational complexity is involved and truncation errors are avoided. This is done in the light of the generalized sampling theory by using the oversampling technique: more samples than strictly necessary are used. For a suitable choice of the sampling period, a necessary and sufficient condition is given in terms of the Kronecker canonical form of a matrix pencil. Comparing with other characterizations in the mathematical literature, the given here has an important advantage: it can be reliable computed by using the GUPTRI form of the matrix pencil. Finally, a practical method for computing the compactly supported reconstruction functions is given for the important case where the oversampling rate is minimum.