Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Учебники / Физика конденсированных сред

.pdf
Скачиваний:
50
Добавлен:
30.03.2022
Размер:
3.97 Mб
Скачать

§ 7. Тензор электропроводности в квантующем поле 311

Запишем компоненты тока Jμ и Jν , используя определение координат центра ларморовской орбиты и координат относительного движения, которые далее будут пониматься как суммарные величины для всей системы электронов

˙ ˙ ˙

Jx = e(ξ + X), Jy = e(η˙ + Y ).

Для упрощения обозначений и сокращения объема формул вве-

дем так называемое скалярное произведение двух операторов

 

 

Кубо A; B(t)

(см. также формулу (5.82)), которое в кван-

тующем магнитном поле является четной функцией времени, если операторы A и B совпадают (доказательство этого замечательного факта в § 9 настоящей главы):

 

 

β

 

 

 

 

(5.130)

A, B(t)

 

= Sp{AB(t + i λ)ρ0}.

 

 

0

 

 

Тогда, используя выражение (5.129), для σxx(0) на нулевой

частоте получаем

 

 

 

 

 

 

 

0

0

 

σxx(0) = dt1e− t1

Jx(t1), Jx) = dt1e− t1

Jx, Jx(t1)) =

 

 

X˙

 

0

= e2 dt1e− t1

ξ,˙

ξ˙(t1) + ξ,˙

(t1) + X˙ ,

˙ ˙

+ X, X(t1) .

ξ˙(t1) +

(5.131)

Покажем, что все члены, кроме последнего, равны нулю. Для доказательства рассмотрим корреляционную функцию

 

 

 

 

 

 

dt

 

 

 

 

 

 

0

 

0

 

 

 

 

dte− t

ξ,˙ ξ˙(t) =

dte− t

d

ξ,˙

ξ(t) =

 

 

 

=

ξ,˙

 

t→∞

 

 

 

 

 

 

ξ,˙

 

 

ξ(0) + lim e− t

ξ,˙ ξ(t) +

dte− t

ξ(t) . (5.132)

0

312

Глава 5. Теория линейного отклика

 

При получении выражения (5.132) мы выполнили интегри-

рование по частям. Пользуясь принципом ослабления корреляций, согласно которому корреляция между двумя физическими

величинами, взятыми в момент времени t1

и t2 , ослабевает с

увеличением временного интервала t = t1 − t2 , получаем

t→∞

 

t

˙

t

{

˙

0}

{

0} →

 

(5.133)

e

 

ξ, ξ(t) → e

 

 

0.

lim

 

 

 

 

Sp ξρ

Sp ξ(t)ρ

 

Для преобразования последнего выражения в формуле (5.132) воспользуемся теоремой Абеля, согласно которой в термодинамическом пределе справедливо равенство

 

 

 

 

 

 

 

lim dte− t

ξ,˙

ξ(t) =

 

 

+0

0

 

 

 

 

lim

˙

˙

} Sp(t)ρ0} = 0.

(5.134)

= t→∞ ξ, ξ(t) = Sp{ξρ0

Обратимся, наконец, к преобразованию первого слагаемого во второй строчке формулы (5.132). Сделав замену переменной интегрирования τ = λ/β в формуле скалярного произведения (5.130), получим

 

 

0

{

}

 

 

 

0

{

}

 

ξ,˙

ξ =

 

β

1 Sp ξρ˙ 0τ ξρ01−τ

 

=

 

β

1 Sp ξρ01−τ ξρ˙ 0τ

 

=

 

 

 

 

 

 

= −β 1 Sp{ξρ0τ ξρ˙ 01−τ } =

 

 

 

 

 

 

1

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

Sp[ρ0, ξ]} =

 

Sp{[ξ, ξ]ρ0} = 0.

 

(5.135)

 

 

i

i

 

При получении результата (5.135) мы сделали замену переменных τ = 1 − τ и воспользовались формулой Кубо (5.60).

Совершенно аналогично можно доказать, что

 

 

 

 

i

{

}

 

eH

 

0

 

 

 

dte− t

ξ,˙

η˙(t) =

 

ξ,˙

η =

1

Sp [ξ, η]ρ0

 

=

nc

;

(5.136)

 

 

 

 

 

 

0

dte− t

X,˙

ξ˙(t) = 0;

 

 

0

dte− t

Y˙

, η˙(t) = 0. (5.137)

§ 8. Электропроводность при рассеянии на фононах 313

Равенство нулю корреляционных функций в формуле (5.137) следует из того факта, что коммутаторы [X, ξ] и [Y, η] равны нулю.

Теперь можно вернуться к формуле (5.131) и записать результат для σxy (0) и других компонент электропроводности в квантующем магнитном поле. Учитывая результаты (5.134) – (5.137), получаем

 

H

 

 

 

 

 

 

0

 

σxy =

enc

+ e2

dt1e− t1 X˙ (t1), Y˙

;

 

 

 

 

 

 

 

 

 

 

0

 

σxx = e2 dt1e− t1 X˙ , X˙

(t1) ;

 

 

 

0

 

 

 

 

σyy = e2 dt1e− t1 Y˙ , Y˙

(t1) .

(5.138)

(5.139)

(5.140)

Анализируя полученные результаты, нетрудно заметить, что в квантующем магнитном поле диагональные компоненты σxx и σyy отличны от нуля только благодаря процессам рассеяния, поскольку, как следует из уравнений движения для этих величин (5.127), σxx и σyy пропорциональны по меньшей мере квадрату константы взаимодействия электронов с рассеивателями. Недиагональная компонента σxy содержит не зависящий от процессов рассеяния бесстолкновительный вклад enc/H и квадратичную по константе взаимодействия электронов с рассеивателями поправку. Важно отметить, что компоненты тензора электропроводности в квантующем магнитном поле выражаются через корреляционные функции координат центров ларморовских орбит, которые в борновском приближении теории рассеяния могут быть непосредственно вычислены.

§8. Вычисление электропроводности в случае квазиупругого рассеяния на фононах

Рассмотрим вычисление σxx и σyy в случае квазиупругого

ˆ

рассеяния на фононах, полагая, что оператором U в формулах

314

Глава 5. Теория линейного отклика

(5.127) является гамильтониан электрон-фононного взаимодействия (4.76).

Чтобы лучше понимать полученные выше результаты для компонент тензора электропроводности в квантующем магнитном поле, полезно сравнить их с результатами, которые дает метод кинетического уравнения в пределе сильного ω0τp 1 (но неквантующего) магнитного поля. На основании формул (4.118), (4.121) и (4.128) получаем

σxy =

e2n

ω0τp2

 

 

enc

;

(5.141)

m

 

1 + (ω0τp)2

H

 

e2n

τ

 

 

e2n

 

 

 

σxx =

 

 

 

 

 

p

 

 

 

 

.

(5.142)

m

 

1 + (ω0τp)2

02τp

Таким образом, в пределе сильного магнитного поля, если при разложении знаменателя в формуле (5.141) по малому параметру 1/(ω0τp)2 оставить только нулевой член разложения, кинетическое уравнение для недиагональной компоненты электропроводности дает тот же самый бесстолкновительный вклад, что и формула (5.138).

Выражение для диагональной компоненты (5.142) позволяет, по крайней мере формально, определить время релаксации импульса в квантующем магнитном поле. Действительно, сравнивая два выражения (5.142) и (5.139), получаем определение для времени релаксации импульса в квантующем магнитном поле:

1

τp

2

 

 

 

 

n

0

 

=

0

dt1e− t1

X˙ , X˙

(t1) .

(5.143)

 

Проблема вычисления компонент σxx и σxy тензора электропроводности в борновском приближении теории рассеяния по существу сводится к квадратурам, поскольку в этом приближении взаимодействием в статистическом операторе и операторе эволюции можно пренебречь и тогда эти операторы имеют только диагональные матричные элементы на классе собствен-

| ˆ

ных функций ν > оператора H0 . Поскольку корреляционная

§ 8. Электропроводность при рассеянии на фононах 315

функция A, B(t) является четной функцией временного аргумента t , интегрирование по переменной t1 в интеграле(5.143) можно распространить до −∞, сделав пределы симметричными. Тогда интегрирование по времени t1 даст δ -функцию.

Далее квантово-статистическое среднее по электронным и фононным переменным представим в виде

 

)

< ν | X˙

| ν >< ν | X˙

*

 

Sp{X˙ X˙ (t)ρ0} =

 

(t) | ν > s ×

 

νν

 

× fν (1 − fν ) ei/ (εν −εν ).

 

 

 

(5.144)

В формуле (5.144) угловые скобки, помеченные индексом s, обозначают квантово-статистическое усреднение по состояниям рассеивателей. При выводе этой формулы мы воспользовались также статистической теоремой Вика – Блоха – Доминисиса (5.75).

Наконец, можно показать (предлагаем это доказательство провести самостоятельно), что

 

 

β

 

 

 

 

 

 

dte− |t|

Sp{X˙ X˙ (t + i λ)ρ0} = β dte− |t| Sp{X˙ X˙ (t)ρ0}.

−∞

 

 

0

 

 

 

−∞

 

 

(5.145)

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая сказанное, выражение для компоненты σxx мож-

но представить в такой форме:

 

 

 

 

 

 

e

2

 

 

 

 

 

 

 

 

σxx =

 

 

 

d t e− |t|

 

dE f (E) δ(E

Hˆ0)X˙ (1

f (Hˆ0)X˙ (t) .

 

 

 

 

2kБT −∞

 

−∞

)

 

 

*

 

 

 

 

 

 

 

 

 

 

 

 

(5.146)

В этой формуле большие угловые скобки . . .

 

обозначают

) *

квантово-статистическое среднее по фононным переменным и квантово-механическое среднее по одночастичным электронным состояниям.

˙ ˙

Запишем в явном виде операторы X и X(t) . Учитывая, что операторы координаты центра ларморовской орбиты коммути-

316

 

 

Глава 5. Теория линейного отклика

 

 

 

 

руют с гамильтонианом

ˆ

 

, получаем

 

 

 

 

 

 

H0

 

 

 

 

 

 

 

 

 

il2

q

qy Cq bq eiq r − Cq bq+ e−iq r!;

 

 

 

 

X˙

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

il2

ei/ Hˆ0t q

 

 

 

 

 

 

 

 

 

 

X˙

(t) =

 

qy Cq bq eiq r e−iΩq t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

(5.147)

 

 

 

 

 

 

−Cq bq+ e−iq reiΩq t! e−i/ H0t.

 

Подставим эти выражения в выражение (5.146):

 

 

 

 

 

 

e2π

 

 

 

 

 

 

l4

 

 

 

(Nq + 1) ×

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σxx = kБT

−∞

 

 

 

 

 

2 qy2 | Cq |2

 

 

 

 

dE f (E)

 

q

 

 

× Sp δ(E − Hˆ0)eiq r

1 − f (Hˆ0) δ(E − Hˆ0 Ωq )e−iq r

+

 

 

+N Sp δ(E Hˆ )e

 

 

 

 

1 f (Hˆ ) δ(E Hˆ + Ω )e

iq r

! .

 

 

q

 

 

 

0

 

iq r

 

 

 

 

0

 

0

q

!

(5.148)

 

 

 

 

 

 

 

 

 

Шпур в формуле (5.148) означает суммирование по полному набору квантовых чисел ν = {n, px, pz , σ}, характеризующих состояние электрона в квантующем магнитном поле.

Для дальнейшего преобразования выражения (5.148) учтем два тождества

Nq [f (E − Ωq ) − f (E)] = f (E)[1 − f (E − Ωq )];

[Nq + 1][f (E + Ωq ) − f (E)] = f (E)[1 − f (E + Ωq )], (5.149)

которые проверяются непосредственной подстановкой функций распределения. С учетом этих соотношений выражение (5.148) можно переписать в более удобной для дальнейших преобразований форме

 

e22π

 

 

f (E

 

Ωq )

 

f (E) Ωq

 

σxx =

 

dE

l4qy2 |

Cq |2

 

Ωq

 

 

×

 

 

 

 

 

kБT

 

−∞

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×Nq [Nq + 1] Sp

δ(E − Hˆ0)eiq rδ(E − Hˆ0 Ωq )e−iq r!.

(5.150)

При выводе этой формулы при преобразовании второго слагаемого выражения (5.148) сделана замена переменных E + Ωq → E .

§ 8. Электропроводность при рассеянии на фононах 317

Полученное выражение справедливо в случае неупругого рассеяния и квазиупругого рассеяния на фононах. В качестве примера рассмотрим упругое рассеяние на фононах. Поскольку нас интересует только принципиальная сторона вопроса, связанная с методикой вычисления кинетических коэффициентов в магнитном поле, ограничимся наиболее простым случаем ультраквантового предела, когда в переносе заряда участвуют только электроны самой нижней подзоны Ландау с номером n = 0. В этом случае матричные элементы от операторов экспонент в выражении (5.148) легко вычисляются, и мы имеем

|< 0, pz + qz , px + qx | eiq r | 0, pz , px >|2= el2q2 /2; q = qx2 + qy2.

При вычислении шпура по электронным переменным суммы удобно заменить интегралами

 

 

LxLz

 

 

2V

 

 

 

 

 

−∞

 

 

2

(2π )2

d pz d px

(2πl)2

 

d pz .

(5.151)

ν σ

 

 

 

 

 

 

 

 

Для того чтобы получить последний результат, необходимо учесть кратность вырождения электронных состояний по квантовому числу px . Для подсчета этого числа наложим условие цикличности на волновую функцию электрона (5.120) по осям X и Z , т. е. потребуем, чтобы координатам x + Lx , z + Lz и x , z соответствовала одна и та же функция. Если учесть реальный вид волновой функции (5.120), то это требование приводит к условию

p

 

=

2π

n , p

 

=

2π

n ,

x

 

z

 

 

 

 

x

 

 

z

 

 

 

Lx

 

 

Lz

где nx и nz – некоторые целые числа. По оси Y не будем налагать условие цикличности, но потребуем, чтобы решение (5.120) существовало только тогда, когда координата центра ларморовской орбиты y0 находится в области

0 <| y0 |< Ly ,

(5.152)

318 Глава 5. Теория линейного отклика

где Ly – размер образца по оси Y (легко проверить, что y0 – это одна из координат центра ларморовской орбиты и y0 = Y ). Та-

ким образом, максимальное значение координаты центра ларморовской орбиты | y0 |max= Ly . Так как | y0 |= c/(eH)px , нахо-

дим максимальное значение квантового числа px : px = /l2 Ly . Тогда интеграл по px в формуле (5.151) равен /l2 Ly и мы получаем последний результат в этой формуле.

Суммирование по волновому вектору фононов также следует заменить интегрированием в цилиндрической системе координат

 

 

 

 

 

 

 

 

 

V

dq

V

q

 

d q

 

d qz .

(5.153)

 

 

 

(2π)3

(2π)2

 

 

 

 

 

 

 

 

0

 

 

−∞

 

 

q

 

 

 

 

 

 

 

 

 

Учитывая сделанные замечания, выражение для статиче-

ской проводимости σxx

 

в случае квазиупругого рассеяния за-

пишем следующим образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

l

4

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂f (E)

 

 

σxx =

2πe

 

 

 

 

 

q d q

d qz

dE qy2

 

Cq

 

2

 

 

 

(2π)

2

|

|

 

∂E ×

 

 

 

0

−∞ −∞

 

 

 

 

 

kБT

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

d pz δ(ε0pz

 

E)δ(ε0pz+ qz

 

E)e−l

q /2.

(5.154)

 

 

 

 

 

 

 

 

2

 

 

 

 

 

× Ωq (2πl)

−∞

 

 

 

 

 

 

 

 

 

 

 

В целях дальнейшего преобразования выражения (5.154)

рассмотрим интеграл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I =

 

d qz

 

 

d pz qt δ(ε0pz E) δ(ε0pz+ qz

− E).

(5.155)

 

 

 

−∞

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из закона сохранения энергии следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(p

z

+ q )2

(p

)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

=

z

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2m

 

 

2m

 

 

 

 

 

 

 

Поэтому участвующие в рассеянии фононы имеют квазиимпульс qz pz , что позволяет легко оценить продольную со-

ставляющую волнового вектора фононов, участвующих в рас-

сеянии,

 

 

qz

2mkБT

λ1 107см1

 

 

§ 8. Электропроводность при рассеянии на фононах 319

для температур, при которых реально проводится эксперимент. Перпендикулярная составляющая q волнового вектора фононов, участвующих в рассеянии, лимитируется обрезающим

фактором

e−l2q2 /2,

поэтому можно считать, что q 1/l и перпендикулярная составляющая волнового вектора фононов, участвующих в рассеянии, по порядку величины совпадает с обратной магнитной длиной. Поскольку в квантующем поле магнитная длина l меньше длины волны электрона λ , будем считать, что выполняется условие λ l . В силу сделанных выше оценок, следует, что q qz и q q . Таким образом, интегралы, входящие в выражение I (5.155), могут быть достаточно просто вычислены:

 

 

 

 

 

 

2

 

 

 

 

 

ω0

 

 

 

 

 

I =

d qz

 

 

d pz qt δ

pz

 

+

 

E

 

 

 

 

 

 

 

 

 

−∞

 

 

−∞

 

 

 

2m

 

 

 

 

2

 

 

×

 

 

 

p

q

 

2q2

 

 

t

2m

1

 

 

 

 

 

z

 

z

 

z

 

 

 

 

 

 

 

 

 

 

 

 

δ

 

 

 

+

 

= q

 

 

 

 

 

.

(5.156)

 

m

2m

 

 

E − ω0

При выполнении интегрирования дельта-функций полезно

использовать известную формулу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

(5.157)

δ ϕ(x) =

i

δ(x − xi) dx x=xi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где xi – это корни уравнения ϕ(x) = 0 .

Подставляя результат (5.156) в формулу (5.154), получаем следующее выражение для диагональной статической электропроводности

2

2

kБT mC

 

 

 

 

∂f

1

 

 

2

σxx =

e

l

 

 

dE

 

 

dq2 qt+1e−l2q /2.

 

 

 

 

 

 

 

 

 

 

 

 

4π3 4s

−∞

 

 

∂E

E − ω0/2 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.158)

При получении этого результата мы использовали обозначе-

ние

 

 

 

 

 

 

 

 

 

 

E02

 

 

 

 

 

 

 

 

2

 

 

t

 

 

 

 

 

 

 

 

Cq

 

= C q , C =

2ρs

.

 

 

320

Глава 5. Теория линейного отклика

Для упрощения дальнейших вычислений рассмотрим только случай сильновырожденного электронного газа. Тогда интеграл по энергии вычисляется элементарно, если воспользоваться аппроксимацией

∂E∂f δ(E − ζ).

Интеграл по q , очевидно, сводится к гамма-функции. Поэтому дальнейшее интегрирование не представляет проблем, и мы сразу получаем

 

e2l2kБT mC 1

2

 

2

 

 

 

 

 

 

t+3

σxx =

 

 

 

 

 

 

4π3 4s

ζ − ω0/2

l2

 

t + 3

Γ . (5.159)

2

Характерной особенностью полученного результата является наличие расходимости в том случае, когда уровень Ферми пересекает подуровень Ландау. Эта особенность возникает изза того, что в квантующем магнитном поле плотность состояний электронов в пространстве энергий имеет корневую особенность на дне каждой подзоны Ландау. Особенно интересные эффекты в квантующем магнитном поле возникают в двумерном металле, который легко реализуется в полевом транзисторе. Нобелевская премия по физике дважды присуждалась за исследование квантового эффекта Холла: в 1985 г. за открытие этого явления и в 1998 г. за открытие и интерпретацию дробного квантового эффекта Холла. Не имея возможности подробнее остановиться на этой интересной теме, отсылаем читателя к специальной литературе [41].

Задача 5.1

Получить выражение для плотности состояний в пространстве энергий для электронов проводимости в квантующем магнитном поле.

Решение

Ввести понятие плотности состояний в пространстве энергий проще всего, используя соотношение

 

 

 

 

 

→ g(E)dE,

(5.160)

 

n px pz σ

0