Bipartite entanglement entropy in fractional quantum Hall states
Abstract
We present a detailed analysis of bipartite entanglement entropies in fractional quantum Hall (FQH) states, considering both abelian (Laughlin) and nonabelian (MooreRead) states. We derive upper bounds for the entanglement between two subsets of the particles making up the state. We also consider the entanglement between spatial regions supporting a FQH state. Using the latter, we show how the socalled topological entanglement entropy of a FQH state can be extracted from wavefunctions for a limited number of particles.
pacs:
03.67.Mn, 71.10.w, 73.43.fI Introduction
The fractional quantum Hall (FQH) states have long fascinated the condensedmatter community due to their remarkable transport properties and the exotic nature of their quasiparticle excitations. It is in the context of FQH states that the notion of topological order in gapped twodimensional states first arose. WenNiu_PRB () Recently there has been enhanced interest in FQH states with nonabelian statistics, MooreRead_NP91 (); ReadRezayi_PRB99 (); Ardonne_Schoutens99 () due to the possibility of implementing quantum computation schemes topologically protected from decoherence. topologicalquantumcomputing () The unusual features of FQH states have been notoriously difficult to characterize using traditional condensedmatter concepts such as local order parameters and point correlation functions.
In a separate development, recent years have seen growing understanding that entanglement measures borrowed from the discipline of quantum information can be useful in probing global features of quantum manyparticle states. VidalLatorreRicoKitaev_PRL03 (); Entnginspinchains (); Cardy_JStatMech04 (); Eisert_2002 () It is thus natural to ask what features of FQH states can be characterized by entanglement measures. In a recent short report,ourprl07 () three of the present authors have shown that one such entanglement measure, the bipartite entanglement entropy, indeed elucidates the subtle correlations and topological order in the simplest FQH states, those in the socalled Laughlin sequence. The bipartite entanglement entropy is defined by partitioning the system under question into two blocks and , and using the reduced density matrix of one part (e.g., obtained by tracing over degrees of freedom) to calculate the von Neumann entropy .
In Ref. ourprl07, , numerical calculations of the entanglement entropy between two spatial regions allowed us to extract from the Laughlin wavefunctions the socalled topological entanglement entropy, a concept introduced in Ref. PreskillKitaev_PRL06, ; LevinWen_PRL06, . [For brevity we write topological entropy where no confusion can arise.] In addition, we provided results on the entanglement entropy between subsets of the particles making up the state. We showed that such particle entanglements are bounded by expressions that manifest the exclusion statistics in the Laughlin states.
In this article, we present a systematic discussion of bipartite entanglement entropies for FQH states, elaborating on our results in Ref. ourprl07, . In addition to the abelian Laughlin (L) states, we consider a series of nonabelian FQH states: the MooreRead (MR) (or pfaffian) states. MooreRead_NP91 (); ReadRezayi96 (); Read_2000overview () In planar geometry, the respective wavefunctions are given by
with denoting the antisymmetric Pfaffian symbol. We shall here consider these same states in spherical geometry, so as to eliminate boundary effects.
For both series of states, we derive upper bounds for particle entanglement entropies. A marked difference between the Laughlin state and the MooreRead states is that in the latter the leading correlations have a 3body nature, whereas those in the Laughlin states are 2body effects. This difference is nicely manifested in the leading terms of a expansion of the upper bounds , which are given by (see section III below for details)
Another marked difference between the abelian and nonabelian states is in the value for the topological entropy . Comparing the Laughlin and MooreRead states at the same filling fraction we have
the difference being due to the nonabelian nature of the MooreRead states (see section IV for details). In this article we extract values for directly from wavefunctions for a limited number of particles in spherical geometry (up to for the Laughlin state and up to for the MooreRead state), finding values that are consistent with the expected result. These results illustrate how the entanglement entropy can be used in diagnosing the topological order for a FQH state that is only known in the form of wavefunctions for a limited number of particles, as is often the case in numerical studies.
In Section II we give a general discussion of possible bipartite entanglement measures in itinerant many body systems, carefully distinguishing between particle and spatial partitioning schemes. In Section III we present analytical and numerical results for particle entanglement in the Laughlin and MooreRead FQH states. In this, the eigenvalue distribution of the reduced density matrix plays a central role. In subsection III.6 we relate the eigenvalue distribution for to the twoparticle correlation function . In Section IV we discuss spatial partitioning, paying particular attention to the numerical procedure followed in extracting the topological entropy .
While most results in this article are for fermionic FQH states (meaning odd in the Laughlin sequence and even for the MooreRead states), we briefly comment on bosonic states in subsection III.3.
Ii Partition choices on the sphere
The entanglement entropy, being a bipartite measure of entanglement, depends on the particular partitions being considered. Obviously, a manyparticle system can be partitioned in many ways. Rather than asking which partition is the “correct” one, we find it more useful to ask what information one can extract from various kinds of partitioning. Accordingly, we have partitioned both the spatial degrees of freedom, and the particles themselves. We find that both schemes are useful, for revealing distinct features of the manyparticle state.
The kinds of partitioning one is able to study depend on the available degrees of freedom. Our calculations are all performed for FQH states in a spherical geometry. Haldane_FQHE_PRL83 (); ArovasAuerbachHaldane_PRL88 () In this representation the fermions are placed on a sphere containing a magnetic monopole. The magnetic orbitals of the relevant Landau level are then represented as angular momentum orbitals; the total angular momentum is half the number of flux quanta, . The orbitals are labeled either to or to . For particles at fractional filling , one finds the interesting FQH states for , where is a finitesize shift. The Laughlin states appear at while for the MooreRead states . The “filling” acquires the usual meaning only in the thermodynamic limit. The orbitals are each localized around a “circle of latitude” on the sphere, with the orbital localized near one “pole.”
Since the FQH wavefunctions on a sphere are obtained in terms of orbital occupancies, one can either partition orbitals or partition particles. Because of the spatial arrangement of the orbitals, partitioning orbitals is in fact equivalent to partitioning spatial regions. The difference between spatial and particle partitioning has not been stressed in the literature because the most common systems studied (in the context of entanglement entropies in manyparticle states) are spin models, for which there is no such distinction. In the cases of itinerant particles where there is a difference, the common default scheme has been spatial partitioning. In particular, conformal field theory results on entanglement scaling HolzheyWilczek_NPB95 (); Cardy_JStatMech04 () and the distinction between gapless and gapped states observed in entanglement scaling VidalLatorreRicoKitaev_PRL03 () actually pertain to the blocking of space rather than the particles or spins themselves. The definition of the topological entropy for twodimensional topologically ordered states is also based on the entanglement entropy between spatial blocks. PreskillKitaev_PRL06 (); LevinWen_PRL06 ()
In previous work, ourprl07 () three of the present authors studied the entanglement entropy between subsets of particles making up a Laughlin FQH state. We presented upper bounds and gave an interpretation in terms of exclusion statistics. In this paper we extend these results to the MooreRead states, where the exclusion effects are more intricate. We refer to Refs. Tsinghua_FQHE_PRA02, ; EntngQHE_Lattore, ; Pcle_entanglement_Santachiara, for other studies of particle entanglement properties.
For orbital or spatial partitioning, we define block to be the first orbitals, extending spatially from one pole of the sphere out to some latitude. In the thermodynamic limit, this is equivalent to choosing a diskshaped block within an infinite planar system. In this limit, since each orbital is associated with a wavefunction of the form in usual complex coordinate language, a disk with orbitals has radius .
The spatial arrangement of the orbitals constrains us to either disc or ringshaped spatial regions as the partition. Since the orbital indexes do not give us access to the full twodimensional degrees of freedom, it is not possible to experiment with the various kinds of topologically nontrivial partitions suggested by Preskill and Kitaev, PreskillKitaev_PRL06 () Levin and Wen, LevinWen_PRL06 () and Furukawa and Misguich. Misguich_dimergamma_dec06 () However, as we have reported previously, ourprl07 () the spherical geometry is sufficient to probe the topological entropy of FQH wavefunctions.
Comparing particle and spatial entanglement entropies, we remark that the effect of correlations is opposite between the two cases. For particle partitioning, the maximal entropy is realized for uncorrelated fermions; correlations tend to lower from the fermion bound . Spatial entanglement, on the other hand, is entirely due to correlations. The point is illustrated by considering the Laughlin state, where the fermions are uncorrelated. Here the particle entanglement entropy equals , while the spatial entanglement entropy vanishes.
Iii Entanglement for particle partitioning
In this section we provide close upper bounds to the entropy of entanglement between particles of the state and the remaining particles. We also discuss the particle reduced density matrices that arise in this context.
iii.1 Multiplet structure and particle entropy bounds
For FQH states on a sphere, the particle reduced density matrices commute with the total angular momentum operators and of the selected particles. This implies that the eigenvalues of are organized in a multiplet structure of the corresponding algebra: an eigenvalue for total angular momentum will be fold degenerate.
For fermions, each having angular momentum , the 2particle states have total angular momenta , , , , for integer (halfinteger), giving a total number of states. A naive upper bound to the entanglement entropy is thus
(1) 
Inspecting the explicit structure of the fermionic Laughlin states with , one finds that the eigenvalues corresponding to 2particle states with , , , all vanish. The reason is that the correlations in the Laughlin states are such that particles cannot come too close together. For example, if a first fermion occupies the orbital, localized near the north pole, the Laughlin wavefunction has zero amplitude for finding a second fermion in orbitals , , , . The highest possible value of the angular momentum of the two fermions combined is thus . The remaining number of nonzero eigenvalues is , leading to an improved bound on the entropy
(2) 
with as before. For , the multiplet structures are more complicated and we need to resort to a different method for finding a nontrivial upper bound to the particle entropy. In the next subsection we give a general derivation for both the Laughlin and the MooreRead series of fermionic FQH states.
iii.2 Upper bounds for fermionic states
For fermionic particles, particles in the block, and the total number of orbitals given by , fermionic statistics lead to an obvious upper limit to the entropy
(3) 
In the FQH states the correlations are such that the particles avoid each other and the entropy is further reduced. To obtain a handle on this, one may reason as follows. The model FQH states in the Laughlin and MooreRead series can be characterized as zeroenergy eigenstates of a Hamiltonian penalizing pairs and/or triplets of particles coming to the same position. After tracing out the coordinates for the set, the dependence on those in the set is such that one still has a zeroenergy eigenstate. However, the number of orbitals available to the particles is larger than what is needed to make the model FQH state in the sector, and one instead has a certain number of quasiholes on top of the set model state. The total ground state degeneracy for this situation has been studied in the literature: see Ref. ReadRezayi96, for the Laughlin and MooreRead states and Ref. NAcounting, for the ReadRezayi and ArdonneSchoutens series of nonabelian FQH states.
For the Laughlin states the details are as follows. The particle Laughlin state is realized on a total of Landau orbitals, corresponding to flux quanta. The Laughlin state for particles would need flux quanta; we thus have an excess flux of . With the Laughlin gauge argument this corresponds to the presence of quasiholes over the groundstate. According to Ref. ReadRezayi96, each of the quasiholes has a number of effective orbitals to choose from, with bosonic counting rules (meaning that two or more quasiholes can be in the same effective orbital). This gives a number of quasihole states equal to
leading to the following upper bound to the entropy
(4) 
We remark that this expression has a clear interpretation in terms of exclusion statistics: the counting factor in Eq. (4) gives the number of ways particles can be placed in orbitals, in such a way that a particle placed in a given orbital excludes particles from orbitals with .
In a expansion we find (assuming )
(5)  
The particle entropy reaches a maximum for . For this case our Eq. (4) gives, in the limit of large ,
(6) 
This bound is sharper than a bound recently presented in Ref. EntngQHE_Lattore, , which gives a larger coefficient for the linearin behavior.
For the fermionic MooreRead states at , with , we can reason in a similar way, with now . As for the Laughlin states we have an excess flux of but now the number of quasiholes is twice this number due to the fact that the fundamental quasiholes correspond to half a flux quantum. Thus, . We now take from Ref. ReadRezayi96, the following result for the total quasihole degeneracy
(7) 
This gives us an upper bound as before.
Putting , one easily checks that coincides with for . In a expansion, the leading deviation from is a 3body term at order ,
(8) 
This result nicely illustrates the fact that the leading correlations in the MooreRead state have a 3body character: the wavefunction vanishes if at least three particles come to the same position.
For the leading correlations do have a 2body character, as for the Laughlin states,
(9) 
Inspecting the particle entanglement at and for large, our bound implies that for the MooreRead state
(10) 
This bound is reduced from the Fermi bound , but it is larger than the bound for the (bosonic) Laughlin state, which has asymptotic form . This indicates that, at equal filling , the particles in a MooreRead state are more entangled than those in a Laughlin state.
The quasihole counting rules for the order clustered spinpolarized (ReadRezayi) and spinsinglet (ArdonneSchoutens) states are all known in the literature. NAcounting () They can be used to generalize the upper bounds on particle entanglement entropy given in this subsection to these more intricate nonabelian FQH states.
iii.3 Bosonic quantum Hall states
We briefly comment on the case of bosonic FQH states. The realization that a rapidly rotating Bose gas may eventually enter a regime of bosonic quantum Hall states motivates the theoretical study of the effects of bosonic statistics.
We consider bosonic Laughlin states at filling fraction with . The naive upper bound to the the entropy associated to placing bosons in orbitals is
(11) 
The expression for remains unchanged, giving the following leading correction in a expansion
(12) 
For a bosonic MooreRead state, with filling fraction with , the leading correction becomes
(13) 
In the case the leading correlations have 3body character, leading to the vanishing of the leading correction.
iii.4 Numerical results
In deriving the upper bound we relied on the fact that a certain number of eigenvalues of the reduced density matrix vanish. The bounds would be exact if all nonzero eigenvalues were equal, but since they are not the bounds overestimate the actual values for the entropies.
Fig. 1 plots the eigenvalues for the particle reduced density matrix for particles on a sphere in the Laughlin state, for which the single particle angular momentum is . The horizontal axis represents the degeneracy of the eigenvalues, in descending order. The eigenvalue at , with degeneracy 47, vanishes; the nonzero eigenvalues show some scatter around an asymptotic value. Due to this scatter the entropy is somewhat lower than the upper bound .
An important difference between the Laughlin and the MooreRead states is the absence of vanishing eigenvalues for the 2particle reduced density matrix. The eigenvalue distribution shown in Fig. 2 illustrates this point.
iii.5 Corrections to due to eigenvalue spread
It is interesting to consider in some detail the deviation between the bounds and the actual entropies computed numerically. As mentioned above, this deviation arises from the fact that the nonzero eigenvalues of the reduced density matrices are not all equal.
To estimate the effect on of the spread in the nonzero eigenvalues, we do a rough modeling of the eigenvalue distribution (Fig. 1) of for the Laughlin state. For this case the number of nonzero eigenvalues is . If these nonzero eigenvalues were all equal (to ), the entanglement entropy would have the maximum value , which is the predicted upper bound (4). We now take into account the deviations from , guided by Fig. 1, with the following toy distribution: we take out of of the eigenvalues to be equal to , with , while the rest of the eigenvalues are at value , , such that the sum of eigenvalues is unity. Assuming leads to
(14) 
Guided by the eigenvalue distributions in Fig. 1, we assume that is of order ; for concreteness we put equal to the multiplicity of the largest eigenvalue, which is . Taking between 1.2 and 1.5 (as observed for the largest available Laughlin wavefunctions) gives a correction in the entropy with coefficient in the range 0.03 – 0.14.
Fitting the difference to a form gives a coefficient . The vanishing eigenvalues account for , see Eq. (5), and we see that the remaining difference is consistent with the correction due to the spread in the nonzero eigenvalues.
We made similar estimates for a the MooreRead state with up to particles, where the eigenvalues are all nonzero and agrees with . In this case the deviation between data and bound show a dependence with a coefficient of about 0.14.
These considerations are of some general interest, as they make the point that a expansion of particle entanglement entropies are indicative of correlations in a manybody state. In the concrete case studied here, the sizeable value of the leading correction in the Laughlin state indicates strong 2body correlations, while the small value for the MooreRead state indicates the absence of such correlations.
iii.6 Reduced density matrices and correlation functions
Since the particle reduced density matrices are obtained by integrating out all but of the particles, one expects these matrices to be related to the particle correlation functions. In this subsection we study this relation for the case . In particular, the eigenvalue distributions of in Figs. 1, 2, although discrete, are reminiscent of the wellknown twoparticle correlation functions for Laughlin and MooreRead states. We will show that the eigenvalue distributions are in fact very closely related to the correlation functions – the eigenvalue distribution function is a kind of discretized version of .
The twoparticle correlation function is conventionally defined as
(15) 
where is the number of particles and is a density, which is chosen such that .
We express the 2particle reduced density matrix on a sphere in a basis of polar spherical coordinates. Because of the rotational symmetry it should be a function of the angular distance between the two particles. In Appendix A we show that can be written in the form
(16) 
In this expression, is the eigenvalue with multiplicity , corresponding to the total angular momentum of the two particles equal to . The functions are explicitly given in Eq. (32).
Since the distance between two particles is simply equal to , with the radius of the sphere, the 2particle correlation function is directly proportional to . Through Eq. (16) it is expressed as a transfrom from space to space, with basis functions . In Fig. 5, we show some curves for ; they agree with known results.
As illustrated in the inset to Fig. 5, the basis functions have a peak structure, with the position of the peak depending on the total angular momentum . Large values of correspond to small angular distances and vice versa. This is easy to understand from the following classical picture. When two particles with angular momenta have total angular momentum , the corresponding vectors and should point a the same direction. The angular momentum vector points at the position of the particle on a sphere, therefore two particle should be close to each other. On the other hand, if total angular momentum is zero then the angular momentum vectors should point into opposite directions. This means that particles are placed at the opposite sides of the sphere.
The fact that the are localized functions, peaked at values monotonically decreasing with , indicates that curve is simply a continuous form of the versus descending curves of Figs. 1 and 2. The similarity between the discrete and the continuous curves is not accidental.
At small distances . For a Laughlin state the lowest value of is , thus . This behavior is a direct consequence of the vanishing of eigenvalues with the largest multiplicities. For the MooreRead state the lowest value of is because there are no vanishing eigenvalues. Therefore at small distances .
Of course, our observations on the 2particle correlations agree with known results; our main point has been to stress the intimate relation with the eigenvalue distribution of the 2particle reduced density matrices.
Iv Spatial entanglement and topological entropy
We now turn to dividing the Landau level orbitals into two blocks and calculating the entropy of entanglement between them. In a previous publication, ourprl07 () we used this scheme to extract the topological entropy of the Laughlin state. Here, we mainly focus on the () MooreRead state. After reviewing the topological entropy and the total quantum dimension , especially in the context of the MooreRead state (IV.1), we detail some issues with taking the thermodynamic limit necessary for extracting from numerical data (IV.2), and then present our numerical results (IV.3). We also present observations on the spectral structure of the reduced density matrices (IV.4).
iv.1 Topological entropy for the MooreRead state
For spatial partitioning of manyparticle states, the general rule (“area law”) is that the entanglement entropy scales as the size of the boundary between the and blocks. Srednicki_PRL93 () Subtle information about the nature of the manyparticle state can be provided by the presence or absence of logarithmic corrections, values of coefficients, or subleading terms in this basic relationship. For topologically ordered states in two dimensions, the following theorem has been presented recently PreskillKitaev_PRL06 (); LevinWen_PRL06 () concerning the scaling of entanglement entropy between spatial partitions. If is the length of the boundary between the two blocks, the entanglement entropy scales as . As usual the scaling law applies to situations where is large and the total system is infinite. The subleading term is called the topological entanglement entropy. The striking result of Refs. PreskillKitaev_PRL06, ; LevinWen_PRL06, has been that it can be expressed as the logarithm of a quantity known as the total quantum dimension of the topological field theory describing the topological order of the state. The total quantum dimension is given by
(17) 
where the ’s are the quantum dimensions of the individual sectors making up the topological field theory. These quantum dimensions are set by fusion rules of the fundamental anyons in the field theory.
The topological field theory for a Laughlin state has a fundamental anyon (of fractional charge ), which generates abelian sectors. The quantum dimension is unity in all sectors so that for the Laughlin state. For this gives .
For states with nonabelian quasiparticles, the situation is more interesting
because some anyon sectors contribute . Details for some examples have
been provided in Refs. PreskillKitaev_PRL06, ; FendleyFisherNayak_JStatPhys07, . In particular, for the MooreRead
state, there are six sectors (two each of quasiparticles denoted by ,
, ) which contribute , ,
, leading to and . The
nonabelian nature shows up in the fact that is larger than
, six being the degeneracy of the MooreRead state on the
torus.
A compact general expression for the total quantum dimension for a ReadRezayi state with order clustering and at filling fraction is (see also Ref. FendleyFisherNayak_JStatPhys07, )
(18) 
It includes the Laughlin states (, ) and the MooreRead states (, ) as special cases. For a general spinsinglet nonabelian FQH state with order clustering and filling fraction , the result is
(19) 
giving for the paired spinsinglet state (, ) at .
In general, the dependence of these expressions for total quantum dimensions is linked to the ground state degeneracy in torus geometry. Denoting the latter by we have the relation
(20) 
The conformal field theories underlying the states at are of WessZuminoWitten type ( for the RR states and for the AS series) and the quantities can be expressed in the modular matrix for these WZW models.
iv.2 from sphere calculations: extrapolations
First of all, we note that, since our degrees of freedom are ordered essentially onedimensionally, we cannot use one of the twodimensional schemes proposed previously PreskillKitaev_PRL06 (); LevinWen_PRL06 (); Misguich_dimergamma_dec06 () in which an appropriate addition/subtraction of the entanglement entropies of several regions cancels the boundary parts of the entropy () leaving the subleading term . With the orbital degrees of freedom on a sphere, we can choose only regions corresponding to disks, concentric rings, and combinations thereof. Any combination of entropies of disk and ringlike regions that cancels out the boundary terms also unfortunately cancels out the term.
We are thus led to using directly the scaling law, . Our choice of block as the first orbitals, extending spatially from one pole out to some latitude, corresponds to a diskshaped block only in the thermodynamic limit. The block area is proportional to the square of while its boundary is proportional to ; these are equivalent only in the same limit. One way to numerically access the thermodynamic limit is to take the entanglement entropy of orbitals with the rest, for accessible wavefunctions of various sizes , and then take the limit. The versus points thus obtained should then follow a linear curve at large , whose vertical intercept gives the topological entropy. Results following this procedure were provided for the Laughlin state in our earlier paper; ourprl07 () here we will focus on the MooreRead state. The extrapolation of values to the thermodynamic limit is a tricky issue. We therefore discuss the extrapolation in some detail here, providing some general results.
We are interested in the function , where . We have access to at several integer values of , and would like to estimate . For each dataset that we have access to (each ; both Laughlin and MooreRead), we note the following: the versus values form a monotonic curve and this curve gets flatter (slope magnitude decreases) with decreasing . Two examples can be seen in the inset to Fig. 6. In other words, the first and second derivatives of the function have the same sign and neither derivative changes sign.
Motivated by the above observations, we provide the following result. Assuming only that the signs of the first two derivatives of the function are the same and that the signs remain unchanged until , we have:

The value corresponding to the smallest value of the available is a strict lower (upper) bound for if the is negative (positive).

The intercept found by connecting the corresponding to the smallest two values (, ), namely
is a strict upper (lower) bound if is negative (positive) and is positive (negative).
The limits and thus obtained give us conservative bounds for the required entanglement entropies in the thermodynamic limit, . To obtain a sharper extrapolation, one can use various polynomial extrapolations and take the average, as done in our earlier work. ourprl07 () Here, we improve the extrapolation by using the extrapolation algorithm of Bulirsh and Stoer (BST algorithm), based on rational polynomial fraction approximations. xtrap_BulirshStoer64 (); xtrap_HenkelSchuetz88 ()
2.0  
3.46601  5.21183  4.69023  4.77647 
3.80408  4.92355  4.74756  
4.01821  4.84476  
4.16743  
3.0002  
3.46601  4.56058  4.51976  4.55722 
3.80408  4.53643  4.59865  
4.01821  4.56761  
4.16743  
2.5926  
3.46601  4.74638  4.58027  4.65108 
3.80408  4.65109  4.65108  
4.01821  4.65108  
4.16743 
The BST procedure involves successive transformations of the original dataset, leading to transformed datasets successively smaller by one element. If the free parameter () of the algorithm is chosen properly, the successive sets will be more and more convergent and will eventually converge to the value. The parameter is chosen to optimize convergence; how precisely this is done is a nontrivial issue, especially for large datasets.
In the present case where the initial dataset consists of bipartite entanglement entropies, it is particularly desirable to be able to use a small number of initial values, corresponding to the largest available . This is because the symmetry makes the data meaningless for extrapolation when gets close to . In other words, the values used should not be too close to the peak at the midpoint of the respective versus curve (e.g., Figure 1 of Ref. ourprl07, ). This restricts us to the largest few values.
In Table 1, we demonstrate some possible choices of starting from four values. The top example uses , which is a common choice used for extrapolationbased numerical integration algorithms. In the second example, is tuned to give the “most converged” set after the first iteration, by choosing to minimize the standard deviation of the second column. We find that this procedure does not always lead to great stability for the subsequent iteration results.
In the last example, we have tuned to give the same last value for the sequences obtained after the first and the second iterations. This also guarantees the third (final) iteration to give the same value (Table 1), which indicates that the convergence is very good. This procedure has the added advantage that actually only the last three points of the initial dataset are used to determine the values tuned to be equal. We thus get a nicely converged estimate based only on three initial points, which is highly desirable as explained above. This is therefore our method of choice for tuning the free parameter in applying the BST algorithm for our problem.
As a side note, we remark that if the largest available entanglement entropies are not known exactly (e.g., if they are calculated from approximate wavefunctions, or using the desnity matrix renormalization group (DMRG) technique), it might be more prudent to use a prescription for choosing that actually uses more than the minimal (three) numbers that we have used here. In the cases reported here and in our earlier work, ourprl07 () we have entanglement entropies calculated from numerically exact wavefunctions, so we have no precision issues to worry about when choosing the number of values to use for the extrapolation.
Once the extrapolation value is (uniquely) determined by the algorithm outlined above, we need an estimate of the uncertainty. We can use the bounds (,) derived earlier to obtain a conservative error estimate:
where is the extrapolated value obtained by the BST algorithm.
iv.3 Numerical results
MooreRead state. Fig. 6 shows results of numerical calculations for the MooreRead state. We used exact wavefunctions up to particles. These wavefunctions were obtained by diagonalizing in an Hilbert space spanned by the “squeezed states” Hal2 (). After numerically obtaining the entanglement entropies from these wavefunctions, we obtain estimates and uncertainties for the extrapolations by the procedure outlined in the previous subsection. The resulting data are plotted in Fig. 6.
The linear versus behavior is expected only for large ; however our large points have the greatest uncertainty. For estimating the topological entropy, we therefore make linear fits after discarding 0 to 5 of the smallest points and/or 0 to 2 of the largest points. This results in estimates of (magnitude of the vertical intercept) scattered between 0.85 and 1.35. The error propagated into our estimate from our extrapolation uncertainties is , larger than that obtained from this scatter. With all this we arrive at the result , quite consistent with the expected value of .
Laughlin state. We used the welldefined procedure of the previous subsection to revisit our previous estimate of the topological entropy for the Laughlin state.ourprl07 () To get the extrapolated , we now use the BST estimates rather than doing several polynomial fits. Dropping 0 to 4 of the smallest points and 0 to 2 of the lrgest points leads to , consistent with the previously reported estimate () and with the expected value . The error estimate reported in Ref. ourprl07, only took into account this variation, due to dropping various number of points. There is also some error propagated from the extrapolation uncertainty. Using the conservative uncertainty estimate proposed in the previous subsection gives us a more conservative and more rigorous error estimate, .
iv.4 Eigenvalue distribution for reduced density matrix
In Fig. 7, we show the largest eigenvalues of reduced density matrices obtained by orbital or spatial partitioning. The eigenvalues are ordered according to decreasing magnitude and plotted on a log scale; the resulting curves are roughly linear, suggesting a roughly exponential decay of the eigenvalue distribution function.
It is interesting to note the complete dissimilarity of this eigenvalue spectrum compared to the particle partitioning case discussed earlier, e.g., Fig. 1. It would also be interesting to put our observations in the context of the spectra of reduced density matrices of manybody systems in general. Reduced density matrices for spatially connected blocks have been studied previously in the context of the convergence of the DMRG algorithm; an overview is available in section IIIB of Schollwöck’s DMRG review. Schollwoeck_DMRGreview05 () From our numeric data, It is difficult to say whether or not the decay of the eigenvalue distributions is slower than exponential.
V Concluding remarks
We have presented a detailed study of the entanglement entropy in abelian and nonabelian quantum Hall states, taking a paradigmatic example of each, the Laughlin state and the MooreRead state.
For entanglement between subsets of particles, we have demonstrated the effects of particleparticle correlations in the deviation of the entanglement entropies from an upper bound set by fermionic statistics only. We presented a close upper bounds for particle entanglement entropies for both bosonic and fermionic states. For a particular case , i.e. entanglement of the one half of the system with the other, the presented upper bound is notably lower than a similar estimate in Ref. EntngQHE_Lattore, . We showed that the distribution of the eigenvalue spectrum of the twoparticle density matrix can be directly related to a twoparticle correlation function . This suggests that the eigenvalue distribution of particle reduced density matrices can be related to particle correlation functions.
For entanglement between spatial regions, indexed by orbitals on a spherical geometry, we have provided a thorough discussion of our efforts to optimize the procedure for extracting the socalled “topological entropy” from finitesize data. Using this procedure, we calculated the topological entropy for the MooreRead state, and found our numerical result to be consistent with the expected value .
Our results on Laughlin and MooreRead states indicate that the computation of topological entanglement entropies can be used for diagnosing topological order in FQH states that are only available for a limited number of particles. Clearly, the practical use of this method hinges on the accuracy by which a value for can be extracted. In this context, it could be particularly interesting to use entropies that are generated in DMRG studies of FQH states.
It is interesting to note the differences between the particle and spatial partitioning cases. In the first case, correlation effects decrease the particle entanglement entropy . In the latter case, correlation effects increase the spatial entanglement entropy; for example, a Wigner crystal would have very high entanglement entropy between two spatial partitions. These differences are particularly dramatically manifested in the difference between the spectral structure of the reduced density matrices for the two cases. In the particle entanglement case, the eigenvalues tend to be of the same order, with relatively small fluctuations around some average value. In the spatial case, the eigenvalue distribution function has an exponentiallike form, and only a small fraction of the eigenvalues contribute to the entanglement entropy.
Acknowledgements.
We thank U. Schollwöck for drawing our attention to the Bulirsh & Stoer (BST) algorithm for extrapolations and we acknowledge discussions with E. Ardonne (on eq. (19)), P. Calabrese. We thank A. Feiguin for discussions and providing us with DMRG data for largesize Laughlin wave functions. The research of OZ and KS is financially supported by the Stichting voor Fundamenteel Onderzoek der Materie (FOM) of the Netherlands. EHR is supported by US DOE under contract DEFG0302ER45981. KS and EHR acknowledge the Institute for Pure and Applied Mathematics (IPAM) at Los Angeles, where part of this work was done.Appendix A 2particle reduced density matrix and correlation functions
In this appendix we derive Eq. (16). The twoparticle reduced density matrix may be written as
(21) 
where is a state of two particles with total angular momentum and projection . The rotational symmetry of the sphere ( must be invariant under rotations) indicates that the eigenvalues do not depend on the projection , i.e, are grouped in multiplets; .
To get the correlation function, we need to compute the density matrix in a basis of polar spherical coordinates,
(22) 
where and are coordinates of particles on a sphere. Rotational symmetry requires that depends only on the angular distance between particles:
(23) 
Now the unnormalized wavefunction of two particles with total angular momentum can be written as Haldane_FQHE_PRL83 ()
(24) 
where is the angular momentum of each particle, , , and , fix the center of mass of the two particles. The same wavefunction can also be written as
(25) 
where is a rotation operator. Let us write , which is a familiar coherent state. Using the identity
(26) 
one obtains
(27) 
The normalization of is equal to Fano ()
(28) 
We need to compute
(29) 
Defining and , one finds
(30) 
One can see that is normalized as
(31) 
Introducing the function
(32) 
we finally write the 2body correlator as
(33) 
Appendix B Entanglement entropies; numerical results with orbital partitioning
1  0.636514  0.673012  0.682908  0.686962  0.689009  0.690186  0.690923  0.691416 

2  1.09861  1.15777  1.23519  1.24961  1.25796  1.2663  1.27095  1.27385 
3  1.09861  1.49971  1.65513  1.7002  1.72587  1.74591  1.75859  1.76716 
4  1.76712  2.0355  2.11077  2.15836  2.19257  2.21444  2.22959  
5  1.88152  2.31807  2.44895  2.52664  2.58123  2.61639  2.64089  
6  2.48295  2.70298  2.82429  2.90754  2.96113  2.99855  
7  2.54282  2.88348  3.06178  3.18053  3.257  3.31045  
8  2.98996  3.24156  3.40254  3.50653  3.57927  
9  3.02558  3.36728  3.57792  3.71402  3.80922  
10  3.44165  3.71089  3.88358  4.00428  
11  3.46601  3.80408  4.01821  4.16743  
12  3.85949  4.12054  4.30135  
13  3.87774  4.19219  4.408 
1  0.682908  0.673012  0.666278  0.661563  0.65811  0.655482  0.653418  0.651757 

2  1.00424  1.05492  1.07339  1.0822  1.08707  1.09006  1.09202  1.09337 
3  1.27703  1.3944  1.44547  1.46998  1.48491  1.49466  1.50147  1.50648 
4  1.27703  1.59387  1.70596  1.76262  1.79672  1.81933  1.8353  1.84717 
5  1.66184  1.87354  1.97796  2.03976  2.08068  2.10962  2.13113  
6  1.9523  2.1202  2.21872  2.28347  2.32919  2.36314  
7  1.9523  2.20091  2.3448  2.43838  2.50425  2.55307  
8  2.22768  2.42555  2.55286  2.64196  2.70783  
9  2.46451  2.63139  2.74713  2.83233  
10  2.46451  2.67761  2.82363  2.93049  
11  2.69293  2.87356  3.00467  
12  2.89817  3.05653 
In tables 2 and 3 we list orbitalpartitioning entanglement entropies calculated using the numerical wavefunctions. For each wavefunction (each colmun), the entanglement entropies are only listed up to their maximum value, typeset in bold, because the values after this are determined by the symmetry . The extrapolation procedure of Sec. IV.2 involves extrapolating each row of numbers to .
Note that, in each column (for a particular ), the omitted part after the midpoint (in bold) is a decreasing function of and thus does not give useful information about the thermodynamic limit of . In the extrapolation, it is therefore important to avoid values from these parts of the table. We therefore restrict ourselves to values which in tables 2 and 3 are to the right of (i.e., above) the diagonal line through the midpoint numbers typeset in bold.
References
 (1) X. G. Wen and Q. Niu, Phys. Rev. B 41, 9377 (1990); X. G. Wen, Phys. Rev. B 41, 12838 (1990); X. G. Wen, Phys. Rev. B 44, 2664 (1991).
 (2) G. Moore and N. Read, Nucl. Phys. B360, 362 (1991).
 (3) N. Read and E. Rezayi, Phys. Rev. B 59, 8084 (1999).
 (4) E. Ardonne and K. Schoutens, Phys. Rev. Lett. 82, 5096 (1999).
 (5) A. Yu. Kitaev, Annals Phys. 303, 2 (2003); M. Freedman, M. Larsen, and Z. Wang, Commun. Math. Phys. 227, 605 (2002); N. E. Bonesteel, L. Hormozi, G. Zikos, S.H. Simon, Phys. Rev. Lett. 95, 140503 (2005); S. Das Sarma, M. Freedman, C. Nayak, Phys. Rev. Lett. 94, 166802 (2005).
 (6) G. Vidal, J. I. Latorre, E. Rico and A. Kitaev, Phys. Rev. Lett. 90, 227902 (2003).
 (7) J. I. Latorre, E. Rico, and G. Vidal, Quantum Information and Computation, 4, 48 (2004); V.E. Korepin, Phys. Rev. Lett. 92, 096402 (2004); V. Popkov and M. Salerno, Phys. Rev. A 71, 012301 (2005).
 (8) P. Calabrese and J. Cardy, J. Stat. Mech. 0406, 002 (2004).
 (9) K. Audenaert, J. Eisert, M. B. Plenio, R. F. Werner, Phys. Rev. A 66, 042327 (2002).
 (10) M. Haque, O. Zozulya, and K. Schoutens, Phys. Rev. Lett. 98, 060401 (2007).
 (11) A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 110404 (2006).
 (12) M. Levin and X. G. Wen, Phys. Rev. Lett. 96, 110405 (2006).
 (13) N. Read and E. Rezayi, Phys. Rev. B 54, 16864 (1996).
 (14) N. Read, Physica 298B, 121 (2001).
 (15) C. Holzhey, F. Larsen, and F. Wilczek, Nucl. Phys. B 424, 443 (1995).
 (16) B. Zeng, H. Zhai, and Z. Xu, Phys. Rev. A 66, 042324 (2002).
 (17) S. Iblisdir, J. I. Latorre, and R. Orús, Phys. Rev. Lett. 98, 060402 (2007).
 (18) R. Santachiara, F. Stauffer, D. Cabra, condmat/0610402; P. Calabrese and M. Mintchev, condmat/0703117.
 (19) S. Furukawa and G. Misguich, condmat/0612227.
 (20) F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983).
 (21) D. P. Arovas, A. Auerbach, and F. D. M. Haldane, Phys. Rev. Lett. 60, 531 (1988).
 (22) V. Gurarie and E. Rezayi, Phys. Rev. B 61, 5473 (2000); E. Ardonne, N. Read, E. Rezayi and K.Schoutens, Nucl. Phys. B 607, 549 (2001); E. Ardonne, J. Phys. A 35, 447 (2002).
 (23) M. Srednicki, Phys. Rev. Lett. 71, 666 (1993); M. M. Wolf, F. Verstraete, M. B. Hastings, J. I. Cirac, arXiv: 0704.3906; M. B. Plenio, J. Eisert, J. Dreißig, M. Cramer, Phys. Rev. Lett. 94, 060503 (2005).
 (24) P. Fendley, M. P. A. Fisher, and C. Nayak, J. Stat. Phys. 126, 1111 (2007).
 (25) R. Bulirsh and J. Stoer, Numer. Math. 6, 413 (1964).
 (26) M. Henkel and G. Schütz, J. Phys. A 21, 2617 (1988).
 (27) F.D.M. Haldane, unpublished.
 (28) U. Schollwöck, Rev. Mod. Phys. 77, 259 (2005).
 (29) G. Fano, F. Ortolani, E. Colombo, Phys. Rev. B 34, 2670, (1986).