# Entanglement Measures for Quasi-Two-Dimensional Fractional Quantum Hall States

###### Abstract

We theoretically examine entanglement in fractional quantum hall states, explicitly taking into account and emphasizing the quasi-two-dimensional nature of experimental quantum Hall systems. In particular, we study the entanglement entropy and the entanglement spectrum as a function of the finite layer thickness of the quasi-two-dimensional system for a number of filling fractions in the lowest and the second Landau levels: = 1/3, 7/3, 1/2, and 5/2. We observe that the entanglement measures are dependent on which Landau level the electrons fractionally occupy, and find that filling fractions 1/3 and 7/3, which are considered to be Laughlin states, weaken with in the lowest Landau level (=1/3) and strengthen with in the second Landau level (=7/3). For the enigmatic even-denominator state, we find that entanglement in the ground state is consistent with that of the non-Abelian Moore-Read Pfaffian state at an optimal thickness . We also find that the single-layer system is not a fractional quantum Hall state consistent with the experimental observation. In general, our theoretical findings based on entanglement considerations are completely consistent with the results based on wavefunction overlap calculations.

###### pacs:

73.43.-f, 71.10.Pm## I Introduction

The discovery of the fractional quantum hall effect (FQHE) in 1982Tsui et al. (1982) has proven to be one of the most significant experimental findings in all of physics. The incompressible quantum fluid that manifests in a two-dimensional (2D) electron system in the presence of a strong perpendicular magnetic field at low temperatures can not be explained by the “conventional” Landau theory of phase transitionsLaughlin (1983). Instead this unique phase is an example of an emergent topological state of matter, where the charged quasi-particle excitations are governed by anyonic rather than fermionic or bosonic statisticsHalperin (1984); Leinaas and Myrheim (1977); Wilczek (1982); Arovas et al. (1984). The unique nature of FQHE becomes obvious with the realization that the phenomenon occurs only in the truly strongly interacting limit of vanishing (or extremely small) kinetic energy with the noninteracting ground state having a macroscopically large degeneracy. Although the FQHE was discovered decades agoDas Sarma and Pinczuk (1997); Prange and Girvin (1990); Jain (2007), many questions still linger to this day. The most notable of these is the nature of the FQH state observed at electron filling fraction . First discovered in 1987Willett et al. (1987), the state is (so far) the only exception to the famous “odd denominator” rule for monolayer FQH systems given by the Laughlin ansatzLaughlin (1983) and the more general “composite fermion” theoryJain (1989, 2007). (The odd denominator rule assumes FQHE in a monolayer system. Bilayer FQH systems, however, have several even denominator states that are well-understood. See, for example Refs. Suen et al., 1992; Eisenstein et al., 1992; He et al., 1993). Currently, the leading theoretical candidate for the state is the Moore-Read (MR) Pfaffian model state introduced by Moore and Read in 1990Moore and Read (1991). Recently it has been noted that the particle-hole conjugate to the MR Pfaffian, the so-called anti-Pfaffian Lee et al. (2007); Levin et al. (2007), is also a viable candidate for the FQHE at . An exceptional feature of these ansatz is the prediction of anyonic quasi-particles with non-Abelian braiding statistics Moore and Read (1991); Nayak and Wilczek (1996); Bonderson et al. (2011). This prediction is attractive given the current interest in topological quantum computing Das Sarma et al. (2005); Nayak et al. (2008), but the true nature of the FQH state is, by no means, a settled questionStorni et al. (2010); Storni and Morf (2011); Wójs (2009). Directly probing the topological nature of FQH states has proven to be both a theoretical and an experimental challenge, making definitive and uncontroversial verification of the MR Pfaffian ansatz elusive. However, recent developments in the field of quantum information have shown that measures of entanglement are useful tools in examining the global features of many-body strongly correlated quantum statesAmico et al. (2008).

The true nature of the FQHE state is one of the most prominent open questions in condensed matter physics, and is a primary motivation for the current study. Recent experimental studies have explored this mysterious state and give some weight to the Moore-Read theory that is believed to explain it. One of these studies is the recent experiment by Venkatachalam et. al. Venkatachalam et al. (2010) that measured the charge of localized excitations in the state to be as predicted by the MR theoryMoore and Read (1991). These results are consistent with previous studies by Radu et. al. and Dolev et. al. that used shot noise to investigate the local charge M. Dolev and Mahalu (2008); Radu et al. (2008); Dolev et al. (2010). In addition to the experiments noted above, Willett et. al. Willett et al. (2009, 2009) have seen evidence of quasiparticle interference oscillations that support the existence of charge excitations at . Another recent study is work performed by Bid et. al.Bid et al. (2010) that experimentally observed the theorized neutral mode of the state consistent with the MR theoryMilovanović and Read (1996); Overbosch and Wen (2008). Although these developments point to the MR theory as the likely candidate for the state, they are not sufficient to unambiguously establish the existence of non-Abelian anyons. Also, it should be noted that the MR theory predicts a spin polarized state at , but recent experimental work Rhone et al. (2011); Stern et al. (2010) suggest that this state may be unpolarized in some cases. Although we do not consider spin in this work, we attempt to leverage recent developments in quantum information theory to observe topological features predicted by the Laughlin and MR theory in numerically obtained FQHE ground states and examine how these features change under a realistic change (in particular, the finite thickness effect with varying thickness) in the effective interaction.

These developments in quantum information theory provide the context for our current work. In this article we provide a detailed numerical theoretical study of entanglement for FQHE states incorporating the dependence on the quasi-2D layer thickness of the transverse dimension (i.e the finite layer thickness effectMacDonald and Aers (1984); Zhang and Das Sarma (1986); He et al. (1990); Ortalano et al. (1997); Park and Jain (1998); Peterson et al. (2008a, b)). We emphasize physics of the finite-thickness effect in this article partly because it provides a qualitative understanding of the FQHE in higher Landau levels as shown in Ref. Peterson et al., 2008a, b. Indeed, the orbital Landau level (LL) dependence of the FQHE is not completely understood Kumar et al. (2010). The theory behind the MR Pfaffian model state, for example, makes no distinction between the half-filled second orbital Landau level (SLL) (i.e., ) and the half-filled lowest orbital Landau level (LLL) (), but no FQHE has so far been observed at in monolayer systems. Furthermore FQH states in the SLL are relatively rare and generally “weaker” (i.e., requiring lower temperatures and higher sample mobilities to experimentally observe due to their relatively small energy gaps on the order of 0.5 K or less) compared to FQH states in the LLLWillett et al. (1987); Pan et al. (1999); Eisenstein et al. (2002); Xia et al. (2004); Csáthy et al. (2005); Choi et al. (2008); Pan et al. (2008). Some light was shed on this phenomenon of fragility of the FQHE in the SLL in Ref. Peterson et al., 2008a, b where it is shown that the finite-thickness effect is qualitatively dependent on Landau level. In particular, it is shown that a finite, non-zero, layer thickness, , helps stabilize the FQHE in the SLL, whereas in the LLL, finite tends to “weaken” the state. Also, for the state, there appears to be an optimal thickness where the MR Pfaffian is “strongest.” Thus, there seems to be a close connection between finite thickness and the SLL FQHE, which we explore in this work by calculating the thickness dependent entanglement properties of the FQHE. The results given in Ref. Peterson et al., 2008a, b are primarily based on calculations of the overlap between numerically obtained ground states and model FQHE wavefunctions foo () — in particular, the Laughlin and MR Pfaffian model states. Although the overlap is a powerful theoretical tool — in fact, the wide acceptance of the Laughlin wavefunction as the appropriate description of the experimentally observed odd-denominator FQHE is arguably based on overlap results — it is not always definitive and can be misleading in some cases. An example of this involves the state that has been shown to have a large overlap with both the Jain composite fermion wavefunction as well as the so-called “Gaffnian” wavefunction, even though these two states have different underlying topological order belonging to different universality classes Simon et al. (2007); Regnault et al. (2009). Therefore, we seek alternative theoretical tools to probe the Landau level dependence of the FQHE through the finite-thickness effect, which should inevitably give insight into the nature of the enigmatic state. In this article, we examine bipartite entanglement as an alternative measure (in particular the entanglement entropy and the entanglement spectrum) and study its dependence on finite layer thickness, , in quasi-2D FQH states. One reason for our study of the thickness-dependent FQHE entanglement is that the thickness parameter enables a continuous tuning of the Hamiltonian, lending to a continuous variation in the entanglement allowing a comprehensive systematic study.

Bipartite entanglement measures are tools designed to quantify the extent to which degrees of freedom are entangled in a bipartitioned systemPlenio and Virmani (2007). The most straight forward of these is the entanglement entropy (EE), defined as the the von-Neuman entropy

(1) |

of the reduced density matrix for state in a Fock space that has been partitioned into two parts. The EE has proven to be a very powerful tool in examining quantum correlations in interacting many-body systemsAmico et al. (2008). In particular, the scaling of the entanglement entropy with system size has been shown to follow certain “area” laws that can identify quantum phase transitions in certain casesEisert et al. (2010). Also, the EE can be used to extract the “topological entanglement entropy” which is an indicator of topological order in the systemKitaev and Preskill (2006); Levin and Wen (2006) (i.e. states in the same topological class will have the same topological entanglement entropy). Given these developments, the EE and the topological EE appear to be attractive tools to probe FQH states. However, as we discuss below in more detail, obtaining a precise quantitative estimate of the topological entropy in FQH states requires a technical procedure prone to introducing significant errors Zozulya et al. (2007); Iblisdir et al. (2007); Läuchli et al. (2010a). In our study, we do not attempt to make such an estimate. Instead we focus on FQHE states with the primary objective of observing the qualitative behavior of the EE as a function of finite layer thickness, , in the LLL and SLL. Our goal is to obtain the qualitative dependence of the EE in FQHE as a function of layer thickness .

We also investigate the “entanglement spectrum” of quasi-2D FQH states in this study. Introduced by Li and HaldaneLi and Haldane (2008), the entanglement spectrum (ES) provides more information than the entanglement entropy alone. In particular, Li and Haldane conjectured that there is a direct correspondence between the low-lying eigenvalues of the operator and the edge modes of the system and, thus, ES can be used to determine characteristics of the underlying conformal field theory (CFT) of the corresponding FQH ground state. As long as these “CFT-like” states are well-separated from the “generic non-CFT-like” states in the entanglement spectrum of the ground state in the thermodynamic limit, then it is conjectured that the identified CFT does, indeed, describe the state. In other words, states described the by the same CFT (i.e. in the same universality class) will have the same low-lying structure in their respective entanglement spectra. Numerical studies Hermanns et al. (2010); Chandran et al. (2011); Zhao et al. (2011); Sterdyniak et al. (2011); Papić et al. (2011) and some recent analytical results Qi et al. (2011) support the Li and Haldane conjecture. In our study, we carefully examine the entanglement spectrum of finite sized FQHE states by quantifying the separations between the suspected low-lying CFT and non-CFT-like states (“entanglement gaps”) as a function of finite-layer thickness. The entanglement gaps serve as a semi-quantitative measure of how well the state in a finite system fits with the universality class described by the CFT. By doing this we discover that the entanglement gaps follow trends qualitatively similar to the EE, leading to similar conclusions obtained in the overlap calculations given in Ref. Peterson et al., 2008a, b (i.e., that the finite-thickness effect strengthens the FQH states in the SLL). We also examine the entanglement spectrum in the so-called “conformal limit” Thomale et al. (2010). The aim of the conformal limit is to remove finite size artifacts from the geometry of the system (in our case, the sphere), allowing for a concrete definition of a full entanglement gap. We find that the conformal limit does result in a full entanglement gap in most cases studied, but not in all situations. The lack of a entanglement gap in these exceptional cases, however, is likely due to our choice of planar pseudo-potentials in our workPapic (2011); further work along this line would be necessary to fully understand these situations.

The structure of this article is as follows: In section II we describe our methods for numerically obtaining the exact ground state of the Coulombic Hamiltonian for FQH ground states at filling fractions , 7/3, 1/2, and 5/2 and our chosen model of the finite thickness of the quasi-2D system as well as the model wavefunctions (Laughlin and MR Pfaffian) to which we compare the numerically obtained ground state. Further, we describe how our entanglement measures are defined and calculated. In section III we provide our results for the entanglement entropy (III.1), entanglement gaps in the entanglement spectrum (III.2), entanglement spectra in the conformal limit (III.3) and density of states calculations of entanglement spectra (III.4) for FQH ground states at the filling fractions listed above. We provide our conclusions in section IV. Finally in the appendix, we discuss our choice of planar Haldane pseudopotentials over spherical pseudopotentials and examine the implications of this choice by comparing entanglement spectra in FQHE ground states at zero thickness obtained with either choice of pseudopotentials.

## Ii Method

We begin by considering a quasi-2D geometry where spinless electrons are confined in the x-y plane with layer thickness along the z-axis and an external magnetic field also along the z-axis. In the non-interacting case, the presence of the magnetic field quantizes the electron energy levels into highly degenerate Landau levels (LL) with energies where and is a non-negative integer that defines the LL index. The degeneracy of each LL per unit area is given by , where is the magnetic length and defines a length scale for the problem. The filling factor is given by where is the particle density per unit area. When we include spin, the degeneracy of each Landau level is doubled. If we assume the electrons are polarized by the magnetic field and the LL’s are filled sequentially by spin up and spin down electrons, then corresponds to states in the lowest Landau level (LLL) with LL index and corresponds to the second Landau level (SLL) with LL index . Thus, the 1/2-filled LLL corresponds to or 3/2 and the 1/2-filled SLL corresponds to or 7/2. The identifications of the 1/3-filled filling factors is obtained similarly.

The presence of electron-electron interactions clearly complicates this picture, but if we assume the electrons are confined to a single LL (i.e., there is no LL mixing), then the kinetic energy is a constant that can be ignored, giving the effective Hamiltonian:

(2) |

where is the total number of particles and is the distance between particles measured in units of the magnetic length . Note that our assumption of no LL mixing might not be a very good approximation Bishara and Nayak (2009); Rezayi and Simon (2011); Wójs et al. (2010) in all cases especially when considering the FQH state at where it has been observed at “low” magnetic fields where the LL mixing parameter ( characterizes the strength of the Coulomb interaction where is the electron charge and is the dielectric constant of the host semiconductor). Furthermore, the Hamiltonian we are using is completely particle-hole symmetric and apparently the ground state does not spontaneously break particle-hole symmetry Peterson et al. (2008c). This means that all of our conclusions about the MR Pfaffian state apply equally to the anti-Pfaffian state since the two ansatz are particle-hole conjugates. LL mixing, as shown recently Bishara and Nayak (2009); Rezayi and Simon (2011); Wójs et al. (2010), can induce three-body terms which explicitly break particle-hole symmetry leading to a possible preference toward the MR Pfaffian or anti-Pfaffian, however, in this work we do not consider such terms and can make no distinction.

This Hamiltonian (Eq. (2)) can be parametrized in terms of the relative angular momentum between two particles by the Haldane “pseudopotential” expansionHaldane (1983):

(3) |

where projects onto states with relative angular momentum and are the Haldane pseudopotentials for a given relative angular momentum and Landau level index . These pseudopotentials are the energies of a pair of particles with relative angular momentum confined to the -th LL. For spinless fermions, only for odd enter the Hamiltonian due to the Pauli exclusion principle. In a planar geometry, the Haldane pseudopotentials are given in terms of the Fourier transform of the interaction potential by

(4) |

where are Laguerre polynomials. This expansion allows us to work strictly in the Hilbert space of the lowest Landau level since all necessary information about electrons confined in higher Landau levels are contained in the pseudopotentials.

The model wavefunctions (Laughlin and MR Pfaffian states) can be obtained by diagonalizing certain “hardcore” Hamiltonians. The Laughlin wavefunction

(5) |

where is the electron coordinate in the complex plane, and is the zero-energy ground state of a special case of Eq. (3). For filling fraction , odd, this Hamiltonian is given by:

(6) |

This two-body “hardcore” potential penalizes any state where two particles have a relative angular momentum smaller than . The MR Pfaffian wavefunction

(7) |

is the exact zero-energy ground state of a three body Hamiltonian projecting onto electron triplets instead of pairs Greiter et al. (1991).

In an ideal, strictly two dimensional system, the electron-electron interaction is given by the 2D Coulomb potential, . The finite extent of an experimental quantum Hall system in the perpendicular direction will alter the ideal 2D interaction, yielding an effective quasi-2D electron-electron interaction. There are several models for the effect of the finite layer thicknessAndo et al. (1982); Stern and Das Sarma (1984); Sarma and Mason (1985); Zhang and Das Sarma (1986), however, these models all provide similar qualitative results Peterson et al. (2008a, b). Therefore, we will focus on one particular model, the infinite square well potential. In this model, we average the three dimensional Coulomb potential over the single-particle ground state of an infinite square well (i.e., ) in the perpendicular dimension, yielding the effective interaction potential

(8) | |||||

Combining Eq. (8) with Eq. (4) gives us effective pseudopotentials as a function of finite thickness , Landau level , and relative angular momentum .

We diagonalize the FQH Hamiltonians (one for each and LL index ) in the spherical geometryHaldane (1983) where electrons are confined to the surface of a sphere. Although we use this geometry, we use the pseudopotentials obtained from the infinite planar geometry (Eq.(4)) since the finite layer thickness effect is more conveniently modeled in this case. Furthermore, the pseudopotentials in the spherical geometry equal those in the planar geometry as the thermodynamic limit is approached (as the spherical radius is taken to infinity) and it can be argued that they provide a better approximation to the thermodynamic limit (this is discussed in detail in Ref. Peterson et al., 2008b and in the appendix). In the spherical geometry, the perpendicular (radial) magnetic field is provided by a magnetic monopole of strength, , quantized in half-integer units, placed at the center of the sphere. The eigenvalues of the squared magnitude and z-component of the angular momentum, and respectively, are good quantum numbers for the single particle wavefunctions, where is related to the LL index by the constraint, and is constrained such that . Thus the degeneracy of a LL with index is given by . The filling factor for a LL is defined in the thermodynamic limit . The uniform ground state has total angular momentum, and therefore, can be obtained in the Hilbert subspace where the total -component of angular momentum, .

In this study, we consider the FQH ground states at the Laughlin filling fractions and with particle number , 7, and 8 and the even-denominator filling fractions and with particle number and . We restrict ourselves to these relatively modest system sizes in order to investigate a large number of FQH ground states for various values of the finite thickness with reasonable computing resources. Although the Hilbert space for particle number at the half fillings is not prohibitively large, this system is also aliased with and could, therefore, yield ambiguous results. Since we are largely concerned with the qualitative features of the finite-thickness effect, these system sizes are adequate.

We calculate the entanglement entropy (EE) and the entanglement spectrum (ES) of FQH ground states by dividing the sphere into two regions. In particular we write our Fock space as the tensor product of two subspaces with respective basis states and where is the occupation number of the Landau orbital with angular momentum, and is the total number of Landau level orbitals. For all cases we choose our partitions when dividing our Fock space into and such that for the number of single particle orbitals even, and for odd, . Geometrically this is equivalent to dividing the sphere along a line of latitude (see Fig. 1).

We calculate three quantities of interest:

(i) The bipartite entanglement entropy EE of the ground state is given by where is the reduced density matrix defined by .

(ii) The entanglement spectrum is obtained from the eigenvalues of the reduced density matrix, , by the simple relation . This is equivalent to finding the Schmidt-decomposition of the matrix where , and . Since the quantum numbers for angular momentum and particle number in each region, , , are constrained such that and , the reduced density matrix is block diagonal with and being good quantum numbers for the eigenstates of . Therefore and are good labels for the corresponding entanglement spectrum.

(iii) The “conformal limit” of the entanglement spectrum. Recently, Thomale et al. Thomale et al. (2010) have introduced the conformal limit when calculating the entanglement spectrum of spherical FQH ground states. This limit is obtained by expressing the ground state in terms of a special choice of unnormalized basis states. The normalized single particle wavefunction on a sphere with angular momentum is given byHaldane (1983); Fano et al. (1986)

where and with and the usual spherical coordinates. In the conformal limit, we “unnormalize” the single particle wavefunctions by removing the prefactor in Eq. (II) such that the wavefunctions take the simple form . This procedure is an attempt to remove the finite size effects inherent in these calculations by basically removing the “length” in the problem. With the ground state redefined in this new basis, the entanglement spectrum is calculated as described above.

## Iii Results

### iii.1 Entanglement Entropy

We now report numerical results for the entanglement entropy (EE) of quasi-2D FQH ground states as a function of the finite layer thickness for FQH states in the LLL ( and 1/2) and the SLL ( and ). As mentioned above we choose the partition to be as close to the equator of the sphere as possible to minimize finite size effects.

The results for EE for the Coulomb ground state at filling fractions and are shown in Fig. 2 as a function of finite layer thickness . For comparison, the EE of the corresponding Laughlin model wavefunction is also shown as a independent horizontal line. In each of the figures, we see that in the LLL, the EE is near that of the Laughlin model wavefunction at and rises slightly as a function of . In contrast, the EE in the SLL is large compared to that of the Laughlin model at , but decreases as a function of and evidently reaches an asymptotic value. The qualitative behavior is independent of system size. If we consider for both the LLL and SLL filling fractions and speculate that is a qualitative measure of how far removed the ground state is from the Laughlin model state, then we see that the LL dependence of as a function of behaves qualitatively similar to that of the overlap between the ground state and the model wavefunction as reported in Ref. MacDonald and Aers, 1984; Zhang and Das Sarma, 1986; Peterson et al., 2008a, b. In particular, the ground state in the LLL is a “strong” FQHE state (i.e., is small) at and gradually becomes “weaker” for increasing (albeit only slightly), whereas in the SLL, the ground state is initially weak at but gets stronger with increasing (i.e., decreases). Thus, the EE for these cases qualitatively and semi-quantitatively captures how well the states are “Laughlin-like” as a function of in similar manner to the overlap.

Our operational definition of “weak” and “strong” depends on how close the EE of the Coulomb state is to the model state which, in this case, is the Laughlin state. In the SLL, becomes closer to the for the Laughlin state but, as mentioned above, appears to saturate at some asymptotic value that is still nearly . In contrast, the EE in the LLL is almost identical to that of the Laughlin state. We conjecture as to the reason for this difference between the EE in the SLL Coulomb ground state compared to the Laughlin state and the difference between the EE in the LLL as compared to the SLL: (i) it is possible that the FQHE at 7/3 is not described by the Laughlin state and is instead described by a state in a different topological universality class such as those given by Read and Rezayi Read and Rezayi (1999) and Bonderson and Slingerland Bonderson and Slingerland (2008), (ii) perhaps composite fermion interactions, which are thought Töke et al. (2005) to be more relevant in higher LLs, are producing this difference in and the Laughlin state, (iii) perhaps the 7/3 FQHE state is in fact a Laughlin state but our model system is leaving out realistic effects such as LL mixing which are crucial to its success.

Fig. 3 gives results for the EE of FQH ground states with even denominator filling fractions and as a function of finite layer thickness . Also shown in the figure is the EE of the Moore-Read Pfaffian state for comparison. In the LLL (), the EE has a weak minima as a function of , in contrast to the Laughlin fractions (this minimum is difficult to discern on our scale). The location of this minima changes with , suggesting a finite size effect, but the qualitative behavior is similar in both cases. In the SLL (), the EE has a very pronounced minima that approaches the EE of the MR Pfaffian model state for and crosses it for . This suggest that the FQH states becomes more MR Pfaffian-like at near an optimal . However, this optimal also changes with . Similarly to the Laughlin fractions, this LL dependence in EE as a function of is also qualitatively similar to that seen in the overlap between the FQH ground states and MR Pfaffian reported in Ref. Peterson et al., 2008a, b. These results suggest that is not particularly well-described by the MR Pfaffian, whereas is better described by the MR Pfaffian model state at finite thickness.

We note that recently, entanglement entropy in the SLL including finite thickness effects has been investigated Friedman and Levine (2008, 2010). However, the previous study did calculations using the torus geometry, in contrast to our spherical geometry, and attempted to isolate and calculate the so-called topological term of the entanglement entropy. The entanglement entropy can be essentially divided into two pieces

(9) |

where is the linear length of the boundary dividing the system into parts A and B (in our case it would be circumference of the sphere where we made our cut). The term is non-topological in origin and the term is the topological entropy and for the Laughlin and MR Pfaffian state can be calculated analytically Zozulya et al. (2007): for the Laughlin state and for the MR Pfaffian state at . Extracting the topological entropy from a Coulomb Hamiltonian requires one to numerically calculate the exact ground state for many different systems sizes and system cuts and perform a thermodynamic extrapolation. This is a labor intensive procedure that inherently induces statistical errors. Such numerical extrapolation, without some strong theoretical guidance about the finite-size behavior of the system, is often unreliable in estimating quantities in the thermodynamic limit.

The conclusion of Refs. Friedman and Levine, 2008, 2010 was that the topological entropy of the ground states of the LLL or SLL Coulomb Hamiltonians was consistent with associated model states (we note, however, that in Ref. Friedman and Levine, 2008 it was concluded that was more consistent with the Read-Rezayi state Read and Rezayi (1999) instead of the Laughlin state). However, they also included finite thickness in the form of an infinite square well potential and, interestingly, found that there was not much difference between the EE and the topological entropy with or without finite thickness included. We, however, clearly see a finite thickness effect on the EE. It is possible that this difference in the two studies (our present study and Refs. Friedman and Levine, 2008, 2010) is due to the different geometry used in the calculations (sphere vs. torus) but we find this scenario unlikely since most quantities of interest produce consistent results in the two geometries Rezayi and Haldane (2000); Peterson et al. (2008a, b). Such a comparison between geometries (torus vs sphere) was shown in Ref. Läuchli et al., 2010b to give similar results for the entanglement spectra of Laughlin states, supporting our suspicion. Moreover results given in Ref. Läuchli et al., 2010a suggests that the the extrapolation procedure performed in Refs. Friedman and Levine, 2008, 2010 may have been inappropriate for the torus. More work is clearly necessary to understand the difference between the results in spherical and toroidal geometries, particularly in the presence of the realistic finite thickness effects.

Before moving on to entanglement spectra we briefly discuss how our results compare to the previous overlap calculations done in Refs. Peterson et al., 2008a, b. The right panels in Figs. 2 and 3 gives one minus the percentage error in the entanglement entropy, . In Ref. Peterson et al., 2008b it is found that the overlap between the Laughlin state and the Coulomb ground state at 1/3-filling in the LLL and SLL is approximately at and is reduced monotonically to at in the LLL and is at and has a maximum of for in the SLL. These overlap trends are very consistent with what we have seen previously in EE. For the 1/2-filled LLL and SLL we find Peterson et al. (2008a, b) the overlap is relatively constant in the LLL at and in the SLL it is at and has a maximum value of nearly at . Again, one minus the percentage error in the entanglement entropy tracks the behavior in the overlap to a remarkable degree. Perhaps this is not a surprise since if the overlap is close to one then the EE (which is a particular combination of ) should also be nearly identical to the EE of the model state .

### iii.2 Entanglement Spectrum

In the previous section, we saw that the entanglement entropy (and in particular, ) as a function of behaves qualitatively similarly to the overlap Peterson et al. (2008a, b). For the half-filled case, increasing makes the calculated closer to the MR Pfaffian state for the SLL () in a rather dramatic way while increasing has very little effect on the in the LLL (), i.e., using the entanglement entropy as a measure we see that the MR Pfaffian is stabilized by finite thickness. For the 1/3-filled case we find that increasing drives away from the Laughlin value in the LLL () and closer to the Laughlin value in the SLL (), however, as in the previous overlap investigations, the value of the entanglement entropy for the case never gets as close to the Laughlin value as the 5/2 entanglement entropy gets to the MR Pfaffian. As discussed above, this could be a hint that something is missing from our understanding of the physics for the FQHE at .

To gain a deeper understanding of entanglement, we now turn our attention to the finite layer thickness dependence of the entanglement spectrum (ES), which as discussed earlier, provides more information than the EE alone. To calculate the ES, we partition the sphere the same as was done for the EE. We follow the convention established by Li and HaldaneLi and Haldane (2008) and restrict ourselves to the part of the ES where the number of particles in the partition, , is the same as that of the “root” configuration for the corresponding Laughlin or Moore-Read Pfaffian model wavefunctionBernevig and Haldane (2008); Li and Haldane (2008) for a given partition size . The “root” configurations describe the occupancy of LL orbitals for MR Pfaffian and Laughlin model states in the thermodynamic limit. Root configurations with a maximum -component of angular momentum, and their corresponding quantum numbers, and , are given in Table 1 for different filling fractions and partition sizes.

FQH state | root config. | |||
---|---|---|---|---|

Laughlin 1/3,7/3 | 8 | ‘10010010’ | 13.5 | 3 |

10 | ‘1001001001’ | 18 | 4 | |

11 | ‘10010010010’ | 24 | 4 | |

MR Pfaffian 1/2, 5/2 | 7 | ‘1100110’ | 16 | 4 |

9 | ‘110011001’ | 24.5 | 5 |

In order to obtain a general qualitative picture of how the ES changes as a function of the finite layer thickness, we calculate the “entanglement gaps” in each ES and plot it as a function of . An entanglement gapLi and Haldane (2008) is defined as the difference between the low-lying “CFT levels” (i.e., those levels displaying the low lying CFT counting structure) and the higher non-CFT levels for a given value of in the spectrum. According to the Li and Haldane conjecture, the state has an underlying CFT if the entanglement gaps are finite in the thermodynamic limit. However only the entanglement gaps at relatively “small” values of are relevant due to finite size effects, where is the total -component of angular momentum of the root configuration. The finite number of LL orbitals limits the number of possible “edge excitations.” Therefore only a few levels are expected to have the same counting structure as the CFT edge modes. The “depth” (i.e., the max ) at which the counting structure in the ES is consistent with the CFT edge modes is dependent on the system size, . It has been conjectured Hermanns et al. (2010) that these “finite size” effects of the entanglement spectra contain information such as the generalized statistics of the underlying FQH state, however, we will not consider such conjectures in this work.

We can determine a suitable cutoff for by examining when the level counting in the ES of the model states deviate from the expected counting in the thermodynamic limit. To illustrate this finite size cutoff, we give the ES of the Laughlin state in Fig. 4. For the Laughlin state, the multiplicity of CFT levels is given by where is the partition function of the integer . The first 7 values of , starting with are 1, 1, 2, 3, 5, 7, and 11. In Fig. 4, we see that for and 7, the level counting begins to deviate from at = 4, and for , the deviation begins at . Thus, for our study we will focus on the entanglement gaps for , 2 and 3 for the and 7 Laughlin systems, and for , we also examine the entanglement gap at .

We determine the finite size cutoff for the entanglement gaps of the half-filled FQH states in a similar manner, which we now illustrate. The ES for the MR Pfaffian model states are shown in Fig. 5. The counting rules for the MR Pfaffian model state depend on where the partition is made, which correspond to choosing one of the three sectors of the corresponding CFTLi and Haldane (2008). For the case of , the partition along the equator is equivalent to the partition in Li and Haldane’s nomenclature (i.e., a cut between two unoccupied orbitals in the root configuration). The CFT level counting for the first 4 Virasoro levels of this partition are 1, 1, 3, and 5. The counting in the MR Pfaffian ES given in Fig. 5 with deviates from this structure at . For , the partition along the equator corresponds to (i.e., a cut between two occupied orbitals of the root configuration), which has a CFT level counting of 1, 2, 4, and 7 for the first 4 Virasoro levels. Examining the ES in the figure for , we see this spectrum also deviates from the expected counting at . Thus for the half-filled FQH states we examine in this study, we concern ourselves only with the entanglement gaps up to .

The entanglement gaps, which we denote as for , are calculated by finding the difference between the largest suspected low-lying CFT level and the next highest level at the given value of in the ES of the numerically obtained Coulomb ground states for varying . We also calculate the minimal gap between CFT and non-CFT levels for , which we denote as where is the cutoff described above. The minimal gap gives us a qualitative measure of how well separated, overall, the CFT levels are from the generic non-CFT levels. The suspected low-lying CFT levels are identified by the expected counting described above. An example of this procedure is shown in Fig. 6, which shows the entanglement gaps in the ES of the ground state for and . Note that throughout this work, when presenting figures showing ES, we color code the suspected CFT levels with red (gray) diamonds connected by a red (gray) dash and all other ES levels with a black dash. The suspected CFT levels are chosen by calculating the ES for the model state (be it the Laughlin or the MR Pfaffian) and noting how many ES levels there are for each . Then, when we consider the ES for the Coulomb Hamiltonians, we identify suspected CFT levels (and color code them) as the lowest ES levels for each .

Entanglement gaps as a function of finite layer thickness for the Laughlin filling fraction are shown in Fig. 7. The entanglement gaps are slightly decreasing with for all cases, indicating that the states are weakening. These trends are similar to those observed in the EE at . Note that the minimal gap for and is initially small and becomes zero for . This may indicate that the FQH state collapses at a finite thickness, as has been shown in previous worksZhang and Das Sarma (1986); He et al. (1990) (the previous works showed the FQHE to collapse at very large ). However, this effect is not seen in the minimal entanglement gap for . This “even-odd” finite size effect is likely due in part to a trade-off between the finite size cutoff and the number of available orbitals. Indeed, the finite size cutoff is the same for and , but the larger Hilbert space for the case allows for more “edge excitations” that strengthen each entanglement gap, not just the minimal gap, compared to . In all cases, however, the overall trends in the entanglement gaps (i.e., slight decrease with ) are qualitatively similar to those seen in the EE (in particular, ) and the overlap in Refs. Peterson et al., 2008a, b.

To illustrate this overall trend in the entanglement gaps for , we provide the ES of the ground states in Fig. 8 for = 0, 2 and 6. We have marked the levels that are consistent with the counting found in the ES of the Laughlin model state shown in Fig. 4 for all values of and indicate our chosen finite size cutoff. We see that qualitatively, the ES is largely insensitive to finite . Moreover, on the right of the finite size cutoff, except for the largest CFT state at , the low-lying CFT levels are well-separated from the higher energy generic levels.

We now examine the case when in comparison. In Fig. 9 are the entanglement gaps as a function of finite layer thickness. For , the root entanglement gap, , generally increases with . and each have a weak, local maxima near and is actually decreasing with . Moreover the minimal entanglement gap is zero throughout. The entanglement gaps for are are each monotonically increasing with , similarly to in the case. The minimal gap, which is initially zero, opens at and then gradually increases with in this case. The case when shows trends similar to the case. Here , , and increase with , has a local maxima near , and decreases with . The minimal gap for is zero throughout. Again we see an “even-odd” finite size effect in the entanglement gaps as was seen with . However, in this case, we have entanglement gaps that increase, decrease, or have a weak maxima as a function of . This is in contrast to the case where all entanglement gaps follow the same trend with finite . The different trends in the entanglement gaps may suggest that the topological signature of the state differs from that of the Laughlin state.

Some illustrative examples of ES at are given in Fig. 10 with and = 0, 4 and 6. The given ES appear to have structure similar to that seen in the case, however, we see that for , the higher energy “CFT” states are virtually indistinguishable from the “generic states”. This “blending” appears to get worse for larger . Again, these results may suggest that the Laughlin model state is not an accurate description for the state.

Results on the entanglement gaps for the even denominator filling fraction are shown in Fig. 11. Here we see that the entanglement gaps slightly decrease with and behave similarly to the EE at this filling fraction. Also note that for , the minimal gap is small and decreases with , while for , the minimal gap is zero throughout. As mentioned earlier, there has been no definitive experimental observation of FQHE at in monolayer systems consistent with our calculations. The full ES of the ground states with are given in Fig. 12 for = 0, 2 and 6 and . Qualitatively, we see that the ES is largely insensitive to the finite-thickness effect. Moreover, the largest suspected CFT level for is well separated from the other CFT levels and appears to be more consistent with the generic levels. Again, this suggests that is not described by the MR Pfaffian wavefunction.

Fig. 13 shows the entanglement gaps at filling fraction . For each entanglement gap peaks at a certain value for . In particular peaks near ; the other gaps peak near . We also see peaks in the entanglement gaps for the case when . Here, the gaps gradually rise to a local maxima near and then slowly decay for increasing . Note that the gaps in this case are generally smaller compared to those observed for . These results may suggest that there is a slight difference in the finite-size effect on the different MR Pfaffian CFT sectors. However, these results are qualitatively similar to the EE results and the results on the overlap in Refs. Peterson et al., 2008b (i.e., the MR Pfaffian signature of the state is strengthened by the finite size effect).

We also provide the ES of the state for in Fig. 14 for , 4, and 6. Here, we see the ES “opens” at , giving a larger separation between the CFT and generic levels in the spectrum compared to = 0 and 6. Again, these results suggests the is, indeed, described by the MR Pfaffian wavefunction, and this description is more stable at finite thickness.

In summary, the entanglement gaps in the ES have a similar dependence on finite thickness as the EE, leading to similar conclusions. However finite size effects prevent us from making definitive statements. In the next section, we attempt to alleviate this problem using the conformal limit.

### iii.3 The Conformal Limit

In the previous section, we used entanglement gaps in the ES to evaluate the “strength” of a state as a function of and we were able to confirm the MR Pfaffian signature of the state and distinguish it from the (lack of) signature of the state. However, we have intentionally ignored a significant part of the ES in order to avoid finite size effects, i.e., we focused on the region of the ES with small (see Fig. 6). We determined the size of this region by examining where the ES of the finite sized MR Pfaffian and Laughlin model states deviate from the conjectured structure in the thermodynamic limit (i.e., the edge state level counting given by the suspected CFT). Thus, we have confirmed the MR Pfaffian and Laughlin signatures only to a certain extent because, in fact, there is not an actual entanglement gap in the ES.

It has been conjectured that the full entanglement spectrum of the finite sized model states contain information on the topological signature of the FQH state.Hermanns et al. (2010); Chandran et al. (2011) Thus, all states in the ES can be used to identify the topological quantities. With this in mind, we now examine the entanglement spectrum of quasi-2D FQH states in the “conformal limit” (CL), which reportedly allows the use of the entire spectrum to examine the state by unambiguously defining a full entanglement gap. As discussed briefly above, and at length by Thomale et al. in Ref. Thomale et al., 2010, the CL works by removing finite size effects due to the curvature of the sphere and gives an ES with a “full” unambiguous entanglement gap in the spectrum for topologically ordered states. Thus the presence of an entanglement gap in the conformal limit is conjectured to be a sign of topological order. A demonstration of an ES before and after the CL is given in Fig. 15. After taking the CL of an ES (CLES), we determine the “minimal gap” by taking the difference between the highest suspected CFT level and the lowest generic level in the entire spectrum. The suspected CFT levels are determined by comparing the CLES to that of the model state with the assumption that all levels in the ES of the model state are CFT levels. For comparison, we examine the entanglement gaps for each value of and define the “average gap” as the average of the individual entanglement gaps. We also define the “maximum gap” as the maximum of the entanglement gaps. Individual gaps that are near infinite (i.e., no levels above the highest CFT level) are ignored. The minimal gap, the average gap, and the maximum gap are calculated for each CLES as a function of the finite layer thickness, .

CLES entanglement gaps as a function of finite layer thickness for (LLL) are shown in Fig. 16. For and 7, the entanglement gap measures decreases with but remains finite throughout. This behavior is qualitatively similar to the ES gaps for small , as well as the EE results, suggesting a weakening of the Laughlin state. The fact that the minimal entanglement gap in the ES for the case (Fig. 7) differs from the minimal gap in the CLES may indicate that the closing of the gap in the ES is due to finite size effects related to the curvature of the geometry rather than the limited number of LL orbitals. However, the minimal gap for the case where seems anomalous. Although the average and maximum gaps follow similar qualitative trends, the minimal gap is at or near 0 for all values of , including . How to interpret this result is unclear since there is a general consensus that the Laughlin state does, indeed, model the state. We can shed some light on this anomaly by examining the CLES of the FQH states directly. Fig. 17 shows the CLES for the , FQH state at finite thicknesses , 2 and 6. The suspected CFT levels are marked in each plot.

We examine CLES in the SLL case () in Fig. 18. In general, the behavior of each gap measure differs with varying . The minimal gap appears fragile and virtually disappears for larger . The average gap has two local maxima in for . Only one of the local maxima in the average gap is preserved when we look at the case, and for , the average gap fluctuates. The maximum gap, in general, increases with increasing but has a notable peak near for . The inconsistency in these results may suggest, from the ES and EE results, that the Laughlin model state is not a suitable model for the state, or other ignored effects are needed for the Coulomb state to be adequately described by the Laughlin state.

Fig. 19 shows the CLES results for , at finite thickness = 0, 4, and 6. We note that for each value of , there is very little separation between the suspected CFT levels and the generic levels. Indeed, if the suspected CFT levels were not marked, there is no clear entanglement gap across the whole spectrum. However, there does appear to be structure in the CLES for small values of (i.e., near the “root” configuration). What this may imply about the topological signature of the state is not clear.

Results for the CLES gap measures in the even denominator FQH state are shown in Fig. 20. We see that the minimal gap is nonzero and gradually increases with for . However, for the minimal gap is zero throughout. The maximum and average gaps decrease with for . For , the average gap has several local maxima, while the maximum gap decreases then suddenly becomes constant with with two sharp peaks. Given our results on the EE and the ES for this state, the inconsistency between the and in the CLES gap measures may suggest that the MR Pfaffian model state is not a suitable model for . We also provide the CLES of the ground states in Fig. 21 for and , 4, and 6. Qualitatively, the CLES do not change very much as a function of , and there is no clear separation between the CFT and generic levels. This, again, suggests that there is no FQH state at this filling fraction.

The CLES gap measures for in the SLL are given in Fig. 22. For , the minimal gap has a very pronounced peak near . The average and maximum gaps, however, have local minima near where the minimal gap is maximum. These “cusps” are a result of level crossings. For , the minimal gap is initially zero, but becomes finite for non-zero and peaks near . The average and maximum gap in this case have similar shapes with a peak near . These results are qualitatively similar to the results of the ES, EE, and the overlap in Refs. Peterson et al., 2008a, b. Moreover, the difference between and may suggest that finite thickness affects each sector of the suspected CFT differently, but larger system sizes are necessary to verify this. In the CLES plots shown in Fig. 23 for and , 4 and 6 respectively, we see the entanglement gap between CFT and generic levels “open” at finite compared to and 6. These results are consistent with results observed with the EE and the ES, indicating that MR Pfaffian signature strengthens with a finite .

In summary, taking the conformal limit of the entanglement spectra provides us with a full entanglement gap in most cases with a finite thickness dependence that is qualitatively similar to the results on the EE. The notable exceptions are the which has little or no entanglement gap consistent with experimental observations, the case at that is not expected given results with other system sizes, and the case that is consistent with the results on the EE and ES suggesting that other physics besides the Laughlin state alone is needed to explain this FQHE. The case with and , however, is inconsistent with most theory and experiment, but when we examine the spectra directly, there are a few “spurious” states that cross an otherwise full gap. The origin of these “spurious” states are related to our use of planar Haldane pseudopotentials rather than spherical pseudopotentials and is discussed in the appendix. However we do not expect this choice to alter the topological features of the state. Therefore, this result may suggest that a “full” quantitative entanglement gap is not necessary to identify a topological state. In the next section, we introduce the concept of a entanglement spectral density of states where a qualitative, “soft” gap may be identified in such cases.

### iii.4 Entanglement Spectrum Density of States

In the entanglement results presented above, we require a model state wavefunction for comparison in order to systematically define and calculate the entanglement gaps. These methods have the obvious disadvantage of requiring an ansatz for comparison. In the conformal limit case, we assume the low-lying states in the entanglement spectrum should have the exact counting as seen in the model entanglement spectrum. This assumption may be premature since other finite size effects may cause the counting to deviate, even after taking the conformal limit especially at the largest , see Fig. 17. With this in mind, we attempt to obtain a general qualitative sense for how the entanglement spectra vary with finite layer thickness by extending the analogy with “energy levels” a bit further by calculating the “density of entanglement spectral states.” With the density of states, we can qualitatively look for entanglement gaps without relying on a model state for comparison. Also, we may be able to detect “soft” gaps where a small number of states may be present within an otherwise prominent gap between two peaks in the density of states. Thus in this section we briefly examine this extension by providing results for the density of states (DOS) of the entanglement spectrum, both with and without the conformal limit, as a function of finite layer thickness, .

The plots shown in Fig. 24 give the density of states of the ES and CLES for as a function of finite layer thickness for . In the DOS for the ES before taking the CL, we see sparse low lying states that are separated from a denser cloud of higher states by a series of gaps. These low lying states are the CFT states from the Li and Haldane conjecture. The states appear, largely to be insensitive to the finite layer thickness. Turning to the DOS for the CLES, a clear gap is much more evident for the cases. Here the higher-energy states appear to makeup a wide, low-density band that is well-separated from a low, dense band of states by a gap that decreases with . This case seems to illustrate the effect of using the CL.

For comparison, we provide the DOS results for in Fig. 25. Here there also appears to be low-lying CFT states in the ES below a high density region of higher energy states. In the DOS of the CLES, a “clear” gap does not appear. But the low lying band in this case does appear qualitatively similar to the case. One may possibly associate a “soft” gap in this case, where a few states appear to be present between two somewhat distinct regions in the DOS. This “soft” gap is qualitatively similar to the “clear” gap in the case and it does appear to slightly decrease as a function of . However, it is difficult to distinguish this “soft” gap from the other small gaps in the spectrum.

In Fig. 26 we provide DOS plots for the FQH state for as a function of . In this case, the ES is especially sensitive to finite layer thickness. However, we still see a series of small gaps separating thin, dense bands at lower energies. After taking the CL, a clear gap at finite (non-zero) thickness has a definite peak corresponding to a level crossing. Below the gap, there appears to be some band crossings as is varied.

In summary, the DOS of the entanglement spectra (with and without the conformal limit) gives us a general qualitative picture of how the ES evolve with a varying parameter (i.e., the finite layer thickness in our case). Thus we expect the DOS of the ES to be a good initial cursory tool in examining topological states with varying parameters.

## Iv Conclusions

In this work we study entanglement in finite sized, quasi-2D FQH states via the entanglement entropy and the entanglement spectrum as a function of the finite layer thickness of the transverse dimension in a realistic FQH system and compare them to the entanglement signatures of the Laughlin and MR Pfaffian model states. For the Laughlin filling fractions, we find that the EE increases (decreases) with finite layer thickness for () in the LLL (SLL) with increasing (decreasing) deviation from the EE of the Laughlin model state. However the EE in the SLL reaches an asymptotic value larger than the EE of the Laughlin state, possibly suggesting the state is modeled by different physics than the Laughlin state. Similar behavior is also seen in the entanglement gaps of the ES for the Laughlin filling fractions. Here we find that the entanglement gaps decrease with finite layer thickness for the Laughlin filling fractions in the LLL. But in the SLL, the behavior of the entanglement gaps depend on the “depth” of the gap. These results suggest that the Laughlin FQH states “weaken” with increasing thickness in the LLL, which is consistent with previous work on quasi-2D FQH states Peterson et al. (2008a, b), but in the SLL, other physics beyond just the Laughlin state alone is needed to describe the FQH state. The LL dependence of the finite thickness effect at half-filling differs slightly. The EE of the state in the LLL is largely insensitive to the finite layer thickness in contrast to that of the SLL state where the EE has a local minima that approaches the EE of the MR Pfaffian at finite . This qualitative behavior is also seen in the entanglement gaps of the ES for half-filled LLs. For in the SLL, we see peaks (local maxima) in the entanglement gaps at finite thickness, suggesting the is more “MR Pfaffian-like” at an optimal thickness, which, again is consistent with previous workPeterson et al. (2008a, b) and strongly suggests the state is, indeed, MR Pfaffian. In contrast the entanglement gaps of the ES for the state suggest that it is not modeled by the MR Pfaffian. Thus, the entanglement gaps in the ES allows us to differentiate the and states, which we could not definitively establish with the EE or the overlap calculations. Of course, we must be cautious with these results since the calculated entanglement gaps made use of only a few “Virasoro levels” in the low-lying CFT due to finite-size effects. Assuming the Li and Haldane conjecture to be correct, we can only say that we have observed the Laughlin and MR Pfaffian signatures up to a few “Virasoro levels”.

We also investigate the conformal limit of the entanglement spectrum which is conjectured to remove curvature in the spectrum due to finite size effects and allow the use of the entire spectrum to determine the topological signature of the state. Our results on the conformal limit, however are inconsistent between varying system sizes and are difficult to interpret. This appears to be due to our choice of using planar pseudopotentials rather than spherical pseudopotentials in obtaining the FQHE ground states. In the appendix we examine this choice by comparing the entanglement spectra of ground states obtained by using either spherical or planar pseudopotentials at and observe that the conformal limit can be affected by components of the ground state that have exponentially small contributions and, therefore, are sensitive to minor details in the interaction (such as the difference between planar and spherical pseudopotentials). Thus the presence of the entanglement gap in the conformal limit is sensitive to certain details in the effective interaction that may not be relevant in determining the topological features of the state. Further work using much larger system sizes would be necessary to resolve this issue which is well beyond the scope of the current work.

We have also introduced the notion of entanglement density of states as a method for examining the idea of an entanglement gap without an explicit reliance on a model wavefunction. Although, far from definitive, the entanglement DOS suggests itself as a powerful tool to determine the topological nature of a particular ground state. Our detailed numerical study establishes the entanglement DOS to be a useful quantity underlying topological FQHE particularly in the context of finite size numerical calculations.

It is interesting to observe that the entanglement measures give results similar to those obtained with overlaps in Refs. Peterson et al., 2008a, b. Whereas the overlap is a simple measure of how well a numerically obtained ground state matches a particular model state (e.g. the Laughlin state or the MR Pfaffian state), the entanglement measures (in particular, the ES) is a more general measure of how well a state fits a suspected CFT (i.e. universality class) that describes the model state. Therefore, it can be said that these results confirm the conclusions in Refs. Peterson et al., 2008a, b in a more general sense in respect to the Laughlin and MR Pfaffian CFTs. However, we must be cautious in this generalization given that we have only observed the Laughlin and MR Pfaffian signature up to a few “Virasoro levels” and different theories can result in the same low-level structure in the ES Turner et al. (2009). More work is necessary to understand how well entanglement measures can definitely identify universality classes in finite systems.

In interpreting our results and conclusions, one may wonder about the importance of finite size effects on our numerical diagonalization. The possible limitations associated with finite-size effects are of course always present in any exact diagonalization study of any FQHE system, and the possibility that some of the conclusions are affected by finite size effects can never be ruled out even if the calculations are carried out on systems much larger than what we use in this work, since in the end any statement about an experimental system based on calculations performed on few-particle systems is always subject to an extrapolation to the thermodynamic limit. We believe that all our conclusions regarding the importance of finite quasi-2D thickness effect on the FQHE entanglement spectra are valid independent of the rather modest size of our finite system diagonalization study because earlier work Peterson et al. (2008a, b) clearly established, when compared with calculations Storni et al. (2010); Feiguin et al. (2008) carried out on much larger systems, that the system size we use in this work, namely N=8, is certainly adequate in making qualitatively correct conclusions about the SLL FQHE. Our goal in this paper has been to study as many FQHE states as feasible as a function of the quasi-2D layer thickness in depth, thus necessarily (due to the computational time restrictions) limiting our system size to N=8 which should be adequate. Nevertheless, we feel that future work should explore larger system size diagonalization in order to study the finite-thickness effect on the entanglement spectra of various FQHE states.

We also note that there are various alternatives to the infinite square well effective potential in examining the finite thickness effect. The Zhang-Das Sarma potential is likely the most well-known and oft-used alternative in this regard. We choose to focus on the infinite square well instead of the Zhang-Das Sarma potential because the infinite square well is more closely aligned to realistic quasi-2D systems. Moreover many observables calculated with the Zhang-Das Sarma potential besides the overlap have been shown to give similar results compared to the infinite square well. Therefore we expect the two potentials to give similar results with entanglement measures as well, but since investigating the nature of these potentials is not a goal of the present study, we leave explicit verification of this assertion for a future work.

It should be noted that the concept of entanglement spectra (or for that matter, the entanglement entropy itself) has no direct experimental or observational implications since it cannot be directly measured. The concept is useful conceptually and theoretically in ascertaining the quantum topological nature of a particular interacting Hamiltonian, and in that sense its experimental consequences are indirect since the topological nature of a system has obvious experimental consequences. We have investigated in this work the utility of the concept of the entanglement spectra in ascertaining the underlying topological nature of realistic FQHE states as a function of the quasi-2D layer thickness, finding that the entanglement spectrum provides results consistent with what has earlier been established in the literature based on the wavefunction overlap studies. Our very detailed study also indicates that at this stage of theoretical development, the entanglement approach is perhaps no more predictive in providing experimental implications of various FQHE states than what is already available in the literature based on the direct wavefunction overlap studies. Further work would be necessary to see if the entanglement approach has some specific advantages in predicting experimental properties of FQHE states not already apparent in wavefunction-based analyses.

In conclusion, we have extended the concept of topological entanglement spectra and entanglement gaps to finite-thickness FQH systems by calculating the FQHE topological properties systemically as a function of finite thickness of the quasi-2D systems, establishing in the process that the FQHE entanglement measures calculated as a function of system thickness are completely consistent with the results obtained earlier in the literature using wavefunction overlap calculations. While our work establishes various entanglement measures as important theoretical quantities classifying FQHE, more work will be necessary to understand the finite size aspects of entanglement spectra and entanglement gaps in the context of realistic fractional quantum Hall systems. Although it is gratifying that the qualitative conclusions of our entanglement-measure-based results in this work are completely consistent with earlier FQHE results obtained on the basis of wavefunction overlap calculations, it remains to be seen whether the entanglement-measure based probes have more predictive power regarding the nature of FQHE than the wavefunction-overlap based probes or it is simply a deeper way of looking at the same physics with no obvious additional implications for the experimental occurrence of FQHE.

## Appendix A Planar vs spherical pseudopotentials at

The analysis presented above is based on ground state wavefunctions obtained by diagonalizing the quasi-2D Coulomb potential in a spherical geometry. However the Haldane pseudopotentials used to construct the Hamiltonian are derived from a infinite planar geometry rather than a spherical geometry. We choose to use planar rather than spherical psuedopotentials because (i) the effective Coulomb potential in a quasi-2D system is more naturally obtained in the infinite planar geometry and (ii) we expect the spherical and planar pseudopotentials to be indistinguishable in the thermodynamic limit. Moreover given the mostly qualitative nature involved in studying entanglement spectra, we expected this choice to make little difference in the results. Nevertheless, there are cases under study where this choice matters. The goal of this appendix is to highlight some of these cases. We show that for , the low energy spectrum in the entanglement spectra are qualitatively similar between ground states obtained from either spherical or planar pseudopotentials, but higher energy spectra can can differ in some cases. This difference does not change the qualitative conclusions drawn from the low energy spectra, but when we consider the conformal limit which looks for a full entanglement gap, the difference can lead to different conclusions (in particular, the case of with ). We leave the comparison of cases with and larger for future work.

6 | 1/3 | 0.9988 | 0.9964 | 0.9921 |
---|---|---|---|---|

6 | 7/3 | 0.9480 | 0.5285 | 0.7369 |

7 | 1/3 | 0.9999 | 0.9964 | 0.9952 |

7 | 7/3 | 0.8648 | 0.6071 | 0.8737 |

8 | 1/3 | 0.9996 | 0.9954 | 0.9954 |

8 | 7/3 | 0.9675 | 0.5719 | 0.7441 |

8 | 1/2 | 0.9978 | 0.9213 | 0.8953 |

8 | 5/2 | 0.9688 | 0.8674 | 0.9639 |

10 | 1/2 | 0 | 0.8891 | 0 |

10 | 5/2 | 0.9720 | 0.8376 | 0.9342 |

In Table 2 we provide several overlap calculations between exact ground states at obtained using either spherical or planar pseudopotentials. In column 3 of the table, we see that the overlap between the ground states from the spherical and planar cases is generally high. The notable exception is the case when and where the overlap is . In this case the ground state obtained with the planar pseudopotentials possesses a different symmetry compared to the ground state of the spherical case, which leads to a vanishing overlap. Excluding these, columns 4 and 5 of the table show that the overlap between the spherical and planar ground states with the model Laughlin or MR Pfaffian states are qualitatively similar.

We now turn our attention to the entanglement spectra and how they may differ with choice of pseudopotentials. ES (with and without the conformal limit) for the exact ground state of the FQHE state at and using spherical and planar pseudopotentials are given in fig. 27. In the figure, we see that the ES with planar pseudopotentials (Fig. 27 a) is qualitatively similar to the spectra obtained with spherical pseudopotentials (Fig. 27 b). The same can also be said with the ES in the conformal limit between the planar case (Fig. 27 c) and the spherical case (Fig. 27 d). Thus, given the results in Fig. 27, we would expect that the choice of pseudopotenials makes little difference in obtaining a qualitative understanding of the ES in this case.

Fig. 28 compares the ES of the FQHE state at filling fraction with . In this case we see that in the ES before the conformal limit (Fig. 28 a and b), the low energy spectra are qualitatively similar between the planar and spherical cases. The higher energy spectra in the ES, however, show notable differences with the planar case having a few CFT levels at much higher energy compared to the spherical case. In the conformal limit, these higher energy CFT levels lead to a vanishing entanglement gap in the conformal limit for the planar case (Fig. 28 c) compared to the spherical case (Fig. 28 d) where there is a full entanglement gap. These are the same “spurious” levels identified earlier in section III.3. These results suggest that the vanishing minimal gap seen in Fig. 16 is due to our choice of planar rather than spherical Haldane pseudopotentials. This may seem surprising given the large overlaps seen in Table 2. However the states associated with the the higher energy CFT levels have exponentially small contributions to the ground state wavefunction, and thus contribute little to the overlap. We might also expect these states to be more sensitive to certain quantitative details of the potential that do not affect the qualitative picture of the FQHE ground state (e.g. values of for “large” ). Thus when taking the conformal limit, the choice of pseudopotential may matter in some cases in order to observe a full entanglement gap. But a qualitative understanding can still be gleaned from the planar case since there does appear to be two distinct regions in the CLES that we can identify, at least qualitatively, as CFT and non-CFT levels.

We now compare the spherical and planar pseudopotentials in the SLL with . Fig. 29 gives the ES (29 a and b) and CLES (29 c and d) for the FQHE state obtained with either planar or spherical pseudopotentials with . The planar and spherical cases are qualitatively similar in both the ES and CLES and both suggest that the Laughlin wavefunction may not describe this state, as discussed in sections III.1, III.2 and III.3.

Results for the even denominator filling fraction with are given in Fig. 30. In this case, the planar results (28 a and c) differ considerably from that of the spherical case (28 b and d). This is not surprising since the overlap between these two states given in table 2 vanishes. However, it appears that neither state is consistent with the MR Pfaffian.

Comparison of FQHE ground states obtained with planar and spherical pseudopotentials for the state with is given in Fig. 31. Similar to the case, the low energy spectra in the ES (31 a and b) are qualitatively similar between the two cases. The higher energy levels in the spectra do differ, but the CLES (31 c and d) does appear to give the same qualitative picture. Recall that in section III.3, the minimal gap for this case becomes non-zero only at finite for the planar case. We would expect a similar result to occur using the spherical psuedopotentials. Verification of this is left for a future work.

###### Acknowledgements.

This work is supported by Microsoft Q and DARPA QUEST.## 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).
- Halperin (1984) B. I. Halperin, Phys. Rev. Lett. 52, 1583 (1984).
- Leinaas and Myrheim (1977) J. Leinaas and J. Myrheim, Il Nuovo Cimento B (1971-1996) 37, 1 (1977), 10.1007/BF02727953.
- Wilczek (1982) F. Wilczek, Phys. Rev. Lett. 49, 957 (1982).
- Arovas et al. (1984) D. Arovas, J. R. Schrieffer, and F. Wilczek, Phys. Rev. Lett. 53, 722 (1984).
- Das Sarma and Pinczuk (1997) S. Das Sarma and A. Pinczuk, Perspectives in quantum Hall effects: novel quantum liquids in low-dimensional semiconductor structures (Wiley, 1997).
- Prange and Girvin (1990) R. Prange and S. Girvin, The Quantum Hall Effect (Springer-Verlag, New York, 1990).
- Jain (2007) J. Jain, Composite fermions (Cambridge University Press, 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).
- Jain (1989) J. K. Jain, Phys. Rev. Lett. 63, 199 (1989).
- Suen et al. (1992) Y. W. Suen, L. W. Engel, M. B. Santos, M. Shayegan, and D. C. Tsui, Phys. Rev. Lett. 68, 1379 (1992).
- Eisenstein et al. (1992) J. P. Eisenstein, G. S. Boebinger, L. N. Pfeiffer, K. W. West, and S. He, Phys. Rev. Lett. 68, 1383 (1992).
- He et al. (1993) S. He, S. Das Sarma, and X. C. Xie, Phys. Rev. B 47, 4394 (1993).
- Moore and Read (1991) G. Moore and N. Read, Nuclear Physics B 360, 362 (1991).
- Lee et al. (2007) S.-S. Lee, S. Ryu, C. Nayak, and M. P. A. Fisher, Phys. Rev. Lett. 99, 236807 (2007).
- Levin et al. (2007) M. Levin, B. I. Halperin, and B. Rosenow, Phys. Rev. Lett. 99, 236806 (2007).
- Nayak and Wilczek (1996) C. Nayak and F. Wilczek, Nuclear Physics B 479, 529 (1996).
- Bonderson et al. (2011) P. Bonderson, V. Gurarie, and C. Nayak, Phys. Rev. B 83, 075303 (2011).
- Das Sarma et al. (2005) S. Das Sarma, M. Freedman, and C. Nayak, Phys. Rev. Lett. 94, 166802 (2005).
- Nayak et al. (2008) C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008).
- Storni et al. (2010) M. Storni, R. H. Morf, and S. Das Sarma, Phys. Rev. Lett. 104, 076803 (2010).
- Storni and Morf (2011) M. Storni and R. H. Morf, ArXiv e-prints (2011), eprint 1101.5290.
- Wójs (2009) A. Wójs, Phys. Rev. B 80, 041104 (2009).
- Amico et al. (2008) L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. Phys. 80, 517 (2008).
- Venkatachalam et al. (2010) V. Venkatachalam, A. Yacoby, L. Pfeiffer, and K. West, Nature 469, 182 (2010).
- M. Dolev and Mahalu (2008) V. U. A. S. M. Dolev, M. Heiblum and D. Mahalu, Nature 452, 829 (2008).
- Radu et al. (2008) I. P. Radu, J. B. Miller, C. M. Marcus, M. A. Kastner, L. N. Pfeiffer, and K. W. West, Science 320, 899 (2008).
- Dolev et al. (2010) M. Dolev, Y. Gross, Y. C. Chung, M. Heiblum, V. Umansky, and D. Mahalu, Phys. Rev. B 81, 161303 (2010).
- Willett et al. (2009) R. L. Willett, L. N. Pfeiffer, and K. W. West, Proceedings of the National Academy of Sciences 106, 8853 (2009).
- Willett et al. (2009) R. L. Willett, L. N. Pfeiffer, and K. W. West, ArXiv e-prints (2009), eprint 0911.0345.
- Bid et al. (2010) A. Bid, N. Ofek, H. Inoue, M. Heiblum, C. Kane, V. Umansky, and D. Mahalu, Nature 466, 585 (2010).
- Milovanović and Read (1996) M. Milovanović and N. Read, Phys. Rev. B 53, 13559 (1996).
- Overbosch and Wen (2008) B. J. Overbosch and X.-G. Wen, ArXiv e-prints (2008), eprint 0804.2087.
- Rhone et al. (2011) T. D. Rhone, J. Yan, Y. Gallais, A. Pinczuk, L. Pfeiffer, and K. West, Phys. Rev. Lett. 106, 196805 (2011).
- Stern et al. (2010) M. Stern, P. Plochocka, V. Umansky, D. K. Maude, M. Potemski, and I. Bar-Joseph, Phys. Rev. Lett. 105, 096801 (2010).
- MacDonald and Aers (1984) A. H. MacDonald and G. C. Aers, Phys. Rev. B 29, 5976 (1984).
- Zhang and Das Sarma (1986) F. C. Zhang and S. Das Sarma, Phys. Rev. B 33, 2903 (1986).
- He et al. (1990) S. He, F. C. Zhang, X. C. Xie, and S. Das Sarma, Phys. Rev. B 42, 11376 (1990).
- Ortalano et al. (1997) M. W. Ortalano, S. He, and S. Das Sarma, Phys. Rev. B 55, 7702 (1997).
- Park and Jain (1998) K. Park and J. K. Jain, Phys. Rev. Lett. 81, 4200 (1998).
- 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).
- 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).
- Pan et al. (1999) W. Pan, J.-S. Xia, V. Shvarts, D. E. Adams, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer, K. W. Baldwin, and K. W. West, Phys. Rev. Lett. 83, 3530 (1999).
- Eisenstein et al. (2002) J. P. Eisenstein, K. B. Cooper, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 88, 076801 (2002).
- 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).
- Csáthy et al. (2005) G. A. Csáthy, J. S. Xia, C. L. Vicente, E. D. Adams, N. S. Sullivan, H. L. Stormer, D. C. Tsui, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 94, 146801 (2005).
- 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).
- (51) Refs. Peterson et al., 2008a, b also showed the emergence of the so-called three-fold topological degeneracy which is a signature of the MR Pfaffian state on the torus geometry with the inclusion of finite thickness.
- Simon et al. (2007) S. H. Simon, E. H. Rezayi, N. R. Cooper, and I. Berdnikov, Phys. Rev. B 75, 075317 (2007).
- Regnault et al. (2009) N. Regnault, B. A. Bernevig, and F. D. M. Haldane, Phys. Rev. Lett. 103, 016801 (2009).
- Plenio and Virmani (2007) M. B. Plenio and S. Virmani, Quant. Inf. Comp. 7, 1 (2007).
- Eisert et al. (2010) J. Eisert, M. Cramer, and M. B. Plenio, Rev. Mod. Phys. 82, 277 (2010).
- Kitaev and Preskill (2006) A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 (2006).
- Levin and Wen (2006) M. Levin and X.-G. Wen, Phys. Rev. Lett. 96, 110405 (2006).
- Zozulya et al. (2007) O. S. Zozulya, M. Haque, K. Schoutens, and E. H. Rezayi, Phys. Rev. B 76, 125310 (2007).
- Iblisdir et al. (2007) S. Iblisdir, J. I. Latorre, and R. Orús, Phys. Rev. Lett. 98, 060402 (2007).
- Läuchli et al. (2010a) A. M. Läuchli, E. J. Bergholtz, and M. Haque, New Journal of Physics 12, 075004 (2010a).
- Li and Haldane (2008) H. Li and F. D. M. Haldane, Phys. Rev. Lett. 101, 010504 (2008).
- Hermanns et al. (2010) M. Hermanns, A. Chandran, N. Regnault, and B. A. Bernevig, ArXiv e-prints (2010), eprint 1009.4199.
- Chandran et al. (2011) A. Chandran, M. Hermanns, N. Regnault, and B. A. Bernevig, ArXiv e-prints (2011), eprint 1102.2218.
- Zhao et al. (2011) J. Zhao, D. N. Sheng, and F. D. M. Haldane, Phys. Rev. B 83, 195135 (2011).
- Sterdyniak et al. (2011) A. Sterdyniak, N. Regnault, and B. A. Bernevig, Phys. Rev. Lett. 106, 100405 (2011).
- Papić et al. (2011) Z. Papić, B. A. Bernevig, and N. Regnault, Phys. Rev. Lett. 106, 056801 (2011).
- Qi et al. (2011) X. Qi, H. Katsura, and A. W. W. Ludwig, ArXiv e-prints (2011), eprint 1103.5437.
- Thomale et al. (2010) R. Thomale, A. Sterdyniak, N. Regnault, and B. A. Bernevig, Phys. Rev. Lett. 104, 180502 (2010).
- Papic (2011) Z. Papic (2011), private communication.
- Bishara and Nayak (2009) W. Bishara and C. Nayak, Phys. Rev. B 80, 121302 (2009).
- Rezayi and Simon (2011) E. H. Rezayi and S. H. Simon, Phys. Rev. Lett. 106, 116801 (2011).
- Wójs et al. (2010) A. Wójs, C. Tőke, and J. K. Jain, Phys. Rev. Lett. 105, 096802 (2010).
- Peterson et al. (2008c) M. R. Peterson, K. Park, and S. Das Sarma, Phys. Rev. Lett. 101, 156803 (2008c).
- Haldane (1983) F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983).
- Greiter et al. (1991) M. Greiter, X.-G. Wen, and F. Wilczek, Phys. Rev. Lett. 66, 3205 (1991).
- Ando et al. (1982) T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 54, 437 (1982).
- Stern and Das Sarma (1984) F. Stern and S. Das Sarma, Phys. Rev. B 30, 840 (1984).
- Sarma and Mason (1985) S. D. Sarma and B. A. Mason, Annals of Physics 163, 78 (1985), ISSN 0003-4916.
- Fano et al. (1986) G. Fano, F. Ortolani, and E. Colombo, Phys. Rev. B 34, 2670 (1986).
- Read and Rezayi (1999) N. Read and E. Rezayi, Phys. Rev. B 59, 8084 (1999).
- Bonderson and Slingerland (2008) P. Bonderson and J. K. Slingerland, Phys. Rev. B 78, 125323 (2008).
- Töke et al. (2005) C. Töke, M. R. Peterson, G. S. Jeon, and J. K. Jain, Phys. Rev. B 72, 125315 (2005).
- Friedman and Levine (2008) B. A. Friedman and G. C. Levine, Phys. Rev. B 78, 035320 (2008).
- Friedman and Levine (2010) B. A. Friedman and G. C. Levine, Int. J. Mod. Phys. B 24, 4707 (2010).
- Rezayi and Haldane (2000) E. H. Rezayi and F. D. M. Haldane, Phys. Rev. Lett. 84, 4685 (2000).
- Läuchli et al. (2010b) A. M. Läuchli, E. J. Bergholtz, J. Suorsa, and M. Haque, Phys. Rev. Lett. 104, 156404 (2010b).
- Bernevig and Haldane (2008) B. A. Bernevig and F. D. M. Haldane, Phys. Rev. Lett. 100, 246802 (2008).
- Turner et al. (2009) A. M. Turner, Y. Zhang, and A. Vishwanath, ArXiv e-prints (2009), eprint 0909.3119.
- Feiguin et al. (2008) A. E. Feiguin, E. Rezayi, C. Nayak, and S. Das Sarma, Phys. Rev. Lett. 100, 166803 (2008).