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

Fractional quantum Hall effect at ν=2+4/9\nu=2+4/9

Ajit C. Balram1 and A. Wójs2 1Institute of Mathematical Sciences, HBNI, CIT Campus, Chennai 600113, India 2Department of Theoretical Physics, Wrocław University of Science and Technology, 50-370 Wrocław, Poland
Abstract

Motivated by two independent experiments revealing a resistance minimum at the Landau level (LL) filling factor ν=2+4/9\nu=2+4/9, characteristic of the fractional quantum Hall effect (FQHE) and suggesting electron condensation into a yet unknown quantum liquid, we propose that this state likely belongs in a parton sequence, put forth recently to understand the emergence of FQHE at ν=2+6/13\nu=2+6/13. While the ν=2+4/9\nu=2+4/9 state proposed here directly follows three simpler parton states, all known to occur in the second LL, it is topologically distinct from the Jain composite fermion (CF) state which occurs at the same ν=4/9\nu=4/9 filling of the lowest LL. We predict experimentally measurable properties of the 4/94/9 parton state that can reveal its underlying topological structure and definitively distinguish it from the 4/94/9 Jain CF state.

pacs:
73.43-f, 71.10.Pm

The fractional quantum Hall effect (FQHE) Tsui et al. (1982); Laughlin (1983) forms a paradigm in our understanding of strongly correlated quantum phases of matter. Of particular interest among the panoply of FQHE phases are the ones observed in the second Landau level (SLL) of ordinary semiconductors such as GaAs. These have attracted widespread attention because of the possibility that the excitations of these phases obey non-Abelian braiding statistics, which could potentially be utilized in carrying out fault-tolerant topological quantum computation Kitaev (2003); Nayak et al. (2008).

FQHE was first observed at filling factor ν=1/3\nu=1/3 Tsui et al. (1982) in the lowest LL (LLL) and was explained by Laughlin using his eponymous wave function Laughlin (1983). Soon a whole zoo of fractions were observed, primarily along the sequence n/(2pn±1)n/(2pn\pm 1) (nn and pp are positive integers) and its particle-hole conjugate Eisenstein et al. (1990). These FQHE states can be understood as arising from the integer quantum Hall effect (IQHE) of composite fermions (CFs) Jain (1989a, 2007), which are bounds states of electron and an even number (2p2p) of quantized vortices. The theory of weakly interacting CFs captures almost all of the observed FQHE phenomenology in the LLL.

In comparison to the FQHE in the LLL, the FQH states in the SLL are fewer in number and are more fragile Willett et al. (1987). Moreover, the nature of many FQH states in the SLL is dramatically different from their LLL counterparts. In particular, one of the strongest FQH states in the SLL occurs at ν=2+1/2\nu=2+1/2 Willett et al. (1987), whereas at the corresponding half-filled LLL a compressible state is observed. A key breakthrough in the field of FQHE came about by a proposal of Moore and Read Moore and Read (1991), who posited a “Pfaffian” wave function to describe the state at ν=5/2\nu=5/2. Subsequently, it was understood that the Pfaffian wave function can be interpreted as a pp-wave paired state of composite fermions Read and Green (2000). The excitations of the Pfaffian state are Majorana fermions which, owing to their non-Abelian braiding properties, could form building blocks of a topological quantum computer Kitaev (2003); Nayak et al. (2008).

The nature of the state at 1/31/3 state in the SLL, although believed to be Laughlin-like, has been under intense debate d’Ambrumenil and Reynolds (1988); Balram et al. (2013); Johri et al. (2014); Peterson et al. (2015); Kleinbaum et al. (2015); Jeong et al. (2017); Balram et al. (2020). FQHE has been observed at ν=2+2/5\nu=2+2/5 Xia et al. (2004); Choi et al. (2008); Pan et al. (2008); Kumar et al. (2010) but is widely believed to be of a parafermionic nature, unlike the Abelian LLL Jain CF state at ν=2/5\nu=2/5 Read and Rezayi (1999); Rezayi and Read (2009); Wójs (2009); Bonderson et al. (2012); Sreejith et al. (2013); Zhu et al. (2015); Mong et al. (2017); Pakrouski et al. (2016). Furthermore, as yet there is no conclusive experimental evidence of FQHE at the next three members of the n/(2n+1)n/(2n+1) Jain sequence, namely ν=3/7\nu=3/7, 4/94/9, and 5/115/11, in the SLL Shingla et al. (2018), though some features of FQHE have been reported in the literature at some of these fillings Pan et al. (2008); Choi et al. (2008); Reichl et al. (2014) (A non-Abelian parton state at 3/73/7 was constructed in Ref. Faugno et al. (2020) and shown to be feasible in the SLL.). However, FQHE has been observed at ν=2+6/13\nu=2+6/13 Kumar et al. (2010). These observations collectively point to the fact that the nature of FQHE states in the LLL and SLL are different from each other, and a description of FQHE states in the SLL likely entails going beyond the framework of noninteracting composite fermions.

Recently, candidate “parton” states have been constructed to describe the FQHE at many of the experimentally observed fillings in the SLL Balram et al. (2018a, b, 2019). The parton theory Jain (1989b) produces model incompressible states beyond the Laughlin Laughlin (1983) and the CF theory Jain (1989a). The “n¯2¯111\bar{n}\bar{2}111” parton states for n=1,2,3n=1,2,3 give a good description of the SLL FQHE states observed at ν=2+2/3,2+1/2\nu=2+2/3,2+1/2 and 2+6/132+6/13 Balram et al. (2018a, b). In particular, the n¯2¯111\bar{n}\bar{2}111 sequence underscores the unusual stability of 6/136/13 in the SLL.

In this work, we consider the next member of the n¯2¯111\bar{n}\bar{2}111 parton sequence, namely 4¯2¯111\bar{4}\bar{2}111 and consider its feasibility at ν=2+4/9\nu=2+4/9. While FQHE at 2+4/92+4/9 has not been established conclusively, indications for it in the form of a minimum in the longitudinal resistance have already been seen in experiments Choi et al. (2008); Reichl et al. (2014). Future experiments on high-quality samples could likely establish FQHE at this filling by an observation of a well-quantized plateau in the Hall resistance. We show that the 4¯2¯111\bar{4}\bar{2}111 parton state gives a good description of the exact SLL Coulomb ground state seen in numerics. Furthermore, the parton state is energetically favorable compared to the 411411 Jain CF state at 2+4/92+4/9 in the thermodynamic limit. Therefore, if FQHE is established at ν=2+4/9\nu=2+4/9, it is highly likely to be distinct from its LLL counterpart at ν=4/9\nu=4/9, which is well described by the 411411 Jain CF state. The parton state is topologically distinct from the Jain CF state and we make predictions for experimentally measurable quantities that can unambiguously distinguish the two. In particular, the parton state supports counter-propagating edge modes that do not occur in the Jain CF state.

The parton theory put forth by Jain Jain (1989b) constructs FQHE states as a product of IQH states. The essential idea is to break each electron into fictitious objects called partons, place the partons into incompressible IQH states, and recover the final state by fusing the partons back into the physical electrons. The NN-electron wave function for the ll-parton Jain state, denoted as “n1n2nln_{1}n_{2}\cdots n_{l}”, is given by

Ψνn1n2nl=𝒫LLLμ=1lΦnμ({zj}).\Psi^{n_{1}n_{2}\cdots n_{l}}_{\nu}=\mathcal{P}_{\rm LLL}\prod_{\mu=1}^{l}\Phi_{n_{\mu}}(\{z_{j}\}). (1)

Here zj=xjiyjz_{j}=x_{j}-iy_{j} is the two-dimensional coordinate of electron jj which is parametrized as a complex number, μ\mu denotes the parton species and 𝒫LLL\mathcal{P}_{\rm LLL} implements projection into the LLL. Each parton species is exposed to the external magnetic field and occupies the same area, which fixes the charge of the μ\mu parton species to qμ=eν/nμq_{\mu}=-e\nu/n_{\mu}, with μqμ=e\sum_{\mu}q_{\mu}=-e, where e-e is the charge of the electron. The state Φnμ\Phi_{n_{\mu}} is the Slater determinant IQH wave function of NN electrons filling the lowest |nμ||n_{\mu}| LLs. We allow for nμ<0n_{\mu}<0 and negative values are denoted by n¯\bar{n}, with Φn¯=Φ|n|=[Φn]\Phi_{\bar{n}}=\Phi_{-|n|}=[\Phi_{n}]^{*}. The parton state of Eq. (1) occurs at a filling factor of ν=(μnμ1)1\nu=(\sum_{\mu}n_{\mu}^{-1})^{-1} and has a shift Wen and Zee (1992) of 𝒮=μnμ\mathcal{S}=\sum_{\mu}n_{\mu} in the spherical geometry.

The Laughlin state Laughlin (1983) is a “1111\cdots” parton state. The parton state n11n11 (n¯11\bar{n}11) with wave function Ψn/(2n+1)CF=𝒫LLLΦnΦ12\Psi^{\rm CF}_{n/(2n+1)}=\mathcal{P}_{\rm LLL}\Phi_{n}\Phi^{2}_{1} (Ψn/(2n1)CF=𝒫LLL[Φn]Φ12\Psi^{\rm CF}_{n/(2n-1)}=\mathcal{P}_{\rm LLL}[\Phi_{n}]^{*}\Phi^{2}_{1}) correspond to the Jain CF states Jain (1989a). Recently, it has been shown that parton states of the form “221,221\cdots,” which are not composite fermion states, could be viable candidates to describe certain FQH states observed in the LLL in wide quantum wells Faugno et al. (2019) and in LLs of graphene Wu et al. (2017); Kim et al. (2019).

The motivation for considering the 4¯2¯111\bar{4}\bar{2}111 parton state stems from the recent application of parton theory to capture states in the SLL Balram et al. (2018a, b, 2019). Consider the family of parton states described by the wave function

Ψν=2n/(5n2)n¯2¯111=𝒫LLLΦn¯Φ2¯Φ13Ψn/(2n1)CFΨ2/3CFΦ1.\Psi^{\bar{n}\bar{2}111}_{\nu=2n/(5n-2)}=\mathcal{P}_{\rm LLL}\Phi_{\bar{n}}\Phi_{\bar{2}}\Phi_{1}^{3}\sim\frac{\Psi^{\rm CF}_{n/(2n-1)}\Psi^{\rm CF}_{2/3}}{\Phi_{1}}. (2)

The \sim sign in Eq. (2) indicates that the states written at both sides of the sign differ slightly in the details of the projection. We expect that such details do not change the topological nature of the states Balram and Jain (2016). Throughout this text, for the n¯2¯111\bar{n}\bar{2}111 parton states we use the form given on the rightmost side of Eq. (2). A nice feature of the parton wave functions stated in Eq. (2) is that they can be evaluated for large system sizes which allow a reliable extrapolation of their thermodynamic energies. One can construct the above parton states for large system sizes because the constituent Jain CF states can be evaluated for hundreds of electrons using the Jain-Kamilla method of projection Jain and Kamilla (1997); Möller and Simon (2005); Jain (2007); Davenport and Simon (2012); Balram et al. (2015a). The Jain CF states in this work are evaluated using the Jain-Kamilla method.

The n=1n=1 member, namely 1¯2¯111\bar{1}\bar{2}111, is likely topologically equivalent to the ν=2/3\nu=2/3 2¯11\bar{2}11 Jain CF state Balram and Jain (2016). The n=2n=2 member 2¯2¯111\bar{2}\bar{2}111 has a good overlap with the exact SLL Coulomb ground state at ν=5/2\nu=5/2 Balram et al. (2018a). The n=3n=3 state, 3¯2¯111\bar{3}\bar{2}111, gives a good description of the Coulomb ground state at ν=2+6/13\nu=2+6/13 Balram et al. (2018b). In the Supplemental Material (SM) SM we provide further evidence in favor of the feasibility of the 3¯2¯111\bar{3}\bar{2}111 parton state to describe the ν=2+6/13\nu=2+6/13 FQHE. We shall consider the n=4n=4 member of this sequence which occurs at filling factor 4/94/9.

Although there is no definitive observation of FQHE at 4/94/9 in the SLL, signatures of incompressibility have been seen at ν=2+4/9\nu=2+4/9 and its particle-hole conjugate at ν=2+5/9\nu=2+5/9 Choi et al. (2008); Reichl et al. (2014). FQHE at ν=2+4/9\nu=2+4/9 is likely swamped by a bubble phase Shingla et al. (2018); however, it is likely that with improvements in the sample quality or for some interaction parameters close to that of the SLL Coulomb one, FQHE will ultimately be observed at 2+4/92+4/9.

For all our calculations, we deploy Haldane’s spherical geometry Haldane (1983), in which NN electrons reside on the surface of a sphere in the presence of a radial magnetic field generated by a monopole of strength 2Q(hc/e)2Q(hc/e) located at the center of the sphere. FQHE ground states occur at flux values 2Q=ν1N𝒮2Q=\nu^{-1}N-\mathcal{S}, where 𝒮\mathcal{S} is a rational number called the shift, which is useful in characterizing the topological nature of the FQHE state Wen and Zee (1992). All FQHE ground states are uniform on the sphere and thus have total orbital angular momentum L=0L=0. The parton states Ψν=2n/(5n2)n¯2¯111\Psi^{\bar{n}\bar{2}111}_{\nu=2n/(5n-2)} of Eq. (2) satisfy the flux-particle relationship 2Q=[(5n2)/(2n)]N(1n)2Q=[(5n-2)/(2n)]N-(1-n); i.e., their filling factors are ν=2n/(5n2)\nu=2n/(5n-2) and their shifts are 𝒮=1n\mathcal{S}=1-n. Of particular interest to us in this work is the 4¯2¯111\bar{4}\bar{2}111 parton state which has a shift of 𝒮=3\mathcal{S}=-3. This parton state is topologically distinct from the 411411 Jain CF state which also occurs at ν=4/9\nu=4/9 but has a shift of 𝒮=6\mathcal{S}=6. We assume a single-component system and neglect the effects of LL mixing and disorder. Under these assumptions, states related by particle-hole conjugation are considered on the same footing.

Throughout this work, we shall write wave functions in the LLL, which is where they are easily evaluated, even though they might apply to states occurring in the SLL. Haldane Haldane (1983) showed that the physics of the SLL can be simulated in the LLL by using an effective interaction that has the same set of Haldane pseudopotentials in the LLL as the Coulomb interaction has in the SLL. In this work, we have used the form of the effective interaction described in Ref. Shi et al. (2008) to simulate the physics of the SLL in the LLL.

Let us begin by testing the viability of the 4¯2¯111\bar{4}\bar{2}111 parton state for ν=2+4/9\nu=2+4/9 FQHE. In Fig. 1 we compare the energies of the 4¯2¯111\bar{4}\bar{2}111 parton and the 411411 Jain CF states at ν=4/9\nu=4/9 in the LLL and the SLL. In the LLL, as expected, we find that the Jain CF state has lower energy than the parton state. However, in the SLL we find that the 4¯2¯111\bar{4}\bar{2}111 parton state is energetically more favorable compared to the Jain CF state. For the sake of completeness, we have also investigated the competition between the parton and Jain CF states in the n=1n=1 LL of monolayer graphene. The effective interaction we use to simulate the physics of the n=1n=1 LL of monolayer graphene in the LLL is described in Ref. Balram et al. (2015b). We find the 411411 Jain CF state has lower energy here, consistent with the fact that experimentally observed FQHE states in the n=1n=1 LL of monolayer graphene are well described by the CF paradigm Amet et al. (2015); Balram et al. (2015b); Zeng et al. (2019). Results for n=0n=0 LL of graphene are identical to those in the LLL of GaAs under our working assumptions of neglecting effects of finite width and LL mixing.

Refer to caption
Figure 1: (color online) Thermodynamic extrapolations of the Coulomb per-particle energies for the 411411 Jain composite fermion (CF) state and the 4¯2¯111\bar{4}\bar{2}111 parton state. The left-hand, middle, and right-hand panels show energies for ν=4/9\nu=4/9 in the n=0n=0 LL, n=1n=1 LL of GaAs, and in the n=1n=1 LL of monolayer graphene respectively. The extrapolated energies, obtained from a linear fit in 1/N1/N, are quoted in units of e2/(ϵ)e^{2}/(\epsilon\ell) on the plot (error in the fit is indicated in the parentheses). The energies include contributions of electron-background and background-background interactions and are density corrected Morf et al. (1986). The LLL Coulomb energy for the ν=4/9\nu=4/9 Jain CF state has been reproduced from Ref. Balram and Jain (2017).

Next, we turn to comparisons of the parton state with the exact SLL Coulomb ground state. The smallest system accessible to exact diagonalization (ED) is that of N=16N=16 electrons at a flux of 2Q=392Q=39 which has a Hilbert space dimension of 7×1087\times 10^{8}. We have evaluated the ground state for this system with the truncated pseudopotentials from the disk geometry, which differ slightly from the spherical pseudopotentials but are known to provide a more reliable extrapolation to the thermodynamic limit Peterson et al. (2008a, b). The exact SLL Coulomb ground obtained by using the truncated disk pseudopotentials has L=0L=0. In Fig. 2 we compare the pair-correlation function foo of this exact SLL Coulomb ground state with that of the 4¯2¯111\bar{4}\bar{2}111 parton state. Both these pair-correlation functions show oscillations that decay at long distances, which is a typical characteristic of incompressible states Kamilla et al. (1997); Balram et al. (2015c). Moreover, the two pair-correlation functions are in reasonable agreement with each other. For completeness, we have also evaluated the exact LLL Coulomb ground state for the same system. The overlaps of the LLL and SLL Coulomb ground states obtained using the disk pseudopotentials is 0.36630.3663 and their pair-correlation functions are also very different from each other SM which indicate that the nature of the ground state in the two LLs are very different.

Refer to caption
Figure 2: (color online) The pair correlation function g(r)g(r) as a function of the arc distance for the exact second Landau level Coulomb ground state obtained using the disk pseudopotentials for width w=0w=0 and w=3w=3\ell, and the 4¯2¯111\bar{4}\bar{2}111 parton state of Eq. (2) for N=16N=16 electrons at a flux of 2Q=392Q=39.

Currently, we do not have a reliable estimate of the thermodynamic values of the gaps predicted by our parton ansatz. However, we can extract the charge and neutral gaps for N=16N=16 particles from exact diagonalization of the SLL Coulomb interaction at the parton shift. The charge gap here is defined as [(2Q=40)+(2Q=38)2(2Q=39)]/4[\mathcal{E}(2Q=40)+\mathcal{E}(2Q=38)-2\mathcal{E}(2Q=39)]/4, where (2Q)\mathcal{E}(2Q) is the exact ground state energy at flux 2Q2Q and the factor of 44 in the denominator accounts for the fact that the addition of a single flux quantum in the parton state produces four fundamental quasiholes. The neutral gap is defined as the difference between the two lowest exact energies at the flux of 2Q=392Q=39. The charge and neutral gaps for N=16N=16, evaluated using exact diagonalization with the disk pseudopotentials, are 0.009e2/ϵ0.009~{}e^{2}/\epsilon\ell and 0.005e2/ϵ0.005~{}e^{2}/\epsilon\ell respectively, where =c/(eB)\ell=\sqrt{\hbar c/(eB)} is the magnetic length and ϵ\epsilon is the dielectric constant of the background host material.

We next consider the effect of finite width on the system, which we model by taking the transverse wave function to be the ground state for an infinite square quantum well of width ww (see SM SM for details). For the disk pseudopotentials, we find that the ground state for 1616 electrons has L=0L=0 for (at least) w5w\leq 5\ell SM . Moreover, the pair-correlation function of the exact ground state agrees well with that of the parton state for the entire range of widths considered in this work (see Fig. 2 and Ref. SM ). Furthermore, we find that the system has robust charge and neutral gaps for all the widths considered. We note that the L=0L=0 ground state is delicate. In particular, the exact SLL Coulomb ground state obtained using the spherical pseudopotentials has L=2L=2. However, the overlap between the lowest energy L=0L=0 state obtained using the spherical pseudopotentials and the ground state obtained using the disk pseudopotentials is 0.96920.9692, which indicates that these two states are close to each other. Encouragingly, with the spherical pseudopotentials, as the quantum well width is increased the ground state turns uniform in the range w(0.5,1]w\in(0.5,1]\ell and stays uniform for w[1,10]w\in[1,10]\ell SM . These results indicate that finite thickness enhances the stability of the parton state.

Now that we have made a case for the plausibility of the 4¯2¯111\bar{4}\bar{2}111 parton state to occur in the SLL, we shall turn to deduce the experimental consequences of this parton ansatz. An additional particle in the factor Φ4¯\Phi_{\bar{4}} has charge e/9e/9, whereas that in the factors Φ2¯\Phi_{\bar{2}} and Φ1\Phi_{1} has a charge 2e/92e/9 and 4e/9-4e/9 respectively. All the quasiparticles of the 4¯2¯111\bar{4}\bar{2}111 parton state obey Abelian braid statistics Wen (1991). The 4/94/9 Jain CF state is also an Abelian state and hosts quasiparticles of charge e/9-e/9 and 4e/9-4e/9.

Next, to infer other topological consequences of the 4¯2¯111\bar{4}\bar{2}111 ansatz, we consider the low-energy effective theory of its edge, which is described by the Lagrangian density WEN (1991); Wen (1992, 1995); Moore and Wen (1998):

=14πKIJϵαβγaαIβaγJ12πϵαβγtIAαβaγI.\mathcal{L}=-\frac{1}{4\pi}K_{\rm IJ}\epsilon^{\alpha\beta\gamma}a_{\alpha}^{\rm I}\partial_{\beta}a_{\gamma}^{\rm J}-\frac{1}{2\pi}\epsilon^{\alpha\beta\gamma}t_{\rm I}A_{\alpha}\partial_{\beta}a_{\gamma}^{\rm I}. (3)

Here ϵαβγ\epsilon^{\alpha\beta\gamma} is the fully anti-symmetric Levi-Civita tensor, AA is the vector potential corresponding to the external electromagnetic field, aa is the internal gauge field, and we have used the Einstein’s convention of summing repeated indices. The integer-valued symmetric KK matrix and the charge vector tt of Eq. (3) for the parton state are given by (see SM SM for a derivation)

K=(2110112101112010002111111),t=( 0 0 0 0 1).K=\begin{pmatrix}-2&-1&-1&0&1\\ -1&-2&-1&0&1\\ -1&-1&-2&0&1\\ 0&0&0&-2&1\\ 1&1&1&1&1\\ \end{pmatrix},\quad t=\begin{pmatrix}\,0\,\\ \,0\,\\ \,0\,\\ \,0\,\\ \,1\,\\ \end{pmatrix}. (4)

The above KK matrix has four negative and one positive eigenvalues and thus the 4¯2¯111\bar{4}\bar{2}111 state hosts four upstream and one downstream edge modes. A naive counting suggests that there are a total of nine edge states for the 4¯2¯111\bar{4}\bar{2}111 ansatz: four from the factor Φ4¯\Phi_{\bar{4}}, two from Φ2¯\Phi_{\bar{2}}, and one from each factor of Φ1\Phi_{1}. However, these edges states are not all independent since the density variations of the five partons must be identified. This results in four constraints and leads to five edge states consistent with that ascertained from the above KK matrix.

Assuming full equilibration of the edge modes, the thermal Hall conductance κxy\kappa_{xy} at temperatures much smaller than the gap takes a quantized value proportional to the chiral central charge cc_{-}, which is defined as the difference in the number of downstream and upstream modes: κxy=c[π2kB2/(3h)]T\kappa_{xy}=c_{-}[\pi^{2}k_{\rm B}^{2}/(3h)]T Kane and Fisher (1997). For 4¯2¯111\bar{4}\bar{2}111 ansatz, we thus predict a thermal Hall conductance of κxy=3[π2kB2/(3h)]T\kappa_{xy}=-3[\pi^{2}k_{\rm B}^{2}/(3h)]T. The Hall viscosity of the 4¯2¯111\bar{4}\bar{2}111 state is also expected to be quantized Read (2009): ηH=ρ0𝒮/4\eta_{\rm H}=\hbar\rho_{0}\mathcal{S}/4, where ρ0=(4/9)/(2π2)\rho_{0}=(4/9)/(2\pi\ell^{2}) is the electron density and 𝒮=3\mathcal{S}=-3 is the shift of the parton state. The ground-state degeneracy of the parton state on a topologically nontrivial manifold with genus gg is |Det(K)|g=18g|{\rm Det}(K)|^{g}=18^{g}. Besides the nn¯111n\bar{n}111 parton states Balram et al. (2020), the 4¯2¯111\bar{4}\bar{2}111 ansatz provides another example of a fully spin polarized Abelian FQH state at ν=a/b\nu=a/b (with a,ba,b coprime), which has a ground state degeneracy on the torus that is greater than bb.

The 4/94/9 Jain CF state is described by the 4×44\times 4 KK matrix K=2C4+𝕀4K=2C_{4}+\mathbb{I}_{4}, where CkC_{k} is the k×kk\times k matrix of all ones and 𝕀k\mathbb{I}_{k} is the k×kk\times k identity matrix, and charge vector t=(1,1,1,1)Tt=(1,1,1,1)^{T}. In contrast to the 4¯2¯111\bar{4}\bar{2}111 state, assuming the absence of edge reconstruction, the 4/94/9 Jain CF state has four downstream edge states and no upstream modes. The 4/94/9 Jain CF state thus has a thermal Hall conductance of κxy=4[π2kB2/(3h)]T\kappa_{xy}=4[\pi^{2}k_{\rm B}^{2}/(3h)]T. Moreover, the Hall viscosity of the 4/94/9 Jain CF state is given by ηH=(3/2)ρ0\eta_{\rm H}=(3/2)\hbar\rho_{0}, corresponding to shift 𝒮=6{\cal S}=6. On a manifold of genus gg, the 4/94/9 Jain CF state has a degeneracy of 9g9^{g}.

The presence of upstream neutral modes can be detected in shot noise experiments Bid et al. (2010); Dolev et al. (2011); Gross et al. (2012); Inoue et al. (2014). Recently, thermal Hall measurements have been carried out at several filling factors in the lowest as well as the second LL Banerjee et al. (2017, 2018); Srivastav et al. (2019). These experiments can be used to test the predictions of the parton theory and therefore can unambiguously distinguish between the topological nature of the 4/9 states in the SLL and the LLL. In particular, including the contributions of the filled LLLs of spin up and spin down, the thermal Hall conductance of the 4¯2¯111\bar{4}\bar{2}111 state in the SLL is [π2kB2/(3h)]T-[\pi^{2}k_{\rm B}^{2}/(3h)]T which is different from what one would expect from the 4/94/9 Jain CF state in the SLL, which has κxy=6[π2kB2/(3h)]T\kappa_{xy}=6[\pi^{2}k_{\rm B}^{2}/(3h)]T.

In summary, we have considered the viability of the “4¯2¯111\bar{4}\bar{2}111” parton state for FQHE at ν=2+4/9\nu=2+4/9 where the first signs of incompressibility in the form of minimum in longitudinal resistance have already been observed experimentally Choi et al. (2008); Reichl et al. (2014). Interestingly, if FQHE eventually stabilizes at this filling factor, then it is likely to be topologically different from its LLL counterpart at ν=4/9\nu=4/9, which is described by a Jain CF state. We also proposed experimental measurements that can reveal the underlying topological structure of the parton state and decisively distinguish it from the 4/94/9 state occurring in the lowest Landau level.

Acknowledgements.
We acknowledge useful discussions with Maissam Barkeshli, Jainendra K. Jain, Mark S. Rudner, Dam T. Son, and Bo Yang. This work was supported by the Polish NCN Grant No. 2014/14/A/ST3/00654 (A. W.). We thank Wrocław Centre for Networking and Supercomputing and Academic Computer Centre CYFRONET, both parts of PL-Grid Infrastructure. Some of the numerical calculations reported in this work were carried out on the Nandadevi supercomputer, which is maintained and supported by the Institute of Mathematical Science’s High-Performance Computing Center. This research was supported in part by the International Centre for Theoretical Sciences (ICTS) during a visit for the program Novel Phases of Quantum Matter (code: ICTS/topmatter2019/12)

References

  • Tsui et al. (1982) D. C. Tsui, H. L. Stormer, and A. C. Gossard, Phys. Rev. Lett. 48, 1559 (1982).
  • Laughlin (1983) R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983).
  • Kitaev (2003) A. Kitaev, Annals of Physics 303, 2 (2003).
  • Nayak et al. (2008) C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008).
  • Eisenstein et al. (1990) J. P. Eisenstein, H. L. Stormer, L. N. Pfeiffer, and K. W. West, Phys. Rev. B 41, 7910 (1990).
  • Jain (1989a) J. K. Jain, Phys. Rev. Lett. 63, 199 (1989a).
  • Jain (2007) J. K. Jain, Composite Fermions (Cambridge University Press, New York, US, 2007).
  • Willett et al. (1987) R. Willett, J. P. Eisenstein, H. L. Störmer, D. C. Tsui, A. C. Gossard, and J. H. English, Phys. Rev. Lett. 59, 1776 (1987).
  • Moore and Read (1991) G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991).
  • Read and Green (2000) N. Read and D. Green, Phys. Rev. B 61, 10267 (2000).
  • d’Ambrumenil and Reynolds (1988) N. d’Ambrumenil and A. M. Reynolds, Journal of Physics C: Solid State Physics 21, 119 (1988).
  • Balram et al. (2013) A. C. Balram, Y.-H. Wu, G. J. Sreejith, A. Wójs, and J. K. Jain, Phys. Rev. Lett. 110, 186801 (2013).
  • Johri et al. (2014) S. Johri, Z. Papić, R. N. Bhatt, and P. Schmitteckert, Phys. Rev. B 89, 115124 (2014).
  • Peterson et al. (2015) M. R. Peterson, Y.-L. Wu, M. Cheng, M. Barkeshli, Z. Wang, and S. Das Sarma, Phys. Rev. B 92, 035103 (2015).
  • Kleinbaum et al. (2015) E. Kleinbaum, A. Kumar, L. N. Pfeiffer, K. W. West, and G. A. Csáthy, Phys. Rev. Lett. 114, 076801 (2015).
  • Jeong et al. (2017) J.-S. Jeong, H. Lu, K. H. Lee, K. Hashimoto, S. B. Chung, and K. Park, Phys. Rev. B 96, 125148 (2017).
  • Balram et al. (2020) A. C. Balram, J. K. Jain, and M. Barkeshli, Phys. Rev. Research 2, 013349 (2020).
  • Xia et al. (2004) J. S. Xia, W. Pan, C. L. Vicente, E. D. Adams, N. S. Sullivan, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer, K. W. Baldwin, and K. W. West, Phys. Rev. Lett. 93, 176809 (2004).
  • Choi et al. (2008) H. C. Choi, W. Kang, S. Das Sarma, L. N. Pfeiffer, and K. W. West, Phys. Rev. B 77, 081301 (2008).
  • Pan et al. (2008) W. Pan, J. S. Xia, H. L. Stormer, D. C. Tsui, C. Vicente, E. D. Adams, N. S. Sullivan, L. N. Pfeiffer, K. W. Baldwin, and K. W. West, Phys. Rev. B 77, 075307 (2008).
  • Kumar et al. (2010) A. Kumar, G. A. Csáthy, M. J. Manfra, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 105, 246808 (2010).
  • Read and Rezayi (1999) N. Read and E. Rezayi, Phys. Rev. B 59, 8084 (1999).
  • Rezayi and Read (2009) E. H. Rezayi and N. Read, Phys. Rev. B 79, 075306 (2009).
  • Wójs (2009) A. Wójs, Phys. Rev. B 80, 041104 (2009).
  • Bonderson et al. (2012) P. Bonderson, A. E. Feiguin, G. Möller, and J. K. Slingerland, Phys. Rev. Lett. 108, 036806 (2012).
  • Sreejith et al. (2013) G. J. Sreejith, Y.-H. Wu, A. Wójs, and J. K. Jain, Phys. Rev. B 87, 245125 (2013).
  • Zhu et al. (2015) W. Zhu, S. S. Gong, F. D. M. Haldane, and D. N. Sheng, Phys. Rev. Lett. 115, 126805 (2015).
  • Mong et al. (2017) R. S. K. Mong, M. P. Zaletel, F. Pollmann, and Z. Papić, Phys. Rev. B 95, 115136 (2017).
  • Pakrouski et al. (2016) K. Pakrouski, M. Troyer, Y.-L. Wu, S. Das Sarma, and M. R. Peterson, Phys. Rev. B 94, 075108 (2016).
  • Shingla et al. (2018) V. Shingla, E. Kleinbaum, A. Kumar, L. N. Pfeiffer, K. W. West, and G. A. Csáthy, Phys. Rev. B 97, 241105 (2018).
  • Reichl et al. (2014) C. Reichl, J. Chen, S. Baer, C. Rössler, T. Ihn, K. Ensslin, W. Dietsche, and W. Wegscheider, New Journal of Physics 16, 023014 (2014).
  • Faugno et al. (2020) W. N. Faugno, J. K. Jain, and A. C. Balram, Phys. Rev. Research 2, 033223 (2020).
  • Balram et al. (2018a) A. C. Balram, M. Barkeshli, and M. S. Rudner, Phys. Rev. B 98, 035127 (2018a).
  • Balram et al. (2018b) A. C. Balram, S. Mukherjee, K. Park, M. Barkeshli, M. S. Rudner, and J. K. Jain, Phys. Rev. Lett. 121, 186601 (2018b).
  • Balram et al. (2019) A. C. Balram, M. Barkeshli, and M. S. Rudner, Phys. Rev. B 99, 241108 (2019).
  • Jain (1989b) J. K. Jain, Phys. Rev. B 40, 8079 (1989b).
  • Wen and Zee (1992) X. G. Wen and A. Zee, Phys. Rev. Lett. 69, 953 (1992).
  • Faugno et al. (2019) W. N. Faugno, A. C. Balram, M. Barkeshli, and J. K. Jain, Phys. Rev. Lett. 123, 016802 (2019).
  • Wu et al. (2017) Y. Wu, T. Shi, and J. K. Jain, Nano Letters 17, 4643 (2017), pMID: 28649831, eprint http://dx.doi.org/10.1021/acs.nanolett.7b01080.
  • Kim et al. (2019) Y. Kim, A. C. Balram, T. Taniguchi, K. Watanabe, J. K. Jain, and J. H. Smet, Nature Physics 15, 154 (2019).
  • Balram and Jain (2016) A. C. Balram and J. K. Jain, Phys. Rev. B 93, 235152 (2016).
  • Jain and Kamilla (1997) J. K. Jain and R. K. Kamilla, Phys. Rev. B 55, R4895 (1997).
  • Möller and Simon (2005) G. Möller and S. H. Simon, Phys. Rev. B 72, 045344 (2005).
  • Davenport and Simon (2012) S. C. Davenport and S. H. Simon, Phys. Rev. B 85, 245303 (2012).
  • Balram et al. (2015a) A. C. Balram, C. Töke, A. Wójs, and J. K. Jain, Phys. Rev. B 92, 075410 (2015a).
  • (46) See Supplemental Material (SM) accompanying this paper for (i) a derivation of the low-energy effective edge theory of the 4¯2¯111\bar{4}\bar{2}111 parton state and (ii) a detailed study of the effect of finite width at ν=2+4/9\nu=2+4/9 and (iii) a comparison of candidate states with exact Coulomb ground states at other filling factors, namely ν=7/3,5/2\nu=7/3,~{}5/2 and 2+6/132+6/13, in the second Landau level. The SM includes Refs. Morf (1998); Wójs (2009); Hutasoit et al. (2017); Zuo et al. (2020).
  • Haldane (1983) F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983).
  • Shi et al. (2008) C. Shi, S. Jolad, N. Regnault, and J. K. Jain, Phys. Rev. B 77, 155127 (2008).
  • Balram et al. (2015b) A. C. Balram, C. Tőke, A. Wójs, and J. K. Jain, Phys. Rev. B 92, 205120 (2015b).
  • Amet et al. (2015) F. Amet, A. J. Bestwick, J. R. Williams, L. Balicas, K. Watanabe, T. Taniguchi, and D. Goldhaber-Gordon, Nat. Commun. 6, 5838 (2015).
  • Zeng et al. (2019) Y. Zeng, J. I. A. Li, S. A. Dietrich, O. M. Ghosh, K. Watanabe, T. Taniguchi, J. Hone, and C. R. Dean, Phys. Rev. Lett. 122, 137701 (2019).
  • Morf et al. (1986) R. Morf, N. d’Ambrumenil, and B. I. Halperin, Phys. Rev. B 34, 3037 (1986).
  • Balram and Jain (2017) A. C. Balram and J. K. Jain, Phys. Rev. B 96, 235102 (2017).
  • Peterson et al. (2008a) M. R. Peterson, T. Jolicoeur, and S. Das Sarma, Phys. Rev. Lett. 101, 016807 (2008a).
  • Peterson et al. (2008b) M. R. Peterson, T. Jolicoeur, and S. Das Sarma, Phys. Rev. B 78, 155308 (2008b).
  • (56) We compare the pair-correlation functions since, due to technical reasons, it has not yet been possible to calculate the overlap of the parton state with the state obtained from exact diagonalization for the large system size considered in this work.
  • Kamilla et al. (1997) R. K. Kamilla, J. K. Jain, and S. M. Girvin, Phys. Rev. B 56, 12411 (1997).
  • Balram et al. (2015c) A. C. Balram, C. Tőke, and J. K. Jain, Phys. Rev. Lett. 115, 186805 (2015c).
  • Wen (1991) X. G. Wen, Phys. Rev. Lett. 66, 802 (1991).
  • WEN (1991) X. WEN, Modern Physics Letters B 05, 39 (1991).
  • Wen (1992) X.-G. Wen, International Journal of Modern Physics B 06, 1711 (1992).
  • Wen (1995) X.-G. Wen, Advances in Physics 44, 405 (1995).
  • Moore and Wen (1998) J. E. Moore and X.-G. Wen, Phys. Rev. B 57, 10138 (1998).
  • Kane and Fisher (1997) C. L. Kane and M. P. A. Fisher, Phys. Rev. B 55, 15832 (1997).
  • Read (2009) N. Read, Phys. Rev. B 79, 045308 (2009).
  • Bid et al. (2010) A. Bid, N. Ofek, H. Inoue, M. Heiblum, C. L. Kane, V. Umansky, and D. Mahalu, Nature 466, 585 (2010).
  • Dolev et al. (2011) M. Dolev, Y. Gross, R. Sabo, I. Gurman, M. Heiblum, V. Umansky, and D. Mahalu, Phys. Rev. Lett. 107, 036805 (2011).
  • Gross et al. (2012) Y. Gross, M. Dolev, M. Heiblum, V. Umansky, and D. Mahalu, Phys. Rev. Lett. 108, 226801 (2012).
  • Inoue et al. (2014) H. Inoue, A. Grivnin, Y. Ronen, M. Heiblum, V. Umansky, and D. Mahalu, Nature Communications 5, 4067 (2014).
  • Banerjee et al. (2017) M. Banerjee, M. Heiblum, A. Rosenblatt, Y. Oreg, D. E. Feldman, A. Stern, and V. Umansky, Nature 545, 75+ (2017).
  • Banerjee et al. (2018) M. Banerjee, M. Heiblum, V. Umansky, D. E. Feldman, Y. Oreg, and A. Stern, Nature 559, 205 (2018).
  • Srivastav et al. (2019) S. K. Srivastav, M. R. Sahu, K. Watanabe, T. Taniguchi, S. Banerjee, and A. Das, Science Advances 5 (2019), eprint https://advances.sciencemag.org/content/5/7/eaaw5798.full.pdf.
  • Morf (1998) R. H. Morf, Phys. Rev. Lett. 80, 1505 (1998).
  • Hutasoit et al. (2017) J. A. Hutasoit, A. C. Balram, S. Mukherjee, Y.-H. Wu, S. S. Mandal, A. Wójs, V. Cheianov, and J. K. Jain, Phys. Rev. B 95, 125302 (2017).
  • Kuśmierz and Wójs (2018) B. Kuśmierz and A. Wójs, Phys. Rev. B 97, 245125 (2018).
  • Scarola et al. (2002) V. W. Scarola, S.-Y. Lee, and J. K. Jain, Phys. Rev. B 66, 155320 (2002).
  • Pakrouski et al. (2015) K. Pakrouski, M. R. Peterson, T. Jolicoeur, V. W. Scarola, C. Nayak, and M. Troyer, Phys. Rev. X 5, 021004 (2015).
  • Zuo et al. (2020) Z.-W. Zuo, A. C. Balram, S. Pu, J. Zhao, A. Wójs, and J. K. Jain, arXiv e-prints arXiv:2002.09356 (2020), eprint 2002.09356.

Supplemental Material for “Fractional Quantum Hall Effect at ν=2+4/9\nu=2+4/9

In this Supplemental Material (SM) we provide

  • a derivation of the low-energy effective edge theory of the 4¯2¯111\bar{4}\bar{2}111 parton state (Sec. I),

  • a detailed study of the effect of finite width at ν=2+4/9\nu=2+4/9 (Sec. II), and

  • a comparison of the Laughlin, Pfaffian and the 3¯2¯111\bar{3}\bar{2}111 parton states with the exact Coulomb ground state at filling factors ν=7/3,5/2\nu=7/3,~{}5/2 and 2+6/132+6/13 respectively (Sec. III).

I Derivation of the low-energy effective theory of the parton 4¯2¯111\bar{4}\bar{2}111 ν=4/9\nu=4/9 edge

To derive the low-energy effective theory of the parton 4/94/9 edge we closely follow the procedure outlined in the Supplemental Material of Ref. Balram et al. (2018b). We are interested in the n=4n=4 member of the n¯2¯111\bar{n}\bar{2}111 parton sequence. We shall first discuss the case of general nn and then specialize to the case of n=4n=4. The unprojected wave function of the n¯2¯111\bar{n}\bar{2}111 parton state can be re-written as

Ψ2n/(5n2)n¯2¯111=[Φn][Φ2]Ψ1/3L,\Psi^{\bar{n}\bar{2}111}_{2n/(5n-2)}=[\Phi_{n}]^{*}[\Phi_{2}]^{*}\Psi^{L}_{1/3}, (S1)

where Φn\Phi_{n} is the Slater determinant state of nn filled Landau levels and Ψ1/3LΦ13\Psi^{L}_{1/3}\equiv\Phi^{3}_{1} is the ν=1/3\nu=1/3 Laughlin state Laughlin (1983). This state can be expressed in terms of partons =f1f2f3\wp=f_{1}f_{2}f_{3}, where the fif_{i}’s are fermionic partons in the following mean-field states:

  • f1f_{1} is in a ν=n\nu=-n integer quantum Hall (IQH) state

  • f2f_{2} is in a ν=2\nu=-2 IQH state

  • f3f_{3} is in a ν=1/3\nu=1/3 Laughlin state.

The charges of these partons are Jain (1989b) q1=2e/(5n2)q_{1}=2e/(5n-2), q2=ne/(5n2)q_{2}=ne/(5n-2) and q3=6ne/(5n2)q_{3}=-6ne/(5n-2), where e-e is the electron charge. For n=2n=2, the parton state of Eq. (S1) describes a non-Abelian state at ν=1/2\nu=1/2 Balram et al. (2018a). From here on we restrict to n2n\neq 2, wherein the parton state of Eq. (S1) describes an Abelian state with a residual U(1)×U(1)U(1)\times U(1) gauge symmetry associated with the transformations:

f1eiθ1f1,f2eiθ1+iθ2f2,f3eiθ2f3.f_{1}\rightarrow e^{i\theta_{1}}f_{1},~{}f_{2}\rightarrow e^{-i\theta_{1}+i\theta_{2}}f_{2},~{}f_{3}\rightarrow e^{-i\theta_{2}}f_{3}. (S2)

Therefore, we have two internal emergent U(1)U(1) gauge fields, denoted by hμh_{\mu} and gμg_{\mu}, associated with the above transformations. The low-energy effective field theory for this parton mean-field state is described by the Lagrangian density (here and henceforth we have set e=1e=1 for convenience) Balram et al. (2018b):

\displaystyle\mathcal{L} =\displaystyle= 14πi=1nα(i)α(i)+12πi=1n(h+q1A)α(i)\displaystyle\frac{1}{4\pi}\sum_{i=1}^{n}\alpha^{(i)}\partial\alpha^{(i)}+\frac{1}{2\pi}\sum_{i=1}^{n}(h+q_{1}A)\partial\alpha^{(i)} (S3)
+\displaystyle+ 14πj=12β(j)β(j)+12πj=12(gh+q2A)β(j)\displaystyle\frac{1}{4\pi}\sum_{j=1}^{2}\beta^{(j)}\partial\beta^{(j)}+\frac{1}{2\pi}\sum_{j=1}^{2}(g-h+q_{2}A)\partial\beta^{(j)}
\displaystyle- 34πγγ+12π(g+q3A)γ,\displaystyle\frac{3}{4\pi}\gamma\partial\gamma+\frac{1}{2\pi}(-g+q_{3}A)\partial\gamma,

where AA is the external physical electromagnetic vector potential, and α(i)\alpha^{(i)}, β(j)\beta^{(j)} and γ\gamma are U(1)U(1) gauge fields describing the current fluctuations of the IQH and Laughlin states. Furthermore, we have used the short-hand notation ααϵμνλαμναλ\alpha\partial\alpha\equiv\epsilon^{\mu\nu\lambda}\alpha_{\mu}\partial_{\nu}\alpha_{\lambda}, where ϵμνλ\epsilon^{\mu\nu\lambda} is the fully anti-symmetric Levi-Civita tensor and repeated indices are summed over.

This Chern-Simons theory can be further simplified by integrating out the internal gauge fields hh and gg. Doing so we obtain the following two constraints Balram et al. (2018b)

i=1nα(i)=j=12β(j)+c,\sum_{i=1}^{n}\alpha^{(i)}=\sum_{j=1}^{2}\beta^{(j)}+c, (S4)

and

γ=j=12β(j)+d,\gamma=\sum_{j=1}^{2}\beta^{(j)}+d, (S5)

where cc and dd are U(1)U(1) gauge fields that satisfy ϵμνλνcλ=0\epsilon^{\mu\nu\lambda}\partial_{\nu}c_{\lambda}=0 and ϵμνλνdλ=0\epsilon^{\mu\nu\lambda}\partial_{\nu}d_{\lambda}=0. Note that when substituting Eq. (S4) and (S5) into Eq. (S3), all terms involving the gauge fields cc and dd vanish. Thus, we end up with a simplified U(1)n+1U(1)^{n+1} Chern-Simons theory which can be described by an integer valued symmetric (n+1)×(n+1)(n+1)\times(n+1) KK matrix Balram et al. (2018b).

Using the constraints of Eqs. (S4) and (S5) we can eliminate α(n)\alpha^{(n)} and β(2)\beta^{(2)} by noting that

α(n)=γα(1)α(2)α(n1)+cd,β(2)=γβ(1)d.\alpha^{(n)}=\gamma-\alpha^{(1)}-\alpha^{(2)}-\cdots-\alpha^{(n-1)}+c-d,~{}\beta^{(2)}=\gamma-\beta^{(1)}-d. (S6)

Substituting these back into the Lagrangian density given in Eq. (S3) and using the fact that all terms involving the gauge fields cc and dd vanish, we obtain the following simplified Lagrangian density:

\displaystyle\mathcal{L} =\displaystyle= 14π(α(1)α(1)+α(2)α(2)++α(n1)α(n1)+(γα(1)α(2)α(n1))(γα(1)α(2)α(n1)))\displaystyle\frac{1}{4\pi}\left(\alpha^{(1)}\partial\alpha^{(1)}+\alpha^{(2)}\partial\alpha^{(2)}+\cdots+\alpha^{(n-1)}\partial\alpha^{(n-1)}+(\gamma-\alpha^{(1)}-\alpha^{(2)}-\cdots-\alpha^{(n-1)})\partial(\gamma-\alpha^{(1)}-\alpha^{(2)}-\cdots-\alpha^{(n-1)})\right)
+\displaystyle+ 14π(β(1)β(1)+(γβ(1))(γβ(1)))\displaystyle\frac{1}{4\pi}\left(\beta^{(1)}\partial\beta^{(1)}+(\gamma-\beta^{(1)})\partial(\gamma-\beta^{(1)})\right)
\displaystyle- 34πγγ+12πAγ.\displaystyle\frac{3}{4\pi}\gamma\partial\gamma+\frac{1}{2\pi}A\partial\gamma.

By defining a new set of gauge fields:

(a1,a2,,an1,an,an+1)=(α(1),α(2),,α(n1),β(1),γ),(a^{1},a^{2},\cdots,a^{n-1},a^{n},a^{n+1})=(\alpha^{(1)},\alpha^{(2)},\cdots,\alpha^{(n-1)},\beta^{(1)},\gamma), (S8)

we can write the Lagrangian in the standard form WEN (1991); Wen (1992, 1995); Moore and Wen (1998):

=14πKIJaIaJ+12πtIAaI.\mathcal{L}=-\frac{1}{4\pi}K_{\rm IJ}a^{\rm I}\partial a^{\rm J}+\frac{1}{2\pi}t^{\rm I}A\partial a^{\rm I}. (S9)

Here the charge vector is t=(0,0,,0,1)Tt=(0,0,\cdots,0,1)^{\rm T} and the KK matrix is given by

K=(210101120100211111).K=\begin{pmatrix}-2&\cdots&-1&0&1\\ \vdots&\ddots&\cdots&0&1\\ -1&\vdots&-2&0&1\\ 0&\cdots&0&-2&1\\ 1&\cdots&1&1&1\\ \end{pmatrix}. (S10)

The KK-matrix can be specified succinctly as

Ki,j={2,1i=jn1,i=n+1j=n+11,1i<jn11,1j<in10,otherwsieK_{i,j}=\begin{cases}-2,&1\leq i=j\leq n\\ 1,&i=n+1\lor j=n+1\\ -1,&1\leq i<j\leq n-1\\ -1,&1\leq j<i\leq n-1\\ 0,&{\rm otherwsie}\end{cases} (S11)

The filling factor is given by Wen (1995):

ν=tTK1t=Kn+1,n+11=2n/(5n2),\nu=t^{\rm T}\cdot K^{-1}\cdot t=K^{-1}_{n+1,n+1}=2n/(5n-2), (S12)

as anticipated from the microscopic wave function. The ground-state degeneracy of the n¯2¯111\bar{n}\bar{2}111 parton state on a manifold with nontrivial genus gg is Wen (1995):

ground-state degeneracy=|Det(K)|g=(5n2)g.\text{ground-state degeneracy}=|{\rm Det}(K)|^{g}=(5n-2)^{g}. (S13)

An interesting point to note is that when n=2mn=2m is even, then the filling factor ν=2n/(5n2)=2m/(5m1)\nu=2n/(5n-2)=2m/(5m-1). For these set of fillings, the ground-state degeneracy is (5n2)g=2g(5m1)g(5n-2)^{g}=2^{g}(5m-1)^{g}. Therefore, for even n2n\neq 2 we get a single component Abelian state at ν=a/b\nu=a/b (with a,ba,b coprime) which has a ground-state degeneracy on the torus (=2b=2b) which is greater than bb (see Ref. Balram et al. (2020) for another example of such a state).

The KK matrix of Eq. (S10) has one positive and nn negative eigenvalues which indicates that the n¯2¯111\bar{n}\bar{2}111 state hosts one downstream and nn upstream edge modes which gives it a chiral central charge of 1n1-n. The particle-hole conjugate of the n¯2¯111\bar{n}\bar{2}111 state, which occurs at filling factor 12n/(5n2)=(3n2)/(5n2)1-2n/(5n-2)=(3n-2)/(5n-2), can be viewed as a fully chiral state with nn downstream modes which results in a chiral central charge of nn. The charges and the braiding statistics of the quasiparticles of the n¯2¯111\bar{n}\bar{2}111 state can be ascertained from its KK matrix following the work of Ref. Wen (1995).

I.1 Coupling to curvature and shift Wen and Zee (1992)

To compute the shift 𝒮\mathcal{S} Wen and Zee (1992) of the n¯2¯111\bar{n}\bar{2}111 state on the sphere from its topological field theory, we need to couple it to the curvature. For a state filling the nthn^{\rm th} Landau level and described by the U(1)U(1) gauge field Γ\Gamma, upon including the coupling to curvature, the effective Lagrangian density is modified by the addition of the following term Wen and Zee (1992):

δ=12πsωΓ,\delta\mathcal{L}=\frac{1}{2\pi}s\omega\partial\Gamma, (S14)

where ω\omega is the spin connection and the spin is s=(n1/2)s=(n-1/2). Including the coupling to curvature for the n¯2¯111\bar{n}\bar{2}111 parton state amounts to adding the following additional terms in the Lagrangian density of Eq. (S3):

δ=12πi=1n(i1/2)ωα(i)12πj=12(j1/2)ωβ(j)+12π32ωγ.\delta\mathcal{L}=-\frac{1}{2\pi}\sum_{i=1}^{n}(i-1/2)\omega\partial\alpha^{(i)}-\frac{1}{2\pi}\sum_{j=1}^{2}(j-1/2)\omega\partial\beta^{(j)}+\frac{1}{2\pi}\frac{3}{2}\omega\partial\gamma. (S15)

Using the constraints of Eqs. (S4) and (S5) we again eliminate α(n)\alpha^{(n)} and β(2)\beta^{(2)} to end up with the following additional term in the Lagrangian density of Eq. (I) which describes coupling of the n¯2¯111\bar{n}\bar{2}111 parton state to the curvature:

δ\displaystyle\delta\mathcal{L} =\displaystyle= 12π((1/2)ωα(1)+(3/2)ωα(2)++((2n3)/2)ωα(n2)+((2n1)/2)ω(γα(1)α(2)α(n1)))\displaystyle-\frac{1}{2\pi}\left((1/2)\omega\partial\alpha^{(1)}+(3/2)\omega\partial\alpha^{(2)}+\cdots+((2n-3)/2)\omega\partial\alpha^{(n-2)}+((2n-1)/2)\omega\partial(\gamma-\alpha^{(1)}-\alpha^{(2)}-\cdots-\alpha^{(n-1)})\right) (S16)
\displaystyle- 12π((1/2)ωβ(1)+(3/2)ω(γβ(1)))+12π32ωγ\displaystyle\frac{1}{2\pi}\left((1/2)\omega\partial\beta^{(1)}+(3/2)\omega\partial(\gamma-\beta^{(1)})\right)+\frac{1}{2\pi}\frac{3}{2}\omega\partial\gamma
=\displaystyle= 12π((n1)ωα(1)+(n2)ωα(2)++ωα(n1)+ωβ(1)((2n1)/2)ωγ)=12π𝔰IωaI.\displaystyle\frac{1}{2\pi}\left((n-1)\omega\partial\alpha^{(1)}+(n-2)\omega\partial\alpha^{(2)}+\cdots+\omega\partial\alpha^{(n-1)}+\omega\partial\beta^{(1)}-((2n-1)/2)\omega\partial\gamma\right)=\frac{1}{2\pi}\mathfrak{s}^{\rm I}\omega\partial a^{\rm I}.

Here we have defined the spin vector 𝔰=((n1),(n2),,1,1,((2n1)/2))T\mathfrak{s}=((n-1),(n-2),\cdots,1,1,-((2n-1)/2))^{\rm T}. The shift 𝒮\mathcal{S} on the sphere is given by Wen (1995):

𝒮=2νtTK1𝔰=1n,\mathcal{S}=\frac{2}{\nu}t^{\rm T}\cdot K^{-1}\cdot\mathfrak{s}=1-n, (S17)

which is consistent with the value ascertained from the microscopic wave function.

I.2 Specializing to the n=4n=4 case for the ν=4/9\nu=4/9 state

For the case of n=4n=4 the charge vector is t=(0,0,0,0,1)Tt=(0,0,0,0,1)^{\rm T} and the KK-matrix following Eq. (S10) is given by

K=(2110112101112010002111111).K=\begin{pmatrix}-2&-1&-1&0&1\\ -1&-2&-1&0&1\\ -1&-1&-2&0&1\\ 0&0&0&-2&1\\ 1&1&1&1&1\\ \end{pmatrix}. (S18)

The filling fraction is:

ν=tTK1t=K5,51=4/9,\nu=t^{\rm T}\cdot K^{-1}\cdot t=K^{-1}_{5,5}=4/9, (S19)

as expected. The ground state-degeneracy on a manifold with genus gg is:

ground-state degeneracy=|Det(K)|g=18g.\text{ground-state degeneracy}=|{\rm Det}(K)|^{g}=18^{g}. (S20)

This KK matrix has one positive and four negative eigenvalues which indicates that the 4¯2¯111\bar{4}\bar{2}111 state hosts one downstream and four upstream edge modes. The charges of the quasiparticles and their braiding statistics can be read off from the KK matrix Wen (1995).

The examples of the 4¯2¯111\bar{4}\bar{2}111 state considered in the present work and the nn¯111n\bar{n}111 states considered in Ref. Balram et al. (2020) show that in general a single-component Abelian fractional quantum Hall state at filling factor p/qp/q (with p,qp,q coprime) has the following properties:

  • Its ground state degeneracy on the torus is kqkq, where kk is a natural number.

  • The minimally charged quasiparticle has a charge (e)/(mq)(-e)/(mq), where mm is a natural number, and (e)(-e) is the charge of the electron.

In general, mkm\neq k. The 4¯2¯111\bar{4}\bar{2}111 state, which occurs at ν=4/9\nu=4/9, has k=2k=2, and m=1m=1. The nn¯111n\bar{n}111 states, which occur at ν=1/3\nu=1/3 for all nn, have k=n2k=n^{2} and m=nm=n Balram et al. (2020). Although these examples suggest that kk is a multiple of mm, we do not know if there exists a relation between kk and mm for a generic single-component Abelian state.

II Effect of finite width at ν=2+4/9\nu=2+4/9

We model the effect of finite thickness on the system by considering an infinite square well of width ww. The disk pseudopotentials of the Coulomb interaction, which are modified by the finite width ww, in the Landau level indexed by nn are given by

Vm\displaystyle V_{m} =\displaystyle= 0q𝑑qV(q)[Ln(q222)]2Lm(q22)eq22,\displaystyle\int_{0}^{\infty}qdq~{}V(q)\left[L_{n}\left(\frac{q^{2}\ell^{2}}{2}\right)\right]^{2}L_{m}(q^{2}\ell^{2})e^{-q^{2}\ell^{2}},
V(q)\displaystyle V(q) =\displaystyle= e2ϵqw/2w/2𝑑z1w/2w/2𝑑z2|ξ(z1)|2|ξ(z1)|2eq|z1z2|,\displaystyle\frac{e^{2}}{\epsilon q}\int^{w/2}_{-w/2}dz_{1}\int^{w/2}_{-w/2}dz_{2}~{}|\xi(z_{1})|^{2}|\xi(z_{1})|^{2}e^{-q|z_{1}-z_{2}|},
ξ(z)\displaystyle\xi(z) =\displaystyle= 2wcos(πzw),\displaystyle\sqrt{\frac{2}{w}}\cos\left(\frac{\pi z}{w}\right), (S21)

where ξ(z)\xi(z) is the wave function in the transverse direction and Ln(x)L_{n}(x) is the nnth order Laguerre polynomial. Analogously, one can also calculate the spherical pseudopotentials for a system with finite width. We shall present results for N=16N=16 electrons at the 4¯2¯111\bar{4}\bar{2}111 flux 2Q=392Q=39 obtained on the spherical geometry using both the disk and spherical pseudopotentials.

Refer to caption
Figure S1: (color online) Exact second Landau level Coulomb spectra N=16N=16 electrons at a flux of 2Q=392Q=39 in the spherical geometry evaluated using the a) disk pseudopotentials, b) disk pseudopotentials for width w=3w=3\ell, c) spherical pseudopotentials and d) spherical pseudopotentials for width w=3w=3\ell.

The exact second Landau level Coulomb spectra obtained using the disk and spherical pseudopotentials for w=0w=0 and w=3w=3\ell for N=16N=16 electrons at a flux 2Q=392Q=39 is shown in Fig. S1. The ground state, as well as the structure of the low-energy spectra obtained using the disk and spherical pseudopotentials, are quite different from each other, which indicates that finite-size effects are quite strong in the second Landau level. However, we find that the overlap between the lowest energy L=0L=0 state obtained using the spherical pseudopotentials and the ground state obtained using the disk pseudopotentials is 0.97\geq 0.97 for w[0,5]w\in[0,5]\ell. This shows that while the choice of pseudopotentials may affect selecting the global ground state in the SLL, the lowest energy L=0L=0 state itself is insensitive to the choice.

Using the disk pseudopotentials, we find that the ground state is uniform for w[0,5]w\in[0,5\ell]. Moreover, the pair-correlation function g(r)g(r) of the Coulomb ground state in the SLL for different widths have a visible but small difference indicating that they are fairly close to one another. Furthermore, the pair correlation function of the exact ground state obtained using the disk pseudopotentials agrees well with that of the 4¯2¯111\bar{4}\bar{2}111 parton state for the entire range of widths considered. With increasing width, the Coulomb g(r)g(r) looks more similar to the parton g(r)g(r) indicating that finite thickness enhances the stability of the parton state (see Fig. 2 of the main text). With the spherical pseudopotentials, the ground state has L=2L=2 for w=0w=0 and w=0.5w=0.5\ell. Encouragingly, as the quantum well width is increased, the ground state crosses from L=2L=2 to L=0L=0 in the range w(0.5,1]w\in(0.5,1]\ell and stays uniform for w[1,10]w\in[1,10]\ell. These results indicate that finite thickness aids in stabilizing the 4¯2¯111\bar{4}\bar{2}111 parton state in the second Landau level. For comparison in Fig. S2) we also show g(r)g(r) for the exact LLL Coulomb ground state which indicates that the ground state in the lowest two LLs for this system are very different from each other.

Refer to caption
Figure S2: (color online) The pair correlation function g(r)g(r) as a function of the arc distance for the exact second Landau level Coulomb ground state obtained using the disk pseudopotentials, and the 4¯2¯111\bar{4}\bar{2}111 state for N=16N=16 electrons at a flux of 2Q=392Q=39. For comparison we also show g(r)g(r) for the exact LLL Coulomb ground state for the same system.

We have also evaluated the charge gap as a function of finite width using both the disk and spherical pseudopotentials. To obtain the charge gap, we calculate the average charging energy 𝒞\mathcal{C} of the full LL per pair:

𝒞=li,jli,j|V|j,i(2l+1)2,\mathcal{C}=\sum_{-l\leq i,j\leq l}\frac{\langle i,j|V|j,i\rangle}{(2l+1)^{2}}, (S22)

where i,j|V|j,i\langle i,j|V|j,i\rangle is the direct matrix element of the interaction VV, l=Q+nl=Q+n is the shell angular momentum and nn is the LL index. Using the spherical symmetry of a filled LL, 𝒞\mathcal{C} can be simplified to

𝒞=ljli,j|V|j,i(2l+1),\mathcal{C}=\sum_{-l\leq j\leq l}\frac{\langle i,j|V|j,i\rangle}{(2l+1)}, (S23)

where ii can be chosen to be any of the orbitals, e.g., the one with m=lm=-l or m=lm=l. The charging energies associated with a model pseudopotential Hamiltonian with a single nonzero pseudopotential H=δm,jH=\delta_{m,j} can be expressed analytically as:

𝒞j=(4l2j+1)(2l+1)2.\mathcal{C}_{j}=\frac{(4l-2j+1)}{(2l+1)^{2}}. (S24)

Note that the 𝒞j\mathcal{C}_{j}’s are normalized, i.e., j=02l𝒞j=1\sum_{j=0}^{2l}\mathcal{C}_{j}=1. The charging energy for an arbitrary interaction specified by the set of pseudopotentials {Vm}\{V_{m}\} is therefore:

𝒞=m=02lVm(4l2m+1)(2l+1)2.\mathcal{C}=\sum_{m=0}^{2l}V_{m}\frac{(4l-2m+1)}{(2l+1)^{2}}. (S25)

For the spherical pseudopotentials the average charging energy for the 1/r1/r Coulomb interaction is equal to the inverse radius, i.e.,

𝒞sphere(1r)=e2ϵR=e2ϵQ=e2ϵln.\mathcal{C}^{\rm sphere}\left(\frac{1}{r}\right)=\frac{e^{2}}{\epsilon R}=\frac{e^{2}}{\epsilon\ell\sqrt{Q}}=\frac{e^{2}}{\epsilon\ell\sqrt{l-n}}. (S26)

In particular, the charging energy for our system of interest at 2l=392l=39 for w=0w=0 is e2/(ϵ39/2)e^{2}/(\epsilon\ell\sqrt{39/2}) and e2/(ϵ37/2)e^{2}/(\epsilon\ell\sqrt{37/2}) in the LLL and SLL respectively. In particular, the charging energy for our system of interest at 2l=392l=39 for w=0w=0 is e2/(ϵ39/2)e^{2}/(\epsilon\ell\sqrt{39/2}) and e2/(ϵ37/2)e^{2}/(\epsilon\ell\sqrt{37/2}) in the LLL and SLL respectively. Note that for the Coulomb interaction V(r)=1/rV(r)=1/r, Eq. (S26) leads to a sum rule for the Coulomb pseudopotentials {VmCoulomb}\{V^{\rm Coulomb}_{m}\}’s, namely m=02lVmCoulomb(4l2m+1)=(2l+1)2/R\sum_{m=0}^{2l}V^{\rm Coulomb}_{m}(4l-2m+1)=(2l+1)^{2}/R. For the planar disk pseudopotentials the average charging energy for the 1/r1/r Coulomb interaction in the lowest two LLs, in terms of the shell angular momentum ll, is given by

𝒞disk,LLL(1r)\displaystyle\mathcal{C}^{\rm disk,~{}LLL}\left(\frac{1}{r}\right) =\displaystyle= (3+4(2l))Γ[2l+3/2]3(2l+1)2Γ[2l+1]e2ϵ\displaystyle\frac{(3+4(2l))\Gamma[2l+3/2]}{3(2l+1)^{2}\Gamma[2l+1]}\frac{e^{2}}{\epsilon\ell} (S27)
𝒞disk,SLL(1r)\displaystyle\mathcal{C}^{\rm disk,~{}SLL}\left(\frac{1}{r}\right) =\displaystyle= (256(2l)3+20(2l)33)Γ[2l1/2]192(2l+1)2Γ[2l+1]e2ϵ,\displaystyle\frac{(256(2l)^{3}+20(2l)-33)\Gamma[2l-1/2]}{192(2l+1)^{2}\Gamma[2l+1]}\frac{e^{2}}{\epsilon\ell},

where Γ[x]\Gamma[x] is the Gamma function.

The value of the charging energy of the finite-width interaction has to be numerically evaluated since no closed-form analytic expression for the finite-width pseudopotentials is known. We find that the 𝒞disk(1/r)<𝒞sphere(1/r)\mathcal{C}^{\rm disk}(1/r)<\mathcal{C}^{\rm sphere}(1/r) since the disk pseudopotentials decrease more quickly with relative angular momentum mm than the spherical pseudopotentials. Throughout this work we consider only fully spin-polarized electrons, therefore, the electron-electron interaction energies only depend on the odd pseudopotentials. However, the average charging energy, which depends on the direct matrix element of the interaction, is a function of both the even and odd pseudopotentials. Note that it is customary to quote the value of 2l2l itself as the flux 2Q2Q since the numerics on the SLL are simulated using effective pseudopotentials in the LLL (and in the LLL 2l=2Q2l=2Q).

Using the average charging energies we compute the charged gaps as follows:

(2Q1)\displaystyle\mathcal{E}(2Q-1) =\displaystyle= E(2Q1)(N2(nqeq)2)𝒞(2Q1)2,\displaystyle E(2Q-1)-(N^{2}-(n_{q}e_{q})^{2})\frac{\mathcal{C}(2Q-1)}{2},
(2Q)\displaystyle\mathcal{E}(2Q) =\displaystyle= E(2Q)(N2)𝒞(2Q)2,\displaystyle E(2Q)-(N^{2})\frac{\mathcal{C}(2Q)}{2},
(2Q+1)\displaystyle\mathcal{E}(2Q+1) =\displaystyle= E(2Q+1)(N2(nqeq)2)𝒞(2Q+1)2,\displaystyle E(2Q+1)-(N^{2}-(n_{q}e_{q})^{2})\frac{\mathcal{C}(2Q+1)}{2},
Δcharge\displaystyle\Delta^{\rm charge} =\displaystyle= (2Q1)+(2Q+1)2(2Q)nq,\displaystyle\frac{\mathcal{E}(2Q-1)+\mathcal{E}(2Q+1)-2\mathcal{E}(2Q)}{n_{q}}, (S28)

where E(2Q)E(2Q) are ground state energies obtained from exact diagonalization of NN electrons at flux 2Q2Q, nq=4n_{q}=4 is the number of quasiholes or quasiparticles produced per flux quantum and eq=1/9e_{q}=1/9 is the magnitude charge of the quasihole or quasielectron in units of the electronic charge. The N2N^{2} term includes the background contribution, and the (nqeq)2(n_{q}e_{q})^{2} term accounts for the fact that the background is different when some charge is accumulated in form of quasiholes or quasielectrons. For the spherical pseudopotentials for w=0w=0 and w=0.5w=0.5\ell we take the lowest energy L=0L=0 state as the ground state to evaluate the charge gap. The neutral gap is defined as the difference between the two lowest energy states at the flux 2Q2Q corresponding to the ground state. We multiply the energies by a factor of (2l2n)ν/N\sqrt{(2l-2n)\nu/N} that corrects for the deviation of the electron density of a finite system from its NN\rightarrow\infty value before extrapolating them to the thermodynamic limit Morf et al. (1986).

We find that both the charge and neutral gaps evaluated using the disk pseudopotentials first increase as the well-width is increased and then slowly decreases as the well-width is further increased. In the entire range of widths we considered, the charge and neutral gaps obtained from the disk pseudopotentials are of the order of 0.01e2/(ϵ)0.01~{}e^{2}/(\epsilon\ell) and 0.005e2/(ϵ)0.005~{}e^{2}/(\epsilon\ell) respectively (see Fig. S3). With the spherical pseudopotentials, the charge gap evaluated using the lowest energy L=0L=0 state is negative for w=0w=0 and w=0.5w=0.5\ell but in Fig. S3 we have shown only its magnitude to follow the conventional way of indicating phase transitions. In the range from w=1w=1\ell to w=5w=5\ell, we find that both the charge and neutral gaps evaluated using the spherical pseudopotentials also monotonically increase as the well-width is increased These results corroborate the fact that the finite thickness of the quantum well enhances the stability of the 2+4/92+4/9 FQHE.

Refer to caption
Figure S3: (color online) Charge and neutral gaps for N=16N=16 electrons at a flux of 2Q=392Q=39 evaluated in the spherical geometry using the disk and spherical pseudopotentials for various widths ww.

For comparison, in Figs. S4 and S5 we show the neutral and charge gaps obtained using the disk and spherical pseudopotentials for various widths at the Laughlin (extrapolated from N=8,9,,15N=8,9,\cdots,15) and Pfaffian (extrapolated from N=8,10,,20N=8,10,\cdots,20) fluxes in the SLL (see also Refs. Morf (1998); Wójs (2009)). For completeness, on the same plots, we have also included the corresponding LLL gaps. The ground state at the half-filled LLL is a compressible Fermi liquid state of composite fermions. Therefore, the gap in the LLL at the Pfaffian flux is irregular as a function of 1/N1/N (the gap is only large when the flux and number of electrons are aliased with a Jain CF state), so we have not extrapolated it to the thermodynamic limit. We find that the gaps for the largest systems at 7/37/3 and 5/25/2 are a factor of 2-4 times larger than those for the 1616-particle 2+4/92+4/9 state. This indicates that the 2+4/92+4/9 state is more fragile as compared to the experimentally observed 7/37/3 and 5/25/2 states.

Refer to caption
Figure S4: (color online) Neutral gaps in the second Landau level (blue) evaluated in the spherical geometry using the disk and spherical pseudopotentials for various widths ww at the 1/31/3 Laughlin (top panels) and 1/21/2 Pfaffian (bottom panels) fluxes. Panels a) and d) [b) and e)] show linear extrapolation of the neutral gap as a function of 1/N1/N for w=0w=0 [w=3w=3\ell]. The extrapolated gaps are shown in panels c) and f). For completeness, we have also shown the corresponding lowest Landau level gaps (red).
Refer to caption
Figure S5: (color online) Charge gaps in the second Landau level (blue) evaluated in the spherical geometry using the disk and spherical pseudopotentials for various widths ww at the 1/31/3 Laughlin (top panels) and 1/21/2 Pfaffian (bottom panels) fluxes. Upon the insertion of a single flux quantum, the 1/31/3 Laughlin state produces a single quasihole while the 1/21/2 Pfaffian state produces two quasiholes. Thus we use nq=1n_{q}=1 and nq=2n_{q}=2 for the 1/31/3 Laughlin and 1/21/2 Pfaffian states respectively [see Eq. (S28)]. Panels a) and d) [b) and e)] show linear extrapolation of the charge gap as a function of 1/N1/N for w=0w=0 [w=3w=3\ell]. The extrapolated gaps are shown in panels c) and f). For completeness, we have also shown the corresponding lowest Landau level gaps (red).

III Comparison between some candidate and exact ground states in the second Landau level

The wave function for the 4¯2¯111\bar{4}\bar{2}111 parton state is most easily evaluated in the LLL in the first-quantized coordinate space (as opposed to the Hilbert space of Fock states). Therefore, to evaluate its second LL Coulomb energy, we use the effective interaction of Ref. Shi et al. (2008) to simulate the physics of the second LL in the LLL. This effective interaction has nearly the same Haldane pseudopotentials as that of the Coulomb interaction in the second LL. Including the contribution of the positively charged background, the per-particle density corrected Morf et al. (1986) energy of the 4¯2¯111\bar{4}\bar{2}111 parton state for the effective interaction for the system of N=16N=16 particles is 0.3839(1)-0.3839(1) (the number in the parenthesis is the Monte Carlo uncertainty in the energy estimate) while the exact energy of the lowest-lying L=0L=0 state is 0.3858-0.3858, both in Coulomb units of e2/(ϵ)e^{2}/(\epsilon\ell). The level of agreement (within 0.49%0.49\%) between these two numbers is comparable with that of other trial states in the second LL Hutasoit et al. (2017); Balram et al. (2018b). The rest of the section is devoted to showing results for other candidate states and their comparison with the exact second LL Coulomb ground states.

In Tables S1 and S2 we compare the overlap and energies of the Coulomb ground state in the two lowest Landau levels with the Laughlin and Pfaffian state respectively. For the largest system of N=15N=15 electrons considered in this work, the energy of the Laughlin state differs from the exact second LL Coulomb energy by about 1.1%1.1\%. Similarly, the energy of the Pfaffian state of N=20N=20 electrons is within 0.6%0.6\% of the exact second LL Coulomb ground state. These numbers are of the same order of magnitude as those of the 4¯2¯111\bar{4}\bar{2}111 parton state shown above. For completeness, in Figs. S6 and S7 we show the thermodynamic extrapolation and the extrapolated per-particle Coulomb energies of the Laughlin and Pfaffian states in the lowest and second Landau level obtained using the spherical and disk pseudopotentials at various finite widths.

Refer to caption
Figure S6: (color online) Thermodynamic extrapolations of the Coulomb per-particle energies for the ν=1/3\nu=1/3 Laughlin state in the lowest and second Landau level obtained using the spherical and disk pseudopotentials at various finite widths. The extrapolated energies are shown in panel c) are obtained for the spherical pseudopotentials from quadratic fits in 1/N1/N.
Refer to caption
Figure S7: (color online) Thermodynamic extrapolations of the Coulomb per-particle energies for the ν=1/2\nu=1/2 Pfaffian state in the lowest and second Landau level obtained using the spherical and disk pseudopotentials at various finite widths. The extrapolated energies shown in panel c) are obtained for the spherical pseudopotentials from quadratic fits in 1/N1/N.
Refer to caption
Figure S8: (color online) Overlap of the exact Laughlin state (left panel) and Pfaffian state (right panel) with the exact Coulomb ground state evaluated in the spherical geometry using the disk and spherical pseudopotentials for various widths ww for the largest system sizes considered in this work.
NN 2Q2Q |ψ1/31LL|ψ1/3L||\langle\psi_{1/3}^{\rm 1LL}|\psi_{1/3}^{L}\rangle| ψ1/31LL|H1LL|ψ1/31LL\langle\psi_{1/3}^{\rm 1LL}|H^{\rm 1LL}|\psi_{1/3}^{\rm 1LL}\rangle ψ1/3L|H1LL|ψ1/3L\langle\psi_{1/3}^{L}|H^{\rm 1LL}|\psi_{1/3}^{L}\rangle |ψ1/30LL|ψ1/3L||\langle\psi_{1/3}^{\rm 0LL}|\psi_{1/3}^{L}\rangle| ψ1/30LL|H0LL|ψ1/30LL\langle\psi_{1/3}^{\rm 0LL}|H^{\rm 0LL}|\psi_{1/3}^{\rm 0LL}\rangle ψ1/3L|H0LL|ψ1/3L\langle\psi_{1/3}^{L}|H^{\rm 0LL}|\psi_{1/3}^{L}\rangle |ψ1/31LL|ψ1/30LL||\langle\psi_{1/3}^{\rm 1LL}|\psi_{1/3}^{\rm 0LL}\rangle|
6 15 0.52848 -0.33306 -0.32604 0.99645 -0.41095 -0.41075 0.58814
7 18 0.60705 -0.33017 -0.32611 0.99636 -0.41082 -0.41062 0.64941
8 21 0.57197 -0.33071 -0.32615 0.99540 -0.41074 -0.41052 0.62836
9 24 0.47941 -0.33101 -0.32617 0.99405 -0.41068 -0.41044 0.53871
10 27 0.54000 -0.33043 -0.32619 0.99295 -0.41063 -0.41038 0.60592
11 30 0.70303 -0.32954 -0.32620 0.99217 -0.41058 -0.41032 0.77057
12 33 0.50296 -0.33039 -0.32621 0.99092 -0.41055 -0.41028 0.57650
13 36 0.54445 -0.33001 -0.32622 0.98977 -0.41052 -0.41024 0.62202
14 39 0.57706 -0.32984 -0.32623 0.98871 -0.41049 -0.41020 0.66098
15 42 0.52978 -0.32991 -0.32623 0.98763 -0.41047 -0.41018 0.61475
Table S1: Table shows the overlap of the exact Coulomb ground state |ψ1/31LL|\psi_{1/3}^{\rm 1LL}\rangle for NN electrons at the Laughlin flux 2Q=3N32Q=3N-3 in the second Landau level with the Laughlin state |ψ1/3L|\psi_{1/3}^{L}\rangle. We also show the per-particle density corrected second Landau level Coulomb energies (which includes the contribution of the positively charged background) of the exact and Laughlin states in Coulomb units of e2/(ϵ)e^{2}/(\epsilon\ell). For completeness, we have also included analogous numbers in the lowest Landau level. Some of these results were previously quoted in Refs. d’Ambrumenil and Reynolds (1988); Balram et al. (2013); Kuśmierz and Wójs (2018).
NN 2Q2Q |ψ1/21LL|ψ1/2MR||\langle\psi_{1/2}^{\rm 1LL}|\psi_{1/2}^{\rm MR}\rangle| ψ1/21LL|H1LL|ψ1/21LL\langle\psi_{1/2}^{\rm 1LL}|H^{\rm 1LL}|\psi_{1/2}^{\rm 1LL}\rangle ψ1/2MR|H1LL|ψ1/2MR\langle\psi_{1/2}^{\rm MR}|H^{\rm 1LL}|\psi_{1/2}^{\rm MR}\rangle |ψ1/20LL|ψ1/2MR||\langle\psi_{1/2}^{\rm 0LL}|\psi_{1/2}^{\rm MR}\rangle| ψ1/20LL|H0LL|ψ1/20LL\langle\psi_{1/2}^{\rm 0LL}|H^{\rm 0LL}|\psi_{1/2}^{\rm 0LL}\rangle ψ1/2MR|H0LL|ψ1/2MR\langle\psi_{1/2}^{\rm MR}|H^{\rm 0LL}|\psi_{1/2}^{\rm MR}\rangle |ψ1/21LL|ψ1/20LL||\langle\psi_{1/2}^{\rm 1LL}|\psi_{1/2}^{\rm 0LL}\rangle|
6 9 0.81624 -0.37528 -0.37186 0.92280 -0.46734 -0.46275 0.53064
8 13 0.86739 -0.37195 -0.36966 0.92130 -0.46175 -0.45921 0.61067
10 17 0.83764 -0.36968 -0.36795 0.90537 -0.46239 -0.45931 0.55247
12 21 0.81939 -0.36900 -0.36678 0.65520 -0.46651 -0.45922 0.24462
14 25 0.69345 -0.36850 -0.36621 0.72226 -0.46367 -0.45875 0.28091
16 29 0.77945 -0.36767 -0.36574 0.74587 -0.46297 -0.45854 0.33699
18 33 0.67656 -0.36764 -0.36536 0.63554 -0.46400 -0.45845 0.20746
20 37 0.67364 -0.36723 -0.36508 0.37029 -0.46621 -0.45834 0.08814
Table S2: Table shows the overlap of the exact Coulomb ground state |ψ1/21LL|\psi_{1/2}^{\rm 1LL}\rangle for NN electrons at the Moore-Read Pfaffian flux 2Q=2N32Q=2N-3 in the second Landau level with the Moore-Read Pfaffian state |ψ1/2MR|\psi_{1/2}^{\rm MR}\rangle. We also show the per-particle density corrected second Landau level Coulomb energies (which includes the contribution of the positively charged background) of the exact and Moore-Read Pfaffian states in Coulomb units of e2/(ϵ)e^{2}/(\epsilon\ell). For completeness, we have also included analogous numbers in the lowest Landau level, where overlaps and energies are shown for the lowest energy L=0L=0 state ( indicates cases where the lowest energy L=0L=0 state is also the ground state in the lowest Landau level). Some of these results were previously quoted in Refs. Morf (1998); Scarola et al. (2002); Pakrouski et al. (2015); Kuśmierz and Wójs (2018).

III.1 Updating results of Ref. Balram et al. (2018b)

The 3¯2¯111\bar{3}\bar{2}111 parton state, which is the predecessor of the 4¯2¯111\bar{4}\bar{2}111 parton state in the n¯2¯111\bar{n}\bar{2}111 parton sequence, was posited as a candidate Balram et al. (2018b) to describe the experimentally observed FQHE ν=2+6/13\nu=2+6/13 FQHE Kumar et al. (2010). In Ref. Balram et al. (2018b), the 3¯2¯111\bar{3}\bar{2}111 parton state was compared against the exact SLL Coulomb ground state for only a single system of N=12N=12 electrons at a flux of 2Q=282Q=28. At the time of publication, the next system size, that of N=18N=18 electrons at a flux of 2Q=412Q=41, which has a Hilbert space dimension of 3.5×1093.5\times 10^{9}, was not accessible to exact diagonalization. We have now been able to evaluate the exact SLL Coulomb ground state for this system size, which, to the best of our knowledge, is the largest system size for which the ground state of an FQHE Hamiltonian has been obtained. Improvements in the matrix times vector multiplication at each Lanczos iteration, which have made it possible to efficiently diagonalize such large matrices, have been described in Ref. Zuo et al. (2020).

The exact SLL Coulomb ground state for N=18N=18 electrons at a flux of 2Q=412Q=41 has L=0L=0 that is consistent with realizing a uniform state. This system has a neutral gap of 0.0096 e2/(ϵ)e^{2}/(\epsilon\ell). In Fig. S9 we compare the pair-correlation function of the exact SLL Coulomb ground state with that of the 3¯2¯111\bar{3}\bar{2}111 parton state. The pair-correlation functions of both the 3¯2¯111\bar{3}\bar{2}111 parton state as well as the exact SLL Coulomb ground state show oscillations that decay at long distances, which is a typical feature of an incompressible state Kamilla et al. (1997); Balram et al. (2015c). Furthermore, g(r)g(r) of the exact SLL Coulomb ground state and the parton trial state match reasonably well with each other. Although the agreement is not conclusive, it is on par with that of other candidates in the SLL d’Ambrumenil and Reynolds (1988); Morf (1998); Scarola et al. (2002); Balram et al. (2013); Pakrouski et al. (2015); Hutasoit et al. (2017); Kuśmierz and Wójs (2018). For example, in Figs. S10 and S11 we show a comparison of the g(r)g(r)’s of the Laughlin state Laughlin (1983) with that of the ν=7/3\nu=7/3 ground state and of the Pfaffian state Moore and Read (1991) with that of the ν=5/2\nu=5/2 ground state for the largest systems considered in this work. Thus, we conclude that the results for the N=18N=18 electrons at flux 2Q=412Q=41 provide further evidence for the feasibility of the 3¯2¯111\bar{3}\bar{2}111 parton state to describe the 2+6/132+6/13 FQHE.

Refer to caption
Figure S9: (color online) The pair correlation function g(r)g(r) as a function of the arc distance for the exact second Landau level Coulomb ground state, and the 3¯2¯111\bar{3}\bar{2}111 parton candidate state for N=18N=18 electrons at a flux of 2Q=412Q=41. For comparison we also show g(r)g(r) for the exact LLL Coulomb ground state as well as that of the ν=2/5\nu=2/5 Jain CF state for the same system size.
Refer to caption
Figure S10: (color online) The pair correlation function g(r)g(r) as a function of the arc distance for the exact second Landau level Coulomb ground state, and the Laughlin state for N=15N=15 electrons at a flux of 2Q=422Q=42. For comparision we also show g(r)g(r) for the exact LLL Coulomb ground state for the same system.
Refer to caption
Figure S11: (color online) The pair correlation function g(r)g(r) as a function of the arc distance for the exact second Landau level Coulomb ground state, and the Pfaffian state for N=20N=20 electrons at a flux of 2Q=372Q=37. For comparison we also show g(r)g(r) for the exact LLL Coulomb ground state as well as that of the ν=4/7\nu=4/7 Jain CF state for the same system size.

Including the contribution of the positively charged background, the per-particle density corrected Morf et al. (1986) energy of the 3¯2¯111\bar{3}\bar{2}111 parton state for the effective interaction that we use to simulate the physics of the SLL in the LLL for the system of N=18N=18 particles at flux 2Q=412Q=41 is 0.3855(1)-0.3855(1) while the exact energy of the ground state is 0.3879-0.3879, both in Coulomb units of e2/(ϵ)e^{2}/(\epsilon\ell). Again, the level of agreement (within 0.64%0.64\%) between these two numbers is comparable with that of other candidate wave functions in the second LL Hutasoit et al. (2017); Balram et al. (2018b).

For completeness, we have also evaluated the exact LLL Coulomb ground state for the N=18N=18 system at the flux of 2Q=412Q=41. The overlaps of the LLL and SLL Coulomb ground states is 0.08720.0872 and their g(r)g(r)’s are also very different from each other (see Fig. S9) which suggests that the ground state in the two LLs are quite different from each other for this system. We note that the system of N=18N=18 electrons at a flux of 2Q=412Q=41 aliases with the ν=2/5\nu=2/5 Jain CF state. As expected, we find excellent agreement between the g(r)g(r)’s of the LLL Coulomb ground state and the ν=2/5\nu=2/5 Jain CF state for this system (see Fig. S9). Microscopically, the Jain CF states provide a very accurate representation of the exact LLL Coulomb ground states whereas the candidate states in the SLL provide only an approximate representation of the exact SLL Coulomb ground state. The neutral gap for the system of N=18N=18 electrons at the flux 2Q=412Q=41 in the LLL is 0.0487 e2/(ϵ)e^{2}/(\epsilon\ell) (five times larger than the corresponding SLL neutral gap quoted above) indicating that the 2/5 Jain CF state is more robust than the 6/136/13 parton state.