Creator

Date

2001-04-02

Description

Comment: 9 pages, 2 figures

A mapping is developed between the quantum Hall plateau transition and

two-dimensional self-interacting lattice polymers. This mapping is exact in the

classical percolation limit of the plateau transition, and diffusive behavior

at the critical energy is shown to be related to the critical exponents of a

class of chiral polymers at the $\theta$-point. The exact critical exponents of

the chiral polymer model on the honeycomb lattice are found, verifying that

this model is in the same universality class as a previously solved model of

polymers on the Manhattan lattice. The mapping is obtained by averaging

analytically over the local random potentials in a previously studied lattice

model for the classical plateau transition. This average generates a weight on

chiral polymers associated with the classical localization length exponent $\nu

= 4/3$. We discuss the differences between the classical and quantum

transitions in the context of polymer models and use numerical results on

higher-moment scaling laws at the quantum transition to constrain possible

polymer descriptions. Some properties of the polymer models are verified by

transfer matrix and Monte Carlo studies.

two-dimensional self-interacting lattice polymers. This mapping is exact in the

classical percolation limit of the plateau transition, and diffusive behavior

at the critical energy is shown to be related to the critical exponents of a

class of chiral polymers at the $\theta$-point. The exact critical exponents of

the chiral polymer model on the honeycomb lattice are found, verifying that

this model is in the same universality class as a previously solved model of

polymers on the Manhattan lattice. The mapping is obtained by averaging

analytically over the local random potentials in a previously studied lattice

model for the classical plateau transition. This average generates a weight on

chiral polymers associated with the classical localization length exponent $\nu

= 4/3$. We discuss the differences between the classical and quantum

transitions in the context of polymer models and use numerical results on

higher-moment scaling laws at the quantum transition to constrain possible

polymer descriptions. Some properties of the polymer models are verified by

transfer matrix and Monte Carlo studies.

Type

Identifier

doi:10.1103/PhysRevB.65.035307

Phys. Rev. B 65, 035307 (2002)

Database

Link to record

Show preview

Hide preview

ar
X

iv :c

on d-

m at

/0 10

40 33

v3 [

co nd

-m at.

me s-h

all ]

4 O ct

20 01

Connecting polymers to the quantum Hall plateau transition

Joel E. Moore Bell Labs Lucent Technologies, 600 Mountain Avenue, Murray Hill, NJ 07974

(February 1, 2008)

A mapping is developed between the quantum Hall plateau transition and two-dimensional self- interacting lattice polymers. This mapping is exact in the classical percolation limit of the plateau transition, and diffusive behavior at the critical energy is shown to be related to the critical exponents of a class of chiral polymers at the θ-point. The exact critical exponents of the chiral polymer model on the honeycomb lattice are found, verifying that this model is in the same universality class as a previously solved model of polymers on the Manhattan lattice. The mapping is obtained by averaging analytically over the local random potentials in a previously studied lattice model for the classical plateau transition. This average generates a weight on chiral polymers associated with the classical localization length exponent ν = 4/3. We discuss the differences between the classical and quantum transitions in the context of polymer models and use numerical results on higher-moment scaling laws at the quantum transition to constrain possible polymer descriptions. Some properties of the polymer models are verified by transfer matrix and Monte Carlo studies.

PACS numbers: 73.43.Nq, 72.15.Rn, 61.41.+e

I. INTRODUCTION

The quantum Hall plateau transition is of great inter- est because it links Anderson localization and the quan- tum Hall effect (QHE), two of the fundamental phenom- ena of condensed matter physics. Noninteracting elec- trons moving in two dimensions in a random potential form energy eigenstates which do not extend to infinity but are exponentially localized in finite regions of the plane. The situation changes drastically in a magnetic field: there are then extended states at discrete critical energiesEc separated by the cyclotron energy ~ωc. These extended states are remnants of the Landau bands at zero disorder. At other energies the electron eigenstates are localized, and the localization length near a critical en- ergy scales according to

ξ(E) = ξ0

( Ec

E − Ec

)ν . (1.1)

An understanding of this behavior, which determines the passage from one quantum Hall plateau to another, is essential to the explanation of the integer QHE1. De- spite progress in finding an effective theory for this tran- sition2,3, the scaling law (1.1) and other universal prop- erties of the transition have still not been obtained ana- lytically. The clearest picture for the scaling (1.1) remains the

connection between hulls of percolation clusters and clas- sical electron trajectories in a strong magnetic field and random potential4. Classical percolation was recently shown to describe correctly some aspects of the spin quantum Hall transition5, but the ordinary quantum Hall transition is known from numerical studies6,7 not to lie in the percolation universality class. Considerable effort has been devoted to how quantum effects modify the per- colation picture, and while there is now an understand- ing via numerics of the essential ingredients required to

model the transition7, analytic progress on generalizing percolation has been quite limited8. This paper devel- ops a mapping between the plateau transition and the physics of self-interacting two-dimensional lattice poly- mers. Many statistical properties are known for such polymers because they in turn can be mapped to sim- ple magnetic systems with known critical properties9,10. The final section discusses how quantum effects might be incorporated in a polymer description. However, nu- merically determined scaling laws for higher moments in the quantum case, which result from strong interference of quantum paths, seem to rule out a simple polymer mapping. The properties of a lattice version of the classical

plateau transition are mapped after disorder averaging onto polymers at the θ-point, which is a tricritical point separating the collapsed and extended phases of polymers with an attractive short-ranged self-interaction. The polymer model for the classical transition gives a useful complementary picture to the percolation description, in which some facts, such as diffusion at the critical energy, are more easily recovered. A side benefit is that some new results on polymers come out naturally from the mapping to the plateau transition. The polymer mapping provides an alternate connec-

tion between the classical limit of the plateau transition and (classical) percolation4,11,12, as ring polymers at the θ-point are equivalent to percolation hulls13,14. One re- sult of this paper is that directly mapping the plateau transition to polymers without the intermediate step of percolation gives many more relations. The exponent νθ which governs the typical size of a polymer (for a polymer of N units, 〈R2〉 ∼ N2νθ , where 〈〉 denote averages over the ensemble of polymers), and µ and γ which determine essentially the number ∼ µNNγ−1 of polymers of length N , can all be connected to the plateau transition. The localization length exponent in (1.1) is ν = 43 for

1

classical percolation, while numerical studies6,7 for the quantum case predict ν = 2.35±0.05 for the lowest Lan- dau level (LLL), consistent with experiments15. In this paper ν will be studied via the subdiffusive propagation of electrons in a magnetic field and quenched random po- tential, which is now reviewed. Recently it was shown by Sinova, Meden, and Girvin16 that the localization length exponent ν appears in the energy-integrated correlation function Π(x, t) ≡ 〈〈ρ¯(0, 0)ρ¯(x, t)〉〉, where ρ¯ is the LLL- projected electron density operator and 〈〈〉〉 denote disor- der averaging. (The discussion in this paper is generally restricted to the LLL, although ν is believed to be uni- versal.) The Fourier transform was verified numerically to have the scaling form

ωImΠ(q, ω) = ω 1

2ν f(q2/ω) (1.2)

in the limit q, ω → 0 with q2 ≪ ω. This scaling form can be understood as resulting from a simple form for Π(q, ω) in this limit,

Π(q, ω) ∼ 1 ω − iD(ω)q2 , (1.3)

with the frequency-dependent diffusion constant D(ω) ∝ ω

1

2ν . The result (1.2) depends on the assumption that only

at isolated critical energies Ec are there extended states. It can be understood from the following argument, which is somewhat different from that in16. Electrons at en- ergy E with localization length ξ(E) move diffusively over short times but cross over to localized behavior once t ≥ ξ(E)2/D0. The diffusion constant D0 should have a finite limit as E → Ec since the conductivity σxx is fi- nite at the transition, and can be approximated by this limiting value in the scaling limit. So for a particle at the origin at t = 0 (where it projects onto eigenstates of different energies),

〈x2(t)〉 ≈ ∫ ∞ ξ(E)=

√ D0t

dE ρξ2(E) +

∫ ξ(E)=√D0t Ec

dE ρD0t.

(1.4)

where ρ is the density of states near Ec. This results in subdiffusive behavior:

D˜ ≡ d dt 〈x2(t)〉 = ρEcD0

( ξ0√ D0

) 1 ν

t− 1

2ν , (1.5)

which corresponds to (1.2) with f(x) ∝ x for x≪ 1. Our starting point to obtain the anomalous diffusion

(1.2) is a single electron moving either classically or quantum-mechanically in the x-y plane in a random po- tential V (x) and strong constant magnetic field Bzˆ. The classical coarse-grained equation of motion

Bx˙i = −ǫij∂jV (x) (1.6)

can also be obtained12 by taking a certain limit of the Li- ouvillian formalism. The lowest-Landau-level projected electron density operator in this limit becomes a clas- sical distribution function of particles moving according to (1.6). Of course, the equation of motion (1.6) can also be derived simply from classical physics: a single electron moving in constant electric and magnetic fields with E < B has average velocity EB c along the direction E×B. Since the direction of motion is always perpendic- ular to ∇V , the particle moves along a constant-energy contour of the potential. The picture underlying network models7 of the transition is that electron propagation is nearly classical except near a saddle point of the poten- tial, where quantum tunneling becomes significant. The first part of this paper shows that a discrete-time

lattice version of (1.6), known to have the correct (per- colative) critical scaling for the classical limit, maps after disorder averaging onto a model of two-dimensional in- teracting polymers on the same lattice. Although the lattice is useful to derive the mapping, the critical prop- erties related by the mapping are universal and hence lattice-independent. In the remainder of the introduc- tion, we outline the lattice model of (1.6) and some basic properties of interacting polymers, then summarize the main results. In order to establish the connection between polymers

and motion along level surfaces, we use a lattice model due to Gurarie and Zee12. The particle is taken to have constant velocity along level surfaces: a nonzero mean velocity at criticality was found numerically in11,12 for similar models, and fixing the particle velocity does not alter the critical scaling. Particles move on the edges of the honeycomb lattice of Fig. 1, where each hexagonal face has an associated random potential energy. Except in section II, the energy E of a particle starting at vertex A will be taken to be the average of the three neigh- boring potentials V1, V2, V3, instead of an independent quantity as in12. The energy E is constant along the particle trajectory. The particle’s first step is chosen so that the potential to the left is larger than the particle energy E, which is larger than the potential to the right. In successive steps, there is always a choice between two directions aside from the direction by which the parti- cle entered, and only one of these choices will satisfy the condition that the energy to the left (right) be greater (less) than E. For each realization of the random poten- tials and each starting point, there is a unique locus of the particle after N steps. The connection to the classi- cal localization exponent ν = 43 is that the mean square displacement after N steps is found to show subdiffusive behavior:

〈〈R2(N)〉〉 ∼ N1− 12ν (1.7) in accord with (1.5). The value 1− 12ν ≈ 0.62 was found by numerical simulation12 of (1.6), compared to the pre- dicted value 58 = 0.625. Now the connection between motion along level sur-

faces (1.6) and percolation hulls is quickly reviewed. Con-

2

1V V3

B V2

V4A

FIG. 1. Sample self-avoiding walk from A to B of 20 steps, with 19 neighboring hexagons. The dotted edges are self-contacts: the edge between V3 and V4 is an antiparal- lel self-contact, while that between V1 and V2 is a parallel self-contact. This path is not allowed classically since the walk passes V2 both on the right and on the left.

sider the level-surface motion on the hexagonal lattice with white-noise potentials (i.e., the disorder correlation length λD is less than the lattice spacing). Then all po- tentials on the left of a trajectory are higher than the trajectory energy E, which is higher than all potentials on the right. Now let all faces with energy higher than E be “colored”, while those with energy lower than E remain uncolored. The trajectory is then a hull separat- ing colored faces from uncolored ones, and the properties of such hulls are a standard problem in percolation. Al- though the choice of lattice affects such properties as the density of colored faces at the critical point, critical expo- nents are universal (independent of the lattice). Similarly a lattice model is used here to establish the mapping to a polymer problem, but the universal polymer properties γ and ν discussed below do not depend on this lattice. The usual way to study this type of lattice model11,12

is by summing numerically over paths at fixed particle energy. Here the path will be held fixed for the inte- gration over random potentials: the goal is to assign a weight to each path according to the fraction of the space of random potentials for which that path is the particle trajectory. For convenience, the potentials are assumed to be uniformly distributed on [−1, 1], so the critical en- ergy is E0 = 0. After N steps the particle has either moved along a self-avoiding walk (SAW) of length N , or else has looped and begun retracing previous steps. There is a constraint of “no parallel self-contacts” (the terminology is explained in section II) on allowed SAWs resulting from the restriction that no hexagon can be passed on both the left and right (Fig. 1). The probability P (a, b,N) to reach b from a after N

steps can be written as a sum over an ensemble of paths including both closed self-avoiding polygons (SAPs) as well as open SAWs with no parallel self-contacts. Writing W i for the weight of closed or open curve i, the disorder-

averaged probability to be at b after N steps starting from a is (H is the number of different hexagons visited by the SAW or SAP)

P (a, b,N) ∝ ∑

SAPs i through a and b, l= length of SAP, q= steps froma to b

δN mod l,qW i

+ ∑

SAWs j of lengthN from a to b

no ‖ self−contacts

W j . (1.8)

In the above each SAP should actually be summed twice, once with distance q and once with distance l − q. The δ-function in the SAW part ensures that the particle lo- cation after N steps is b. Sections II and III carry out the disorder average to calculate the weightsW i,j exactly for cases of interest. The weights turn out to have a natural interpretation in terms of self-interacting polymers. In section II we demonstrate that the trajectories of

the lattice model at the critical energy E are related to a chiral polymer model, whose exact critical properties are found. The property of diffusion at the critical energy is shown to follow from the critical exponents γ = 67 and

νθ = 4 7 of the polymer problem. Then in section III we

modify the model so that the energy of a trajectory is a function of initial position rather than an independent quantity, and study the localization exponent ν. An in- equality is derived which connects ν to the exponents of the associated polymer problem. Finally, in section IV we discuss modifications resulting from quantum-mechanical effects, which are most clearly apparent in scaling laws for higher moments of the particle distribution function.

II. TRAJECTORIES AT THE CRITICAL

ENERGY

At the critical energy, the particle motion is diffusive: 〈〈R2(N)〉〉 ∼ N for long timesN . This section shows that the conditional probability P (a, b,N,E0) for the particle moves from a to b after N steps, given that the particle energy is the critical energy E0, is related to critical prop- erties of self-interacting polymers at the θ-point. Then in the following section the same mapping will be shown to give information about trajectories at other energies, and hence about ν. We note in passing that in the percola- tion picture, diffusion at the critical energy is somewhat surprising. A particle on the hull of the infinite cluster moves superdiffusively, while a particle on a finite cluster has only bounded motion: the diffusive motion obtained after averaging over initial position essentially interpo- lates between these two limits. For a path PAB at the critical energy E0, the proba-

bility that PAB is the trajectory in a random potential realization is proportional to 2−HL2−HR = 2−H . Here HL (HR) is the number of hexagons passed on the left (right) by the path, and H = HL +HR: the probability

3

2−H comes about because each hexagon i with potential Vi is as likely to have Vi > E0 as to have Vi < E0. Then the ensemble (1.8) becomes

P (a, b,N,E0) = ∑

SAPs through a and b, l= length of SAP, q= steps from a to b

δN mod l,q 2H

+ ∑

SAWsof lengthN from a to b

no ‖ self−contacts

1

2H . (2.1)

The connection to self-interacting polymers appears because the number of hexagons visited by an SAW is related to the number of self-contacts of the SAW. A self-contact is a point where the SAW is within one edge of intersecting itself. Counting hexagons in lieu of self-contacts gives rise to the famous θ′ model13,14 of a two-dimensional self-interacting polymer. The num- ber of hexagons visited by an SAW of length N is H = N + 1−N2 − 2N3, where N2 and N3 are the numbers of hexagons visited twice and thrice by the SAW. Checking possible paths on the lattice shows thatH = N+1−I−I ′, where I is the number of self-contacts and I ′ the num- ber of a certain type of next-nearest-neighbor contacts. The effects of I ′ are not believed to alter the universality class of the model13,14 and will be ignored. The weight 2−H = 2−N−1+I thus corresponds to the grand-canonical ensemble for polymers at chemical potential µ = − log 2 and with an attractive interaction energy βU = − log 2 for each self-contact. We call a self-contact parallel (an- tiparallel) if, once a direction is defined along the poly- mer, the two sections of polymer in contact have the same (opposite) direction. For a long polymer, almost all self- contacts are antiparallel, as might be expected since par- allel self-contacts are a boundary effect, in the weak sense that a closed polymer has none. There are three phases of the θ′ model for a two-

dimensional self-interacting polymer. At high temper- ature, the statistical properties are those of the non- interacting SAW, and the mean radius of gyration is R ∼ N3/4. At low temperature, the polymer is in a collapsed phase with R ∼ N1/2. There is a tricritical point, called the θ-point, separating these two behaviors, with R ∼ N4/7. The importance of the θ-point for the plateau transition is that the weight 2−H corresponds ex- actly to the θ-point on a honeycomb lattice. The chirality constraint will be shown to change the scaling and give the same universal properties as the solvable Manhattan lattice θ-point. Now the diffusion at the critical energy can be obtained

from (1.8). Considering for the moment only the SAW term in (1.8), the mean particle displacement after N steps is

〈〈R2(N)〉〉 = ∑ b

P (a, b,N)(xb − xa)2

= ∑

SAWsof lengthN

R2SAW 2H

∼ µNNγ−1+2νθ . (2.2)

Here we have introduced the standard polymer exponents γ and νθ, defined through

∑ SAWsof lengthN

1

2H ∼ µNNγ−1

∑ SAWsof lengthN

R2 SAW

2H∑ SAWsof lengthN

1 2H

∼ N2νθ . (2.3)

For ordinary polymers (no chirality constraint) at θ, µ = 1, γ = 87 , and νθ =

4 7 . The effect of the chirality

constraint is clearly to reduce γ, since some polymers are forbidden: in fact we now show that γ = 67 with the chirality constraint (µ and νθ are unchanged), so that 〈〈R2(N)〉〉 ∝ N in (2.2) and motion is diffusive at the critical energy. The critical properties of chiral polymers at θ are ac-

tually related in a very simple way to those of ordinary polymers at θ. As shown at the end of this section, the transfer matrix for L chiral polymers on a cylinder of finite circumference N hexagons has the same leading eigenvalue as the transfer matrix of 2L nonchiral poly- mers on the same cylinder. This means that the critical exponents for the chiral model can be deduced from the known values for the nonchiral model. The “watermelon” exponents xL

10,14 are defined from the correlation functions Gn,L(a, b) of L mutually avoid- ing SAW’s from a to b at criticality: Gn,L(a, b) ∝ |a − b|−2xL(n). The natural generalization for the chiral case is that the L self-avoiding walks have no parallel self- contacts. Then the computation below of the transfer matrix at criticality (µ = 1) on finite strips shows that the chiral exponents x˜L are identical to the nonchiral ex- ponents x2L for twice as many connectors. The nonchiral values xL = (L

2 − 1)/12 known from the Coulomb-gas technique10 then determine all the x˜L. Now we connect the watermelon exponents to physical

properties such as γ and νθ. First, the size exponent νθ = (2 − x2)−1 is unchanged by the chirality constraint because it takes the same value for ring polymers as for linear polymers, and ring polymers are unaffected by the chirality constraint. The exponent γ is given by νθ(d − 2x˜1) = νθ(d − 2x2) = 67 . Note that the ring exponent α is14 also 67 so the two terms of (1.8) scale with the same power of N , as required for consistency. The diffusion result 〈〈R2(N)〉〉 ∝ N might seem almost

coincidental. However, it follows directly from the rela- tionship x˜1 = x2 between the chiral exponent with one leg and the nonchiral exponent with two legs:

γ + 2νθ − 1 = 4− 2x2 2− x2 − 1 = 1. (2.4)

The “mysterious cancellation of exponents”12 which yields diffusion in the percolation picture is relatively

4

FIG. 2. A directed walk on the Manhattan lattice. Note that all allowed self-contacts on this lattice are antiparallel.

simple in the polymer picture, and does not depend on the specific value of x2. The result x˜L = x2L for chiral polymers is exactly

the same as for polymers at θ on the Manhattan directed lattice (Fig. 2), which by construction has no parallel self- contacts17. Hence we learn that the detailed structure of the Manhattan lattice is in some sense irrelevant: it is the short-ranged constraint of no parallel self-contacts which determines the universality class. Another piece of information about polymers follows from the beautiful

result of Cardy18 for the conductivity σxx = √

3e2

4h at the critical energy. This fixes the lattice diffusion constant through the Einstein relation11,12, and therefore predicts a value for the combination of prefactors in (2.2). The chiral polymer model discussed here is just one

point of a two-parameter family of models with antipar- allel self-contacts weighted by some real number w and parallel self-contacts by some possibly different number v. Then w = v gives ordinary two-dimensional self- interacting polymers, while v = 0 gives the chiral poly- mer ensemble. One expects a “coiled” polymer phase for w = 0 and v → ∞, different from the conventional col- lapsed polymer phase. The full phase diagram of these models in the (w, v) plane is a rich subject; the corre- sponding problem defined in terms of parallel and an- tiparallel self-contacts, rather than hexagons, has been investigated numerically on the square lattice19,20. A conformal field theory approach suggests the possibility of continuously varying γ between the chiral polymer and ordinary polymer θ-points21. We remark in passing that the lattice θ′ model defined here in terms of hexagon weights (w, z) has a number of advantages for this prob- lem: the critical point w = 12 is known exactly, and the exact relation discussed below between transfer matri- ces suggests that this model may be solvable by vertex methods. The remainder of this section establishes the connec-

tion x˜L = x2L between the watermelon exponents of chi- ral polymers on the hexagonal lattice and those of nonchi- ral polymers, and can be skipped by nonspecialists. A powerful method to obtain properties of polymer models

is by using conformal invariance to analyze the results of finite-size transfer matrix calculations22,23. Since there are reviews of the technique23,24, some minor details will be omitted. The goal will be to find the transfer matrix for L poly-

mers on a cylinder of circumference h hexagons. The scaling dimensions x˜L can be derived from the finite- size correlation length ξL,h, which is determined by the largest eigenvalue λL,h of the transfer matrix:

ξL,h = − 1 logλL,h

. (2.5)

The x˜L are then estimated for successively larger cylin- ders using

x˜L,h = − (

2√ 3

) h

2π logλL,h. (2.6)

The geometrical factor 2√ 3 comes from the hexagonal lat-

tice dimensions and is 1 for a square lattice. The transfer matrix acts on “configurations” of hor-

izontal edges. A configuration consists of the state of all the horizontal edges, plus information on which oppo- sitely directed edges are paired (originated from the same loop), plus information on which hexagons between hor- izontal edges have been passed on the left or right. The entries in the matrix sum over the different possible states of the vertical edges which can link two configurations. Each entry is weighted by a factor w for each hexagon passed on one side, v for each hexagon passed on both sides, and µ for each edge. For the critical point of chiral polymers, w = 12 , v = 0, and µ = 1; setting v = w =

1 2

gives the ordinary θ-point, while taking w = v = 1 and µ = µc = (2 +

√ 2)−1 gives the critical point of noninter-

acting SAWs. Note that parallel self-contacts can only occur for one

polymer in the cylindrical geometry when the polymer winds around the cylinder. As a result the surface crit- ical exponents, which follow from the transfer matrix on the strip (closed boundary conditions) rather than on the cylinder (periodic boundary conditions), are un- modified from the nonchiral case. Upon conformal map- ping from the cylinder back to the plane, polymers which wrap around the cylinder become polymers which wrap around the origin of the plane, and closed boundary con- ditions correspond to a branch cut which polymers can- not cross. The equivalence of surface exponents to those of the nonchiral θ-point was previously obtained for the Manhattan lattice25. Table I gives the estimated scaling dimensions x˜1 and

x˜2 from cylinders of various sizes. The first few can be done by hand, while the larger matrices are done by com- puter. The leading eigenvalues are exactly the same as for those of twice as many nonchiral polymers. This con- nection is in retrospect not too surprising, since the con- dition of no parallel self-contacts for a polymer from A to B in the chiral case is exactly the condition that another

5

h x˜1,h = x2,h xextrap. x˜2,h = x4,h xextrap.

2 0.254768 3 0.259127 0.23788 4 0.256212 0.24859 1.52861 1.2814 5 0.254221 0.24911 1.39654 1.26398 6 0.253007 1.34308 1.25603 7 0.25224 1.31513 8 1.2984

∞ 1

4 = 0.25 5

4 = 1.25

TABLE I. Results of transfer matrix calculations for one and two chiral polymers (identical to results for two and four nonchiral polymers), on cylinders of circumference h hexagons. The largest eigenvalue of the transfer matrix is related to xL,h through (2.6). The convergence to the pre- dicted values x˜1 =

1

4 and x˜2 =

5

4 is seen to be quite rapid.

The extrapolated values are obtained by using three consec- utive values of xh to fix the constants in xh = c1 + c2h

−c3 , then taking c1 as an estimate of x∞.

polymer can be added from A to B in the nonchiral case. The subleading eigenvalues can differ, however, so there may not be a simple equivalence between states of L chi- ral polymers and 2L nonchiral polymers. The critical properties of the nonchiral model follow from Coulomb gas results for the O(n) model10, so we have

xL = L2 − 1 12

, x˜L = 4L2 − 1

12 . (2.7)

The chiral exponents x˜L are the same as those of the θ-point on the Manhattan lattice26. There are Monte Carlo results for another hexagonal lattice model believed to lie in the Manhattan universality class, the “smart kinetic growth SAW”27, which are consistent with the above values.

III. THE CLASSICAL LOCALIZATION

EXPONENT

When the particle energy E moves away from the crit- ical energy, the trajectories become less extended and the mean distance from the origin after N time steps is reduced. In this section, we take the lattice disorder average in order to express P (a, b,N), the probability that after N steps the particle has moved from a to b, as a weighted sum over linear and ring polymers. Simple properties of the weight function then yield an inequality connecting the localization exponent ν to polymer expo- nents γ and νθ. As in the preceding section, we fix an open or closed

curve on the lattice and ask what fraction of poten- tial realizations make this curve the correct trajectory.

The particle energy E is also varied in order to find the energy-integrated diffusion constant, and hence ν. The weight of a curve is determined by the numbers of hexagons touched by the curve to the left and right: the requirement for a path to be the correct trajectory is that all the hexagons to the right lie above the particle energy. The probability that a path is the particle trajectory is a function of the number of different hexagons visited by the path. The requirement is that all the hexagons to the immediate left have energies larger than the particle energy, while those to the immediate right have energies smaller than the particle energy. The HL hexagons on the left must have higher energies than the HR on the right, which is true for (HL+HRHL )

−1 of potentials. Fur- thermore, the particle energy must lie in the window of width ∼ (HL + HR)−1 between the lowest potential on the left and the highest potential on the right. So the weight of an allowed path PAB is

W (PAB) ∝ 1 (HL +HR)(

HL+HR HL

) ≈ HL!HR!

(HL +HR + 1)! .

(3.1)

The same result is obtained by integrating the probability (1+E2 )

HL(1−E2 ) HR over particle energy E to obtain a beta

function. Now we can again connect the expression (1.8) and the

weight (3.1) to known properties of polymers. For fixed H = HL+HR, the weight is minimized if HL = HR, and

using Stirling’s approximation is then W (H) ≈ √

2pi 2HH1/2

. The probability to get from a to b after N steps thus satisfies

P (a, b,N) ≥ ∑

SAPs througha and b, l= length of SAP, q=steps froma to b

δN mod l,q 2HH1/2

+ ∑

SAWsof lengthN froma to b

no ‖ self−contacts

1

2HH1/2 , (3.2)

up to a possible numerical constant. The typical number of hexagons H scales linearly in N to sufficient accuracy thatH1/2 can be replaced byN1/2 (this is verified numer- ically by Monte Carlo simulations, and if false would re- quire an unexpected multifractality at the θ-point). Sum- ming over final positions b to find the mean squared dis- placement then gives

〈〈R2(N)〉〉 ≥ Nγ−1+2νθ− 12 . (3.3) Then from equation (1.7) we obtain an inequality con- necting the localization exponent ν for the plateau tran- sition to polymer exponents νθ and γ:

1− 1 2ν ≥ γ + 2νθ − 3

2 . (3.4)

For the chiral polymer model, the resulting prediction is ν ≥ 1, which is satisfied by the actual value ν = 43 .

6

The usual nonchiral polymer exponents at θ would pre- dict ν ≥ 73 (cf. section IV), so again it is seen that the chirality constraint is essential. The fact that the lower bound is not reached shows that even in the limit of long paths, the number of hexagons to the left and right of the path cannot be assumed equal in calculating the weight (3.1). At the critical energy (section II), hexagons to the left and right contribute equally and this difference is ir- relevant, but away from the critical energy the difference affects the scaling. The exact value ν = 43 is derived in the percolation

picture4 from the equivalence of closed trajectories at energy E to percolation hulls with p−pc ∝ E−E0, where pc is the critical probability for percolation and E0 is the critical energy. Such percolation hulls28 have average size ξ ∼ (p − pc)− 43 . We remark in passing that the value ν = 43 can be understood in the polymer context from the fact the crossover exponent of the tricritical θ-point is φ = 37 (this value was first obtained using the connection

to percolation14): then N−φ ∼ (p − pc) and ξ ∼ (p − pc)

−νθ/φ = (p − pc)− 43 . The reasons for stressing the inequality (3.4) here rather than the exact result are that the inequality follows immediately from the classical path weight on polymers and can be used to gain information on higher moments (Section IV). The main result of this section is that the exact path

weight induced by averaging over disorder and particle energy can be calculated for the classical lattice model. This weight yields the inequality (3.4) connecting statis- tics of chiral polymers to the critical exponent ν of the classical plateau transition. The focus of the next section will be whether a similar relationship to polymers exists for the quantum plateau transition.

IV. HIGHER MOMENTS AND THE QUANTUM

TRANSITION

The quantum Hall plateau transition shows several qualitative similarities with the semiclassical limit stud- ied in the preceding sections of this paper. Both the quantum transition and the semiclassical limit have power-law delocalization at the critical energy and a fi- nite critical conductivity. However, the quantum tran- sition has proved much more difficult to describe theo- retically, and remains a major open problem. A natural question is whether any generalization of the polymer mapping developed for the classical limit would serve as a useful approach for the quantum case. The goal of this section is to show that the diffusive behavior usu- ally assumed to exist up to the localization length in the quantum case would be inconsistent with almost any such generalization, and present a numerical method and pre- liminary results to verify this assumption. We focus on one generalization in particular (to nonchiral polymers at θ, for reasons described below) for conciseness.

The scaling laws of moments of the particle distribu- tion function

〈R2n(t)〉 ∼ tα(n) (4.1)

demonstrate an essential difference between the classi- cal limit and the conventional picture of the quantum case. In this section 〈〉 indicates averaging over parti- cle energy and random potentials, while 〈〉E indicates averaging over random potentials at fixed particle en- ergy E. If particle motion is essentially diffusive on short length scales in the quantum case, then higher scaling laws beyond n = 1 in (4.1) do not contain additional in- formation. As discussed in the previous section, the mean squared displacement 〈R2n(t)〉 ∼ t1− 12ν contains the lo- calization exponent ν. This formula was obtained from the assumption that 〈R2(t)〉E increases linearly in time at each energy until the localization length is reached (R2 ≈ ξ(E)2), then saturates. If the particle motion is truly diffusive up to the local-

ization length, then 〈R2n(t)〉E ∼ tn until the localization length is reached, and

〈R2n(t)〉 ∼ tn− 12ν , (4.2)

with the localization length exponent ν ≈ 2.35 ± 0.056. So if motion in the quantum case is diffusive up to the localization length, there are no nontrivial exponents to be found in higher moments of the particle displacement. Higher moments in the classical case show nontrivial

scaling, and consequently highly extended trajectories are much more common in the classical case than the quantum case, even though the localization length di- verges more rapidly near Ec for the quantum case. For walks at the critical energy, it follows from the results of Section II that

〈R2n(t)〉Ec ∼ t2nνθ+γ−1 = t(8n−1)/7. (4.3)

Hence, although the mean square displacement does in- crease linearly with time, higher moments have nontrivial power laws because the particle trajectories are not ran- dom walks but instead have “memory,” as required for the absence of self-intersections. The higher moments are more extended than they would be for simple diffu- sive motion (random walks). A similar result holds for the classical case even when

the average is extended to include the particle energy. The inequality (3.4) derived in Section III between poly- mer exponents at the chiral θ-point and the scaling of the second moment α(1) = 1− 12ν from (4.1):

1− 1 2ν ≥ 2νθ + γ − 3

2 =

1

2 . (4.4)

For the classical transition with ν = 43 , the left side is 5 8 and the inequality is satisfied. Similarly for higher

moments 〈R2n(t)〉 ∼ tαc(n)

7

αc(n) ≥ 2nνθ + γ − 3 2 =

1

2 +

8(n− 1) 7

. (4.5)

Hence for sufficiently large n the classical scaling ex- ponents αc(n) are larger than the quantum exponents α(n) = n − 12ν , if the motion in the quantum case is indeed diffusive. In the remainder of this section, we consider the ques-

tion of how quantum interference keeps the quantum case from being related to a polymer ensemble in the same way as the classical case. One polymer ensemble in particular is attractive for the quantum case because the ensemble is similar to the classical one and the value ν = 73 appears in this ensemble, but this connection predicts nondiffusive motion up to the localization length. Monte Carlo nu- merics are used to test the assumption of diffusive motion in the quantum case; if verified this assumption would rule out a simple connection to polymers. The approach of the previous section was to generate a

positive weight on electron paths by averaging over disor- der with the electron path held constant. In lattice mod- els for the quantum transition, it should be possible to at- tribute a positive weight to each path on the lattice, and then these weights may be connected to some polymer problem, presumably different from the chiral ensemble discussed above for the classical limit. The first state- ment (that there is an assignment of weights) is some- what trivial from a mathematical point of view: there are many paths on the lattice between any pair of points, and hence given any positive probabilities to reach different points on the lattice after N steps, there is some assign- ment of positive path weights which results in the given probabilities. The difficult question is whether there is an assignment of weights which is physically meaningful and related to some local two-dimensional theory, as in the classical case. The next paragraphs define a real (not necessarily positive) path weight; this weight is positive in the absence of interference, and the leading interfer- ence contribution to this path weight from “cooperons” vanishes, though higher contributions do not. For fixed disorder, different paths W iAB from A to B

contribute to the amplitude, and the probability PAB to get fromA to B includes both diagonal terms |W iAB|2 and cross terms. Here and in the following we assume a dis- cretized model for the quantum case, similar to the lattice model introduced previously for the classical limit. We start by considering two paths which do not cross: then the disorder average generates a random phase which cancels the cross terms, leaving only the (positive) di- agonal terms. The remaining question is what occurs for intersecting paths; this is known to be the case of interest for the quantum effects causing weak localization. If the cross terms did vanish, then in the discretized

model where λD is effectively zero, PAB would be a sum over (not necessarily self-avoiding) paths with some pos- itive weight, the “quantum path weight” (QPW). A real weight can be defined even if cross terms are present by adding to the direct contribution of each path half of all

its disorder-averaged cross terms with other paths. The QPW picture can break down if, for energies near the critical energy, the cross terms become large enough to drive the weight negative for long paths. However, the leading interference corrections vanish upon disorder av- eraging, and a finite strength of interference is required to drive the path weight negative, so it is seemingly pos- sible that a positive QPW exists for the paths near the critical energy which determine ν. This motivates the conjecture, tested in the remainder of this section, that the universal large-length-scale properties of the plateau transition may be related to those of some classical gener- alized polymer model (i.e., a sum over paths with positive weights), in similar fashion to the relationships found in sections II and III between the classical percolation limit and the chiral polymer model. The remaining step is to determine whether any uni-

versality class of classical polymers can reproduce the weights which follow from disorder averaging in the quan- tum case. It seems worthwhile to identify possibilities, since exact results have been obtained for many two- dimensional polymer models by Coulomb gas and CFT techniques. The QPW should give nearly the classical weight (3.1) to paths which are classically allowed or in- clude a small number of quantum tunneling events, but should not allow of order N tunneling events for an N - step path since then the motion is simply diffusive even away from the critical energy. A speculative possibil- ity for the quantum transition is the nonchiral polymer ensemble at θ. The expectation that quantum mechan- ics should allow some unfavorable steps (but fewer than ∼ N) matches the fact that a typical polymer in the nonchiral ensemble has some parallel self-contacts, but fewer than of order N . A simple argument that the num- ber of parallel self-contacts is subextensive (< N) is that a ring polymer has no parallel self-contacts, so that par- allel self-contacts are in some sense a boundary property. A surprise is that the value ν = 73 (which has attracted

attention as the simplest rational consistent with numer- ics) appears from exponents of the nonchiral ensemble. The inequality (3.4) connecting the localization exponent ν to nonchiral polymer exponents predicts, since γ = 87 ,

1− 1 2ν ≥ γ + 2νθ − 3

2 =

11

14 , (4.6)

or ν ≥ 73 . Hence the value 73 appears in the critical prop- erties of a polymer ensemble closely related to the poly- mer ensemble describing the classical plateau transition. The the lower bound is realized if paths have asymptot- ically the same number of hexagons to the left as to the right (HR ∼ HL), which should be a better approxima- tion for the less convoluted paths in the quantum case. As seen below, however, this inequality predicts scaling laws for higher moments which appear to be ruled out numerically in the quantum case. Hence quantum inter- ference seems to be relevant at the quantum transition even beyond the level of changing path weights.

8

We note in passing that the desired property of dif- fusion at the critical energy does not have any simple interpretation as a statement about the polymer ensem- ble, since it is only after integrating over particle energy that the polymer weights may appear. This situation is familiar from the Liouvillian approach to the transi- tion16,31, where the transition is mapped onto a different problem which contains the exponent ν but not the crit- ical conductivity. Note that the previous appearance of the value 73 in a semiclassical average over a single per-

colation trajectory29 does not clearly relate to a critical point, when the electron is delocalized over multiple tra- jectories. Now we discuss how numerics can test whether the appearance of this value in this polymer ensemble is just a numerical coincidence. There are numerically testable consequences which

can be used to check whether nonchiral polymers at θ are indeed related to the quantum case. The polymer problem predicts various moments of the disorder- and energy-averaged displacement: it was shown in the pre- vious paragraph that the mean squared displacement 〈R2(N)〉 ∼ N11/14, so that ν = 73 . Similar predic- tions follow for higher moments, such as 〈R4(N)〉 ∼ N4ν+γ−3/2 = N27/14, or α(2) = 2714 . This can be com- pared to the null hypothesis of diffusion up to the local- ization length, which predicts 〈R4(N)〉 ∼ N1.78. We have performed Monte Carlo simulations with up

to 1800 states in the lowest Landau level to track the evo- lution of a localized wave packet in a disordered potential (the method is similar to that of30). The error bars are larger for the fourth moment than for the second because finite-size effects are more pronounced on the extended paths which dominate the fourth moment, but it appears that α(2) = 1.8±0.1, which if correct is sufficient to rule out α(2) = 27/14 ≈ 1.92. With larger system sizes, it should be straightforward to confirm the assumption of diffusive motion up to the localization length. Then a polymer description would have to have ν = 1/2 (either dense polymers or random walks), but no appropriate ensemble is obvious. It seems more likely that strong quantum interference prevents a physically meaningful assignment of path weights in the quantum case. To summarize, this section discussed differences be-

tween the classical and quantum transitions which be- come apparent in higher moments of the particle dis- tribution function. Numerics seem to support the pic- ture of diffusion up to the localization length and rule out the simplest polymer model for the quantum case. In closing, we mention briefly connections between the polymer models discussed in this paper and conformal field theory (CFT) approaches to the transition. The low-temperature O(n) phase also appears in a large-N expansion of the disorder-averaged Liouvillian theory, which is similar to a partially supersymmetric complex O(N) model, with N → 1 the physical limit31. Recent work32,33 on the critical point of Dirac fermions in a non- abelian random vector potential found a c = −2 dense

polymer problem hidden in the critical theory for sev- eral different types of disorder. Since upon adding addi- tional disorder (random mass and chemical potential) the abelian version of this critical point flows to the plateau transition fixed point34, the appearance of polymer sub- algebras may be generic to this class of random critical points. The author wishes to thank H. Saleur and S. Girvin

for helpful suggestions.

1 R. B. Laughlin, Phys. Rev. B 23, 5632 (1981). 2 A. M. M. Pruisken, Nucl. Phys. B 235, 227 (1984). 3 M. Zirnbauer, J. Math. Phys. 38, 2007 (1997). 4 S. Trugman, Phys. Rev. B 27, 7539 (1983). 5 I. A. Gruzberg, A. W. W. Ludwig, and N. Read, Phys. Rev. Lett. 82, 4524 (1999).

6 B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 7 J. T. Chalker and P. D. Coddington, J. Phys. C 21 2665 (1988).

8 S. Das Sarma, in Perspectives in Quantum Hall Effects, ed. S. Das Sarma and A. Pinczuk (Wiley, New York, 1997).

9 P. G. de Gennes, Phys. Lett. A38, 339 (1972). 10 B. Nienhuis, Phys. Rev. Lett. 49, 1063 (1982); B. Nienhuis, in Phase Transitions and Critical Phenomena, ed. C. Domb and J. L. Lebowitz (Academic, London, 1986), vol. 11.

11 F. Evers, Phys. Rev. E 55, 2321 (1997). 12 V. Gurarie and A. Zee, Int. J. Mod. Phys. B 15, 1225 (2001).

13 A. Coniglio et al., Phys. Rev. B 35, 3617 (1987). 14 B. Duplantier and H. Saleur, Phys. Rev. Lett. 59, 539 (1987).

15 H. P. Wei, D. C. Tsui, M. A. Paalanen and A. M. M. Pruisken, Phys. Rev. Lett. 61, 1294 (1988); H. P. Wei, S. Y. Lin, D. C. Tsui, and A. M. M. Pruisken, Phys. Rev. B 45, 3926 (1992).

16 J. Sinova, V. Meden, and S. M. Girvin, Phys. Rev. B 62, 2008 (2000).

17 R. M. Bradley, Phys. Rev. A 41, 914 (1990). 18 J. Cardy, Phys. Rev. Lett. 84, 3507 (2000). 19 A. Trovato and F. Seno, Phys. Rev. E 56, 131 (1997). 20 T. Prellberg and B. Drossel, Phys. Rev. E 57, 2045 (1998). 21 J. L. Cardy, Nucl. Phys. B 419, 411 (1994). 22 B. Derrida, J. Phys. A 14, L5 (1981). 23 H. Saleur, J. Phys. A 20, 455 (1987). 24 C. Vanderzande, Lattice Models of Polymers (Cambridge, 1998).

25 S. L. M. de Queiroz and J. M. Yeomans, J. Phys. A 24, L933 (1991).

26 T. Prellberg and A. Owczarek, J. Phys. A 27, 1811 (1994). 27 D. Bennett-Wood, A. L. Owczarek, and T. Prellberg, Phys- ica A 206, 283 (1994).

28 M. den Nijs, J. Phys. A 12, 1857 (1979). 29 G. V. Mil’nikov and I. M. Sokolov, JETP Lett. 48, 536 (1988).

9

30 S. Boldyrev and V. Gurarie, cond-mat/0009203 (2000). 31 J. E. Moore, A. Zee, and J. Sinova, Phys. Rev. Lett. 87, 8704-8707 (2001).

32 A. LeClair, cond-mat/0011413 (2000). 33 M. J. Bhaseen, J.-S. Caux, I. I. Kogan, and A. M. Tsvelik,

cond-mat/00012240 (2000). 34 A. W. W. Ludwig, M. P. A. Fisher, R. Shankar, and G. Grinstein, Phys. Rev. B 50 7526 (1994).

10

iv :c

on d-

m at

/0 10

40 33

v3 [

co nd

-m at.

me s-h

all ]

4 O ct

20 01

Connecting polymers to the quantum Hall plateau transition

Joel E. Moore Bell Labs Lucent Technologies, 600 Mountain Avenue, Murray Hill, NJ 07974

(February 1, 2008)

A mapping is developed between the quantum Hall plateau transition and two-dimensional self- interacting lattice polymers. This mapping is exact in the classical percolation limit of the plateau transition, and diffusive behavior at the critical energy is shown to be related to the critical exponents of a class of chiral polymers at the θ-point. The exact critical exponents of the chiral polymer model on the honeycomb lattice are found, verifying that this model is in the same universality class as a previously solved model of polymers on the Manhattan lattice. The mapping is obtained by averaging analytically over the local random potentials in a previously studied lattice model for the classical plateau transition. This average generates a weight on chiral polymers associated with the classical localization length exponent ν = 4/3. We discuss the differences between the classical and quantum transitions in the context of polymer models and use numerical results on higher-moment scaling laws at the quantum transition to constrain possible polymer descriptions. Some properties of the polymer models are verified by transfer matrix and Monte Carlo studies.

PACS numbers: 73.43.Nq, 72.15.Rn, 61.41.+e

I. INTRODUCTION

The quantum Hall plateau transition is of great inter- est because it links Anderson localization and the quan- tum Hall effect (QHE), two of the fundamental phenom- ena of condensed matter physics. Noninteracting elec- trons moving in two dimensions in a random potential form energy eigenstates which do not extend to infinity but are exponentially localized in finite regions of the plane. The situation changes drastically in a magnetic field: there are then extended states at discrete critical energiesEc separated by the cyclotron energy ~ωc. These extended states are remnants of the Landau bands at zero disorder. At other energies the electron eigenstates are localized, and the localization length near a critical en- ergy scales according to

ξ(E) = ξ0

( Ec

E − Ec

)ν . (1.1)

An understanding of this behavior, which determines the passage from one quantum Hall plateau to another, is essential to the explanation of the integer QHE1. De- spite progress in finding an effective theory for this tran- sition2,3, the scaling law (1.1) and other universal prop- erties of the transition have still not been obtained ana- lytically. The clearest picture for the scaling (1.1) remains the

connection between hulls of percolation clusters and clas- sical electron trajectories in a strong magnetic field and random potential4. Classical percolation was recently shown to describe correctly some aspects of the spin quantum Hall transition5, but the ordinary quantum Hall transition is known from numerical studies6,7 not to lie in the percolation universality class. Considerable effort has been devoted to how quantum effects modify the per- colation picture, and while there is now an understand- ing via numerics of the essential ingredients required to

model the transition7, analytic progress on generalizing percolation has been quite limited8. This paper devel- ops a mapping between the plateau transition and the physics of self-interacting two-dimensional lattice poly- mers. Many statistical properties are known for such polymers because they in turn can be mapped to sim- ple magnetic systems with known critical properties9,10. The final section discusses how quantum effects might be incorporated in a polymer description. However, nu- merically determined scaling laws for higher moments in the quantum case, which result from strong interference of quantum paths, seem to rule out a simple polymer mapping. The properties of a lattice version of the classical

plateau transition are mapped after disorder averaging onto polymers at the θ-point, which is a tricritical point separating the collapsed and extended phases of polymers with an attractive short-ranged self-interaction. The polymer model for the classical transition gives a useful complementary picture to the percolation description, in which some facts, such as diffusion at the critical energy, are more easily recovered. A side benefit is that some new results on polymers come out naturally from the mapping to the plateau transition. The polymer mapping provides an alternate connec-

tion between the classical limit of the plateau transition and (classical) percolation4,11,12, as ring polymers at the θ-point are equivalent to percolation hulls13,14. One re- sult of this paper is that directly mapping the plateau transition to polymers without the intermediate step of percolation gives many more relations. The exponent νθ which governs the typical size of a polymer (for a polymer of N units, 〈R2〉 ∼ N2νθ , where 〈〉 denote averages over the ensemble of polymers), and µ and γ which determine essentially the number ∼ µNNγ−1 of polymers of length N , can all be connected to the plateau transition. The localization length exponent in (1.1) is ν = 43 for

1

classical percolation, while numerical studies6,7 for the quantum case predict ν = 2.35±0.05 for the lowest Lan- dau level (LLL), consistent with experiments15. In this paper ν will be studied via the subdiffusive propagation of electrons in a magnetic field and quenched random po- tential, which is now reviewed. Recently it was shown by Sinova, Meden, and Girvin16 that the localization length exponent ν appears in the energy-integrated correlation function Π(x, t) ≡ 〈〈ρ¯(0, 0)ρ¯(x, t)〉〉, where ρ¯ is the LLL- projected electron density operator and 〈〈〉〉 denote disor- der averaging. (The discussion in this paper is generally restricted to the LLL, although ν is believed to be uni- versal.) The Fourier transform was verified numerically to have the scaling form

ωImΠ(q, ω) = ω 1

2ν f(q2/ω) (1.2)

in the limit q, ω → 0 with q2 ≪ ω. This scaling form can be understood as resulting from a simple form for Π(q, ω) in this limit,

Π(q, ω) ∼ 1 ω − iD(ω)q2 , (1.3)

with the frequency-dependent diffusion constant D(ω) ∝ ω

1

2ν . The result (1.2) depends on the assumption that only

at isolated critical energies Ec are there extended states. It can be understood from the following argument, which is somewhat different from that in16. Electrons at en- ergy E with localization length ξ(E) move diffusively over short times but cross over to localized behavior once t ≥ ξ(E)2/D0. The diffusion constant D0 should have a finite limit as E → Ec since the conductivity σxx is fi- nite at the transition, and can be approximated by this limiting value in the scaling limit. So for a particle at the origin at t = 0 (where it projects onto eigenstates of different energies),

〈x2(t)〉 ≈ ∫ ∞ ξ(E)=

√ D0t

dE ρξ2(E) +

∫ ξ(E)=√D0t Ec

dE ρD0t.

(1.4)

where ρ is the density of states near Ec. This results in subdiffusive behavior:

D˜ ≡ d dt 〈x2(t)〉 = ρEcD0

( ξ0√ D0

) 1 ν

t− 1

2ν , (1.5)

which corresponds to (1.2) with f(x) ∝ x for x≪ 1. Our starting point to obtain the anomalous diffusion

(1.2) is a single electron moving either classically or quantum-mechanically in the x-y plane in a random po- tential V (x) and strong constant magnetic field Bzˆ. The classical coarse-grained equation of motion

Bx˙i = −ǫij∂jV (x) (1.6)

can also be obtained12 by taking a certain limit of the Li- ouvillian formalism. The lowest-Landau-level projected electron density operator in this limit becomes a clas- sical distribution function of particles moving according to (1.6). Of course, the equation of motion (1.6) can also be derived simply from classical physics: a single electron moving in constant electric and magnetic fields with E < B has average velocity EB c along the direction E×B. Since the direction of motion is always perpendic- ular to ∇V , the particle moves along a constant-energy contour of the potential. The picture underlying network models7 of the transition is that electron propagation is nearly classical except near a saddle point of the poten- tial, where quantum tunneling becomes significant. The first part of this paper shows that a discrete-time

lattice version of (1.6), known to have the correct (per- colative) critical scaling for the classical limit, maps after disorder averaging onto a model of two-dimensional in- teracting polymers on the same lattice. Although the lattice is useful to derive the mapping, the critical prop- erties related by the mapping are universal and hence lattice-independent. In the remainder of the introduc- tion, we outline the lattice model of (1.6) and some basic properties of interacting polymers, then summarize the main results. In order to establish the connection between polymers

and motion along level surfaces, we use a lattice model due to Gurarie and Zee12. The particle is taken to have constant velocity along level surfaces: a nonzero mean velocity at criticality was found numerically in11,12 for similar models, and fixing the particle velocity does not alter the critical scaling. Particles move on the edges of the honeycomb lattice of Fig. 1, where each hexagonal face has an associated random potential energy. Except in section II, the energy E of a particle starting at vertex A will be taken to be the average of the three neigh- boring potentials V1, V2, V3, instead of an independent quantity as in12. The energy E is constant along the particle trajectory. The particle’s first step is chosen so that the potential to the left is larger than the particle energy E, which is larger than the potential to the right. In successive steps, there is always a choice between two directions aside from the direction by which the parti- cle entered, and only one of these choices will satisfy the condition that the energy to the left (right) be greater (less) than E. For each realization of the random poten- tials and each starting point, there is a unique locus of the particle after N steps. The connection to the classi- cal localization exponent ν = 43 is that the mean square displacement after N steps is found to show subdiffusive behavior:

〈〈R2(N)〉〉 ∼ N1− 12ν (1.7) in accord with (1.5). The value 1− 12ν ≈ 0.62 was found by numerical simulation12 of (1.6), compared to the pre- dicted value 58 = 0.625. Now the connection between motion along level sur-

faces (1.6) and percolation hulls is quickly reviewed. Con-

2

1V V3

B V2

V4A

FIG. 1. Sample self-avoiding walk from A to B of 20 steps, with 19 neighboring hexagons. The dotted edges are self-contacts: the edge between V3 and V4 is an antiparal- lel self-contact, while that between V1 and V2 is a parallel self-contact. This path is not allowed classically since the walk passes V2 both on the right and on the left.

sider the level-surface motion on the hexagonal lattice with white-noise potentials (i.e., the disorder correlation length λD is less than the lattice spacing). Then all po- tentials on the left of a trajectory are higher than the trajectory energy E, which is higher than all potentials on the right. Now let all faces with energy higher than E be “colored”, while those with energy lower than E remain uncolored. The trajectory is then a hull separat- ing colored faces from uncolored ones, and the properties of such hulls are a standard problem in percolation. Al- though the choice of lattice affects such properties as the density of colored faces at the critical point, critical expo- nents are universal (independent of the lattice). Similarly a lattice model is used here to establish the mapping to a polymer problem, but the universal polymer properties γ and ν discussed below do not depend on this lattice. The usual way to study this type of lattice model11,12

is by summing numerically over paths at fixed particle energy. Here the path will be held fixed for the inte- gration over random potentials: the goal is to assign a weight to each path according to the fraction of the space of random potentials for which that path is the particle trajectory. For convenience, the potentials are assumed to be uniformly distributed on [−1, 1], so the critical en- ergy is E0 = 0. After N steps the particle has either moved along a self-avoiding walk (SAW) of length N , or else has looped and begun retracing previous steps. There is a constraint of “no parallel self-contacts” (the terminology is explained in section II) on allowed SAWs resulting from the restriction that no hexagon can be passed on both the left and right (Fig. 1). The probability P (a, b,N) to reach b from a after N

steps can be written as a sum over an ensemble of paths including both closed self-avoiding polygons (SAPs) as well as open SAWs with no parallel self-contacts. Writing W i for the weight of closed or open curve i, the disorder-

averaged probability to be at b after N steps starting from a is (H is the number of different hexagons visited by the SAW or SAP)

P (a, b,N) ∝ ∑

SAPs i through a and b, l= length of SAP, q= steps froma to b

δN mod l,qW i

+ ∑

SAWs j of lengthN from a to b

no ‖ self−contacts

W j . (1.8)

In the above each SAP should actually be summed twice, once with distance q and once with distance l − q. The δ-function in the SAW part ensures that the particle lo- cation after N steps is b. Sections II and III carry out the disorder average to calculate the weightsW i,j exactly for cases of interest. The weights turn out to have a natural interpretation in terms of self-interacting polymers. In section II we demonstrate that the trajectories of

the lattice model at the critical energy E are related to a chiral polymer model, whose exact critical properties are found. The property of diffusion at the critical energy is shown to follow from the critical exponents γ = 67 and

νθ = 4 7 of the polymer problem. Then in section III we

modify the model so that the energy of a trajectory is a function of initial position rather than an independent quantity, and study the localization exponent ν. An in- equality is derived which connects ν to the exponents of the associated polymer problem. Finally, in section IV we discuss modifications resulting from quantum-mechanical effects, which are most clearly apparent in scaling laws for higher moments of the particle distribution function.

II. TRAJECTORIES AT THE CRITICAL

ENERGY

At the critical energy, the particle motion is diffusive: 〈〈R2(N)〉〉 ∼ N for long timesN . This section shows that the conditional probability P (a, b,N,E0) for the particle moves from a to b after N steps, given that the particle energy is the critical energy E0, is related to critical prop- erties of self-interacting polymers at the θ-point. Then in the following section the same mapping will be shown to give information about trajectories at other energies, and hence about ν. We note in passing that in the percola- tion picture, diffusion at the critical energy is somewhat surprising. A particle on the hull of the infinite cluster moves superdiffusively, while a particle on a finite cluster has only bounded motion: the diffusive motion obtained after averaging over initial position essentially interpo- lates between these two limits. For a path PAB at the critical energy E0, the proba-

bility that PAB is the trajectory in a random potential realization is proportional to 2−HL2−HR = 2−H . Here HL (HR) is the number of hexagons passed on the left (right) by the path, and H = HL +HR: the probability

3

2−H comes about because each hexagon i with potential Vi is as likely to have Vi > E0 as to have Vi < E0. Then the ensemble (1.8) becomes

P (a, b,N,E0) = ∑

SAPs through a and b, l= length of SAP, q= steps from a to b

δN mod l,q 2H

+ ∑

SAWsof lengthN from a to b

no ‖ self−contacts

1

2H . (2.1)

The connection to self-interacting polymers appears because the number of hexagons visited by an SAW is related to the number of self-contacts of the SAW. A self-contact is a point where the SAW is within one edge of intersecting itself. Counting hexagons in lieu of self-contacts gives rise to the famous θ′ model13,14 of a two-dimensional self-interacting polymer. The num- ber of hexagons visited by an SAW of length N is H = N + 1−N2 − 2N3, where N2 and N3 are the numbers of hexagons visited twice and thrice by the SAW. Checking possible paths on the lattice shows thatH = N+1−I−I ′, where I is the number of self-contacts and I ′ the num- ber of a certain type of next-nearest-neighbor contacts. The effects of I ′ are not believed to alter the universality class of the model13,14 and will be ignored. The weight 2−H = 2−N−1+I thus corresponds to the grand-canonical ensemble for polymers at chemical potential µ = − log 2 and with an attractive interaction energy βU = − log 2 for each self-contact. We call a self-contact parallel (an- tiparallel) if, once a direction is defined along the poly- mer, the two sections of polymer in contact have the same (opposite) direction. For a long polymer, almost all self- contacts are antiparallel, as might be expected since par- allel self-contacts are a boundary effect, in the weak sense that a closed polymer has none. There are three phases of the θ′ model for a two-

dimensional self-interacting polymer. At high temper- ature, the statistical properties are those of the non- interacting SAW, and the mean radius of gyration is R ∼ N3/4. At low temperature, the polymer is in a collapsed phase with R ∼ N1/2. There is a tricritical point, called the θ-point, separating these two behaviors, with R ∼ N4/7. The importance of the θ-point for the plateau transition is that the weight 2−H corresponds ex- actly to the θ-point on a honeycomb lattice. The chirality constraint will be shown to change the scaling and give the same universal properties as the solvable Manhattan lattice θ-point. Now the diffusion at the critical energy can be obtained

from (1.8). Considering for the moment only the SAW term in (1.8), the mean particle displacement after N steps is

〈〈R2(N)〉〉 = ∑ b

P (a, b,N)(xb − xa)2

= ∑

SAWsof lengthN

R2SAW 2H

∼ µNNγ−1+2νθ . (2.2)

Here we have introduced the standard polymer exponents γ and νθ, defined through

∑ SAWsof lengthN

1

2H ∼ µNNγ−1

∑ SAWsof lengthN

R2 SAW

2H∑ SAWsof lengthN

1 2H

∼ N2νθ . (2.3)

For ordinary polymers (no chirality constraint) at θ, µ = 1, γ = 87 , and νθ =

4 7 . The effect of the chirality

constraint is clearly to reduce γ, since some polymers are forbidden: in fact we now show that γ = 67 with the chirality constraint (µ and νθ are unchanged), so that 〈〈R2(N)〉〉 ∝ N in (2.2) and motion is diffusive at the critical energy. The critical properties of chiral polymers at θ are ac-

tually related in a very simple way to those of ordinary polymers at θ. As shown at the end of this section, the transfer matrix for L chiral polymers on a cylinder of finite circumference N hexagons has the same leading eigenvalue as the transfer matrix of 2L nonchiral poly- mers on the same cylinder. This means that the critical exponents for the chiral model can be deduced from the known values for the nonchiral model. The “watermelon” exponents xL

10,14 are defined from the correlation functions Gn,L(a, b) of L mutually avoid- ing SAW’s from a to b at criticality: Gn,L(a, b) ∝ |a − b|−2xL(n). The natural generalization for the chiral case is that the L self-avoiding walks have no parallel self- contacts. Then the computation below of the transfer matrix at criticality (µ = 1) on finite strips shows that the chiral exponents x˜L are identical to the nonchiral ex- ponents x2L for twice as many connectors. The nonchiral values xL = (L

2 − 1)/12 known from the Coulomb-gas technique10 then determine all the x˜L. Now we connect the watermelon exponents to physical

properties such as γ and νθ. First, the size exponent νθ = (2 − x2)−1 is unchanged by the chirality constraint because it takes the same value for ring polymers as for linear polymers, and ring polymers are unaffected by the chirality constraint. The exponent γ is given by νθ(d − 2x˜1) = νθ(d − 2x2) = 67 . Note that the ring exponent α is14 also 67 so the two terms of (1.8) scale with the same power of N , as required for consistency. The diffusion result 〈〈R2(N)〉〉 ∝ N might seem almost

coincidental. However, it follows directly from the rela- tionship x˜1 = x2 between the chiral exponent with one leg and the nonchiral exponent with two legs:

γ + 2νθ − 1 = 4− 2x2 2− x2 − 1 = 1. (2.4)

The “mysterious cancellation of exponents”12 which yields diffusion in the percolation picture is relatively

4

FIG. 2. A directed walk on the Manhattan lattice. Note that all allowed self-contacts on this lattice are antiparallel.

simple in the polymer picture, and does not depend on the specific value of x2. The result x˜L = x2L for chiral polymers is exactly

the same as for polymers at θ on the Manhattan directed lattice (Fig. 2), which by construction has no parallel self- contacts17. Hence we learn that the detailed structure of the Manhattan lattice is in some sense irrelevant: it is the short-ranged constraint of no parallel self-contacts which determines the universality class. Another piece of information about polymers follows from the beautiful

result of Cardy18 for the conductivity σxx = √

3e2

4h at the critical energy. This fixes the lattice diffusion constant through the Einstein relation11,12, and therefore predicts a value for the combination of prefactors in (2.2). The chiral polymer model discussed here is just one

point of a two-parameter family of models with antipar- allel self-contacts weighted by some real number w and parallel self-contacts by some possibly different number v. Then w = v gives ordinary two-dimensional self- interacting polymers, while v = 0 gives the chiral poly- mer ensemble. One expects a “coiled” polymer phase for w = 0 and v → ∞, different from the conventional col- lapsed polymer phase. The full phase diagram of these models in the (w, v) plane is a rich subject; the corre- sponding problem defined in terms of parallel and an- tiparallel self-contacts, rather than hexagons, has been investigated numerically on the square lattice19,20. A conformal field theory approach suggests the possibility of continuously varying γ between the chiral polymer and ordinary polymer θ-points21. We remark in passing that the lattice θ′ model defined here in terms of hexagon weights (w, z) has a number of advantages for this prob- lem: the critical point w = 12 is known exactly, and the exact relation discussed below between transfer matri- ces suggests that this model may be solvable by vertex methods. The remainder of this section establishes the connec-

tion x˜L = x2L between the watermelon exponents of chi- ral polymers on the hexagonal lattice and those of nonchi- ral polymers, and can be skipped by nonspecialists. A powerful method to obtain properties of polymer models

is by using conformal invariance to analyze the results of finite-size transfer matrix calculations22,23. Since there are reviews of the technique23,24, some minor details will be omitted. The goal will be to find the transfer matrix for L poly-

mers on a cylinder of circumference h hexagons. The scaling dimensions x˜L can be derived from the finite- size correlation length ξL,h, which is determined by the largest eigenvalue λL,h of the transfer matrix:

ξL,h = − 1 logλL,h

. (2.5)

The x˜L are then estimated for successively larger cylin- ders using

x˜L,h = − (

2√ 3

) h

2π logλL,h. (2.6)

The geometrical factor 2√ 3 comes from the hexagonal lat-

tice dimensions and is 1 for a square lattice. The transfer matrix acts on “configurations” of hor-

izontal edges. A configuration consists of the state of all the horizontal edges, plus information on which oppo- sitely directed edges are paired (originated from the same loop), plus information on which hexagons between hor- izontal edges have been passed on the left or right. The entries in the matrix sum over the different possible states of the vertical edges which can link two configurations. Each entry is weighted by a factor w for each hexagon passed on one side, v for each hexagon passed on both sides, and µ for each edge. For the critical point of chiral polymers, w = 12 , v = 0, and µ = 1; setting v = w =

1 2

gives the ordinary θ-point, while taking w = v = 1 and µ = µc = (2 +

√ 2)−1 gives the critical point of noninter-

acting SAWs. Note that parallel self-contacts can only occur for one

polymer in the cylindrical geometry when the polymer winds around the cylinder. As a result the surface crit- ical exponents, which follow from the transfer matrix on the strip (closed boundary conditions) rather than on the cylinder (periodic boundary conditions), are un- modified from the nonchiral case. Upon conformal map- ping from the cylinder back to the plane, polymers which wrap around the cylinder become polymers which wrap around the origin of the plane, and closed boundary con- ditions correspond to a branch cut which polymers can- not cross. The equivalence of surface exponents to those of the nonchiral θ-point was previously obtained for the Manhattan lattice25. Table I gives the estimated scaling dimensions x˜1 and

x˜2 from cylinders of various sizes. The first few can be done by hand, while the larger matrices are done by com- puter. The leading eigenvalues are exactly the same as for those of twice as many nonchiral polymers. This con- nection is in retrospect not too surprising, since the con- dition of no parallel self-contacts for a polymer from A to B in the chiral case is exactly the condition that another

5

h x˜1,h = x2,h xextrap. x˜2,h = x4,h xextrap.

2 0.254768 3 0.259127 0.23788 4 0.256212 0.24859 1.52861 1.2814 5 0.254221 0.24911 1.39654 1.26398 6 0.253007 1.34308 1.25603 7 0.25224 1.31513 8 1.2984

∞ 1

4 = 0.25 5

4 = 1.25

TABLE I. Results of transfer matrix calculations for one and two chiral polymers (identical to results for two and four nonchiral polymers), on cylinders of circumference h hexagons. The largest eigenvalue of the transfer matrix is related to xL,h through (2.6). The convergence to the pre- dicted values x˜1 =

1

4 and x˜2 =

5

4 is seen to be quite rapid.

The extrapolated values are obtained by using three consec- utive values of xh to fix the constants in xh = c1 + c2h

−c3 , then taking c1 as an estimate of x∞.

polymer can be added from A to B in the nonchiral case. The subleading eigenvalues can differ, however, so there may not be a simple equivalence between states of L chi- ral polymers and 2L nonchiral polymers. The critical properties of the nonchiral model follow from Coulomb gas results for the O(n) model10, so we have

xL = L2 − 1 12

, x˜L = 4L2 − 1

12 . (2.7)

The chiral exponents x˜L are the same as those of the θ-point on the Manhattan lattice26. There are Monte Carlo results for another hexagonal lattice model believed to lie in the Manhattan universality class, the “smart kinetic growth SAW”27, which are consistent with the above values.

III. THE CLASSICAL LOCALIZATION

EXPONENT

When the particle energy E moves away from the crit- ical energy, the trajectories become less extended and the mean distance from the origin after N time steps is reduced. In this section, we take the lattice disorder average in order to express P (a, b,N), the probability that after N steps the particle has moved from a to b, as a weighted sum over linear and ring polymers. Simple properties of the weight function then yield an inequality connecting the localization exponent ν to polymer expo- nents γ and νθ. As in the preceding section, we fix an open or closed

curve on the lattice and ask what fraction of poten- tial realizations make this curve the correct trajectory.

The particle energy E is also varied in order to find the energy-integrated diffusion constant, and hence ν. The weight of a curve is determined by the numbers of hexagons touched by the curve to the left and right: the requirement for a path to be the correct trajectory is that all the hexagons to the right lie above the particle energy. The probability that a path is the particle trajectory is a function of the number of different hexagons visited by the path. The requirement is that all the hexagons to the immediate left have energies larger than the particle energy, while those to the immediate right have energies smaller than the particle energy. The HL hexagons on the left must have higher energies than the HR on the right, which is true for (HL+HRHL )

−1 of potentials. Fur- thermore, the particle energy must lie in the window of width ∼ (HL + HR)−1 between the lowest potential on the left and the highest potential on the right. So the weight of an allowed path PAB is

W (PAB) ∝ 1 (HL +HR)(

HL+HR HL

) ≈ HL!HR!

(HL +HR + 1)! .

(3.1)

The same result is obtained by integrating the probability (1+E2 )

HL(1−E2 ) HR over particle energy E to obtain a beta

function. Now we can again connect the expression (1.8) and the

weight (3.1) to known properties of polymers. For fixed H = HL+HR, the weight is minimized if HL = HR, and

using Stirling’s approximation is then W (H) ≈ √

2pi 2HH1/2

. The probability to get from a to b after N steps thus satisfies

P (a, b,N) ≥ ∑

SAPs througha and b, l= length of SAP, q=steps froma to b

δN mod l,q 2HH1/2

+ ∑

SAWsof lengthN froma to b

no ‖ self−contacts

1

2HH1/2 , (3.2)

up to a possible numerical constant. The typical number of hexagons H scales linearly in N to sufficient accuracy thatH1/2 can be replaced byN1/2 (this is verified numer- ically by Monte Carlo simulations, and if false would re- quire an unexpected multifractality at the θ-point). Sum- ming over final positions b to find the mean squared dis- placement then gives

〈〈R2(N)〉〉 ≥ Nγ−1+2νθ− 12 . (3.3) Then from equation (1.7) we obtain an inequality con- necting the localization exponent ν for the plateau tran- sition to polymer exponents νθ and γ:

1− 1 2ν ≥ γ + 2νθ − 3

2 . (3.4)

For the chiral polymer model, the resulting prediction is ν ≥ 1, which is satisfied by the actual value ν = 43 .

6

The usual nonchiral polymer exponents at θ would pre- dict ν ≥ 73 (cf. section IV), so again it is seen that the chirality constraint is essential. The fact that the lower bound is not reached shows that even in the limit of long paths, the number of hexagons to the left and right of the path cannot be assumed equal in calculating the weight (3.1). At the critical energy (section II), hexagons to the left and right contribute equally and this difference is ir- relevant, but away from the critical energy the difference affects the scaling. The exact value ν = 43 is derived in the percolation

picture4 from the equivalence of closed trajectories at energy E to percolation hulls with p−pc ∝ E−E0, where pc is the critical probability for percolation and E0 is the critical energy. Such percolation hulls28 have average size ξ ∼ (p − pc)− 43 . We remark in passing that the value ν = 43 can be understood in the polymer context from the fact the crossover exponent of the tricritical θ-point is φ = 37 (this value was first obtained using the connection

to percolation14): then N−φ ∼ (p − pc) and ξ ∼ (p − pc)

−νθ/φ = (p − pc)− 43 . The reasons for stressing the inequality (3.4) here rather than the exact result are that the inequality follows immediately from the classical path weight on polymers and can be used to gain information on higher moments (Section IV). The main result of this section is that the exact path

weight induced by averaging over disorder and particle energy can be calculated for the classical lattice model. This weight yields the inequality (3.4) connecting statis- tics of chiral polymers to the critical exponent ν of the classical plateau transition. The focus of the next section will be whether a similar relationship to polymers exists for the quantum plateau transition.

IV. HIGHER MOMENTS AND THE QUANTUM

TRANSITION

The quantum Hall plateau transition shows several qualitative similarities with the semiclassical limit stud- ied in the preceding sections of this paper. Both the quantum transition and the semiclassical limit have power-law delocalization at the critical energy and a fi- nite critical conductivity. However, the quantum tran- sition has proved much more difficult to describe theo- retically, and remains a major open problem. A natural question is whether any generalization of the polymer mapping developed for the classical limit would serve as a useful approach for the quantum case. The goal of this section is to show that the diffusive behavior usu- ally assumed to exist up to the localization length in the quantum case would be inconsistent with almost any such generalization, and present a numerical method and pre- liminary results to verify this assumption. We focus on one generalization in particular (to nonchiral polymers at θ, for reasons described below) for conciseness.

The scaling laws of moments of the particle distribu- tion function

〈R2n(t)〉 ∼ tα(n) (4.1)

demonstrate an essential difference between the classi- cal limit and the conventional picture of the quantum case. In this section 〈〉 indicates averaging over parti- cle energy and random potentials, while 〈〉E indicates averaging over random potentials at fixed particle en- ergy E. If particle motion is essentially diffusive on short length scales in the quantum case, then higher scaling laws beyond n = 1 in (4.1) do not contain additional in- formation. As discussed in the previous section, the mean squared displacement 〈R2n(t)〉 ∼ t1− 12ν contains the lo- calization exponent ν. This formula was obtained from the assumption that 〈R2(t)〉E increases linearly in time at each energy until the localization length is reached (R2 ≈ ξ(E)2), then saturates. If the particle motion is truly diffusive up to the local-

ization length, then 〈R2n(t)〉E ∼ tn until the localization length is reached, and

〈R2n(t)〉 ∼ tn− 12ν , (4.2)

with the localization length exponent ν ≈ 2.35 ± 0.056. So if motion in the quantum case is diffusive up to the localization length, there are no nontrivial exponents to be found in higher moments of the particle displacement. Higher moments in the classical case show nontrivial

scaling, and consequently highly extended trajectories are much more common in the classical case than the quantum case, even though the localization length di- verges more rapidly near Ec for the quantum case. For walks at the critical energy, it follows from the results of Section II that

〈R2n(t)〉Ec ∼ t2nνθ+γ−1 = t(8n−1)/7. (4.3)

Hence, although the mean square displacement does in- crease linearly with time, higher moments have nontrivial power laws because the particle trajectories are not ran- dom walks but instead have “memory,” as required for the absence of self-intersections. The higher moments are more extended than they would be for simple diffu- sive motion (random walks). A similar result holds for the classical case even when

the average is extended to include the particle energy. The inequality (3.4) derived in Section III between poly- mer exponents at the chiral θ-point and the scaling of the second moment α(1) = 1− 12ν from (4.1):

1− 1 2ν ≥ 2νθ + γ − 3

2 =

1

2 . (4.4)

For the classical transition with ν = 43 , the left side is 5 8 and the inequality is satisfied. Similarly for higher

moments 〈R2n(t)〉 ∼ tαc(n)

7

αc(n) ≥ 2nνθ + γ − 3 2 =

1

2 +

8(n− 1) 7

. (4.5)

Hence for sufficiently large n the classical scaling ex- ponents αc(n) are larger than the quantum exponents α(n) = n − 12ν , if the motion in the quantum case is indeed diffusive. In the remainder of this section, we consider the ques-

tion of how quantum interference keeps the quantum case from being related to a polymer ensemble in the same way as the classical case. One polymer ensemble in particular is attractive for the quantum case because the ensemble is similar to the classical one and the value ν = 73 appears in this ensemble, but this connection predicts nondiffusive motion up to the localization length. Monte Carlo nu- merics are used to test the assumption of diffusive motion in the quantum case; if verified this assumption would rule out a simple connection to polymers. The approach of the previous section was to generate a

positive weight on electron paths by averaging over disor- der with the electron path held constant. In lattice mod- els for the quantum transition, it should be possible to at- tribute a positive weight to each path on the lattice, and then these weights may be connected to some polymer problem, presumably different from the chiral ensemble discussed above for the classical limit. The first state- ment (that there is an assignment of weights) is some- what trivial from a mathematical point of view: there are many paths on the lattice between any pair of points, and hence given any positive probabilities to reach different points on the lattice after N steps, there is some assign- ment of positive path weights which results in the given probabilities. The difficult question is whether there is an assignment of weights which is physically meaningful and related to some local two-dimensional theory, as in the classical case. The next paragraphs define a real (not necessarily positive) path weight; this weight is positive in the absence of interference, and the leading interfer- ence contribution to this path weight from “cooperons” vanishes, though higher contributions do not. For fixed disorder, different paths W iAB from A to B

contribute to the amplitude, and the probability PAB to get fromA to B includes both diagonal terms |W iAB|2 and cross terms. Here and in the following we assume a dis- cretized model for the quantum case, similar to the lattice model introduced previously for the classical limit. We start by considering two paths which do not cross: then the disorder average generates a random phase which cancels the cross terms, leaving only the (positive) di- agonal terms. The remaining question is what occurs for intersecting paths; this is known to be the case of interest for the quantum effects causing weak localization. If the cross terms did vanish, then in the discretized

model where λD is effectively zero, PAB would be a sum over (not necessarily self-avoiding) paths with some pos- itive weight, the “quantum path weight” (QPW). A real weight can be defined even if cross terms are present by adding to the direct contribution of each path half of all

its disorder-averaged cross terms with other paths. The QPW picture can break down if, for energies near the critical energy, the cross terms become large enough to drive the weight negative for long paths. However, the leading interference corrections vanish upon disorder av- eraging, and a finite strength of interference is required to drive the path weight negative, so it is seemingly pos- sible that a positive QPW exists for the paths near the critical energy which determine ν. This motivates the conjecture, tested in the remainder of this section, that the universal large-length-scale properties of the plateau transition may be related to those of some classical gener- alized polymer model (i.e., a sum over paths with positive weights), in similar fashion to the relationships found in sections II and III between the classical percolation limit and the chiral polymer model. The remaining step is to determine whether any uni-

versality class of classical polymers can reproduce the weights which follow from disorder averaging in the quan- tum case. It seems worthwhile to identify possibilities, since exact results have been obtained for many two- dimensional polymer models by Coulomb gas and CFT techniques. The QPW should give nearly the classical weight (3.1) to paths which are classically allowed or in- clude a small number of quantum tunneling events, but should not allow of order N tunneling events for an N - step path since then the motion is simply diffusive even away from the critical energy. A speculative possibil- ity for the quantum transition is the nonchiral polymer ensemble at θ. The expectation that quantum mechan- ics should allow some unfavorable steps (but fewer than ∼ N) matches the fact that a typical polymer in the nonchiral ensemble has some parallel self-contacts, but fewer than of order N . A simple argument that the num- ber of parallel self-contacts is subextensive (< N) is that a ring polymer has no parallel self-contacts, so that par- allel self-contacts are in some sense a boundary property. A surprise is that the value ν = 73 (which has attracted

attention as the simplest rational consistent with numer- ics) appears from exponents of the nonchiral ensemble. The inequality (3.4) connecting the localization exponent ν to nonchiral polymer exponents predicts, since γ = 87 ,

1− 1 2ν ≥ γ + 2νθ − 3

2 =

11

14 , (4.6)

or ν ≥ 73 . Hence the value 73 appears in the critical prop- erties of a polymer ensemble closely related to the poly- mer ensemble describing the classical plateau transition. The the lower bound is realized if paths have asymptot- ically the same number of hexagons to the left as to the right (HR ∼ HL), which should be a better approxima- tion for the less convoluted paths in the quantum case. As seen below, however, this inequality predicts scaling laws for higher moments which appear to be ruled out numerically in the quantum case. Hence quantum inter- ference seems to be relevant at the quantum transition even beyond the level of changing path weights.

8

We note in passing that the desired property of dif- fusion at the critical energy does not have any simple interpretation as a statement about the polymer ensem- ble, since it is only after integrating over particle energy that the polymer weights may appear. This situation is familiar from the Liouvillian approach to the transi- tion16,31, where the transition is mapped onto a different problem which contains the exponent ν but not the crit- ical conductivity. Note that the previous appearance of the value 73 in a semiclassical average over a single per-

colation trajectory29 does not clearly relate to a critical point, when the electron is delocalized over multiple tra- jectories. Now we discuss how numerics can test whether the appearance of this value in this polymer ensemble is just a numerical coincidence. There are numerically testable consequences which

can be used to check whether nonchiral polymers at θ are indeed related to the quantum case. The polymer problem predicts various moments of the disorder- and energy-averaged displacement: it was shown in the pre- vious paragraph that the mean squared displacement 〈R2(N)〉 ∼ N11/14, so that ν = 73 . Similar predic- tions follow for higher moments, such as 〈R4(N)〉 ∼ N4ν+γ−3/2 = N27/14, or α(2) = 2714 . This can be com- pared to the null hypothesis of diffusion up to the local- ization length, which predicts 〈R4(N)〉 ∼ N1.78. We have performed Monte Carlo simulations with up

to 1800 states in the lowest Landau level to track the evo- lution of a localized wave packet in a disordered potential (the method is similar to that of30). The error bars are larger for the fourth moment than for the second because finite-size effects are more pronounced on the extended paths which dominate the fourth moment, but it appears that α(2) = 1.8±0.1, which if correct is sufficient to rule out α(2) = 27/14 ≈ 1.92. With larger system sizes, it should be straightforward to confirm the assumption of diffusive motion up to the localization length. Then a polymer description would have to have ν = 1/2 (either dense polymers or random walks), but no appropriate ensemble is obvious. It seems more likely that strong quantum interference prevents a physically meaningful assignment of path weights in the quantum case. To summarize, this section discussed differences be-

tween the classical and quantum transitions which be- come apparent in higher moments of the particle dis- tribution function. Numerics seem to support the pic- ture of diffusion up to the localization length and rule out the simplest polymer model for the quantum case. In closing, we mention briefly connections between the polymer models discussed in this paper and conformal field theory (CFT) approaches to the transition. The low-temperature O(n) phase also appears in a large-N expansion of the disorder-averaged Liouvillian theory, which is similar to a partially supersymmetric complex O(N) model, with N → 1 the physical limit31. Recent work32,33 on the critical point of Dirac fermions in a non- abelian random vector potential found a c = −2 dense

polymer problem hidden in the critical theory for sev- eral different types of disorder. Since upon adding addi- tional disorder (random mass and chemical potential) the abelian version of this critical point flows to the plateau transition fixed point34, the appearance of polymer sub- algebras may be generic to this class of random critical points. The author wishes to thank H. Saleur and S. Girvin

for helpful suggestions.

1 R. B. Laughlin, Phys. Rev. B 23, 5632 (1981). 2 A. M. M. Pruisken, Nucl. Phys. B 235, 227 (1984). 3 M. Zirnbauer, J. Math. Phys. 38, 2007 (1997). 4 S. Trugman, Phys. Rev. B 27, 7539 (1983). 5 I. A. Gruzberg, A. W. W. Ludwig, and N. Read, Phys. Rev. Lett. 82, 4524 (1999).

6 B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 7 J. T. Chalker and P. D. Coddington, J. Phys. C 21 2665 (1988).

8 S. Das Sarma, in Perspectives in Quantum Hall Effects, ed. S. Das Sarma and A. Pinczuk (Wiley, New York, 1997).

9 P. G. de Gennes, Phys. Lett. A38, 339 (1972). 10 B. Nienhuis, Phys. Rev. Lett. 49, 1063 (1982); B. Nienhuis, in Phase Transitions and Critical Phenomena, ed. C. Domb and J. L. Lebowitz (Academic, London, 1986), vol. 11.

11 F. Evers, Phys. Rev. E 55, 2321 (1997). 12 V. Gurarie and A. Zee, Int. J. Mod. Phys. B 15, 1225 (2001).

13 A. Coniglio et al., Phys. Rev. B 35, 3617 (1987). 14 B. Duplantier and H. Saleur, Phys. Rev. Lett. 59, 539 (1987).

15 H. P. Wei, D. C. Tsui, M. A. Paalanen and A. M. M. Pruisken, Phys. Rev. Lett. 61, 1294 (1988); H. P. Wei, S. Y. Lin, D. C. Tsui, and A. M. M. Pruisken, Phys. Rev. B 45, 3926 (1992).

16 J. Sinova, V. Meden, and S. M. Girvin, Phys. Rev. B 62, 2008 (2000).

17 R. M. Bradley, Phys. Rev. A 41, 914 (1990). 18 J. Cardy, Phys. Rev. Lett. 84, 3507 (2000). 19 A. Trovato and F. Seno, Phys. Rev. E 56, 131 (1997). 20 T. Prellberg and B. Drossel, Phys. Rev. E 57, 2045 (1998). 21 J. L. Cardy, Nucl. Phys. B 419, 411 (1994). 22 B. Derrida, J. Phys. A 14, L5 (1981). 23 H. Saleur, J. Phys. A 20, 455 (1987). 24 C. Vanderzande, Lattice Models of Polymers (Cambridge, 1998).

25 S. L. M. de Queiroz and J. M. Yeomans, J. Phys. A 24, L933 (1991).

26 T. Prellberg and A. Owczarek, J. Phys. A 27, 1811 (1994). 27 D. Bennett-Wood, A. L. Owczarek, and T. Prellberg, Phys- ica A 206, 283 (1994).

28 M. den Nijs, J. Phys. A 12, 1857 (1979). 29 G. V. Mil’nikov and I. M. Sokolov, JETP Lett. 48, 536 (1988).

9

30 S. Boldyrev and V. Gurarie, cond-mat/0009203 (2000). 31 J. E. Moore, A. Zee, and J. Sinova, Phys. Rev. Lett. 87, 8704-8707 (2001).

32 A. LeClair, cond-mat/0011413 (2000). 33 M. J. Bhaseen, J.-S. Caux, I. I. Kogan, and A. M. Tsvelik,

cond-mat/00012240 (2000). 34 A. W. W. Ludwig, M. P. A. Fisher, R. Shankar, and G. Grinstein, Phys. Rev. B 50 7526 (1994).

10

Comments