one of the most important questions in condensed matter physics relates to how we should characterize and classify the many different quantum phases of matter .
a large part of the story is the theory of symmetry breaking , which tells us that we should classify various phases based on the symmetries of the ground state wave function . yet with the discovery of the fractional quantum hall ( fqh ) states@xcite came also the understanding that there are many distinct and fascinating quantum phases of matter , called topological phases , whose characterization has nothing at all to do with symmetry .
how should we characterize and systematically classify the different possible topological phases that may occur in a fqh system ?
let us run through a few obvious possibilities .
we know that the fqh states contain topology - dependent and topologically stable degenerate ground states , which allow us to introduce the concept of topological order in fqh states.@xcite such topology - dependent degenerate ground states suggest that the low energy theories describing the fqh states are topological gauge field theories in 2 + 1 dimensions , called chern - simons theories.@xcite so one possibility is that we may try to classify the different fqh phases by classifying all of the different possible chern - simons theories .
but this is not a satisfactory approach for non - abelian fqh states because we do not have a good way of knowing which chern - simons theories can possibly correspond to a physical system made of electrons and which can not .
another possibility comes through the connection between fqh wavefunctions and conformal field theory ( cft ) .
it was discovered around 1990 that correlation functions in certain two - dimensional conformal field theories may serve as good model wavefunctions for fqh states.@xcite thus perhaps we may classify fqh states by classifying all of the different conformal field theories that may be used to construct a `` valid '' fqh wavefunction . while the connection between fqh wave functions and cft correlators has been extremely fruitful in both constructing new model wavefunctions and understanding their topological properties , the cft approach is incomplete in that there does not exist a complete classification of all possible conformal field theories that may be used to construct fqh wavefunctions .
this precludes the possibility of systematically classifying fqh states .
the cft approach also obscures the essential physics and the essential properties of the conformal field theory that allow one to obtain amenable fqh wavefunctions .    in an attempt to obtain a systematic classification of fqh states ,
it was shown recently that a wide class of single - component fqh states and their topological excitations can be classified by their pattern of zeros , which describe the way ideal fqh wavefunctions go to zero when various clusters of particles are brought together.@xcite this analysis led to the discovery of some new non - abelian fqh states whose corresponding cft has not yet been identified .
it also helped elucidate the role of cft in constructing fqh wavefunctions .
the cft encodes the way the wavefunction goes to zero as various clusters of electrons are brought together .
the order of these zeros must satisfy certain conditions and the solutions to these conditions correspond to particular cfts .
thus in classifying and characterizing fqh states , one can bypass the cft altogether and proceed directly to classifying the different allowed pattern of zeros and subsequently obtaining the topological properties of the quasiparticles from the pattern of zeros .
this construction can then even be thought of as a classification of the allowed cfts that can be used to construct fqh states .
furthermore , these considerations give way to a natural notion of which pattern of zeros solutions are simpler than other ones . in this sense , then , one can see that the moore - read pfaffian quantum hall state@xcite is the `` simplest '' non - abelian generalization of the laughlin state .    in this paper
, we generalize the pattern - of - zeros classification to multilayer fqh wave functions .
such a systematic classification leads us to the construction of non - abelian multilayer fqh states and also helps identify the simplest non - abelian generalizations of the halperin@xcite bilayer fqh states .
therefore , in the same way that the pfaffian fqh state is the simplest non - abelian single - layer state and so one of the first non - abelian states expected to be realized experimentally , our construction allows us to identify the simplest non - abelian bilayer states and therefore some of the ones that deserve further consideration in future experimental and numerical work .
we would like to point out that the `` 1d charge - density - wave '' characterization of single - component fqh states@xcite is closely related to the pattern - of - zeros approach .
our pattern - of - zeros approach for multilayer fqh states can also be viewed as a way to generalize the `` 1d charge - density - wave '' characterization to multicomponent cases.@xcite    this paper is organized as follows .
we begin by describing the ideal fqh wave functions that we can characterize by the pattern of zeros and their relation to symmetric holomorphic polynomials . in section [ pozcharacterization ]
, we describe the different ways to characterize the pattern of zeros . in section
[ conditions ] we find the conditions that the pattern of zeros must satisfy in order to describe valid fqh wave functions . in section [ idealhamiltonians ] ,
we sketch how one may begin to construct ideal hamiltonians whose ground states will be fqh wave functions with a given pattern of zeros . after a brief summary of the pattern of zeros data and conditions in section [ summary ]
, we explain in section [ cftsection ] the relation between the pattern - of - zeros approach and the cft approach to fqh wave functions . in section [ examples ] , we describe some example solutions of this systematic classification of multilayer fqh wave functions , which yields many non - abelian multilayer states . in section [ discussion ]
we discuss some of the simplest of these non - abelian bilayer states that may be relevant for experiments on two - component quantum hall systems and that warrant further numerical study .
the ground state wave function of a two - dimensional system of electrons in the lowest landau level can be written in the form @xmath4 where @xmath5 , @xmath6 are the coordinates of the @xmath7 electron , and @xmath8 is a holomorphic function of @xmath9 .
since the electrons obey fermi statistics , @xmath10 is anti - symmetric under interchange of any two coordinates @xmath9 and @xmath11 when all of the electrons are identical . in many physical situations ,
the electrons may be distinguished by various quantum numbers , such as a spin index ( when the zeeman energy is not too high ) , a layer index ( in a multilayer two - dimensional electron system ) , a valley index ( such as in graphene and sige heterostructures ) , _ etc . in such cases ,
the ground state wave function in the lowest landau level may instead be written in the form @xmath12 where @xmath13 is a flavor index and @xmath14 is the number of different flavors .
@xmath10 is then necessarily antisymmetric only under interchange of @xmath15 and @xmath16 for any @xmath17 and @xmath18 . given any antisymmetric polynomial @xmath19 , we can uniquely construct a symmetric polynomial : @xmath20 @xmath21 will also be a polynomial because @xmath22 must vanish when any two identical particles approach each other .
thus the above division by the factor @xmath23 will never produce any poles in the resulting function .
_    therefore to classify fqh phases of electrons , we can restrict our attention mainly to symmetric polynomials @xmath24 , where @xmath10 is invariant under the interchange of @xmath15 and @xmath16 for any @xmath17 and @xmath18 , but not necessarily invariant under the interchange of @xmath15 and @xmath25 if @xmath26 . in this paper
we will often refer to @xmath27 as a layer index . in the following
, @xmath10 will always refer to such a symmetric multilayer polynomial .
we will introduce data , such as @xmath28 , @xmath29 , and @xmath30 , to characterize and classify bosonic fqh states ( symmetric polynomials ) @xmath24 . from the above discussion
, we see that the same set of data also characterizes fermionic fqh states whose wave functions are given by @xmath31 .
progress in understanding the various possible topological phases in the fqh effect has occurred largely because of the discovery of various kinds of model wave functions and the ideal hamiltonians that yield these wavefunctions as their ground states . in the pattern of zeros
approach , we classify all of the possible ideal wave functions for fqh phases . in this section
we will explain what we mean by `` ideal . ''    for each topological phase in the fqh system
, we want to have a certain representative wave function that captures the topological properties of the phase .
the prototypical example of such an ideal wavefunction is the laughlin wave function at filling fraction @xmath32 : @xmath33 at every particle in @xmath34 , there are @xmath29 zeros , and there are no off - particle zeros .
the exact ground state of the microscopic hamiltonian with the coloumb interaction does not have this simple property .
nevertheless , the ideal laughlin wavefunction captures all of the essential topological properties of these phases .
the laughlin wavefunctions are also ground states of certain ideal hamiltonians that contain interaction potentials that are constructed only from @xmath35-functions and their derivatives .
for example , if the interaction potential is of the form @xmath36 then the following wavefunction is an exact ground state , with vanishing total potential energy : @xmath37 on the other hand , if the interaction potential between two electrons is of the form @xmath38 with @xmath39 and @xmath40 ,
then the zero - energy ground state will be @xmath41    more complicated ground states can be obtained through more complicated interactions .
for example , consider the following three - body interaction between electrons : @xmath42 where @xmath43 is the total symmetrization operator between @xmath44 , @xmath45 , and @xmath46 .
an exact zero - energy ground state of this interaction is the pfaffian wavefunction at @xmath47 : @xmath48 where @xmath49 is the total antisymmetrization operator between @xmath50 .
since the pfaffian and laughlin wavefunctions are exact zero - energy ground states of interaction potentials that are constructed only from @xmath35-functions and their derivatives , these wavefunctions can be characterized completely by their pattern of zeros , i.e. by the order of the zeros in the wavefunction as different numbers of particles are taken together .
these wave functions have certain ideal properties ; for example , the zeros that are bound to a particle lie precisely at the location of that particle , not slightly away from it . in this paper , we will classify ideal fqh wavefunctions , which are wave functions that are exact zero energy ground states of such ideal hamiltonians .
the existence of these ideal hamiltonians is crucial .
we can write down any arbitrary complex function of @xmath51 variables , but we can not know whether it corresponds to a valid topological phase of matter in the limit @xmath52 unless we also know that it is the ground state of some local , gapped hamiltonian .
if this is the case , then it is conceivable that there exists some physical situation in which the low energy effective interactions between electrons yields a phase of matter that is in the same universality class as the ideal wave functions that we constructed .
unfortunately , judging whether a many - particle interacting hamiltonian is gapped in the thermodynamic limit is in general an intractable problem .
so here we limit ourselves to classifying ideal wave functions , which at least we believe can be realized as ground states of local hamiltonians .
whether the corresponding hamiltonian can be gapped in the thermodynamic limit is a question that we must attempt to answer along the way .
the ideal quasihole wavefunctions for these states are also zero energy eigenstates of their corresponding ideal hamiltonians .
the true ground state is distinguished by being the unique translationally invariant state with lowest angular momentum .
thus we can also classify topologically distinct quasiholes by their pattern of zeros , i.e. by the order of the zeros in the quasihole wavefunction as different numbers of electrons are taken to the quasihole . while the quasihole wavefunctions are also zero energy states , the quasiparticles will be gapped and will have the same topological properties as the quasiholes .
therefore in the following , when we discuss the pattern of zeros of the quasiparticles , we are referring to the pattern of zeros of the quasihole wave functions .
not all fqh phases have such ideal wave functions .
the hierarchy states and the composite fermion states , for example , do not have ideal wave functions . while these phases do have their own model wave functions , they are not ideal in the sense that they can not be directly described by their pattern of zeros or written as a correlation function of conformal primary operators in a cft
. therefore the pattern - of - zeros construction does not directly classify these phases .
we will however discuss how they are related to the ideal multilayer fqh wave functions in section [ sec : n11 ] .
the spirit of the pattern - of - zeros approach is to consider bringing together @xmath53 particles of type @xmath27 , for @xmath13 , and asking how @xmath10 goes to zero under such a procedure .
the order of the zero will be denoted @xmath30 , where @xmath54 . in the following
we will more precisely define @xmath30 and discuss some different yet equivalent ways of characterizing the pattern of zeros .
this discussion is a straightforward generalization of the discussion in the single - layer case.@xcite      consider a set of @xmath53 coordinates of each type @xmath27 , and set @xmath55 .
define @xmath30 as the minimal power of @xmath56 in the polynomial @xmath10 .
this means that if we set @xmath57 and we take @xmath58 , then @xmath59 where @xmath60 is a polynomial in @xmath61 and the remaining coordinates @xmath62 and @xmath63 .
we refer to @xmath62 as the coordinate of an @xmath64-cluster .
we assume that @xmath30 is independent of the choice of @xmath62 , which must be the case for translationally invariant wavefunctions .
we also assume that @xmath30 is independent of the choice of @xmath65 and that different polynomials @xmath66 obtained from different choices of @xmath67 are linearly dependent .
this is the assumption of _ unique fusion .
_    we can immediately deduce some basic properties of @xmath30 .
since @xmath10 has no poles , it is clear that @xmath68 .
since @xmath10 must be single - valued under rotating @xmath69 in the complex plane by an angle @xmath70 , @xmath30 must be an integer .
let @xmath71 be the minimal power of @xmath72 ; that is , @xmath73 .
a translationally invariant @xmath10 will have @xmath74 , otherwise it will vanish everywhere .
thus , for a translationally invariant polynomial , @xmath30 is a nonnegative integer that characterizes the order of zero that results when the size of an @xmath64-cluster goes to zero .      in the previous section
, we introduced the derived polynomials @xmath76 . as a consequence of the unique fusion condition , these polynomials are actually independent of @xmath77 .
we may consider more general derived polynomials by bringing together other sets of coordinates in @xmath78 to obtain @xmath79 .
then we may consider bringing together an @xmath64-cluster and a @xmath80-cluster : @xmath81 thus , @xmath82 characterizes the order of the zeros in the derived polynomials as a cluster of @xmath64 electrons are brought close to a cluster of @xmath83 electrons .
the unique - fusion condition assumes that the derived polynomials obtained from different ways of fusion are always linearly dependent .
the fact that @xmath10 is a single - valued , symmetric polynomial implies @xmath84    we can deduce a relation between @xmath82 and @xmath30 as follows . the order of the zero obtained by creating an @xmath85-cluster is @xmath86 .
one way of creating such a cluster is by first creating an @xmath64-cluster , then creating a @xmath83-cluster , and finally bringing together the two clusters to create an @xmath87-cluster .
the order of zero in this case will be @xmath88 . thus @xmath82 can be obtained from @xmath89 through the formula @xmath90 since @xmath91 , where recall @xmath92 is the unit vector in the @xmath27 direction , we also have @xmath93 from this recursion relation and from the fact that @xmath91
, we may obtain @xmath30 from the sequence @xmath82 .
therefore we may equivalently label the pattern - of - zeros data using @xmath89 or @xmath94 .
the integer @xmath30 has the following meaning .
a polynomial with @xmath53 particles of @xmath95 kind ( @xmath53 particles in the @xmath95 layer ) has a total order of @xmath30 . in other words
the total angular momentum of the quantum hall droplet is @xmath30 if the droplet has @xmath96 particles in the @xmath95 layer .
if we remove one particle from the @xmath95 layer , then the total angular momentum of the quantum hall droplet will be reduced to @xmath97 .
thus we can interpret @xmath98 as the angular momentum of the highest occupied orbital in the @xmath95 layer for a quantum hall droplet with @xmath99 particles in the @xmath100 layer .
the @xmath14-dimensional sequence of vectors @xmath101 will be called the sequence of highest occupied orbitals ( hoo ) .
we see that @xmath102 makes sense only when @xmath103 .
we will set @xmath104 when it does not make sense . from
, we also see that there is one - to - one correspondence between the sequence @xmath30 and @xmath105 .
thus we can also use @xmath106 to characterize the pattern of zeros in the wave function .
a fqh wave function @xmath107 defined on a sphere forms a representation of @xmath108 . in such a case , @xmath109 represents the stereographic projection onto the plane of a point on the sphere .
a single particle in the lowest landau level can fill any of the @xmath110 orbitals , so the representation of @xmath108 formed in this case is the one with angular momentum @xmath111 .
the su(2 ) lie algebra is generated by @xmath112 in the multilayer case , then , the angular momentum of a particle of type @xmath27 will be @xmath113 where @xmath114 is the total number of flux quanta through the sphere seen by the particles in the @xmath95 layer . note that here we allow the numbers of flux quanta in different layers to be different .
the total angular momentum of an @xmath64 cluster in the @xmath109 direction will be the eigenvalue of the operator @xmath115 the operator @xmath116 counts the total power of a polynomial .
since the minimum total power of @xmath117 is @xmath30 , the minimum total angular momentum of an @xmath64-cluster is given by @xmath118 .
this means that the @xmath64-cluster carries an angular momentum of @xmath119 where @xmath120 and @xmath121 .
we will use this relation later to construct ideal hamiltonians and to place conditions on the pattern of zeros for when they can correspond to rotationally invariant wave functions on the sphere .
for the pattern of zeros to describe a valid fqh wavefunction , it must satisfy certain consistency conditions .
we already encountered several such conditions above .
for instance , we found that @xmath30 is a nonnegative integer , @xmath122 , and @xmath123 is even . in the following we find additional conditions that the pattern of zeros must satisfy .
one of the most important conditions on the wave function is simply the condition that the wave function have no poles .
this condition is remarkably restrictive on the allowed pattern - of - zeros sequences .
consider a derived polynomial @xmath124 and fix all coordinates but @xmath62 , thus viewing it as a complex function @xmath125 .
@xmath125 has zeros at isolated points , but no poles anywhere .
some of the zeros are located at @xmath126 , @xmath127 , .
those zeros are called on - particle zeros .
the rest of the zeros are called off - particle zeros .
if we imagine taking @xmath62 around @xmath126 without enclosing any off - particle zeros , then @xmath128 will pick up a phase @xmath129 . similarly ,
if we take @xmath62 around @xmath127 without enclosing any off - particle zeros , then @xmath128 will pick up a phase @xmath130 .
now consider taking @xmath131 .
under such a process , some nearby off - particle zeros will also be taken to @xmath127 . therefore , if we take @xmath62 around a contour that encloses both @xmath126 and @xmath127 in the limit that @xmath132 , the complex function @xmath128 must change by a phase that is greater than or equal to @xmath133 .
the phase can never be less than this amount because that would require the existence of off - particle poles that get taken to @xmath127 in order to diminish the strength of the on - particle zeros . by definition ,
the phase change of @xmath128 under the above procedure is @xmath134 .
therefore , the condition that the wavefunction have no poles leads directly to the following concavity condition on the integers @xmath82 : @xmath135 in cases where all of the zeros are located on the particles and there are no off - particle zeros , the above inequality is saturated .
this occurs in the laughlin states @xmath136 , and their multilayer abelian generalizations , the halperin states @xmath137    in the following we will rewrite the concave condition as @xmath138 @xmath139      the cluster condition is a way to associate some kind of grading to the polynomials that is physically meaningful . let @xmath140 for @xmath13 be a set of vectors that generate an @xmath14-dimensional lattice , where @xmath14 is , as before , the number of flavors of particles ( or the number of layers ) .
the cluster condition states that the concave condition is saturated , _
i.e. @xmath141 if either @xmath64 , @xmath83 , or @xmath142 lie on the lattice generated by @xmath143
. that is , if either @xmath64 , @xmath83 , or @xmath142 can be written as a linear combination with integer coefficients of the vectors @xmath144 .
this means that a derived polynomial containing a @xmath145 cluster is non - zero unless @xmath146 coincides with the coordinates of another cluster ; viewed as a function of the single variable @xmath146 , it has no off - particle zeros .
a consequence of this is that if all of the particles are fused to form @xmath147-clusters , then the resulting derived polynomial has the laughlin - halperin form ( see ( [ halperin ] ) ) and there are no off - particle zeros . _
the single - layer read - rezayi @xmath148 parafermion wave functions satisfy an @xmath28-cluster condition and they are exact ground states of hamiltonians with @xmath149-body interactions . for a fixed filling fraction , as @xmath28 increases , the number of topologically distinct quasiparticles , the ground state degeneracy on higher genus surfaces and the complexity of interactions necessary to realize the state all increase .
this suggests that the energy gap typically decreases with increasing @xmath28 .
wave functions that do not obey a cluster condition can be thought of as having infinite @xmath28 and are not expected to correspond to gapped phases .
this intuition also comes from the cft approach to fqh wave functions ; infinite @xmath28 corresponds to an irrational conformal field theory , which does not yield a finite number of quasiparticles and a finite ground state degeneracy on the torus . in the multilayer case
, we may use the volume of the unit cell spanned by @xmath140 as one way to measure the complexity of a given fqh state .
the cluster condition is extremely powerful and simplifying because it allows us to determine the entire pattern - of - zeros sequence from knowledge of a `` small '' number of them . to see how this works ,
first observe using ( [ clustercond ] ) that @xmath150 where we have defined the matrix @xmath151 .
so for any vector @xmath152 , where @xmath153 is an integer and @xmath154 , we have : @xmath155 the above equations imply @xmath156 where we have also defined the matrix @xmath157 .    in terms of the sequence @xmath158
, this implies that for @xmath152 , where @xmath153 is an integer , @xmath159 therefore , all of the integers @xmath30 are specified by the points @xmath160 within the unit cell spanned by @xmath144 .    in terms of the hoo squence @xmath106
, we have : @xmath161 finally , note that since @xmath162 is an even integer , we have : @xmath163      the density profiles of the single - particle states in the lowest landau level , @xmath164 , are in the shape of a ring , with a peak at a radius @xmath165 , where @xmath166 is the magnetic length .
such a wave function has an angular momentum @xmath29 .
when many of these orbitals are filled by particles , the total wavefunction will describe a uniform , rotationally symmetric state that goes to zero at a radius @xmath167 , where @xmath168 refers to the filled orbital with maximum angular momentum .
therefore , a given quantum hall wavefunction will describe a qh droplet of area @xmath169 , and @xmath168 is given by the maximum power of @xmath44 ( or @xmath9 for any other fixed @xmath17 ) .
@xmath168 is also equal to the number of flux quanta , @xmath170 .
an important constraint on the multilayer quantum hall wavefunctions is that they must describe systems in which each layer occupies the same area , up to small corrections .
the requirement that each layer occupies exactly the same area amounts to the requirement that each layer has exactly the same number of flux quanta , @xmath171 .
however , it is reasonable to include states in which different layers occupy equal areas only up to @xmath172 corrections , where @xmath173 is the number of electrons .
such a requirement of approximately equal area layers is summarized in the following equation : @xmath174 we wish to see how this condition translates into a condition on the pattern of zeros .
the conditions are slightly different depending on whether we ultimately want to characterize gapped fqh phases of fermions or bosons .
if we are interested in fermionic phases , we require that @xmath175 be a valid fqh wave function of fermions , which does not require that @xmath176 be a valid fqh wave function of bosons . in what follows we
will explicitly analyze the bosonic case , where we require @xmath177 to be a valid multilayer fqh wave function of bosons .
@xmath114 is equal to the maximal power of @xmath72 ; for the boson wavefunction @xmath10 , this is given by @xmath178 where recall @xmath179 is the number of particles of type @xmath27 .
using the cluster condition , we find @xmath180 where we have set @xmath181 .
requiring ( [ equalareacond1 ] ) , we obtain the following condition on the pattern of zeros : @xmath182 this can be seen most easily by ignoring the @xmath172 terms in ( [ nphieqn1 ] ) , taking @xmath183 , and inverting the result to obtain @xmath184 , which must be nonnegative . from this analysis
, we learn that if @xmath185 is not invertible , then the pattern of zeros can not fix the ratio of particles @xmath186 in the different layers .
therefore the corresponding fqh state has a gapless mode corresponding to the relative density fluctuations between the different layers .    as a simple example of this analysis , consider the @xmath187 halperin bilayer state , which is known to have a gapless density mode and for which @xmath188 , which is not invertible . a macroscopic number of particles can freely go from one layer to the other without changing the area of the quantum hall droplets , signalling the existence of a gapless relative density mode
.    inverting ( [ nphieqn1 ] ) yields @xmath189 from which we can read off the filling fraction in each layer : @xmath190 the total filling fraction is the sum of the filling fraction of each layer : @xmath191 .    for fermions , ( [ nphipoz ] )
is modified to @xmath192 due to the extra factor @xmath193 in @xmath194 .
note that @xmath195 still describes the pattern of zeros of the symmetric polynomial @xmath10 .
the result for fermions is therefore @xmath196 where @xmath197 is the @xmath198 identity matrix .
if @xmath199 is not invertible in the fermionic case , then there are gapless relative density modes , which is why the filling fraction in each layer becomes undefined .
consider a multilayer quantum state with @xmath179 particles in the @xmath95 layer .
we want to put the quantum state on a sphere with @xmath200 flux quanta in the @xmath95 layer .
we would like to know for which set of @xmath200 can the quantum hall state completely fill the sphere ? naively , one may expect @xmath200 and @xmath179 are related by the filling fraction in each layer @xmath201 .
however the precise relation between the number of flux quanta and the number of electrons includes a shift , @xmath202 where @xmath43 is of order 1 in the large @xmath173 limit ( see eqn .
( [ nphieqn2 ] ) ) .    more precisely , completely filling the sphere
means that the quantum hall state is rotationally invariant with zero total angular momentum .
using , we find that , for a bosonic fqh state characterized by @xmath30 with @xmath179 particles and @xmath203 flux quanta in the @xmath95 layer , the maximum total angular momentum is given by @xmath204 where @xmath205 and @xmath206 .
for a fermionic fqh state characterized by @xmath30 with @xmath179 particles and @xmath203 flux quanta in the @xmath95 layer , the maximum total angular momentum is given by @xmath207 in the above , @xmath208 must satisfy @xmath209 in order for the wave function to fit into each layer .    completely filling the sphere
requires that @xmath210 and @xmath211 .
we see that @xmath212 and @xmath208 must satisfy @xmath213 which implies ( for both bosons and fermions ) @xmath214 if a given @xmath212 does not satisfy , then the corresponding quantum hall state ( with @xmath179 particles on the @xmath95 layer ) can not completely fill the sphere . for @xmath212
that satisfies , the corresponding quantum hall state can completely fill the sphere and has zero total angular momentum if @xmath208 is given by .
can generally be satisfied only if @xmath212 lies on the lattice spanned by @xmath144 .
we would like to remark that it is easy to have different numbers of flux quanta on different layers in numerical calculations .
the pattern of @xmath215 where the quantum hall state has zero total angular momentum on the sphere can be used as a fingerprint to identify different quantum hall states through numerical calculations ( for examples , see tables [ fluxnotable23 ] and [ fluxnotable45 ] ) .
the analysis of the single - layer case in has suggested an additional condition : @xmath217 there , it was found that allowing @xmath218 allows for certain pattern - of - zeros sequences that either do not correspond to single - valued wavefunctions ( such as the square root of the pfaffian ) or could not correspond to translationally invariant wavefunctions .
it was suggested that one way to rule out such possibilities is to impose ( [ delta3even ] ) .
how should this condition be generalized to the multilayer situation ?
one natural generalization is to impose @xmath219 for all @xmath64 , @xmath83 , and @xmath142 .
however , we find that this condition is too restrictive .
it rules out certain known fqh wavefunctions , such as the @xmath220 non - abelian spin singlet states.@xcite the need to relax this condition while still having it remain compatible with the single - layer situation suggests that we should impose ( [ delta3even ] ) only for choices of @xmath64 , @xmath83 , @xmath142 that are collinear through the origin .    while it was found that allowing @xmath221 allows for pattern - of - zeros sequences that do not seem to correspond to valid translationally invariant , single - valued wavefunctions , there are known cases of cfts with @xmath221 that do seem to yield translationally invariant , single - valued wavefunctions .
one such example is the so - called gaffnian wavefunction , which has @xmath222 and which can be constructed using the minimal model cft @xmath223.@xcite this cft however is non - unitary .
it has been suggested that fqh wavefunctions constructed using non - unitary cfts correspond to gapless phases;@xcite whether this is always necessarily the case is currently an important open question in fqh theory .
these considerations suggest that in order to restrict ourselves to pattern - of - zeros sequences that have a corresponding unitary cft , we should impose @xmath224 for those @xmath64 , @xmath83 , and @xmath142 that are collinear through the origin . in our search for pattern - of - zeros solutions , we will impose this condition and analyze the resulting states . the precise connection , if any , between this condition and valid fqh wavefunctions that correspond to unitary cfts remains to be clarified
given a pattern - of - zeros sequence , it is important to be able to construct a local , gapped hamiltonian whose ground state wave function has the given pattern of zeros . if this is possible , then we know that the corresponding pattern - of - zeros sequence describes a topological phase of matter .
whether this particular phase is realized in an experiment then depends on the particular types of low energy effective interactions between the electrons in the unfilled landau levels .
we can go about constructing such a hamiltonian by noticing that on a sphere , the integers @xmath30 are directly related to the angular momentum of the @xmath64-cluster . for an electron system on a sphere with @xmath225 flux quanta for
the @xmath27th layer , an electron of type @xmath27 will carry an angular momentum @xmath226 . for an @xmath227-cluster ,
the maximum angular momentum is therefore @xmath228 .
however , for a polynomial @xmath229 described by a pattern of zeros @xmath89 , the maximum allowed angular momentum of the @xmath64 cluster is only @xmath230 .
the pattern of zeros forbids the appearance of angular momentum @xmath231 for any @xmath227-clusters in @xmath232 .
such a condition can be easily enforced by writing the hamiltonian as a sum of projection operators , @xmath233 .
let @xmath234 be a projection operator that acts on the @xmath64-cluster hilbert space .
@xmath233 projects onto the subspace of @xmath227-clusters that have total angular momenta greater than @xmath235 .
now consider the hamiltonian @xmath236 where @xmath237 sums over all of the @xmath227-clusters for a fixed @xmath64 .
the wave function described by @xmath195 will clearly be a zero - energy ground state of the above @xmath238 .
in many cases , there is only one unique ground state wave function with minimal total angular momentum , but in general there can be many independent polynomials with the same pattern of zeros .
in such a situation , the hamiltonian would need to be modified further to select for a particular polynomial with the given pattern of zeros .    in order for the above hamiltonian to be local ,
@xmath239 must be limited to a small , finite number of @xmath64-clusters .
but as a result , we can not guarantee in general that the ground state wave functions will all be described by the sequence @xmath89 for every @xmath64 , or even that they will obey the cluster condition .
in many of the known cases , such as the laughlin , moore - read , and read - rezayi wave functions , the sum over @xmath64-clusters can indeed be terminated after the first few clusters while still yielding a unique zero energy ground state wave function with minimal angular momentum which is described by @xmath195 .
the above construction for @xmath240 should therefore be viewed as a starting point for constructing an ideal hamiltonian that is local , gapped , and whose unique minimal angular momentum ground state wave function is described by @xmath158 . in some of the simplest cases ,
we know that this construction suffices .
we have found that the polynomials @xmath241 , @xmath242 , that may correspond to stable fqh states are described by the following data : @xmath243 where @xmath28 and @xmath29 are @xmath198 matrices with integer entries that satisfy @xmath244 the above implies that @xmath185 is a symmetric matrix .
furthermore , for the pattern of zeros to fix the relative densities of particles in each layer , we have @xmath245 otherwise , there are gapless relative density fluctuations .
we also have @xmath246 the total filling fraction is @xmath247 .
note that the @xmath158 need to be specified only for values of @xmath64 that are contained in the unit cell spanned by @xmath248 , where @xmath249 corresponds to the @xmath95 row of the matrix @xmath28 : @xmath250 .
using the cluster condition , we can determine @xmath30 for all @xmath64 simply from @xmath29 , the fact that the @xmath91 , and from the values of @xmath251 for which @xmath83 lies in the unit cell spanned by @xmath144 : @xmath252 where @xmath145 .
@xmath30 must satisfy : @xmath253 where @xmath254 @xmath255
finally we impose @xmath256 for those @xmath64 , @xmath83 , and @xmath142 that are collinear through the origin .
consider two polynomials of @xmath257-cluster form , @xmath10 and @xmath258 , and consider their product : @xmath259 .
the pattern of zeros of @xmath260 is the sum of the pattern of zeros of @xmath10 and @xmath258 : @xmath261 similarly , the data in terms of @xmath29 and @xmath262 are also additive .
note that all of the conditions that we impose on the pattern - of - zeros are linear .
thus the pattern of zeros of @xmath260 is also valid .
however , notice that the condition for filling the sphere is not linear .
there may be two fqh wavefunctions @xmath10 and @xmath258 that can fill the sphere but whose product @xmath260 can not fill the sphere .
thus , we can divide the pattern - of - zeros solutions into primitive and non - primitive solutions .
primitive solutions are solutions that can not be written as a sum of two other solutions .
the pattern - of - zeros approach is closely related to the conformal field theory approach to constructing fqh wavefunctions . in the cft approach ,
the symmetric polynomial @xmath263 that describes a multilayer fqh state can be written as a correlation function of a set of electron operators @xmath264 in a cft:@xcite @xmath265 the operators @xmath264 are written in the form : @xmath266 where @xmath267 is a vertex operator in a @xmath268 cft .
it has scaling dimension @xmath269 .
@xmath270 is a simple current operator ; that is , it satisfies the following fusion relation : @xmath271 this abelian fusion rule is the cft version of the unique - fusion condition .
the cluster condition implies that @xmath272 satisfies @xmath273 where @xmath250 .
an @xmath64-cluster of electrons will be described by the operator @xmath274 thus we see that the cluster condition implies that an @xmath249 cluster is described by a vertex operator @xmath275 .
if all of the particles are grouped into @xmath28-clusters , then the corresponding derived polynomial will be just a correlation function of vertex operators in a @xmath268 theory , which will have no off - particle zeros and will be of the laughlin - halperin form .
let us denote the scaling dimension of the operator @xmath276 as @xmath277 where @xmath278 is the scaling dimension of the simple current @xmath279 and @xmath280 is the scaling dimension of the vertex operator .
notice that since @xmath281 , the simple current scaling dimensions satisfy @xmath282 .
the scaling dimension of the gaussian part is given by @xmath283 the pattern of zeros are related to the scaling dimensions through the relation @xmath284 this allows us to obtain the scaling dimensions from the pattern of zeros . using the cluster condition , some algebra shows that @xmath285 and
so ( [ scalingdimga ] ) becomes @xmath286 the scaling dimensions of the simple - current part can also be determined from the pattern of zeros by using the fact that @xmath287 , @xmath288 , and applying ( [ pozsc ] ) iteratively .
this yields : @xmath289 multiplying both sides by @xmath290 gives @xmath291 in a similar manner , one can obtain @xmath292 which determines @xmath278 in terms of the pattern - of - zeros .    note that the correlation function of the gaussian part is , leaving the background charge implicit , @xmath293 thus the fqh wave function is of the form @xmath294 where @xmath295 @xmath296 arises from the correlation function of the simple current sector and is the `` non - abelian part '' of the wavefunction .    in this manner , each pattern - of - zeros solution corresponds to the current algebra of a rational cft . the connection between the pattern of zeros approach and the cft approach can be thought of in the following way .
the pattern of zeros describes the essential properties of the cft that yield valid fqh wavefunctions .
so in order to classify ideal fqh wavefunctions , one can bypass the cft altogether and go directly to the heart of the matter : characterizing the allowed pattern - of - zeros solutions . furthermore , since each pattern - of - zeros solution corrresponds to a cft , the pattern - of - zeros classification can be viewed as a classification of the allowed cfts that can be used to construct fqh wave functions .    in this formulation
, the pattern of zeros classifies all those ideal fqh wavefunctions that can be formulated as a correlation function of _ conformal primary fields .
there are also many fqh wavefunctions , such as the hierarchy states and the jain series , that can not be written in this way .
these wavefunctions are outside of the pattern - of - zeros classification .
however they may be closely related to the pattern - of - zeros construction .
we comment on this connection elsewhere . _      using ( [ scalingdimsc ] )
, we can derive a formula for @xmath30 in terms of @xmath297 , @xmath28 , and @xmath29 : @xmath298 thus there is a one - to - one correspondence between the simple - current scaling dimensions and the sequence @xmath195 .
this means that there is yet another way to label the pattern of zeros .
earlier , we found that one convenient labeling of the pattern of zeros is with an @xmath198 matrix @xmath28 , an @xmath198 matrix @xmath29 , and the value of the non - negative integers @xmath30 for @xmath64 lying inside the unit cell spanned by the rows of @xmath28 .
an alternative , equivalent labeling of the pattern of zeros is by specifying the following data : @xmath299 for @xmath64 lying inside the unit cell spanned by @xmath28 .
this labeling is convenient because it makes close contact with the corresponding cft description ; if @xmath300 for some @xmath64 , then the cft has a non - trivial simple - current structure and therefore generally also has some form of non - abelian statistics .      in the single - layer case , many of the pattern - of - zeros solutions were found to be closely related to the @xmath302 parafermion cfts of zamalodchikov and fateev.@xcite what are the relevant cfts to expect in the multilayer case ?
the answer is that some of the corresponding cfts in the multilayer case will be closely related to the @xmath303 parafermion cfts that were constructed by gepner , where @xmath0 is a simple affine lie algebra at level @xmath304 and @xmath305 is the rank of the lie algebra @xmath306.@xcite the case @xmath307 is equivalent to the @xmath302 parafermion cfts of zamalodchikov and fateev .
the simple - current algebra of the @xmath301 parafermion cft has the following structure .
for every element @xmath308 of the @xmath305-dimensional root lattice of @xmath306 , associate a simple current operator @xmath309 .
the simple - currents will have the fusion rules @xmath310 furthermore , @xmath311 if @xmath312 is an element of @xmath304 times the long root lattice of @xmath306 .
the scaling dimension of @xmath309 is given by @xmath313 where @xmath314 is an integer equal to the minimum number of roots from which @xmath315 is composed .
the inner product @xmath316 is defined with respect to the quadratic form matrix of @xmath306 .
if we are considering quantum hall states with @xmath14 layers , then we would expect to see the appearance of these parafermion cfts with rank @xmath317 .
therefore in the bilayer case , one class of states that we expect to see should be related to @xmath318 parafermion cfts where @xmath306 is a simple lie algebra of rank 2 .
there are only three simple lie algebras of rank 2 : @xmath319 , @xmath320 , and @xmath321 . of these ,
only @xmath319 is simply laced , so the long root lattice is the same as the root lattice .
this means that for the pattern - of - zeros solutions that correspond to @xmath322 , we expect to see solutions that correspond to @xmath323 cfts . in appendix [ su3para ]
we will describe the @xmath220 cft in more detail .
the parafermion cfts for @xmath324 and @xmath325 , on the other hand , are more complicated because the long root lattice is different from the root lattice . for example , @xmath326 cfts will generically be relevant for @xmath327 while @xmath328 cfts will generically be relevant for @xmath329 .
in this section we examine explicitly several simple bilayer pattern of zeros solutions .
we first fix the cluster structure @xmath28 to have a simple form .
then we try to find all the solutions @xmath158 that satisfy the conditions listed in section [ summary ] .
the simplest non - abelian states can be obtained from the simplest cluster structures @xmath330 and @xmath331 .
these are the simplest non - abelian generalizations of the bilayer abelian laughlin - halperin states .    note that by definition the ordering of the rows in @xmath28 is arbitrary ; we choose it so that @xmath332 .
interchanging the layers yields the same physical system but corresponds to interchanging @xmath333 with @xmath334 and @xmath335 with @xmath336 , so two matrices @xmath28 and @xmath337 that are related by such an interchange are regarded as equivalent .
to list the solutions for those simplest cases , we may use some known cfts to construct the appropriate simple - current algebra that corresponds to the pattern - of - zeros solutions . using this known cft
, we can then write the wavefunction explicitly . however , the wave function that we write down may not be unique in some cases ; there may be several independent polynomials that have the same pattern of zeros .
this corresponds to there being several distinct cfts whose simple - current algebra possesses the same pattern of zeros . in the following examples
, we will make this choice when necessary so that we can explicitly write down a wave function with a certain pattern of zeros .
these states are all abelian and correspond to the halperin @xmath339 states .
the unit cell spanned by the @xmath249 consists only of the points @xmath340 and @xmath341 ; by translation invariance @xmath342 . thus these states are completely characterized by the matrix @xmath29 and are of the form @xmath343 the @xmath29 matrix here is exactly the @xmath344 matrix that describes all abelian fqh states.@xcite we also have @xmath345 for all of these states  the cft has no simple - current part and corresponds to a @xmath346 gaussian cft .
although the single - layer hierarchy states , such as the @xmath347 hierarchy state , do not have ideal single - layer wave functions , there are ideal multilayer states that have the same topological orders as those single - layer hierarchy states .
for example , the topological order in the @xmath347 hierarchy state is described by the @xmath344-matrix @xmath348 ( in the symmetric basis).@xcite such a topological order can be represented by the ideal bilayer state with @xmath338 and @xmath349 .
so although the pattern - of - zeros construction does not directly classify those single - layer hierarchy states , their topological orders can still be described by the pattern - of - zeros approach .
this choice of @xmath28 requires that the electron operators in the cft must take the form @xmath351 where @xmath352 and @xmath272 has scaling dimension that is integer or half - integer .
this latter fact can be obtained from the condition @xmath353 .
one general way of constructing such a simple - current operator is by expressing it as a set of majorana fermions from several copies of the ising cft : @xmath354 where @xmath355 is the majorana fermion from the @xmath17th ising cft .
such an operator has scaling dimension @xmath356 and gives rise to the following fqh wavefunction : @xmath357 @xmath358 represents the coordinates in both layers : @xmath359 @xmath360 is defined in ( [ phiga ] ) .
note that the simple - current algebra in this case implies @xmath361 .
that is , these states can all also be viewed as satisfying the cluster condition for @xmath362 , but with a different choice @xmath29 . for every pattern - of - zeros solution found here
, there is an equivalent one for @xmath363 .
later , we will list the solutions with the @xmath364 cluster structure .
some of those solutions actually have @xmath365 .
[ [ n - left - protectbeginmatrix-1 - 1 - 0 - 1-protectendmatrix - right ] ]      here , the electron operators must be of the form @xmath367 the fact that @xmath368 however also forces @xmath369 .
thus in fact the corresponding cfts do not have a simple - current part ; they all correspond to a @xmath346 guassian cft .
all of these states are therefore abelian and correspond to the halperin @xmath339 states .
thus , all pattern - of - zeros solutions where @xmath370 actually also satisfy the cluster condition for @xmath371 .    using the cluster condition and the fact that @xmath372 , it is easy to see that the pattern of zeros for these states are completely characterized by the matrix @xmath29 and therefore that these states are identical to the halperin states .    [
[ n - left - protectbeginmatrix-1 - 0 - 0 - 2-protectendmatrix - right ] ]      in this case , the electron operators in the cft must take the form @xmath374 where @xmath352 . from @xmath375
, it follows that the scaling dimension of @xmath272 is integer or half - integer .
this means that we can in general write it as a product of @xmath376 majorana fermion operators from @xmath376 independent copies of the ising cft .
thus the pattern - of - zeros solution for this choice of @xmath28 includes only two classes of states : the abelian halperin states and the following version of the pfaffian states : @xmath377 such a state spontaneously breaks the discrete @xmath378 symmetry associated with interchanging the two layers .      for this choice of @xmath28 ,
the corresponding cfts have two simple current operators , @xmath380 and @xmath381 , which each square to the identity : @xmath382 .
thus there are a total of three distinct primary currents : @xmath381 , @xmath380 , and @xmath383 .
the pattern of zeros can be fully specified by specifying the scaling dimensions of these simple current operators , @xmath384 , @xmath385 , and @xmath386 , and the matrix @xmath29 , which specifies the @xmath346 part of the electron operators in the cft .
applying @xmath387 for this choice of @xmath28 shows that the simple - current scaling dimensions are all integer or half integer : @xmath388 .
one general way of encoding these fusion relations and the associated scaling dimensions is to write the electron operators in the cft using simple - current operators that consist of multiple copies of the @xmath220 parafermion cft .
the most general choice for the electron operators can be written in the form @xmath389 where @xmath390 is a simple current from the @xmath376th copy of the @xmath220 parafermion cft and @xmath391 and @xmath392 refer to either @xmath340 , @xmath341 , or @xmath393 .
some explicit forms for such correlators of simple - current operators in the @xmath220 parafermion cft were discussed in .
computing these correlation functions provides one way  not necessarily unique  of constructing a wave function with the desired pattern of zeros .    these pattern - of - zeros solutions
can naturally be grouped into two distinct classes . in the first class ,
@xmath394 and @xmath395 contain only one kind of simple - current , either @xmath380 , @xmath381 , or @xmath396 , from each copy of the @xmath220 cft .
for example , if @xmath397 appears in @xmath394 , then @xmath395 can not contain @xmath398 or @xmath399 .
in such a situation , we can think of @xmath400 as being the majorana fermion from the ising cft .
this means that these states can be written in terms of multiple copies of the ising cft .
this class of multilayer states can be written by choosing @xmath401 where @xmath402 is the majorana fermion from the @xmath376th copy of the ising cft .
the wave function for this class of states is therefore :    @xmath403    @xmath358 represents the coordinates in both layers : @xmath359    this is the simplest generalization of the single - layer non - abelian states to a class of non - abelian bilayer states : the interlayer pfaffian states .
the simplest version of this , with @xmath404 , is closely related to ( and identical to for certain choices of @xmath29 ) the spin - charge separated non - abelian spin singlet wavefunction proposed in .
the second class of states can not be written in terms of multiple copies of the ising cft
the full @xmath220 parafermion cft is necessary .
the first class of states , which can be written only using the ising cft , have the property that their pattern of zeros satisfies @xmath405 for all choices of @xmath64 , @xmath83 , and @xmath142 .
the second class of states , for which the full @xmath220 cft is necessary , satisfies @xmath406 for certain choices of @xmath64 , @xmath83 , and @xmath142 .
let us compare the simple - current algebra of the @xmath220 parafermion cft to what one would obtain using two copies of the ising cft .
if we used two copies of the ising cft , we could have @xmath407 with scaling dimension @xmath408 , @xmath409 with scaling dimension @xmath408 , and @xmath410 with scaling dimension @xmath411 .
this satisfies @xmath412 . on the other hand , in the @xmath220 parafermion cft
, the only difference is that @xmath396 also has scaling dimension @xmath408 .
thus in this latter theory , two fermions combine to give another fermion .
this yields @xmath413 .
the fact that there are valid single - valued translationally invariant fqh wavefunctions that arise from unitary cfts and that have @xmath406 for certain choices of @xmath64 , @xmath83 , and @xmath142 suggests ( see section [ delta3evensec ] ) that we should impose @xmath414 not in general but only if @xmath64 , @xmath83 , and @xmath142 are collinear through the origin .
let us examine the pattern of zeros for a few of the simplest examples of these non - abelian bilayer states .
there is a fermionic @xmath1 state with @xmath415 this is the pattern of zeros for the interlayer pfaffian state , which is of the form @xmath416 we use the notation @xmath417 there are also fermionic states at @xmath418 and @xmath419 .
these have the following pattern of zeros : @xmath420 @xmath421 the state at @xmath419 is the non - abelian spin singlet state that was proposed in .
note once again that the pattern of zeros @xmath29 and @xmath158 refer to the pattern of zeros of the symmetric polynomial , @xmath422 .
in single - layer quantum hall samples , a quantum hall plateau is seen at @xmath423 , but not at @xmath424 .
the reason is that even though in all of these cases there is a single half - filled landau level , the existence of the two filled extra landau levels modifies the effective interactions between the electrons in the unfilled level . in the @xmath423 case
, numerical calculations suggest that the these effective interactions are modified in such a way that a non - abelian quantum hall state may be realized .
experiments on multicomponent quantum hall systems should be able to probe an even wider variety of regimes with distinct effective interactions .
for example , for a two - component fqh system , we can study systems in which the spin degree of freedom is present , two - dimensional electron systems with two quantum wells , wide single - layer systems in which the electrons spontaneously form a double - layer system due to coulomb repulsion , or systems in which there may be two valleys for the free quasiparticle spectrum ( such as in graphene or sige heterostructures ) , etc .
in many of these cases , experimentalists can also tune to some extent the degree of correlation between the two components .
for example , in double layer systems , application of a parallel magnetic field can tune the tunneling and correlation between the layers .
there may also be some degree of tunability in the relative densities between the two components in addition to being able to probe fqh states with different numbers of filled landau levels . with this greatly increased amount of variability and tunability in the effective interactions between electrons in the unfilled landau levels , it is possible that a non - abelian state can be realized in a two - component quantum hall system .
since the pattern of zeros provides a systematic classification and characterization of a wide variety of quantum hall states , it provides us with a general sense of how all of the non - abelian bilayer states are related and which ones are simpler than other ones .
just as we know that the single - layer pfaffian quantum hall state is the simplest non - abelian generalization of the laughlin states , we can determine the simplest non - abelian generalization of the halperin bilayer states and therefore single out some of the possibilities that may be experimentally viable .    in
, we have given an overview of some of the simplest non - abelian bilayer states that we find and that occur at filling fractions at which experiments on two - component fqh systems have already observed incompressible states . here we briefly summarize that discussion and supplement details of the calculations of various topological properties of the candidate states .
experiments have so far observed fqh plateaus in two - component systems at @xmath1 , @xmath2 , @xmath3 , @xmath425 , @xmath426 , @xmath427 , @xmath428 , etc .
@xcite in some cases , these plateaus have been observed in both bilayer and spin - unpolarized single - layer systems , while in others , the plateau has only been observed in one of them .
at all of these filling fractions , there exists also one ( or several ) candidate abelian phase(s ) ; in most cases , it is assumed that these plateaus are described by one of the abelian phases .
however , the pattern - of - zeros construction also yields many simple non - abelian states at these filling fractions . in some situations
, we expect the non - abelian states to be good candidate states .
there are at least two dimensionless quantities that are important determining factors for which fqh state is realized .
the first parameter is @xmath429 , where @xmath430 is the potential for interlayer repulsion and @xmath431 is the potential for intralayer repulsion .
the second parameter is @xmath432 , where @xmath433 is the interlayer hopping amplitude . in the limit @xmath434 and @xmath435
, the system will be a fqh state that consists of two independent single - layer fqh states in each layer . in the limit @xmath436 and @xmath434 ,
a single - layer fqh state may be observed .
but if we keep @xmath435 and increase @xmath437 from @xmath438 , then the fqh state formed by two independent single - layer fqh states in each layer must undergo a phase transition into either a compressible phase or a new incompressible state . in the latter case ,
an abelian hierarchy state ( such as a bilayer composite fermion state ) may form , which would in most cases be a state described by a @xmath439 or more complicated @xmath344-matrix and would have four or more edge modes .
the other possibility is that a non - abelian two - component state may form .
the pattern - of - zeros construction yields non - abelian two - component wave functions that have zeros when particles from different layers approach each other , indicating that they can accomodate situations in which @xmath440 . additionally , these states generally have less than @xmath441 edge modes ; if we use the number of edge modes as a measure of the complexity of the state , then the non - abelian states are simpler and may therefore be realized experimentally .    at @xmath1 , experiments on wide single quantum wells have observed a phase transition from a bilayer to single - layer state while experiments on single - layer systems have seen a phase transition from a spin - polarized to a spin - unpolarized state . in the limit @xmath442 and @xmath435
, the system should be in the @xmath443 state . as @xmath437 is increased while @xmath435 , one possibility is the @xmath444 state .
this wave function appears unphysical , because it has higher order zeros as particles from different layers approach each other than particles from within the same layer .
another wave function , which has the same topological order as @xmath444 , is a spin - singlet composite fermion state.@xcite there are two other plausible non - abelian states in this situation .
one is the following interlayer pfaffian state ( see ( [ interpfpoz ] ) and ( [ interpfwf ] ) ): @xmath445 the other is the following intralayer pfaffian state @xmath446 which has even higher order zeros as particles from different layers approach each other .
@xmath447 has a cluster structure @xmath448 while @xmath449 has a cluster structure @xmath450 . at @xmath451
there are also two single - layer possibilities that may be realized as @xmath452 is increased .
these are the particle - hole conjugate of the @xmath453 laughlin state and the @xmath454 parafermion read - rezayi state .    .
[ proposed ] proposed explanations for incompressible states at experimentally relevant filling fractions , @xmath1 , @xmath2 , @xmath3 , and @xmath428 , in two - component fqh systems .
the bilayer composite fermion state @xmath455 @xcite refers to the state @xmath456 , where @xmath457 is a single layer composite fermion state at filling fraction @xmath458 . for @xmath459 ,
we have taken the single layer @xmath460 state to be the particle - hole conjugate of the laughlin state .
@xmath461 indicates that there are @xmath462 right - moving edge modes and @xmath463 left - moving edge modes .
see appendix [ calculations ] for details of how to calculate the number of edge modes and the shift @xmath43 . [
cols="^,<,^,^",options="header " , ]     in table [ proposed ] we summarize some of the filling fractions at which incompressible states have been experimentally observed in two - component fqh systems . for each filling fraction
we list some of the proposed wave functions that may characterize the topological order of those phases , the number of edge modes , and their respective shifts on the sphere .
we list the qasiparticles , their electric charges , and their scaling dimensions for the interlayer pfaffian state at @xmath1 in table [ ipfqpops ] . in table
[ qpdata ] , we list the quasiparticles with the minimal electric charge and their scaling dimensions @xmath464 for the non - abelian fqh states discussed in this paper [ see eqns .
( [ interpfwf ] ) , ( [ 45poz ] ) , ( [ 47poz ] ) , ( [ su3states ] ) , ( [ na14 ] ) ] .
those minimally charged quasiparticles may dominate interedge tunneling and give rise to the following i - v curve : @xmath465 in the @xmath466 limit .    in summary
, we find many simple non - abelian bilayer states that occur at experimentally observed filling fractions . for certain effective interactions among the electrons in the unfilled landau levels , these states may be more favorable than their abelian counterparts . in these cases , the non - abelian states have larger interlayer correlations and therefore may be energetically more favorable in situations in which the interlayer repulsion is comparable to the intralayer repulsion .      a useful tool for identifying fqh states in numerical studies of exact diagonalization on finite systems on a sphere is to look at what values of the shift , @xmath467 , a ground state with zero total angular momentum is found .
this then limits the possibilities of which topological phase is realized in the system to those that have that particular value of the shift .
similarly , in such numerical studies of multilayer systems , one can look for the different sets @xmath468 that yield a ground state with zero total angular momentum .
each topological phase will have its own list of @xmath469 that let it fill the sphere ; analyzing this can be a useful way of determining which topological phase is obtained numerically . in section [ shiftrisphere ] , we found conditions that @xmath212 and @xmath470 should satisfy for the fqh state to fill the sphere .     & & + & & 2 & 4 & 6 & 8 & 10 + & 2 & ( 3,3 ) & ( 7,5 ) & ( 11 , 7 ) & ( 15,9 ) & ( 19,11 ) + & 4 & ( 5,7 ) & ( 9,9 ) & ( 13,11 ) & ( 17,13 ) & ( 21,15 ) + & 6 & ( 7,11 ) & ( 11,13 ) & ( 15,15 ) & ( 19,17 ) & ( 23 , 19 ) + & 8 & ( 9,15 ) & ( 13,17 ) & ( 17,19 ) & ( 21,21 ) & ( 25,23 ) + & 10 & ( 11,19 ) & ( 15,21 ) & ( 19 , 23 ) & ( 23,25 ) & ( 27,27 ) +     & & + & & 2 & 4 & 6 & 8 & 10 + & 2 & ( 2,2 ) & ( 6,3 ) & ( 10 , 4 ) & ( 14,5 ) & ( 18,6 ) + & 4 & ( 3,6 ) & ( 7,7 ) & ( 11,8 ) & ( 15,9 ) & ( 19,10 ) + & 6 & ( 4,10 ) & ( 8,11 ) & ( 12,12 ) & ( 6 , 13 ) & ( 20 , 14 ) + & 8 & ( 5,14 ) & ( 9,15 ) & ( 13,6 ) & ( 17,17 ) & ( 21,18 ) + & 10 & ( 6,18 ) & ( 10,19 ) & ( 14 , 20 ) & ( 18,21 ) & ( 22,22 ) +    for states that have a cluster structure @xmath471 , we find that the condition ( [ nss ] ) becomes trivial as long as @xmath181 , where @xmath153 is an integer .
this means that as long as @xmath472 and @xmath473 are even and @xmath474 , @xmath475 satisfy ( [ snsne ] ) , then these states can fill the sphere . in this case
, we find that ( [ snsne ] ) reduces to the form @xmath476 where @xmath477 is a @xmath478 matrix and @xmath43 is the shift , which can be calculated using eqn .
( [ shift ] ) . the states that we have been considering are of the form @xmath479 , for which @xmath480 .
tables [ fluxnotable23 ] and [ fluxnotable45 ] lists some examples .
in this paper , we generalized the pattern - of - zeros characterization and classification of fqh states to multicomponent cases .
we found that the topological orders in a multicomponent fqh state can be characterized by the following data : a matrix @xmath28 that describes the cluster structure , a matrix @xmath29 and a sequence @xmath158 that describes the pattern of zeros .
our pattern - of - zeros characterization gives us a general quantitative view on a large class of abelian and non - abelian bilayer fqh states , which allow us to determine which states are simpler than other states .
we find some simplest non - abelian generalizations of the laughlin - halperin abelian bilayer states .
those simple non - abelian states may describe some of the bilayer / spin - unpolarized fqh states observed in experiments and numerical calculations .
this research is partially supported by nsf grant no .
in the single - component pattern - of - zeros description , there is an occupation number characterization that is a useful way to understand both the ground states and the quasiparticles in fqh states .
the generalization to multilayer states does not appear to be quite as simple or useful , but for the sake of completeness we will analyze it below .
a convenient set of single - particle basis states for particles in the lowest landau level are monomials of the form @xmath481 , for integer @xmath29 .
thus , a basis for symmetrized wavefunctions of @xmath179 particles of type @xmath27 is given by : @xmath482 where @xmath483 is a permutation of the particles of type @xmath27 , @xmath484 is an integer , and @xmath485 is a vector whose component @xmath486 is the number of particles of type @xmath27 that occupy the @xmath487 orbit .
the polynomial @xmath10 that we are interested in can be expanded in terms of these basis states as : @xmath488    now we may ask what kind of boson occupations @xmath489 will be present in the sum ( [ basisexpansion ] ) for a polynomial @xmath10 with a given pattern - of - zeros @xmath158 . to answer this question , let us set @xmath490 in @xmath491 . since @xmath10 is nonzero when @xmath490 due to translation invariance , there must be a boson occupation @xmath492 in the above sum that contains at least one boson occupying the @xmath493 orbital .
that is , there is a term in the above sum with @xmath494 . now
, suppose that we bring a second particle of the same type , @xmath495 to 0 .
the minimal power of @xmath495 in @xmath496 is @xmath497 : @xmath498 thus , among those @xmath489 that have at least one boson of type @xmath411 occupying the @xmath493 orbital , there must also be an @xmath489 that contains a second boson of type @xmath411 occupying the @xmath499 orbital where @xmath500 .
next , assume that two bosons occupy the @xmath501 and @xmath499 orbital , and bring a third particle of type @xmath411 to 0 ; the minimal power of @xmath502 is @xmath503 : @xmath504 thus , among those @xmath492 that have two type @xmath411 bosons occupying the @xmath505 and @xmath506 orbitals , there is a third boson of type @xmath411 occupying the @xmath507 orbital .
continuing in this way , we see that there must be a type @xmath411 boson occupying the orbitals @xmath508 for @xmath509 . after taking all the type @xmath411 particles to 0
, we may begin to take the type 2 particles to zero , one by one , thus obtaining that there must be a type @xmath510 boson occupying the orbitals @xmath511 . continuing this argument for bosons of every type
, we find that there must be a term in the sum ( [ basisexpansion ] ) with occupation number described by the above sequence of @xmath102 s .
however , in the above argument , we chose a particular sequence in which to take various particles to zero .
we first took all of the type 1 particles to zero one - by - one , and then all of the type 2 ones , and so on .
but we could just as well have made the argument with any sequence .
suppose that after taking @xmath17 particles to the origin , there is a @xmath512-cluster at the origin .
thus @xmath513 is a sequence that describes the order in which we take particles to the origin until all particles are at the origin . for every such sequence ,
we may make the above argument and argue that if @xmath514 , then there must be a type @xmath27 boson occupying the orbital @xmath515 .
if we enumerate all the different sequences @xmath516 by an integer @xmath437 , then by considering each @xmath437 , we see that there must be a term in the sum ( [ basisexpansion ] ) with occupation number @xmath517 .
@xmath518 would be the number of @xmath17 , along the sequence @xmath516 , for which @xmath519 .
notice that @xmath102 must be non - zero .
thus we have the following important condition on @xmath30 : @xmath520    the analysis above can be thought of in the following way .
consider an @xmath14-dimensional lattice @xmath521 , where @xmath14 is the number of layers .
at every site @xmath64 of this lattice ( @xmath522 ) we can associate the nonnegative integer @xmath30 . on each link ( @xmath64 , @xmath523 ) , of the lattice
we may also associate an integer @xmath524 .
now consider any directed path from the origin to @xmath212 ( @xmath179 is the number of particles of type @xmath27 ) , in which the sum of the coordinates of every point on the path is one larger than the sum of the coordinates of the point preceding it , and enumerate the set of these paths by @xmath437 .
to each such path we associate an occupation number sequence @xmath525 , where @xmath518 is the number of links along the path @xmath437 whose @xmath526 .
if @xmath10 has a pattern - of - zeros @xmath158 , then its basis expansion ( [ basisexpansion ] ) must contain a term with occupation number @xmath517 .
thus we may rewrite ( [ basisexpansion ] ) as @xmath527 the two sequences @xmath158 and @xmath528 contain the same information and are one - to - one labellings of each other .
however , @xmath528 is redundant in the sense that it does not need to be specified for every @xmath437 in order to reconstruct @xmath106 . the @xmath529 that appear in the second sum
characterize the subleading terms that appear when coordinates are brought together ; thus those @xmath530 correspond to sequences @xmath531 where @xmath532 .    in the single - layer case
, @xmath533 naturally defined an occupation number sequence @xmath534 , which also described the fqh state in the thin - cylinder limit . in the multi - component generalization
, we have @xmath535 , which seems to admit no simple generalization of the above occupation number sequence .
instead , one has such occupation number distributions for a large number of sequences which we enumerated above by @xmath437 .
we have not analyzed on general grounds which particular sequences contribute the most weight to the wavefunction in the thin cylinder limit .
some of the simplest non - abelian bilayer states are closely related to @xmath220 parafermion cft .
this cft has central charge @xmath536 and has three simple currents , @xmath537 , @xmath538 , and @xmath539 , all of which square to the identity and have scaling dimension @xmath408 .
there are four other primary fields , which are associated with the fundamental representation of @xmath319 .
their scaling dimensions are listed in table [ su3paraf ] .
the fusion rules for these fields all follow from the following fusion rule : @xmath540     cft operator & scaling dimension + @xmath541 & @xmath542 + @xmath543 & @xmath542 + @xmath544 & @xmath545 + @xmath546 & @xmath542 +
the total number of edge modes is equal to the central charge of the corresponding cft for the states that are described by the pattern of zeros . for the hierarchy states ,
the number of edge modes is given by the rank of the @xmath344-matrix .
furthermore , in the latter case , the number of right ( left ) -moving edge modes is given by the number positive ( negative ) eigenvalues of the @xmath344-matrix .
the interlayer pfaffian states are described by a cft that consists of the ising cft , with @xmath547 , and two scalar boson cfts , each with @xmath548 .
thus the number of edge modes for the interlayer pfaffian is @xmath549 .
the central charge of the @xmath220 parafermion cft is @xmath536 ; the two - component fqh states based on this are described by the @xmath220 theory and two scalar bosons , for a total of @xmath551 edge modes .    for the @xmath552 states ,
the @xmath344-matrix is @xmath553 .
these states have 2 edge modes ; if @xmath554 , all edge modes move in the same direction ; if @xmath555 , then there is one right - moving and one left - moving edge mode .    for the states @xmath556 , the @xmath344 matrix is @xmath557 where @xmath558 is the @xmath344-matrix in the hierarchical basis of the state @xmath559 . for the @xmath560 state ,
@xmath558 in the hierarchical basis is @xmath561 for the @xmath562 state , @xmath558 in the hierarchical basis is @xmath563 for the @xmath451 p - h conjugate of the @xmath564 laughlin state , @xmath558 in the hierarchical basis is @xmath565 using this , we find that the p - h conjugate of the @xmath564 laughlin state has edge modes @xmath566 , @xmath567 has edge modes @xmath568 , and @xmath569 has @xmath570 .        now consider the bilayer composite fermion state @xmath556 : @xmath574 let @xmath575 be the maximum power of @xmath44 in @xmath576 .
it satisfies : @xmath577 where @xmath578 is the shift of the state @xmath559 .
the factor @xmath579 increases the power of @xmath44 by @xmath580 .
thus the maximum power of @xmath44 in @xmath581 is @xmath474 : @xmath582 in our cases , @xmath583 , and the number of flux quanta is the same in each layer , so @xmath584 @xmath585 is the filling fraction in one layer .
thus we see that the shift of @xmath581 is also @xmath578 .    for @xmath459
, we take @xmath586 to be the particle - hole conjugate of the @xmath564 laughlin state .
@xmath586 has shift @xmath587 .
thus @xmath459 also has shift @xmath572 .
@xmath588 has shift @xmath441 , because @xmath589 has shift @xmath441 .
the electron operators for the @xmath220 fqh states that we discuss are of the form : @xmath593 where @xmath594 , @xmath408 , and @xmath595 for the @xmath418 , @xmath3 and @xmath425 states , respectively .
the quasiparticle operators with minimal total charge are of the form : @xmath596 and have scaling dimension @xmath597 .
the total charge of the quasiparticle is @xmath598 , @xmath599 , and @xmath600 for the @xmath418 , @xmath3 and @xmath425 states , respectively .
@xmath601 , @xmath411 , and @xmath599 , respectively , for these states .