Good Wannier bases in Hilbert modules associated to topological insulators
Matthias Ludewig, Guo Chuan Thiang

TL;DR
This paper investigates the existence of smooth, well-localized Wannier functions in spectral subspaces of operators relevant to topological insulators, linking their existence to $K$-theoretic invariants of associated Hilbert modules.
Contribution
It establishes a general criterion connecting Wannier basis existence to the freeness of Hilbert modules over group $C^*$-algebras, highlighting the role of $K$-theory in topological insulators.
Findings
Existence of Wannier bases is equivalent to Hilbert module freeness.
$K$-theoretic invariants classify topological phases.
Provides criteria for Wannier basis construction in complex systems.
Abstract
For a large class of physically relevant operators on a manifold with discrete group action, we prove general results on the (non-)existence of a basis of smooth well-localised Wannier functions for their spectral subspaces. This turns out to be equivalent to the freeness of a certain Hilbert module over the group -algebra canonically associated to the spectral subspace. This brings into play -theoretic methods and justifies their importance as invariants of topological insulators in physics.
| acting on Euclidean | Non-abelian acting on Riemannian manifold |
|---|---|
| Unit cell | Fundamental domain |
| Continuous functions on quasi-momentum space / Brillouin torus | Group -algebra |
| Bloch–Floquet transform | , Eq. (15) |
| Continuous sections of Bloch bundle over | Hilbert -module inside |
| Hamiltonian with Band structure | -elliptic operator on |
| Smooth functions on | The subalgebra |
| Finite-rank spectral subbundle of Bloch bundle (“spectrally isolated bands”) | as a finitely generated projective submodule of |
| Topologically (resp. smoothly) trivialialisable bands | Free Hilbert -module (resp. pre-Hilbert -module) |
| Chern/topological -theory class of bands | Noncommutative Chern class/ operator -theory class of |
| Orthonormal Wannier basis functions /“atomic limit” | Free (pre)-Hilbert module generators |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Good Wannier bases in Hilbert modules associated to topological insulators
Matthias Ludewig
Guo Chuan Thiang
Abstract
For a large class of physically relevant operators on a manifold with discrete group action, we prove general results on the (non-)existence of a basis of well-localised Wannier functions for their spectral subspaces. This turns out to be equivalent to the freeness of a certain Hilbert module over the group -algebra canonically associated to the spectral subspace. This brings into play -theoretic methods and justifies their importance as invariants of topological insulators in physics.
1 Introduction
In solid state physics, one often studies the Schrödinger equation on with a potential which is periodic under the action of a lattice of translations preserving an underlying crystalline structure (the atomic positions, say). The spectral subspace corresponding to the spectrum of the Hamiltonian operator lying between some spectral gaps is then invariant under the translations. A (composite) Wannier basis for comprises a set of wavefunctions and their translates by , such that the set
[TABLE]
is an orthonormal basis for , thus identifying as copies of the regular representation of sitting inside . Wannier bases are convenient for expanding the “effective” electron states in a spectral subspace of physical interest, and it is often desirable to choose the to be as localised as possible, so that one may reasonably think of the as “atomic orbitals” localised at the atomic positions labelled by , see Fig. 1 for an illustration.
Wannier basis construction usually proceeds via the Bloch–Floquet decomposition over the character space of the translation symmetry group (reviewed in §7.1). From this vantage point, much effort has been devoted to proving the existence, for arbitrary , of “good” Wannier bases with, say, exponential decay [31, 12, 38, 8, 33, 37, 13, 14]. Remarkably, in , there is a topological obstruction — the Chern class of the so-called Bloch bundle over [8] — which persists even if we relax the the decay condition significantly (cf. Remark 3.5). The non-existence of good Wannier bases, or “atomic limits” [7], is a paradigmatic feature of so-called topological insulators. When a boundary is subsequently introduced, the slow decay rate of bad Wannier bases means that there is no meaningful length scale to decouple the “bulk” degrees of freedom from the “boundary” ones. This is one way to see why a bulk-boundary correspondence [6, 18, 28, 32, 35, 25, 40] should be expected in topological insulators.
How much of these insights depend on having an abelian symmetry group and the classical Fourier transform? In this paper, we abstract the salient features of the “good Wannier basis existence problem”, and demonstrate that the above topology/localization dichotomy [8, 33, 37, 14] holds much more generally, for nonabelian symmetry groups, and also for projective symmetry group representations (whether the group is abelian or not) which occur for quantum Hall Hamiltonians (see §7.4).
Basic Setup. Let be a complete connected Riemannian manifold with an effective, cocompact, properly discontinuous, isometric action of a countable group . Let be an complex -equivariant vector bundle over with a -invariant fiber metric.
The relevant space of fields is then the space of square-integrable sections111If is just a trivialised line bundle, we will simply write , omitting reference to . of , which admits a right action of the reduced group -algebra (see §2) in a canonical way. Now given a Hamiltonian , i.e. a self-adjoint operator acting on sections of , we assume that we are given a compact subset of the spectrum of , which is separated from the rest of the spectrum by spectral gaps. A typical example has the spectrum of bounded from below, and the so-called Fermi level lying in a spectral gap — this describes an insulator. The spectral subspace for energies below the Fermi level is the subspace of occupied energy states, and is of particular physical interest in determining various material properties.
The idea is now to consider the spectral subspace as a module over the reduced group -algebra . However, in general, this space is too large to be finitely generated and projective. In the case , where , continuous functions on the Brillouin torus, this parallels the fact that is identified with the space of square-integrable (not necessarily continuous) sections of the Bloch bundle; clearly this is not finitely generated as a module over and moreover does not carry any topological information.
To remedy the situation, we use a construction of Roe [42, pp. 243] which identifies a dense subspace , which is a Hilbert module over , in other words, admits a -valued scalar product. The construction is such that in the abelian case , the intersection is precisely the space of continuous sections of the Bloch bundle, which is finitely generated and projective as a -module, by the theorem of Serre–Swan.
We then prove the following result.
Theorem 1.1**.**
Assume in addition that has polynomial growth. Let be a self-adjoint equivariant differential operator acting on sections of , which is either of Laplace type, or first order elliptic. Suppose that is a compact subset of the spectrum which is separated from the rest of the spectrum and let be the corresponding spectral subspace. Then the subspace
[TABLE]
is a finitely generated, projective -module. Moreover, the following assertions are equivalent.
- ()
* is a free -module of rank .* 2. ()
There exist functions such that the set , , is an orthonormal basis of .
*Here denotes the space of smooth sections of which decay faster than any polynomial, together with their derivatives, c.f. Def. 5.1 below. *
In fact, we also prove a stronger, quantitative version of the theorem above, generalizing the results of Kuchment [33] (see also [14]): Namely, we show that for any , the existence of a module such that is free of rank , is equivalent to the existence of a tight frame of consisting of “good” Wannier functions and their -translates; c.f. §6.
The above results do not require to be abelian, so they apply, for instance, to the theory of crystalline topological insulators in solid state physics, for which , the Euclidean space, and is a generally non-abelian crystallographic group (with trivial). Any such group has polynomial growth; we remark that by Gromov’s theorem [22], groups of polynomial growth are precisely those having a nilpotent subgroup of finite index. Our results show that the -theory of , which is a proposed way to classify crystalline topological insulators [19, 46], presents computable obstructions to the existence of good Wannier bases (in the sense of being in the Schwartz class). Furthermore, these topological insulator inspired ideas extend to the non-Euclidean setting, as we exhibit in an example in §7.5. Let us mention that -algebra methods have previously been applied to the mathematics of topological insulators in (implicitly) Euclidean settings, see e.g. [4, 6, 18, 28, 32, 40], and independently of the Wannier basis construction problem.
We remark that we could also directly consider as a module over the group von-Neumann-algebra , which may seem simpler than considering the intersection necessary in the -case. However, we lose topological information: For example, in the case , we will always have as -modules, corresponding to the fact that the Bloch bundle aways has a square-integrable (non-continuous) trivialization.
After fleshing out the constructions and proof of the main theorem in §2-§6, we formulate the general “good Wannier basis existence problem” (Problem 1) in the final section §7, and apply our theory to compute explicitly the presence/absence of obstructions in several physically interesting new examples. For the convenience of readers who are not acquainted with noncommutative geometry ideas, we include Table 1 to convert some familiar ideas in usual Bloch theory to the setup of non-commutative geometry.
2 Algebras of rapidly decaying sequences
Let be a finitely generated group. Its group algebra is the set of all finite formal linear combinations of the group elements . It has a -involution given by
[TABLE]
The group algebra acts on the Hilbert space
[TABLE]
by translation (i.e. left regular representation),
[TABLE]
Multiplication by is bounded, hence one obtains a representation of on the bounded operators , which is in fact a -representation. The reduced group -algebra is then the -algebra obtained by completing with respect to the operator norm. Denoting by the element which is one at the unit of and zero otherwise, the map
[TABLE]
provides a continuous embedding of into , and further into the space of all bounded, -indexed sequences. In particular, elements of are bounded -indexed sequences.
Let be a finite generating set of with and , and let
[TABLE]
be the corresponding length function. Now for any auxiliary Banach space , we define
[TABLE]
the space of rapidly decreasing sequences indexed by with values in .
Remark 2.1**.**
Above, we wrote for the pointwise norms of the -valued sequence , and we will also write for the sequence of pointwise -norms. When , we simply write , which can be thought of as the space of as “smooth functions on the dual of ”. For example, comprises the Fourier coefficient sequences of smooth functions on the dual circle .
The space is topologized by the increasing sequence of Hilbert space norms
[TABLE]
which turn into a Fréchet space. One easily checks that the definition of is independent of the choice of finite generating set ; in particular, any such choice gives rise to an equivalent set of norms (3).
Lemma 2.2**.**
Assume that has polynomial growth, meaning that
[TABLE]
*for some large enough. Then for any -algebra , the space is a Fréchet algebra in a natural way, which is continuously included in the tensor product222For two -algebras and , their algebraic tensor product sits inside and its completion therein is the spatial tensor product, see Appendix T.5 of [48] for details of this and other possible -algebra tensor products. In our case where has polynomial growth and is hence amenable, its group -algebra is nuclear. So for any -algebra , there is actually a unique -algebra tensor product . For details, see [9] Example 2.6.6, Thm. 2.6.8 and Thm. 3.8.7. . *
Proof*.*
We first show that there exist such that
[TABLE]
for all , which implies that is continuously included in . Namely, for (the algebraic tensor product), we have by the triangle inequality that
[TABLE]
ence by the Cauchy–Schwarz inequality,
[TABLE]
with the second -independent factor finite by (4). Since is dense in and in , this shows the claim.
Clearly, is an algebra in the obvious way. To show that the multiplication is continuous, we will show that for any , there exist such that
[TABLE]
for all . It suffices to consider the case . Then
[TABLE]
Using the triangle inequality
[TABLE]
we obtain
[TABLE]
where are the respective pointwise -norms of ; moreover, we used the estimate (5). Since the norms (3) are increasing in strength, this implies (6) with .
Remark 2.3**.**
More generally, a group is said to have property (RD), if there exist such that for all , which implies . The lemma above shows that groups of polynomial growth have property (RD). Conversely, it is known that an amenable group has property (RD) if and only if it is of polynomial growth [27]. The proof above also shows that also under the assumption (RD), is a Fréchet algebra for any -algebra , i.e. that we have estimates of the form (6).
A subalgebra of a -algebra is called spectral, if for any that is invertible in , we already have . Moreover, we say that is closed under holomorphic functional calculus when for any function that is is holomorphic on an open neighborhood of the spectrum of (in ).
Proposition 2.4**.**
*Suppose that has polynomial growth. Then for any unital -algebra , the dense subalgebra of is spectral and closed under holomorphic functional calculus. *
In order to prove this proposition, we use the following criterion of Ji, Thm. 1.2 of [26]. In fact, Ji proves a version of Prop. 2.4 for the special case where , a finite-dimensional algebra, but under the weaker condition that only satisfies property (RD). We adapt his proof to show that under the stronger condition of polynomial growth, one can allow to be arbitrary.
Proposition 2.5**.**
Let be a unital -algebra and be a closed subalgebra containing the unit of . Let
[TABLE]
*be a closed, unbounded derivation. Then is a subalgebra of , and if is dense in , then it is spectral in and closed under holomorphic functional calculus. *
We remark that the theorem of Ji does not include the additional statement that is closed under holomorphic functional calculus. However, this follows easily from an additional argument by Valette, see [47], below Prop. 8.12.
Proof* (of Prop. 2.4).*
Let for some Hilbert space . We will use Prop. 2.5 for and . To get our hands on a derivation , we define the densely defined, unbounded operator
[TABLE]
given by pointwise multiplication with the length function , where we set , making densely defined. Now define
[TABLE]
Hence for any , the operator extends by continuity to an element of , again denoted by ; the operator is then given by
[TABLE]
It is straightforward to show by induction that for any , we have the formula
[TABLE]
for and (note that this is a well-defined, not necessarily square-summable, -indexed sequence, as for fixed , the sum over is in fact finite). From the triangle inequality (7), we obtain for , that
[TABLE]
Here denotes the element of (as a bounded -indexed sequence) obtained by taking the pointwise -norm of ; similarly for . This shows that for each , hence the intersection of all is dense. Moreover, since has polynomial growth, Lemma 2.2 implies that there exist such that
[TABLE]
This shows that
[TABLE]
which implies that , for any .
To see that is closed, let converge to some . We claim that this means in particular that in for every . Namely, given , denote by the element of which is at the unit element of and zero otherwise. Then we have
[TABLE]
which converges to zero as , proving the claim. Having established this fact, formula (8) implies that for all ,
[TABLE]
pointwise, in the sense that the coefficients converge in for each ; here is some -indexed sequence, not necessarily square-summable. But now if converges to , then in particular, converges pointwise to , for all . Hence in fact for each . We obtain that , hence is densely defined, and moreover, is bounded, since is bounded. This shows , and , hence is a closed operator.
With a view on the result of Ji cited above, it remains to show that
[TABLE]
in order to show that is spectral in , and it is this step where we need the stronger condition of polynomial growth instead of just property (RD). Namely, let for any . Given , denote by the element of which is at the unit element of and zero otherwise. Then as , we have
[TABLE]
in fact
[TABLE]
for some constant , depending only on . In particular, we have that , for any . Since has polynomial growth, this implies that .
3 The Hilbert module of a group action
In this section, we review the notion of (Pre-)Hilbert modules over -algebras, and introduce the basic examples relevant to our paper.
Let be a pre--algebra, by which we mean a -closed subalgebra of a -algebra. By a pre-Hilbert -module, we mean a right -module together with an -valued inner product,
[TABLE]
which is -antilinear in the first argument, -linear in the second argument, and satisfies the identities
[TABLE]
for all and . Moreover, we require that lies in the cone of positive elements of , for each ; here we say that an element of is positive if it is positive in the completion of . Using the -property of the norm of , one deduces the Cauchy–Schwarz inequality
[TABLE]
which implies that the positive homogeneous functional
[TABLE]
satisfies the triangle inequality, hence defines a norm on . If is complete with respect to the norm (11), it is called a Hilbert -module over or simply a Hilbert -module. Otherwise, starting from a pre-Hilbert -module , one can “simultaneously complete” and to obtain a Hilbert -module, see pp. 4–5 of [34].
Example 3.1** (Standard Hilbert -module).**
Let be a separable Hilbert space and a -algebra. The interior tensor product is defined as follows (for general facts regarding the interior product, which in this case happens to coincide with the exterior product, see e.g. [34, Chapter 4]). On the algebraic tensor product of and , define the inner product
[TABLE]
We then denote by the completion of the algebraic tensor product with respect to the norm given by (11) in terms of this inner product. For , we obtain the standard (countably-generated) Hilbert -module discussed in [48, p. 237]. If we take , we obtain the standard finitely generated free module .
Example 3.2**.**
Let be a finitely generated group having polynomial growth or, more generally, property (RD), so that the space is a -subalgebra of , c.f. Prop. 2.2 and Remark 2.3. Given a separable Hilbert space , the space as defined in Eq. (2) is a right -module in a natural way. Moreover, with the inner product
[TABLE]
becomes a pre-Hilbert -module. This inner product coincides with the one defined in (12) above for . The inner product indeed takes values in : For any , we have
[TABLE]
hence
[TABLE]
using the estimate (6). We also have , because
[TABLE]
using estimate (5) together with (14). Of course, the inner product (13) on coincides with the restriction of the inner product of from Example 3.1; this shows that (13) has all required properties of a -valued inner product.
For our last example, we assume the Basic Setup from the introduction; in other words, let be a complete Riemannian manifold and let be a group acting properly discontinuously on by isometries, such that is compact. Moreover, let be a -equivariant vector bundle on with a -equivariant fiber metric. Denote by the space of square-integrable sections of . The data induce an isometric right action of on by pullback, explicitly,
[TABLE]
for , , . With respect to a choice of fundamental domain for the action of , we obtain a -equivariant isometry
[TABLE]
where the tensor product is a tensor product of Hilbert spaces; here acts on by right multiplication. Notice that we need the action to be effective in order to ensure surjectivity of . Of course, these right -actions extend in the obvious fashion to right -module actions of the group algebra .
Example 3.3** (The module ).**
We now describe a construction due to Roe, cf. pp. 243 of [42] of a Hilbert -module associated to the action of the group on . As a space, is a dense subset of .
On , we consider the pairing
[TABLE]
a priori taking values in the space of bounded sequences indexed by . Now starting with compactly supported , we see that with positive (in ), and that the pairing satisfies (9), so becomes a pre-Hilbert -module (Lemma 2.1 of [42]). We then define
[TABLE]
One shows that the right action of preserves and that the bracket, restricted to , takes values in . This turns into a Hilbert -module.
We remark that if is the element which is one at the unit element and zero otherwise, then
[TABLE]
hence which shows that in fact, the completion (17) can be realized as a subspace of . Therefore, we always stipulate .
Given a fundamental domain , we see that defined in (15) restricts to a vector space isomorphism
[TABLE]
Moreover, the bracket (16) can be equivalently written as
[TABLE]
which shows that maps the bracket of to the bracket (12). Extending this observation by continuity, one has the following result.
Proposition 3.4**.**
After choosing a fundamental domain for the -action, restricts to an isomorphism of Hilbert -modules
[TABLE]
*where on the right hand side, we have the Hilbert -module from Example 3.1. *
Remark 3.5**.**
It is easy to see that if has the decay condition
[TABLE]
with respect to some fundamental domain , then . Condition Eq. (18) was considered in [33], for , and trivial , in connection with Chern class obstructions to localised Wannier bases (see §7 for more on Wannier bases). Also in this abelian setting with , another decay condition
[TABLE]
was considered in [37] for the related problem of finding the optimal decay of a (composite) Wannier basis for a possibly topologically non-trivial Bloch bundle. Amongst other results, it was found that any could always be achieved whether or not the Bloch bundle is topologically trivial; the threshold case , corresponding to -trivialisability of the Bloch bundle, is attainable by some Wannier basis exactly when the Bloch bundle is topologically trivialisable. Unfortunately, it seems that there is no simple relationship between the decay condition (19) and membership in our .
3.1 Comparison with literature
Roe constructed the Hilbert -module of Example 3.3 for the mildly less general case where is a trivial line bundle, but with allowed to be a general proper metric space with a proper cocompact isometric -action [42]. One motivation of [42] was to connect coarse geometry to the analytic assembly map in the arena of the Baum–Connes conjectures; in particular, Lemma 2.3 of [42] provides an identification of the compact operators on this Hilbert -module with the equivariant Roe algebra (completion in of the -equivariant locally compact, finite propagation operators). In a subsequent work of the authors [35], the algebra plays a prominent role in proving, via coarse index methods, the robustness of edge states for Chern topological insulators of generic shapes .
One important role of the bundle is to allow for gauge (associated) bundles, especially in the context of electromagnetism (with as local gauge group) and quantum Hall Hamiltonians/magnetic Schrödinger operators. Here, the -action on is allowed to be lifted to a projective unitary representation on , as encoded by a -valued group 2-cocycle . For free, properly discontinuous, cocompact -actions on a Riemannian manifold , Gruber (§V.B of [23]) made constructions almost identical to Example 3.3, modified slightly to the projective-linear context, to produce inside a standard Hilbert -module over the twisted group algebra . Note that is generally noncommutative, even if is an abelian group such as (see §7.4). The paper [23] uses Hilbert -module language to study noncommutative Bloch theory (also studied in [10, 36]), with the intention to understanding some spectral properties of abstract -elliptic operators on Hilbert -modules, with a general unital -algebra.
The commutative Bloch theory setting also applies to under a so-called rational flux condition — for such rational cocycles , the algebra is Morita equivalent to . For example, the discrete magnetic translations on corresponding to rational fluxes per unit cell will generate such a . A suitable superlattice (a subgroup of whose fundamental domain has integer flux though it) will act linearly rather than projective linearly, and commutative Bloch theory can be carried out with respect to this superlattice. In this context, the Bloch–Floquet transform converts into a Hilbert bundle over , and topological aspects of its spectral subbundles (sometimes called Bloch bundles) as modules over is studied in [15], especially §7 therein. The availability of a bona fide topological space whose bundle theory mirrors the module theory of is very useful for physics applications, including finding the strongest statement of the so-called localisation-topology dichotomy/correspondence (e.g. [37]).
Our focus, however, is on the setting in which no such topological space exists for or , yet we can state a reasonable form of the localisation-noncommutative topology dichotomy/correspondence (e.g. Theorem 1.1, and examples of applications in §7).
4 Finitely generated projective pre-Hilbert modules
For a unital -algebra , it is known, cf. Thm. 15.3.8 of [48] that there is a correspondence between finitely generated projective (f.g.p.) -modules and projections in , the algebra of compact adjointable module maps on the Hilbert -module , c.f. Example 3.1. Explicitly, this means that given a projection , one has that the Hilbert -submodule is isomorphic to for some .
For the dense subalgebra , we can achieve a similar result for pre-Hilbert -modules.
Lemma 4.1**.**
*Let be a finitely generated group of polynomial growth. Let be a projection in , where is a separable Hilbert space. Then is isomorphic as a pre-Hilbert -module to a f.g.p. -module. *
Proof*.*
By Prop. 2.4, is spectral and closed under holomorphic functional calculus in for as well as , the -algebra of compact operators on with a unit adjoined. The proof is then essentially an adaptation of Lemma 15.4.1 of [48].
We may assume that and write . Let denote the standard rank projection in , tensored with the identity in . As is an approximate identity for , after writing for the truncation (no longer a projection in general but still self-adjoint), we can arrange for
[TABLE]
Furthermore, is positive with , so that
[TABLE]
Note that can be regarded as an element of . If is in the spectrum of (which can be taken in or in since the first is spectral in the latter), then by the above calculation, . This implies that
[TABLE]
so has a spectral gap at . The function which is 0 to the left of and 1 to the right is therefore holomorphic in a neighbourhood of the spectrum of , so we obtain a projection which satisfies by construction (that follows because is closed under holomorphic functional calculus). Moreover, we have
[TABLE]
by the choice of . Therefore
[TABLE]
and a standard construction (e.g. Prop. 5.2.6 in [48]) gives
[TABLE]
implementing , where is the unital -algebra given by adjoining a unit to . Then , so that is invertible in . In fact, we have , as by Prop. 2.4, is spectral in . Note that gives and also , so with being the unitary in the polar decomposition of , which again exists by Prop. 2.4, we obtain a unitary equivalence in .
The unitary constructed above gives an isomorphism
[TABLE]
Note that since , we have , hence . Clearly, commutes with right multiplication by , hence is a module map. Moreover, since is unitary, we have
[TABLE]
hence preserves the inner product. This finishes the proof.
Of course, Lemma 4.1 also implies that is unitarily isomorphic as a Hilbert -module to a direct summand of . A priori however, the minimal required could be smaller in the -module case than in the -module case. The next lemma shows that this is not the case.
Lemma 4.2**.**
*Let be a projection in as in Lemma 4.1. If is isomorphic as a Hilbert -module to a direct summand of for some , then is isomorphic as a pre-Hilbert -module to a direct summand of for the same . *
In particular, is (isomorphic to) a pre-Hilbert -module of rank if and only if is (isomorphic to) a free Hilbert -module of rank .
Proof*.*
Lemma 4.1 gives and for some , with a projection in . Now suppose a unitary implements a further reduction
[TABLE]
Choose a unitary close to so that is a projection in close to .
On the other hand, pick (not necessarily a projection) close to such that its spectrum avoids the line . As in the proof of Lemma 4.1, holomorphic functional calculus gives a projection close to .
Now the projections and are in and are both close to , hence close to each other; therefore, they are unitarily equivalent in (take the polar decomposition of the invertible ). This gives a unitary such that
[TABLE]
which finishes the proof.
The above results combine to the following corollary.
Corollary 4.3**.**
Let be a projection in as in Lemma 4.1, and let be the corresponding -module. Then is a finitely generated and projective Hilbert -module. Moreover, for each , there exists a Hilbert -module such that
[TABLE]
if and only if there exists a pre-Hilbert -module such that
[TABLE]
**
5 Rapidly decaying functions and admissible operators
Let be a complete Riemannian manifold and let be a Hermitian vector bundle over . In this section, we consider the following function spaces inside , whose definitions refer to a choice of basepoint , but are in fact independent of it. Later, we will relate these function spaces under the assumption of polynomial growth.
Definition 5.1**.**
For some fixed , we define
[TABLE]
The first is the space of -sections of having rapid decay, while the second one is the space of -valued Schwartz functions on .
If we are in our Basic Setup from the introduction, which we will assume from now on, the Milnor–Švarc lemma applies [16, Thm. 8.37], so must be finitely generated, and, as a metric space with its word metric, be quasi-isometric to . In other words, there exist constants such that
[TABLE]
for all and all . Of course, these constants (as well as the length function) depend on the choice of finite generating set for (c.f. Section 2). We moreover assume that has polynomial volume growth, meaning that there exists such that for some (equivalently, for any) , we have
[TABLE]
for some constant and all , where is the ball around of radius . The estimates (20) imply that this condition is equivalent to the condition that has polynomial growth in the sense of (4). The estimate (21) now implies the inclusions
[TABLE]
in fact, the second inclusion follows from the following lemma.
Lemma 5.2**.**
For any bounded fundamental domain, the map defined in (15) restricts to an isomorphism of pre-Hilbert -modules
[TABLE]
*where . *
Proof*.*
For , let be the set of all such that is not included in . Then for and any , we have the estimate
[TABLE]
Together with (20), this shows that if , then .
To see the converse, note that for , we have , hence (20) implies that
[TABLE]
here we assume that . Using this, we obtain for any and any
[TABLE]
If now , this sum converges, showing that .
Wrapping up the results from above, with respect to the choice of a bounded fundamental domain, we have the following diagram, where each of the horizontal arrows, given by the map , is an isomorphism preserving the respective (pre-)Hilbert module structures.
[TABLE]
Remark 5.3**.**
Under the condition of polynomial growth, is in fact a nuclear space [27, Thm. 3.1.7], hence the algebraic tensor product with any other Banach space has a unique tensor product topology. In particular, the tensor product is unambiguously defined and equal to H^{\infty}\bigl{(}\Gamma,L^{2}({\mathcal{F}},E)\bigr{)}. However, we will not need this fact.
Definition 5.4**.**
We say that an operator is admissible if it has a smooth integral kernel which is rapidly decaying in the sense that for and each , there exists a constant such that
[TABLE]
for all .
Lemma 5.5**.**
*If is an admissible operator and , then . *
Proof*.*
Let be the integral kernel of and let . We have
[TABLE]
which is absolutely convergent by the polynomial growth condition on and the decay of . Now for any , we have
[TABLE]
By the triangle inequality, we have for any that
[TABLE]
which is integrable with respect to for large enough, since has polynomial growth. Thus for any and , is bounded independent from , which was what we needed to show.
Proposition 5.6**.**
Let be an admissible -invariant projection. Then is a f.g.p. Hilbert submodule of . Moreover, given , there exists a f.g.p. Hilbert -module such that
[TABLE]
if and only if there exists there exists a f.g.p. pre-Hilbert -module such that
[TABLE]
**
Proof*.*
With a view on (23), the result follows from Corollary 4.3, if we can show that
[TABLE]
for some bounded fundamental domain . Corollary 4.3 states that the existence of a f.g.p. Hilbert -module as above implies the existence of a f.g.p. pre-Hilbert -module such that and vice versa.
However, if for , then with . Hence, since is a projection, , which by Lemma 5.5 is contained in , as . In other words, we have
[TABLE]
To show (25), let be the integral kernel of . Fix , and let . We have
[TABLE]
where we used that the integral kernel of is -invariant, for all , . Hence
[TABLE]
In other words, we have
[TABLE]
which is a compact operator, as it has a bounded integral kernel. For the operator norm of , we have the estimate
[TABLE]
for any , where we used that is admissible, and that is bounded, together with (20). This follows from the fact that is admissible and that is quasi-isometric to . This implies (25) since has polynomial growth.
6 The main theorem
Let be as in the Basic Setup of the introduction. Thm. 1.1 from the introduction is a special case of the following result, Thm. 6.1, combined with Thm. 6.3 further below.
Recall that a tight frame of a Hilbert space is a collection of (possibly linearly dependent) elements of such that for all , we have
[TABLE]
Theorem 6.1**.**
Suppose that has polynomial growth and let be a -invariant, self-adjoint (unbounded) operator on such that is admissible for each Schwartz function . Let be a compact subset of the spectrum of which is separated from the rest of the spectrum and let be the corresponding spectral subspace. Set
[TABLE]
- (i)
* is a f.g.p. Hilbert -module.* 2. (ii)
Suppose that there exists a f.g.p. Hilbert -module such that . Then there exist such that these functions and their translates form a tight frame of . Moreover, if , then we can arrange for this tight frame to be an orthonormal basis. 3. (iii)
Conversely, suppose that there exist that together with their -translates form a tight frame of . Then there exists a f.g.p. Hilbert -module such that , and if the together with their translates are linearly independent, then is a free Hilbert -module.
**
Remark 6.2**.**
Statement (iii) above is false if are not constrained to be in . Indeed, there are examples where there exists an orthonormal basis of consisting of and their translates, without being free, see the case in §7.3. However, Thm. 6.1 states that for a basis comprising and their translates, if the do satisfy the mild decay property of being contained in (c.f. Remark 3.5), then is free, and one can further choose the from the Schwartz space .
Proof*.*
Let be the spectral projection in associated to the subset . Since is separated from the rest of the spectrum, we can write using functional calculus, where is a compactly supported smooth function such that for and whenever . Thus is admissible in the sense of Def. 5.4, and by Prop. 5.6,
[TABLE]
is finitely generated and projective, which is claim (i).
To show claim (ii), let be a f.g.p. Hilbert -module such that . By Prop. 5.6, there exists a pre-Hilbert -module and an isomorphism of pre-Hilbert -modules
[TABLE]
By Lemma 5.5, we have . Define by , where is the projection onto the first factor and is the -th unit vector. Since (being an isomorphism of pre-Hilbert -modules) preserves the inner products, the vectors satisfy
[TABLE]
Composing the -valued inner product with the standard trace on gives a scalar product on , and by (26), the vectors , , , form an orthonormal basis of the completion with respect to this inner product.
The inner product on obtained this way coincides with the restriction of the standard inner product on , and since is dense in , the completion of is . Because the orthogonal projection of an orthonormal basis to a subspace forms a tight frame of the subspace, this shows that the sections , , form a tight frame of . Clearly, if , then the are linearly independent, as then , the are linearly independent and is a vector space isomorphism.
To show (iii), let be such that , , , forms a tight frame of . From the characterization of tight frames [24], there exists a Hilbert space and an orthonormal basis , , of such that . Setting
[TABLE]
for defines a right representation of on . This action preserves ; in fact, if we let , then . Setting
[TABLE]
defines a -valued inner product on the subspace consisting of finite linear combinations of the , , . Letting be the completion of with respect to the norm induced from the inner product and the norm on defines a Hilbert--module such that , via the obvious isomorphism of Hilbert--modules sending to . If the were in fact linearly independent, then , so that .
The following result shows that there are many examples of operators for which Thm. 6.1 applies.
Theorem 6.3**.**
*Suppose that has polynomial growth. Let be a self-adjoint, -invariant differential operator acting on sections of and assume that either is elliptic of order one; or that is of order two and of Laplace type. Then for each Schwartz function , the operator is admissible in the sense of Def. 5.4. *
Here by a Laplace type operator, we mean a second order operator such that in local coordinates on , it is given by
[TABLE]
where is the inverse matrix of the coefficient matrix of the Riemannian metric on and , are certain sections of the endomorphism bundle of . Note that the self-adjointness requirement on poses some additional restrictions on the lower order coefficients and .
Proof*.*
The argument is similar to the one in [11]. We will show that for all , , there exists a constant such that
[TABLE]
whenever is a smooth section of with support in a compact set . Here denotes the set of with and for an open subset , is the Sobolev space of sections in with square-integrable weak derivatives up to order . The result then follows from the results in [17]: The estimates (27) imply that is a quasi-local smoothing operator (c.f. Def. 2.14 ibid.); Sobolev embedding together with the results of Section 2.4 in [17] shows that is admissible in the sense of Def. 5.4 above.
First consider the case that is of order one. In that case, the wave operator has finite propagation, meaning that there exists a constant such that whenever has support in a compact set , then has support in . We use the formula
[TABLE]
for Schwartz functions , where denotes the Fourier transform of . Given , let . Then by elliptic estimates, we have (if one defines norms suitably). Let be the function that is equal to zero on and identically one on the complement. Now
[TABLE]
by the finite propagation of . This gives the estimate (27), since the operator is uniformly bounded independent of and is rapidly decaying.
If is order two and of Laplace type, the spectrum of is bounded below. By possibly replacing by , we may assume that the spectrum of is bounded below by . We may then take to be an even function, in which case formula (28), applied for instead of yields
[TABLE]
here is the Fourier transform of . The point is now that is the solution operator to the wave equation, , which again has finite propagation speed [45, 2]. The proof is then similar to the argument before.
7 Application to good Wannier basis existence problem
We apply our results to the old problem of constructing well-localised Wannier bases in solid state physics. Most existing results on Wannier bases apply to the basic case where is a lattice of translations acting on affine Euclidean space . Fourier transform identifies with , the continuous functions on
[TABLE]
the Brillouin zone/torus in physics, and one studies -invariant Hamiltonian operators acting on the Hilbert space of quantum mechanical wavefunctions.
We will recast the idea of Wannier bases in noncommutative topology/geometry language, so that our results become applicable, e.g. to all crystallographic at once, and also in certain non-Euclidean settings as illustrated by our final example. It is instructive, however, to first recall the basic notions in the commutative case , which proceeds via classical Bloch theory.
7.1 Commutative Bloch–Floquet transform and Wannier bases
The fundamental domain for the action on is an affine torus (not to be confused with the Brillouin torus . The Bloch–Floquet transform is a unitary map , cf. Eq. (15) with , explicitly defined by
[TABLE]
This sum of -shifted versions of weighted by the phase factor is often called a Bloch sum. If we replace in the Bloch sum by , , we get the Bloch wave condition
[TABLE]
Thus each fixed quasimomentum , the function extends to a -quasiperiodic Bloch “wavefunction” on , albeit not normalisable over but only over . Equivalently, we write where each is the line bundle obtained by quotienting by , . Thus in total (see §D.3 of [19]), there is a Hilbert bundle with fibres
[TABLE]
and is an -section of . The action of translation by is represented unitarily on these sections by pointwise multiplication by the continuous function . An inversion formula
[TABLE]
holds, recovering as a “superposition” of Bloch wavefunctions . In reverse, given a section , its Wannier function is the inverse Bloch–Floquet transform (29) of .
Often, a certain finite-rank subbundle of is of interest, e.g. if the spectrum of has band structure, the spectral subspace for spectra lying between some given spectral gaps is a -invariant subspace obtainable (after taking Bloch–Floquet transform) as the -sections of a locally trivial subbundle , called a Bloch bundle [39, 33, 19]. If there are bands (so has rank ), one can always find orthonormal measurable sections for ; then each gives rise to a corresponding Wannier wavefunction such that the translates are mutually orthonormal [33]. So we obtain an orthonormal Wannier basis , , for the spectral subspace of interest.
As mentioned in the introduction, when , the Chern class of obstructs choosing the to be continuous, thereby obstructing the existence of a Wannier basis comprising exponentially decaying wavefunctions [8], or even much more mildly decaying ones [33]. This failure has been turned into a triumph in recent years, because of the experimental discovery and burgeoning theoretical interest in these topological insulators in physics, usually characterised exactly by the topological nontriviality of as detected, e.g. by -theory classes [29, 19, 46].
The same obstruction can occur for a range of decay conditions on the Wannier basis. From a physical perspective, a working definition of a topological insulator is one for which does not admit an “atomic limit” [7], which we can think of as the nonexistence of localised Wannier bases for . For this, a weak decay condition is preferred, so as to argue that a topological insulator necessarily has very delocalised Wannier bases. On the other hand, when the obstruction is not present, we would like to be able to choose Wannier wavefunctions which are as regular (smooth) and/or localised as possible. Let us remark that the rather extreme condition of compactly supported Wannier wavefunctions was considered in [41] in connection with algebraic -theory obstructions, cf. our construction of as a pre-Hilbert -module in §3.
For our purposes, we say that a (Wannier) wavefunction is “good” if it belongs to the Schwartz class , cf. Definition 5.1.
7.2 Existence of good Wannier basis: nonabelian symmetry groups
Quite generally, we can ask the nonabelian- analogue of the good Wannier basis existence problem:
Problem 1**.**
Given a spectral subspace of an admissible -invariant operator on (as defined in §5-6), does there exist a good Wannier basis, i.e. a set which together with their translates , form an orthonormal basis for ?
When has polynomial growth, our Thm. 6.1 says that the dense subspace of “not-so-poorly decaying” functions forms a f.g.p. Hilbert -module. Furthermore, if is freely generated by , we can even choose these to be “good”, i.e. in . Thus, Thm. 6.1 answers Problem 1 in the affirmative, for arbitrary , if all f.g.p. modules over are free. A simple example where this occurs is (see Fig. 1), and a new wallpaper group example is given in the next subsection (see Fig. 3). Generically, there can be f.g.p. -modules which are not free, in which case our Thm. 6.1 answers Problem 1 in the negative — we can only achieve a tight frame if we want to be in or better.
The semigroup structure of f.g.p. -modules is difficult to understand in general, so a tractable first step is to compute . For crystallographic groups, one can appeal to the Baum–Connes conjecture, or use algebraic topology methods after converting to a twisted equivariant -theory group of as in [19]. For physical application, a reasonable justification for restricting the search to stably (non-)free f.g.p. modules can be made, following [29]: may be supplemented by well-localised inner atomic shells (free -modules) which were not accounted for in the specification of . Thus the -theory of provides physically meaningful invariants which label topological phases given a group of symmetries [46].
7.3 Crystallographic group examples
In this subsection, will be a trivial bundle over Euclidean and we omit reference to it. The -invariant Hamiltonians whose spectral subspaces we are considering are assumed to satisfy the generic conditions of Thm. 6.3, so the conclusions of the main theorem 6.1 apply.
Case . The crucial difference between and is that all f.g.p. modules over are free, whereas there are topologically nontrivial bundles over , thus non-free f.g.p. modules over . A vector bundle over (f.g.p. -module) is trivial iff its first Chern class vanishes. Furthermore, the reduced -theory of coincides with in this case, so a nonvanishing reduced -theory class detects not just the failure of a f.g.p. -module to be stably free, but its failure to be free.
The first Chern class obstruction, , leads to “Chern topological insulators” in , for which there is no good Wannier basis for . Note that we can measurably trivialise a topologically non-trivial Bloch bundle, and obtain a bad Wannier basis for , e.g. §3 of [33], as illustrated in Fig. 2. For , there are further -theory/Chern class obstructions to being free/stably free.
Case . An interesting non-abelian example is , where the second copy of acts on the first by the nontrivial automorphism of reflection. In crystallography, this group appears as the 2D crystallographic space group (a “wallpaper group”) , and can be realised as a group of isometries of 2D Euclidean space . Since is torsion free, the action is free and the fundamental domain is a manifold diffeomorphic to the Klein bottle, as illustrated in Fig. 3.
Proposition 7.1**.**
*Any finitely generated projective module over is isomorphic to the free module for some unique . *
Proof*.*
The Klein bottle is a , and invoking the Baum–Connes isomorphism [3] and low-dimension of , we compute that
[TABLE]
The isomorphism takes the generator to the class of the identity projection in , or equivalently, the free rank-1 module, cf. §7 of [47] or §5.3 of [21]. Thus every f.g.p. -module is stably isomorphic to with its -theory class. Now, we note that has a canonical faithful finite trace which extracts the coefficient at the identity, and that this trace is invariant under the reflection action of the second in . Furthermore, has (topological) stable rank 1 in the sense that invertibles are dense in , cf. Prop. 1.7 of [43]. So Rieffel’s Thm. 10.8 in [43] applies, saying that every stably free f.g.p. -module is actually already free.
This nice property of allows us to invoke Thm. 1.1, to answer the good Wannier basis existence problem in the affirmative:
Corollary 7.2**.**
*Let the crystallographic group act on the Euclidean plane as above, and let be a Hamiltonian satisfying the conditions of Thm. 1.1. Then for compact separated part , the module from Thm. 6.1 is a free -module (of rank say), so that there exists a good Wannier basis , , , for . *
To our knowledge, our technique is the first one that can achieve a result such as Corollary 7.2 for nonabelian .
Case . Maybe the simplest non-abelian crystallographic group is , with acting on as . The group acts on with the (lift of the) generator of effecting reflection about some origin. We may compute that
[TABLE]
where denotes the representation ring of , and is an extra projective -module which we will construct later (Fig. 5). This computation should be contrasted with
[TABLE]
The part is generated by f.g.p. -modules induced from irreducible -representations as follows. There is a natural action of on . Let be a finite-dimensional representation of and endow with the diagonal action. Then is a “free -module” in the equivariant -theory sense described in §11.2 of [5], and f.g.p. -modules are direct summands of such . There is a natural way to turn f.g.p. -modules into f.g.p. -modules, giving a correspondence of the equivariant and , see §11.7 of [5]. If we let be the regular representation of (which contains one copy each of the trivial and sign representations), then the above construction recovers as the basic free -module, and exhibits it as a direct sum of and , see Fig. 4. There is another summand of , illustrated in Fig. 5.
Alternatively, in the Fourier transformed picture, where has the dual “flip” action of character-conjugation (momentum-reversal in physics). Then one computes , with generators explicitly realised by equivariant line bundles, cf. Lemma 4.3 of [20], which are physically the Bloch bundles corresponding to , and .
7.4 Magnetic translations and twisted group algebras
Let be contractible (vanishing first and second cohomology also suffices), and let be a -invariant closed 2-form on . For , might arise as the curvature form of a magnetic field perpendicular to a surface. Pick a 1-form (connection/vector potential) such that , and an origin . For , let be a function on satisfying , such that for all . For instance, take . Let denote multiplication by the phase function . Then the (left) magnetic translations satisfy , where satisfies the group 2-cocycle condition,
[TABLE]
Here, we may verify each is indeed independent of , and is just a unimodular number.
Thus these magnetic translations furnish a projective unitary representation of on with cocycle , and one can verify that they commute with and thus with the magnetic Laplacian . A -periodic potential may also be added.
On , the left regular representation can be twisted by , by taking , obtaining the left -regular representation. We also have the right -regular representation, , where denotes the dual cocycle to . Using the cocycle identity, one sees that these projective representations commute with each other. We can construct the twisted group algebra of finite linear combinations of , as well as the norm completion in . Example 3.3 may be modified accordingly to obtain, inside , a right Hilbert -module isomorphic to (cf. Proposition 3.4) via the map of Eq. (15). Up to this point, the discussion in this Subsection can also be found in §B of [23] with minor differences in conventions.
The good Wannier basis existence problem (Problem 1) carries over to this -twisted setting. Namely, one considers spectral subspaces of admissible Hamiltonians acting on which are -invariant, while the notion of Wannier basis just needs a modification of “translates” by “-translates.” A Wannier basis is thus some number of copies of the left -regular representation (realised inside ).
If has polynomial growth, we can also construct the Fréchet algebras of rapidly decaying sequences as in §2, since the relevant estimates are unaffected by the scalings. In particular, Lemma 2.2 and Proposition 2 carry over in an identical way. Then the abstract results of §4 and §5 follow, as does the Main theorem 6.1. Thus we deduce that the solution to Problem 1 in the -twisted setting is solved by checking whether
[TABLE]
is a free f.g.p. Hilbert -module.
Let us mention that for the Euclidean plane, a lattice of translations, and a constant magnetic field, the algebra is just the famous noncommutative torus, and there exist non-trivial projective modules over — the detection of such modules via -theory and noncommutative Chern characters is of great importance in the quantum Hall effect [4] modelled typically by magnetic Schrödinger operators. Furthermore, the Baum–Connes conjectures also admit twisted versions, facilitating the computation of , in particular whether any non stably-free projective modules exist or not. If a spectral subspace (e.g. of a quantum Hall Hamiltonian) has being such a non-free module, then it cannot admit a good Wannier basis. Our analysis corroborates the physics heuristic that magnetic fields are usually responsible for lack of localization of the Wannier basis functions.
7.5 A non-Euclidean example
Consider the Heisenberg Lie group manifold
[TABLE]
Its Lie algebra elements, i.e. tangent vectors at the identity, extend to left invariant vector fields on . An inner product on the Lie algebra similary extends to a left-invariant Riemannian metric on . Thus is topologically , but geometrically very different. Physically, such a non-Euclidean geometry could model a uniform density of screw dislocations along the -direction [30, 25]. By restricting to , one obtains the the discrete subgroup , which of course acts on isometrically by translations. Note that there is a central extension
[TABLE]
and decomposing over the character space of the central subgroup furnishes as a continuous field of noncommutative tori. It is known that , see e.g. [1], with one of the generators being the class of the trivial projection, and that the -theory class of a general projection may be computed via parings with cyclic cocycles [25]. So we may compute in principle whether a f.g.p. -module is stably free or not.
Proposition 7.3**.**
*If a f.g.p. module over is stably free, it is even a free module. *
Proof*.*
We require, for a -algebra , the notions of general stable rank, gsr(), and connected stable rank, csr(), as defined in [43]. They are related by gsr()csr(), see pp. 328 of [43], and the computation csr() was carried out in [44], so we have gsr( or . As in the proof of Thm. 10.8 of [43], to show that gsr(, it suffices to observe that has a faithful finite trace. Then by Corollary 10.7 of [43], gsr() implies the claim of this Proposition.
Note that is nilpotent, thus of polynomial growth, and acts freely and isometrically on with quotient/fundamental domain a compact nilmanifold. Thus we may apply our Main Thm. 1.1 to conclude that when a spectral subspace of a -invariant Hamiltonian (satisfying the generic conditions of that theorem) has with -theory class , we can already conclude that a good Wannier basis , , for exists.
Acknowledgements. We would like to thank Vito Zenobi, Varghese Mathai, Giuseppe De Nittis, Gianluca Panati, and Domenico Monaco for helpful discussion. The first-named author was supported by Australian Research Council Discovery Project grant FL170100020, under Chief Investigator and Australian Laureate Fellow Mathai Varghese, and the second-named author by Australian Research Council Discovery Early Career Researcher Award grant DE170100149 and Discovery Projects DP200100729.
Data availability statement. Data sharing is not applicable to this article as no new data were created or analyzed in this study.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] J. Anderson, W. Paschke.: The rotation algebra. Houston J. Math. 15 (1) 1–26 (1989)
- 2[2] C. Bär, N. Ginoux, F. Pfäeffle.: Wave Equations on Lorentzian Manifolds and Quantization (ESI Lectures in Mathematics and Physics), EMS Publishing House, 2007.
- 3[3] P. Baum, A. Connes, N. Higson.: Classifying space for proper actions and K 𝐾 K -theory of group C ∗ superscript 𝐶 C^{*} -algebras. C ∗ superscript 𝐶 C^{*} -algebras: 1943–1993 (San Antonio, TX, 1993), Contemp. Math. 167 , Amer. Math. Soc. (1994), 240–291.
- 4[4] J. Bellissard.: K 𝐾 K -theory of C ∗ superscript 𝐶 C^{*} -algebras in solid state physics. Statistical mechanics and field theory: mathematical aspects, Springer, pp. 99–156, 1986
- 5[5] B. Blackadar.: K 𝐾 K -theory for operator algebras. Math. Sci. Res. Inst. Publ. vol. 5, Cambridge University Press, 1998.
- 6[6] C. Bourne, A. Carey, A. Rennie.: The bulk-edge correspondence for the quantum Hall effect in Kasparov theory. Lett. Math. Phys. 105 (9) 1253–1273 (2015)
- 7[7] B. Bradlyn, L. Elcoro, J. Cano, M.G. Vergniory, Z. Wang, C. Felser, M.I. Aroyo, B.A. Bernevig: Topological quantum chemistry. Nature 547 (7663) 298 (2017)
- 8[8] C. Brouder, G. Panati, M. Calandra, C. Mourougane, N. Marzari.: Exponential Localization of Wannier Functions in Insulators. Phys. Rev. Lett. 98 046402 (2007)
