Singular flat bands

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 nonsingular 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 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.


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 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.

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 [1][2][3][4][5][6]. In condensed matter physics, flat bands are considered ideal to study many-body phenomena such as the ferromagnetism [7][8][9][10][11], superconductivity [12][13][14][15][16], and Wigner crystal formation [17][18][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 [20][21][22][23][24][25][26].
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 view point of the topological band theory, the flat band is simply considered trivial due to 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 [78,79]. 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 2nπ (trivial) due to the band-flatness condition. [80]. Therefore, the flat band cannot exhibit the conventional bulk-boundary correspondence [81][82][83][84][85][86][87] such as the chiral edge modes from the nonzero Chern number in Chern bands [88][89][90][91], the Fermi arcs between two topological monopole charges with opposite signs in the Weyl semimetal [92], and helical surface modes from Z 2 bulk topological invariant of topological insulators [93][94][95][96][97].
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 ψ 1 and ψ 2 are normalized quantum states [80,[98][99][100]. The quantum distance d, which takes a value between 0 and 1, measures how close the two states ψ 1 and ψ 2 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 k 1 and k 2 is usually becomes zero when |k 1 − k 2 | → 0 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. [79,80]. 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 bandcrossing point is a bulk invariant measuring the strength of the singularity at the band crossing point. [80]. 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 [63,79]. Moreover, the Landau levels of the singular flat band with a quadratic band-crossing show anomlous behavior, which is characterized by the maximum quantum distance [80]. 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 (or distance). In solids, Berry curvature and the quantum metric take the roles of the curvature and metric in the geometric description of the Bloch wave functions respectively, in momentum space. While the curvature part has been studied exhaustively in past decades related to the semiclassical formulation of electron dynamics [101] and the topological classifications of materials [102], physical implications of the metric or distance part have been much less investigated although there have been several studies on the physical roles of them in the current noise [103], the orbital magnetic susceptibility [104], the superfluid weight of a nearly flat Chern band [105], etc. To achieve a unified understanding of the geometric nature of solids, further studies on the relations between the quantum metric or distance and physical quantities in various condensed matter phenomena have to be perfomed. 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 Sec. 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 Sec. 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 Sec. 4, we show how to define the strength of the singularity of the singular flat band via the maximum quantum distance or the psedospin structure. In Sec. 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 the divergence of the orbital magnetic susceptibility.
Unit cell (c) Then, we show several candidate realistic materials, where those results can be observed. Finally, in Sec. 6, we give conclusions and outlook of the research on the singular flat band systems.
2 Flat band and singuarity

Compact localized state
The group velocity, the first derivative of an energy dispersion with respect to crystal momentum, represents the mobiility 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 [1,4,79,[106][107][108][109]. As a representative example, the compact localized state of the kagome lattice with the nearest neighbor hopping processes is shown in Fig. 1(a). 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. 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. 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. 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 [79]. If the hopping processes are allowed within a finite range, each matrix element of the Bloch Hamiltonian H(k) of the tight binding model of such system is in the form of the finite sum of Bloch phases, which is given by where d is the dimension of the system described by d primitive vectors ā i , k i = k · ā i , n i is an integer ranging from p lm i to q lm i , and h lm (n 1 , · · · , n d ) is a complex coefficient of the Bloch phase e in1k1 · · · e in d k d . In a simple manner, we denote the Bloch phas as e ik·R = e in1k1 · · · e in d k d , where R = (n 1 , · · · , n d ) is the lattice vector. Then, the eigenvector v fb for the flat band with energy ǫ fb is obtained by solving coupled homogeneous equations given by [H(k)−ǫ fb ]v fb (k) = 0. Since all the matrix elements of H(k) − ǫ fb are also in the form of the finite sum of Bloch phases, one can always find a unnormalized solution for the v fb (k) 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 |ψ fb (k) is the Bloch eigenstate of the flat band, c χ is the normalization constant for |χ R . Here, α k is a mixing coefficient of the Bloch eigenstates with different k's, in the form of a multiplication between |v fb (k)| and an arbitrary finite sum of Bloch phases. Note that the Bloch eigenstate is in the form |ψ fb (k) = N −1/2 R ′ Q j=1 e ik·R ′ṽ fb (k)| j |R ′ , j , where N is the total number of the unit cells in the system,ṽ fb (k)| j is the j-th component of the normalized eigenvector of the flat band given byṽ fb = v fb /|v fb (k)|, |R ′ , j represents the j-th orbital in the unit cell located at r = R ′ , and Q is the size of the Bloch Hamiltonian. As a result, the amplitude of the j-th orbital of |χ R in the unit cell at R ′ is obtained as Since α k v fb /|v fb (k)| is in the form of the finite sum of Bloch phases, the value of R ′ , j|χ R vanishes if (R − R ′ ) is out of the range of the lattice vectors in the Bloch phases in α k v fb /|v fb (k)|. This implies that |χ R is a compact localized state, and one can have various kinds of the compact localized state depending on the choice of α k . 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 ā 1 = a(1, 0), ā 2 = a(−1/2, √ 3/2), and ā 3 = −ā 1 −ā 2 . If the unnormalized eigenvector for the flat band of this model is chosen to be v kagome fb = (e ik1 −1, 1− e −ik2 , e −ik2 −e ik1 ) T , a compact localized state centered at r = 0 is obatined from (4) as This is the compact localized state with the smallest size in the kagome lattice model as plotted in Fig. 1(a). The j-th element of the column vector on the right-hand side of (6) represents the amplitude of the compact localized state at the site-j, which is indicated by three colors in Fig. 1(a). For example, the first element −δ R ,(−1,0) + δ R ,(0,0) of this vector means that the amplitudes at the site-1 in the unit cells located at R = (−1, 0) and R = (0, 0) are the coefficients of δ R ,(−1,0) and δ R ,(0,0) , namely -1 and 1 respectively. In fact, these coeffecients can be inferred from the unnormalized form of the Bloch eigenstate v kagome fb because there is an one-to-one correspondence between the coefficients of e −ik·R ′ in v kagome fb and that of δ R ,R ′ in (6). In general, one can easily derive the compact localized state from the unnormalized Bloch eigenstate in this way.

Singular and non-singular flat bands
One fundamental question on the compact localized state is that can one find a complete set of N compact localized states to span the whole flat band? [79,106] We need to answer this question to understand the localized nature of the flat band completely. It seems quite obvious that N compact localized states centered at N different unit cells are linearly independent. Here N 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 N compact localized states fail to form a complete set sometimes [106]. 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. [79] This singular point can only be generated by the bandcrossing 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. [78]. A flat band with the singularity is called a singular flat band, and otherwise, it is called a non-singular flat band [79].
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 (k = 0). We check whether the corresponding Bloch eigenstate is singular or non-singular at the band-crossing point as follows. The normalized where |v kagome This implies thatṽ kagome fb 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 Fig. 1(b). In this model, only the nearest neighbor hopping processes with amplitude 1/2 are allowed between the sites in the same layer, and the inter-layer hopping parameters corresponding to the red and green dashed lines are −2 and 1/2 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ṽ BSL fb = (1, 1) T / √ 2 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 [79].
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 [79]. To show this, we denote N number of different compact localized states (labeled by j) by which have different shapes and center positions depending on the choice of α (j) k and R j respectively. Then, a determinant diagnosing their linear independence or dependence is given by because the Bloch wave functions |ψ fb (k) 's with different momenta in (9) form a complete set. As noted in Sec. 2.1, to ensure the compact localization of |χ R j , α (j) k should be proportional to |v fb (k)|, where v fb (k) 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 v fb (k) is zero at a crystal momentum, namely the Bloch eigenstatẽ v fb (k) is singular at that momentum, D also becomes zero implying that N compact localized states are always linearly dependent of each other. As one can note from (9), in this case, the Bloch wave function corresponding to the singular momentum is always absent in constructing the compact localized state because α (j) k is zero at that momentum. Namely, N compact localized states of a singular flat band are made of a smaller (less than N ) number of Bloch wave functions, which implies that they cannot be linearly independent. For example, in the kagome lattice, the N compact localized states obtained by translating the minimal compact localized state in (6) 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 v fb (k) for the flat band, one can obtain non-zero D by choosing α (j) k = |v fb (k)| and {R j } to be N distinct lattice vectors. Therefore, one can find a complete set of N 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 ā i direction by N i , so that the total number of the unit cells in the system is given by N = d i=1 N i . Then, the crystal momentum discretized under periodic boundary condition along each reciprocal lattice vector is k i = (2π/N i )m, where m is an integer. Therefore, even though the continuum expression of the band dispersion shows a singular band-crossing at a certain crystal momentum k 0 , one can find a complete set of N compact localized states if k 0 cannot be represented by the form k i = (2π/N i )m for the given system size N i . 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 k 0 = (π, π). Note that this k 0 is allowed only when both N x and N y are even numbers. If we choose a compact localized state for this model as shown in Fig. 1(c), any linear combination of the N = N x N y translated copies of this compact localized state is nonzero for odd values of N x and N y . This implies that these compact localized states form a complete set. On the other hand, if N x and N y are even integers, the same staggered sum of them vanishes. In this case, these compact localized states are linearly dependent on each other. Note that if the singular band-crossing is at k = 0, we are free from this finite size issue in general.

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 noncompact 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 Fig. 1(a) [106]. They are extended at least along one spatial direction while compactly localized along other directions. 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. 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. In the kagome lattice case, it was noted that N − 1 number of compact localized states can be linearly independent of each other [106]. Therefore, the existence of these two additional non-contractible loops states is consistent with the fact that there are N + 1 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. 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. 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.
In general, it was shown that such non-contractible loop states are guaranteed to exist in a two dimensional singular flat band [79]. If the band crossing point is at k = k 0 , 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 the equation (9). As a result, any kind of eigenmodes including the Bloch wave function at k = k 0 is linearly independent of the compact localized states. While the Bloch wave function at k = k 0 is extended along both directions ā 1 and ā 2 , 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 H(k 1 , k 2 ) be the Hamiltonian under consideration. If we fix k 1 to k 0,1 , where k 0 = (k 0,1 , k 0,2 ), we obtain an effective one dimensional flat band Hamiltonian H(k 0,1 , k 2 ) with momentum k 2 . Since any one dimensional flat bad is a non-singular flat band as noted previously, one can always find a compact localized state for this effective Hamiltonian localized along ā 2 direction, which includes the Bloch wave function at k = k 0 too. This is the non-contractible loop state of the full Hamiltonian H(k 1 , k 2 ) because it is still extended along ā 1 direction with momentum k 0,1 . In the similar way, one can find another non-contractible loop state extended along ā 2 direction too.
In three dimensional case, one can have different kind of partially extended states with nontrivial real space topology, called the non-contractible planar states (NPSs) [79]. If we denote a singular flat band Hamiltonian for this case by H(k 1 , k 2 , k 3 ), one should fix two components of k = (k 1 , k 2 , k 3 ) to the corresponding components of the singular momentum k 0 = (k 0,1 , k 0,2 , k 0,3 ) to obtain an effective one dimensional flat band Hamiltonian in this case. Let us assume that those two fixed momenta are k 1 = k 0,1 and k 2 = k 0,2 , for convenience. Then, we can construct a planar eigenstate, namely a NPS, which is extended along two directions ā 1 and ā 2 , and compactly localized along ā 3 . By fixing other pairs of momentum components, one can obtain other NPSs along different directions.

Open boundary system
Robust boundary mode (a) (b) 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 R 1 to R 9 ). If one subtract a compact localized state centered at R 2 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 [79]

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 [79]. The robust boundary mode of the kagome lattice is shown in Fig. 2(a) 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 Fig. 2(c) 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 [88][89][90][91], Z 2 index in topological insulators [93][94][95][96][97], the Zak phase of the one dimensional insulators with reflection symmetry [86,[110][111][112], 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. Therefore, one cannot detect the robust boundary mode by electronic spectroscopy such as the angle-resolved photoemission spectroscopy (ARPES) [113]. 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. 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. It is worth noting that the compact localized states of a singluar flat bands are linearly dependent so that their linear combination becomes zero under periodic boundary condition, whereas the same linear combination generates a robust boundary mode under open boundary condition. 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 few number of compact localized This figure is adapted from [79] 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 Fig. 2(b). Let us note that the presence of robust boundary modes is confirmed in a recent experiment in photonic lattices as discussed in Sec. 3.

Flatness enforced band-crossing
Conventionally, a band-crossing is protected by crystalline symmetries or topological charge [114]. 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 [79]. 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 σ x , σ y , σ z are Pauli matrices, and σ 0 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 H k yields a flat band and a parabolic band touching each other at k = 0. One can note that the flat band is a singular flat band only when both t 1 and t 4 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 H k can be always transformed toH k = (t ′ 1 k 2 x + t ′ 2 k x k y + t ′ 3 k 2 y )(σ z + σ 0 ). SinceH ′ k has no off-diagonal terms, a perturbation of the form H ′ = δλ(σ z − σ 0 ) 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 δλ = 0. This is summarized in Fig. 3 (b).
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 H ′ = m x σ x + m y σ y + m z σ z + m 0 σ 0 is added to the system. From the flatness conditions, detH k = 0 and det H k + H ′ = 0, we always obtain t 4 = 0 if one of the mass parameters is nonzero. This contradicts with the fact that t 4 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 Fig. 3(c). 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 H ′ = δ(λ 2 + λ 7 ), where λ i 's are the Gell-Mann matrices, to H kagome (k), 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 [79].
3 Observation of the real-space topology in photonic lattices

Flat bands in photonic lattices
One can fabricate a two dimensional periodic array of optical waveguides, socalled the photonic lattice [115][116][117][118][119][120][121][122][123][124][125][126][127][128][129][130], in a three dimensional medium by using the femtosecond laser writing technique [131,132]. A light beam propagating in this system satisfies a wave equation given by within the paraxial approximaition [116,121,122], where ψ(x, y, z) describes the electric field envelope of the light propagating along z-axis, n 0 and n(x, y) are the refractive indices outside and inside the waveguide of the medium, k 0 is the wavenumber in vacuum, and ∇ 2 ⊥ = ∂ 2 x + ∂ 2 y . This can be regarded as a 2D Schrödinger equation by substituting z with time, and n(x, y) with the periodic potential. Namely, the 2D cross section of the medium perpendicular to the z-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 [57][58][59][60] and kagome lattice models [61][62][63] have been studied by using the photonic lattice experimental set-ups [64][65][66][67][68][69].
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 z-direction from z = 0 to z = L z , an incident light beam, whose amplitudes and phases correspond to those of the compact localized state, enters into the waveguides at z = 0. As the light beam propagates along z-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  [79] compact localized state remains 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 (z = L z ), is the same incident light beam only multiplied by a phase factor, and this is how one demonstrate the existence of a flat band in the photonic lattice experiments.

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 [60]. 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 simple open boundary condition, and we need a more delicate decoration of the boundaries for the line state to become an eigenstate [63]. 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 Fig. 4(a) was realized in the photonic lattice as shown in Fig. 4(b). Then, as in Fig. 4(c) and (d), 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 Fig. 4(e), 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 Fig. 4(f). 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.

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 [133]. 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 [92]. 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 C is a closed path enclosing the band-crossing point and |v k is the cell-periodic part of the Bloch wave function [134,135]. Here, the singularity of the Bloch wave function at the touching point also plays a key role when gamma is nonzero modulo 2π 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 [80]. First, the generic singular flat band with a quadratic band crossing does not have a quantized Berry phase (22) because all the Pauli matrices are used in the generic Hamiltonian in (14). Even though the Berry phase is quantized due to some symmetries such as C 2z T and P T , where C 2z , P , and T are two-fold rotation about the z-axis, inversion, time-reversal symmetries, respectively, the Berry phase of the singular flat band reads γ = 2πn with integer n. 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 [80]. From (1), the quantum distance between two Bloch states at momenta k 1 and k 2 is given by d 2 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 k 1 and k 2 are momenta close to the band-crossing point. By definition, the strength of the singularity d max is valued from 0 to 1. In the case of the non-singular flat band, the Bloch eigenmode of (14) is always independent of momentum, so that the quantum distance is zero between an arbitrary pair of momenta. As a result, d max = 0 for the non-singular flat band. On the other hand, in the case of the singular flat band, d max remains finite even if we consider momenta extremely close to the band-crossing point. Note that the eigenvector of the Hamiltonian (14) is independent of k, which is the magnitude of the momentum displacement with respect to the band crossing point, because all the dependencies on k can be factored out from this quadratic Hamiltonian. The general form of the maximum quantum distance of the Hamiltonian (14) is evaluated as where m 1 and m 2 are the minimum and maximum effective masses of the quadratic band touching with the flat band. The strength of the singularity, d max , also represents how strong the interband coupling between the flat band and parabolic band as explained below by a counting argument [80]. First, the number of independent parameters of the Hamiltonian (14) is four because five parameters among nine parameters in (14) are determined by others through the flatness conditions from (16) to (20). Since we need three parameters to describe the dispersion relation of the quadratic band of the form E quad (k) = a 1 k 2 x +a 2 k x k y +a 3 k 2 y , 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 (14) can be represented by four parameters M xx , M xy , M yy , and is the mass tensor for the parabolic band. This implies that d max 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 [79]. The eigenvector of the flat band of this model is evaluated as v fb = 1 where θ is the polar angle with respect to the band-crossing point at k = 0. Then, the quantum distance between two eigenvectors at θ 1 and θ 2 is obtained as This gives us the maximum quantum distance d max = 1 when θ 1 − θ 2 = π/2, which is consistent with the formula (24) because the band parameters for H eff kagome are given by t 1 = b 1 = −t 3 = b 3 = 1/2, t 4 = −1, 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.

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 |v k is the eigenvector of the flat band and f α (k) is the coefficient of σ α in (14) given by f x (k) = t 6 k 2 y , f y (k) = t 4 k x k y + t 5 k 2 y , and f z (k) = t 1 k 2 x + t 2 k x k y + t 3 k 2 y [80]. One can represent the pseudospin on a surface of a unit sphere because f x (k), f y (k), and f z (k) are all in quadratic order of k = |k| 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 Fig. 5. Since the eigenvector of (14) is independent of k, 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 Fig. 5(a). 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 Fig. 5(b) and (c). 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 k 1 and k 2 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 ∆θ max = max∆θ(k 1 , k 2 ) as the strength of the singularity of the singular flat band instead of the maximum quantum distance d max . 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 (∆θ max = π or d max = 1), 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 HamiltonianH eff Then,H eff kagome is described by f x (k) = (−k 2 x + k 2 y )/2, f y (k) = −k x k y , f z (k) = 0, and f 0 (k) = (k 2 x + k 2 y )/2. As a result, the pseudospin formula is given by which shows double winding around the band-crossing point, and ∆ max = π.
5 Anomalous Landau levels of singular flat bands

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, S 0 (ε) is the area of a closed path at the energy ε on a band, B is the magnetic field, e is the electric charge, h = 2π is the Planck constant, n is an integer corresponding to the Landau level index, and γ ε,B is the quantum correction due to Berry phase, magnetic susceptibility, and higher order magnetic responses [134,[136][137][138][139]. Only the energy values of ε satisfying the above equation are allowed as eigenenergies, and they are quantized because the Landau level index n 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 v F , we have S 0 (ε) = π(ε/ v F ) 2 and γ ε,B = π. This leads to the well-known Landau level quantization ε = v F (2e Bn) 1/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 S 0 (ε) 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 S 0 (ε) is not singlevalued. What we can naively infer from the fact that the effective mass m * is infinite in the flat band is that the cyclotron energy ( ω c = eB/m * ) 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 [80].
To consider the effect of the magnetic field on the flat band, the transformations k x → (a + a † )/ √ 2l B and k y → i(a − a † )/ √ 2l B are applied to the quadratic band-crossing flat band Hamiltonian in (14), where a and a † are the annihilation and creation operators, respectively satisfying [a, a † ] = 1. To ensure the Hermiticity, k x k y is symmetrized by replacing it with (k x k y + k y k x )/2. Then, each element of (14) becomes an operator. This Hamiltonian is solved by applying a trial wave function of the form where |u n is the normalized eigenstate of the harmonic oscillator satisfying a|u n = √ n|u n−1 and a † |u n = √ n + 1|u n+1 , and C n and D n are complex coefficients.
The generic structure of the Landau levels of the flat band model (14) obtained by the procedure described above is shown in Fig. 6. 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 equation (31) 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 [80]. In the case shown in Fig. 6(a), 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 d max of the Bloch eigenstate before applying magnetic field. More specifically, as shown in Fig. 6(e), the ratio ∆/ ω c is a simple monotonic function of d max given by where ∆ is the Landau level spreading and ω c = eB/ √ m 1 m 2 is the cyclotron frequency of the parabolic band. It is important to note that ∆/ ω c is only a function of d max regardless of an arbitrary choice of band parameters of the model Hamiltonian in (14). 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 d max . 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 d max share the same value of ∆/ ω c . 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, was shown to be proportional to d max in the weakly singular limit (d max ≪ 1).
While the explicit form of ξ(x) is unavailable, its approximate form, which fits the numerical data accurately, was found to be where c = 0.041, r = 2.344, and s = 0.487. This was obtained from several analytic results such as ξ(d 2 max → 0) = (3/40)d 2 max and ξ(d 2 max → 1) = (2 √ 3 − 3)/4. By using this formula, one can extract the geometric information, namely the maximum quantum distance d max , 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 n close to the flat band's energy as shown in Fig. 6(d). Here, the energy of the flat band is assume to be zero without loss of generality. As an example, in the case of the maximally singular (d max = 1) model in (25), the explicit formula for the Landau levels corresponding to the singular flat band is given by where l B = /eB is the magnetic length. Then, close to the flat band's energy (E = 0), namely for n ≫ 1, the above formula is approximated to The same behavior appears in the weakly singular cases (d max ≪ 1) too. One of the simplest example of the weakly singular flat band models is obtained by the band parameters t 1 = b 1 = t 3 , t 2 = b 2 = t 5 = 0, t 6 = t 4 4t 2 1 + t 2 4 /(2t 1 ), and b 3 = t 2 1 + t 2 4 (4t 2 1 + t 2 4 )/(4t 2 1 ), where t 1 and t 4 are free parameters with t 4 ≪ t 1 . In this case, we have d max ≈ t 4 /(2t 1 ). 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

Lattice models
The results in Sec. 5.2 obtained from the continuum Hamiltonian (14) were successfully applied to several lattice models such as the kagome and the checkerboard lattice models [80]. Such lattice models are described by the tight binding Hamiltonians of the form  where t ij is the hopping parameter between the i-th and j-th sites, and c i represents the annihilation operator at the i-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 φ 0 = h/e is the magnetic flux quantum [140]. Denoting the magnetic flux penetrating the unit cell before applying the magnetic field by φ, we only consider the cases, where φ/φ 0 = p/q with natural numbers p and q.
The Landau levels of the kagome lattice with only the nearest neighbor hopping processes are shown in Fig. 7(b), and those of the checkerboard lattice model with the nearest and the next nearest hopping processes are plotted in Fig. 7(f). 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 (33), the Landau level spreading should be ∆ = ω c (2 where ω c is the Landau level spacing of the parabolic band, because the singular flat bands of the kagome and checkerboard lattices yield d max = 1 around their band-crossing points. As shown in Fig. 7(d) and (h), this continuum results (solid lines) are consistent with the numerical data obtained from the lattice models (circular symbols). Moreover, the 1/n behavior of the Landau levels of the continuum model is clearly shown in the Landau levels of the lattice models too as plotted in Fig. 7(c) and (g). However, this trend breaks down for large n'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 Q times larger than the original unit cell, where Q is a natural number determined by the flux number φ/φ 0 . Then, to preserve the total number of states before and after the application of the magnetic field, we have Q magnetic bands, corresponding to the Landau levels, per each original energy band. In those two lattice models under the magnetic field described by φ/φ 0 = 1/q, one can find a gauge choice where the magnetic unit cell is q times larger than the original unit cell, and we have q magnetic bands corresponding to their flat band. It was shown that only q − 1 magnetic bands are developed in the band gap region below the flat band, and one remnant magnetic band is above it.

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 t 1 and t 2 , and the SOC strength by λ respectively. Then, we have a nearly flat band when the later two processes are perturbatively small, namely t 2 , λ ≪ t 1 . There are several candidate materials with kagome-like structures hosting nearly flat bands belonging to different regimes in the parameter space of t 2 /t 1 and λ/t 1 as follows. First, the carbon networks such as the cyclic graphene, cyclic graphyne, and cyclic graphdiyne, belong to the case where λ ≈ 0 because the SOC can be neglected due to the low atomic number of carbon [11,18,80]. In this case, the gap opening at the band-corssing point is almost negligible while the flat band obtains a small dispersion due to the finite value of t 2 /t 1 . Second, the nearly commensurate-charge density wave (NC-CDW) phase of 1T-TaS 2 corresponds to the regime t 2 /t 1 ≈ 0 because the second nearest neighbor hopping processes as well as the longer range hopping parameters are exponentially suppressed due to the large insulating domains [141]. Third, in the trans-Au-THTAP, both t 2 /t 1 and λ are non-negligible although they are small enough to have a nearly flat band [142]. 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 a 0 = − √ 3π(t 1 /2 + 2t 2 ), and a 1 = π 2 (t 2 1 + 8t 2 2 ) [80]. Here, the zeroth Landau level E + 0 (E − 0 ) corresponds to the lowest (highest) Landau level of the flat band touching with an electron-like (hole-like) parabolic band. The electronlike (hole-like) parabolic band can be obtained by the negative (positive) value of t 1 . This formula can be applied to all the three cases discussed in the 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 band width of the nearly flat band. As a result, under very weak magnetic field, the band width 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 band width, the Landau level spreading of the nearly flat band can be clearly observed.

Diverging orbital magnetic susceptibility
At zero temperature, the orbital magnetic susceptibility is defined by where E tot 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 √ B. Since the Landau level degeneracy is proportional to B, the total energy below the Dirac point shows B √ B dependence, which results in 1/ √ B divergence in the weak magnetic field limit. Therefore, the peculiar divergence of the orbital magnetic susceptibility of graphene is due to the √ B dependence of the Landau levels [143].
On the other hand, the Landau levels of the singular flat band are proportional to B as in the conventional cases. Instead, they show the intriguing 1/n dependence [80]. 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 φ/φ 0 = 1/q with q ≫ 1. This form of the total energy leads to the logarithmically divergent paramagnetism, namely χ orb ∝ − ln B + const. 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 [80]. In this case, the singular flat band becomes an insulator due to the gap opening at the bandcrossing 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.

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 nonsingular flat band and singular flat band, and the quantum distance plays a central role [80]. 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 [80], an onsite disorder potential given by V disorder = ρ i ǫ i c † i c i is considered, where i is the site index, ǫ i is a random number ranging from −0.5 to 0.5, 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 Γ = ( 2 ω c /2πτ 0 ) 1/2 ∝ √ B, where τ 0 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 ω c is the cyclotron energy of the quadratic band. Note that d max represents the strength of the inter-band coupling as discussed in Sec. 5.2. While the √ B 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 (∼ √ B) of the parabolic band, which appears in Eq.(43) 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 [80] that the critical magnetic field and the disorder strength are within the experimentally accessible range.

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 metric or 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 [27,28] and FeSn [30], 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 [80]. 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 groundwork to figure out the underlying mechanism of the strongly correlated physics in flat band system, such as the superconducivity and magnetism.

Disclosure statement
No potential conflict of interest was reported by the authors.