This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Kähler geometry and Chern insulators: Relations between topology and the quantum metric

Bruno Mera Instituto de Telecomunicações, 1049-001 Lisboa, Portugal Departmento de Física, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001 Lisboa, Portugal Departmento de Matemática, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001 Lisboa, Portugal    Tomoki Ozawa Advanced Institute for Materials Research (WPI-AIMR), Tohoku University, Sendai 980-8577, Japan
Abstract

We study Chern insulators from the point of view of Kähler geometry, i.e. the geometry of smooth manifolds equipped with a compatible triple consisting of a symplectic form, an integrable almost complex structure and a Riemannian metric. The Fermi projector, i.e. the projector onto the occupied bands, provides a map to a Kähler manifold. The quantum metric and Berry curvature of the occupied bands are then related to the Riemannian metric and symplectic form, respectively, on the target space of quantum states. We find that the minimal volume of a parameter space with respect to the quantum metric is π|𝒞|\pi|\mathcal{C}|, where 𝒞\mathcal{C} is the first Chern number. We determine the conditions under which the minimal volume is achieved both for the Brillouin zone and the twist-angle space. The minimal volume of the Brillouin zone, provided the quantum metric is everywhere non-degenerate, is achieved when the latter is endowed with the structure of a Kähler manifold inherited from the one of the space of quantum states. If the quantum volume of the twist-angle torus is minimal, then both parameter spaces have the structure of a Kähler manifold inherited from the space of quantum states. These conditions turn out to be related to the stability of fractional Chern insulators. For two-band systems, the volume of the Brillouin zone is naturally minimal provided the Berry curvature is everywhere non-negative or nonpositive, and we additionally show how the latter, which in this case is proportional to the quantum volume form, necessarily has zeros due to topological constraints.

I Introduction

The notion of topological phases has drastically changed our understanding of gapped phases of matter. There is much to learn about a band insulator beyond the assertion that it has a gap separating the valence bands, which form the occupied bands at zero temperature, from the conduction bands. In the particular case of two spatial dimensions, in the absence of other symmetries, the occupied bands may have a non-trivial topological “twist” determining what is called a Chern insulator. This topological twist is not merely a mathematical observation and it is manifested in the topological response of the system as known in integer and fractional quantum Hall effects, in which the transverse Hall conductivity is proportional to the first Chern number, a topological invariant of the vector bundle of occupied states over the two-dimensional Brillouin zone. Finer properties of insulating states refer not just to the topology, but also to their geometry [1, 2]. Of particular interest is the momentum-space quantum metric, related to the overlap of Bloch states at nearby quasi-momenta [1, 2, 3]. As it turns out, the integral of this metric over the Brillouin zone [4] is associated to the spread functional and the localization tensor of the material [5, 6, 7, 8, 3]. There are numerous proposals in the literature for extracting the quantum metric [9, 10, 11, 12, 13, 14, 15]. A proposal [4], making use of the fluctuation-dissipation theorem, shows that the localization tensor can also be measured through spectroscopy in synthetic quantum matter, such as ultracold atomic gases or trapped ions. The momentum-space quantum metric has recently been observed in ultracold gases [16] and in a microcavity hosting exciton-polariton modes [17].

The connection of microscopic geometric quantities, such as the Brillouin zone Berry curvature and quantum metric, to macroscopic geometric quantities associated with ground state properties of the band insulator is given, at zero temperature, by integrating the former quantities over the Brillouin zone. Equivalently, one can consider the system in finite size and take periodic boundary conditions twisted by phases. As one varies these phases, one obtains a family of many-particle ground states over the so-called twist-angle space and from which we can extract the macroscopic geometric quantities mentioned above. In two spatial dimensions, the Berry curvature in twist-angle space yields the famous Thouless–Kohomoto–Nightingale–den Nijs result for the Hall conductivity [18, 19, 20] and the quantum metric is related to the localization tensor [4, 20]. The anisotropy of the localization tensor is measured by a modular parameter τ\tau in the upper half-plane which describes a complex structure in twist-angle space [20].

The study of the geometry of band insulators can also be used to understand if the material can host stable fractional topological phases [21, 22]. In particular, for a Chern band to have an algebra of projected density operators which is isomorphic to the WW_{\infty} algebra found by Girvin, MacDonald, and Platzman for the fractional quantum Hall effect [23], certain geometric constraints are naturally found which enforce a compatibility relation between the momentum-space quantum metric and the Berry curvature of the band [21, 22]. This compatibility relation between the quantum metric and the Berry curvature is, as noted in Ref. [22], the same compatibility required, in an oriented surface, between the Riemannian metric and a symplectic form for these to endow the surface with a Kähler structure.

Motivated by the above, we provide in this paper a detailed study of the geometry of Chern insulators from the perspective of Kähler geometry. In particular, we find that the minimal quantum volume, as determined by the Riemannian volume form induced by the quantum metric, is given by π|𝒞|\pi|\mathcal{C}|, where 𝒞\mathcal{C} is the first Chern number of the occupied band. Furthermore, we determine the conditions of achieving the minimal quantum volume both in the Brillouin zone and in the twist-angle space. The minimal volume of the Brillouin zone, provided the quantum metric is everywhere non-degenerate, is achieved when the latter is endowed with the structure of a Kähler manifold induced from the one of the space of quantum states. If the quantum volume of the twist-angle torus is minimal, then both tori have the structure of a Kähler manifold inherited from the space of quantum states. These conditions are related to the geometric stability conditions for fractional Chern insulators presented in Refs. [21, 22], measuring the deviation from the lowest Landau-level physics. For two-band systems, the quantum volumes are naturally minimal provided the Berry curvature is everywhere non-negative or nonpositive, and we additionally show how the latter, which in this case is proportional to the quantum volume form, necessarily has zeros due to topological constraints. The results of this paper are accompanied, complemented, and supported by those of Ref. [24], in which physical implications of the mathematical structures we find in this paper are studied and verified with explicit models.

This paper is organized as follows. In Sec. II, we explore the geometry of Chern insulators from the point of view of Kähler geometry. In Sec. III, we present our main results, namely Theorems 12 and 3, together with some corollaries and remarks. In Sec. IV, we provide a physical interpretation for the results of Theorem 1 in terms of the many-particle ground state of the band insulator. Finally in Sec. V we draw the conclusions. In Appendix A, the reader can find detailed mathematical derivations of the results presented in Sec. III.

II The geometry of Chern insulators

In a tight-binding description of an insulator in two spatial dimensions, in the presence of translation invariance, the Hamiltonian is described by an n×nn\times n matrix H(k)H(\textnormal{{k}}) which is smooth as a function of quasimomentum k in the Brillouin zone BZ2\textnormal{\text{BZ}}^{2}, topologically, a two-torus. The number nn specifies the internal degrees of freedom. The presence of a gap below the Fermi level EFE_{F} allows one to define the Fermi projector P(k)=Θ(EFH(k))P(\textnormal{{k}})=\Theta(E_{F}-H(\textnormal{{k}})), where Θ()\Theta(\cdot) is the Heaviside step function, which describes the ground state of the insulator and, up to adiabatic deformation, it determines its topological properties. In fact, the smoothness and the gap condition imply that we have a rank rr Hermitian vector bundle EBZ2E\to\textnormal{\text{BZ}}^{2}, where rr is the number of occupied bands and is the constant rank of the Fermi projector, whose fiber over kBZ2\textnormal{{k}}\in\textnormal{\text{BZ}}^{2} is simply Im[P(k)]n\text{Im}\left[P(\textnormal{{k}})\right]\subset\mathbb{C}^{n}. The smooth family of orthogonal projectors {P(k)}kBZ2\{P(\textnormal{{k}})\}_{\textnormal{{k}}\in\textnormal{\text{BZ}}^{2}} can be seen as a classifying map P:BZ2Grk(n)P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{k}(\mathbb{C}^{n}) for the vector bundle EBZ2E\to\textnormal{\text{BZ}}^{2}, where we interpret the Grassmannian manifold of rr-dimensional subspaces of n\mathbb{C}^{n}, denoted Grr(n)\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}), as the set of orthogonal projectors of rank rr. The Chern insulator is determined by the reduced complex K-theory class [E]K~0(BZ2)[E]\in\widetilde{K}^{0}(\textnormal{\text{BZ}}^{2}), which is determined by the idempotent PP [25, 26, 27].

The space Grr(n)\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) has the natural structure of a Kähler manifold (see, for instance, Refs. [28, 29] for more details), i.e., it is a symplectic manifold together with a compatible complex structure JJ giving it the structure of a complex manifold. By a symplectic manifold structure we mean that it is a smooth manifold together with a non-degenerate closed two-form ω\omega – the symplectic form. The almost complex structure JJ is a linear map on the tangent spaces that squares to minus the identity. Morally, it acts as multiplication by 1\sqrt{-1}. The fact that Grr(n)\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) is a complex manifold means that we can find local complex coordinates zjz^{j} and, with respect to these,

J(zj)=izj and J(z¯j)=iz¯j.\displaystyle J\left(\frac{\partial}{\partial z^{j}}\right)=i\frac{\partial}{\partial z^{j}}\text{ and }J\left(\frac{\partial}{\partial\overline{z}^{j}}\right)=-i\frac{\partial}{\partial\overline{z}^{j}}. (1)

Finally, the fact that ω\omega and JJ are compatible means that g=ω(,J)g=\omega(\cdot,J\cdot) is a Riemannian metric. It is said in this case that the complex structure JJ is ω\omega-compatible. In a Kähler manifold, the symplectic form is also known as the Kähler form. In the particular case of the Grassmannian, the compatible triple of structures stem from the usual Fubini-Study Kähler structure on the complex projective space Pm\mathbb{C}P^{m}, with m=(nr)1m=\binom{n}{r}-1, through the Plücker embedding ι:Grr(n)(Λrn)Pm\iota:\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n})\hookrightarrow\mathbb{P}(\Lambda^{r}\mathbb{C}^{n})\cong\mathbb{C}P^{m} (see Ref. [30]). The latter map sends Im(P)\text{Im}(P), or equivalently, PP to the line spanned by a Slater determinant made of a linear basis for Im(P)\text{Im}(P). For this reason we write these structures as ωFS,JFS\omega_{FS},J_{FS}, and gFSg_{FS}.

One can write these three structures in a familiar form by using orthogonal projectors. To do this, it is useful to recall that, since a projector satisfies P2=PP^{2}=P, if Q=IPQ=I-P is the orthogonal complement projector, where II is the identity matrix, then we have the relation PdPP=QdPQ=0PdPP=QdPQ=0 and hence the differential dPdP satisfies dP=QdPP+PdPQdP=QdPP+PdPQ. Furthermore, since P=PP^{\dagger}=P for orthogonal projectors, we see that dPdP is completely determined by the triangular matrix QdPPQdPP [in an orthonormal basis where P=diag(Ir,0nr)P=\text{diag}(I_{r},0_{n-r}) it is represented by a lower triangular matrix whose entries are one-forms].

The symplectic and Riemannian structures are given by

ωFS\displaystyle\omega_{FS} =i2Tr((QdPP)QdPP)=i2Tr(PdPdP),\displaystyle=-\frac{i}{2}\textnormal{\text{Tr}}\left(\left(QdPP\right)^{\dagger}\wedge QdPP\right)=-\frac{i}{2}\textnormal{\text{Tr}}\left(PdP\wedge dP\right),
gFS\displaystyle g_{FS} =Tr((QdPP)QdPP)=Tr(PdPdP),\displaystyle=\textnormal{\text{Tr}}\left(\left(QdPP\right)^{\dagger}QdPP\right)=\textnormal{\text{Tr}}\left(PdPdP\right), (2)

and the complex structure is defined by

(QdPP)JFS=iQdPP and (PdPQ)JFS=iPdPQ,\displaystyle\left(QdPP\right)\circ J_{FS}\!=\!iQdPP\!\text{ and }\!\left(PdPQ\right)\circ J_{FS}\!=\!-iPdPQ, (3)

where the second equation is essentially the Hermitian conjugate of the first since JFSJ_{FS} is an \mathbb{R}-linear operator and, therefore, it is completely specified by the first equation. It is not hard to see that, indeed, ωFS(,JFS)=gFS\omega_{FS}(\cdot,J_{FS}\cdot)=g_{FS}. A remarkable property of this structure is that it is invariant under the action of U(n)\text{U}(n) by conjugation PUPUP\mapsto UPU^{\dagger}, for UU(n)U\in\text{U}(n), so all the structures are completely determined by knowing them at one particular projector. Before moving on, it is important to note that ωFS\omega_{FS} is related to the curvature two-form of a connection on a line bundle over Grr(n)\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}). Namely, over Grr(n)\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) we can consider the vector bundle whose fiber over an orthogonal projector PP is simply Im(P)\text{Im}(P) — the tautological bundle E0Grr(n)E_{0}\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}). This rank rr-vector bundle is naturally a subbundle of the trivial bundle θn=Grr(n)×n\theta^{n}=\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n})\times\mathbb{C}^{n}. The Berry connection is then simply given by projection of the exterior derivative. The curvature of this connection is determined by the matrix of two-forms

Θ=iPdPdPP.\displaystyle\Theta=iPdP\wedge dPP. (4)

We can then consider the top exterior power of this bundle, ΛrE0Grr(n)\Lambda^{r}E_{0}\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}), which is a line bundle whose fiber at PP is spanned by a Slater determinant

v1vr, with span{vi}i=1r=Im(P).\displaystyle v_{1}\wedge\dots\wedge v_{r},\text{ with }\text{span}_{\mathbb{C}}\{v_{i}\}_{i=1}^{r}=\text{Im}(P). (5)

The connection induced on this line bundle from the Berry connection on the tautological vector bundle has curvature two-form

Ω0=Tr(Θ)=iTr(PdPdPP)=iTr(PdPdP).\displaystyle\Omega_{0}=\textnormal{\text{Tr}}\left(\Theta\right)=i\textnormal{\text{Tr}}\left(PdP\wedge dPP\right)=i\textnormal{\text{Tr}}\left(PdP\wedge dP\right). (6)

The relation between ωFS\omega_{FS} and Ω0\Omega_{0} is then clear: ωFS=(1/2)Ω0\omega_{FS}=-(1/2)\Omega_{0}. In the physics literature, the quantity χ=gFS+iωFS=gFSi(Ω0/2)\chi=g_{FS}+i\omega_{FS}=g_{FS}-i(\Omega_{0}/2) is usually known as the quantum geometric tensor.

Geometric quantities related to gFSg_{FS} and ωFS\omega_{FS} appear naturally in the context of insulators because we have a map P:BZ2Grr(n)P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}), which provides the means to pull back the rich geometry in the target to the Brillouin zone. The pullback by PP of the tautological vector bundle, i.e., the vector bundle whose fiber at kBZ2\textnormal{{k}}\in\textnormal{\text{BZ}}^{2} is Im[P(k)]\text{Im}\left[P(\textnormal{{k}})\right], is simply the occupied band vector bundle PE0=EBZ2P^{*}E_{0}=E\to\textnormal{\text{BZ}}^{2}. The pullback of the rrth exterior power of the tautological vector bundle is a line bundle whose fiber at k is spanned by a Slater determinant of a basis for Im[P(k)]\text{Im}\left[P(\textnormal{{k}})\right]. This bundle specifies the ground state completely as, physically, the ground state is obtained precisely by occupying the rr bands below the Fermi level EFE_{F}.

The Berry curvature for ΛrE=PΛrE0\Lambda^{r}E=P^{*}\Lambda^{r}E_{0} is

Ω=PΩ0=ii,j=12Tr(P(k)Pki(k)Pkj(k))dkidkj.\displaystyle\Omega=P^{*}\Omega_{0}=i\sum_{i,j=1}^{2}\textnormal{\text{Tr}}\left(P(\textnormal{{k}})\frac{\partial P}{\partial k_{i}}(\textnormal{{k}})\frac{\partial P}{\partial k_{j}}(\textnormal{{k}})\right)dk_{i}\wedge dk_{j}. (7)

Meanwhile, the pullback of the metric gFSg_{FS} by PP, is the so-called quantum metric

g=PgFS=i,j=12Tr(P(k)Pki(k)Pkj(k))dkidkj.\displaystyle g=P^{*}g_{FS}=\sum_{i,j=1}^{2}\textnormal{\text{Tr}}\left(P(\textnormal{{k}})\frac{\partial P}{\partial k_{i}}(\textnormal{{k}})\frac{\partial P}{\partial k_{j}}(\textnormal{{k}})\right)dk_{i}dk_{j}. (8)

Following the discussion in Ref. [20], we can consider threading a flux through a finite system of size N×NN\times N with periodic boundary conditions, i.e., in a torus. This means that the fermions acquire phases eiθie^{i\theta_{i}}, i=1,2i=1,2, as they are adiabatically moved around the fundamental cycles of position space. The angles 𝜽=(θ1,θ2)\bm{\theta}=(\theta_{1},\theta_{2}) are referred to as the twist-angles and they span the twist-angle torus Tθ2T^{2}_{\theta}. In momentum space, these angles can be accounted for by sampling the matrix H(k)H(\textnormal{{k}}) at points k=(2π/N)m+𝜽/N\textnormal{{k}}=(2\pi/N)\textnormal{{m}}+\bm{\theta}/N, with m{0,,N1}2\textnormal{{m}}\in\{0,\dots,N-1\}^{2}. The ground state of the system, as we vary 𝜽\bm{\theta}, produces a smooth Hermitian line bundle Tθ2\mathcal{L}\to T^{2}_{\theta}, whose Berry curvature has the form

Ω~=1N2mΩ12(2πmN+𝜽N)dθ1dθ2,\displaystyle\widetilde{\Omega}=\frac{1}{N^{2}}\sum_{\textnormal{{m}}}\Omega_{12}\left(\frac{2\pi\textnormal{{m}}}{N}+\frac{\bm{\theta}}{N}\right)d\theta_{1}\wedge d\theta_{2}, (9)

where Ω12=Ω(/k1,/k2)\Omega_{12}=\Omega(\partial/\partial k_{1},\partial/\partial k_{2}) is the single component of the Berry curvature in Eq. (7). In the thermodynamic limit, NN\to\infty, we get

Ω~(BZ2d2k(2π)2Ω12(k))dθ1dθ2=𝒞2πdθ1dθ2,\displaystyle\widetilde{\Omega}\to\left(\int_{\textnormal{\text{BZ}}^{2}}\frac{d^{2}k}{(2\pi)^{2}}\Omega_{12}(\textnormal{{k}})\right)d\theta_{1}\wedge d\theta_{2}=\frac{\mathcal{C}}{2\pi}d\theta_{1}\wedge d\theta_{2}, (10)

where 𝒞=BZ2Ω/(2π)\mathcal{C}=\int_{\textnormal{\text{BZ}}^{2}}\Omega/(2\pi) is the first Chern number of the occupied Bloch bundle EBZ2E\to\textnormal{\text{BZ}}^{2}. The relation to linear response theory is then

e2Ω~=σHalldθ1dθ2,\displaystyle e^{2}\widetilde{\Omega}=\sigma_{\text{Hall}}d\theta_{1}\wedge d\theta_{2}, (11)

where σHall=e2𝒞/(2π)\sigma_{\text{Hall}}=e^{2}\mathcal{C}/(2\pi) is the Hall conductivity of the insulator and ee is the charge of the fermions. We find that the Berry curvature is, in the thermodynamic limit NN\to\infty, flat. Observe additionally that

𝒞=BZ2Ω2π=Tθ2Ω~2π,\displaystyle\mathcal{C}=\int_{\textnormal{\text{BZ}}^{2}}\frac{\Omega}{2\pi}=\int_{T^{2}_{\theta}}\frac{\widetilde{\Omega}}{2\pi}, (12)

and, as a consequence, the line bundle Tθ2\mathcal{L}\to T^{2}_{\theta} is isomorphic, in the smooth category, to the line bundle ΛrEBZ2\Lambda^{r}E\to\textnormal{\text{BZ}}^{2}, where we identify BZ2Tθ2\textnormal{\text{BZ}}^{2}\cong T^{2}_{\theta}, through k=(k1,k2)(θ1(k),θ2(k))=(k1,k2)\textnormal{{k}}=(k_{1},k_{2})\mapsto(\theta^{1}(\textnormal{{k}}),\theta^{2}(\textnormal{{k}}))=(k^{1},k^{2}). In a similar fashion, we can determine the quantum metric in twist-angle space, yielding,

g~=1N2i,j=12mgij(2πmN+𝜽N)dθidθj,\displaystyle\widetilde{g}=\frac{1}{N^{2}}\sum_{i,j=1}^{2}\sum_{\textnormal{{m}}}g_{ij}\left(\frac{2\pi\textnormal{{m}}}{N}+\frac{\bm{\theta}}{N}\right)d\theta_{i}d\theta_{j}, (13)

where gij=g(/ki,/kj)g_{ij}=g(\partial/\partial k_{i},\partial/\partial k_{j}) are the components of the quantum metric in Eq. (8). In the same way that the Berry curvature is flat, also the quantum metric is flat in the thermodynamic limit NN\to\infty

g~=i,j=12(BZ2d2k(2π)2gij(k))dθidθj.\displaystyle\widetilde{g}=\sum_{i,j=1}^{2}\left(\int_{\textnormal{\text{BZ}}^{2}}\frac{d^{2}k}{(2\pi)^{2}}g_{ij}(\textnormal{{k}})\right)d\theta_{i}d\theta_{j}. (14)

This metric also appears in the properties of the insulator in question, namely, it is related to the so-called localization tensor. The fact that we have a flat metric tensor, allowed us in Ref. [20], to extract a modular parameter τ\tau\in\mathcal{H}, which determines a complex structure j~\widetilde{j} on Tθ2T^{2}_{\theta}. This complex structure is defined by a rotation of 90 degrees in the standard orientation of the torus with respect to the quantum metric g~\widetilde{g}. This complex structure turned out to be related to the anisotropy of the localization tensor.

In two dimensions having a metric and an orientation specifies a complex structure. Indeed, provided gg is non-degenerate, we get, by a rotation of 90 degrees in the standard orientation of BZ2\textnormal{\text{BZ}}^{2} with respect to gg, an almost complex structure over it which we call jj. We can get an explicit expression for jj as follows. At a given point k of the Brillouin zone we can always choose an orthonormal frame e1,e2e_{1},e_{2}, consistent with the orientation (i.e., the determinant of the matrix relating e1,e2e_{1},e_{2} to /k1,/k2\partial/\partial k_{1},\partial/\partial k_{2} is positive), for the quantum metric gg, which is a Riemannian metric, so that

g(ei,ej)=δij,i,j=1,2.\displaystyle g(e_{i},e_{j})=\delta_{ij},\ i,j=1,2.

Then at that point j(e1)=e2j(e_{1})=e_{2} and j(e2)=e1j(e_{2})=-e_{1}. Hence, it follows that

g(j(ei),ej)=εij,i,j=1,2,\displaystyle g\left(j(e_{i}),e_{j}\right)=\varepsilon_{ij},\ i,j=1,2, (15)

where εij\varepsilon_{ij} is the Levi-Civita symbol in two dimensions. Equivalently, g(j,)=e1e2g(j\cdot,\cdot)=e^{1}\wedge e^{2} is the volume form associated with gg, where e1,e2e^{1},e^{2} is the dual basis of e1,e2e_{1},e_{2}. We will now give an expression for jj in the original (k1,k2)(k_{1},k_{2}) coordinates. Let us denote by gijg_{ij} and jjij^{i}_{j}, 1i,j21\leq i,j\leq 2, the components of gg and jj in the (k1,k2)(k_{1},k_{2}) coordinates, respectively. From g(j,)=e1e2=det(g)dk1dk2g(j\cdot,\cdot)=e^{1}\wedge e^{2}=\sqrt{\det(g)}dk_{1}\wedge dk_{2}, where det(g)=det(gij)\det(g)=\det(g_{ij}), we get,

k=12gkjjik=det(g)εij,i,j=1,2.\displaystyle\sum_{k=1}^{2}g_{kj}j^{k}_{i}=\sqrt{\det(g)}\varepsilon_{ij},\ i,j=1,2. (16)

It is also useful to write the equivalent form

gij=k=12det(g)εikjjk,i,j=1,2.\displaystyle g_{ij}=\sum_{k=1}^{2}\sqrt{\det(g)}\varepsilon_{ik}j^{k}_{j},\ i,j=1,2. (17)

Either way, it follows that

jji=1det(g)k=12εikgkj,i,j=1,2,\displaystyle j^{i}_{j}=-\frac{1}{\sqrt{\det(g)}}\sum_{k=1}^{2}\varepsilon^{ik}g_{kj},\ i,j=1,2, (18)

or, in matrix form,

j=1det(g)[g12g22g11g12].\displaystyle j=\frac{1}{\sqrt{\det(g)}}\left[\begin{array}[]{cc}-g_{12}&-g_{22}\\ g_{11}&g_{12}\end{array}\right]. (21)

By construction, jj does indeed satisfy the requirements of an almost complex structure, namely that j2=Ij^{2}=-I, provided gg is non-degenerate. The expression presented above for jj is global, i.e., we can use it at any point over the Brillouin zone. However, being able to find local coordinates x1,x2x_{1},x_{2} in a neighborhood of any point kBZ2\textnormal{{k}}\in\textnormal{\text{BZ}}^{2} such that jj assumes the canonical form j(/x1)=/x2j(\partial/\partial x_{1})=\partial/\partial x_{2} and j(/x2)=/x1j(\partial/\partial x_{2})=-\partial/\partial x_{1} is a much more subtle question, namely that of integrability of jj, and it implies solving a partial differential equation. If such local coordinates exist, then z=x1+ix2z=x_{1}+ix_{2} will be a local complex coordinate since j(/z)=i/zj(\partial/\partial z)=i\partial/\partial z and j(/z¯)=i/z¯j(\partial/\partial\overline{z})=-i\partial/\partial\overline{z}. What we did above by choosing a orthonormal basis at a point (which is a linear combination of the natural tangent vectors /k1\partial/\partial k_{1} and /k2\partial/\partial k_{2}) was to make jj look canonical at a specific point. Indeed we can always choose local coordinates xi(k1,k2)x_{i}(k_{1},k_{2}), i=1,2i=1,2, such that /xi\partial/\partial x_{i}, i=1,2i=1,2, become, at a specific point, an orthonormal frame for gg but away from that point that will, in general, not happen. Finding local coordinates such that this happens everywhere, up to a conformal factor, is equivalent to solving the Beltrami equation for isothermal coordinates. Hence, solving the Beltrami equation gives us local holomorphic coordinates in two dimensions. Being able to find local holomorphic coordinates for which jj assumes the canonical form of multiplication by ii is a very special feature of oriented surfaces. Indeed, in two dimensions any almost complex structure jj is integrable and it endows the surface with the structure of a complex manifold [31]. One is then naturally lead to ask under which conditions the geometry induced by PP on BZ2\textnormal{\text{BZ}}^{2}, given by the triple (PωFS,j,PgFS)=(ω,j,g)(P^{*}\omega_{FS},j,P^{*}g_{FS})=(\omega,j,g), is compatible and gives rise to a Kähler structure. In this case, the map P:BZ2Grr(n)P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) is a Kähler map in the sense that the pullback of the triple (ωFS,JFS,gFS)(\omega_{FS},J_{FS},g_{FS}) is well-defined as a compatible triple (ω=PωFS,j=PJFS,g=PgFS)(\omega=P^{*}\omega_{FS},j=P^{*}J_{FS},g=P^{*}g_{FS}) giving BZ2\textnormal{\text{BZ}}^{2} the structure of a Kähler manifold. The meaning of j=PJFSj=P^{*}J_{FS} needs further justification. In the case that g=PgFSg=P^{*}g_{FS} is a Riemannian metric the map PP has to be an immersion, i.e., the differential dPdP is full-rank at every point of the Brillouin zone. In that case dPdP provides an isomorphism of the tangent spaces of the Brillouin zone to their images. The condition that ω(,j)\omega(\cdot,j\cdot), with ω\omega non-degenerate, then turns out to imply dPj=JFSdPdP\circ j=J_{FS}\circ dP, i.e., the map is holomorphic — this is what we mean by j=PJFSj=P^{*}J_{FS}. We will see that when (ω,j,g)(\omega,j,g) form a compatible triple, the Chern number is necessarily negative. We point out that, under the condition of non-degeneracy of gg, the Riemannian volume form det(g)dk1dk2\sqrt{\det(g)}dk_{1}\wedge dk_{2} (notice that the orientation of BZ2\textnormal{\text{BZ}}^{2} is explicitly used to define it), where det(g)\det(g) is the determinant of the matrix representing the metric in periodic coordinates, is symplectic, i.e., it is non-degenerate and closed. As a consequence, if we consider the triple of structures (det(g)dk1dk2,j,g)(\sqrt{\det(g)}dk_{1}\wedge dk_{2},j,g), we easily see that it is compatible, cf. Eq. (17), and it gives BZ2\textnormal{\text{BZ}}^{2} a Kähler structure. However, this triple does not generally coincide with the triple (ω,j,g)(\omega,j,g). We will see that only when the map is Kähler these two triples coincide. One can then ask what are the physical consequences of this condition. Furthermore, we could also ask whether (ω~=Ω~/2,j~,g~)(\widetilde{\omega}=-\widetilde{\Omega}/2,\widetilde{j},\widetilde{g}) gives Tθ2T^{2}_{\theta} the structure of a Kähler manifold and what are the physical consequences.

As we will see, there will also be cases where (ω,j,g)(\omega,-j,g) form a compatible triple. Namely, the almost complex structure is defined with respect to the orientation opposite to the natural orientation specified by the ordered pair of coordinates (k1,k2)(k_{1},k_{2}). This situation corresponds to the case where the Chern number is positive.

We note that the relation between holomorphicity of the wave function, related to the holomorphicity of the map to the projective space (see Sec. IV), and the Berry curvature being a Kähler form has been studied in the different context of the fractional quantum Hall effect in curved background geometry [32], where the relevant metric there is that of the spatial sample as opposed to the quantum metric.

III Main results

In this section we answer the questions formulated in the end of previous section. We refer the reader to Appendix A for the proofs.

We first note that the condition of P:BZ2Grr(n)P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) to be an immersion is equivalent to g=PgFSg=P^{*}g_{FS} being non-degenerate [det(g)0\det(g)\neq 0 everywhere], which we use frequently in the discussion below.

Next, we observe that for P:BZ2Grr(n)P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) to be Kähler it is enough for PP to be an immersion and to be holomorphic, i.e., to satisfy

QPz¯=0,\displaystyle Q\frac{\partial P}{\partial\overline{z}}=0, (22)

in local complex coordinates, see Proposition 2 and its proof in Appendix A. We then have the following result:

Theorem 1.

Suppose P:BZ2Grr(n)P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) is an immersion, then the Cauchy-Schwarz like inequality is saturated

det(gij(k))=Ω12(k)2, for all kBZ2,\displaystyle\sqrt{\det(g_{ij}(\textnormal{{k}}))}=\frac{\mp\Omega_{12}(\textnormal{{k}})}{2},\text{ for all }\textnormal{{k}}\in\textnormal{\text{BZ}}^{2},

if and only if the map PP is Kähler with respect to the triple of structures (ω,±j,g)(\omega,\pm j,g).

We will use the notation volg(BZ2)\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2}) and volg~(Tθ2)\textnormal{\text{vol}}_{\widetilde{g}}(T^{2}_{\theta}) to denote the integrals

volg(BZ2)=BZ2det(g)d2k\displaystyle\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2})=\int_{\textnormal{\text{BZ}}^{2}}\sqrt{\det(g)}d^{2}k

and

volg~(Tθ2)=Tθ2det(g~)d2θ,\displaystyle\textnormal{\text{vol}}_{\widetilde{g}}(T^{2}_{\theta})=\int_{T^{2}_{\theta}}\sqrt{\det(\widetilde{g})}d^{2}\theta,

where det(g)\det(g) and det(g~)\det(\widetilde{g}) are the determinants of the matrices representing gg and g~\widetilde{g} in the periodic coordinates k=(k1,k2)\textnormal{{k}}=(k_{1},k_{2}) and 𝜽=(θ1,θ2)\bm{\theta}=(\theta_{1},\theta_{2}), respectively. In the case where the quantum metrics are non-degenerate the quantities volg(BZ2)\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2}) and volg~(Tθ2)\textnormal{\text{vol}}_{\widetilde{g}}(T^{2}_{\theta}) correspond to the volumes of BZ2\textnormal{\text{BZ}}^{2} and Tθ2T^{2}_{\theta}, respectively, as measured by the quantum metrics. For this reason, we will refer to them as the quantum volumes and the top-forms det(g)dk1dk2\sqrt{\det(g)}dk_{1}\wedge dk_{2} and det(g~)dθ1dθ2\sqrt{\det(\widetilde{g})}d\theta_{1}\wedge d\theta_{2} as the quantum volume-forms of the Brillouin zone and twist-angle space, respectively. Furthermore we will use the diffeomorphism BZ2Tθ2\textnormal{\text{BZ}}^{2}\cong T^{2}_{\theta} given by k=(k1,k2)(θ1(k),θ2(k))=(k1,k2)\textnormal{{k}}=(k_{1},k_{2})\mapsto(\theta^{1}(\textnormal{{k}}),\theta^{2}(\textnormal{{k}}))=(k^{1},k^{2}) to identify the two tori.

We have the following result and associated corollaries.

Theorem 2.

The following inequalities hold

π|𝒞|volg(BZ2)volg~(Tθ2),\displaystyle\pi|\mathcal{C}|\leq\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2})\leq\textnormal{\text{vol}}_{\widetilde{g}}(T^{2}_{\theta}),

where the left-hand side inequality is saturated if and only if det(g)=|Ω12|2\sqrt{\det(g)}=\frac{|\Omega_{12}|}{2} and provided Ω12\Omega_{12} does not change sign over BZ2\textnormal{\text{BZ}}^{2}, and the right-hand side inequality is satisfied if and only if, for every kBZ2\textnormal{{k}}\in\textnormal{\text{BZ}}^{2},

gij(k)=e2f(k)g~ij(θ=k), 1i,j2,\displaystyle g_{ij}(\textnormal{{k}})=e^{2f(\textnormal{{k}})}\widetilde{g}_{ij}(\theta=\textnormal{{k}}),\ 1\leq i,j\leq 2,

for some function fC(BZ2)f\in C^{\infty}(\textnormal{\text{BZ}}^{2}), i.e., gg is related to g~\widetilde{g} by a Weyl rescaling, implying that, provided PP is an immersion, they share the same complex structure jj (and hence the same modular parameter).

Observe that the topological invariant π|𝒞|\pi|\mathcal{C}| appears as the lower bound for both quantum volumes volg(BZ2)\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2}) and volg~(Tθ2)\textnormal{\text{vol}}_{\widetilde{g}}(T^{2}_{\theta}). In other words, the first Chern number of the occupied Bloch bundle determines the minimal quantum volume. As a consequence of Theorem 2, we have the following corollaries.

Corollary 1.

If PP is an immersion, and thus det(g)0\det(g)\neq 0 everywhere, then volg(BZ2)=π|𝒞|\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2})=\pi|\mathcal{C}| if and only if PP is Kähler by Theorem 1.

Corollary 2.

If g~\widetilde{g} is non-degenerate and volg~(Tθ2)=π|𝒞|\textnormal{\text{vol}}_{\widetilde{g}}(T^{2}_{\theta})=\pi|\mathcal{C}| then (12Ω~,g~,±j~)(-\frac{1}{2}\widetilde{\Omega},\widetilde{g},\pm\widetilde{j}) is a flat Kähler structure. The converse is also true. In both cases, we also have volg(BZ2)=π|𝒞|\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2})=\pi|\mathcal{C}|. If in the first implication we have, in addition, that PP is an immersion, implying gg in non-degenerate, then PP is a Kähler map for the same complex structure. This implies that the same Weyl rescaling that relates gg and g~\widetilde{g} in this case, relates also Ω\Omega and Ω~\widetilde{\Omega}, i.e.,

Ω=e2fΩ~.\displaystyle\Omega=e^{2f}\widetilde{\Omega}.

Because the integral of the curvatures yields 𝒞\mathcal{C}, we have BZ2d2k(2π)2e2f(k)=1\int_{\textnormal{\text{BZ}}^{2}}\frac{d^{2}k}{(2\pi)^{2}}e^{2f(\textnormal{{k}})}=1, which is consistent with the findings in the proof of Theorem 2.

Remark 1.

Corollaries 1 and 2 are, in the mathematics literature, presented as a corollary of Wirtinger’s inequality (see Proposition 3 and its proof in Appendix A) which states that every complex submanifold of a Kähler manifold is volume minimizing in its homology class.

Corollary 3.

If any of the volumes volg(BZ2)<π\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2})<\pi or volg~(Tθ2)<π\textnormal{\text{vol}}_{\widetilde{g}}(T^{2}_{\theta})<\pi we have 𝒞=0\mathcal{C}=0.

Corollary 3 is perhaps one of the most dramatic consequences in the sense that the volumes measured by the quantum metrics, which are purely geometric quantities, allow us to infer about non-trivial topological properties of the insulator.

The results presented place strong constraints in the geometry of Chern bands, which is relevant to the study of the stability of fractional topological insulators. Consider now, for the sake of simplicity of the following discussion, the case of a single band, i.e., r=1r=1. In the works of Refs. [21, 22], the stability of topological phases arising in fractionally filled Chern insulators is studied. The saturation of the inequality in Proposition 3 together with the Fubini-Study metric being uniform throughout the Brillouin zone is then found as a criterion for the algebra of projected operators to be isomorphic to the one in the fractional quantum Hall effect — the WW_{\infty} algebra found by Girvin, MacDonald and Platzman [23]. Observe that this condition corresponds to the saturation of both inequalities in Theorem 2 and, provided non degeneracy of the metric, it is the situation described in Corollary 2, supplemented with the additional condition that the conformal factor is trivial to have constancy of the metric in momentum space. Furthermore, in Ref. [22] the saturation of the inequality of Proposition 3 for every kBZ2\textnormal{{k}}\in\textnormal{\text{BZ}}^{2} was identified as having a Kähler structure in the Brillouin zone. The saturation of the inequality is not enough to have a Kähler structure as one needs the immersion condition or, equivalently, det(g)0\det(g)\neq 0 everywhere. We remark that the holomorphicity condition QP/z¯=0Q\partial P/\partial\bar{z}=0 found in Proposition 2 is the momentum space form of the real space equation Qα(xiy)Pα=0Q_{\alpha}(x^{\prime}-iy^{\prime})P_{\alpha}=0 found in Ref. [21], for the case when both inequalities in Theorem 2 are satisfied and the complex structure in BZ2\textnormal{\text{BZ}}^{2} is the same as the one in Tθ2T^{2}_{\theta} and completely determined by a single (multi-valued) complex coordinate of the form z=k1±τk2z=k_{1}\pm\tau k_{2}, for some modular parameter τ\tau\in\mathcal{H} and the sign takes into account a possible orientation flip with respect to the standard one in the Brillouin zone. In Ref. [33] and in related works [34, 35, 36], a Kähler structure where the associated complex structure is given by the flat isotropic one, i.e., determined by z=k1+τk2z=k_{1}+\tau k_{2}, with τ=i\tau=i, plays an important role in the construction of ideal fractional Chern insulator models.

Finally, we present a few results which we noticed as a consequence of analyzing the particular case of two-band models.

Theorem 3.

For two-dimensional two-band models, there must exist a point in the Brillouin zone where det(g)=0\det(g)=0. In other words, the map P:BZ2Gr1(2)=P1S2P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{1}(\mathbb{C}^{2})=\mathbb{C}P^{1}\cong S^{2} cannot be an immersion.

The next proposition is known in the literature and, as far as we know, the first to notice it were Yu-Quan Ma et al[37, 38].

Proposition 1.

For two-dimensional two-band models the inequality of Proposition 3 is saturated det(g)=|Ω12|/2\sqrt{\det(g)}=|\Omega_{12}|/2.

Looking at the proof of Proposition 3 in the Appendix, this last result follows from the fact that the orthogonal complement projector Q(k)=IP(k)Q(\textnormal{{k}})=I-P(\textnormal{{k}}) has rank 11 in this case and saturation of the inequality is automatic. As a consequence of this proposition we have the two following corollaries.

Corollary 4.

For two-dimensional two-band models, there must exist a point in the Brillouin zone where the Berry curvature vanishes.

We point out that in Ref. [33] it was noted that it is impossible to have a uniform Berry curvature for two-dimensional two-band models – a manifestation of Theorem 3 and Corollary 4.

Corollary 5.

For two-dimensional two-band models, if Ω12\Omega_{12} or Ω12-\Omega_{12} is everywhere nonnegative then volg(BZ2)=π|𝒞|\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2})=\pi|\mathcal{C}|.

This implies that the quantum volume volg(BZ2)\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2}) directly tells us about the topology of quantum states in this particular case.

IV Physical interpretation of the Kähler map condition

In this section, we provide a physical interpretation of the Kähler map condition appearing in Theorem 1, namely, that the Fermi projector P(k)=Θ(EFH(k))P(\textnormal{{k}})=\Theta(E_{F}-H(\textnormal{{k}})), besides providing an immersion of the Brillouin zone BZ2\textnormal{\text{BZ}}^{2} to the Grassmannian Grr(n)\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}), satisfies

QPz¯=0,\displaystyle Q\frac{\partial P}{\partial\overline{z}}=0, (23)

in local complex coordinates zz coming from the complex structure jj (or j-j, depending on the sign of the first Chern number, see the proof of Theorem 1 in Appendix A) in the Brillouin zone. We will phrase this condition in terms of the many-particle ground state of the band insulator, which is obtained by filling the Bloch bands below the Fermi level. Recall that this state is, at the formal level, a Slater determinant over the momenta k in the Brillouin zone, where each momentum piece is itself a Slater determinant determined by the rr-dimensional subspace Im[P(k)]n\text{Im}\left[P(\textnormal{{k}})\right]\subset\mathbb{C}^{n}, where nn is the total number of bands in the system [see Eq. (5)]. We may construct a linear independent set of local Bloch wave functions |v1,k,.,|vr,k|v_{1,\textnormal{{k}}}\rangle,....,|v_{r,\textnormal{{k}}}\rangle spanning Im[P(k)]n\text{Im}\left[P(\textnormal{{k}})\right]\subset\mathbb{C}^{n}, depending smoothly on k, for k in a small open set, which we can write in terms of the canonical basis of n\mathbb{C}^{n},

|vi,k=j=1nvij(k)|j,i=1,,r.\displaystyle|v_{i,\textnormal{{k}}}\rangle=\sum_{j=1}^{n}v_{i}^{j}(\textnormal{{k}})|j\rangle,\ i=1,...,r. (24)

Associated to these vectors, we have operators which create them in the many-particle Fock space:

ξik=j=1nvij(k)cjk,i=1,,r,\displaystyle\xi_{i\textnormal{{k}}}^{\dagger}=\sum_{j=1}^{n}v_{i}^{j}(\textnormal{{k}})c_{j\textnormal{{k}}}^{\dagger},\ i=1,...,r, (25)

where cjkc_{j\textnormal{{k}}}^{\dagger} denotes the original fermion creation operator, which creates a fermion with momentum k and internal degree of freedom jj, with j=1,..,nj=1,..,n. The k\textnormal{{k}}-piece of the many-particle ground state is then determined by

i=1rξik|0,\displaystyle\prod_{i=1}^{r}\xi_{i\textnormal{{k}}}^{\dagger}|0\rangle, (26)

where |0|0\rangle denotes the vacuum. This piece is only determined up to a gauge transformation. Namely, we are allowed to perform linear combinations among the states |vi,k|v_{i,\textnormal{{k}}}\rangle, i=1,,ri=1,...,r, which do not change Im[P(k)]\text{Im}\left[P(\textnormal{{k}})\right]. This amounts to taking the n×rn\times r matrix v(k)=[vij(k)]1ir,1jnv(\textnormal{{k}})=\left[v_{i}^{j}(\textnormal{{k}})\right]_{1\leq i\leq r,1\leq j\leq n} and multiplying it on the right by a k×kk\times k invertible matrix S(k)S(\textnormal{{k}}):

v(k)v(k)S(k),\displaystyle v(\textnormal{{k}})\longmapsto v(\textnormal{{k}})S(\textnormal{{k}}), (27)

and this can be done locally, in a neighbourhood of the considered momentum, in a smooth way. The effect of this gauge transformation is to transform the k piece of the many-particle ground state as

i=1rξik|0det(S(k))i=1rξik|0,\displaystyle\prod_{i=1}^{r}\xi_{i\textnormal{{k}}}^{\dagger}|0\rangle\longmapsto\det\left(S(\textnormal{{k}})\right)\prod_{i=1}^{r}\xi_{i\textnormal{{k}}}^{\dagger}|0\rangle, (28)

which clearly does not change the state. It is useful to write the projector P(k)P(\textnormal{{k}}) in terms of the matrix v(k)v(\textnormal{{k}}). To do this, note that from v(k)v(\textnormal{{k}}) we can build a unitary gauge, i.e., a gauge u(k)=[uij(k)]1ir,1jnu(\textnormal{{k}})=[u_{i}^{j}(\textnormal{{k}})]_{1\leq i\leq r,1\leq j\leq n} for which the rr columns of the new matrix form an orthonormal basis of Im[P(k)]\text{Im}\left[P(\textnormal{{k}})\right]. In matrix notation, this is equivalent to the requirement u(k)u(k)=Iru^{\dagger}(\textnormal{{k}})u(\textnormal{{k}})=I_{r}. Since this is a gauge transformation, this can be achieved by multiplying on the right by a matrix S(k)S(\textnormal{{k}}) and we have

u(k)=v(k)S(k), such that u(k)u(k)=Ir.\displaystyle u(\textnormal{{k}})=v(\textnormal{{k}})S(\textnormal{{k}}),\text{ such that }u^{\dagger}(\textnormal{{k}})u(\textnormal{{k}})=I_{r}. (29)

The unitarity condition above implies that the gauge transformation S(k)S(\textnormal{{k}}) has to satisfy

SvvS=Ir(SS)1=vv,\displaystyle S^{\dagger}v^{\dagger}vS=I_{r}\implies(SS^{\dagger})^{-1}=v^{\dagger}v, (30)

where here and below, we drop the k dependence from the expressions just to make the formulas less cumbersome to read. If we find one SS satisfying this equation, obviously, multiplying it on the right by a unitary r×rr\times r matrix will give another valid solution, which is just a manifestation of the fact that all we have done with this choice (coming from the Hilbert space inner product structure) was to reduce the full complex linear gauge group GL(r;)\text{GL}(r;\mathbb{C}) to the unitary gauge group U(r)GL(r;)\text{U}(r)\subset\text{GL}(r;\mathbb{C}). We can then write

P=uu=vSSv=v(vv)1v.\displaystyle P=uu^{\dagger}=vSS^{\dagger}v^{\dagger}=v(v^{\dagger}v)^{-1}v^{\dagger}. (31)

If we recall that the Berry connection, denoted \nabla, is determined by the projection of the exterior derivative, we can determine the one-form connection coefficients A=[Aij]1i,jrA=[A_{i}^{j}]_{1\leq i,j\leq r} in the gauge provided by vv:

v=Pdv=v((vv)1vdv)A=(vv)1vdv.\displaystyle\nabla v=Pdv=v\left((v^{\dagger}v)^{-1}v^{\dagger}dv\right)\implies A=(v^{\dagger}v)^{-1}v^{\dagger}dv. (32)

The reader should compare the above expression to the case in which the gauge is unitary, denoted uu, for which the expression of the connection coefficients, because uu=Iru^{\dagger}u=I_{r}, reduces to the familiar form

A\displaystyle A =(uu)1udu=udu\displaystyle=(u^{\dagger}u)^{-1}u^{\dagger}du=u^{\dagger}du
=[uiduj]1i,jr=[uik|d|ujk]1i,jr,\displaystyle=\left[u_{i}^{\dagger}du_{j}\right]_{1\leq i,j\leq r}=\left[\langle u_{i\textnormal{{k}}}|d|u_{j\textnormal{{k}}}\rangle\right]_{1\leq i,j\leq r}, (33)

where uiu_{i}, i=1,,ri=1,\dots,r, denote the columns of the n×rn\times r matrix uu which provide a basis of Bloch wave functions |uik|u_{i\textnormal{{k}}}\rangle, i=1,,ri=1,\dots,r, spanning the rr Bloch bands below the Fermi level at momentum k, namely Im(P(k))\text{Im}(P(\textnormal{{k}})).

Now equipped with the notion of holomorphic local coordinate, we have preferred local gauges which are holomorphic, namely, since can write

A=A(0,1)+A(1,0),\displaystyle A=A^{(0,1)}+A^{(1,0)}, (34)

where A(0,1)A^{(0,1)} only contains dz¯d\bar{z}, while A(1,0)A^{(1,0)} only contains dzdz. A holomorphic gauge is one in which A(0,1)=0A^{(0,1)}=0. By applying a gauge transformation, vvSv\mapsto vS, we have

AA=S1AS+S1dS.\displaystyle A\longmapsto A^{\prime}=S^{-1}AS+S^{-1}dS. (35)

A holomorphic gauge can be found by solving,

A(0,1)=S1A(0,1)S+S1Sz¯dz¯=0.\displaystyle A^{\prime(0,1)}=S^{-1}A^{(0,1)}S+S^{-1}\frac{\partial S}{\partial\bar{z}}d\bar{z}=0. (36)

In complex dimension one this equation has no obstruction to integrability, because there are no (2,0)(2,0)-forms, and one can always find a holomorphic gauge (see Proposition 3.7 of Ref. [39]). It is then clear the reason for calling it a holomorphic gauge: if we have any two holomorphic gauge choices, vv and vv^{\prime}, defined on the same local neighbourhood, then the gauge transformation SS relating them, v=vSv^{\prime}=vS, is seen to be holomorphic by Eq. (36) (set A(0,1)=A(0,1)=0A^{(0,1)}=A^{\prime(0,1)}=0). Essentially, by choosing local holomorphic gauges, what we are doing is to equip our occupied Bloch vector bundle EBZ2E\to\textnormal{\text{BZ}}^{2} with the structure of a holomorphic vector bundle where the Cauchy-Riemann operator is simply the (0,1)(0,1) part of the Berry connection. Since the Berry connection preserves the Hermitian inner product, i.e., parallel transport is unitary, the Berry connection becomes the Chern connection of the Hermitian holomorphic vector bundle EBZ2E\to\textnormal{\text{BZ}}^{2}(see Proposition 4.9 of [39]).

In particular, in a holomorphic gauge vv, we have

A(0,1)=(vv)1vvz¯dz¯=0vvz¯=0,\displaystyle A^{(0,1)}=\left(v^{\dagger}v\right)^{-1}v^{\dagger}\frac{\partial v}{\partial\bar{z}}d\bar{z}=0\implies v^{\dagger}\frac{\partial v}{\partial\bar{z}}=0, (37)

since vvv^{\dagger}v is invertible. Now the condition QP/z¯=0Q\partial P/\partial\bar{z}=0 means that actually, locally, we can choose a holomorphic gauge and that in such gauge v=v(z)v=v(z), i.e., the local Bloch wave functions themselves can be chosen to be holomorphic functions. To see this, we will write the condition explicitly. We begin by taking vv to be a holomorphic gauge, meaning vv/z¯=0v^{\dagger}\partial v/\partial\overline{z}=0 and we will write PP in terms of vv. Also, it will be more convenient to work with the equivalent condition Q(P/z¯)P=0Q(\partial P/\partial\bar{z})P=0 which follows due to PP being a projection operator (since Q=IPQ=I-P and QdQQ=0QdQQ=0 we have QdP=QdQ=QdQP=QdPPQdP=-QdQ=-QdQP=QdPP). We then have,

QPz¯P=Qz¯(v(vv)1v)P\displaystyle Q\frac{\partial P}{\partial\bar{z}}P=Q\frac{\partial}{\partial\bar{z}}\left(v(v^{\dagger}v)^{-1}v^{\dagger}\right)P
=Q[vz¯(vv)1vv(vv)1z¯(vv)(vv)1v\displaystyle=Q\Big{[}\frac{\partial v}{\partial\bar{z}}(v^{\dagger}v)^{-1}v^{\dagger}-v(v^{\dagger}v)^{-1}\frac{\partial}{\partial\bar{z}}\left(v^{\dagger}v\right)(v^{\dagger}v)^{-1}v^{\dagger}
+v(vv)1vz¯]P\displaystyle+v(v^{\dagger}v)^{-1}\frac{\partial v^{\dagger}}{\partial\bar{z}}\Big{]}P
=Q[vz¯(vv)1vv(vv)1vz¯v(vv)1v\displaystyle=Q\Big{[}\frac{\partial v}{\partial\bar{z}}(v^{\dagger}v)^{-1}v^{\dagger}-v(v^{\dagger}v)^{-1}\frac{\partial v^{\dagger}}{\partial\bar{z}}v(v^{\dagger}v)^{-1}v^{\dagger}
v(vv)1vvz¯v(vv)1v+v(vv)1vz¯]P\displaystyle-v(v^{\dagger}v)^{-1}v^{\dagger}\frac{\partial v}{\partial\bar{z}}v(v^{\dagger}v)^{-1}v^{\dagger}+v(v^{\dagger}v)^{-1}\frac{\partial v^{\dagger}}{\partial\bar{z}}\Big{]}P
=Q[vz¯(vv)1vv(vv)1vz¯v(vv)1v\displaystyle=Q\Big{[}\frac{\partial v}{\partial\bar{z}}(v^{\dagger}v)^{-1}v^{\dagger}-v(v^{\dagger}v)^{-1}\frac{\partial v^{\dagger}}{\partial\bar{z}}v(v^{\dagger}v)^{-1}v^{\dagger}
+v(vv)1vz¯]P\displaystyle+v(v^{\dagger}v)^{-1}\frac{\partial v^{\dagger}}{\partial\bar{z}}\Big{]}P
=Q[vz¯(vv)1vv(vv)1vz¯\displaystyle=Q\Big{[}\frac{\partial v}{\partial\bar{z}}(v^{\dagger}v)^{-1}v^{\dagger}-v(v^{\dagger}v)^{-1}\frac{\partial v^{\dagger}}{\partial\bar{z}}
+v(vv)1vz¯]P=Q[vz¯(vv)1v]P=0,\displaystyle+v(v^{\dagger}v)^{-1}\frac{\partial v^{\dagger}}{\partial\bar{z}}\Big{]}P=Q\left[\frac{\partial v}{\partial\bar{z}}(v^{\dagger}v)^{-1}v^{\dagger}\right]P=0, (38)

where we used the holomorphic gauge condition and the fact that P=v(vv)1vP=v(v^{\dagger}v)^{-1}v^{\dagger} and it is a projection operator, i.e., P2=PP^{2}=P. Since Pv=v(vv)1vv=vPv=v(v^{\dagger}v)^{-1}v^{\dagger}v=v it follows that vP=vv^{\dagger}P=v^{\dagger}, so we can write the previous condition as

Qvz¯(vv)1v=0.\displaystyle Q\frac{\partial v}{\partial\bar{z}}(v^{\dagger}v)^{-1}v^{\dagger}=0. (39)

Multiplying on the right by vv, we get

Qvz¯=0.\displaystyle Q\frac{\partial v}{\partial\bar{z}}=0. (40)

Now since from the holomorphic gauge constraint we have A(0,1)=0z¯v=Pv/z¯=0A^{(0,1)}=0\Leftrightarrow\nabla_{\frac{\partial}{\partial\bar{z}}}v=P\partial v/\partial\bar{z}=0, it follows that vv satisfies v/z¯=0\partial v/\partial\bar{z}=0, i.e. it is holomorphic. We conclude that the ground state of the band insulator can be locally described by local holomorphic Bloch wave functions |v1k,,|vrk|v_{1\textnormal{{k}}}\rangle,\dots,|v_{r\textnormal{{k}}}\rangle. This should be compared to what happens in the integer quantum Hall effect in which the ground state wave function is, in the symmetric gauge, a Gaussian factor times a holomorphic Slater determinant of single-particle states taken from the lowest Landau level. The Gaussian factor in the integer Hall effect appears because the symmetric gauge is unitary, but one could also work with a holomorphic gauge, in which case that factor does not appear and, effectively, the single-particle states from the lowest Landau level in the plane are just holomorphic functions which are square integrable with respect to an L2L^{2}-inner product. In analogy with what happens in the lowest Landau level, here the local wave functions can be chosen to effectively only depend on half of the degrees of freedom, meaning they only depend on the zz coordinate and not on z¯\bar{z}. An alternative way of understanding this feature for the integer quantum Hall effect is to use the Landau gauge, in which the lowest Landau level is identified with the Hilbert space of square integrable functions in the momentum variable with respect to which translation symmetry is preserved and, hence, it corresponds to the Hilbert space of a particle in dimension 11. The sign of the first Chern number 𝒞\mathcal{C} associated with the Fermi projector, by determining whether the map is Kähler with respect to jj or j-j (see the proof of Theorem 1 in Appendix A), controls the notion of orientation or chirality of the wave function in momentum space, for if we locally write z=|z|eiθz=|z|e^{i\theta} any holomorphic function will wind in the positive orientation of the zz plane, i.e., in the positive θ\theta direction. Alternatively and in a complementary way, we can also connect this notion of chirality to the sign of the first Chern number by considering parallel transport with respect to the Berry connection around loops. The fact that PωFSP^{*}\omega_{FS} is a symplectic form when the map is an immersion has as a consequence that PωFSP^{*}\omega_{FS} is a volume form on the Brillouin zone. This volume form (because it is nowhere vanishing) determines an orientation on the Brillouin zone and, because the map is Kähler, it is also the one induced by the complex structure jj (or j-j) which makes the map P:BZ2Grk(n)P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{k}(\mathbb{C}^{n}) holomorphic. The fact that PωFS=Ω/2P^{*}\omega_{FS}=-\Omega/2, where Ω\Omega is the Berry curvature, then tells us that in this orientation the Chern number is negative (see also the proof of Theorem 1 in Appendix A). Suppose for a moment that we are considering a single band, i.e. r=1r=1, for the sake of simplicity (equivalently, consider the top exterior power line bundle ΛrEBZ2\Lambda^{r}E\to\textnormal{\text{BZ}}^{2}, whose fibers describe Slater determinants of vectors in the corresponding fibers of the occupied Bloch vector bundle EBZ2E\to\textnormal{\text{BZ}}^{2}). If γ\gamma is a closed loop which is the boundary of some region ΣBZ2\Sigma\subset\textnormal{\text{BZ}}^{2}, with Σ\Sigma having the same orientation as BZ2\textnormal{\text{BZ}}^{2}, we will have, by Stokes’ theorem,

exp(iγA)=exp(iΣΩ)=exp(2iΣPωFS),\displaystyle\exp\left(i\int_{\gamma}A\right)=\exp\left(i\int_{\Sigma}\Omega\right)=\exp\left(-2i\int_{\Sigma}P^{*}\omega_{FS}\right), (41)

and because PωFSP^{*}\omega_{FS} is the volume form,

ΣPωFS=Σω>0\displaystyle\int_{\Sigma}P^{*}\omega_{FS}=\int_{\Sigma}\omega>0 (42)

so that the argument of the parallel transport phase, ΣΩ\int_{\Sigma}\Omega, which equals, in magnitude, twice the symplectic volume of Σ\Sigma measured with respect to ω\omega, or equivalently, due to PP being Kähler, twice the volume of Σ\Sigma measured with respect to the quantum metric gg, will necessarily have negative sign. We notice also that the statement is gauge invariant and is valid for every loop γ\gamma enclosing some arbitrary region ΣBZ2\Sigma\subset\textnormal{\text{BZ}}^{2} in the same orientation as the Brillouin zone.

The phenomenon of chirality described above is also similar to what happens in the lowest Landau level, in which the sign of the magnetic field also controls the chirality of the wave function. The fact that the geometrical degrees of freedom are cut in half is, intuitively, and, in fact, formally due to Theorem 1, also consistent with the saturation of the Cauchy-Schwarz like inequality in Theorem 1 in that this band insulator will have less geometrical degrees of freedom in its description because the quantum metric and the Berry curvature are no longer independent quantities.

We point out that a closely related result has been obtained also by Lee et al. [33]; they have shown, in the context of a single band, that the (local) Bloch wave function is a holomorphic function with respect to the isotropic flat complex structure, τ=i\tau=i, if and only if the inequality tr(g)=2det(g)|Ω12|\mathrm{tr}(g)=2\sqrt{\det(g)}\geq|\Omega_{12}| is saturated. Our result holds more generally for a general complex structure given by Eq. (21), also including a general number of bands rr. Besides, the results presented here re-express this local holomorphicity condition (expressed in terms of local Bloch wave functions) as a global holomorphicity condition for the map to the Grassmannian Grr(n)\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) induced by the Fermi projector P(k)=Θ(EFH(k))P(\textnormal{{k}})=\Theta(E_{F}-H(\textnormal{{k}})), which, provided the additional immersion condition, det(g)0\det(g)\neq 0 everywhere in BZ2\textnormal{\text{BZ}}^{2}, is translated into a Kähler map condition.

V Conclusions

In this paper we have studied the geometry of Chern insulators using the natural Kähler geometry of the space of orthogonal projectors of a given rank, the latter being associated with the number of occupied bands. In the process, we have established important results, namely Theorem 1, which renders the saturation of the Cauchy-Schwarz-like inequality equivalent to a Kähler map condition for the Fermi projector, and Theorem 2, which determines the conditions of minimal quantum volume both in the Brillouin zone and in the twist-angle space. The minimal quantum volume is, in both cases, precisely given by the topological invariant π|𝒞|\pi|\mathcal{C}|, where 𝒞\mathcal{C} is the first Chern number of the occupied bands. Supplemented with a non-degeneracy condition, the minimal volume conditions amount to these tori being endowed with the structure of a Kähler manifold inherited from the Kähler geometry of the space of quantum states. Strikingly, for the case of a single Chern band, when both quantum volumes are equal and minimal, together with the additional condition that the conformal factor relating the metrics is trivial, we get that the algebra generated by the projected density operators is isomorphic to the WW_{\infty} algebra in the fractional Hall effect [21]. Hence, this can also be used as a necessary criterion for the stability of fractional topological insulators.

We have related the Kähler map condition appearing in Theorem 1 to the possibility of writing local holomorphic Bloch functions, thus obtaining a description of the many-particle ground state of the insulator which is quite close to the one in the integer quantum Hall effect in which, apart from an overall Gaussian factor, we have a Slater determinant of holomorphic functions. The sign of the Chern number is seen to control the “chirality” of the wave function just like the sign of the magnetic field does in the case of the quantum Hall effect.

A remarkable consequence of Theorem 2 is Corollary 3 which states that if any of the quantum volumes is smaller than π\pi then we immediately have that 𝒞=0\mathcal{C}=0, i.e., we have a trivial Chern insulator. In the accompanying paper [24], it is also shown with several tight-binding models that even when 𝒞0\mathcal{C}\neq 0, the quantum volume often gives a good estimate of the Chern number. This shows how the quantum volumes, intrinsically geometric quantities, are tied to the non-trivial topology of the Chern insulator.

Additionally, we found that, for two-band systems, it is impossible to satisfy the non-degeneracy condition due to topological constraints placed by the fundamental group of the two-torus. However, provided the Berry curvature is everywhere non-negative or nonpositive, the quantum volume will automatically reach the minimal value π|𝒞|\pi|\mathcal{C}|, the reason being due to the minimal dimension of the total number of internal degrees of freedom.

Acknowledgements.
B.M. acknowledges very stimulating discussions with J. P. Nunes and J. M. Mourão. B.M. and T.O. acknowledge the anonymous referees for their insightful comments which lead to substantial improvement of this article. TO acknowledges support from JSPS KAKENHI Grant Number JP20H01845, JST PRESTO Grant Number JPMJPR19L2, JST CREST Grant Number JPMJCR19T1, and RIKEN iTHEMS. BM acknowledges the support from SQIG – Security and Quantum Information Group, the Instituto de Telecomunicações (IT) Research Unit, Ref. UIDB/50008/2020, funded by Fundação para a Ciência e a Tecnologia (FCT), European funds, namely, H2020 project SPARTA, as well as projects QuantMining POCI-01-0145-FEDER-031826 and PREDICT PTDC/CCI-CIF/29877/2017.

Appendix A Detailed mathematical proofs

Proposition 2.

If P:ΣGrr(n)P:\Sigma\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) is an immersion of a Riemann surface to the Grassmannian which is holomorphic, i.e., it satisfies, in every local complex zz coordinate on the surface,

QPz¯=0,\displaystyle Q\frac{\partial P}{\partial\bar{z}}=0,

then it is a Kähler map with respect to the usual Kähler structure on Grr(n)\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) given by (ωFS,JFS,gFS)(\omega_{FS},J_{FS},g_{FS}). The converse is also true, i.e., if P:ΣGrr(n)P:\Sigma\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) is Kähler then it is an holomorphic immersion of Σ\Sigma.

Proof.

Since PP is an immersion, we will immediately have that g=PgFSg=P^{*}g_{FS} is a Riemannian metric. Furthermore, the two-form ω=PωFS\omega=P^{*}\omega_{FS} is automatically closed. It remains to show that if PP is holomorphic then ω(,j)=g\omega(\cdot,j\cdot)=g and also that ω\omega is non-degenerate. Recall that PP being holomorphic means that the differential of PP intertwines the two complex structures (see, for example, Ref. [30]), i.e.,

JFSdP=dPj.\displaystyle J_{FS}\circ dP=dP\circ j. (43)

If we take a complex coordinate zz associated to jj over Σ\Sigma, so that

j(z)=iz and j(z¯)=iz¯,\displaystyle j\left(\frac{\partial}{\partial z}\right)=i\frac{\partial}{\partial z}\text{ and }j\left(\frac{\partial}{\partial\overline{z}}\right)=-i\frac{\partial}{\partial\overline{z}}, (44)

we have

JFS(Pz)=iPz.\displaystyle J_{FS}\left(\frac{\partial P}{\partial z}\right)=i\frac{\partial P}{\partial z}. (45)

Observing that dP=QdPP+PdPQdP=QdPP+PdPQ, that (QdPP)JFS=iQdPP\left(QdPP\right)\circ J_{FS}=iQdPP and that (PdPQ)JFS=iPdPQ\left(PdPQ\right)\circ J_{FS}=-iPdPQ, we get,

iQPzPiPPzQ=iQPzP+iPPzQ.\displaystyle iQ\frac{\partial P}{\partial z}P-iP\frac{\partial P}{\partial z}Q=iQ\frac{\partial P}{\partial z}P+iP\frac{\partial P}{\partial z}Q. (46)

This then implies the equation

PPzQ=0,\displaystyle P\frac{\partial P}{\partial z}Q=0, (47)

or, equivalently, the Hermitian conjugate

QPz¯P=QPz¯=0,\displaystyle Q\frac{\partial P}{\partial\overline{z}}P=Q\frac{\partial P}{\partial\overline{z}}=0, (48)

where the last equality follows from QdPQ=Q(dQ)Q=0QdPQ=Q(-dQ)Q=0. Observe that

PωFS=i2Tr(PdPdP)\displaystyle P^{*}\omega_{FS}=-\frac{i}{2}\textnormal{\text{Tr}}\left(PdP\wedge dP\right)
=i2(Tr(PPz¯Pz)Tr(PPzPz¯))dz¯dz\displaystyle=-\frac{i}{2}\left(\textnormal{\text{Tr}}\left(P\frac{\partial P}{\partial\bar{z}}\frac{\partial P}{\partial z}\right)-\textnormal{\text{Tr}}\left(P\frac{\partial P}{\partial z}\frac{\partial P}{\partial\bar{z}}\right)\right)d\bar{z}\wedge dz
=i2(Tr(PPz¯QPz)Tr(PPzQPz¯))dz¯dz\displaystyle=-\frac{i}{2}\left(\textnormal{\text{Tr}}\left(P\frac{\partial P}{\partial\bar{z}}Q\frac{\partial P}{\partial z}\right)-\textnormal{\text{Tr}}\left(P\frac{\partial P}{\partial z}Q\frac{\partial P}{\partial\bar{z}}\right)\right)d\bar{z}\wedge dz
=i2Tr(PPz¯QPz)dz¯dz=i2Tr(PPz¯Pz)dzdz¯,\displaystyle=-\frac{i}{2}\textnormal{\text{Tr}}\left(P\frac{\partial P}{\partial\bar{z}}Q\frac{\partial P}{\partial z}\right)d\bar{z}\wedge dz=\frac{i}{2}\textnormal{\text{Tr}}\left(P\frac{\partial P}{\partial\bar{z}}\frac{\partial P}{\partial z}\right)dz\wedge d\bar{z},

where we used PdPP=0PdPP=0 and QdQQ=0QdQQ=0 which imply that PdP=PdPQ=dPQ=dQQPdP=PdPQ=dPQ=-dQQ and similarly dPP=QdPP=QdPdPP=QdPP=QdP. Also, by the same reasoning

PgFS\displaystyle P^{*}g_{FS} =Tr(PdPdP)=Tr(PPz¯Pz)dz¯dz.\displaystyle=\textnormal{\text{Tr}}\left(PdPdP\right)=\textnormal{\text{Tr}}\left(P\frac{\partial P}{\partial\bar{z}}\frac{\partial P}{\partial z}\right)d\bar{z}dz.

Using Eq. (44) we get

i2dzdz¯(,j)\displaystyle\frac{i}{2}dz\wedge d\bar{z}(\cdot,j\cdot) =i2(dzdz¯dz¯dz)(,j)\displaystyle=\frac{i}{2}\left(dz\otimes d\bar{z}-d\bar{z}\otimes dz\right)(\cdot,j\cdot)
=12(dzdz¯+dz¯dz)=dzdz¯.\displaystyle=\frac{1}{2}\left(dz\otimes d\bar{z}+d\bar{z}\otimes dz\right)=dzd\bar{z}.

Thus, we conclude

ω(,j)=PωFS(,j)=PgFS=g.\displaystyle\omega(\cdot,j\cdot)=P^{*}\omega_{FS}(\cdot,j\cdot)=P^{*}g_{FS}=g.

Finally, ω\omega is non-degenerate because gg is so and they are both determined by the same function, namely Tr(PPz¯Pz)\textnormal{\text{Tr}}\left(P\frac{\partial P}{\partial\overline{z}}\frac{\partial P}{\partial z}\right).

The converse is also true, for if P:ΣGrr(n)P:\Sigma\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) is Kähler it automatically is an immersion since that is implied by gg being non-degenerate and it has to be holomorphic in order to be Kähler. Hence, the proof is concluded. ∎

Proposition 3 (Cauchy-Schwarz-type inequality and relation to the Wirtinger inequality).

We have the inequality

det[gij(k)]|Ω12(k)|2,\displaystyle\sqrt{\det[g_{ij}(\textnormal{{k}})]}\geq\frac{|\Omega_{12}(\textnormal{{k}})|}{2},

with the equality, for a given kBZ2\textnormal{{k}}\in\textnormal{\text{BZ}}^{2}, given by at least one of the following conditions

  • (i)

    Q(k)Pk1(k)=0Q(\textnormal{{k}})\frac{\partial P}{\partial k_{1}}(\textnormal{{k}})=0,

  • (ii)

    Q(k)Pk2(k)=0Q(\textnormal{{k}})\frac{\partial P}{\partial k_{2}}(\textnormal{{k}})=0,

  • (iii)

    λ{0}:Q(k)Pk1(k)=λQ(k)Pk2(k)\exists\lambda\in\mathbb{C}-\{0\}:Q(\textnormal{{k}})\frac{\partial P}{\partial k_{1}}(\textnormal{{k}})=\lambda Q(\textnormal{{k}})\frac{\partial P}{\partial k_{2}}(\textnormal{{k}}).

The above result is a manifestation of the so-called Wirtinger inequality for general Kähler manifolds [40]. The Wirtinger inequality states that on a Kähler manifold, the kkth exterior power of the Kähler form, when evaluated on a simple (decomposable) 2kk-vector of unit volume, is bounded above by k!k!. The case at hand follows with k=1k=1.

Proof.

The inequality has been first noted by Roy [21], for the case r=1r=1. A physicist-friendly proof, for r1r\geq 1, can be found in the accompanying paper [24]. Below we give an alternative proof using the Wirtinger inequality. Before doing that, we give a few additional remarks here to the proof of Ref. [24] because the conditions for the saturation of the inequality are stated here in terms of projectors. We first observe that conditions (i), (ii), (iii) are statements about the image of the differential of PP at k because dkP(/ki)d_{\textnormal{{k}}}P(\partial/\partial k_{i}) is given by Q(k)(P/ki)(k)+H.c.Q(\textnormal{{k}})\left(\partial P/\partial k_{i}\right)(\textnormal{{k}})+\text{H.c.}, i=1,2i=1,2. Moreover, the fact that the Plücker embedding is an isometry allows us to reduce the proof of the proposition to the case of r=1r=1, i.e., for the complex projective space (Λrn)Pm=Gr1(m+1)\mathbb{P}(\Lambda^{r}\mathbb{C}^{n})\cong\mathbb{C}P^{m}=\textnormal{\text{Gr}}_{1}(\mathbb{C}^{m+1}) with m=(nr)1m=\binom{n}{r}-1. Effectively, this amounts to, for each kBZ2\textnormal{{k}}\in\textnormal{\text{BZ}}^{2}, replacing the initial rank rr projector by the rank 11 projector associated to the wedge product (Slater determinant) associated to an orthonormal basis of the image of the former. In this case, we have a rank 11 projector locally specified by a (normalized) Bloch wave function P(k)=|ψ(k)ψ(k)|P(\textnormal{{k}})=|\psi(\textnormal{{k}})\rangle\langle\psi(\textnormal{{k}})|. It follows from the proof in Ref. [24] that the inequality is saturated if and only if Q(k)/k1|ψ(k)=0Q(\textnormal{{k}})\partial/\partial k_{1}|\psi(\textnormal{{k}})\rangle=0 or Q(k)/k2|ψ(k)=0Q(\textnormal{{k}})\partial/\partial k_{2}|\psi(\textnormal{{k}})\rangle=0 or Q(k)/k1|ψ(k)=λQ(k)/k2|ψ(k)=0Q(\textnormal{{k}})\partial/\partial k_{1}|\psi(\textnormal{{k}})\rangle=\lambda Q(\textnormal{{k}})\partial/\partial k_{2}|\psi(\textnormal{{k}})\rangle=0 for some λ\lambda\in\mathbb{C}. Additionally, observing that

Q(k)Pki(k)=Q(k)(ki|ψ(k))ψ(k)|, for i=1,2,\displaystyle Q(\textnormal{{k}})\frac{\partial P}{\partial k_{i}}(\textnormal{{k}})=Q(\textnormal{{k}})\left(\frac{\partial}{\partial k_{i}}|\psi(\textnormal{{k}})\rangle\right)\langle\psi(\textnormal{{k}})|,\text{ for }i=1,2, (49)

because Q(k)|ψ(k)=0Q(\textnormal{{k}})|\psi(\textnormal{{k}})\rangle=0, the proof is complete.

Now we consider the relation with the Wirtinger inequality, which provides an alternative proof. For k=1k=1, in our particular setting, this means that if e1,e2e_{1},e_{2} are orthonormal with respect to gFSg_{FS}, we have

|ωFS(e1,e2)|1\displaystyle|\omega_{FS}(e_{1},e_{2})|\leq 1 (50)

If e1e_{1} and e2e_{2} form an orthonormal basis for the image of the differential dPdP for some momentum k, then,

ei=j=12aijdP(kj),i=1,2\displaystyle e_{i}=\sum_{j=1}^{2}a_{i}^{j}dP\left(\frac{\partial}{\partial k_{j}}\right),i=1,2 (51)

for some invertible matrix a=[aij]1i,j2a=[a_{i}^{j}]_{1\leq i,j\leq 2}. Then,

ωFS(e1,e2)\displaystyle\omega_{FS}(e_{1},e_{2}) =det(aij)ωFS(dP(k1),dP(k1))\displaystyle=\det(a_{i}^{j})\omega_{FS}(dP\left(\frac{\partial}{\partial k_{1}}\right),dP\left(\frac{\partial}{\partial k_{1}}\right))
=det(aij)(PωFS)12=12Ω12det(aij).\displaystyle=\det(a_{i}^{j})\left(P^{*}\omega_{FS}\right)_{12}=-\frac{1}{2}\Omega_{12}\det(a_{i}^{j}).

Next, observe that since gFS(ei,ej)=δijg_{FS}(e_{i},e_{j})=\delta_{ij}, 1i,j21\leq i,j\leq 2, we have that

gij=(PgFS)(ki,kj)=k,l=12(a1)ik(a1)jlδkl.\displaystyle g_{ij}=(P^{*}g_{FS})\left(\frac{\partial}{\partial k_{i}},\frac{\partial}{\partial k_{j}}\right)=\sum_{k,l=1}^{2}(a^{-1})^{k}_{i}(a^{-1})^{l}_{j}\delta_{kl}.

As a consequence, det(gij)=det(a1)2\det(g_{ij})=\det(a^{-1})^{2}. Thus, the Wirtinger inequality implies

|ωFS(e1,e2)|=|12Ω12det(aij)|1\displaystyle|\omega_{FS}(e_{1},e_{2})|=\left|-\frac{1}{2}\Omega_{12}{\det}(a_{i}^{j})\right|\leq 1 (52)

or

det(g)12|Ω12|,\displaystyle\sqrt{\det(g)}\geq\frac{1}{2}|\Omega_{12}|, (53)

as claimed. Concerning the saturation of the inequality, two cases are to be considered in light of Wirtinger’s inequality. If at a given point kBZ2\textnormal{{k}}\in\textnormal{\text{BZ}}^{2} the mapping is not an immersion the dimension of image of the differential of PP at k is smaller than 22. This implies that det[gij(k)]=Ω12(k)=0\sqrt{\det[g_{ij}(\textnormal{{k}})]}=\Omega_{12}(\textnormal{{k}})=0 at that point because we can not build an orthonormal basis with two elements in the image. This covers the cases (i) and (ii), and also (iii) provided λ\lambda is real since they imply that the differential of PP at k is not injective (the vanishing of i=12aiQ(k)P/ki(k)\sum_{i=1}^{2}a^{i}Q(\textnormal{{k}})\partial P/\partial k_{i}(\textnormal{{k}}) is equivalent to the vanishing of the image under the differential of PP at k of the tangent vector i=12ai/kiTkBZ2\sum_{i=1}^{2}a^{i}\partial/\partial k_{i}\in T_{\textnormal{{k}}}\textnormal{\text{BZ}}^{2} with aia^{i}\in\mathbb{R}, i=1,2i=1,2). We are left with the non-degenerate case, i.e., when the image of the differential of PP at k is two-dimensional, in which we can apply Wirtinger’s inequality and consider the situation in which it is saturated. We will show in the following that this is equivalent to the case (iii) with Im(λ)0\text{Im}\left(\lambda\right)\neq 0. Because (ωFS,JFS,gFS)(\omega_{FS},J_{FS},g_{FS}) is a compatible triple, the saturation of Wirtinger’s inequality holds iff e2=±JFS(e1)e_{2}=\pm J_{FS}(e_{1}). Indeed, if e2=±JFS(e1)e_{2}=\pm J_{FS}(e_{1}) we have |ωFS(e1,e2)|=gFS(e1,e1)=1|\omega_{FS}(e_{1},e_{2})|=g_{FS}(e_{1},e_{1})=1. Conversely if |ωFS(e1,e2)|=1|\omega_{FS}(e_{1},e_{2})|=1, we have

|ωFS(e1,e2)|\displaystyle|\omega_{FS}(e_{1},e_{2})| =|g(JFS(e1),e2)|\displaystyle=|g(J_{FS}(e_{1}),e_{2})|
=1\displaystyle=1
=gFS(e1,e1)=gFS(e2,e2).\displaystyle=g_{FS}(e_{1},e_{1})=g_{FS}(e_{2},e_{2}). (54)

Furthermore, by skew-symmetry, ωFS(ei,ei)=0=gFS(ei,J(ei))\omega_{FS}(e_{i},e_{i})=0=g_{FS}(e_{i},J(e_{i})), i=1,2i=1,2 and, by compatibility, gFS(JFS(e1),JFS(e1))=ωFS(e1,JFS(e1))=g(e1,e1)g_{FS}(J_{FS}(e_{1}),J_{FS}(e_{1}))=\omega_{FS}(e_{1},J_{FS}(e_{1}))=g(e_{1},e_{1}). It follows that JFS(e1)J_{FS}(e_{1}) is orthogonal to e1e_{1} and it has norm 11. From the saturation of the Cauchy-Schwarz inequality for the metric gFSg_{FS}

1=gFS(e2,e2)gFS(JFS(e1),JFS(e1))\displaystyle 1\!=\!g_{FS}(e_{2},e_{2})g_{FS}(J_{FS}(e_{1}),J_{FS}(e_{1})) |gFS(e2,JFS(e1))|2\displaystyle\geq|g_{FS}(e_{2},J_{FS}(e_{1}))|^{2}
=|ωFS(e1,e2)|2\displaystyle=|\omega_{FS}(e_{1},e_{2})|^{2}
=1,\displaystyle=1, (55)

it follows that J(e1)=±e2J(e_{1})=\pm e_{2}. Now if e1,e2e_{1},e_{2} form an orthonormal basis of the image of dPdP at k as in Eq. (51) this in turn implies that

j=12a1j(iQPkjiPkjQ)=±j=12a2j(QPkj+PkjQ).\displaystyle\sum_{j=1}^{2}a_{1}^{j}\left(iQ\frac{\partial P}{\partial k_{j}}-i\frac{\partial P}{\partial k_{j}}Q\right)=\pm\sum_{j=1}^{2}a_{2}^{j}\left(Q\frac{\partial P}{\partial k_{j}}+\frac{\partial P}{\partial k_{j}}Q\right). (56)

Due to PP and QQ being orthogonal to each other it follows that

j=12a1jiQPkj=±j=12a2jQPkj,\displaystyle\sum_{j=1}^{2}a_{1}^{j}iQ\frac{\partial P}{\partial k_{j}}=\pm\sum_{j=1}^{2}a_{2}^{j}Q\frac{\partial P}{\partial k_{j}}, (57)

which further implies

(ia11a21)QPk1=(ia21±a22)QPk2.\displaystyle\left(ia_{1}^{1}\mp a_{2}^{1}\right)Q\frac{\partial P}{\partial k_{1}}=\left(-ia_{2}^{1}\pm a_{2}^{2}\right)Q\frac{\partial P}{\partial k_{2}}. (58)

Now computing

λ=(ia21±a22)(ia11a21)\displaystyle\lambda=\frac{(-ia^{1}_{2}\pm a^{2}_{2})}{(ia^{1}_{1}\mp a^{1}_{2})} =(ia21±a22)(ia11a21)((a21)2+(a11)2)\displaystyle=\frac{(-ia^{1}_{2}\pm a^{2}_{2})(-ia^{1}_{1}\mp a^{1}_{2})}{\left((a^{1}_{2})^{2}+(a^{1}_{1})^{2}\right)}
=a11a21a21a22+i(a11a22±a21a21)((a21)2+(a11)2)\displaystyle=\frac{-a^{1}_{1}a^{1}_{2}-a^{1}_{2}a^{2}_{2}+i\left(\mp a^{1}_{1}a^{2}_{2}\pm a^{1}_{2}a^{1}_{2}\right)}{\left((a^{1}_{2})^{2}+(a^{1}_{1})^{2}\right)}
=a11a21a21a22±idet(aji)((a21)2+(a11)2),\displaystyle=\frac{-a^{1}_{1}a^{1}_{2}-a^{1}_{2}a^{2}_{2}\pm i\det(a^{i}_{j})}{\left((a^{1}_{2})^{2}+(a^{1}_{1})^{2}\right)}, (59)

which, because det(aji)0\det(a^{i}_{j})\neq 0 and (a21)2+(a11)20(a^{1}_{2})^{2}+(a^{1}_{1})^{2}\neq 0 (since e20e_{2}\neq 0), implies condition (iii) with Im(λ)0\mbox{Im}(\lambda)\neq 0. This concludes our proof. ∎

Theorem 1.

Suppose P:BZ2Grr(n)P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) is an immersion, then the Cauchy-Schwarz like inequality is saturated

det(gij(k))=Ω12(k)2, for all kBZ2,\displaystyle\sqrt{\det(g_{ij}(\textnormal{{k}}))}=\frac{\mp\Omega_{12}(\textnormal{{k}})}{2},\text{ for all }\textnormal{{k}}\in\textnormal{\text{BZ}}^{2},

if and only if the map PP is Kähler with respect to the triple of structures (ω,±j,g)(\omega,\pm j,g).

Proof.

By assumption, PP is an immersion so the equality

det(gij(k))=|Ω12(k)|2,\displaystyle\sqrt{\det(g_{ij}(\textnormal{{k}}))}=\frac{|\Omega_{12}(\textnormal{{k}})|}{2},

implies that condition (iii) of Proposition 3 holds pointwise, i.e., there exists a function λC(BZ2)\lambda\in C^{\infty}(\textnormal{\text{BZ}}^{2}) such that

Q(k)Pk1(k)=λ(k)Q(k)Pk2(k), for all kBZ2.\displaystyle Q(\textnormal{{k}})\frac{\partial P}{\partial k_{1}}(\textnormal{{k}})=\lambda(\textnormal{{k}})Q(\textnormal{{k}})\frac{\partial P}{\partial k_{2}}(\textnormal{{k}}),\text{ for all }\textnormal{{k}}\in\textnormal{\text{BZ}}^{2}.

Observe furthermore that by the fact that PP is an immersion we have that gg is non-degenerate and hence Im[λ(k)]0\mbox{Im}[\lambda(\textnormal{{k}})]\neq 0 for all k\textnormal{{k}}\in\mathbb{R}. Otherwise the quantities QP/kiQ\partial P/\partial k_{i}, i=1,2i=1,2, which represent the images of the tangent vectors /ki\partial/\partial k_{i}, i=,1,,2i=,1,,2, would be linearly dependent.

Observe that if we introduce a complex variable zz such that

z¯=12(k1λk2),\displaystyle\frac{\partial}{\partial\overline{z}}=\frac{1}{2}\left(\frac{\partial}{\partial k_{1}}-\lambda\frac{\partial}{\partial k_{2}}\right),

the statement is

QPz¯=0.\displaystyle Q\frac{\partial P}{\partial\overline{z}}=0.

To find zz we would have to solve,

k1z¯=12 and k2z¯=λ2,\displaystyle\frac{\partial k_{1}}{\partial\bar{z}}=\frac{1}{2}\text{ and }\frac{\partial k_{2}}{\partial\bar{z}}=-\frac{\lambda}{2},
k1z=12 and k2z=λ¯2.\displaystyle\frac{\partial k_{1}}{\partial z}=\frac{1}{2}\text{ and }\frac{\partial k_{2}}{\partial z}=-\frac{\overline{\lambda}}{2}.

Equivalently, since

(12[11λ¯λ])1=1iIm(λ)[λ1λ¯1],\displaystyle\left(\frac{1}{2}\left[\begin{array}[]{cc}1&1\\ -\bar{\lambda}&-\lambda\end{array}\right]\right)^{-1}=\frac{1}{i\text{Im}(\lambda)}\left[\begin{array}[]{cc}\lambda&1\\ -\bar{\lambda}&-1\end{array}\right],

well defined since Im(λ)0\mbox{Im}(\lambda)\neq 0, we have

dz=1iIm(λ)(λdk1+dk2) and dz¯=1iIm(λ)(λ¯dk1+dk2).\displaystyle dz\!=\!\frac{1}{i\text{Im}(\lambda)}\!\left(\lambda dk_{1}+dk_{2}\right)\!\text{ and }d\overline{z}\!=\!-\frac{1}{i\text{Im}(\lambda)}\!\left(\overline{\lambda}dk_{1}+dk_{2}\right)\!.

The resulting partial differential equations are locally integrable (see, for instance, Example 5 of Chapter 1 in Ref. [31] or Theorem 4.16 of Ref. [41]) because by solving them we are finding local isothermal coordinates, in two dimensions, with respect to PgFSP^{*}g_{FS}. These local systems of coordinates then glue together holomorphically giving the Brillouin zone the structure of a complex manifold. We now refer to Proposition 2 and one side of the implication is proved. Note, however, that the new local complex coordinates zz can induce the opposite orientation with respect to the usual one induced by the (multi-valued) complex coordinate k1+ik2k_{1}+ik_{2}. This is the case when the 1st Chern number is positive. In this case, the map is Kähler with respect to the complex structure one gets by rotation of 9090 degrees as determined by gg and the orientation opposite to the standard orientation of the Brillouin zone; namely the almost complex structure is j-j, with respect to jj as defined in Sec. II. This follows from the fact that the function Tr(PPz¯Pz)\textnormal{\text{Tr}}\left(P\frac{\partial P}{\partial\overline{z}}\frac{\partial P}{\partial z}\right) determining PωFS=PΩ0/2=Ω/2P^{*}\omega_{FS}=-P^{*}\Omega_{0}/2=-\Omega/2 in the local holomorphic coordinate zz is always positive (because it also determines the metric) and so the integral in an orientation consistent with it will be positive.

Let us now consider the converse, i.e., if P:BZ2Grr(n)P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{r}(\mathbb{C}^{n}) is Kähler (which automatically implies that PP is an immersion) with respect to the triple of structures (ω,±j,g)(\omega,\pm j,g), then the inequality

det(gij(k))=|Ω12(k)|2, for all kBZ2,\displaystyle\sqrt{\det(g_{ij}(\textnormal{{k}}))}=\frac{|\Omega_{12}(\textnormal{{k}})|}{2},\text{ for all }\textnormal{{k}}\in\textnormal{\text{BZ}}^{2},

is saturated. The reason is that the Kähler condition ω(,±j)=g(,)\omega(\cdot,\pm j\cdot)=g(\cdot,\cdot), with ω=PωFS\omega=P^{*}\omega_{FS} and g=PgFSg=P^{*}g_{FS} is locally given by

j=1212Ωij(k)jkj(k)=gik(k),\displaystyle\mp\sum_{j=1}^{2}\frac{1}{2}\Omega_{ij}(\textnormal{{k}})j^{j}_{k}(\textnormal{{k}})=g_{ik}(\textnormal{{k}}),

with jji(k)j^{i}_{j}(\textnormal{{k}}) being the matrix elements of the complex structure jj, i.e.,

j(ki)=j=12jijkj,i=1,2.\displaystyle j\left(\frac{\partial}{\partial k_{i}}\right)=\sum_{j=1}^{2}j^{j}_{i}\frac{\partial}{\partial k_{j}},\ i=1,2.

In matrix form

±ωj=12Ωj=g.\displaystyle\pm\omega j=\mp\frac{1}{2}\Omega j=g.

Taking determinants, we get

14Ω122det(j)=14|Ω12|2det(j)=det(g).\displaystyle\frac{1}{4}\Omega_{12}^{2}\det(j)=\frac{1}{4}|\Omega_{12}|^{2}\det(j)=\det(g).

Since the right-hand side has positive determinant it follows that det(j)>1\det(j)>1 and hence, because j2=Ij^{2}=-I we have det(j)2=1\det(j)^{2}=1 so det(j)=1\det(j)=1. Hence we get the equality

12|Ω12|=det(g),\displaystyle\frac{1}{2}|\Omega_{12}|=\sqrt{\det(g)},

as desired. The theorem is, thus, proved. ∎

Theorem 2.

The following inequalities hold

π|𝒞|volg(BZ2)volg~(Tθ2),\displaystyle\pi|\mathcal{C}|\leq\mbox{vol}_{g}(\textnormal{\text{BZ}}^{2})\leq\mbox{vol}_{\widetilde{g}}(T^{2}_{\theta}),

where the left-hand side inequality is saturated iff, for every kBZ2\textnormal{{k}}\in\textnormal{\text{BZ}}^{2}, det(g)=|Ω12|2\sqrt{\det(g)}=\frac{|\Omega_{12}|}{2} and provided Ω12\Omega_{12} does not change sign over BZ2\textnormal{\text{BZ}}^{2}, and the right-hand side inequality is satisfied iff

gij(k)=e2f(k)g~ij(θ=k), 1i,j2,\displaystyle g_{ij}(\textnormal{{k}})=e^{2f(\textnormal{{k}})}\widetilde{g}_{ij}(\theta=\textnormal{{k}}),\ 1\leq i,j\leq 2,

for some function fC(BZ2)f\in C^{\infty}(\textnormal{\text{BZ}}^{2}), i.e., gg is related to g~\widetilde{g} by a Weyl rescaling, implying that, provided PP is an immersion, they share the same complex structure jj (and hence the same modular parameter).

Proof.

The left-hand side inequality is a consequence of Proposition 3 and the definitions of the 1st Chern number and the quantum volume. The right-hand side inequality follows from an L2L^{2}-Cauchy-Schwarz inequality; a detailed derivation can be found in the accompanying paper [24]. Here we will just prove the last additional statement. The saturation of the inequality volg(BZ2)volg~(Tθ2)\textnormal{\text{vol}}_{g}(\textnormal{\text{BZ}}^{2})\leq\textnormal{\text{vol}}_{\widetilde{g}}(T^{2}_{\theta}) is shown in [24] to be equivalent to the equality

gij(k)=g11(k)cij,for every kBZ2,\displaystyle g_{ij}(\textnormal{{k}})=g_{11}(\textnormal{{k}})c_{ij},\text{for every }\textnormal{{k}}\in\textnormal{\text{BZ}}^{2},

for some constants cijc_{ij}, 1i,j21\leq i,j\leq 2. This means that gg and g~\widetilde{g}, defined by Eq. (14), differ by a Weyl rescaling – hence, provided PP is an immersion so that gg is non-degenerate, they share the same complex structure. It then follows that

g~ij(θ)=BZ2d2k(2π)2gij(k)=(BZ2d2k(2π)2g11(k))cij.\displaystyle\widetilde{g}_{ij}(\theta)=\int_{\textnormal{\text{BZ}}^{2}}\frac{d^{2}k}{(2\pi)^{2}}g_{ij}(\textnormal{{k}})=\left(\int_{\textnormal{\text{BZ}}^{2}}\frac{d^{2}k}{(2\pi)^{2}}g_{11}(\textnormal{{k}})\right)c_{ij}.

The result follows immediately, for

e2f(k)=g11(k)BZ2d2k(2π)2g11(k).\displaystyle e^{2f(\textnormal{{k}})}=\frac{g_{11}(\textnormal{{k}})}{\int_{\textnormal{\text{BZ}}^{2}}\frac{d^{2}k}{(2\pi)^{2}}g_{11}(\textnormal{{k}})}.

Theorem 3.

For two-dimensional two-band models, there must exist a point in the Brillouin zone where det(g)=0\det(g)=0. In other words, the map P:BZ2Gr1(2)=P1S2P:\textnormal{\text{BZ}}^{2}\to\textnormal{\text{Gr}}_{1}(\mathbb{C}^{2})=\mathbb{C}P^{1}\cong S^{2} cannot be an immersion.

Proof.

In this particular case, since Gr1(2)=P1S2\textnormal{\text{Gr}}_{1}(\mathbb{C}^{2})=\mathbb{C}P^{1}\cong S^{2} is the Bloch sphere, we get a map from a torus to the two-sphere, BZ2kn(k)S2\textnormal{\text{BZ}}^{2}\ni\textnormal{{k}}\mapsto n(\textnormal{{k}})\in S^{2}. If this map were an immersion, then by the inverse function theorem, it would be a local diffeomorphism. As a consequence it would be a covering map. Now Proposition 1.3.1. of Ref. [42], tells us that if we have a covering map then the induced group homomorphism in fundamental groups is injective. The fundamental group of the Brillouin torus BZ2\textnormal{\text{BZ}}^{2} is the \mathbb{Z}\oplus\mathbb{Z} Abelian group generated by the non-contractible loops associated with the k1k_{1} and k2k_{2} directions. Meanwhile, that of the sphere is the trivial group since any loop is contractible to a point. As a consequence, there can be no injective group homomorphism between these two groups. Hence, we find that there can be no immersion map from the torus to the sphere and, as a conclusion, the quantum metric must be degenerate, i.e. non-invertible, somewhere in the Brillouin zone. ∎

References

  • Resta [2011] R. Resta, The insulating state of matter: a geometrical theory, The European Physical Journal B 79, 121 (2011).
  • Resta [2018] R. Resta, Theory of the insulating state, La Rivista del Nuovo Cimento 41, 463 (2018).
  • Matsuura and Ryu [2010] S. Matsuura and S. Ryu, Momentum space metric, nonlocal operator, and topological insulators, Physical Review B 82, 245113 (2010).
  • Ozawa and Goldman [2019] T. Ozawa and N. Goldman, Probing localization and quantum geometry by spectroscopy, Phys. Rev. Research 1, 032019 (2019).
  • Marzari and Vanderbilt [1997] N. Marzari and D. Vanderbilt, Maximally localized generalized Wannier functions for composite energy bands, Physical review B 56, 12847 (1997).
  • Resta [1998] R. Resta, Quantum-mechanical position operator in extended systems, Physical Review Letters 80, 1800 (1998).
  • Resta and Sorella [1999] R. Resta and S. Sorella, Electron localization in the insulating state, Physical Review Letters 82, 370 (1999).
  • Souza et al. [2000] I. Souza, T. Wilkens, and R. M. Martin, Polarization and localization in insulators: Generating function approach, Physical Review B 62, 1666 (2000).
  • Kolodrubetz et al. [2013] M. Kolodrubetz, V. Gritsev, and A. Polkovnikov, Classifying and measuring geometry of a quantum ground state manifold, Phys. Rev. B 88, 064304 (2013).
  • Neupert et al. [2013] T. Neupert, C. Chamon, and C. Mudry, Measuring the quantum geometry of Bloch bands with current noise, Phys. Rev. B 87, 245103 (2013).
  • Ozawa [2018] T. Ozawa, Steady-state Hall response and quantum geometry of driven-dissipative lattices, Phys. Rev. B 97, 041108(R) (2018).
  • Ozawa and Goldman [2018] T. Ozawa and N. Goldman, Extracting the quantum metric tensor through periodic driving, Phys. Rev. B 97, 201117(R) (2018).
  • Bleu et al. [2018a] O. Bleu, D. D. Solnyshkov, and G. Malpuech, Measuring the quantum geometric tensor in two-dimensional photonic and exciton-polariton systems, Phys. Rev. B 97, 195422 (2018a).
  • Bleu et al. [2018b] O. Bleu, G. Malpuech, Y. Gao, and D. D. Solnyshkov, Effective theory of nonadiabatic quantum evolution based on the quantum geometric tensor, Phys. Rev. Lett. 121, 020401 (2018b).
  • Klees et al. [2020] R. L. Klees, G. Rastelli, J. C. Cuevas, and W. Belzig, Microwave spectroscopy reveals the quantum geometric tensor of topological Josephson matter, Phys. Rev. Lett. 124, 197002 (2020).
  • Asteria et al. [2019] L. Asteria, D. T. Tran, T. Ozawa, M. Tarnowski, B. S. Rem, N. Fläschner, K. Sengstock, N. Goldman, and C. Weitenberg, Measuring quantized circular dichroism in ultracold topological matter, Nature physics 15, 449 (2019).
  • Gianfrate et al. [2020] A. Gianfrate, O. Bleu, L. Dominici, V. Ardizzone, M. De Giorgi, D. Ballarini, G. Lerario, K. West, L. Pfeiffer, D. Solnyshkov, et al., Measurement of the quantum geometric tensor and of the anomalous Hall drift, Nature 578, 381 (2020).
  • Witten [2016] E. Witten, Fermion path integrals and topological phases, Reviews of Modern Physics 88, 035001 (2016).
  • Kudo et al. [2019] K. Kudo, H. Watanabe, T. Kariyado, and Y. Hatsugai, Many-body Chern number without integration, Phys. Rev. Lett. 122, 146601 (2019).
  • Mera [2020a] B. Mera, Localization anisotropy and complex geometry in two-dimensional insulators, Phys. Rev. B 101, 115128 (2020a).
  • Roy [2014] R. Roy, Band geometry of fractional topological insulators, Phys. Rev. B 90, 165139 (2014).
  • Jackson et al. [2015] T. S. Jackson, G. Möller, and R. Roy, Geometric stability of topological lattice phases, Nature Communications 610.1038/ncomms9629 (2015).
  • Girvin et al. [1986] S. M. Girvin, A. H. MacDonald, and P. M. Platzman, Magneto-roton theory of collective excitations in the fractional quantum Hall effect, Phys. Rev. B 33, 2481 (1986).
  • Ozawa and Mera [2021] T. Ozawa and B. Mera, Relations between topology and the quantum metric for Chern insulators, Phys. Rev. B 104, 045103 (2021).
  • Kitaev [2009] A. Kitaev, Periodic table for topological insulators and superconductors, in AIP Conference Proceedings, Vol. 1134 (AIP, 2009) pp. 22–30.
  • Park [2008] E. Park, Complex topological K-theory, Vol. 111 (Cambridge University Press, 2008).
  • Mera [2020b] B. Mera, The product of two independent Su-Schrieffer-Heeger chains yields a two-dimensional Chern insulator, Phys. Rev. B 102, 155150 (2020b).
  • da Silva [2001] A. C. da Silva, Lectures on symplectic geometry, Vol. 3575 (Springer, 2001).
  • Nakahara [2003] M. Nakahara, Geometry, Topology and Physics, 2nd Edition (CRC Press, 2003).
  • Huybrechts [2005] D. Huybrechts, Complex geometry: an introduction (Springer Science & Business Media, 2005).
  • Chern [1967] S.-S. Chern, Complex manifolds without potential theory, 15 (Springer, 1967).
  • Nguyen et al. [2017] D. X. Nguyen, T. Can, and A. Gromov, Particle-hole duality in the lowest Landau level, Phys. Rev. Lett. 118, 206602 (2017).
  • Lee et al. [2017] C. H. Lee, M. Claassen, and R. Thomale, Band structure engineering of ideal fractional Chern insulators, Phys. Rev. B 96, 165150 (2017).
  • Lee et al. [2013] C. H. Lee, R. Thomale, and X.-L. Qi, Pseudopotential formalism for fractional Chern insulators, Phys. Rev. B 88, 035101 (2013).
  • Lee et al. [2015] C. H. Lee, Z. Papić, and R. Thomale, Geometric construction of quantum Hall clustering hamiltonians, Phys. Rev. X 5, 041003 (2015).
  • Claassen et al. [2015] M. Claassen, C. H. Lee, R. Thomale, X.-L. Qi, and T. P. Devereaux, Position-momentum duality and fractional quantum hall effect in Chern insulators, Phys. Rev. Lett. 114, 236802 (2015).
  • Ma et al. [2013] Y.-Q. Ma, S.-J. Gu, S. Chen, H. Fan, and W.-M. Liu, The Euler number of Bloch states manifold and the quantum phases in gapped fermionic systems, EPL (Europhysics Letters) 103, 10008 (2013).
  • Ma [2020] Y.-Q. Ma, Euler characteristic number of the energy band and the reason for its non-integer values, arXiv:2001.05946  (2020).
  • Kobayashi [1987] S. Kobayashi, Differential geometry of complex vector bundles (Princeton University Press, 1987).
  • Bangert et al. [2009] V. Bangert, M. G. Katz, S. Shnider, and S. Weinberger, E7E_{7}, Wirtinger inequalities, Cayley 44-form, and homotopy, Duke Mathematical Journal 146, 35 (2009).
  • McDuff and Salamon [2017] D. McDuff and D. Salamon, Introduction to symplectic topology (Oxford University Press, 2017).
  • Hatcher [2002] A. Hatcher, Algebraic topology (Cambridge Univ. Press, Cambridge, 2002).