Relaxation and Edge Reconstruction in Integer Quantum Hall Systems
Relaxation and Edge Reconstruction in Integer Quantum Hall Systems
This content has been downloaded from IOPscience. Please scroll down to see the full text.
(http://iopscience.iop.org/1367-2630/14/10/105009)
View the table of contents for this issue, or go to the journal homepage for more
Download details:
IP Address: 84.25.233.197
This content was downloaded on 31/08/2015 at 09:16
4
Author to whom any correspondence should be addressed.
Content from this work may be used under the terms of the Creative Commons Attribution-NonCommercial-
ShareAlike 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title
of the work, journal citation and DOI.
New Journal of Physics 14 (2012) 105009
1367-2630/12/105009+14$33.00 © IOP Publishing Ltd and Deutsche Physikalische Gesellschaft
2
Contents
1. Introduction 2
2. Unreconstructed edge 3
2.1. Three-body scattering formalism . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.2. Results for the unreconstructed edge . . . . . . . . . . . . . . . . . . . . . . . 6
3. Spin reconstruction 7
4. Charge reconstruction 8
5. Conclusions 10
Acknowledgments 11
Appendix A. Calculation of the Coulomb matrix element Vq (k1 − k2 ) 11
Appendix B. Calculation of the three-body matrix element T123
10 20 30 12
References 13
1. Introduction
The kinetic properties of one-dimensional (1D) quantum systems are an active area of current
research [1, 2]. What makes the field exciting is that many-particle physics is drastically
different in one spatial dimension. This difference is already evidenced in basic nonequilibrium
properties such as the microscopic mechanisms of relaxation. Within the scope of Fermi-liquid
theory, relaxation processes in higher dimensions proceed by pair collisions of electrons, which
provide an efficient mechanism for the relaxation of initial nonequilibrium states. In contrast,
conservation of energy and momentum strongly restricts scattering in one spatial dimension
so that pair collisions necessarily result in zero-momentum exchange or an interchange of
the momenta of the colliding particles. Neither process causes relaxation. This poses the
fundamental question of the microscopic origin of relaxation in 1D systems. Notably, the
absence of relaxation by pair collisions, which holds true regardless of the strength of
interaction, has received experimental support [3].
The question of equilibration emerges in a diverse set of 1D many-body systems. These
include energy and momentum-resolved tunneling experiments with nanoscale quantum wires
[3, 4], quench dynamics of cold atomic gases [5, 6] as well as energy-spectroscopy experiments
on quantum Hall edge states driven out of equilibrium [7–10]. This paper is motivated by the
latter experiments, which were carried out in a high-mobility two-dimensional electron system
at Landau-level filling factor ν = 2. This system hosts two co-propagating edge states, which can
be driven out of equilibrium by inter-edge tunneling in the vicinity of quantum point contacts.
This generates a nonequilibrium distribution of electron energy (in the sense of electronic edge
transport) downstream from the contact, which is monitored as a function of the propagation
distance by means of a quantum-dot-based energy spectrometer. The experiments show that
the initial nonequilibrium distribution relaxes to a stationary form that is close to the thermal
distribution but with an effective temperature and chemical potential.
Edge-state equilibration was also probed in experiments at Landau-level filling factor
ν = 1 [11]. Heat is carried unidirectionally by the single chiral edge mode as confirmed by
thermopower measurements along the edge. These experiments found that hot electrons injected
locally into the edge cool down while propagating along the edge. It is worth emphasizing that
2. Unreconstructed edge
A confinement potential Vc (x) that is sharp on the scale of the Coulomb interaction
(i.e. Vc0 e2 /(κl 2B ), where κ is the dielectric constant and l B denotes the magnetic length)
remains stable against interaction-induced reconstructions and the electron dispersion ε(k)
can be obtained approximately from the noninteracting Schrödinger equation [22]. A generic
electronic dispersion of an unreconstructed edge is sketched in figure 1(a), exhibiting a
confinement-induced bending of the Landau levels near the edge of the sample.
In the limit of high magnetic fields (Vc0 ωc /l B ), the electron states near the edge can be
described by the lowest Landau-level wave functions
√
ψ X (x, y) = (Ll B π )−1/2 e−iky e−(x−X ) /2l B
2 2
(1)
in the Landau gauge. Here, k = X/l 2B and L denotes the length of the sample edge (taken along
the y direction). The defining feature of the unreconstructed edge is the sharp zero-temperature
occupation function νσ (X ) = 2(−X ) of Landau-level states with the guiding center X when
the Zeeman splitting εZ is negligible (see figure 1(b)).
The single-particle dispersion near the Fermi energy (corresponding to momentum kF ) is
controlled by the confinement potential and can be approximated as
ε(k) = vF (k − kF ) + (k − kF )2 /2m c . (2)
The dispersion is parametrized through the edge velocity vF = Vc0l 2B at the Fermi energy and
the curvature 1/m c = Vc00l 4B . Note that these parameters become maximal for an infinitely sharp
edge for which vF ∼ ωcl B and 1/m c ∼ 1/m [22]. Note, however, that a description in terms of
the wave functions in equation (1) is no longer valid in this extreme limit.
The finite curvature of the dispersion implies that at least three particles are required
for an energy and momentum conserving relaxation process. Relaxation of a high-energy
electron (labeled by i = 1 in figure 1(a)) is possible by scattering two electrons (labeled i = 2, 3
in figure 1(a)) near the Fermi energy. Indeed, due to the curvature of the dispersion near
the Fermi energy, exciting electron i = 2 from the Fermi energy requires more energy than
scattering electron i = 3 deeper into the Fermi sea. Clearly, this relaxation process relies on
finite temperature, and the typical energy transfers for electrons i = 2, 3 at the Fermi energy
are of the order of T . Quantitatively, this process can relax the hot particle with excess
energy ε ≈ vF (k1 − kF ) = vF 1k by q1 = q3 (k20 − k3 )/1k, where qi = ki 0 − ki is the momentum
transferred to particle i in the collision. Note that q1 q3 so that relaxation occurs in many
small steps vF q1 ∼ T 2 /ε.
For Landau-level filling factor ν = 2, these considerations apply when the Zeeman splitting
is small compared to the temperature. In the opposite limit εZ T , the curvature of the
New Journal of Physics 14 (2012) 105009 (http://www.njp.org/)
5
dispersion implies that the Fermi momenta and hence the Fermi velocities differ for the two
spin directions. In this case, relaxation is dominated by processes in which the electrons i = 1
and 2 have opposite spins and thus different Fermi momenta kF j and Fermi velocities v j with
j = 1, 2. To include a finite Zeeman splitting at Landau-level filling factor ν = 2 as well as for
later convenience, it is thus beneficial to consider a modified dispersion
v (k − kF j ), k ≈ kF j ,
ε(k) = j (3)
v1 (k − k1 ) + ε, k ≈ k1 ,
which is linearized in the vicinity of each of the three particles, including the hot particle with
velocity v1 and momentum k1 . This captures the behavior in the regime of strong Zeeman
splitting εZ T , which we will focus on in the following. Nevertheless, we can also recover the
results for the quadratic dispersion and weak Zeeman splitting εZ T by identifying v2 − v3
with the typical velocity difference T /(vF m c ) due to the curvature of the dispersion.
Using the dispersion in equation (3), energy and momentum conservation leads to
v2 − v3
q1 = q3 . (4)
v1 − v2
The velocity difference v2 − v3 = εZ /(v2 m c ) is controlled by the Zeeman splitting, which we
assume to be small compared to the excitation energy ε such that (v2 − v3 ) (v1 − v2 ).
Since for large Zeeman energy the particles at k2 and k3 have opposite spins, there is no
exchange contribution (remember that exchange is appreciable for small momentum transfers
only) and equation (9) fully determines the three-body matrix element. The corresponding
energy relaxation rate can then be obtained by power counting, which yields
2 4
1 e (m c v12 )4 T 4
∼ m c (v2 − v3 )2
. (10)
τE κv2 ε8
Here 1/m c = Vc00l 4B , (v2 − v3 ) = εZ /(v2 m c ) and we also used 1k = m c (v1 − v2 ) = ε/v1 . In
obtaining equation (10), a factor of L/(v1 − v2 ) emerges from eliminating the energy δ-function
in equation (5), each summation over the remaining k2 , k3 , q3 contributes a phase space factor
of ∼ T /v2 and the weighting factor v1 q1 /ε takes the form v1 (v2 − v3 )T /[v2 (v1 − v2 )ε]. Finally,
we have to account for the competition between excitation (q1 > 0) and relaxation (q1 < 0) of
are no longer necessarily of opposite spin. T1123 0 30 20 can then be obtained from equation (9) by
replacing q3 → k20 − k3 , which does not change the power counting argument.
It is therefore possible to combine both cases by setting v2 − v3 = max{εZ , T }/(v2 m c ). In
the case of v1 ≈ v2 = Vc0l 2B , it is then possible to rewrite equation (10) as
4 4
max{εZ , T }2 Vc0 2 e2
1 T
∼ , (11)
τE Vc00l B2
Vc00l B κl B ε8
which applies for the regime Vc0l 2B ε(Vc00l 2B /Vc0 ) e2 /κ.
For a later comparison of the relaxation rates before and after edge reconstruction, it will
be useful to consider the unreconstructed case as the v1 v2 ∼ e2 /κ limit of equation (10)
(which can be applied for ε (e2 /κl B )2 /(Vc00l 2B )). Formally, this regime leaves the condition for
applicability of the Taylor expansion of the confinement potential that defines 1/m c = Vc00 (µ)l 4B
and would lead to another inverse mass Vc00 (µ + ε)l 4B for curvature effects at energies of the order
of ε. Distinguishing these different masses does not, however, lead to qualitative changes in the
results, and for brevity of the presentation we assume a quadratic confinement potential over the
energy interval [µ, µ + ε]. We can then rewrite equation (10) as
max{εZ , T } 2 00 2 T 4
1
∼ Vc l B , (12)
τE(u) e2 /(κl B ) ε
where we used that in this regime, m c v12 ∼ . The crossover between equations (11) and (12)
can be obtained at their limits of applicability by setting ε = e2 Vc0 /(κl 2B Vc00 ) and Vc0 = e2 /(κl 2B ).
Note that for a spin polarized edge, equations (10)–(12) apply only if the Coulomb
interaction is not screened for momenta of the order of T /v2 . For a screened short-range
interaction (T /v2 1/λ), the Pauli principle then leads to a suppression of the energy relaxation
rate by an additional factor of (T k/v2 )4 1 [27].
3. Spin reconstruction
The edge reconstruction in quantum Hall systems results from competition between the
Coulomb interaction and the confinement potential. Spin reconstruction at ν = 2 takes place
when the confinement potential Vc varies sufficiently slowly so that Vc0 < e2 /κl 2B and can be
understood at the level of the Hartree–Fock approximation [29–32]. Once the slope of the
confinement potential becomes weaker than that of the repulsive Hartree potential VH , it is
favorable to deposit charges outside the edge. This can be done without paying exchange energy
by a relative shift of the Fermi momenta of spin-up and spin-down particles, as depicted in
figure 2. In the absence of Zeeman splitting, εZ = 0, this is a second-order phase transition with
spontaneous breaking of the spin symmetry. Then, the distance between the two Fermi momenta
varies as kF2 − kF3 ∝ (|VH0 | − Vc0 )1/2 , eventually saturating at ∼ 1/l B [29]. For finite Zeeman
splitting εZ , the spin symmetry is lifted by the Zeeman field and the transition is smeared on the
scale of kF2 − kF3 ∼ εZ /v2 .
Spin reconstruction leads to characteristic changes in the single-particle dispersion that
develops an ‘eye structure’ (see figure 2(a)). Important for the relaxation dynamics is the
increase of v2 − v3 = (kF2 − kF3 )/m c , which enhances the typical energy transferred per step
of relaxation (see equation (4)).
For truly long-range interactions, the Hartree–Fock approximation predicts a logarithmic
singularity ∼ e2 /j ln(|k − kF |l B ) of the particle velocity at the Fermi energy, which is however
cut off in the presence of screening, say by a nearby gate electrode. The Fermi velocity is thus
still of the order of v2 , v3 ∼ e2 /κ for typical choices of the screening length.
Even with spin reconstruction, the relaxation of hot particles can be described within the
model dispersion of equation (3). We consider the case when the hot particle (not shown in
figure 2) is injected well outside the energy window e2 /(κl B ) of the reconstructed region. This is
compatible with the condition for validity of a perturbative expansion, which reduces to v1 v2
for the case of the Fermi velocity determined by the interaction.
The energy relaxation rate 1/τE(s) can now be derived in the same way as for the
unreconstructed edge, and consequently, equation (10) applies also to spin reconstructed edges.
The crucial difference is that the velocity difference v2 − v3 is now strongly enhanced by
the spin reconstruction, taking values up to v2 − v3 ∼ 1/(m cl B ). Comparing the rates before
[v2 − v3 ∼ max{εZ , T }/(m c v2 )] and well after spin reconstruction, we find an enhancement of
the relaxation rates given by
2
e /(κl B )
2
1 1
(s)
∼ . (13)
τE max{εZ , T } τE(u)
4. Charge reconstruction
For confinement potentials that vary even more smoothly, changing by e2 /κl B over a region
w > l B , charge reconstruction may occur such that some of the electrons at the edge are pushed
away from the bulk by a length of the order of l B [31, 32]. It leads to nonmonotonic behavior
of the dispersion with momentum and the creation of two additional counter-propagating5 edge
modes, as depicted in figure 3.
A minimal model for charge reconstruction considers filling factor ν = 1 within the
Hartree–Fock approximation [31]. It is convenient to formally model the confinement potential
by a positive background charge which is distributed spatially as if it was occupying the lowest
Landau-level wave functions ψ X with occupation numbers νc (X ) = 2(−X ). The advantage
of this model is that such a confinement potential exactly cancels the Hartree potential
of the electrons for an unreconstructed edge. In this case, the electron occupation of the
unreconstructed edge is stabilized by the (attractive) exchange potential.
The reconstruction transition can then be modeled by changing the abrupt drop of νc (X )
into a linear decrease over a length w. For the unreconstructed electron occupations, this leads
to negative (at X < 0) and positive (at X > 0) excess charges, causing a dipole field that favors
separating electrons from the bulk. Once this dipole field overcomes the exchange potential,
a charge reconstruction transition takes place. Within the Hartree–Fock approximation, this
happens for w ∼ 8l B . Due to the particle–hole symmetric choice of the confinement potential
around X = 0, the width and distance of the additional stripe from the bulk electron droplet
both take the same value b. Moreover, the transition is of first order in the sense that b changes
abruptly at the transition from zero to a value of the order of l B .
Note that the same mechanism induces new (weaker) effective dipole fields at each of the
three Fermi points as the edge becomes yet smoother. Thus, increasing w even further causes
additional stripes to appear, eventually approaching the limit of a compressible edge which is
expected for w l B [23]. In the following we will focus on w & l B , remaining well outside the
compressible limit.
Energy relaxation in the charge reconstructed case can also be captured by the dispersion
(3) when setting v3 < 0 and choosing the particle i = 2 to lie in one of the co-propagating
branches6 . The three Fermi velocities of the charge reconstructed edge are essentially
determined by the variation of the exchange potential, which is short range such that b & l B
already approximates the bulk edge (b → ∞) behavior. Consequently, the magnitudes of the
5
Note that the conductance may stay at the ν = 1 level due to impurity-induced elastic backscattering of electrons
at the Fermi level. At higher energies, inelastic scattering leads to a shorter mean free path and we can therefore
omit the disorder-induced backscattering.
6
Processes where the particle i = 2 is on the inner or outer co-propagating branch contribute at the same order
in the power counting argument. The choice depicted in figure 3 is slightly favored because of a larger Coulomb
matrix element due to the shorter distance to the hot particle.
5. Conclusions
Acknowledgments
We acknowledge financial support from the Deutsche Forschungsgemeinschaft under SPP 1538
(FUB) and by DOE under contract no. DEFG02-08ER46482 (Yale) and by the US–Israel BSF
under grant no. 2010366 (Yale). AL acknowledges support from Michigan State University.
Within this section we provide all essential details needed for the derivation of equation (7)
presented in the main text. We assume that the edge is smooth enough that we can approximate
the electron wave functions by those of the bulk. We start from the interaction matrix element
in real space
e2 /κ (1) (2)
(1) (2)
VdX (X, X ) = ψ X +dX , ψ X 0 −dX (1)
0 ψ , ψX 0 . (A.1)
|r − r(2) | X
In the following, we will measure all length scales in units of magnetic length l B . In these units
the guiding center coordinate directly translates to momenta. With the lowest Landau-level wave
functions of equation (1), we then find that
√
2 Z − 1x 2 +1y 2 /λ
e e
e−i1y dX e− 2 (x−X ) e− 2 (x−X −dX )
1 2 1 2
VdX (X, X 0 ) = dx dy d1x d1y p
πκ L 2
1x 2 + 1y 2
× e− 2 (x+1x−X ) e− 2 (x+1x−X +dX ) ,
1 0 2 1 0 2
(A.2)
√
where we used the screened Coulomb potential which carries an extra factor e− 1x +1y /λ
2 2
with λ being the distance to a screening gate. The integration over 1y gives 2K 0 (|1x dX |)
in the case when dX 1/k, where K 0 is the Bessel function of imaginary argument. If,
however, dX 1/λ the integral is cut off and the result changes to 2K 0 (|1x/λ|). We will derive
New Journal of Physics 14 (2012) 105009 (http://www.njp.org/)
12
results for the dX 1/λ case and keep in mind appropriate changes for the other limit. After
y integration, which gives a factor of L, we obtain the intermediate step
2e2 − 1 dX 2
Z
dx d1x K 0 (|1x dX |)e− 2 (2x−X −X +1x) e− 2 (X −X +dX +1x) .
1 0 2 1 0 2
VdX (X, X ) =
0
e 2 (A.3)
πκ L
√
Performing now the Gaussian integral over x, which gives a factor of π/2, followed by using
the Landau gauge to replace guiding center coordinates by momenta, one arrives at
r
2 e2 −q 2 /2
Z
dξ e− 2 (k+q+ξ ) K 0 (|ξ q|),
1 2
Vq (k) = e (A.4)
π κ
where we used the short-hand notation X − X 0 = kl 2B = k (with l B = 1). Note that Vq (k) =
Vq (−k − 2q) and is therefore not symmetric, which plays an important role. We see immediately
that scattering processes with large momentum transfer q 1 are exponentially suppressed. We
therefore concentrate on the opposite limit of q 1 when the exponential prefactor e−q /2 can
2
be set to unity.
Let us study limiting cases of equation (A.4). In the√
case when k 1, one can approximate
the exponential under the integral by the delta-function 2πδ(k + q + ξ ) and thus obtain
2e2
Vq (k) = K 0 (|q(k + q)|). (A.5)
κ
Using the asymptotic form of the Bessel function and restoring units of l B , one recovers the
second limit in equation (7).
In the other limiting case when k 1, one can approximate the exponential under the
integral of equation (A.4) by e−ξ /2 and then complete integration exactly with the result
2
e2
Vq (k) = K 0 (q 2 /4). (A.6)
κ
With the logarithmic accuracy at small q this translates into the first limit of equation (7).
In general, the three-particle scattering amplitude h10 20 30 |V G 0 V |123ic contains six terms: one
direct and five exchange contributions [28]. As explained in the text, we need only the former
one, which reads explicitly [28]
δ6,6 0 Vk30 −k3 (k3 − k2 )Vk10 −k1 (k20 − k10 ) Vk10 −k1 (k1 − k3 )Vk20 −k2 (k30 − k20 )
123
T10 20 30 = 2 +
L εk3 + εk2 − εk30 − εk2 +k3 −k30 εk1 + εk3 − εk10 − εk3 +k1 −k10
Vk20 −k2 (k2 − k1 )Vk30 −k3 (k10 − k30 ) Vk20 −k2 (k2 − k3 )Vk10 −k1 (k30 − k10 )
+ +
εk2 + εk1 − εk20 − εk1 +k2 −k20 εk2 + εk3 − εk20 − εk3 +k2 −k20
Vk10 −k1 (k1 − k2 )Vk30 −k3 (k20 − k30 ) Vk30 −k3 (k3 − k1 )Vk20 −k2 (k10 − k20 )
+ + , (B.1)
εk1 + εk2 − εk10 − εk2 +k1 −k10 εk3 + εk1 − εk30 − εk1 +k3 −k30
where the spin structure is δ6,6 0 = δσ1 ,σ10 δσ2 ,σ20 δσ3 ,σ30 and the Coulomb matrix element Vq (k) was
derived in appendix A. Now using the dispersion relation from equation (3) and the constraint
References
[1] Deshpande V V, Bockrath M, Glazman L I and Yacoby A 2010 Nature 464 209
[2] Imambekov A, Schmidt T L and Glazman L I 2011 arXiv:1110.1374
Imambekov A, Schmidt T L and Glazman L I 2011 Rev. Mod. Phys. submitted
[3] Barak G, Steinberg H, Pfeiffer L N, West K W, Glazman L, von Oppen F and Yacoby A 2010 Nature Phys.
6 489
[4] Chen Y-F, Dirks T, Al-Zoubi G, Birge N O and Mason N 2009 Phys. Rev. Lett. 102 036804
[5] Kinoshita T, Wenger T and Weiss D S 2006 Nature 440 900
[6] Hofferberth S, Lesanovsky I, Fischer B, Schumm T and Schmiedmayer J 2007 Nature 449 324
[7] Altimiras C, le Sueur H, Gennser U, Cavanna A, Mailly D and Pierre F 2010 Nature Phys. 6 34
[8] le Sueur H, Altimiras C, Gennser U, Cavanna A, Mailly D and Pierre F 2010 Phys. Rev. Lett. 105 056803
[9] Altimiras C, le Sueur H, Gennser U, Cavanna A, Mailly D and Pierre F 2010 Phys. Rev. Lett. 105 226804
[10] Paradiso N, Heun S, Roddaro S, Sorba L, Beltram F and Biasiol G 2011 Phys. Rev. B 84 235318
[11] Granger G, Eisenstein J P and Reno J L 2009 Phys. Rev. Lett. 102 086803
[12] Wen X G 1990 Phys. Rev. Lett. 64 2206
Wen X G 1991 Phys. Rev. B 43 11025
[13] Kane C L and Fisher M P A 1994 Phys. Rev. B 51 13449
Kane C L and Fisher M P A 1995 Phys. Rev. B 52 17393
Kane C L and Fisher M P A 1997 Phys. Rev. B 55 15832