Quasi-one-dimensional Hall physics in the Harper-Hofstadter-Mott model

We study the ground-state phase diagram of the strongly interacting Harper-Hofstadter-Mott model at quarter flux on a quasi-one-dimensional lattice consisting of a single magnetic flux quantum in $y$-direction. In addition to superfluid phases with various density patterns, the ground-state phase diagram features quasi-one-dimensional analogues of fractional quantum Hall phases at fillings $\nu=1/2$ and $3/2$, where the latter is only found thanks to the hopping anisotropy and the quasi-one-dimensional geometry. At integer fillings - where in the full two-dimensional system the ground-state is expected to be gapless - we observe gapped non-degenerate ground-states: At $\nu=1$ it shows an odd 'fermionic' Hall conductance, while the Hall response at $\nu=2$ consists of the transverse transport of a single particle-hole pair, resulting in a net zero Hall conductance. The results are obtained by exact diagonalization and in the reciprocal mean-field approximation.


Introduction
The prospect of realizing and measuring topologically non-trivial bosonic phases remains an intriguing and important challenge of condensed matter physics. The bosonic statistics can lead to different topological properties than the ones observed in fermionic systems [1][2][3][4][5][6]. Furthermore, as bosons condense in the absence of interactions, topologically non-trivial bosonic phases are inherently many-body in nature, as an interaction is needed to introduce a gap. From an experimental point of view, bosonic atoms are easier to control in cold atom experiments with optical lattices [7][8][9][10][11], making them prime candidates in the search for interacting topological properties.
While both ED and DMRG have provided great insight in such models, they rely on finite system sizes. In two dimensions, these may be too small for ED whereas DMRG has difficulties converging for gapless phases due to rapid entanglement growth; i.e., there is a preference for gapped, low entanglement phases. In order to obtain the full phase diagram, we see therefore a need to develop new methods, such as the recently proposed reciprocal cluster mean-field (RCMF) method, which has also been applied to the two-dimensional HHMm [22], and which tends to favor the opposite: it is defined in the thermodynamic limit and favors condensed phases (although it certainly can find topologically non-trivial phases as we will see). A systematic comparison of the two approaches therefore provides a promising path towards understanding the ground-state properties of the

HHMm model on a cylinder
The Hamiltonian of the HHMm can be written as U n n n e h . c . 2 1 , 1 x y x y x y x y y x y x y x y x y x y x y where the coordinates (x, y) parameterize a system of size L x ×L y , with periodic boundary conditions. The operators b x y , ( †) are the annihilation (creation) operators at site (x, y), and n b b x y x y x y , , , = † is the occupation number operator. The hopping amplitudes in xand y-direction are t x and t y , respectively, Φ is the flux through each plaquette, U is the strength of the on-site interaction, and μ is the chemical potential. Throughout our analysis we focus on the hard-core boson limit, U  ¥, and the flux Φ=π/2. We are left with two dimensionless parameters, the hopping anisotropy t t t x y max -( ) , and the chemical potential μ/t max , where t t t max , x y max = { }. Furthermore, we focus on the cylinder-geometry in the thermodynamic limit, with L x  ¥ and L y =4. In the case Φ=π/2 the magnetic unit-cell is of size 1×4 for the Landau gauge used in equation (1). This choice of the lattice size therefore makes the cylinder quasi-one-dimensional, as only one unit-cell is present in the ydirection, and different ground-state phases than those of the fully two-dimensional HHMm [2, 4, 18-20, 22, 35-37] can be expected. Note that even though the magnetic unit cell can be chosen to consist of less than 4 sites in y-direction, e.g. of size 2×2 [38], the particles would still need to perform a (gauge-invariant) loop around 4 plaquettes in y-direction in order to pick up a trivial phase of 2π. Having just two plaquettes in y-direction on the other hand, would correspond to just half of the flux quantum 2π. In this sense our approach to quasi-one-dimensionality is different to the one of [26,39,40], where there is just one plaquette in y-direction on a ladder geometry.
The Hamiltonian (1) can be block-diagonalized in the hard-core boson limit with the Fourier transform This Hamiltonian can be further written in terms of the h-vector [22], which measures the Hall response-see appendix A for details. Note that the Hamiltonian (1) in the hard-core boson case is invariant-up to a constant-under the charge-conjugation transformation . 3 x y x y , , This implies that the ground-states at positive and negative chemical potentials (or equivalently at densities n and 1−n), are related by the holes taking on the role of particles and the Hall response changing sign (σ xy a−σ xy ) [22,41]. In the following we will therefore restrict ourselves to the case of μ0 (i.e. n1/2), with the phases at positive chemical potentials related to the ones discussed in this work by equation (3).

RCMF theory
The main method we employ for the analysis of the model is RCMF [22], whose results we will further support by ED. It is defined in the thermodynamic limit and variationally approaches both condensed and uncondensed phases in models with multiorbital unit-cells and non-trivial dispersions, while preserving the translational symmetry of the lattice. Within this method, a Hamiltonian such as (1) in the thermodynamic limit is decoupled into a set of identical clusters with size C x ×C y (here C x =C y =4) through a combination of momentum coarse-graining [42] and the mean-field decoupling approximation [43]. is a generalized N Φ ×N Φ dispersion of the unit cell and k, q are momenta in x, y direction, respectively. We start by dividing the lattice of size L x ×L y into L L C C x y x y clusters of size C x ×C y . As is done also in the dynamical cluster approximation [42], we decompose the position coordinates as x y X x Y y , , = + + ( ) (˜˜)and the quasi-momenta as where (X, Y), (K, Q) are the intra-cluster components, and x y k q , , , (˜˜) (˜˜) are the inter-cluster components. The dispersion can now be rewritten as , and the inter-cluster term capturing hopping processes between different clusters with the remainder of the dispersion K k Q q , In the next step we perform the mean-field decoupling in the inter-cluster term, expanding the bosonic k q ,0 ,0 is the condensate order parameter in momentum space. In the second equality of equation (7), we have assumed that the clusters are large enough, such that the minima of the non-interacting dispersion can be reproduced by the quasi-momenta of the C x ×C y clusters (i.e. n C m C K 2 ,2 n m x y , p p = ( )). In this case, for local translational-invariant interactions, we can safely assume to have condensation only in the cluster momenta (i.e. at k q 0 = =˜).
in the inter-cluster couplings, we find the effective cluster-local RCMF Hamiltonian, Here, t R R , ¢ is the hopping-amplitude of the coarse-grained dispersion from equation (6) which can be solved iteratively. For a more detailed derivation of the method and technical details, see [22].
It should be noted that, as the approximation neglects quadratic fluctuations between particles on neighboring clusters (which we know matter at this level of accuracy, see [44,45]), RCMF tends to systematically overestimate the condensate order parameter at the expense of gapped phases. This can probably be improved, but is left for future work. Note, however, that RCMF reduces to ED (up to a renormalization of the parameters originating from coarse graining in momentum space) where it finds gapped phases: it is thus well equipped to find non-trivial phases and only approaches the thermodynamic limit in a non-conventional but mathematically sound way. Since conventional ED lacks U(1)-symmetry-breaking, regions of the phase diagram where both methods agree are therefore a strong indication of trustworthiness.

Results
The ground-state phase diagram of the quasi-one-dimensional model computed with RCMF (see section 3) is shown in figure 1 in terms of the hopping anisotropy t t t x y max -( ) and the chemical potential μ/t max . In the limit t t t 1 x y max -= ( ) (i.e. t y =0), the system consists of four decoupled infinite chains each exhibiting an uncondensed one-dimensional superfluid (or band-insulating (BI)) ground-state. At t t t 1 . t x = 0) instead, the cylinder reduces to an infinite set of decoupled 4-site rings which-depending on the filling -can exhibit gapped phases. The phase diagram therefore strongly differs at negative anisotropies from the one observed in the fully two-dimensional model [22]. The fact that the phases discussed in the following are in part entirely related to the quasi-one-dimensional geometry of the lattice is further evidenced by the fact that the situation changes drastically as soon as L y is changed from 4 to 8 with two unit cells in the y-direction, as discussed in appendix B, where the L y =8 results are much closer to the fully two-dimensional results than the ones for L y =4.
Before discussing the different phases in more detail, let us illustrate the role of sign asymmetry in the hopping anisotropy as well as perform a benchmarking of RCMF against ED by studying the density as a function of chemical potential, shown in figure 2. If μ is sufficiently negative the model is in a trivial BI phase with zero particles for either anisotropy. For hopping anisotropy t t t 0.35 , RCMF always finds a condensed ground-state except in a narrow region at density n=1/2. The ED results (which can only exhibit plateaus at densities commensurable with the system size and are computed by comparing the grand-canonical energies of the respective particle-number sectors) tend towards the continuous RCMF results as L x is increased, as all observed plateaus quickly shrink with system size (again excluding the narrow region at density n=1/2, which remains gapped in either method). For hopping anisotropy t t t 0.35 figure 2(b)), the agreement between RCMF and ED increases further, as the RCMF ground-state now also shows plateaus at integer (ν=1, 2, i.e. n=1/4, 1/2) and fractional (ν=1/2, 3/2, i.e. n=1/8, 3/8) fillings. The fractional plateaus, which look almost system size independent in ED, are however smaller in RCMF. In the limit of low densities however, we are able to extend the RCMF cluster from a size of 4×4 to 8×4 (see appendix C). As can be seen in figures 5(a) and C1 for the plateau at ν=1/2 this increases the accuracy substantially, finding excellent agreement with ED.

Symmetry-broken phases
The ground-state away from integer (ν=1, 2) and half-integer (ν=1/2, 3/2) filling is always condensed, exhibiting different density and condensate modulations as a function of chemical potential and hopping anisotropy. The three resulting phases we will discuss in the following exhibit first order phase transitions.
At positive anisotropy a large part of the phase diagram is occupied by the vertically striped superfluid (VSF) where both the total and the condensate density exhibit vertical stripes (see figure 3(a)). With increasing negative anisotropy, these patterns are rotated by π/2 and the superfluid becomes horizontally striped (HSF). Finally, in the vicinity of filling ν=1/2 the ground-state exhibits a checkerboard superfluid pattern, apparently due to doping mechanisms on top of the degenerate ground-state (CS, see figure 3(b)). ) . The observed phases are: band insulator (BI, light blue), gapped phases at integer fillings ν=1, 2 (Int., green), quasi-one-dimensional analogs of fractional quantum Hall phases at fillings ν=1/2, 3/2 (fQH, dark gray), gapless liquid (liquid, pink), and superfluid phases with different patterns: vertically striped (VSF, white), horizontally striped (HSF, brown) or with checkerboard order (CS, dark blue). The phases in the purely one-dimensional case are an uncondensed superfluid at t t t 1 , and gapped decoupled 4-site rings at t t t 1 x y max -=-( ) (except for the band-insulator at density n = 0). In the dashed region around zero chemical potential the results are inconclusive whether the phase is gapped (Int.) or gapless (liquid). The fractional plateau at ν=1/2 was computed using an 8×4 cluster, unlike the results for higher chemical potentials/ densities which were computed using a 4×4 cluster, see also appendix C. The plateau at ν=3/2 is most likely larger than the one found here (see figure 5(c)), however both RCMF and ED are inconclusive on the exact phase boundaries within the accessible cluster-/system-sizes. Comparison between results computed with RCMF (black) and ED on system sizes 4×4 (red) and 8×4 (blue) with periodic boundary conditions. Particle density n as a function of chemical potential μ/t max at hopping anisotropy t t t 0.35

U(1)-symmetry preserving phases
In the absence of a finite condensate order parameter R f , the Hamiltonian (8) reduces to a finite U(1)symmetry-preserving 4×4 torus with periodic boundaries. We can therefore turn to ED in order to further analyze the properties of U(1)-symmetry preserving phases (the only difference is a rescaling of the energy unit).
In order to extrapolate to the thermodynamic limit L x  ¥ we apply the twisted boundary conditions This allows us to estimate the many-body gap as E min , 11 where ε GS and ε 1 are the energies of the ground-state(s) and the first excited state, respectively. In addition, by analyzing the behavior of the quasi-one-dimensional h-vector [22] as a function of the twisting angle θ x we can extrapolate the topological properties of the system. If the h-vector shows a closed-loop as a function of θ x , and the ground-state does not mix with the excited states such that (11) stays finite, it implies that the many-body ground-state is adiabatically translated by a single site in y-direction during one charge-pump cycle, resulting in different quantized Hall conductances depending on the filling. For details on twisted boundaries and the hvector, see appendix A.

Integer filing
At (low enough) negative hopping anisotropy we observe non-degenerate ground-states at integer fillings ν=1, and 2 (see 'Int.' in figure 1). In figures 4(b) and (d) we show the many-body gap computed with ED for these fillings as a function of hopping anisotropy. The gap computed with (11) on a 4×4 system, as well as the scaling when going from L x =4 to L x =8 clearly indicate a gapped ground-state. Furthermore, as shown in figures 5(b) and (d), the single-particle gap computed with ED agrees perfectly with the region where we observe the gapped integer filling phases in RCMF. For ν=2 a known candidate for gapped non-trivial many-body states in the HHMm is the biQH phase with transverse conductance σ xy =2 [19,28,38,[46][47][48]. While this phase has been found in the HHMm with hardcore bosons at lower fluxes [19], the case of Φ=π/2 is special, as filling ν=2 corresponds to a density of n=1/2, where hard-core bosons show a  -symmetry (as discussed in section 2). We checked numerically that this symmetry is not spontaneously broken, by computing the overlap of the ED ground-state with its  -transform (i.e. the complex-conjugated ground-state after a particle-hole transform), yielding always 1. The Hall conductance of this ground-state can therefore only be σ xy =0.
At filling ν=1, a biQH phase can only be expected in the presence of two different bosonic species [49] (i.e. filling ν=1+1). Another candidate for gapped phases is the non-abelian Moore-Read state [50,51], which however is characterized by a degenerate ground-state, not observed here.
As discussed in appendix A, the h-vector winds once around the origin as a function of the twisting angle for both fillings. For ν=2 this implies the quantized transverse transport of a single particle-hole pair during one charge-pumping cycle, resulting in a total Hall conductance of σ xy =0, consistent with the charge conjugationsymmetry of (3). For ν=1, on the other hand, this implies the transport of a single particle, and thereby a Hall conductance of σ xy =1, as would be observed in a fermionic integer quantum Hall effect [52].
Especially the latter phase may seem surprising, as such a bosonic phase with odd Hall conductance is expected to be either gapless, or show intrinsic topological order and thereby a degenerate ground-state [28,29,50,51]. However, while in the two-dimensional lattice phases with f=0 at these fillings are found to be always gapless liquids [22], in the quasi-one-dimensional setup these are connected to the limit of decoupled

4-site rings at t x =0
, where hard-core bosons behave as free fermions gapped by the finite size of the rings, providing an intuitive explanation for the 'fermionic' behavior of the ν=1 phase. In fact, the liquid phase of the two-dimensional model [22]-just as metallic Fermi-liquid-like phases of hard-core bosons predicted in the lowest Landau level [53][54][55]-shows an average Hall response of 1 which is however not protected against disorder due to its gapless nature. Here, it appears that this phase is gapped through the finite-size in y-direction. Eventually, when going to positive hopping anisotropies, the gaps computed with (11) and the strongly sizedependent scaling of the gaps with periodic boundaries shown in figures 4(b) and (d) indicate a gapless groundstate for both fillings. These phases are equivalent to the anisotropic gapless liquid observed in the fully twodimensional model [22]. In fact, at positive anisotropies, those phases are connected to the limit of decoupled infinite chains at t y =0, where the hard-core bosons are in a superfluid ground-state.

Fractional filling
In the region of negative hopping anisotropies we observe commensurate phases with at fillings ν=1/2, and 3/ 2. For both fillings the system is characterized by a two-fold degenerate ground-state, and is gapped for negative (and low positive) anisotropies, as evidenced by the ED results shown in figures 4(a) and (c).
As for the integer filling phases, the h-vector shows a closed loop as a function of the twisting angle (see appendix A) indicating Hall conductances of σ xy =1/2, and 3/2, respectively. In the two-dimensional case a fQH phase at filling ν=1/2 has already been observed in previous ED [2,4,18], variational Gutzwiller mean field [37], and DMRG [19,20] studies, while the one at ν=3/2 has not been observed in the two-dimensional limit. As for the integer phases, which are only gapped on the cylinder, the fractional phases therefore appear to be much more stable on a quasi-one-dimensional geometry. The observed phases on the cylinder differ from their isotropic counterpart in the two-dimensional continuum by a weak charge density wave order (see ) computed using twisted boundaries on a 4×4 system (black dashed), periodic boundaries on a 4×4 system (red), and periodic boundaries on a 8×4 system (blue) for fillings ν=1/2 (a), ν=1 (b), ν=3/2 (c), and ν=2 (d).
appendix A). However, while for the integer filling phases the many-body gap closes for twisting in the ydirection and yx xy s s ¹ is not quantized, the fractional phases stay gapped for twisting in both directions, and therefore show a truly two-dimensional Hall response σ yx =σ xy . We conclude that these phases are quasi-onedimensional analogs of fQH phases, which continuously develop into their two-dimensional counterparts as the circumference of the cylinder is increased [23][24][25][26]. A similar case has been discussed in the 1+1-dimensional two-leg ladder with an additional external parameter in [26].
While showing the same general trend, the single particle gaps predicted by ED for the fQH phases in figures 5(a) and (c) are significantly larger than the ones predicted by RCMF using a 4×4 cluster, which tends to underestimate long-range-entangled phases. However, when employing an 8×4 cluster for low densities (see appendix C), the results on the ν=1/2 plateau show excellent agreement with ED. Unfortunately, for the ν=3/2 plateau such a cluster size is out of reach for RCMF, and also ED appears to not be converged with system size, such that we can only conclude that the phase boundaries of the fractional phase at ν=3/2 lie in between the ones found with ED and RCMF.
In the isotropic two-dimensional system, the fractional phase at filling ν=1/2 predicted by other methods [2,4,[18][19][20]37] is found to be slightly metastable in RCMF [22]. However, seeing the scaling when increasing C x to 8, we expect that an 8×8 cluster would be needed to fully capture the ground-state behavior at ν=1/2 in the two-dimensional case. Finally, we note that both our RCMF results and our ED results do not point to a gapped phase at filling ν=2/3 claimed in [19] (see also appendix D). Figure 5. Regions of the ground-state phase diagram computed with RCMF where we observe non-trivial phases with f=0 (gray), compared with the single particle gaps measured with ED using periodic boundaries on system sizes 4×4 (red) and 8×4 (blue) for fillings ν=1/2 (a), ν=1 (b), ν=3/2 (c), and ν=2 (d). In panel (a) the RCMF plateau is computed using an 8×4 cluster, while the RCMF results using a 4×4 cluster (which is employed for all other fillings) are shown as a black dashed line.

Conclusion
In this work we studied the ground-state properties of the HHMm model with quarter-flux per plaquette and hard-core bosons on a quasi-one-dimensional lattice consisting of a single flux quantum along the y-direction.
We found that the stability of topologically non-trivial phases is significantly enhanced by the quasi-onedimensional geometry. The ground-state phase diagram features quasi-one-dimensional analogs of fractional quantum Hall phases at fillings ν=1/2 and ν=3/2, where the latter was not observed in the two-dimensional system. We further observed new gapped non-degenerate ground-states at fillings ν=1 (characterized by an odd 'fermionic' Hall conductance of σ xy =1) and at filling ν=2-with total zero Hall conductance, but characterized by the quantized transverse transport of a single particle-hole pair as Hall response.
By systematically comparing RCMF and ED, which approach the thermodynamical limit from opposite sides (the first method favors gapless, the second one gapped phases), we are able to give conclusive quantitative answers on the phase boundaries of gapped phases.
These unconventional integer filling phases-which do not exist in the two-dimensional case-illustrate the peculiarity of quasi-one-dimensional geometries in topological systems. The increased stability of the gapped phases in this setup-which could be realized by mapping the finite spatial direction y onto internal degrees of freedom-could facilitate the experimental search for topologically non-trivial bosonic phases.
By examining the spectrum on a finite system, twisted boundaries offer a reasonable tool to compensate for finite-size effects in ED calculations regarding the robustness of spectral gaps. In our calculations on the 4×4 torus, the y-direction is treated exactly (as L y =4). We can therefore only introduce a twisted boundary angle in the x-direction, θ x , defined as Figures A1 and A2 show the dependency of the spectrum on the twisting angle θ x for fillings ν=1/2, 1, 3/2 and 2 at t t t 0.35 ) and t t t 0.35 , respectively. The six lowest eigenvalues are shown, as computed with ED on the 4×4 lattice. When the spectrum mixes, the many-body ground-state on the lattice in the thermodynamic limit can be assumed to be gapless, see equation (11).
In figures A1 and A2 we see that the ground-state at fractional fillings is 2-fold degenerate. Furthermore, we observe how the ground-states mix with the excited ones as a function of θ x at t t t 0.35 ) and fillings ν=1 and ν=3/2, while the other cases shown (i.e., ν=1/2 and ν=2 in figure A2 and all negative values of t t t 0.35 figure A1) appear to be gapped. Let us now rewrite the Hamiltonian (2) into the form used in [22] suited to introduce the h á ñ -vector. To this end, we define . This means that the many-body ground-state is adiabatically translated by a single site in y-direction. For a single charge-pump process ( : 0 2 x q p  ) and filling ν (i.e. total number of particles N=ν L x for L y =4) this amounts to the transverse transport of N=ν 4 particles by one site, or equivalently ν particles through the full system (i.e. by L y =4 sites) in y-direction.
The winding shown in figure A3 for ν=1/2 therefore implies a Hall conductance of σ xy =1/2. Similarly, ν=1 and ν=3/2 correspond to σ xy =1, and σ xy =3/2, respectively. The winding in the ν=2 case is special, due to the charge-conjugation symmetry at n=1/2 (see equation (3)). For each particle being transported there is also a hole being transported, resulting in a total quantized transport of a particle-hole pair, and zero Hall conductivity.
Unlike at fractional fillings, for the integer filling phases-as their many-body gap is a direct consequence of the anisotropic geometry-the Hall conductance is highly anisotropic, in the sense that it is quantized for charge-pump processes in the x-direction only. As the ground-state mixes with the excited states as a function of twisting in y-direction (as discussed in [22]), yx xy s s ¹ is not quantized, emphasizing the quasi-onedimensional nature of these phases. It should be noted that if the y-direction is mapped onto four different species/internal states, such that d k d k    In order to fully classify the gapped phases at fractional filling, we turn to the density-density correlations in x-direction n n n n x y y x y y , 0, , 0, á ñ-á ñá ñ | | . As shown in figure A4 for the two degenerate ground-states at filling ν=1/2 and t x =t y for y=0, these correlations indicate a charge density wave order: after decreasing for short distances they quickly saturate to a small finite value as a function of x and appear to be converged in the system size L x for L x =10 and 12. A similar behavior is also observed for other values of y and at filling ν=3/2. Unlike their two-dimensional fQH counterparts in the continuum, these fractional phases are therefore characterized by a weak charge density wave order. The Hall response, on the other hand, is fully isotropic σ yx =σ xy =ν, just as in the two-dimensional case.

Appendix B. Quasi-one-dimensional versus two-dimensional geometry
The quantity which differentiates between the different lattice geometries in RCMF is the coarse-grained dispersion on the 4×4 cluster, i.e.
The general two-dimensional system can be written as H H k k = å , where k=(k, q), k is the quasimomentum in x-direction, q the quasi-momentum in y-direction, and In the thermodynamic limit in x-direction, assumed for all geometries (L x  ¥), k 0, 2p Î [ ] is continuous. For L y =4 considered in the main text, q=0, such that the coarse-grained dispersion is given by a one-dimensional integral k d

( )
In figure B1 we compare the three geometries (L y =4, L y =8 and L y  ¥). We show the condensate density n c C C R R 1 2 x y f = å | | for the three different lattices as a function of chemical potential, illustrating how the L y =8 cylinder with just two unit-cells in y-direction already shows much better agreement with the gapless two-dimensional results [22], further indicating the exclusive quasi-one-dimensional nature of the L y =4 cylinder.
Appendix C. RCMF scaling at low densities As mentioned in section 3, for the RCMF approximation to work, it is indispensable that the momentum values of the four minima of the non-interacting dispersion at k=0, π,±π/2 can be represented within the C x ×4 cluster impurity. After the 4×4 cluster employed in this work the next cluster size which is consistent with this approach would be therefore 8×4, which in the absence of particle number conservation is computationally out of reach.
In the case of low density however we can restrict the basis of the 8×4 cluster to the particle number sectors N=0, 1, 2, K, 6, making the search for a symmetry-broken ground-state numerically accessible. The resulting hopping amplitudes with respect to the 4×4 cluster change accordingly to In figure C1 we compare results computed with such a restricted basis using an 8×4 cluster in the vicinity of density n=1/8 (i.e. ν=1/2), where such an approach is still controlled. We see how ED and RCMF scale differently with cluster size, converging towards each other: since ED prefers gapped phases, the fractional plateau at ν=1/2 shrinks with cluster size. As RCMF prefers gapless phases, the fractional plateau increases with cluster size (see also figure 5(a)). Figure B1. Condensate density n c as a function of chemical potential μ/t max for different lattices, i.e. L y =4 (solid, black), L y =8 (dashed, blue) and L y  ¥ (dotted, red). In (a) the hopping anisotropy is t t t 0.35 Jain's sequence [56], a composite fermion approach to the HHMm predicts a series of gapped phases at fillings ν=p/(p+1) for integer p. In this picture the fQH phase at filling ν=1/2 is equivalent to the first state of Jain's sequence at p=1. The bosonic quantum Hall phase at ν=2 with σ xy =2 is equivalent to p=−2. Note that it does not exist for Φ=π/2 with hard-core bosons (the charge-conjugation symmetry of equation (3) imposes σ xy =0 at ν=2) but was measured for lower fluxes in [19].
Reference [19] further reports on a gapped phase at filling ν=2/3 and no hopping anisotropy at flux Φ=π/2, equivalent to the next state in Jain's sequence after the ν=1/2 fQH (i.e. p=2), i.e. Halperin's 211 ( ) state, which we do not observe in RCMF. We therefore turn to the ED spectrum of system sizes L x ×4 with L x =6, 9 and 12 to look for signatures of a gapped phase at filling ν=2/3 in figure D1(a). As can be seen for both L x =6 and L x =12 we do not observe a clear indication of a three-fold degeneracy of the ground-state expected for Halperin's 211 ( ) state. In figure D1(c) the distance between the lowest and the degenerate second and third eigenstate is shown as a function of L x . While it generally decreases (and goes almost to zero at L x =9), so does the gap between the third and the lowest eigenstate (which would be the many-body gap in case the ground-state were truly three-fold degenerate). For comparison, we show in figure D1(b) the spectrum for ν=1/2 for the same system sizes. Here, as expected, the ground-state is two-fold degenerate, and the many-body gap to the first excited state is considerably larger (around 1/4 of the hopping). In line with RCMF, the ED results do not seem to point toward the existence of a gapped state at filling ν=2/3. Figure C1. Comparison of particle density (a) and condensate density (b) as a function of chemical potential for t x =t y computed with ED on a 4×4 torus (black), ED on an 8×4 torus (green), RCMF with a 4×4 cluster impurity (blue) and RCMF with an 8×4 cluster impurity (red). Figure D1. Comparison of the ED spectra at fillings ν=2/3 and ν=1/2. (a) and (b) Eigenstates as a function of momentum sector K at t x =t y and system sizes 6×4 (red) and 12×4 for fillings ν=2/3 (a) and ν=1/2 (b). (c) Many-body gap of the fractional phase at ν=1/2 (red) and gaps at ν=2/3 between the lowest and the first excited state (blue, dashed) and between the lowest and the third excited state (black) as a function of L x .