![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
We review recent progresses in the study of flat band systems, especially focusing on the fundamental physics related to the singularity of the flat band’s Bloch wave functions. We first explain that the flat bands can be classified into two classes: singular and non-singular flat bands, based on the presence or absence of the singularity in the flat band’s Bloch wave functions. The singularity is generated by the band crossing of the flat band with another dispersive band. In the singular flat band, one can find a special kind of eigenmodes, called the non-contractible loop states and the robust boundary modes, which exhibit nontrivial real-space topology. Then, we review the experimental realization of these topological eigenmodes of the flat band in the photonic lattices. While the singularity of the flat band is topologically trivial, we show that the maximum quantum distance around the singularity is a bulk invariant representing the strength of the singularity which protects the robust boundary modes. Finally, we discuss how the maximum quantum distance or the strength of the singularity manifests itself in the anomalous Landau level spreading of the singular flat band when it has a quadratic band-crossing with another band.
Graphical abstract
![](/cms/asset/5cdecd12-1446-4efd-ad14-b2cb3ca526eb/tapx_a_1901606_uf0001_oc.jpg)
1 Introduction
A flat band indicates a type of band structures with constant energy independent of the crystal momentum. The key feature of the flat band is that the charge carriers in it have a zero group velocity and an infinite effective mass. Due to the suppressed kinetic energy, flat band systems have been considered intriguing in various research areas [Citation1–6]. In condensed matter physics, flat bands are considered ideal to study many-body phenomena such as the ferromagnetism [Citation7–11], superconductivity [Citation12–16], and Wigner crystal formation [Citation17–19] because the kinetic energy of the carriers in the flat band is quenched and dominated by the electron–electron interaction. In bosonic systems such as the photonic crystals, the realization of the slow light via a flat band has gotten a great attention because it can be applied to enhance the light–matter interaction [Citation20–26].
Usually, we need a fine-tuning of the system parameters to have a flat band, and it has hampered the experimental realization of such a flat band. However, there have been considerable experimental efforts to synthesize real materials or artificial systems with a nearly flat band. In the condensed matter physics community, nearly flat bands are discovered in various kagome-type materials such as CoSn [Citation27,Citation28], Co3Sn2S2 [Citation29], FeSn [Citation30] Fe3Sn2 [Citation31,Citation32], and YCr6Ge6 [Citation33,Citation34], as well as pyrochlore oxides [Citation10,Citation35]. Also, by using the novel band structure engineering technique, so-called twistronics, many nearly flat band systems have been realized in various Moire superlattice systems composed of misaligned two-dimensional layers, such as the twisted bilayer graphene [Citation36–42], twisted bilayer transition metal dichalcogenides [Citation43,Citation44], and twisted multi-layer silicene [Citation45]. In particular, the twisted bilayer graphene at the magic angle [Citation46] has attracted a great attention due to its unconventional superconductivity [Citation36,Citation47–57] and Mott insulating phases [Citation36]. Furthermore, in artificial systems, such as the photonic lattices [Citation58–71], cold atom systems [Citation72,Citation73], engineered atomic lattices [Citation74], and metamaterials [Citation75–79], one can have more controllability of the system parameters than the fermionic systems, and indeed almost flat bands were realized.
To understand the condensed matter phenomena properly, we usually start from the non-interacting physics of the given system in the band theory level, and then examine the effects of various interactions, such as the electron–electron interaction, based on the knowledge of the non-interacting electronic properties. However, notwithstanding the rising interest in the interaction physics in flat bands discussed above, the nature of flat bands has not been well-understood even in the band theory level until recently.
In the viewpoint of the topological band theory, the flat band is simply considered trivial due to the following reasons. First, if the flat band is isolated from others and hopping processes are allowed within a finite range, the Chern number is always zero because the Bloch wave function is enforced to be analytic all over the Brillouin zone due to the flatness of the energy dispersion [Citation80,Citation81]. Second, even if the flat band is semimetallic by having a band-crossing point with another dispersive band, it cannot be a topological semimetal because the Berry phase obtained along a path enclosing the band-crossing point is usually unquantized and, even if it is quantized, its value is always (trivial) due to the band-flatness condition [Citation82]. Therefore, the flat band cannot exhibit the conventional bulk-boundary correspondence [Citation83–89] such as the chiral edge modes from the nonzero Chern number in Chern bands [Citation90–93], the Fermi arcs between two topological monopole charges with opposite signs in the Weyl semimetal [Citation94], and helical surface modes from z2 bulk topological invariant of topological insulators [Citation95–99].
Although a flat band is topologically trivial, it was recently revealed that it can exhibit intriguing geometric properties. To see what this means, let us consider the Hilbert–Schmidt quantum distance, which is defined as
where and
are normalized quantum states [Citation82,Citation100–102]. The quantum distance
, which takes a value between 0 and 1, measures how close the two states
and
are. In the case of flat band systems, we consider the quantum distance between Bloch eigenstates in momentum space. The quantum distance between two different Bloch states at momenta
and
is usually becomes zero when
because the two states take the identical form in this limit. However, when a flat band has a band crossing with another dispersive band, the Bloch wave functions of the flat band around the band-crossing point develop a nonzero quantum distance even if their momenta are close to each other [Citation81,Citation82]. In this case, we call such a flat band the singular flat band. Interestingly, it is recently shown that the maximum value of the quantum distance among all the possible pairs of the Bloch eigenstates around the band-crossing point is a bulk invariant measuring the strength of the singularity at the band-crossing point [Citation82]. Due to the singularity of the singular flat band, intriguing eigenmodes dubbed the non-contractible loop state and the robust boundary mode appear in systems with open boundaries and exhibit topological robustness in real space [Citation64,Citation81]. Moreover, the Landau levels of the singular flat band with a quadratic band-crossing show anomalous behavior, which is characterized by the maximum quantum distance
[Citation82].
These recent studies clearly demonstrate that the singular flat band is an intriguing new platform where the geometrical properties of Bloch states can be studied. In general, the most crucial notions consisting of geometry are the curvature and metric. While the curvature indicates how much a given manifold is curved, the metric measures the distance between two infinitesimally close points on it. In solids, the Berry curvature and the quantum metric take the roles of the curvature and metric, respectively, in the geometric description of the Bloch wave functions in momentum space. Both the Berry curvature and the quantum metric are naturally defined from the quantum distance defined in (1). By assuming that the two momenta of the Bloch states are extremely close to each other, such that and
, one obtains
where is the quantum geometric tensor [Citation103–106] of the
th band given by
Here, is the periodic part of the Bloch eigenfunction, and
. In (2), the Einstein convention is used for
and
. Then the quantum metric is the real part of
, and the Berry curvature is the imaginary part of
[Citation107,Citation108]. While the curvature part has been studied exhaustively in past decades related to the semiclassical formulation of electron dynamics [Citation109] and the topological classifications of materials [Citation110], physical implications of the metric and quantum distance have been much less investigated although there have been several studies on the physical roles of them in the current noise [Citation111], the orbital magnetic susceptibility [Citation112,Citation113], the steady-state Hall response [Citation114], the superfluid weight of a nearly flat Chern band [Citation115,Citation116], etc. Very recently, it was reported that the quantum metric can be measured experimentally [Citation117,Citation118]. It is worth mentioning that in the case of the singular flat band with a quadratic band crossing, the quantum metric around the crossing point is directly connected to the maximum quantum distance
via a relation
, where
is the polar angle with respect to the crossing point and
is the quantum metric in the polar coordinate [Citation119]. To achieve a unified understanding of the geometric nature of solids, further studies on the role of the quantum metric and quantum distance in various condensed matter phenomena have to be performed. The singular flat band model can be an ideal starting point of such research direction because this model’s band-crossing point is characterized by the quantum distance instead of Berry curvature.
In this article, we pedagogically review the recent progress in the study of singular flat bands. In particular, we focus on the novel geometrical properties arising from the fundamental relation between the singularity of the band crossing and the maximum quantum distance of a singular flat band. In Section 2, we explain the relation between the existence of the singularity and the completeness of compact localized states, characteristic localized eigenstates of the flat band. Here we show that the nonzero maximum quantum distance guarantees the existence of non-local eigenstates, called the non-contractible loop state and the robust boundary mode, which induce intriguing topological features in real space. In Section 3, we briefly introduce the recent progress in the photonic lattice research, and then describe how the novel real-space topology was probed experimentally in the photonic lattices. In Section 4, we show how to define the strength of the singularity of the singular flat band via the maximum quantum distance or the pseudospin structure. In Section 5, we explore how the maximum quantum distance is manifested in the exotic Landau level structure of the singular flat band with a quadratic band-crossing, which will pave the way for the measurement of the quantum distance of solids. Here we also discuss how this phenomenon is related the divergence of the orbital magnetic susceptibility. Then, we show several candidate realistic materials, where those results can be observed. Finally, in Section 6, we give conclusions and outlook of the research on the singular flat band systems.
2 Flat band and singularity
Figure 1. (a) A compact localized state (CLS) and two non-contractible loop states (NLSs) of the nearest neighbor tight-binding model on the kagome lattice. The red dashed box represents the unit cell located at . (b) The bilayer square lattice. Hopping parameters corresponding to the black and green dashed arrows are
, and those for the red arrows are
. (c) A compact localized state and two non-contractible loop states of the Lieb lattice model with only the nearest neighbor hopping processes. Band structures for the models in (a), (b), and (c) are represented in the lower panels (d), (e), and (f), respectively. All the figures are adapted from [Citation81]
![Figure 1. (a) A compact localized state (CLS) and two non-contractible loop states (NLSs) of the nearest neighbor tight-binding model on the kagome lattice. The red dashed box represents the unit cell located at R=(0,0). (b) The bilayer square lattice. Hopping parameters corresponding to the black and green dashed arrows are 1/2, and those for the red arrows are −2. (c) A compact localized state and two non-contractible loop states of the Lieb lattice model with only the nearest neighbor hopping processes. Band structures for the models in (a), (b), and (c) are represented in the lower panels (d), (e), and (f), respectively. All the figures are adapted from [Citation81]](/cms/asset/9b0c782e-8631-41b5-b2d1-948dc098b43f/tapx_a_1901606_f0001_oc.jpg)
2.1 Compact localized state
The group velocity, the first derivative of an energy dispersion with respect to crystal momentum, represents the mobility of the charge carriers in the band. In the flat band, the group velocity is strictly zero for all momenta in the Brillouin zone. This implies that charge carriers in the flat band are immobile or localized. The localized nature of the flat band is explained by the existence of a special localized eigenstate called the compact localized state, whose amplitude is finite only inside a finite region in real-space, and exactly zero outside of it [Citation1,Citation4,Citation81,Citation120–125]. As a representative example, the compact localized state of the kagome lattice with the nearest neighbor hopping processes is shown in ). One can easily check that this compact localized state is an eigenstate of the kagome lattice model by applying the corresponding Hamiltonian operator to it [Citation120]. Due to the destructive interference provided by the special lattice structure of the kagome lattice, which is also related to the geometrically frustrated lattice structure, the amplitudes of the compact localized state do not leak from its boundary, and the compact localized state recovers its original form after the hopping processes [Citation120]. Note that the compact localized state is not unique, and one can have various types of the compact localized states by linear combinations of the smallest compact localized states centered at different positions [Citation81]. Although the compact localized state looks similar to the Wannier function, they are quite different in the sense that (i) compact localized states usually do not form an orthonormal set while Wannier functions do, (ii) a compact localized state is an eigenmode while Wannier function is not usually.
A compact localized state is guaranteed to exist when the system has a flat band and the corresponding Hamiltonian is described by a tight-binding model with a finite hopping range [Citation81]. If the hopping processes are allowed within a finite range, each matrix element of the Bloch Hamiltonian of the tight-binding model of such system is in the form of the finite sum of Bloch phases. That is, the
th component of
is given by
where is the dimension of the system described by
primitive vectors
,
,
is an integer ranging from an integer
to another integer
, and
is a complex coefficient of the Bloch phase
. In a simple manner, we denote the Bloch phase as
, where
is the lattice vector. Then, the eigenvector
for the flat band with energy
is obtained by solving coupled homogeneous equations given by
. Since all the matrix elements of
are also in the form of the finite sum of Bloch phases, one can always find an unnormalized solution for the
in the form of the finite sum of Bloch phases. By using the fact that any arbitrary linear combination of the Bloch wave functions in the flat band is also an eigenstate, one can construct an eigenstate expressed by
where is the Bloch eigenstate of the flat band,
is the normalization constant for
. While
is just an arbitrary smooth function of crystal momentum, we show below that
becomes a compact localized state if
is chosen to be a multiplication between
and an arbitrary finite sum of Bloch phases. Note that the Bloch eigenstate is in the form
, where
is the total number of the unit cells in the system,
is the
th component of the normalized eigenvector of the flat band given by
,
represents the
th orbital in the unit cell located at
, and
is the size of the Bloch Hamiltonian. As a result, the amplitude of the
th orbital of
in the unit cell at
is obtained as
Since is in the form of the finite sum of Bloch phases, the value of
vanishes if
is out of the range of the lattice vectors in the Bloch phases in
. This implies that
is a compact localized state, and one can have various kinds of the compact localized state depending on the choice of
. This is also a systematic scheme to obtain the compact localized states from the Bloch eigenstate of the flat band.
Let us consider the kagome lattice model as an example. Its Bloch Hamiltonian is given by
where ,
, and
. If the unnormalized eigenvector for the flat band of this model is chosen to be
, a compact localized state centered at
is obtained from (6) as
This is the compact localized state with the smallest size in the kagome lattice model as plotted in ). The th element of the column vector on the right-hand side of (8) represents the amplitude of the compact localized state at the site-
, which is indicated by three colors in ). For example, the first element
of this vector means that the amplitudes at the site-1 in the unit cells located at
and
are the coefficients of
and
, namely −1 and 1, respectively. In fact, these coeffecients can be inferred from the unnormalized form of the Bloch eigenstate
because there is a one-to-one correspondence between the coefficients of
in
and that of
in (8). In general, one can easily derive the compact localized state from the unnormalized Bloch eigenstate in this way.
2.2 Singular and non-singular flat bands
One fundamental question on the compact localized state is that can one find a complete set of compact localized states to span the whole flat band? [Citation81,Citation120] We need to answer this question to understand the localized nature of the flat band completely. It seems quite obvious that
compact localized states centered at
different unit cells are linearly independent. Here
indicates the total number of unit cells in the system with a periodic boundary condition. However, it was noted that, under the periodic boundary condition, those
compact localized states fail to form a complete set sometimes [Citation120]. Later, it was rigorously shown that one cannot find a complete set of compact localized states if the Bloch eigenstate of the flat band is discontinuous at a momentum in the Brillouin zone [Citation81]. This singular point can only be generated by the band-crossing of the flat band with another dispersive band, and the Bloch eigenstate in an isolated flat band is always analytic over the whole Brillouin zone. This is consistent with the fact that the Chern number of an isolated flat band is always zero when the hopping range of the model is finite [Citation80]. A flat band with the singularity is called a singular flat band, and otherwise, it is called a non-singular flat band [Citation81].
The flat band of the kagome lattice model in Sec. 2.1 is an example of the singular flat band. First, note that the flat band has a quadratic band-crossing at point (
). We check whether the corresponding Bloch eigenstate is singular or non-singular at the band-crossing point as follows. The normalized eigenvector of
in (7) is given by
where . Note that all the elements of the unnormalized eigenvector
vanish at
point simultaneously. Since the normalization factor
is also zero at
point, every element of the normalized eigenvector
is in the form of zero over zero. Therefore, the value of
depends on the choice of the path approaching
point. For instance, one can show that
and
. This implies that
is discontinuous at
point, and therefore the flat band of the kagome lattice is singular.
Let us consider another example of the flat band with a band-crossing, the bilayer square lattice model described in ). In this model, only the nearest-neighbor hopping processes with amplitude are allowed between the sites in the same layer, and the inter-layer hopping parameters corresponding to the red and green dashed lines are
and
, respectively. The corresponding Hamiltonian is given by
Although the flat band has a band touching at point, it is a non-singular flat band because its eigenvector
is just a constant and therefore analytic for all momenta.
As shown by these two examples, in general, we have to investigate the Bloch eigenstates of a flat band to determine whether the flat band is singular or not. However, we note that a flat band of any one-dimensional system is always non-singular irrespective of the presence or absence of band crossings [Citation81].
The most crucial result of the existence of the singularity in a flat band is that one cannot find a complete set of compact localized states spanning a singular flat band under the periodic boundary condition [Citation81]. Kagome lattice model is a representative example where the compact localized states form an incomplete set [Citation120]. To show this, we denote number of different compact localized states (labeled by
) by
which have different shapes and center positions depending on the choice of and
respectively. Then, a determinant diagnosing their linear independence or dependence is given by
because the Bloch wave functions ’s with different momenta in (11) form a complete set. As noted in Sec. 2.1, to ensure the compact localization of
,
should be proportional to
, where
is a unnormalized eigenvector of the flat band in the form of the finite sum of Bloch phases. Therefore, one can note that
As a result, if is zero at a crystal momentum, namely the normalized Bloch eigenstate
is singular at that momentum,
also becomes zero implying that
compact localized states are always linearly dependent of each other.
An important point here is that the zeros of the unnormalized eigenvector appear only at the band-crossing points in momentum space. Note that we only consider the zeros which cannot be removed by factoring out a function vanishing at that momenta from all the components of the unnormalized eigenvector because any removable zeros are irrelevant to the completeness of the compact localized states. This is justified as follows. For simplicity, let us consider a Hamiltonian with a real-valued flat band eigenvector. In this case, the unnormalized eigenvector can be represented as
where is a real-valued finite sum of Bloch phases because
is an eigenvector of the flat band. We assume that
has an irremovable zero at
. Since any finite sum of Bloch phases is analytic for all momenta,
can be expanded in the vicinity of
such that
, where
and
is the polar angle with respect to
. Then, two the components of the normalized eigenvector around
becomes
and
respectively. Note that one can have different normalized eigenvectors at depending on the polar angle because the ratio between
and
is different from that between
and
due to the irremovability of the zero at
. As
can take an arbitrary value, it may imply that there are infinitely many different eigenvectors at the flat band’s energy depending on
. Such a situation occurs generally when there is a degeneracy. This is because any linear combination of the degenerate eigenvectors also becomes an eigenvector, which is different from the case of non-degenerate point at which the eigenvectors can be defined uniquely up to the overall U(1) phase factor. Therefore, there should be a degeneracy or the band-crossing at the zeros of the unnormalized eigenvector. One can show that the same conclusion applies to the larger Hamiltonians and complex flat band eigenvector. Therefore, if we want to determine whether the flat band is singular or not, we only have to investigate the Bloch wave functions at the band-crossing points.
As one can note from (11), the Bloch wave function corresponding to the singular momentum is always absent in constructing the compact localized state because is zero at that momentum. Namely,
compact localized states of a singular flat band are made of a smaller (less than
) number of Bloch wave functions, which implies that they cannot be linearly independent. For example, in the kagome lattice, the
compact localized states obtained by translating the minimal compact localized state in (8) by all the lattice vectors are known to be linearly dependent of each other because a simple sum of all of them vanishes. On the other hand, if the flat band is non-singular, namely there exists a nonzero
for the flat band, one can obtain non-zero
by choosing
and
to be
distinct lattice vectors. Therefore, one can find a complete set of
compact localized states for a non-singular flat band.
We note that there is a subtle issue related to the finite size effect. Let us denote the number of unit cells along direction by
, so that the total number of the unit cells in the system is given by
. Then, the crystal momentum discretized under periodic boundary condition along each reciprocal lattice vector is
, where
is an integer. Therefore, even though the continuum expression of the band dispersion shows a singular band-crossing at a certain crystal momentum
, one can find a complete set of
compact localized states if
cannot be represented by the form
for the given system size
. Let us consider the Lieb lattice model as an example. Considering only the nearest neighbor hopping processes, the Hamiltonian is given by
The eigenvector of its flat band is obtained as
which is singular at . Note that this
is allowed only when both
and
are even numbers. If we choose a compact localized state for this model as shown in ), any linear combination of the
translated copies of this compact localized state is nonzero for odd values of
and
. This implies that these compact localized states form a complete set. On the other hand, if
and
are even integers, the same staggered sum of them vanishes. In this case, these compact localized states are linearly dependent on each other. A similar situation can be found in
-
model too, whose flat band has two triple degeneracy points [Citation126,Citation127]. Note that if the singular band-crossing is at
, we are free from this finite size issue in general.
2.3 Non-contractible loop and planar states
Since we cannot construct a complete set of basis wave functions for the singular flat band from the compact localized states as discussed in Sec. 2.2, some non-compact eigenstates should be included in addition to the compact localized states to span the flat band completely. In the kagome lattice, as an example, such non-compact modes were found to be the so-called non-contractible loop states, which are shown in ) [Citation120]. They are extended at least along one spatial direction while compactly localized along other directions [Citation120]. One can check that there are two such extended eigenmodes in the kagome lattice under the periodic boundary condition, which are independent of the compact localized states [Citation120]. Any other non-contractible loop states at different positions and shapes can be constructed by a linear combination of the given non-contractible loop states and the compact localized states [Citation120]. In the kagome lattice case, it was noted that number of compact localized states can be linearly independent of each other [Citation120]. Therefore, the existence of these two additional non-contractible loops states is consistent with the fact that there are
number of degenerate Bloch eigenstates at the energy of the flat band, which indicates that the flat band should be touching with another dispersive band at a momentum [Citation120].
The non-contractible loop states exhibit intriguing topological aspects in real space as follows. For a two-dimensional system under periodic boundary condition, the two independent non-contractible loop states form two closed loops encircling the torus geometry along the toroidal and poloidal directions, respectively [Citation120]. Also, they are robust because a non-contractible loop state cannot be cut in the middle by any linear combination between the non-contractible loop state and compact localized states, which merely deforms the shape of the non-contractible loop state [Citation64,Citation81,Citation120].
In general, it was shown that such non-contractible loop states are guaranteed to exist in a two-dimensional singular flat band [Citation81]. If the band-crossing point is at , the Bloch wave function corresponding to this momentum does not participate in constructing a compact localized state by the linear combination of Bloch wave functions as discussed previously via Equationequation (11)
(11)
(11) . As a result, any kind of eigenmodes including the Bloch wave function at
is linearly independent of the compact localized states. While the Bloch wave function at
is extended along both directions
and
, one can make it compactly localized along one of those two directions by a linear combination of it with other Bloch wave functions as follows. Let
be the Hamiltonian under consideration. If we fix
to
, where
, we obtain an effective one-dimensional flat band Hamiltonian
with momentum
. Since any one-dimensional flat band is a non-singular flat band as noted previously, one can always find a compact localized state for this effective Hamiltonian localized along
direction, which includes the Bloch wave function at
too. This is the non-contractible loop state of the full Hamiltonian
because it is still extended along
direction with momentum
. In the similar way, one can find another non-contractible loop state extended along
direction too.
In three-dimensional case, one can have a different kind of partially extended states with nontrivial real-space topology, called the non-contractible planar states (NPSs) [Citation81]. If we denote a singular flat band Hamiltonian for this case by , one should fix two components of
to the corresponding components of the singular momentum
to obtain an effective one-dimensional flat band Hamiltonian in this case. Let us assume that those two fixed momenta are
and
, for convenience. Then, we can construct a planar eigenstate, namely a NPS, which is extended along two directions
and
, and compactly localized along
. By fixing other pairs of momentum components, one can obtain other NPSs along different directions.
Figure 2. (a) The robust boundary mode of the kagome lattice model. This can be obtained by adding all the compact localized states centered at all the unit cells (from to
). If one subtract a compact localized state centered at
to try to cut the robust boundary mode, it cannot be cut but deformed as illustrated in the right-hand side. (b) A boundary eigenmode on the bilayer square lattice, which is obtained by adding compact localized states centered at boundary sites. This boundary mode can be cut by subtracting a compact localized state on the boundary as shown in the right-hand figure. In (c), it is illustrated schematically how the non-contractible loop states are connected to the robust boundary mode by deforming the periodic boundary condition to the open boundary condition. All the figures are from [Citation81]
![Figure 2. (a) The robust boundary mode of the kagome lattice model. This can be obtained by adding all the compact localized states centered at all the unit cells (from R1 to R9). If one subtract a compact localized state centered at R2 to try to cut the robust boundary mode, it cannot be cut but deformed as illustrated in the right-hand side. (b) A boundary eigenmode on the bilayer square lattice, which is obtained by adding compact localized states centered at boundary sites. This boundary mode can be cut by subtracting a compact localized state on the boundary as shown in the right-hand figure. In (c), it is illustrated schematically how the non-contractible loop states are connected to the robust boundary mode by deforming the periodic boundary condition to the open boundary condition. All the figures are from [Citation81]](/cms/asset/dc41bc6b-7b57-4412-a5bd-46d8810eaffd/tapx_a_1901606_f0002_oc.jpg)
2.4 Robust boundary modes
Under the open boundary condition, the singularity of a singular flat band is manifested by the existence of a boundary eigenmode, called the robust boundary mode [Citation81]. The robust boundary mode of the kagome lattice is shown in ) as an example. In fact, the robust boundary mode is closely related to the non-contractible loop state because the non-contractible loop states can transform to the robust boundary mode as we modify the periodic boundary condition to the open boundary condition as follows. We first prepare four non-contractible loop states on a torus geometry of the system respecting the periodic boundary condition, two along the toroidal direction and the other two along the poloidal direction. Then, we transform the torus into a finite system with the open boundary condition by cutting the surface of the torus between non-contractible loop states as illustrated in ). In this process, the four non-contractible loop states lead to the robust boundary mode. Since the appearance of the non-contractible loop states is from the singularity of the singular flat band, the existence of the robust boundary mode is also due to the same singularity.
This is a new kind of bulk-boundary correspondence between the singularity in the bulk system and the boundary mode of the finite system. The conventional bulk-boundary correspondence in topological systems is the correspondence between a bulk topological invariant and the in-gap boundary modes of the finite system. Examples of the topological invariant are the Chern number of the Chern bands or Landau levels [Citation90–93], Z index in topological insulators [Citation95–99], the Zak phase of the one dimensional insulators with reflection symmetry [Citation88,Citation128–130], and so on. The robust boundary mode of the singular flat band is distinguished from the topological boundary modes in the following perspectives. First, the robust boundary mode is not in the bulk band gap but at the same energy of the flat band. When we consider the periodic boundary condition, a linear combination of all the compact localized states can vanish in the case of the singular flat band. On the other hand, under open boundary condition, the same linear combination of compact localized states vanishes in most parts, except at the boundary of the system. Interestingly, the states on the boundary cannot be removed via any linear combination of compact localized states. In fact, the linear combination of the compact localized states, which vanishes under periodic boundary condition, gives the robust boundary modes under open boundary condition. Also, as the robust boundary mode is obtained by the linear combination of the compact localized states, it has the same energy as the compact localized state, which is identical to the flat band’s energy. Therefore, one cannot detect the robust boundary mode by electronic spectroscopy such as the angle-resolved photoemission spectroscopy (ARPES) [Citation131]. Instead, one should use other kinds of experimental tools, which can access the wave function directly. The photonic lattice set-up is appropriate for this because here wave functions can be prepared as well as detected quite freely. Second, the robust boundary mode is protected by the singularity and the flatness of the band while the topological edge states are protected by the topological charge or certain crystal symmetries. Due to the existence of the singularity, one can find a linear combination of compact localized states, which vanishes under periodic boundary condition but gives the robust boundary modes under open boundary condition. Also, because of the flatness, the robust boundary mode can also be an eigenmode as explained above Third, robust boundary mode is linearly dependent of the compact localized states of the bulk, while the topological edge modes are not. Note that compact localized states are linearly independent of each other under the open boundary condition, and form a complete set spanning the flat band as discussed above. Since the robust boundary mode is obtained from a macroscopic number, corresponding to the system size, of the compact localized states, it cannot be cut or destroyed by adding a few number of compact localized states, and this is what we mean by the robustness. Even in the non-singular flat band, one can construct a boundary eigenmode by collecting all the compact localized states spanning the boundaries of the finite system. However, such boundary mode is not the robust boundary mode because it can be cut by just adding a compact localized state near the boundary as shown by the square lattice bilayer example in ). Let us note that the presence of robust boundary modes is confirmed in a recent experiment in photonic lattices as discussed in Section 3.
Figure 3. (a) A two-dimensional flat band with a quadratic band touching (QBT). (b) If the QBT is non-singular, one can gap out the QBT while maintaining the band flatness. When the flat band is shifted upward in the energy space, one obtains a line touching between the flat band and the parabolic band. (c) If the band-crossing point is singular, one should break the flatness of the flat band during the gap-opening process of the QBT point. If the flat band is mixed with the quadratic band, one obtain point touching between two bands. This figure is adapted from [Citation81]
![Figure 3. (a) A two-dimensional flat band with a quadratic band touching (QBT). (b) If the QBT is non-singular, one can gap out the QBT while maintaining the band flatness. When the flat band is shifted upward in the energy space, one obtains a line touching between the flat band and the parabolic band. (c) If the band-crossing point is singular, one should break the flatness of the flat band during the gap-opening process of the QBT point. If the flat band is mixed with the quadratic band, one obtain point touching between two bands. This figure is adapted from [Citation81]](/cms/asset/712cd7eb-d15e-43fc-a876-d3a4e88a9438/tapx_a_1901606_f0003_b.gif)
2.5 Flatness enforced band-crossing
Conventionally, a band-crossing is protected by crystalline symmetries or topological charge [Citation132]. Interestingly, in the case of the singular flat band, the protection mechanism of its singular band-crossing does not belong to any of them. Instead, it was shown that it is protected by just the flatness of the flat band. One can analyze the generic protection mechanism for the band-crossing of the flat band by using a low-energy continuum Hamiltonian [Citation81]. For simplicity, we consider two-dimensional systems with two bands touching each other. If the flat band has a point touching with another band, the band-crossing type should be at least parabolic, because the energy dispersion of the band touching with the flat band becomes not differentiable in the case of the linear band-crossing. The most general quadratic form of the two-band Hamiltonian is given by
where ,
,
are Pauli matrices, and
is the identity matrix. Any quadratic Hamiltonian matrix can be unitarily transformed to the above form. The condition for the existence of a flat band is given by
where we assume that the flat band is at the zero energy without loss of generality. From this we obtain five constraints as follows.
Under these constraints, the Hamiltonian yields a flat band and a parabolic band touching each other at
. One can note that the flat band is a singular flat band only when both
and
are nonzero. The non-singular flat band and singular flat band show distinct features during the gap-opening process as follows.
First, the non-singular flat band can be decoupled from the parabolic band while maintaining the flatness of the energy dispersion. In this case, it was shown that the Hamiltonian matrix can be always transformed to
. Since
has no off-diagonal terms, a perturbation of the form
can control the energy of the flat band without destroying its flatness. As a result, we obtain a flatness-preserved gap opening or a line toughing between two bands when
. This is summarized in .
On the other hand, in the case of the singular flat band, the flat band cannot be detached from the parabolic band while maintaining its flatness when a mass term of the form is added to the system. From the flatness conditions,
and
, we always obtain
if one of the mass parameters is nonzero. This contradicts with the fact that
should be nonzero to have a singular flat band. Therefore, it is impossible to shift the energy of the singular flat band while preserving its flatness, which implies that the singular flat band’s band-crossing is enforced by the band’s flatness. Depending on the choice of the mass parameters, one can either have the gap-opening at the band-crossing point or observe the generation of two Dirac points at different crystal momenta. These results are illustrated in ). Moreover, one can find a mass term, which results in a Chern band separated from other bands. As an example, in the kagome lattice, if we add a gap-opening perturbation
, where
’s are the Gell–Mann matrices, to
, we have a nearly flat band with a nonzero Chern number. This can be a systematic route to obtain a nearly flat Chern band starting from a singular flat band when one uses a simple scheme to design a singular flat band model proposed in [Citation81].
3 Observation of the real-space topology in photonic lattices
3.1 Flat bands in photonic lattices
One can fabricate a two-dimensional periodic array of optical waveguides, so-called the photonic lattice [Citation133–148], in a three-dimensional medium by using the femtosecond laser writing technique [Citation149,Citation150]. A light beam propagating in this system satisfies a wave equation given by
within the paraxial approximation [Citation134,Citation139,Citation140], where describes the electric field envelope of the light propagating along
-axis,
and
are the refractive indices outside and inside the waveguide of the medium,
is the wavenumber in vacuum, and
. This can be regarded as a 2D Schrödinger equation by substituting
with time, and
with the periodic potential. Namely, the 2D cross section of the medium perpendicular to the
-axis corresponds to a 2D lattice model in solid-state physics. In the photonic lattice set-up, one can have more freedom to design lattice models, control band parameters of them, and prepare an initial state compared with the conventional condensed matter experiments. As a result, many flat band systems such as the Lieb [Citation58–61] and kagome lattice models [Citation62–64] have been studied by using the photonic lattice experimental set-ups [Citation65–71].
In the photonic lattice, a flat band is probed indirectly by showing that a compact localized state is an eigenmode because a compact localized state is guaranteed to exist in a flat band system and represents the localized nature of this dispersionless band. If the waveguides embedded in a medium extend along -direction from
to
, an incident light beam, whose amplitudes and phases correspond to those of the compact localized state, enters into the waveguides at
. As the light beam propagates along
-axis, it would also spread along the transverse direction in general due to the coupling between neighboring waveguides. However, the light beam corresponding to the compact localized state remains the same up to an overall phase factor during the propagation because the compact localized state is an exact eigenstate of the flat band system. As a result, what we observe at the other end of the waveguides (
) is the same incident light beam only multiplied by a phase factor, and this is how one demonstrates the existence of a flat band in the photonic lattice experiments.
Figure 4. (a) The kagome lattice in the Corbino geometry. (b) Experimental realization of (a) in the photonic lattice. (c) Observation of the non-contractible loop state. With the out-of-phase configuration, the initial state propagates along -direction without any transverse dissipation. (d) An initial state with the in-phase condition is not an eigenmode. (e) A finite kagome lattice with the open boundary condition. (f) Experimental realization of the robust boundary mode. The robustness of the robust boundary mode is tested by adding (g) a compact localized state, and (h) a defect. These figures are adapted from [Citation64]
![Figure 4. (a) The kagome lattice in the Corbino geometry. (b) Experimental realization of (a) in the photonic lattice. (c) Observation of the non-contractible loop state. With the out-of-phase configuration, the initial state propagates along z-direction without any transverse dissipation. (d) An initial state with the in-phase condition is not an eigenmode. (e) A finite kagome lattice with the open boundary condition. (f) Experimental realization of the robust boundary mode. The robustness of the robust boundary mode is tested by adding (g) a compact localized state, and (h) a defect. These figures are adapted from [Citation64]](/cms/asset/08d7dd50-6967-4bfe-b6cb-c7106a23d035/tapx_a_1901606_f0004_oc.jpg)
3.2 Observing non-contractible loop states and robust boundary modes in photonic lattices
As described in the previous section, non-contractible loop states and robust boundary modes should also be considered together with compact localized states to understand the fundamental properties of a flat band completely. Probing one of those states is equivalent to the observation of the Bloch wave function’s singularity of the flat band. The first experimental attempt to observe non-contractible loop states was done in the photonic Lieb lattice [Citation61]. They considered a finite Lieb lattice with an open boundary condition, and showed that an extended state, called the line state is an eigenmode. Although the line state looks similar to the non-contractible loop state, it is an indirect demonstration of the non-contractible loop state because the stability of the line state relies on the shape of the boundary of the finite system. Moreover, in the case of the kagome lattice, one cannot find a stable line state for such a simple open boundary condition, and we need a more delicate decoration of the boundaries for the line state to become an eigenstate [Citation64]. In principle, it is impossible to observe the non-contractible loop state directly in a finite system because the non-contractible loop state is well-defined and independent of compact localized states only under the periodic boundary condition.
The periodic boundary conditions in two dimensions are usually hard to be implemented in real space simultaneously. However, we need just one periodic boundary condition along one direction to stabilize one of the non-contractible loop states of a flat band. This was realized experimentally by arranging the waveguides of the photonic crystal in an annular disk, called the Corbino disk, where the spatial periodicity is preserved along the azimuthal direction. The kagome lattice in the Corbino geometry illustrated in ) was realized in the photonic lattice as shown in ). Then, as in ) and , only the light beam under out-of-phase condition, corresponding to the non-contractible loop state, propagates without dissipation along the transverse direction, while the light beam under in-phase condition is destroyed after the propagation. The demonstration of the stability of the light beam winding the whole system with out-of-phase configuration is the first direct observation of the real-space topology represented by the non-contractible loop state.
Another topological object manifesting the singularity of the flat band is the robust boundary mode. The robust boundary mode only requires the open boundary condition, which is more realistic than the periodic boundary condition. As shown in ), a finite kagome lattice is synthesized in a photonic lattice. Then, a state spanning the boundary of the system with out-of-phase configuration is found to be an eigenmode corresponding to the robust boundary mode as presented in ). Moreover, the robustness of the robust boundary mode is also checked by showing that the robust boundary mode cannot be cut, although it is slightly deformed, when a defect is placed on the boundary of the system.
4 Strength of the singularity
4.1 Maximum quantum distance
In topologically non-trivial bands, the bulk topological invariant is nonzero due to a singular behavior of the corresponding Bloch wave function in momentum space. For example, a Chern band yields a nonzero Chern number as a bulk topological number because the Bloch wave function of this band cannot be determined uniquely and smoothly over the whole Brillouin zone [Citation151]. In the case of the Weyl semimetal, the nonzero quantized monopole charge of the Weyl point is from the discontinuity of the Bloch wave function at this point [Citation94]. In two-dimensional band-crossings like the Dirac point of graphene, the Berry phase given by
is usually used to determine whether the band-crossing point is topological or not, where is a closed path enclosing the band-crossing point and
is the cell-periodic part of the Bloch wave function [Citation152,Citation153]. Here, the singularity of the Bloch wave function at the touching point also plays a key role when gamma is nonzero modulo
and quantized.
On the other hand, the band-crossing point of the singular flat band is always topologically trivial due to the flatness condition of the band although the Bloch wave function is singular there [Citation82]. First, the generic singular flat band with a quadratic band crossing does not have a quantized Berry phase (27) because all the Pauli matrices are used in the generic Hamiltonian in (19). Even though the Berry phase is quantized due to some symmetries such as and
, where
,
, and
are twofold rotation about the z-axis, inversion, time-reversal symmetries, respectively, the Berry phase of the singular flat band reads
with integer
. These imply that the band-crossing point of the singular flat band generally does not carry a quantized topological charge, and it is considered topologically trivial.
It was recently proposed that the quantum distance can be used to define a bulk number characterizing the singular band-crossing point of the singular flat band, representing the strength of the singularity [Citation82]. From (1), the quantum distance between two Bloch states at momenta and
is given by
. Then, the strength of the singularity is defined by the maximum value of the quantum distance between all the possible pairs of Bloch eigenstates of the flat band around the singular point. Namely,
where and
are momenta close to the band-crossing point. By definition, the strength of the singularity
is valued from 0 to 1. In the case of the non-singular flat band, the Bloch eigenmode of (19) is always independent of momentum, so that the quantum distance is zero between an arbitrary pair of momenta. As a result,
for the non-singular flat band. On the other hand, in the case of the singular flat band,
remains finite even if we consider momenta extremely close to the band-crossing point. Note that the eigenvector of the Hamiltonian (19) is independent of
, which is the magnitude of the momentum displacement with respect to the band-crossing point, because all the dependencies on
can be factored out from this quadratic Hamiltonian. The general form of the maximum quantum distance of the Hamiltonian (19) is evaluated as
where and
are the minimum and maximum effective masses of the quadratic band touching with the flat band.
The strength of the singularity, , also represents how strong the inter-band coupling between the flat band and parabolic band as explained below by a counting argument [Citation82]. First, the number of independent parameters of the Hamiltonian (19) is four because five parameters among nine parameters in (19) are determined by others through the flatness conditions from (21) to (25). Since we need three parameters to describe the dispersion relation of the quadratic band of the form
, and no parameters for the flat band fixed at the zero energy, the remnant one parameter should be related to the inter-band coupling. Indeed, it was shown that the general flat band model (19) can be represented by four parameters
,
,
, and
, where
is the mass tensor for the parabolic band. This implies that
is very natural inter-band coupling parameter of the two-dimensional flat band model with a quadratic band touching.
As an example, let us consider the kagome lattice. Its low energy effective Hamiltonian around the quadratic band-crossing is given by
where the flat band’s energy is shifted to zero [Citation81]. The eigenvector of the flat band of this model is evaluated as
where is the polar angle with respect to the band-crossing point at
. Then, the quantum distance between two eigenvectors at
and
is obtained as
This gives us the maximum quantum distance when
, which is consistent with the formula (29) because the band parameters for
are given by
,
, and zero for other parameters. Since the quantum distance is ranging from 0 to 1 by definition, the kagome lattice’s flat band corresponds to the maximally singular case.
4.2 Maximum pseudospin canting angle
The singular band-crossing point of the flat band can be also characterized by the canting structure of the pseudospin given by
where is the eigenvector of the flat band and
is the coefficient of
in (19) given by
,
, and
[Citation82]. One can represent the pseudospin on a surface of a unit sphere because
,
, and
are all in quadratic order of
so that the size of the pseudospin is constant. Pseudospin texture is drawn on a circular path encircling the band-crossing point as illustrated in . Since the eigenvector of (19) is independent of
, the pseudospin structure is also independent of the size of circle. In the case of the non-singular flat band, the pseudospins are aligned to each other as shown in ). As a result, the pseudospin is well-defined even at the band-crossing point. However, in the case of the singular flat band, we have a canting structure of the pseudospins around the band-crossing point as illustrated in ) and . This implies that the pseudospin cannot be fixed uniquely at the touching point, which is consistent with the fact that the corresponding Bloch eigenstate is not uniquely defined there.
It was shown that the quantum distance between two Bloch states at and
can be represented by the relative angle
between two pseudospins at the same momenta as follows.
which represents the one-to-one correspondence between those two quantities. Therefore, one can use the maximum canting angle as the strength of the singularity of the singular flat band instead of the maximum quantum distance
. According to the maximum canting angle, the value of the strength of the singularity ranges from 0 (non-singular flat band) to
(maximally singular flat band). In general, due to the topological triviality, the pseudospin texture of the singular flat band does not exhibit any winding structure except the maximally singular case (
or
), where we have a winding number 2 of the pseudospin around the band-crossing point. However, the singular flat band can be characterized by the pseudospin canting structure or the finite maximum quantum distance, which can be manifested in reality via Landau level structure as explained in the next section.
Let us consider again the kagome lattice as an example. For convenience, we use a unitary transformed Hamiltonian , where
. Then,
is described by
,
,
, and
. As a result, the pseudospin formula is given by
which shows double winding around the band-crossing point, and .
5 Anomalous Landau levels of singular flat bands
5.1 Onsager’s semiclassical theory
In solids, the band structure of electrons under a magnetic field, namely the Landau level structure, can be understood intuitively as well as quantitatively by using the Onsager’s semiclassical quantization rule described by
where represents energy,
is the area of a closed path at the energy
on a band,
is the magnetic field,
is the electric charge,
is the Planck constant,
is an integer corresponding to the Landau level index, and
is the quantum correction due to Berry phase, magnetic susceptibility, and higher-order magnetic responses [Citation152,Citation154–157]. Only the energy values of
satisfying the above equation are allowed as eigenenergies, and they are quantized because the Landau level index
is discretized. The key idea of this scheme is that one can predict the Landau level structure from the information obtained before applying magnetic field such as the band dispersion and geometric properties of the Bloch wave function. Conversely, one can extract the geometric properties of the Bloch wave function, such as the Berry phase, from the Landau level structure via the above formula. Let us consider graphene as an example. Denoting the Fermi velocity around its Dirac point by
, we have
and
. This leads to the well-known Landau level quantization
.
Figure 6. (a) A flat band and a parabolic band crossing each other at a point. and
represent the minimum and maximum effective masses. Generic shape of the Landau level structure of (b) the non-singular flat band and (c) the singular flat band. In (c),
is the Landau level spreading of the flat band. (d)
dependence of the anomalous Landau levels corresponding to the singular flat band. (e) The ratio
is plotted as a function of
. Diamond symbols are the numerical results from the tight-binding lattice model, and the red solid curve is the fitting function for the data
![Figure 6. (a) A flat band and a parabolic band crossing each other at a point. m1 and m2 represent the minimum and maximum effective masses. Generic shape of the Landau level structure of (b) the non-singular flat band and (c) the singular flat band. In (c), Δ is the Landau level spreading of the flat band. (d) 1/n dependence of the anomalous Landau levels corresponding to the singular flat band. (e) The ratio Δ/ℏωc is plotted as a function of dmax2. Diamond symbols are the numerical results from the tight-binding lattice model, and the red solid curve is the fitting function for the data](/cms/asset/cce032c3-771e-4bea-87b5-2137103e2c6c/tapx_a_1901606_f0006_oc.jpg)
5.2 Landau level spreading of the singular flat band with a quadratic band-crossing
The Onsager’s scheme cannot be applied to the flat bands because is not well-defined in this case. Namely, one can have infinitely many choices of the closed path at the energy of the flat band, and therefore
is not single-valued. What we can naively infer from the fact that the effective mass
is infinite in the flat band is that the cyclotron energy (
) is vanishing and all the Landau levels are trivially developed at the same energy of the flat band. However, it was recently shown that the Landau levels of a flat band are nontrivial and determined by a completely different scheme involving the quantum distance of the Bloch states of the flat band around the singular momentum, where the flat band is touching with another quadratically, as described below [Citation82].
To consider the effect of the magnetic field on the flat band, the transformations and
are applied to the quadratic band-crossing flat band Hamiltonian in (19), where
and
are the annihilation and creation operators, respectively, satisfying
. To ensure the Hermiticity,
is symmetrized by replacing it with
. Then, each element of (19) becomes an operator. This Hamiltonian is solved by applying a trial wave function of the form
where is the normalized eigenstate of the harmonic oscillator satisfying
and
, and
and
are complex coefficients.
The generic structure of the Landau levels of the flat band model (19) obtained by the procedure described above is shown in . The most intriguing feature is that the Landau levels are developed at the energies, where the density of states was zero before applying the magnetic field. This behavior is completely anomalous from the semiclassical point of view because the semiclassical quantization Equationequation (36)(36)
(36) does not have a solution when the energy
is in the band gap region. This phenomenon is called the Landau level spreading of a flat band [Citation82]. In the case shown in ), the energy of the flat band is lower than that of the parabolic band, the energy difference between the flat band and the lowest Landau level defines the magnitude of the Landau level spreading. On the other hand, if the flat band is above the hole-like quadratic band, the Landau level spreading is the energy of the highest Landau level with respect to the flat band’s energy.
While the Landau level spreading appears only in the singular flat bands, its magnitude can be shown to be determined by the maximum quantum distance of the Bloch eigenstate before applying magnetic field. More specifically, as shown in ), the ratio
is a simple monotonic function of
given by
where is the Landau level spreading and
is the cyclotron frequency of the parabolic band. It is important to note that
is only a function of
regardless of an arbitrary choice of band parameters of the model Hamiltonian in (19). This implies that the quantum distance plays the key role in the Landau level structure of the flat band, and the strength of the singularity of the band-crossing point manifests itself as the Landau level spreading. In fact, it was rigorously shown that two Hamiltonians with distinct sets of band parameters are equivalent to each other within a scale factor if they share the same
. The scale factor is the ratio between the cyclotron energies of the parabolic bands of those two Hamiltonians. Therefore, all the singular flat band systems with the same
share the same value of
. Moreover, the Landau level spreading can be understood as a result of the level repulsion between Landau levels originating from the flat band and the parabolic band. Indeed, the coupling strength between Landau levels from the flat and the parabolic bands, which is relevant to the level repulsion, were shown to be proportional to
in the weakly singular limit (
).
While the explicit form of is unavailable, its approximate form, which fits the numerical data accurately, was found to be
where ,
, and
. This was obtained from several analytic results such as
and
. By using this formula, one can extract the geometric information, namely the maximum quantum distance
, of the singular flat band around the band-crossing point from the Landau level structure.
Another intriguing feature of the Landau level structure of the singular flat band with a quadratic band-crossing is that the Landau levels are inversely proportional to the Landau level index close to the flat band’s energy as shown in ). Here, the energy of the flat band is assumed to be zero without loss of generality. As an example, in the case of the maximally singular (
) model in (30), the explicit formula for the Landau levels corresponding to the singular flat band is given by
where is the magnetic length. Then, close to the flat band’s energy (
), namely for
, the above formula is approximated to
The same behavior appears in the weakly singular cases () too. One of the simplest examples of the weakly singular flat band models is obtained by the band parameters
,
,
, and
, where
and
are free parameters with
. In this case, we have
. According to the level repulsion mechanism mentioned in the previous paragraph, it was shown that the approximate formula for the Landau levels corresponding to the flat band is given by
Figure 7. Lattice structures and band structures of the kagome and checkerboard lattices are shown in (a) and (e), respectively. Landau level structures of those models are plotted in (b) and (f). On the right-hand-side panels, the Landau levels around the flat band’s energy are highlighted. In (c) and (g), the dependency of the Landau level energies around the zero energy are exhibited. It is shown that the Landau levels obtained from the tight-binding lattice model (circular symbols) are consistent with the continuum results (solid curves) from (38) in (d) and (h)
![Figure 7. Lattice structures and band structures of the kagome and checkerboard lattices are shown in (a) and (e), respectively. Landau level structures of those models are plotted in (b) and (f). On the right-hand-side panels, the Landau levels around the flat band’s energy are highlighted. In (c) and (g), the 1/n dependency of the Landau level energies around the zero energy are exhibited. It is shown that the Landau levels obtained from the tight-binding lattice model (circular symbols) are consistent with the continuum results (solid curves) from (38) in (d) and (h)](/cms/asset/5895b2fe-87b0-4b16-b79e-2aeb0da45868/tapx_a_1901606_f0007_oc.jpg)
5.3 Lattice models
The results in Sec. 5.2 obtained from the continuum Hamiltonian (19) were successfully applied to several lattice models such as the kagome and the checkerboard lattice models [Citation82]. Such lattice models are described by the tight-binding Hamiltonians of the form
where is the hopping parameter between the
th and
th sites, and
represents the annihilation operator at the
th site. The effect of the magnetic field on this system can be considered by the Peierls substitution given by
where is the vector potential corresponding to the applied magnetic field, and
is the magnetic flux quantum [Citation158]. Denoting the magnetic flux penetrating the unit cell before applying the magnetic field by
, we only consider the cases, where
with natural numbers
and
.
The Landau levels of the kagome lattice with only the nearest-neighbor hopping processes are shown in ), and those of the checkerboard lattice model with the nearest and the next nearest hopping processes are plotted in ). Most importantly, the Landau level spreading phenomenon of the flat band discussed in the continuum model analysis appears in the Landau level structures of both lattice models too. According to the continuum result in (38), the Landau level spreading should be , where
is the Landau level spacing of the parabolic band, because the singular flat bands of the kagome and checkerboard lattices yield
around their band-crossing points. As shown in ) and , this continuum results (solid lines) are consistent with the numerical data obtained from the lattice models (circular symbols). Moreover, the
behavior of the Landau levels of the continuum model is clearly shown in the Landau levels of the lattice models too as plotted in ) and . However, this trend breaks down for large
’s due to the lattice effect. A distinguishing feature of the results of the lattice model is that there are finite number of Landau levels while there are infinite ones in the continuum model. In the lattice model, after the Peierls substitution, we have an enlarged unit cell, called the magnetic unit cell. Let us assume that the magnetic unit cell is
times larger than the original unit cell, where
is a natural number determined by the flux number
. Then, to preserve the total number of states before and after the application of the magnetic field, we have
magnetic bands, corresponding to the Landau levels, per each original energy band. In those two lattice models under the magnetic field described by
, one can find a gauge choice where the magnetic unit cell is
times larger than the original unit cell, and we have
magnetic bands corresponding to their flat band. It was shown that only
magnetic bands are developed in the band gap region below the flat band, and one remnant magnetic band is above it.
5.4 Realistic systems
In real materials, one cannot have a perfectly flat band because the hopping parameters of such systems cannot be fine-tuned as in the toy models treated in the previous subsections. For instance, while a perfectly flat band appears in the kagome lattice when we consider only the nearest neighbor hopping processes, it becomes dispersive as soon as we include long-range hopping processes or the spin–orbit coupling (SOC), which exist in the realistic situations. Let us denote the first and second nearest-neighbor hopping parameters by and
, and the SOC strength by
, respectively. Then, we have a nearly flat band when the later two processes are perturbatively small, namely
. There are several candidate materials with kagome-like structures hosting nearly flat bands belonging to different regimes in the parameter space of
and
as follows. First, the carbon networks such as the cyclic graphene, cyclic graphyne, and cyclic graphdiyne, belong to the case where
because the SOC can be neglected due to the low atomic number of carbon [Citation11,Citation18,Citation82]. In this case, the gap opening at the band-crossing point is almost negligible while the flat band obtains a small dispersion due to the finite value of
. Second, the nearly commensurate-charge density wave (NC-CDW) phase of 1 T-TaS
corresponds to the regime
because the second nearest-neighbor hopping processes as well as the longer range hopping parameters are exponentially suppressed due to the large insulating domains [Citation159]. Third, in the trans-Au-THTAP, both
and
are non-negligible although they are small enough to have a nearly flat band [Citation160]. In the last two cases, the singularity of the unperturbed flat band is completely removed by the energy gap induced by the SOC. However, this singularity is still manifested in the Landau level structure of the nearly flat band as shown in the following formula for the zeroth Landau level of the nearly flat band:
where , and
[Citation82]. Here, the zeroth Landau level
(
) corresponds to the lowest (highest) Landau level of the flat band touching with an electron-like (hole-like) parabolic band. The electron-like (hole-like) parabolic band can be obtained by the negative (positive) value of
. This formula can be applied to all the three cases discussed above, and it was shown that the maximum quantum distance can be extracted successfully from this in the strong magnetic field limit. More specifically, the Landau level spreading of the unperturbed flat band competes with the energy scale of the bandwidth of the nearly flat band. As a result, under a very weak magnetic field, the bandwidth of the nearly flat band dominates the Landau level spreading so that the Onsager’s semiclassical scheme can be successfully applied to obtain the Landau levels. However, if the magnetic field is strong enough so that the Landau level spreading is much larger than the bandwidth, the Landau level spreading of the nearly flat band can be clearly observed.
While we have dealt with electronic systems as realistic systems hosting the anomalous Landau level phenomenon so far, the artificial systems with bosonic states can also be another kind of promising candidates. Indeed, the Landau levels were observed in the acoustic crystal [Citation161], photonic crystal [Citation143], and photon resonator [Citation162], where the magnetic field of the electronic systems is artificially mimicked by the strain-induced gauge field. In this kind of artificial systems, the singular flat band can be realized more ideally than the electronic systems because of their wide controllability. For example, if the kagome lattice is realized in the photonic crystal, one can easily make smaller by increasing the lattice constant, and
can be also varied by controlling the coupling between the medium and light polarization [Citation163–166]. We expect that one might also study the quantum distance origin of the anomalous Landau levels of the singular flat band by applying the artificial magnetic field to the singular flat band in these artificial systems.
5.5 Diverging orbital magnetic susceptibility
At zero temperature, the orbital magnetic susceptibility is defined by
where is the total energy per unit area of the system under the magnetic field. As a representative example, the diamagnetic divergence of the orbital magnetic susceptibility of graphene can be understood from this definition as follows. The Landau levels of graphene are proportional to
. Since the Landau level degeneracy is proportional to
, the total energy below the Dirac point shows
dependence, which results in
divergence in the weak magnetic field limit. Therefore, the peculiar divergence of the orbital magnetic susceptibility of graphene is due to the
dependence of the Landau levels [Citation167].
On the other hand, the Landau levels of the singular flat band are proportional to as in the conventional cases. Instead, they show the intriguing
dependence [Citation82]. As a result, the total energy for a fully occupied singular flat band is given by
where we assume that the parabolic band touching with the singular flat band at the zero energy is electron-like, and the magnetic flux per unit cell is with
. This form of the total energy leads to the logarithmically divergent paramagnetism, namely
. On the contrary, if the parabolic band is hole-like, the singular flat band gives the diverging diamagnetism. This logarithmic divergence of the orbital magnetic susceptibility can be seen even for the partially filled flat band.
The fate of the orbital magnetic susceptibility when the flat band gains a small dispersion due to a finite SOC is also important [Citation82]. In this case, the singular flat band becomes an insulator due to the gap opening at the band-crossing point. When the resulting nearly flat band is fully occupied, namely the Fermi level is in the gap between the flat band and the parabolic band, it was found that the orbital magnetic susceptibility is no more divergent. However, the orbital magnetic susceptibility still remains finite even though there is no state in the gap, and its value could be large for the small enough SOC.
5.6 Landau level broadening due to disorder
When the system is disordered, its Landau levels cannot be described by the Dirac delta function type density of states anymore, and are broadened. It was shown that the Landau level broadening appears differently in the non-singular flat band and singular flat band, and the quantum distance plays a central role [Citation82]. Quantitatively, the Landau level broadening is defined as the energy width of the density of states of the broadened Landau level at its half-maximum. In the recent work [Citation82], an onsite disorder potential given by
is considered, where
is the site index,
is a random number ranging from
to
, and
is the disorder strength. First, in the case of the non-singular flat band, where all the Landau levels corresponding to the flat band are degenerate at the energy of the flat band, the Landau level broadening
is proportional to
, and almost independent of the magnetic field. This is a different behavior from the usual Landau level broadening of parabolic dispersions within the Born approximation given by
, where
is the relaxation time in the absence of the magnetic field. On the other hand, in the case of the singular flat band, all the Landau levels spread into the energy gap are non-degenerate, and the Landau level broadening of each of them should be investigated. In the weakly singular cases, it was found that the Landau level broadening of the singular flat band is due to the inter-band coupling between the singular flat band and the parabolic band touching with it. Indeed, the approximate formula for the Landau level broadening is given by
where is the cyclotron energy of the quadratic band. Note that
represents the strength of the inter-band coupling as discussed in Sec. 5.2. While the
behavior of the Born approximation is recovered, it is not the intrinsic property of the singular flat band. Instead, it arises from the Landau level broadening (
) of the parabolic band, which appears in EquationEquation.(48)
(48)
(48) through the inter-band coupling due to the singularity at the band-crossing point. To observe the Landau level spreading in experiments, it should be much larger than the Landau level broadening and it was shown in [Citation82] that the critical magnetic field and the disorder strength are within the experimentally accessible range.
6 Conclusions and outlook
We have reviewed the theoretical and experimental progress of the fundamental understanding of singular flat band systems. A singular flat band has a quite simple band structure, composed of a flat band crossing with another dispersive band at a momentum. However, the singularity of the flat band’s Bloch states gives unexpected physical outcomes. The presence of partially extended eigenstates, such as the non-contractible loop states and the robust boundary modes. More surprisingly, the singularity is deeply rooted from the quantum distance of Bloch wave functions around the band crossing point. Such intriguing geometric properties of the singular flat band eventually leads to the unusual Landau level spreading, which goes beyond the paradigm of Onsager’s semi-classical theory.
We believe that singular flat bands are an ideal new platform where intriguing geometric properties of Bloch states associated with quantum distance can be investigated. Although it is trivial in view of the conventional band topology, it provides an avenue to explore the geometry of quantum states from a different perspective, beyond the physics of Berry curvature. In this respect, recent experimental discovery of various materials hosting nearly flat bands, such as CoSn [Citation27,Citation28] and FeSn [Citation30], is quite encouraging. We note that while the most actively studied flat band system recently is the twisted bilayer graphene, the nearly flat bands of the twisted bilayer graphene just belong to the non-singular flat band according to the classification scheme in this review [Citation82]. On the other hand, the flat band of CoSn and FeSn, arising from the underlying kagome lattice structure, belongs to the category of singular flat bands. We propose that the flat band materials based on kagome lattice structure would provide new opportunities to reveal the geometric properties related to the nonzero quantum distance predicted in this paper. It is worth noting that, unlike the physics of Berry curvature and Berry phase, which is already well-established, relatively less attention has been paid to the physical significance of the quantum metric and quantum distance. We believe that the study of singular flat bands would pave the way for the complete understanding of geometric properties of Bloch wave functions in condensed matter physics, and it could lay the groundwork to figure out the underlying mechanism of the strongly correlated physics in flat band system, such as the superconductivity and magnetism.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- Sutherland B. Localization of electronic wave functions due to local topology. Phys Rev B. 1986;34:5208.
- Lieb EH. Two theorems on the hubbard model. Phys Rev Lett. 1989;62:1201.
- Arai M, Tokihiro T, Fujiwara T, et al. Strictly localized states on a two-dimensional penrose lattice. Phys Rev B. 1988;38:1621.
- Aoki H, Ando M, Matsumura H. Hofstadter butterflies for flat bands. Phys Rev B. 1996;54:R17296.
- Shima N, Aoki H. Electronic structure of super-honeycomb systems: a peculiar realization of semimetal/semiconductor classes and ferromagnetism. Phys Rev Lett. 1993;71:4389.
- Leykam D, Andreanov A, Flach S. Artificial flat band systems: from lattice models to experiments. Adv Phys X. 2018;3:1473052.
- Mielke A, Tasaki H. Tasaki,“ferromagnetism in the hubbard-model-examples from models with degenerate single-electron ground-states” commun. Comm Math Phys. 1993;158:341–41.
- Tasaki H. From nagaoka’s ferromagnetism to flat-band ferromagnetism and beyond: an introduction to ferromagnetism in the hubbard model. Prog Theoretical Phys. 1998;99:489–548.
- Mielke A. Stability of ferromagnetism in hubbard models with degenerate single-particle ground states. J Phys A. 1999;32:8411.
- Hase I, Yanagisawa T, Aiura Y, et al. Possibility of flat-band ferromagnetism in hole-doped pyrochlore oxides Sn2Nb2O7 and Sn2Ta2O7. Phys Rev Lett. 2018;120:196401.
- You J-Y, Gu B, Su G. flat band and hole-induced ferromagnetism in a novel carbon monolayer. Sci Rep. 2019;9:1–7.
- Volovik GE. The fermi condensate near the saddle point and in the vortex core. JETP Lett. 1994;59:830.
- Yudin D, Hirschmeier D, Hafermann H, et al. Fermi condensation near van hove singularities within the hubbard model on the triangular lattice. Phys Rev Lett. 2014;112:070403.
- Volovik GE. Graphite, graphene, and the flat band superconductivity. JETP Lett. 2018;107:516–517.
- Aoki H. Theoretical possibilities for flat band superconductivity. J Supercond Novel Magn2020;33:2341-2346.
- Kononov A, Endres M, Abulizi G, et al. Flat band superconductivity in type-ii weyl-semimetal wte2 induced by a normal metal contact. arXiv preprint arXiv:2007.04752, 2020.
- Congjun W, Bergman D, Balents L, et al. Flat bands and wigner crystallization in the honeycomb optical lattice. Phys Rev Lett. 2007;99:070401.
- Chen Y, Shenglong X, Xie Y, et al. Ferromagnetism and wigner crystallization in kagome graphene and related structures. Phys Rev B. 2018;98:035135.
- Jaworowski B, Güçlü AD, Kaczmarkiewicz P, et al. Wigner crystallization in topological flat bands. New J Phys. 2018;20:063023.
- Gersen H, Karle TJ, Engelen RJP, et al. Real-space observation of ultraslow light in photonic crystal waveguides. Phys Rev Lett. 2005;94:073903.
- Baba T, Mori D. Slow light engineering in photonic crystals. J Phys D Appl Phys. 2007;40:2659.
- Krauss TF. Slow light in photonic crystal waveguides. J Phys D Appl Phys. 2007;40:2666.
- Settle MD, Engelen RJP, Salib M, et al. Flatband slow light in photonic crystals featuring spatial pulse compression and terahertz bandwidth. Opt Express. 2007;15:219–226.
- Baba T. Slow light in photonic crystals. Nat Photonics. 2008;2:465–473.
- Corcoran B, Monat C, Grillet C, et al. Green light emission in silicon through slow-light enhanced third-harmonic generation in photonic-crystal waveguides. Nat Photonics. 2009;3:206–210.
- Schulz SA, O’Faolain L, Beggs DM, et al. Dispersion engineered slow light in photonic crystals: a comparison. J Opt. 2010;12:104004.
- Kang M, Fang S, Linda Y, et al. Topological flat bands in frustrated kagome lattice cosn. Nat Commun. 2020;11:4004.
- Liu Z, Man L, Wang Q, et al. Orbital-selective dirac fermions and extremely flat bands in frustrated kagome-lattice metal cosn. Nat Commun. 2020;11:1–7.
- Yin J-X, Zhang SS, Chang G, et al. Negative flat band magnetism in a spin–orbit-coupled correlated kagome magnet. Nat Phys. 2019;15:443–448.
- Kang M, Linda Y, Fang S, et al. Dirac fermions and flat bands in the ideal kagome metal fesn. Nat Mater. 2020;19:163–169.
- Linda Y, Kang M, Liu J, et al. Massive dirac fermions in a ferromagnetic kagome metal. Nature. 2018;555:638–642.
- Lin Z, Choi J-H, Zhang Q, et al. Flatbands and emergent ferromagnetic ordering in fe 3 sn 2 kagome lattices. Phys Rev Lett. 2018;121:096401.
- Yang TY, Wan Q, Wang YH, et al. Evidence of orbit-selective electronic kagome lattice with planar flat-band in correlated paramagnetic YCr6Ge6. arXiv preprint arXiv:1906.07140, 2019.
- Wang P, Wang Y, Zhang B, et al. Experimental observation of electronic structures of kagome metal YCr6Ge6. Chin Phys Lett. 2020;37:087102.
- Hase I, Yanagisawa T, Kawashima K. Flat-band in pyrochlore oxides: a first-principles study. Nanomaterials. 2019;9:876.
- Cao Y, Fatemi V, Fang S, et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature. 2018;556:43–50.
- Cao Y, Fatemi V, Demir A, et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature. 2018;556:80–84.
- Yoo H, Engelke R, Carr S, et al. Atomic and electronic reconstruction at the van der waals interface in twisted bilayer graphene. Nat Mater. 2019;18:448–453.
- Xie Y, Lian B, Jäck B, et al. Spectroscopic signatures of many-body correlations in magic-angle twisted bilayer graphene. Nature. 2019;572:101–105.
- Kerelsky A, McGilly LJ, Kennes DM, et al. Maximized electron interactions at the magic angle in twisted bilayer graphene. Nature. 2019;572:95–100.
- Iqbal Bakti Utama M, Koch RJ, Lee K, et al.Visualization of the flat electronic band in twisted bilayer graphene near the magic angle twist. Nat Phys. 2020;17:184-188.
- Ramires A, Lado JL. Electrically tunable gauge fields in tiny-angle twisted bilayer graphene. Phys Rev Lett. 2018;121:146801.
- Wang L, Shih E-M, Ghiotto A, et al.Correlated electronic phases in twisted bilayer transition metal dichalcogenides. Nature Mater Pages. 2020;19:861–866.
- Zhang Z, Wang Y, Watanabe K, et al.Flat bands in twisted bilayer transition metal dichalcogenides. Nat Phys. 2020;16:1093-1096.
- Zhi L, Zhuang J, Wang L, et al. Realization of flat band with possible nontrivial topology in electronic kagome lattice. Sci Adv. 2018;4:eaau4511.
- Kim K, DaSilva A, Huang S, et al. Tunable moiré bands and strong correlations in small-twist-angle bilayer graphene. Proceedings of the National Academy of Sciences, 114( 13):3364–3369, 2017.
- Liu -C-C, Zhang L-D, Chen W-Q, et al. Chiral spin density wave and d+ i d superconductivity in the magic-angle-twisted bilayer graphene. Phys Rev Lett. 2018;121:217001.
- Fengcheng W, MacDonald AH, Martin I. Theory of phonon-mediated superconductivity in twisted bilayer graphene. Phys Rev Lett. 2018;121:257001.
- Po HC, Zou L, Vishwanath A, et al. Origin of mott insulating behavior and superconductivity in twisted bilayer graphene. Phys Rev X. 2018;8:031089.
- Yankowitz M, Chen S, Polshyn H, et al. Tuning superconductivity in twisted bilayer graphene. Science. 2019;363:1059–1064.
- Lian B, Wang Z, Bernevig BA. Twisted bilayer graphene: a phonon-driven superconductor. Phys Rev Lett. 2019;122:257002.
- Roy B, Vladimir J. Unconventional superconductivity in nearly flat bands in twisted bilayer graphene. Phys Rev B. 2019;99:121407.
- Gonzalez J, Stauber T. Kohn-luttinger superconductivity in twisted bilayer graphene. Phys Rev Lett. 2019;122:026801.
- Nunes LHCM, Smith CM. Flat-band superconductivity for tight-binding electrons on a square-octagon lattice. Phys Rev B. 2020;101:224514.
- Balents L, Dean CR, Efetov DK, et al. Superconductivity and strong correlations in moiré flat bands. Nat Phys. 2020;16:725–733.
- Arora HS, Polski R, Zhang Y, et al. Superconductivity in metallic twisted bilayer graphene stabilized by WSe2. Nature. 2020;583:379–384.
- Peri V, Song Z, Andrei Bernevig B, et al. Fragile topology and flat-band superconductivity in the strong-coupling regime. arXiv preprint arXiv:2008. 02288, 2020.
- Guzmán-Silva D, Meja-Cortés C, Bandres MA, et al. Experimental observation of bulk and edge transport in photonic lieb lattices. New J Phys. 2014;16:063061.
- Vicencio RA, Cantillano C, Morales-Inostroza L, et al. Observation of localized states in lieb photonic lattices. Phys Rev Lett. 2015;114:245503.
- Mukherjee S, Spracklen A, Choudhury D, et al. Observation of a localized flat-band state in a photonic lieb lattice. Phys Rev Lett. 2015;114:245504.
- Xia S, Ramachandran A, Xia S, et al. Unconventional flatband line states in photonic lieb lattices. Phys Rev Lett. 2018;121:263902.
- Zong Y, Xia S, Tang L, et al. Observation of localized flat-band states in kagome photonic lattices. Opt Express. 2016;24:8877–8885.
- Xia S, Hu Y, Song D, et al. Demonstration of flat-band image transmission in optically induced lieb photonic lattices. Opt Lett. 2016;41:1435–1438.
- Jina M, Rhim J-W, Tang L, et al. Direct observation of flatband loop states arising from nontrivial real-space topology. Phys Rev Lett. 2020;124:183901.
- Mukherjee S, Thomson RR. Observation of localized flat-band modes in a quasi-one-dimensional photonic rhombic lattice. Opt Lett. 2015;40:5443–5446.
- Weimann S, Morales-Inostroza L, Real B, et al. Transport in sawtooth photonic lattices. Opt Lett. 2016;41:2414–2417.
- Travkin E, Diebel F, Denz C. Compact flat band states in optically induced flatland photonic lattices. Appl Phys Lett. 2017;111:011104.
- Cantillano C, Mukherjee S, Morales-Inostroza L, et al. Observation of localized ground and excited orbitals in graphene photonic ribbons. New J Phys. 2018;20:033028.
- Mukherjee S, Thomson RR. Observation of robust flat-band localization in driven photonic rhombic lattices. Opt Lett. 2017;42:2243–2246.
- Leykam D, Flach S. Perspective: photonic flatbands. Apl Photonics. 2018;3:070901.
- Milićević M, Montambaux G, Ozawa T, et al. Type-iii and tilted dirac cones emerging from flat bands in photonic orbital graphene. Phys Rev X. 2019;9:031010.
- Gyu-Boong J, Guzman J, Thomas CK, et al. Ultracold atoms in a tunable optical kagome lattice. Phys Rev Lett. 2012;108:045305.
- Taie S, Ozawa H, Ichinose T, et al. Coherent driving and freezing of bosonic matter wave in an optical lieb lattice. Sci Adv. 2015;1:e1500854.
- Drost R, Ojanen T, Harju A, et al. Topological states in engineered atomic lattices. Nat Phys. 2017;13:668–671.
- Masumoto N, Kim NY, Byrnes T; Naoyuki Masumoto, Na Young Kim, Tim Byrnes, Kenichiro Kusudo, Andreas Löffler, Sven Höfling, Alfred Forchel, and Yoshihisa Yamamoto. Exciton–polariton condensates with flat bands in a two-dimensional kagome lattice. New J Phys. 2012;14:065002.
- Nakata Y, Okada T, Nakanishi T, et al. Observation of flat band for terahertz spoof plasmons in a metallic kagomé lattice. Phys Rev B. 2012;85:205128.
- Kajiwara S, Urade Y, Nakata Y, et al. Observation of a nonradiative flat band for spoof surface plasmons in a metallic lieb lattice. Phys Rev B. 2016;93:075126.
- Slot MR, Gardenier TS, Jacobse PH. Guido CP van Miert, Sander N Kempkes, Stephan JM Zevenhuizen, Cristiane Morais Smith, Daniel Vanmaekelbergh, and Ingmar Swart. Experimental realization and characterization of an electronic lieb lattice. Nat Phys. 2017;13:672–676.
- Whittaker CE, Cancellieri E, Walker PM, et al. Exciton polaritons in a two-dimensional lieb lattice with spin-orbit coupling. Phys Rev Lett. 2018;120:097401.
- Chen L, Mazaheri T, Seidel A, et al. The impossibility of exactly flat non-trivial chern bands in strictly local periodic tight binding models. J Phys A. 2014;47:152001.
- Rhim J-W, Yang B-J. Classification of flat bands according to the band-crossing singularity of bloch wave functions. Phys Rev B. 2019;99:045107.
- Rhim J-W, Kim K, Yang B-J. Quantum distance and anomalous landau levels of flat bands. Nature. 2020;584:59–63.
- Hatsugai Y. Bulk-edge correspondence in graphene with/without magnetic field: chiral symmetry, dirac fermions and edge states. Solid State Commun. 2009;149:1061–1067.
- Essin AM, Gurarie V. Bulk-boundary correspondence of topological insulators from their respective green’s functions. Phys Rev B. 2011;84:125132.
- Fukui T, Shiozaki K, Fujiwara T, et al. Bulk-edge correspondence for chern topological phases: a viewpoint from a generalized index theorem. J Phys Soc Jpn. 2012;81:114602.
- Graf GM, Porta M. Bulk-edge correspondence for two-dimensional topological insulators. Commun Math Phys. 2013;324:851–895.
- Bourne C, Carey AL, Rennie A. The bulk-edge correspondence for the quantum hall effect in kasparov theory. Lett Math Phys. 2015;105:1253–1273.
- Rhim J-W, Behrends J, Bardarson JH. Bulk-boundary correspondence from the intercellular zak phase. Phys Rev B. 2017;95:035421.
- Rhim J-W, Bardarson JH, Slager R-J. Unified bulk-boundary correspondence for band insulators. Phys Rev B. 2018;97:115143.
- Hatsugai Y. Chern number and edge states in the integer quantum hall effect. Phys Rev Lett. 1993;71:3697.
- Balents L, Fisher MPA. Chiral surface states in the bulk quantum hall effect. Phys Rev Lett. 1996;76:2782.
- Schulz-Baldes H, Kellendonk J, Richter T. Simultaneous quantization of edge and bulk hall conductivity. J Phys A. 2000;33:L27.
- Kellendonk J, Richter T, Schulz-Baldes H. Edge current channels and chern numbers in the integer quantum hall effect. Rev Math Phys. 2002;14:87–119.
- Yang K-Y, Yuan-Ming L, Ran Y. Quantum hall effects in a weyl semimetal: possible application in pyrochlore iridates. Phys Rev B. 2011;84:075129.
- Kane CL, Mele EJ. Quantum spin hall effect in graphene. Phys Rev Lett. 2005;95:226801.
- Kane CL, Mele EJ. Z 2 topological order and the quantum spin hall effect. Phys Rev Lett. 2005;95:146802.
- Andrei Bernevig B, Hughes TL, Zhang S-C. Quantum spin hall effect and topological phase transition in hgte quantum wells. science. 2006;314:1757–1761.
- Liang F, Kane CL, Mele EJ. Topological insulators in three dimensions. Phys Rev Lett. 2007;98:106803.
- Liang F, Kane CL. Topological insulators with inversion symmetry. Phys Rev B. 2007;76:045302.
- Bužek V, Hillery M. Quantum copying: beyond the no-cloning theorem. Phys Rev A. 1996;54:1844.
- Berry MV. Geometric phases in physics. Singapore: World Scientific; 1989.
- Dodonov VV, Man’Ko OV, Man’Ko VI, et al. Hilbert-schmidt distance and non-classicality of states in quantum optics. J Mod Opt. 2000;47:633–654.
- Provost JP, Vallee G. Riemannian structure on manifolds of quantum states. Commun Math Phys. 1980;76:289–301.
- Zanardi P, Giorda P, Cozzini M. Information-theoretic differential geometry of quantum phase transitions. Phys Rev Lett. 2007;99:100603.
- Yu-Quan M, Chen S, Fan H, et al. Abelian and non-abelian quantum geometric tensor. Phys Rev B. 2010;81:245129.
- Srivastava A, Ataç I. Signatures of bloch-band geometry on excitons: nonhydrogenic spectra in transition-metal dichalcogenides. Phys Rev Lett. 2015;115:166802.
- Fubini G. Sulle metriche definite da una forma hermitiana. Opere Selecte. 1904;2:10–20.
- Study E. Kürzeste wege im komplexen gebiet. Mathematische Annalen. 1905;60:321–378.
- Xiao D, Chang M-C, Niu Q. Berry phase effects on electronic properties. Rev Mod Phys. 2010;82:1959.
- Zahid Hasan M, Kane CL. Colloquium: topological insulators. Rev Mod Phys. 2010;82:3045.
- Neupert T, Chamon C, Mudry C. Measuring the quantum geometry of bloch bands with current noise. Phys Rev B. 2013;87:245103.
- Raoux A, Piéchon F, Fuchs J-N, et al. Orbital magnetism in coupled-bands models. Phys Rev B. 2015;91:085120.
- Piéchon F, Raoux A, Fuchs J-N, et al. Geometric orbital susceptibility: quantum metric without berry curvature. Phys Rev B. 2016;94:134423.
- Ozawa T. Steady-state hall response and quantum geometry of driven-dissipative lattices. Phys Rev B. 2018;97:041108.
- Peotta S, Päivi T. Superfluidity in topologically nontrivial flat bands. Nat Commun. 2015;6:1–9.
- Julku A, Peotta S, Vanhala TI, et al. Geometric origin of superfluidity in the lieb-lattice flat band. Phys Rev Lett. 2016;117:045303.
- Gianfrate A, Bleu O, Dominici L, et al. Measurement of the quantum geometric tensor and of the anomalous hall drift. Nature. 2020;578:381–385.
- Min Y, Yang P, Gong M, et al. Experimental measurement of the quantum geometric tensor using coupled qubits in diamond. Natl Sci Rev. 2020;7:254–260.
- Hwang Y, Jung J, Rhim J-W, et al. Wave function geometry of band crossing points in two-dimensions. arXiv preprint arXiv:2012.15122, 2020.
- Bergman DL, Congjun W, Balents L. Band touching from real-space topology in frustrated hopping models. Phys Rev B. 2008;78:125104.
- Maimaiti W, Andreanov A, Park HC, et al. Compact localized states and flat-band generators in one dimension. Phys Rev B. 2017;95:115135.
- Dubail J, Read N. Tensor network trial states for chiral topological phases in two dimensions and a no-go theorem in any dimension. Phys Rev B. 2015;92:205307.
- Read N. Compactly supported wannier functions and algebraic k-theory. Phys Rev B. 2017;95:115309.
- Gligorić G, Maluckov A, Hadžievski L, et al. Nonlinear localized flat-band modes with spin-orbit coupling. Phys Rev B. 2016;94:144302.
- Zegadlo K, Dror N, Hung NV, et al. Single and double linear and nonlinear flatband chains: spectra and modes. Phys Rev E. 2017;96:012204.
- Raoux A, Morigi M, Fuchs J-N, et al. From dia-to paramagnetic orbital susceptibility of massless fermions. Phys Rev Lett. 2014;112:026402.
- Piéchon F, Fuchs JN, Raoux A, et al. Tunable orbital susceptibility in α-t 3 tight-binding models. In Journal of Physics: Conference Series, volume 603, 012001. IOP Publishing, 2015.
- Vanderbilt D, King-Smith RD; David Vanderbilt and RD King-Smith. Electric polarization as a bulk quantity and its relation to surface charge. Phys Rev B. 1993;48:4442.
- Ryu S, Hatsugai Y. Entanglement entropy and the berry phase in the solid state. Phys Rev B. 2006;73:245115.
- Delplace P, Ullmo D, Montambaux G. Zak phase and the existence of edge states in graphene. Phys Rev B. 2011;84:195452.
- Damascelli A, Hussain Z, Shen Z-X. Angle-resolved photoemission studies of the cuprate superconductors. Rev Mod Phys. 2003;75:473.
- Chiu C-K, Teo JCY, Schnyder AP, et al. Classification of topological quantum matter with symmetries. Rev Mod Phys. 2016;88:035005.
- Efremidis NK, Sears S, Christodoulides DN, et al. Discrete solitons in photorefractive optically induced photonic lattices. Phys Rev E. 2002;66:046602.
- Fleischer JW, Segev M, Efremidis NK, et al. Observation of two-dimensional discrete solitons in optically induced nonlinear photonic lattices. Nature. 2003;422:147–150.
- Longhi S, Marangoni M, Lobino M, et al. Observation of dynamic localization in periodically curved waveguide arrays. Phys Rev Lett. 2006;96:243901.
- Schwartz T, Bartal G, Fishman S, et al. Transport and anderson localization in disordered two-dimensional photonic lattices. Nature. 2007;446:52–55.
- Peleg O, Bartal G, Freedman B, et al. Conical diffraction and gap solitons in honeycomb photonic lattices. Phys Rev Lett. 2007;98:103901.
- Lahini Y, Avidan A, Pozzi F, et al. Anderson localization and nonlinearity in one-dimensional disordered photonic lattices. Phys Rev Lett. 2008;100:013906.
- Lederer F, Stegeman GI, Christodoulides DN, et al. Discrete solitons in optics. Phys Rep. 2008;463:1–126.
- Chen Z, Segev M, Christodoulides DN. Optical spatial solitons: historical overview and recent advances. Rep Prog Phys. 2012;75:086401.
- Nixon M, Ronen E, Friesem AA, et al. Observing geometric frustration with thousands of coupled lasers. Phys Rev Lett. 2013;110:184102.
- Segev M, Silberberg Y, Christodoulides DN. Anderson localization of light. Nat Photonics. 2013;7:197–204.
- Rechtsman MC, Zeuner JM, Tünnermann A, et al. Strain-induced pseudomagnetic field and photonic landau levels in dielectric structures. Nat Photonics. 2013;7:153–158.
- Rechtsman MC, Plotnik Y, Zeuner JM, et al. Topological creation and destruction of edge states in photonic graphene. Phys Rev Lett. 2013;111:103901.
- Rechtsman MC, Zeuner JM, Plotnik Y, et al. Photonic floquet topological insulators. Nature. 2013;496:196–200.
- Plotnik Y, Rechtsman MC, Song D, et al. Observation of unconventional edge states in ‘photonic graphene’. Nat Mater. 2014;13:57–62.
- Zeuner JM, Rechtsman MC, Plotnik Y, et al. Observation of a topological transition in the bulk of a non-hermitian system. Phys Rev Lett. 2015;115:040402.
- Hassan AE, Kunst FK, Moritz A, et al. Corner states of light in photonic waveguides. Nat Photonics. 2019;13:697–700.
- Miura Davis K, Miura K, Sugimoto N, et al. Writing waveguides in glass with a femtosecond laser. Opt Lett. 1996;21:1729–1731.
- Pertsch T, Peschel U, Lederer F, et al. Discrete diffraction in two-dimensional arrays of coupled waveguides in silica. Opt Lett. 2004;29:468–470.
- Kohmoto M. Topological invariant and the quantization of the hall conductance. Ann Phys. 1985;160:343–354.
- Mikitik GP, Sharlai YV. Manifestation of berry’s phase in metal physics. Phys Rev Lett. 1999;82:2147.
- Kopelevich Y, Kopelevich Y; Igor A Luk?yanchuk. Phase analysis of quantum oscillations in graphite. Phys Rev Lett. 2004;93:166402.
- Onsager L. Interpretation of the de haas-van alphen effect. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 43( 344):1006–1008, 1952.
- Roth LM. Semiclassical theory of magnetic energy levels and magnetic susceptibility of bloch electrons. Phys Rev. 1966;145:434.
- Gao Y, Niu Q. Zero-field magnetic response functions in landau levels. Proceedings of the National Academy of Sciences, 114( 28):7295–7300, 2017.
- Fuchs J-N, Piéchon F, Montambaux G. Landau levels, response functions and magnetic oscillations from a generalized onsager relation. SciPost Phys. 2018;4:024.
- Hofstadter DR. Energy levels and wave functions of bloch electrons in rational and irrational magnetic fields. Phys Rev B. 1976;14:2239.
- Lee JM, Geng C, Park JW, et al. Stable flatbands, topology, and superconductivity of magic honeycomb networks. Phys Rev Lett. 2020;124:137002.
- Yamada MG, Soejima T, Tsuji N, et al. First-principles design of a half-filled flat band of the kagome lattice in two-dimensional metal-organic frameworks. Phys Rev B. 2016;94:081102.
- Wen X, Qiu C, Yajuan Q, et al. Acoustic landau quantization and quantum-hall-like edge states. Nat Phys. 2019;15:352–356.
- Schine N, Ryou A, Gromov A, et al. Synthetic landau levels for photons. Nature. 2016;534:671–675.
- Bliokh KY, Niv A, Kleiner V, et al. Geometrodynamics of spinning light. Nat Photonics. 2008;2:748.
- Nalitov AV, Malpuech G, Terças H, et al. Spin-orbit coupling and the optical spin hall effect in photonic graphene. Phys Rev Lett. 2015;114:026803.
- Sotto M, Debnath K, Tomita I, et al. Spin-orbit coupling of light in photonic crystal waveguides. Phys Rev A. 2019;99:053845.
- Zhang Z, Liang S, Feng L, et al. Spin-orbit coupling in photonic graphene. Optica. 2020;7:455–462.
- Ominato Y, Koshino M. Orbital magnetism of graphene nanostructures. Solid State Commun. 2013;175:51–61.