# Charge and Spin Transport at the Quantum Hall Edge of Graphene

###### Abstract

Landau level bending near the edge of graphene, described using 2d Dirac equation, provides a microscopic framework for understanding the quantum Hall Effect (QHE) in this material. We review properties of the QHE edge states in graphene, with emphasis on the novel phenomena that arise due to Dirac character of electronic states. A method of mapping out the dispersion of edge states using scanning tunneling probes is proposed. The Zeeman splitting of Landau levels is shown to create a particularly interesting situation around the Dirac point, where it gives rise to counter-circulating modes with opposite spin. These chiral spin modes lead to a rich variety of spin transport phenomena, including spin Hall effect, spin filtering and injection, and electric detection of spin current. The estimated Zeeman spin gap, enhanced by exchange, of a few hundred Kelvin, makes graphene an attractive system for spintronics. Comparison to recent transport measurements near is presented.

## I Introduction

Isolation and gating Novoselov04 () of graphene, a monolayer of graphite, has enabled the observation of interesting transport effects, resulting from Dirac fermion-like character of excitations. In particular, graphene hosts an integer quantum Hall effect (QHE) with unusual plateau structure Novoselov05 (); Zhang05 (). It was found that the QHE plateaus in monolayer are arranged symmetrically around the neutrality point, occurring at filling factors which are half-integer multiples of four, which is the combined spin and valley degeneracy of graphite. This form of QHE came to be known as the anomalous or half-integer QHE.

The simplest framework allowing to understand this behavior of QHE is provided
by the structure of Landau levels (LL)
of a 2d massless Dirac equation ^{0}^{0}footnotemark: 0
which has the particle-hole symmetric spectrum
^{0}^{0}footnotetext:
The result (1) follows from the
well-known representation
of the Schrödinger-Pauli Hamiltonian as a square of a massless Dirac Hamiltonian,
.
The Dirac spectrum (1) is then obtained
by taking a square root of the nonrelativistic Landau level spectrum.

(1) |

with the sign plus for positive , and minus for negative . Here is the magnetic field and is the velocity at the graphite Dirac point or . Due to spin and valley degeneracy (we shall discuss the role of Zeeman splitting below), each of the levels (1) contributes a step of to the quantized Hall conductivity. The particle-hole symmetry of the spectrum (1), with the level positioned at , suggests that the QHE plateaus must occur at , which is indeed what is observed in experiment Novoselov05 (); Zhang05 (). The anomalous QHE in graphene can be understood in a more fundamental way in terms of a quantum anomaly of the zeroth Landau level Gusynin05 (). The special character of the monolayer spectrum (1) is underscored by the difference between the QHE properties observed in graphene monolayer and bilayer systems Novoselov06 ().

One of the most dramatic consequences of the Dirac LL spectrum (1) is the appearance of the new energy scale . The square root dependence leads to a much larger level spacing than that for electrons with the quadratic dispersion conventional for semiconductors. For typical magnetic field of , the separation of the lowest LL () is quite large, , which enables QHE to persist up to room temperature GeimKimAPS2007 ().

To gain insight into the microscopic origin of the anomalous QHE, it is useful to develop the edge-states approach, which provides an intuitive and simple picture of the conventional QHE Halperin (). The edge states for graphene were studied using a numerical solution of the tight-binding model Peres () and also with the help of the Dirac equation Abanin06a (); Brey06 (). It was found that the energy levels (1), valley-degenerate in the graphene bulk, are split near the edge due to valley mixing at the boundary. Interestingly, it turns out that the structure and dispersion of the edge states depend on particular crystallographic orientation of the edge. For the so-called armchair edge, a simple particle-hole symmetric splitting was found Abanin06a (). Particle-like and hole-like states have different sign of energy dispersion, giving rise to counter-propagating modes with opposite chirality. A somewhat more complicated situation occurs near the zigzag edge, where dispersing edge states coexist with an additional dispersionless surface state Peres (). Despite these differences, the armchair, zigzag and other edges have same numbers of dispersing edge states of both chiralities. As we discuss below in Sec.IIII, this ensures the universal half-integer character of QHE in graphene.

One unique aspect of QHE in graphene is that its electronic states, owing to the monolayer character of this material, are fully exposed and, similar to the surface states of 3d materials, can be investigated by scanning tunneling microscopy (STM) probes with atomic resolution. Moreover, some useful information can even be obtained by imaging the top layer of 3d graphite, as demonstrated by recent STM studies of Dirac Landau levels Matsui05 (); Niimi06b () and of electron states near atomically thin edges Kobayashi06 (); Niimi06a (). In Sec.IIII we discuss new possibilities for scanning experiments that arise in graphene. The characteristic spatial scale of the states in a Landau level, given by the magnetic length , is about for the field of . Being large compared to the STM spatial resolution, it allows to image individual electronic states with sub- resolution and, in particular, to study Landau level bending near graphene edge. This bending in fact mimics the edge states momentum dispersion, due to the position-momentum duality of the Landau levels. As we shall see, the STM technique has sufficient resolution to map out the dispersion of QHE edge states.

Another novel feature of graphene is the simultaneous presence of the QHE edge modes of opposite chiralities, propagating in the opposite directions. Being particle-like and hole-like, they occur at the energies and , respectively. For electron density detuned from the neutrality point, , only one of the chiralities contributes to transport. However, as we shall see, near the states of both chiralities can participate in transport, leading to rather unusual transport properties. In particular, in the presence of Zeeman spin splitting of Landau levels, the state features an energy gap in the bulk and, simultaneously, a pair of edge states of opposite chirality and opposite spin polarization Abanin06a (); Fertig06 (). These states carry spin-up and spin-down electrons in the opposite directions along graphene edge, exhibiting quantized spin Hall effect but no charge Hall effect, owing to the particle-hole symmetry at . As we discuss in Sec.IIIIII, these spin-polarized chiral edge states exhibit interesting spin transport phenomena, such as spin filtering and spin injection, whereas the spin Hall effect provides a natural tool for the detection of spin current.

Interestingly, the counter-propagating edge states manifest themselves directly in charge transport. As we discuss in Sec.IVIV, near longitudinal resistivity remains finite, , while the Hall effect, which is absent at the particle-hole symmetry point , appears at nonzero due to conductivity in the bulk. This leads to Hall resistance changing sign at without exhibiting a clear plateau. The bulk conductivity short-circuits the edge transport and suppresses longitudinal resistivity, leading to a prominent peak in near and a plateau in These predictions, as well as the behavior of resistance fluctuations, which are enhanced near zero , are in agreement with experiment Zhang06 (); Abanin07 ().

Looking beyond graphene, interesting massless Dirac fermion states have been predicted a while ago Volkov85 () at interfaces of narrow-gapped HgTe and PbTe semiconductors. In Ref.Volkov85 (), which anticipated many of the features of electronic states in graphene, 2d Dirac states occur in a band-inverted heterojunction plane due to spin orbit interaction in 3d bulk, in the absence of magnetic field.

We note also that there have been interesting predictions of quantized spin Hall effect in certain 2d insulators KaneMele (); Bernevig06 (). In particular, the proposal of Ref.Bernevig06 () can be viewed as a 2d version of the 3d situation discussed in Ref.Volkov85 (). In these schemes an energy gap forms in the bulk due to spin-orbit even in the absence of magnetic field, similar to Ref.Volkov85 (). At the same time, counter-propagating modes carrying opposite spins appear on the edge, which are responsible for the quantized spin Hall current. As was emphasized in Ref.Fu06 (), there are general symmetry requirements, rooted in the time-reversal symmetry, protecting counter-propagating gapless excitations at an insulator boundary. In particular, certain invariants must exist, which are realized as spin component in the case of graphene KaneMele (); Abanin06a (), and were linked to supersymmetry for heterojunction systems Volkov85 ().

We note that the nature of the edge states of Refs.KaneMele (); Bernevig06 () is essentially identical to ours. The only difference, apart from different size of the energy gap, is that Rashba spin-orbit may cause backscattering in our case (see Sec.III) but not in the situation of Refs.Volkov85 (); KaneMele (); Bernevig06 (), while magnetic impurities will cause backscattering in both cases. In particular, our discussion in Sec.IV of and and our proposals to operate and detect spin current can serve as diagnostic tools should these schemes become realized experimentally.

The rest of the paper is organized as follows. In Sec.II we introduce the edge states using the Dirac equation framework. We focus on the two main edge types, armchair and zigzag, however we emphasize the generic aspects that are applicable to other edges. Then we discuss the anomalous QH effect and the possibility of imaging the edge states with scanning tunneling probes. In Sec.III we focus on spin-polarized chiral edge states and related spin transport phenomena. We also comment on spin relaxation mechanisms and present estimates of the spin relaxation time. In Sec.IV we introduce a transport model which accounts for both edge and bulk transport. This allows us to connect the chiral spin-polarized edge picture with recent transport measurements near the neutrality point Zhang06 (); Abanin07 ().

## Ii II. Dirac QHE edge states

Here we analyze electron states near zigzag and armchair edges (Fig.1a), the two most common graphene edge types, using massless Dirac model diVincenzo (). This exercise, which amounts to setting the boundary conditions for the Dirac spinor and solving an appropriate 1d eigenvalue problem, provides a fully microscopic picture of graphene QHE. Below we use this approach to illustrate the interplay between the QHE edge states and surface states for the zigzag edge, and discuss the possibility of imaging the edge states with scanning tunneling probes.

Let us first recall how low-energy graphene excitations are obtained in the tight-binding model diVincenzo () near the Dirac valleys and , located at the two non-equivalent Brillouin zone corners (see Fig.1(b)). There are two linearly independent zero-energy Bloch functions for each of the points , , each residing only on one sublattice ( or ) and vanishing on the other sublattice. Our choice of Bloch functions for the valley is shown in Fig.1(c). The Bloch functions for the valleys and are related by complex conjugation.

The wave function of
low-lying excitations near and ,
is written as a superposition of these four zero-energy Bloch functions
multiplied by slowly
varying envelope functions , ,
, ,
with and being the wave function amplitudes on the sublattice
and .
(Our choice of the signs for the valley is
convenient for treating an armchair boundary,
as we shall see below.)
The envelope functions , , and ,
describe excitations near and , respectively.
The effective low-energy Hamiltonian, obtained by keeping only
lowest-order gradients of and , takes the massless Dirac form
^{1}^{1}1We note that the Hamiltonian used in our
paper Abanin06a () is related to the Hamiltonian
(2)
by variable transformation ,
.:

(2) |

where , .

Landau levels in an external -field, described with the gauge , can be obtained for the states with the dependence from 1d Hamiltonians

(3) |

where and . Here and are measured in the units of and , respectively. The spectrum of , Eq.(3), yields the Dirac Landau levels, Eq.(1), where the eigenstates for the two valleys are given by

(4) | |||

(5) |

Here
are the eigenfunctions of the magnetic oscillator, ,
the normalization factor equals 1 for and for
, with and .
Note that zeroth LL states reside solely on sublattice
for valley and on sublattice for valley
^{2}^{2}2
This property is specific for the zeroth LL,
making the splitting of the LL due to Coulomb interaction
distinctly different from that of other LLs (see Abanin06b (); Goerbig06 ())..

We now analyze how LL spectrum is modified near the armchair edge. We consider graphene sheet in the halfplane with an armchair edge parallel to the axis (see Fig.1(a)). Energy levels near the edge are determined from the Dirac eigenvalue equations , where , and are given by Eq.(3). To analyze this eigenvalue problem, we exclude components and consider eigenvalue equations with spectral parameter for components:

(6) |

The boundary conditions for Eqs.(6) can be obtained from the tight-binding model, which is valid up to the very last row near , by setting the wave function equal zero at the boundary. Since the armchair edge has lattice sites of both and type (see Fig.1(a)), the wave function on both sublattices should vanish at the edge. In terms of the envelope functions , , taken at , this condition translates into

(7) |

We obtain a pair of differential equations (6) on the semi-axis , coupled at the boundary via Eq.(7). To simplify this problem, let us map Eq.(6) for onto the positive semi-axis , by , while keeping on the negative semi-axis , and treat it as an eigenvalue problem in the domain with the wavefunction given by at negative and by at positive . The first boundary condition, Eq.(7), then means that the wavefunction is continuous at , while the second condition implies continuity of the derivative . (This can be seen by expressing in terms of using Eqs.(3).) Thus we obtain a 1d Schrödinger problem in the potential Abanin06a ()

(8) |

defined on the entire axis. After finding the spectrum numerically, we obtain the energy levels of the Dirac fermions as

(9) |

whereby the particle-hole symmetry is restored due to the two possible signs of . The dispersion (9) is illustrated in Fig.2a. We note that the LL double valley degeneracy in the bulk is lifted at the boundary.

The particle-hole symmetric edge states spectrum in Fig.2a instantly explains the “half-integer” Hall quantization in graphene. Indeed, for any electron density with integer in the bulk there is an odd number of the edge modes crossing the Fermi level, which means that the Hall conductivity is quantized as , where the factor two accounts for spin degeracy.

Spin degeracy of the Landau levels is lifted by the Zeeman interaction, which is substantial in graphene,

(10) |

for (compare to in GaAs quantum wells). Zeeman-split edge states, depicted in Fig.2b, have interesting characteristics for electron density near neutrality, . At this density the state in the bulk is spin-polarized, with the Zeeman gap further enhanced by exchange (see below). The two branches of counterpropagating edge states near , carrying opposite spin, have interesting properties that will be discussed in more detail in Sec.III.

We now analyze the zigzag edge, which even at hosts a band of dispersionless zero-energy states bound to the edge Fujita (). We shall refer to these states as surface states, to distinguish them from the dispersing QHE edge states. As we shall see, the surface states contribute to the splitting of LL near the zigzag edge, in agreement with the tight-binding calculations Peres ().

We consider graphene sheet in the half-plane , with its first row consisting of atoms (see Fig.1a). For the states with the dependence , with , , from (2) we obtain 1d Hamiltonians

(11) |

where . Similarly to the armchair case, the spectrum can be found from the eigenvalue equation , where , which should be supplemented with the boundary conditions. For our zigzag edge the wavefunction have must vanish on all sites at . For that both envelope functions have to vanish at the boundary,

(12) |

Excluding components, we obtain two separate eigenvalue problems for the spectral parameter ,

(13) |

where both and satisfy the hard wall boundary conditions (12). The amplitudes on the sublattice can be expressed via ampliudes on the sublattice and eigenenergy ,

(14) |

The eigenvalue problems (13) with the hard-wall boundary conditions (12) are familiar from the theory of edge states in the conventional QHE Halperin (), and their spectrum can be found numerically. The Dirac fermion energy dispersion is shown in Fig.3.

The behavior of LL’s is similar to the armchair case: there are two branches of the edge states, one for each valley, degenerate in the bulk, , which split near the edge. The zeroth LL, however, coexists with the surface state, which makes its behavior rather peculiar and different for the two valleys.

In the valley, which we discuss first, the zeroth LL states reside solely on the sublattice, see Eq.(4), and therefore automatically satisfy the boundary condition . Thus there are zero-energy states for arbitrary values of , of the form . Let us consider the states with far outside the graphene half-plane, . Not too far from the boundary, such states can be approximated by an exponential

(15) |

which is identical to surface state wave function Fujita (). Thus the zeroth LL for valley near the edge transforms into the surface mode. Being dispersionless, this mode does not contribute to the edge transport. The edge state spectrum for the valley is displayed in Fig.3a.

Now let us consider the zeroth LL for the valley . For , we approximate the ground state of the oscillator (13) with the hard-wall boundary condition as

(16) |

The ground state energy is then approximated by

(17) |

where is the effective Hamiltonian for component, Eq.(13), and denotes averaging over the normalized wave function (16).

Evaluating for , when the state (16) has unit norm with exponential accuracy, we obtain

(18) |

From the relation , we find the energies for the two branches of dispersing edge states

(19) |

Plugging this expression in Eq.(14), we obtain the wave function on sublattice for these two branches,

(20) |

which is again the surface state wave function (compare to Eq.(15)). We therefore conclude that for the valley the zeroth Landau level and the surface state mix giving rise to two dispersing edge modes. This is in agreement with the spectrum displayed in Fig.3b. We see that, although the Dirac model is applicable only in a small part of the Brilloin zone, near points and , it provides a description of the states at the zigzag edge, including the surface state, which is in agreement with the results of the tight-binding model of Ref.Peres (). The surface mode in the vicinity of and is given by Eqs.(15),(20).

Interestingly, the and sites contribute equally to the splitting of the zeroth LL, . This is somewhat counterintuitive, since this LL is solely on the sublattice in the bulk, while the surface mode is solely on the sublattice . This equal participation property can be understood as follows. The spinor states with , given by Eqs.(16),(20), are eigenstates of the Dirac Hamiltonian with the boundary condition (12), with the energies . Thus these states are orthogonal, which implies that the integrals of and are equal. We further note that the integral of the square of the component of our edge state wave function (20) over indeed equals one, in agreement with our choice of normalization on the sublattice in Eq.(16).

To sum up, for the zigzag edge, the zeroth LL gives rise to two dispersing edge states for one of the valleys, while for the other valley the zeroth LL morphs into the dispersionless surface mode which does not contribute to the edge current. Therefore, despite the presence of the surface mode, the number of dispersing QHE edge states with and for the zigzag boundary is the same as for the armchair boundary, giving rise to “half-integer” quantization of Hall current.

Finally, we briefly discuss how the edge states in graphene can be investigated using the STM technique Matsui05 (); Niimi06b (); Kobayashi06 (); Niimi06a (). Due to the Landau level momentum-position duality relation, , the edge state momentum dispersion shown in Figs.2,3 translates into the excitation energy dependence on the distance from the edge. The characteristic scale for the latter is set by the magnetic length , which for typical fields is about 50-80 times greater than the spatial resolution of STM instruments on graphite surface. This makes STM technique particularly convenient for this kind of studies.

A link between the edge states dispersion and the position-dependent tunneling spectroscopy can be established as follows. We shall use the solutions for the edge state wave function given above to calculate the local density of states (LDOS) near the zigzag edge (other edge types can be dealt with similarly). For each of the graphene sublattices LDOS is given by

(21) |

where is the distance from the edge, and denotes the set of eigenstates of the and Hamiltonians (11) with the hard-wall boundary condition (12).

Using the eigenfunctions , and the energies found from Eqs.(13),(14) as discussed above, we obtain LDOS for the and sublattices which is displayed in Fig.4. We see that the position-independent Landau level bands, dominating LDOS far from the edge, bend away from near the edge. This bending mimics the edge states momentum dispersion shown in Fig.3. Note, however, that LDOS is nonzero only for , whereas the edge state momentum can be both positive and negative. The spatial width of the bending bands is determined by the width of the eigenfunctions , , which is of the magnetic length scale.

## Iii III. Spin-polarized chiral edge states and spin transport.

As we noted above, at the neutrality point graphene hosts gapless spin-polarized edge states (see Fig.2(b)). The Zeeman energy gap in the bulk, Eq.(10), is enhanced by the Coulomb interaction. A Hartree-Fock estimate of this enhancement Abanin06a () gives

(22) |

where is RPA screening function, and the parameter describes relative strength of Coulomb and exchange correlations. Assuming , i.e. ignoring correlations of electrons with opposite spin, we obtain a spin gap for . Taking into account the substrate dielectric constant, , changes the result only slightly ( for ). This approximation, while pointing at a correct order of magnitude of a few hundred Kelvin, probably somewhat overestimates the spin gap since it ignores correlations and disorder effects.

The chiral spin-polarized edge states offer a unique setting to study spin transport. In particular, the spin-split state may be used to generate and detect spin-polarized currents. This spin transport regime seems attractive due to the large bulk gap and high stability of the edge states. Moreover, increased quality of samples should allow existence of spin polarized edge states even at relatively low magnetic fields.

Since the purpose of this section is mostly illustrative, we will keep our discussion as simple as possible. In particular, we shall ignore transport in the bulk, leaving the discussion of its role for Sec.IV. We also first neglect spin flip backscattering between edge states within one edge. Estimates of the spin flip rate wil be given below, Eq.(24). A general approach, based on the Landauer-Büttiker formalism Buttiker (), which can be used to calculate spin and charge currents at the edge for any configuration of current and voltage leads, was presented in Ref.Abanin06a ().

In this approach, transport is described by a scattering matrix Buttiker (), with the edge states playing the role of scattering channels, and the reservoirs supplying in-states and absorbing out-states. Current in each mode is described by the relation , where is the reservoir chemical potential for given spin projection. We consider the Hall bar geometry with four contacts 1-4 (see Fig.5), where the contacts 1 and 4 serve as current source and drain. For these two contacts we do not assume spin mixing, so that the injected and drained current may be spin polarized. The contacts 2, 3 are voltage probes, which means that they do not drain current from the system. Furthermore, we assume that the probes provide full spin mixing, i.e. chemical potentials of outgoing spin-up and spin-down electrons are equal.

The simplest situation arises when unpolarized current is injected through contact 1. Then the up- and down-spins spatially separate in a symmetric way, flowing along the opposite edges of the bar. This can be interpreted as circulating spin current, and described as spin-Hall effect with quantized spin conductance . No electric voltage will be induced between the voltage probes 2, 3 in this case (zero charge-Hall effect).

This device can be used as a detector of spin polarized current, made possible by the reciprocal of the spin Hall effect, in which the electric Hall voltage is directly related to spin rather than charge current. Suppose the up-spin and down-spin electrons, injected through contact 1, have unequal chemical potentials, . Then the currents flowing into the probes 2 and 3, , after equilibration and spin mixing in the probes, induce voltages . The resulting Hall voltage is directly proportional to spin current. At the same time, an unpolarized current (for which ) flows symmetrically in the upper and lower edges without generating Hall voltage.

Spin transport at also allows to realize spin filter. Suppose that the upper and lower edges of the device in Fig.5 are made asymmetric, which can be achieved, for example, simply by removing probe 3. Then we inject unpolarized current through contact 1. The injected current will be distributed equally between the upper and lower edges in cross section A. In cross section B, however, the net current will be spin polarized due to spin mixing in probe 2. The down-spin current reaching the drain in the upper edge equals while the up-spin current in the lower edge is . Therefore, the total drained current becomes spin polarized. The spin polarized current can be fed into another system (see Fig.5), where it can be detected using Hall probes 5 and 6 as discussed above.

More complicated circuits can be assembled which generate spin currents and detect them elsewhere. Note that the important principle is that as long as backscattering is not allowed, the edge current can travel long distances and the circuit is nonlocal, just as in the integer QHE Buttiker (). In this case the current-voltage relationship is obtained by solving the circuit equations as described in Ref.Abanin06a (). The spatial scale of nonlocality is controlled by spin relaxation which can be due to spin-orbit interaction or due to magnetic impurities near graphene edge.

For simplicity, here we limit the discussion to the effects of spin-orbit. There are two main spin-orbit terms in the graphene Hamiltonian KaneMele (); Min06 (), the so-called intrinsic and Rashba interaction, given by

(23) |

where Pauli matrices act in the sublattice space (Dirac spinor), while act in the valley space, and represent physical spin. Estimates from band calculations Min06 () give and a negligibly small .

To estimate the backscattering rate due to the spin-orbit interaction (23), we note that for an ideal atomically sharp edge the spin-orbit would couple the left and right states with the same momentum, opening a minigap at branch crossing: . However, this momentum-conserving interaction alone cannot backscatter edge states, and we need to take disorder into account. Edges of graphite monolayers have been imaged using STM probes Kobayashi06 (); Niimi06a (), where it was found that typically edge disorder can be viewed as patches of missing atoms of characteristic size . Taking into account the left-right branch mixing by spin-orbit , characterized by small mixing ratio of away from branch crossing, we obtain an estimate of the backscattering mean free path:

(24) |

which gives for typical . The factor accounts for the magnetic field dependence of disorder matrix elements.

The quadratic energy dependence in (24), with spin flip rate having a sharp peak near branch crossing, suggests Abanin06a () the possibility to control backscattering using local gate. By tuning local chemical potential to and from the branch crossing, where the spin flip rate has a sharp peak, Eq.(24), we can induce or suppress backscattering in a controlled way. Spin filtering is achieved by controlling local gates on opposite sides of the Hall bar asymmetrically.

## Iv IV. Edge and bulk transport at .

Spin flip backscattering (24) can be incorporated in the edge transport model, described by coupled equations for particle density in the two spin-polarized modes:

(25) |

where is the backscattering mean free path (24) taken for at the Fermi level, and are compressibilities of the modes. (For brevity, we use 1 and 2 instead of and .) In writing Eqs.(25) we implicitly assume that fast energy relaxation maintains local equilibrium of each of the modes, which is consistent with metallic temperature dependence of transport coefficients Abanin07 ().

In a stationary state, Eqs.(25) have an integral which expresses current conservation at the edge. [In this section we use the units of .] The general solution in the stationary current-carrying state is

(26) |

Taking into account that is the current in one edge, we calculate the total current as

(27) |

To describe the longitudinal resistance in the four-terminsl geometry, one must add potential drop on voltage probes Abanin07 (), which gives , where is the distance between the probes. Comparing to the data for at we estimate Abanin07 () . This mean free path value, which is relatively small on the scale predicted by Eq.(24), can be explained if spin flip processes are dominated by nonintrinsic effects, such as magnetic impurities localized near the edge.

It is crucial that the edge transport model (25) treats both edges of a Hall bar in an identical way, thus predicting zero Hall effect. In order to understand the observed density dependence of Hall coefficient Zhang06 (); Abanin07 (), which changes sign smoothly at without exhibiting a plateau, and of which has a sharp peak at , we need to incorporate transport in the bulk in our model. In the full edge+bulk model, the density dependence of transport coefficients arises from bulk currents short-circuiting edge currents away from . This explains, as we shall now see, the Hall effect, the peak of , the resistance fluctuations near , as well as the behavior of and .

We describe the transport problem in the bulk by the current-field relation, separately for each spin projection:

where are electrochemical potentials for two spin states. We assume that the bulk conductivities , as a function of density , are peaked at the spin-split Landau levels. For simplicity, here we ignore possible valley splitting, in which case the spin up and down Landau levels occur at around the Dirac point. As a simplest model, below we use Gaussians

(28) |

with the parameter describing the width of the levels. The Hall conductivities exhibit plateaus on either side of the peak in . The dependence of on can be modeled with the help of the semicircle relation which often provides a good description of conventional QHE systems semicircle_relation ().

The condition of charge continuity, , gives a 2d Laplace’s equation for the potentials, . This equation must be solved together with the boundary conditions phenomenologically describing bulk-edge coupling:

(29) |

where is a normal vector to the boundary, and represents the edge-bulk leakage current density.

Although a general solution of this problem can be given with the help of Fourier method, here we consider only the case when the potentials are varying slowly on the scale of the bar width , which will suffice for our analysis of a homogeneous current flow. In this case, linearizing in the direction transverse to the bar, we can write Eqs.(29) for both edges of the bar as

(30) |

, where the primed and unprimed quantities denote variables at opposite edges of the bar.

Equations for the edge variables are obtained by adding the bulk-edge leakage term to Eqs.(25), giving

(31) |

along with a similar pair of equations for , at the opposite edge.

The solution of these eight equations, describing uniform current, is of the form , , etc., with the same linear part for all quantities. Using the algebraic structure of this linear system and the symmetry between the edges, we reduce the number of equations from eight to two. First, it is convenient to express the parameters through using Eqs.(31), which gives

(32) |

Writing similar equations for the variables at the opposite edge to express , through , , and substituting the result in Eqs.(30), we obtain four equations for and which have the form

(33) |

where the coefficients in this linear system are defined as

(34) |

The quantities represent the sum of the bulk and edge contributions to Hall conductivity for each spin.

Symmetry between the edges allows to further reduce the number of independent variables. For that we add the first two equations to obtain . Also we note that all potentials can be changed by the same constant that can be chosen so that the new quantities and satisfy . After that Eqs.(33) yield

(35) |

These two equations can be solved to find .

Now we can find the current as a sum of the edge and bulk contributions, , where

and

After expressing through with the help of Eqs.(32) and using the solution of Eqs.(35), we obtain a relation , where

(36) |

The quantities are defined as . The quantity , Eq.(36), replaces in Eq.(27). In the absence of bulk conductivity, , we recover the result for pure edge transport, .

The Hall voltage can be calculated from this solution as , where , are variables at opposite edges. We obtain , where

(37) |

This quantity vanishes at , since and at this point due to particle-hole symmetry.

Transport coefficients, obtained from this model for typical parameter values, are displayed in Fig.6 which reproduces many of the key features of the data (see Fig.1 in Ref.Abanin07 ()). In particular, the peak in is due to edge transport near . The suppression of at finite is due to the bulk conductivity short-circuiting the edge transport. The bulk and edge contributions to transport can be discerned from the double peak structure in in Fig.6. The peaks correspond to the bulk Landau level contributions, Eq.(28), whereas the part of between the peaks, exceeding the superposition of two Gaussians, Eq.(28), is the edge contribution. The Hall resistance is nonzero due to imbalance in for opposite spin polarizations away from . Interestingly, in Fig.6 exhibits no plateau, while calculated from and displays an under-developed plateau-like feature. Overall, this behavior resembles that of the experimentally measured transport coefficients Zhang06 (); Abanin07 ().

Another notable feature of the measured and is enhanced fluctuations near zero . These fluctuations are found to be strong in Ref.Zhang06 (), where changes sign several times near . They are also present, although are not as dramatic, in Ref.Abanin07 (). In the latter case, both and exhibit noisy behavior in the interval near comparable to the peak width. As Ref.Abanin07 () points out, this behavior is consistent with the edge transport model. In the absence of bulk transport, the distribution of potential along the edge depends on the local backscattering rate , whereby Eq.(26) is replaced by

Fluctuations of arise due to its sensitivity to the local value of Fermi energy in the spin-orbit scattering model, Eq.(24), and, similarly, for the magnetic impurity scattering mechanism. Assuming that the random part of is of a white noise character, we obtain strong fluctuations along the edge of magnitude that scales as a square root of the edge length. These fluctuations will contribute equally to the longitudinal and transverse voltage, since they are uncorrelated on the opposite sides of the Hall bar. The absence of fluctuations away from can be understood as a result of bulk conductivity short-circuiting the edge current, which will equilibrate potentials on the opposite sides of the Hall bar.

The above discussion summarizes the results drawn from an attempt to model quantum Hall transport in graphene at by counter-circulating edge states. By taking into account backscattering within one edge as well as conduction in the bulk which short-circuits edge transport away from the neutrality point, this model accounts for the observed behavior of transport coefficients. Still, since no direct evidence for spin polarization has yet been found, more experimental and theoretical work will be needed to confirm the chiral spin-polarized edge picture of the state. If proven to exist in graphene, these states will provide a unique setting to study spin transport as well as other interesting phenomena.

This work is supported by NSF MRSEC Program (DMR 02132802), NSF-NIRT DMR-0304019 (DA, LL), and NSF grant DMR-0517222 (PAL).

## References

- (1) K. S. Novoselov et al., Science, 306, 666 (2004); Proc. Natl. Acad. Sci. USA, 102, 10451 (2005).
- (2) K. S. Novoselov et al., Nature 438, 197 (2005);
- (3) Y. Zhang, Y.-W. Tan, H. L. Stormer and P. Kim, Nature 438, 201 (2005).
- (4) V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 146801 (2005).
- (5) K. S. Novoselov et al., Nature Physics 2, 177 (2006).
- (6) K. S. Novoselov et al., Science 315, 1379 (2007).
- (7) B. I. Halperin, Phys. Rev. B25, 2185 (1982).
- (8) N. M. R. Peres, F. Guinea, A. H. Castro Neto, Phys. Rev. B73, 125411 (2006).
- (9) D. A. Abanin, P. A. Lee and L. S. Levitov, Phys. Rev. Lett. 96, 176803 (2006).
- (10) L. Brey and H. A. Fertig, Phys. Rev. B73, 195408 (2006)
- (11) T. Matsui et al., Phys. Rev. Lett. 94, 226403 (2005)
- (12) Y. Niimi et al., Phys. Rev. Lett. 97, 236804 (2006)
- (13) Y. Kobayashi et al., Phys. Rev. B 71, 193406 (2005).
- (14) Y. Niimi et al., Phys. Rev. B 73, 085421 (2006).
- (15) H. A. Fertig and L. Brey, Phys. Rev. Lett. 97, 116805 (2006).
- (16) Y. Zhang et al., Phys. Rev. Lett., 96, 136806 (2006).
- (17) D. A. Abanin et al., Phys. Rev. Lett. 98, 196806 (2007).
- (18) B. A. Volkov and O. A. Pankratov, Pis’ma Zh. Eksp. Teor. Fiz. 42, 145 (1985) [Engl. transl. JETP Lett. 42, 178 (1985)]; review in: O. A. Pankratov, Semicond. Sci. Technol. 5, S204-S209 (1990).
- (19) C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005).
- (20) B. A. Bernevig, T. Hughes, S.-C. Zhang, Science 314, 1757 (2006).
- (21) D. P. DiVincenzo and E. J. Mele, Phys. Rev. B 29, 1685 (1984).
- (22) L. Fu, C. L. Kane, cond-mat/0606336, unpublished.
- (23) M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, J. Phys. Soc. Jpn. 65, 1920 (1996).
- (24) M. Büttiker, Phys. Rev. B 38, 9375 (1988).
- (25) H. Min et al., Phys. Rev. B74, 165310 (2006).
- (26) A. M. Dykhne and I. M. Ruzin, Phys. Rev. B50, 2369 (1994); S. S. Murzin, M. Weiss, A. G. Jansen, and K. Eberl, Phys. Rev. B66, 233314 (2002)
- (27) M. O. Goerbig, R. Moessner, B. Doucot, Phys. Rev. B74 161407 (2006)
- (28) D. A. Abanin, P. A. Lee, and L. S. Levitov, Phys. Rev. Lett. 98, 156801 (2007).