## Interaction-driven plateau transition between integer and fractional Chern Insulators

Leon Schoonderwoerd, 1 Frank Pollmann, 2, 3 and Gunnar M¨ oller 1

1 Physics and Astronomy, Division of Natural Sciences,

University of Kent, Canterbury CT2 7NZ, United Kingdom

2 Department of Physics, TFK, Technische Universit¨t M¨nchen, a u

James-Franck-Straße 1, D-85748 Garching, Germany

3 Munich Center for Quantum Science and Technology (MCQST), Schellingstraße 4, D-80799 M¨nchen, Germany u

We present numerical evidence of an interaction-driven quantum Hall plateau transition between a | C &gt; | 1 Chern Insulator (CI) and a ν = 1 3 Laughlin state in the Harper-Hofstadter model. / We study the model at flux densities p/q , where the lowest Landau level (LLL) manifold comprises p magnetic sub-bands. For weak interactions, the model realises integer CIs corresponding to filled sub-bands, while strongly interacting candidate states include fractional quantum Hall (FQH) states at LLL filling fractions ν = r/t . These phases may compete at the same particle density when p = . As a concrete example, we numerically explore the physics at flux density t n φ = 3 11, where / we show evidence that a direct transition occurs between a CI and a ν = 1 3 Laughlin state, which / we characterise in terms of its critical, topological and entanglement properties. We also show that strong interactions generically stabilise a ν = 1 3 Laughlin state even when the LLL is split / into multiple bands, and introduce a powerful methodology to extract its topological entanglement entropy by exploiting the scaling of magnetic length with n φ .

Lattice fractional quantum Hall (FQH) effects have received intense recent interest, owing to proposals and subsequent successful generation of topological band structures in cold atomic gases [1-8], and using superlattice structures in graphene systems [9-13]. Similar prospects exist in 2D materials with strong spin-orbit coupling [14-16], although no concrete materials have yet been identified. All of these realizations can provide topological flat bands with an effective magnetic length comparable to the lattice spacing. When combined with strong repulsive interactions, such systems may exhibit fractional Chern insulator (FCI) states: FQH liquids stabilised specifically by a lattice potential [10, 17-19]. While literature on FCIs focuses on models where a single low-lying flat topological band mirrors the physics of the lowest Landau level (LLL) [14-16, 20-22], composite fermion (CF) theory approaches lattice quantum Hall (QH) effects from a more general angle, and predicts incompressible states even when the LLL is split into numerous sub-bands [17, 18, 23, 24].

scription [38]. The critical point is then described by a massless multi-flavour Dirac theory coupled to a gauge field [37, 38, 40-42]. A microscopic realization was found in graphene heterostructures, where transitions between FCI states within the same Jain series are driven by tuning the strength of a periodic potential [40].

QH states are inherently interesting for fundamental physics as they can support exotic fractionalized quasiparticles obeying non-Abelian exchange statistics [25] and could be exploited as a platform for universal topological quantum computation [26]. Transitions between QH states have likewise attracted much attention [2730]. Conventionally, QH plateau transitions are driven by field strength tuning the Landau level filling. Their critical properties can be understood in terms of the percolation of current-carrying modes [29], captured by network models [31] or renormalisation group approaches [32, 33]. Several scenarios for transitions into or between FQH states have also been described in clean systems [34-40]. These are often cast as Chern number changing transitions in the underlying CF [40] or parton de-

This manuscript explores interaction-driven QH plateau transitions in the Harper-Hofstadter model [4345], motivated partly by its recent experimental realization [2-4, 9]. We study flux densities n φ where the LLL is split into several magnetic sub-bands to allow competition between integer sub-band filling and fractional filling of the LLL. We show that repulsive interaction strengths can drive transitions between integer and fractional CIs in the Harper-Hofstadter model. Specifically, we provide numerical evidence suggesting a direct transition between a C = +4 Chern insulator (CI) and a ν = 1 3 Laugh-/ lin state [46], at n φ = 3 11, where the LLL comprises / three sub-bands. These two states are members of different Jain series, as the flux attached to CFs changes. Hence, this interaction-driven plateau transition cannot be described as CF Chern number changing, and the previously considered critical theories [38, 40] do not apply.

Most studies of lattice FQH effects or FCIs have studied regimes where the low-energy manifold of states consists of a single band [19-22, 47-53], simplifying the numerical treatment. Here, we explicitly consider situations where the LLL is split into multiple sub-bands. This situation arises naturally in the Harper-Hofstadter model (Fig. 1), owing to its self-similar bandstructure [45, 54]. The model is described by the fermionic tight-binding Hamiltonian

$$\mathfrak { s } _ { \mathfrak { s } _ { \mathfrak { g } _ { - } } } \hat { \mathcal { H } } _ { H o t } = - t \sum _ { \langle i j \rangle } \left [ e ^ { i \phi _ { i j } } \hat { a } _ { i } ^ { \dagger } \hat { a } _ { j } + \mathrm H. c. \right ] + V \sum _ { \langle i j \rangle } \hat { n } _ { i } \hat { n } _ { j }, \quad ( 1 )$$

with hopping amplitude t , nearest-neighbour repulsion

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-e3c2eae5-67cc-43f8-b9f2-5ef7a8c22fe2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-e1576318-8154-460a-8caf-e312477a47e3.jpg)

-

Figure 1. (a) The Harper-Hofstadter model in the Landau gauge with constant flux per plaquette (center) is 'wrapped' around the cylinder (bottom), and corresponding Matrix Product State representation (top). (b) Single-particle band structure at n φ = 3 11 (lowest 3 bands). /

V , annihilation (creation) operators ˆ (ˆ ), the number a a † operator ˆ = ˆ ˆ and n a a † φ ij = ∫ i j A · d l are Aharonov-Bohm phases on 〈 i, j 〉 . We choose the Landau gauge with vector potential A = xB ˆ y to yield a homogeneous magnetic field strength B with flux density n φ = B/ Φ 0 in units of the magnetic flux quantum Φ 0 = h e .

At rational n φ = p/q , the single-particle Hofstadter spectrum consists of q bands, each with areal density of states n band = 1 /q . The largest gaps arise at particle density n = n φ , and are connected continuously to LLL physics in the limit of n φ → 0 [45, 55]. As a corollary, for general n φ this LLL is split into p bands. All other gaps in the spectrum can be considered as fractal copies of the principal gap in subordinate instances of the self-similar spectrum [45]. Moreover, all bands are topological, carrying a non-zero Chern number which can be inferred from the diophantine relations, n = Cn φ + s [56, 57].

Non-interacting fermions filling the lowest m bands in the spectrum yield a CI state with integer Hall conductivity σ H = e /h 2 ∑ m i =1 C i [56-59]. However, with increasing V &gt; 0, interaction energy starts to dominate over kinetic energy, and can no longer be treated perturbatively on the non-interacting ground-state. Consequently, the system favours larger inter-particle separation at the expense of lifting particles beyond the lowestlying bands. The question now is which magnetic bands contribute significantly to the ground state interacting problem. CF theory of the Harper-Hofstadter model posits a regime dominated by a low-energy manifold of m s = qn s bands located below a sufficiently large singleparticle gap [17-19], enabling the prediction of candidate FCI states with filling fraction in the low-energy manifold ν = n/n s = r/ kCr ( +1), with k, r ∈ Z [19]. For flux densities n φ allowing m s = | kCr +1 bands in the low-| energy manifold, the CI filling the lowest m bands can compete directly with an FCI realised at the same particle density, but involving states in the lowest m s bands. In principle, larger interaction strength could lead to further mixing of higher-lying bands, destabilising the CF FQH state.

Concretely, we consider the simplest example of such a transition between a CI and a ν = 1 3 Laughlin state /

in the C = 1 LLL manifold ( k = 2, r = 1). We require m s = 2 | × 1 × 1 + 1 | = 3 sub-bands in the LLL, thus n φ = 3 /q , or p = m s = 3. We further specialise to an example with n φ = 3 11, where the lowest sub-band has / a Chern number C 1 = 4.

We numerically analyse the competition between the candidate states using infinite Density Matrix Renormalization Group (iDMRG) methods [60, 61], expressing (1) in a Matrix Product Operator (MPO) representation [Fig. 1(b)] [62, 63] for the well-established infinitecylinder lattice geometry [64-68]. Entanglement grows linearly with cylinder circumference, limiting the maximum circumference L y that can be accurately approximated at a given Matrix Product State (MPS) bond dimension χ .

Fig. 2 summarizes the evolution of several observables as a function of interaction strength V , for our target scenario: n φ = 3 11, / L y = 6. For additional data for L y = 7 and L y = 8, see Appendix D. We find multiple indications of a direct CI-FCI phase transition at a critical V c .

The first sign of a phase transition is found in the correlation length ξ . At V c , our data for ξ χ ( ) suggest a divergence with bond dimension [Fig. 2(a)]. This is consistent with a continuous quantum phase transition occurring at V c . The entanglement spectrum (ES) {-log λ i ( k y ) } of the Schmidt decomposition into two half-infinite states | Ψ = 〉 ∑ i,k y λ i ( k y ) | χ L i,k y 〉 ⊗ | χ R i,k y 〉 is shown in Fig. 2(b). At V c , it displays a sudden opening of the lowest gap, as well as other discontinuous changes in higher entanglement energies. In the ES, we define the entanglement gap ∆ as the gap between the lowest two levels within the ξ same quantum number and momentum sector. Fig. 2(c), shows a discontinuity of ∆ ξ at V c , and a tendency of the gap to decrease for V &lt; V c approaching the transition. Thus, the entanglement properties of the system also support a phase transition.

Fig. 2(d) displays the single-particle density matrix ρ (1) ij = 〈 ˆ ˆ c c † j i 〉 calculated within our MPS unit cell of 1 × 6 magnetic unit cells, yielding six k -points in k y -space. In the CI phase, its spectrum is dominated by six low-lying eigenvalues and a seventh nearby (including at least one state at each k -point), while in the FCI phase a gap opens above 18 low-lying states and two further nearby states (including at least three states at each k -point). This matches the expectation that only states in the lowest band contribute significantly to ρ (1) ij in the CI phase, while it involves all states in the LLL (i.e., three bands) for the Laughlin state [69]. As a function of V , ρ (1) ij changes discontinuously between the CI and Laughlin phases at V c , signalling a direct transition. Fig. 2(e) shows the fidelity susceptibility χ F = ( | 〈 ψ V ( ) | ψ V ( + δV ) 〉 | 2 -1) / δV ( ) 2 , which peaks at V c , also supporting a direct phase transition.

Finally, we measure the Hall conductance using flux

λ

log

-

V

Figure 2. Several observables through the transition (dashed line). All results are for L y = 6, n φ = 3 11, / χ = 500 [except 600 ≤ χ ≤ 1000 in (a)]. (a) Correlation length for several values of the bond dimension. (b) Momentum-space entanglement spectrum (colors denote different quantum number sectors). Also shows the entanglement gap. (c) Entanglement gap. (d) Spectrum of the single-particle density matrix. The dashed line indicates the lowest 6 (18) states on the left (right) of the transition. (e) Fidelity susceptibility χ F = ( |〈 ψ V ( ) | ψ V ( + δV ) 〉| 2 -1) / δV ( ) 2 . (f) Charge 〈 ∆ q 〉 transported along the cylinder by the insertion of 3 × 2 π flux.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-d442c7d3-381f-4d77-8071-3e9d5798a079.jpg)

insertion [57, 70]. We calculate the ground state's evolution under an adiabatic insertion of 3 × 2 π external flux through the cylinder, and measure the charge ∆ q thereby transported along the cylinder [66, 68, 71, 72], which yields the Hall conductivity via σ H = ∆ q/ ∆ φ ext . Flux insertion for V = 0 confirms that the system realizes an integer CI state with Hall conductivity proportional to the Chern number of the lowest Harper-Hofstadter band, σ H ∝ C 1 = 4. Fig. 2(f) shows that this picture persists for small V &gt; 0. At interaction strengths V &gt; V c , we find σ H = 1 3, consistent with the prediction for a / ν = 1 3 / Laughlin state. The Hall conductivity σ H ( V ) discontin- uously jumps between these two extremal values. However, near V c we cannot reliably perform the flux pumping procedure. Here, the insertion violates adiabaticity, evidenced by discontinuities in q φ ( ). This is understood by the increase of the correlation length near V c , requiring larger bond dimensions to fully capture the many-body state.

In sum, these multiple signatures clearly demonstrate a direct transition from the σ H = 4 CI state to a σ H = 1 3 / Laughlin state. Near the transition, we see that our iDMRG calculation is not fully converged in the bond dimension χ , highlighted by the continued increase in correlation length with χ in Fig. 2(a). This behaviour would be consistent with a continuous transition. Contrarily, the state's entanglement structure appears to evolve discontinuously at the transition, suggesting a (weakly) first order transition. As iDMRG on cylinders is incapable to unambiguously distinguish the nature of a quantum phase transition in two dimensions, our data could reflect either scenario. Further indications could be found by examining the evolution of the ground state degeneracy in a torus geometry.

While the transition at L y = 6 is sharply defined even at low χ , features at higherL y are more challenging to resolve using DMRG. At L y = 7 a sharp feature in ξ and S emerges only between 2000 ≤ χ ≤ 3000 while we only resolved a broader maximum for L y = 8 up to χ = 2000. Additionally, the ES and spectrum of ρ ij vary more smoothly rather than discontinuously with V when examined at low χ = 600, as shown in Appendix D. However, all three circumferences clearly support distinct phases at low and high V .

To identify possible theories describing the system at criticality, we extract the central charge for a putative conformal field theory from the simultaneous scaling of S E and ξ with χ . These are expected to follow the relation S E ( χ ) ∼ ( c/ 6) log ξ χ ( ) [73], where c is the central charge, under finite entanglement scaling. This procedure was originally developed for 1D systems [74, 75], but has subsequently been applied successfully to 2D systems [76, 77]. We show numerical scaling data for V ∼ V c and L y = 6 7 8 in Fig. 3(a), alongside the linear scaling pre-, , dicted for c = 1 2 and 3 (additionally, see Appendix D , for underlying data).

The excellent fits show a rapid increase of c with L y , which suggests that the critical system may be described by a Fermi Surface (FS) theory. Following [76], and assuming a circular FS with radius k F glyph[lscript] B = √ 2 / 3, we count the expected numbers of 'wires' N w cutting the FS (due to the quantisation of k y = 1 L y (2 πn + φ CF ), with φ CF the emergent flux insertion experienced by CF) for different values of L y (Fig. 3(c) illustrates L y = 8).

As shown in Fig. 3(b), our results conform with an ordinary Fermi liquid (where c = N w ) iff the CF experience periodic (anti-periodic) boundary conditions at L y = 6, 8 ( L y = 7), respectively.

Figure 3. (a) Scaling of S E and ξ with 1000 ≤ χ ≤ 2000 (3000) for L y = 6: blue circles, L y = 7: orange squares (and stars) and L y = 8: green triangles, at V = 1 634 . , 1 5 . , (and 1 55), and 1 45 respectively, and expected linear scaling . . for c = 1, 2, (2) and 3, respectively. Dashed black lines denote linear fits. (b) c L ( y ), as inferred from linear fit to data (blue circles) and predicted for CF Fermi surface with φ CF = 0 (orange crosses) and φ CF = π (green pluses). Error bars denote one standard deviation σ . (c) View of a putative Fermi surface for L y = 8, showing intersections with three (two) lines of allowed k y momentum for φ CF = 0( π ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-bda52921-4326-4382-8b31-c0f68738a1e1.jpg)

We further substantiate the identification and stability of the Laughlin state, when the LLL-manifold has several magnetic sub-bands, i.e., n φ = p/q , p &gt; 1. Firstly, we compute the entanglement entropy S E of bisected cylinders with varying L y . Unlike previous (i)DMRG simulations of the Harper-Hofstadter model, we simultaneously vary n φ , motivated by the computational advantages this brings in the torus geometry [48, 49]. This allows us to combine data for several n φ , providing a significant technical improvement upon earlier works. We obtain a much larger data set at reasonably low computational cost. Importantly, we find that the entanglement entropy depends on L y expressed in units of the magnetic length glyph[lscript] B = √ glyph[planckover2pi1] /eB = a √ q/ 2 πp , so smaller values of n φ allow for larger L y in units of lattice sites. Our data for S E ( L /glyph[lscript] y B ) are shown in Fig. 4, revealing a collapse onto a single line [78]. This data is consistent with the area law prediction of the entanglement entropy, S E = αL y + S topo , with S topo the topological entanglement entropy [79]. Extrapolation from our data yields S topo = -0 551(15), accurately matching the prediction . of S topo = -ln √ 3 ≈ -0 549 for the . ν = 1 3 Laughlin / state [80-82]. The collapsed data include cases where the LLL consists of p = 2 3 , , or 4 bands, confirming that all of these cases consistently yield the topological nature of the Laughlin state, which is further supported by the respective entanglement spectra shown in Appendix C.

Discussion: Given our data supporting a direct QH plateau transition between a CI state and FCI state in the Hofstadter model, we examine possible scenarios explaining the transition. We consider field theory descriptions for a continuous transition. Several works have considered CF flux-attachment transitions between different QHstates [34-40]. Transitions between a superconductor and a ν = 1 2 Laughlin state [38] and between different / states in the same Jain series [40, 42] seem most similar

Figure 4. Entanglement entropy in the Laughlin phase at V = 10 as a function of cylinder circumference L y (in units of the magnetic length glyph[lscript] B ) at χ = 600. Results for cylinders of various 4 ≤ L y ≤ 10 and flux densities n φ = p/q are combined.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-d5736ee2-f241-4a28-8fbf-2100a695fed5.jpg)

to our current setting. Both of these describe the transition as a Chern number changing transition experienced by a single fermionic parton band and the gap closing at multiple Dirac points simultaneously.

In our work, the CI phase does not require a parton construction, and is described by fermions in the lowest Hofstadter band with Chern number C 1 = 4. The FCI phase is described by two parton fields experiencing C = 1 Chern bands to effectuate flux attachment, while the third parton field describes the effective physics of the CFs-also a C = 1 Chern band for the Laughlin state. Hence, to describe the transition, one would have to explain the simultaneous change of the CF band Chern number, and fractionalisation of particles into three partons to give rise to the two additional fields. While effective theories in which this type of fractionalisation must occur at once may exist, a more generic case might display a sequence of several transitions affecting individual parton fields. Furthermore, our numerics indicate that the central charge at the critical point grows systematically with cylinder circumference (cf. Fig. 3). While a theory of multiple Dirac points would predict central charges similar in magnitude to ours, it cannot explain the scaling of c L ( y ). A single transition would thus require a more exotic explanation, such as the emergence of a 2D Fermi surface at the critical point [83]. The construction of theories for the interaction-driven CI-FCI plateau transition forms an exciting field of future study.

In conclusion, we have shown evidence of a direct phase transition between a ν = 1 3 Laughlin state and an in-/ teger CI in the Harper-Hofstadter model at flux density n φ = 3 11. / Additionally, we used simultaneous scaling of cylinder circumference and flux density to obtain smooth data coverage of the cylinder measured in magnetic lengths and obtained an improved estimate for the topological entanglement entropy.

While we have provided a first step towards a theoretical explanation for the CI-FCI phase transition in terms of a 2D Fermi surface, it would be interesting to find a more complete theoretical description. Moreover, while the transition at n φ = 3 11 appears to be a direct /

phase transition, we have found some indications that at n φ = 3 10, the transition occurs via a yet uncharacter-/ ized intermediate phase. Future research could elucidate this latter transition by formulating a field theoretical model of flux attachment in the transition.

The authors would like to thank Nigel R. Cooper for helpful discussions, Steven H. Simon for comments on the manuscript, and Johannes Hauschild for support using TenPy [63]. LS acknowledges funding from a University of Kent Vice Chancellor's Scholarship. FP acknowledges the support of the DFG Research Unit FOR 1807 through grants no. PO 1370/2-1, TRR80, the Nanosystems Initiative Munich (NIM) by the German Excellence Initiative, the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy - EXC-2111-390814868 and the European Research Council (ERC) under the European Unions Horizon 2020 research and innovation program (grant agreement no. 771537). GM acknowledges funding from the Royal Society under grants URF UF120157 and URF R 180004 and from the Engineer-\ \ ing and Physical Sciences Research Council (EPSRC) under grant EP/P022995/1. GM and FP acknowledge the hospitality of the Aspen Center for Physics, which is supported by National Science Foundation grant PHY1607611. GM thanks the Theory of Condensed Matter Physics Group at the Cavendish Laboratory for their hospitality.

- [1] Y.-J. Lin, R. L. Compton, K. Jim´ enez-Garc´ ıa, J. V. Porto, and I. B. Spielman, 'Synthetic magnetic fields for ultracold neutral atoms,' Nature 462 , 628-632 (2009).
- [2] MAidelsburger, M Atala, M Lohse, J T Barreiro, B Paredes, and I Bloch, 'Realization of the Hofstadter Hamiltonian with Ultracold Atoms in Optical Lattices,' Physical Review Letters 111 , 185301 (2013).
- [3] Hirokazu Miyake, Georgios A. Siviloglou, Colin J. Kennedy, William Cody Burton, and Wolfgang Ketterle, 'Realizing the Harper Hamiltonian with Laser-Assisted Tunneling in Optical Lattices,' Physical Review Letters 111 , 185302 (2013).
- [4] M Aidelsburger, M Lohse, C Schweizer, M Atala, J T Barreiro, S Nascimb` ene, N R Cooper, I Bloch, and N Goldman, 'Measuring the Chern number of Hofstadter bands with ultracold bosonic atoms,' Nature Physics 11 , 162-166 (2014).
- [5] Gregor Jotzu, Michael Messer, R´mi Desbuquois, Martin e Lebrat, Thomas Uehlinger, Daniel Greif, and Tilman Esslinger, 'Experimental realization of the topological Haldane model with ultracold fermions,' Nature 515 , 237-240 (2014).
- [6] N.R. Cooper, 'Rapidly rotating atomic gases,' Adv. Phys. 57 , 539-616 (2008).
- [7] Jean Dalibard, Fabrice Gerbier, Gediminas Juzeli¯nas, u and Patrik Ohberg, ' ¨ Colloquium : Artificial gauge potentials for neutral atoms,' Reviews of Modern Physics

## 83 , 1523-1543 (2011).

- [8] N Goldman, G Juzeli¯nas, P Ohberg, u ¨ and I B Spielman, 'Light-induced gauge fields for ultracold atoms,' Reports on Progress in Physics 77 , 126401 (2014).
- [9] C. R. Dean, L. Wang, P. Maher, C. Forsythe, F. Ghahari, Y. Gao, J. Katoch, M. Ishigami, P. Moon, M. Koshino, T. Taniguchi, K. Watanabe, K. L. Shepard, J. Hone, and P. Kim, 'Hofstadter's butterfly and the fractal quantum Hall effect in moir´ superlattices,' Nature e 497 , 598-602 (2013).
- [10] Eric M. Spanton, Alexander A. Zibrov, Haoxin Zhou, Takashi Taniguchi, Kenji Watanabe, Michael P. Zaletel, and Andrea F. Young, 'Observation of fractional Chern insulators in a van der Waals heterostructure,' Science 360 , 62-66 (2018).
- [11] T. Hensgens, U. Mukhopadhyay, P. Barthelemy, R. F. L. Vermeulen, R. N. Schouten, S. Fallahi, G. C. Gardner, C. Reichl, W. Wegscheider, M. J. Manfra, and L. M. K. Vandersypen, 'A capacitance spectroscopy-based platform for realizing gate-defined electronic lattices,' J. Appl. Phys. 124 , 124305 (2018).
- [12] Carlos Forsythe, Xiaodong Zhou, Kenji Watanabe, Takashi Taniguchi, Abhay Pasupathy, Pilkyung Moon, Mikito Koshino, Philip Kim, and Cory R. Dean, 'Band structure engineering of 2D materials using patterned dielectric superlattices,' Nat. Nanotechnol. 13 , 566-571 (2018).
- [13] Christina Knapp, Eric M. Spanton, Andrea F. Young, Chetan Nayak, and Michael P. Zaletel, 'Fractional Chern insulator edges and layer-resolved lattice contacts,' Phys. Rev. B 99 , 081114(R) (2019).
- [14] Evelyn Tang, Jia-Wei Mei, and Xiao-Gang Wen, 'HighTemperature Fractional Quantum Hall States,' Physical Review Letters 106 , 236802 (2011).
- [15] Titus Neupert, Luiz Santos, Claudio Chamon, and Christopher Mudry, 'Fractional quantum hall states at zero magnetic field,' Physical Review Letters 106 , 236804 (2011).
- [16] Kai Sun, Zhengcheng Gu, Hosho Katsura, and S Das Sarma, 'Nearly flatbands with nontrivial topology,' Physical Review Letters 106 , 236803 (2011).
- [17] A. Kol and N. Read, 'Fractional quantum Hall effect in a periodic potential,' Physical Review B 48 , 8890-8898 (1993).
- [18] Gunnar M¨ oller and N. R. Cooper, 'Composite fermion theory for bosonic quantum hall states on lattices,' Physical Review Letters 103 , 105303 (2009).
- [19] Gunnar M¨ oller and N. R. Cooper, 'Fractional Chern Insulators in Harper-Hofstadter Bands with Higher Chern Number,' Physical Review Letters 115 , 126401 (2015).
- [20] N. Regnault and B. A. Bernevig, 'Fractional Chern Insulator,' Physical Review X 1 , 021014 (2011).
- [21] Emil J. Bergholtz and Zhao Liu, 'Topological Flat Band Models and Fractional Chern Insulators,' International Journal of Modern Physics B 27 , 1330017 (2013).
- [22] Siddharth A Parameswaran, Rahul Roy, and Shivaji L Sondhi, 'Fractional quantum Hall physics in topological flat bands,' Comptes Rendus Physique 14 , 816-839 (2013).
- [23] Jainendra K Jain, 'Composite-fermion approach for the fractional quantum Hall effect,' Physical Review Letters 63 , 199-202 (1989).
- [24] Ana Lopez and Eduardo Fradkin, 'Fractional quantum

Hall effect and Chern-Simons gauge theories,' Phys. Rev. B 44 , 5246-5262 (1991).

- [25] Gregory Moore and Nicholas Read, 'Nonabelions in the fractional quantum hall effect,' Nuclear Physics B 360 , 362-396 (1991).
- [26] Chetan Nayak, Steven H. Simon, Ady Stern, Michael Freedman, and Sankar Das Sarma, 'Non-Abelian anyons and topological quantum computation,' Reviews of Modern Physics 80 , 1083-1159 (2008).
- [27] H. P. Wei, D. C. Tsui, M. A. Paalanen, and A. M. M. Pruisken, 'Experiments on Delocalization and University in the Integral Quantum Hall Effect,' Physical Review Letters 61 , 1294-1296 (1988).
- [28] L. Engel, H. P. Wei, D. C. Tsui, and M. Shayegan, 'Critical exponent in the fractional quantum Hall effect,' Surf. Sci. 229 , 13-15 (1990).
- [29] Bodo Huckestein, 'Scaling theory of the integer quantum Hall effect,' Reviews of Modern Physics 67 , 357-396 (1995).
- [30] S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Shahar, 'Continuous quantum phase transitions,' Reviews of Modern Physics 69 , 315-333 (1997).
- [31] J. T. Chalker and P. D. Coddington, 'Percolation, quantum tunnelling and the integer Hall effect,' J. Phys. C 21 , 2665-2679 (1988).
- [32] D.E. Khmelnitskii, 'Quantization of Hall conductivity,' JETP Letters 38 , 552-556 (1983).
- [33] Steven Kivelson, Dung-Hai Lee, and Shou-Cheng Zhang, 'Global phase diagram in the quantum Hall effect,' Physical Review B 46 , 2223-2238 (1992).
- [34] Xiao-Gang Wen and Yong-Shi Wu, 'Transitions between the quantum Hall states and insulators induced by periodic potentials,' Physical Review Letters 70 , 1501-1504 (1993).
- [35] Wei Chen, Matthew P. A. Fisher, and Yong-Shi Wu, 'Mott transition in an anyon gas,' Physical Review B 48 , 13749-13761 (1993).
- [36] Jinwu Ye and Subir Sachdev, 'Coulomb Interactions at Quantum Hall Critical Points of Systems in a Periodic Potential,' Phys. Rev. Lett. 80 , 5409-5412 (1998).
- [37] Tarun Grover and Ashvin Vishwanath, 'Quantum phase transition between integer quantum Hall states of bosons,' Physical Review B 87 , 045129 (2013).
- [38] Maissam Barkeshli and John McGreevy, 'Continuous transition between fractional quantum Hall and superfluid states,' Physical Review B 89 , 235116 (2014).
- [39] M. Barkeshli, N. Y. Yao, and C. R. Laumann, 'Continuous Preparation of a Fractional Chern Insulator,' Physical Review Letters 115 , 026802 (2015).
- [40] Jong Yeon Lee, Chong Wang, Michael P. Zaletel, Ashvin Vishwanath, and Yin-Chen He, 'Emergent multi-flavor QED 3 at the plateau transition between fractional chern insulators: Applications to graphene heterostructures,' Physical Review X 8 , 031015 (2018).
- [41] Yuan-Ming Lu and Dung-Hai Lee, 'Quantum phase transitions between bosonic symmetry-protected topological phases in two dimensions: Emergent QED 3 and anyon superfluid,' Physical Review B 89 , 195143 (2014).
- [42] Ruochen Ma and Yin-Chen He, 'Emergent QCD 3 quantum phase transitions of fractional chern insulators,' Phys. Rev. Research 2 , 033348 (2020).
- [43] P. G. Harper, 'Single Band Motion of Conduction Electrons in a Uniform Magnetic Field,' Proceedings of the Physical Society. Section A 68 , 874-878 (1955).
- [44] M. Ya. Azbel, 'Energy spectrum of a conduction electron in a magnetic field,' Soviet Physics JETP 19 , 634-645 (1964).
- [45] Douglas R. Hofstadter, 'Energy levels and wave functions of Bloch electrons in rational and irrational magnetic fields,' Physical Review B 14 , 2239-2249 (1976).
- [46] R. B. Laughlin, 'Anomalous quantum Hall effect: An incompressible quantum fluid with fractionally charged excitations,' Physical Review Letters 50 , 1395-1398 (1983).
- [47] Rahul Roy, 'Band geometry of fractional topological insulators,' Physical Review B 90 , 165139 (2014).
- [48] David Bauer, T. S. Jackson, and Rahul Roy, 'Quantum geometry and stability of the fractional quantum Hall effect in the Hofstadter model,' Physical Review B 93 , 235133 (2016).
- [49] Bartholomew Andrews and Gunnar M¨ oller, 'Stability of fractional Chern insulators in the effective continuum limit of Harper-Hofstadter bands with Chern number | C &gt; | 1,' Physical Review B 97 , 035159 (2018).
- [50] Thomas Scaffidi and Gunnar M¨ller, 'Adiabatic contino uation of fractional Chern insulators to fractional quantum Hall states,' Physical Review Letters 109 , 246805 (2012).
- [51] T. S. Jackson, Gunnar M¨ oller, and Rahul Roy, 'Geometric stability of topological lattice phases,' Nature Communications 6 , 8629 (2015).
- [52] Bartholomew Andrews, Madhav Mohan, and Titus Neupert, 'Abelian topological order of ν = 2 5 and 3 7 / / fractional quantum Hall states in lattice models,' Phys. Rev. B 103 , 075132 (2021).
- [53] Bartholomew Andrews, Titus Neupert, and Gunnar M¨ oller, 'Stability, phase transitions, and numerical breakdown of fractional Chern insulators in higher Chern bands of the Hofstadter model,' Phys. Rev. B 104 , 125107 (2021).
- [54] Similarly, a parent Chern band could be split into multiple sub-bands by applying periodic superlattice potentials. Care would have to be taken to ensure that the resulting sub-bands do not overlap and gaps open fully.
- [55] G.H. Wannier, 'A Result Not Dependent on Rationality for Bloch Electrons in a Magnetic Field,' Phys. Stat. Sol. B 88 , 757-765 (1978).
- [56] P Streda, 'Theory of quantised Hall conductivity in two dimensions,' J. Phys. C 15 , L717 (1982).
- [57] D. J. Thouless, Mahito Kohmoto, M. P. Nightingale, and M. den Nijs, 'Quantized hall conductance in a twoDimensional periodic potential,' Physical Review Letters 49 , 405-408 (1982).
- [58] Qian Niu, D. J. Thouless, and Yong-Shi Wu, 'Quantized Hall conductance as a topological invariant,' Physical Review B 31 , 3372-3377 (1985).
- [59] D. Osadchy and J. E. Avron, 'Hofstadter butterfly as quantum phase diagram,' Journal of Mathematical Physics 42 , 5665-5671 (2001).
- [60] Steven R. White, 'Density matrix formulation for quantum renormalization groups,' Physical Review Letters 69 , 2863-2866 (1992).
- [61] Ian P. McCulloch, 'Infinite size density matrix renormalization group, revisited,' (2008), arXiv:0804.2509.
- [62] Ulrich Schollw¨ck, 'The density-matrix renormalization o group in the age of matrix product states,' Annals of Physics 326 , 96-192 (2011).
- [63] Johannes Hauschild and Frank Pollmann, 'Efficient nu-

merical simulations with Tensor Networks: Tensor Network Python (TeNPy),' SciPost Physics Lecture Notes 5 (2018), 10.21468/SciPostPhysLectNotes.5.

- [64] Simeng Yan, David A. Huse, and Steven R. White, 'Spin-Liquid Ground State of the S = 1/2 Kagome Heisenberg Antiferromagnet,' Science 332 , 1173-1176 (2011), arXiv:1011.6114.
- [65] Johannes Motruk, Adolfo G Grushin, Fernando De Juan, and Frank Pollmann, 'Interaction-driven phases in the half-filled honeycomb lattice: An infinite density matrix renormalization group study,' Phys. Rev. B 92 , 085147 (2015), arXiv:1505.01674 [cond-mat.str-el].
- [66] Johannes Motruk, Michael P Zaletel, Roger S K Mong, and Frank Pollmann, 'Density matrix renormalization group on a cylinder in mixed real and momentum space,' Physical Review B 93 , 155139 (2016).
- [67] Matthias Gerster, Matteo Rizzi, Pietro Silvi, Marcello Dalmonte, and Simone Montangero, 'Fractional quantum Hall effect in the interacting Hofstadter model via tensor networks,' Physical Review B 96 , 195123 (2017).
- [68] Johannes Motruk and Frank Pollmann, 'Phase transitions and adiabatic preparation of a fractional Chern insulator in a boson cold-atom model,' Physical Review B 96 , 165107 (2017).
- [69] Antoine Sterdyniak, Nicolas Regnault, and Gunnar M¨ oller, 'Particle entanglement spectra for quantum Hall states on lattices,' Physical Review B 86 , 165314 (2012).
- [70] R. B. Laughlin, 'Quantized Hall conductivity in two dimensions,' Physical Review B 23 , 5632-5633 (1981).
- [71] Michael P Zaletel, Roger S K Mong, and Frank Pollmann, 'Flux insertion, entanglement, and quantized responses,' Journal of Statistical Mechanics: Theory and Experiment 2014 , P10007 (2014).
- [72] Adolfo G. Grushin, Johannes Motruk, Michael P. Zaletel, and Frank Pollmann, 'Characterization and stability of a fermionic ν = 1 3 fractional chern insulator,' / Phys. Rev. B 91 , 035136 (2015).
- [73] Pasquale Calabrese and John Cardy, 'Entanglement Entropy and Quantum Field Theory,' Journal of Statistical Mechanics: Theory and Experiment 06002 , 33 (2004), arXiv:hep-th/0405152.
- [74] L. Tagliacozzo, Thiago R. De Oliveira, S. Iblisdir, and J. I. Latorre, 'Scaling of entanglement support for matrix product states,' Phys. Rev. B 78 , 024410 (2008).
- [75] Frank Pollmann, Subroto Mukerjee, Ari M. Turner, and Joel E. Moore, 'Theory of finite-entanglement scaling at one-dimensional quantum critical points,' Physical Review Letters 102 , 255701 (2009).
- [76] S. D. Geraedts, M. P. Zaletel, R. S. K. Mong, M. A. Metlitski, A. Vishwanath, and O. I. Motrunich, 'The half-filled Landau level: The case for Dirac composite fermions,' Science 352 , 197-201 (2016).
- [77] Matthias Gohlke, Ruben Verresen, Roderich Moessner, and Frank Pollmann, 'Dynamics of the KitaevHeisenberg Model,' Physical Review Letters 119 , 157203 (2017).
- [78] We have excluded an outlier for the case of p/q = 2 7 / , L y = 5, because its density profile differed significantly from that of a Laughlin state. We have also excluded results for the highest flux densities, p/q ≥ 3 10 as / experimental results suggest the Laughlin state is destabilized for n φ glyph[greaterorsimilar] 0 4 [84]. .
- [79] Alexei Kitaev and John Preskill, 'Topological Entanglement Entropy,' Physical Review Letters 96 , 110404

## (2006).

- [80] Masudul Haque, Oleksandr Zozulya, and Kareljan Schoutens, 'Entanglement entropy in fermionic laughlin states,' Physical Review Letters 98 , 060401 (2007).
- [81] Andreas M. L¨ auchli, Emil J. Bergholtz, and Masudul Haque, 'Entanglement scaling of fractional quantum Hall states through geometric deformations,' New Journal of Physics 12 , 075004 (2010).
- [82] Michael P Zaletel, Roger S K Mong, and Frank Pollmann, 'Topological characterization of fractional quantum hall ground states from microscopic hamiltonians,' Physical Review Letters 110 , 236801 (2013).
- [83] Shou-Shu Gong, Wayne Zheng, Mac Lee, Yuan-Ming Lu, and D. N. Sheng, 'Chiral spin liquid with spinon fermi surfaces in the spin1 2 triangular heisenberg model,' Phys. Rev. B 100 , 241111 (2019).
- [84] M. Hafezi, A. S. Sørensen, E. Demler, and M. D. Lukin, 'Fractional quantum Hall effect in optical lattices,' Physical Review A 76 , 023613 (2007).
- [85] Frank Pollmann and Ari M Turner, 'Detection of symmetry-protected topological phases in one dimension,' Phys. Rev. B 86 , 125441 (2012).
- [86] L Cincio and Guifr´ e Vidal, 'Characterizing topological order by studying the ground states on an infinite cylinder,' Physical Review Letters 110 , 067208 (2013), arXiv:1208.2623 [cond-mat.str-el].
- [87] Hui Li and F. D. M. Haldane, 'Entanglement Spectrum as a Generalization of Entanglement Entropy: Identification of Topological Order in Non-Abelian Fractional Quantum Hall Effect States,' Phys. Rev. Lett. 101 , 010504 (2008).
- [88] Xiao-Liang Qi, Hosho Katsura, and Andreas W. W. Ludwig, 'General Relationship between the Entanglement Spectrum and the Edge State Spectrum of Topological Quantum States,' Physical Review Letters 108 , 196402 (2012).

Figure 5. Flux insertion for various flux densities n φ . In all cases, L y = 6, V = 10, χ = 600.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-2cfef268-a095-4b69-8c12-475d5af130a3.jpg)

## Appendix A: Details of the Laughlin state in multiple sub-bands

We here present further details to demonstrate the existence and stability of the Laughlin state at flux densities n φ = p/q , p &gt; 1, i.e., where the lowest Landau level (LLL) consists of more than one sub-band. We include cases with p = 2, 3 and 4 sub-bands in the LLL. For each p , we consider several q .

Firstly, Fig. 5 shows the response to a flux pump procedure (as discussed in the main text) for all flux densities, on an L y = 6 cylinder with V = 10. We observe that, although there is some variation in the precise course of the charge transport, in all cases a single charge quantum is transported upon adiabatic insertion of 3 × 2 π flux. Thus, all cases show a Hall response σ H = ∆ q/ ∆ φ ext = 1 3, / which is consistent with a ν = 1 3 Laughlin state. /

Secondly, Fig. 6 shows both the density-density correlation functions as well as the density profile for n φ = 3 / 11 for L y = 10 at V = 10. The density-density correlations reveal a strong correlation hole at small r , followed by rapidly damped oscillations at intermediate r , before the asymptotic value representative of a homogeneous fluid is reached. This is consistent with the Laughlin state, which is a liquid with no long range order, but with strong short-range correlations maximizing typical inter-particle distances. With periodic boundary conditions along ˆ (i.e., on the cylinder geometry), the density y profile is homogeneous (not shown). The inset in Fig. 6 therefore shows the density profile for a strip geometry (i.e., with infinite boundary conditions along ˆ, but open x boundary conditions along ˆ) to demonstrate the effect y of the edges of the system. This density profile shows translational invariance along the cylinder, as well as oscillations around the expected bulk value and overshoot near the edge. While some small finite size effects are visible in the bulk, all these results are consistent with a Laughlin state.

Finally, Fig. 7 shows the scaling of the entanglement entropy with interaction V and bond dimension χ , for

Figure 6. Density-density correlation functions 〈 :ˆ n 0 ,y ˆ n x,y : 〉 along an L y = 10 cylinder. The inset shows the density profile on a strip geometry (open boundary conditions in the ˆ-direction). y

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-27d32a57-6a47-4f0d-94cb-65aa87b75595.jpg)

cylinders with L y = 4 6 7 8 and 10. , , , The singular values s i for any bond in the Matrix Product State give the entanglement entropy as S = -∑ i s 2 i log s 2 i , which we evaluate for a bond that corresponds to a bisection of the cylinder. We note that with increasing L y , it becomes increasingly difficult to resolve the physics at the transition, as evidenced by the strong growth of S χ ( ) with bond dimension. Hence, we cannot fully converge the ground state wave function near the critical point, for L y ≥ 6. In all cases, we observe signatures of a phase transition, where for bigger L y the observed features are less sharp than at smaller ones. Near the transition, there is significant growth of S as function of χ , indicative of a critical point. This implies that DMRG convergence is challenging near the critical point and must be analyzed in terms of finite entanglement scaling. We include additional discussion of the high-bond dimension data for L y = 7 and L y = 8 in appendix D.

## Appendix B: Single-particle density matrix

The single-particle density matrix ρ ij = 〈 ψ c c | † i j | ψ 〉 can be transformed to momentum space using the Fourier transformation of the creation- and annihilation operators. Transforming only the y -coordinate, we have

$$\hat { c } _ { x, y } = \frac { 1 } { \sqrt { L _ { y } } } \sum _ { k _ { y } = 0 } ^ { L _ { y } - 1 } e ^ { i 2 \pi k _ { y } y / L _ { y } } \hat { c } _ { x, k _ { y } }, \quad \ \ ( 2 )$$

and thus

$$\text{ws} _ { \text{os} ^ { - } } \quad \langle c _ { x, k _ { y } } ^ { \dagger } c _ { x ^ { \prime } k _ { y } ^ { \prime } } \rangle = \frac { 1 } { L _ { y } } \sum _ { y y ^ { \prime } } e ^ { i 2 \pi ( k _ { y } y - k _ { y } ^ { \prime } y ^ { \prime } ) / L _ { y } } \left \langle c _ { x y } ^ { \dagger } c _ { x ^ { \prime } y ^ { \prime } } \right \rangle. \ \ ( 3 )$$

We combine this basis transformation with a trivial reordering of the basis to make k y the dominant index, leading to a matrix that is block-diagonal in momentum sectors [Fig. 8(b)]. The spectrum of the matrix remains unchanged throughout this procedure [Fig. 8(a)]. This

Figure 7. Entanglement entropy (Von Neumann) S as function of interaction strength V and bond dimensions χ up to χ = 1000 throughout, and up to χ = 3000 in the transition region, for (a) L y = 4, (b) L y = 6, (c) L y = 7, (d) L y = 8, (e) L y = 10.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-a9671466-686d-454a-8a40-84821f86c64f.jpg)

Figure 8. The single-particle density matrix at L y = 6. (a) The spectrum of the matrix remains unchanged under the transformations described in the main text, and equations 2 and 3. Vertical dashed lines separate every 6 consecutive values. (b) After the final reordering of the indices, the matrix is block-diagonal in k y . Each block contains 11 × 11 values, due to the size (in ˆ of the MUC and MPS unit cell. x

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-25d29260-4c72-4ded-ae55-c3d54a74da71.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-2e6f1427-e3b8-4e7a-af2b-09315d0bd2a8.jpg)

change of basis allows us to classify each eigenvalue by its k y sector, shown in terms of the coloring of symbols in Fig. 2(d) in the main text, and supporting our finding that the number of the most significant, low-lying density-matrix eigenvalues can be mapped to the eigenstates for the lowest 1 /ν single-particle bands.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-a5a437e8-dcaf-4b77-bfec-522b326b29cc.jpg)

-

-

Figure 9. Momentum-resolved entanglement spectrum for a bisected cylinder with L y = 8, χ = 600. Colors denote different MPS quantum number sectors. (a) p/q = 1 3, dashed / lines indicate the semi-degeneracies of the first few levels within the q = 0 sector. (b) p/q = 3 11. / Both spectra show the partition counting of semi-degenerate levels.

## Appendix C: Momentum-resolved entanglement spectra

To further probe the topological nature of the Laughlin state, we compute the momentum-resolved entanglement spectrum (ES) [85, 86] for flux densities that we have studied, examples of which are shown in Fig. 9. Momenta of the Schmidt eigenstates were extracted from the eigenvalues of translations around the cylinder ˆ ( T y a ) | χ R i 〉 = e -i 2 πk /L y y | χ R i 〉 . The calculation proceeds by constructing a mixed transfer matrix of this operator inserted into the ground-state wave function to evaluate the matrix element 〈 χ R i | ˆ ( T y a ) | χ R i 〉 .

Our data clearly show the partition counting associated with the edge modes of the Laughlin state, expected to show up due to the bulk-boundary correspondence relating the ES to the edge modes' energy spectrum [87, 88]. In the lowest-lying charge sector, we find a clear signature of the universal U (1) edge counting of the number of states in consecutive momentum sectors given by the partition counting sequence 1, 1, 2, 3, 5, 7 , . . . . For the charge-neutral sector, our data show reveal the fourfive momentum sectors satisfying this counting of lowlying quasi-degenerate states, before these are no longer clearly distinguishable from the continuum above the entanglement gap. In the adjacent charged charge sectors, the entanglement energies tend to be higher overall, so it becomes more difficult to discern the edge from bulk states.

## Appendix D: Quantum Hall plateau transition at L y = 7 and L y = 8

In the main text, we discuss the plateau transition between a C = +4 Chern insulator and a ν = 1 3 Laugh-/ lin state at n φ = 3 11, at a circumference / L y = 6. At that circumference, the transition appears most clearly. Here, we show additional results for L y = 7 (Fig. 10) and

Figure 10. The plateau transition at L y = 7. Shown are: (a) Correlation length in lattice sites (the circumference log L y ≈ 1 95 falls below the range shown.) . and (b) entanglement entropy for several values of the bond dimension 600 ≤ χ ≤ 1000 throughout, and increasing to χ = 2000 near the transition and up to χ = 3000 between V = 1 5 and . V = 1 6. . (c) Momentum-space ES (colours denote different quantum number sectors). (d) Entanglement gap. (e) Spectrum of the single-particle density matrix. The dashed black and dash-dotted red lines denote the lowest 7 and 21 levels, respectively. (f) Fidelity susceptibility χ F . Subfigures (c) through (f) have χ = 600.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-dc6d01d3-d50f-4600-88bc-98ace166b64c.jpg)

L y = 8 (Fig. 11). The flux insertion procedure has not been repeated for these circumferences due to the large computational cost associated with the required pseudoadiabatic evolution of the wave function. All other observables shown in the main text are included here.

We observe that examined with an MPS of the same bond dimension, the transition appears less sharp at L y = 7 and even less so for L y = 8 than it does at L y = 6. Overall, the transition region is broader than for L y = 6 in both these geometries. However, upon increasing the bond dimension, we see sharper features appear in the entanglement and correlation length at L y = 7.

We first broadly identify the transition region by prob-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-969e724f-a043-425d-9f58-1f5fab49d36d.jpg)

V

Figure 11. The plateau transition at L y = 8. Shown are: (a) Correlation length in lattice sites (the circumference log L y ≈ 2 089 falls below the range shown) and (b) en-. tanglement entropy for several values of the bond dimension 600 ≤ χ ≤ 1000 throughout and up to χ = 2000 near the transition. (c) Momentum-space ES (colours denote different quantum number sectors). (d) Entanglement gap. (e) Spectrum of the single-particle density matrix. The dashed black and dash-dotted red lines indicate the lowest 8 and 24 levels, respectively. (f) Fidelity susceptibility χ F . Subfigures (c) through (f) have χ = 600.

ing for the largest continued increase of the correlation length and entanglement entropy with MPS bond dimension. For L y = 7 [Fig. 10(a-b)], we see that both ξ χ ( ) and S χ ( ) go through a relatively broad plateau for interactions up to V = 1 2 and . V = 1 5, respectively, before . falling off towards the FCI phase. However, on top of this plateau, we find an additional feature linked to a steep rise in both quantities for χ &gt; 1400 and located between V = 1 5 and . V = 1 6. . In order to further investigate this feature, we have run calculations up to a bond dimension of χ = 2000 for the region surrounding the transition, and up to χ = 3000 between V = 1 5 and . V = 1 6. . We have analysed these highχ data to verify

Figure 12. Finite entanglement scaling for the plateau transition at L y = 7, demonstrating the emergence of a linear relation between the entanglement entropy S and the logarithm of the correlation length ξ is consistent with a central charge of c = 2 near the critical point of V c glyph[similarequal] 1 55. . The underlying data include simulations up to bond dimension χ = 3000 at each interaction strength.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00988v2-figure-033bf175-0361-4345-831a-93ba7dfceec0.jpg)

the simultaneous scaling of log( χ ) and S E ( χ ) with the bond dimension, and find a good fit of linear scaling in the throughout the transition region. A detail of these data are is shown in Fig. 12. We see that linear relation corresponding to S E = c/ 6 log( ξ ) provides a good description of the data: The maximum entropy and correlation lengths are seen at V = 1 55, our data point clos-. est to criticality. DMRGconvergence has been noticeably slowed down at this point compared to nearby interaction strengths, and some numerical noise is still visible. In particular, this DMRG run crossed over from a lowerentropy state to a higher entropy one during subsequent increases of the bond dimension. Between χ = 2200 and χ = 2800 the data are fit with a slope of 1 6 × 1 98, while . the χ = 3000 datapoint lies slightly below this line. At the slightly lower interaction of V = 1 5, a linear fit of . S E vs log( ξ ) matches the central charge c = 2 when restricted to bond dimensions below χ = 2000 (data shown in the main text), but then tails off slightly at larger bond dimension. Interestingly, our calculation at the interaction strength V = 1 525 has instead converged to a state . that yields a scaling consistent with c = 1. According to our analysis in terms of a CF Fermi surface state, both c = 1 and c = 2 could describe the critical theory, corresponding to two possible values of 0 or π for the emergent flux felt by composite fermions in the ground state. This emergent flux could be selected randomly depending on the random initial state used in our DMRG runs, so both outcomes are consistent with our proposed scenario for a critical point. We also note that calculations at lower bond dimensions shown in Fig. 10(b) tend to show a dip in the correlation length for V glyph[similarequal] V c , which could be related to the alternative lower central charge c = 1 branch. By contrast, our higher bond dimension calculations overwhelmingly are consistent with c = 2 scaling (cf. Fig. 12), so we believe this reflects the dominant state for the critical behaviour. At larger interaction strength of V = 1 6, we find the slope is noticeably smaller, and is . best fit by a non-integer value of c glyph[similarequal] 1 54, scaling linearly . up to the highest bond dimensions we explored. The non-integer value reflects that this interaction strength is situated outside the critical region. Overall, these findings match the phenomenology for a critical point with central charge c = 2, as we would expect to see continued linear scaling with the bond dimension only precisely at the critical point. With increasing distance from the critical point, we instead expect the growth in entanglement to tail off once the bond dimension is sufficient to accurately capture the finite correlation length of the offcritical system, as seen in these data. Altogether, these features allow us to situate a transition described by a c = 2 central charge in the close vicinity of V = 1 55. .

For L y = 8 [Fig. 11(a-b)], we find that the correlation length goes through a broad maximum which is reached at V = 1 3 to . V = 1 4, depending on the bond dimension . χ . At larger V , we then find a rapid drop of the correlation length, and to the right of the maximum, there is a very strong dependence of the correlation length on χ . Furthermore, the precise location of this drop in ξ , which one may take as an estimate of the numerically determined transition point V c ( χ ), moves towards larger values as the bond dimension is increased. We believe that this shift is due to a competition between the two states for entanglement resources and that it can be explained by noting that the higher-entanglement Chern insulator phase becomes more favourable compared to the fractional Chern phase once the MPS can hold sufficient entanglement to represent it sufficiently accurately. In other words, the DMRG algorithm has a bias towards the FCI phase at low bond dimension, given that the CI state has a higher intrinsic entanglement entropy than the FCI state. At larger V , the correlation length continues to evolve smoothly towards lower values. Unlike the data for L y = 6, or 7, we did not identify any sharp peaks in ξ , but our numerical results are not sufficiently closely spaced on the V -axis to exclude this possibility, and even larger bond dimensions would likely be required to reveal any detailed structure. The entanglement entropy in [Fig. 11(b)] shows very similar behaviour to the correlation length, but the peak is overall slightly broader and values of S start sloping down before undergoing a jump to lower values associated with the FCI phase between V = 1 4 and . V = 1 45. .

The broad behaviour for L y = 7, 8 is similar to what we observed at L y = 6: Within the transition region, both the correlation length ξ χ ( ) and the entanglement entropy S χ ( ) [Fig. 10(a,b), Fig. 11(a,b)] show a continual growth with bond dimension, which is consistent

with a divergence in the bond dimension, as well as with | V -V c | . The region near the transition shows several marked features, with a sharp peak emerging at L y = 7 and the most prominent feature at L y = 8 being a sudden drop in the two quantities towards the right side of the transition. We believe that both are consistent with a continuous quantum phase transition occurring at V c . The scaling of the entanglement entropy with correlation length, shown in Fig. 3 and Fig. 12, further corroborates this view. The data shown in Fig. 3 were taken at the interaction strengths of V L y =6 scale = 1 634, . V L y =7 scale,1 = 1 5 (orange squares) and . V L y =8 scale = 1 45, . respectively. These points are located within the transition region and slightly to the FCI side of an estimated critical interaction strength, where we see clear linear scaling up to χ = 2000. Additionally, we showed the scaling for L y = 7 and V L y =7 scale,2 = 1 55 (orange stars) up to bond . dimension χ = 3000.

The additional data in Fig. 10(c-f) and Fig. 11(c-f) show a number of entanglement spectral properties and the fidelity susceptibility of the ground state wave function for a lower bond dimension of χ = 600. Here, the underlying calculations were performed while semiadiabatically increasing the interaction strength V in order to obtain a consistent allocation of charge quantum numbers between neighbouring interaction strengths (running the DMRG from random initial states can yield different members of the degenerate set of ground states, which have different allocations of charge sectors). We have confirmed that these results are not dependent on the direction of change in V , i.e., that no hysteresis occurs.

Where the momentum-resolved entanglement spectrum exhibited a discontinuous change at L y = 6, for both L y = 7 and L y = 8 we observe a more continuous evolution of the spectrum in the transition region, which also includes some level crossings. For example, in Fig. 10(c) we see the lowest two entanglement eigenstates cross near V = 1 6, and we further clarify that the .

respective states are located in different momentum sectors (while they carry the same charge quantum number), so the crossing is protected by translation symmetry. In Fig. 11(c,e) we note that the entanglement spectrum and single-particle density matrix eigenvalues display a significant change between V = 1 6 and . V = 1 7, and then . remain nearly unchanged for larger values of V . Based on this information, we would locate the transition at a slighly larger value of V c than what we have determined from the entanglement entropy and correlation length.

Sharp discontinuities in the entanglement gap [Fig. 10(d) and Fig. 11(d)] can be ascribed to the gap suddenly being measured from a different ground state (and thus in a different momentum- and quantum number sector), rather than a physical effect. In particular, the level crossing discussed in the previous paragraph explains the sharp feature seen at the corresponding interaction strength in the entanglement gap in Fig. 10(d). The spectra of the single-particle density matrix in [Fig. 10(e) and Fig. 11(e)] show a first gap at a number of eigenstates matching the number of states in the lowest band to the left of the transition, while a second gap emerges adjacent to thrice that number of states to the right of the transition. At L y = 7 the first gap closes smoothly with V and the second one opens gradually above the transition. Al L y = 8, a number of states seem to descend into the first gap on approaching the transition from low V , but the second gap remains at a fairly constant value above the transition. Neither oft these cases show any other prominent gaps, again pointing towards a direct transition from the CI to FCI states. We should caution that these data were taken at bond dimension of χ = 600, and corrections to the higher-lying states are expected at larger bond dimension in the transition region. Results for the fidelity susceptibility χ F [Fig. 10(f) and Fig. 11(f)] again do not demonstrate sharp features, but do suggest a broad region where the ground state varies rapidly with the interaction V .
Leon Schoonderwoerd Frank Pollmann Gunnar Möller
Interaction-driven plateau transition between integer and fractional Chern Insulators
We present numerical evidence of an interaction-driven quantum Hall plateau
transition between a $|C|>1$ Chern Insulator (CI) and a $\nu = 1/3$ Laughlin
state in the Harper-Hofstadter model. We study the model at flux densities
$p/q$, where the lowest Landau level (LLL) manifold comprises $p$ magnetic
sub-bands. For weak interactions, the model realises integer CIs corresponding
to filled sub-bands, while strongly interacting candidate states include
fractional quantum Hall (FQH) states at LLL filling fractions $\nu=r/t$. These
phases may compete at the same particle density when $p=t$. As a concrete
example, we numerically explore the physics at flux density $n_{\phi} = 3/11$,
where we show evidence that a direct transition occurs between a CI and a $\nu
= 1/3$ Laughlin state, which we characterise in terms of its critical,
topological and entanglement properties. We also show that strong interactions
generically stabilise a $\nu = 1/3$ Laughlin state even when the LLL is split
into multiple bands, and introduce a powerful methodology to extract its
topological entanglement entropy by exploiting the scaling of magnetic length
with $n_\phi$.