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

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

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

§ 1. Решение в приближении времени релаксации 191

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

 

 

˙

∂f0

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p pf =

∂εp

 

evE + pf1

 

c [v

× H].

 

(4.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

и по-

При получении результата (4.12) мы учли, что pf0

этому вклад

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

× H] = 0,

 

 

 

 

 

 

 

 

 

 

pf0 c [v

 

 

 

 

 

 

поскольку

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v [v × H] = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя результаты (4.7), (4.10), (4.12) в кинетическое

уравнение (4.6), получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

∂f0

 

εp − ζ

 

 

f1

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

∂εp

 

= τp

 

 

 

 

(4.13)

v eε −

T

T

+ pf1 c

[v × H].

В том случае, когда магнитное поле равно нулю, выражение (4.13) сразу позволяет определить поправку к функции распределения f1 , линейную по градиенту электрохимического потенциала ε = (ϕ + 1/e ζ) (1.14) и градиенту температуры T :

f

 

= τ

p

∂f0

v eε

εp − ζ

T .

(4.14)

 

∂εp

 

 

1

 

 

T

 

Если магнитное поле не равно нулю, то для решения уравнения (4.13) будем искать поправку f1 в виде

f1 =

∂f0

v χ(εp) ,

(4.15)

∂εp

где χ(εp) – неизвестная векторная функция, зависящая только от энергии.

192 Глава 4. Кинетическое уравнение для электронов

Вычислим градиент в импульсном пространстве от функции f1 . Пользуясь определением (4.15), получаем

 

 

 

 

2f0

 

 

 

 

 

 

 

∂f0

 

 

 

 

 

 

pf1 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂εp2

 

 

 

 

∂εp

 

 

 

 

 

 

 

v v χ(εp)

 

 

 

p v χ(εp) ,

 

 

p

v χ(ε

) = v

 

p

χ (ε

) + χ (ε

)

p

v ,

 

 

 

 

 

∂χi(εp)

 

 

 

 

 

 

 

 

 

χ(εp)

(4.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p χi(εp) =

∂εp

 

 

v,

χi(εp) p vi =

m .

Подставляя полученные результаты в последнее слагаемое в правой части (4.13) и учитывая, что члены, пропорциональные вектору скорости v , вклада не дадут, получаем простое выражение

 

e

 

∂f0

 

 

 

 

 

 

 

 

 

 

[v × H] = ∂εp

ω0 v [h × χ(εp)].

(4.17)

pf1 c

При выводе формулы (4.17) мы воспользовались определением частоты ларморовской прецессии ω0 электронов в магнитном поле

ω0

=

eH

,

(4.18)

mc

 

 

 

 

ввели единичный вектор h , ориентированный вдоль направ-

ления вектора индукции магнитного поля (H = h H) , H – модуль вектора индукции магнитного поля, и переставили порядок векторов в векторно-скалярном произведении.

Подставляя результаты (4.15), (4.17) в формулу (4.13) и производя необходимые сокращения, получаем векторное уравне-

ние для определения функции χ(εp)

 

 

 

 

χ(εp) + [a ×

 

 

 

 

(4.19)

χ(εp)] = b;

 

 

 

 

eε −

εp − ζ

 

.

(4.20)

a = ω0τp h;

b = τp

T

T

Для решения уравнения (4.19) умножим его один раз скалярно, а второй раз векторно на вектор a . В результате простых

§ 1. Решение в приближении времени релаксации 193

алгебраических преобразований можно выразить вектор χ(εp)

через векторы a и b :

 

 

 

 

 

 

χ(εp) =

b + a (a b ) [a × b ]

.

(4.21)

 

 

1 + a2

 

 

 

Для получения более удобной структуры решения воспользуемся тождеством

b =

1

a (a b ) [a × [a × b ] ]!

(4.22)

a2

и подставим это выражение для b

тате получаем представление для

через векторы a и b :

в решение (4.21). В резульфункции χ(εp) , выраженое

1

χ(εp) = a2

 

 

 

2

[a

 

 

 

a (a b )

[a × b] + 1/a

 

× [a × b ] ]

.

(4.23)

 

 

1 + a2

 

 

 

 

Наконец, учитывая явный вид векторов a и b и подставляя их значения (4.20) в формулу (4.23), получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω0τp [h

× ε ] + [h ×

[h × ε ] ]

 

 

 

 

 

χ(εp) = p (h ε )h −

 

 

 

 

 

+

 

 

 

 

1 + (ω0τp)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω0

τp [h

 

 

 

 

 

+τ

ζ − εp

(h

T )h

 

×

T ] + [h × [h × T ] ]

.

(4.24)

 

 

p

T

 

 

 

 

1 + (ω0τp)2

 

 

 

Формулы (4.14), (4.15), (4.24) будут использованы в дальнейшем для определения потоков заряда и тепла и вычисления кинетических коэффициентов, определяющих термомагнитные и гальваномагнитные явления в проводящих кристаллах.

194 Глава 4. Кинетическое уравнение для электронов

§ 2. Условия применимости квазиклассического описания электронов в проводящих кристаллах

Записанное в предыдущем разделе кинетическое уравнение (4.6), (4.13) является квазиклассическим. Поскольку хорошо известно, что электроны в кристалле – это квантовые объекты и есть достаточно много убедительных эффектов (например дифракция электронов в кристаллах), в которых квантовые свойства электронов наглядно проявляются, встает вопрос о применимости такого описания. В действительности квазиклассическое описание накладывает некоторые ограничения на условия проведения физического эксперимента, но можно показать, что для большинства реальных ситуаций, в которых производится измерение кинетических явлений в твердых телах, применение квазиклассического описания вполне оправданно. Ниже сформулированы основные условия применимости квазиклассического кинетического уравнения для описания кинетических явлений в проводящих кристаллах при наличии постоянного внешнего магнитного поля и без него.

Эти условия приводят к трем основным ограничениям. Во-первых, длина волны электрона λ должна быть мень-

ше других характерных пространственных масштабов задачи, что позволяет рассматривать электрон как точечный объект. В отсутствие магнитного поля естественным параметром размерности длины является длина свободного пробега l . Поэтому квазиклассическое описание возможно, если

λ l.

Во-вторых, неопределенность в энергии электрона E , которая является следствием квантово-механического принципа неопределенности, должна быть малой по сравнению со средней энергией электрона ε (средняя энергия ε kБT для невырожденного случая и ε ζ в условиях вырождения)

ε, τ0

p,

§ 2. Условия применимости квазиклассического описания 195

где τ0 – характерное время взаимодействия электрона с другими подсистемами кристалла. Поэтому время характерного взаимодействия электронов с рассеивателями τ0 должно быть достаточно большим. В этом случае столкновительное уширение энергетических уровней можно считать пренебрежимо малым и температура будет единственным параметром, хаотизирующим движение носителей заряда. Это условие служит основанием для описания электронной системы на языке функции распределения. Вместе с тем характерное время столкновений должно быть существенно меньше, чем время между двумя последовательными столкновениями, поскольку каждое из них рассматривается как независимый процесс и считается, что после каждого столкновения в системе успевает сформироваться неравновесное распределение. Поэтому в качестве верхней оценки времени τ0 можно взять время между двумя последовательными столкновениями и считать, что τ0 ≈ τp . Это тем более оправданно, поскольку величина τp легко поддается экспериментальной оценке. В этом случае условие

≈ ε

τ0 τp

по существу сводится к условию λ l. В этом легко убедиться, если предыдущее условие записать в виде

v

ε

умножив левую и правую часть неравенства на среднюю скорость электронов v . Тогда в левой части неравенства стоит величина /p ≈ λ , а в правой – длина свободного пробега l = p . Таким образом, первое и второе ограничения в отсутствие магнитного поля приводят к одному и тому же условию λ l.

Еще одно ограничение возникает в том случае, когда электрон находится в области действия внешних силовых полей. В этом случае, если электрон рассматривается как точечный объект, изменение его энергии на длине порядка длины волны де Бройля должно быть много меньше средней энергии электрона. Иначе говоря, если, например, рассматривать движение электрона во внешнем электрическом поле, то должно выполняться условие λ eE ε kБT для невырожденного случая.

196 Глава 4. Кинетическое уравнение для электронов

Это ограничение не является очень существенным, поскольку простые оценки дают ограничение E < 106 В/м , что вполне приемлемо для большинства экспериментальных ситуаций.

Рассмотрим условия применимости кинетического уравнения для описания неравновесных носителей тока в магнитном поле. В этом случае имеется три характерных параметра размерности длины: lH , так называемая «магнитная длина»– характерный размер ларморовской орбиты электрона

 

c

1/2

lH =

 

 

,

eH

длина волны электрона де Бройля

 

λ = 2π

 

 

2

и длина свободного пробега электронов между двумя последовательными столкновениями l .

В этих условиях первым критерием квазиклассического описания является условие

2π c 1/2. 2mε eH

Полагая в этом выражении ε kБT , получаем хорошо известное в литературе условие применимости квазиклассического описания электронов в магнитном поле:

ω0 kБT.

(4.25)

Это условие допускает простую интерпретацию: при квазиклассическом описании расстояния между квантованными уровнями энергии электронов в магнитном поле должны быть малыми, по сравнению со средней энергией теплового движения электронов.

Другое условие применимости кинетического уравнения в магнитном поле также связано с влиянием магнитного поля на орбитальное движение электронов. Расстояние между уровнями Ландау в магнитном поле ω0 должно быть существенно

§ 2. Условия применимости квазиклассического описания 197

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

ω0τp 1.

(4.26)

Условие (4.26) может иметь и другую интерпретацию: для того чтобы квазиклассическое описание было применимо, необходимо, чтобы электрон, двигаясь по циклотронной орбите, между двумя актами рассеяния успевал пройти лишь малую часть периода T круговой траектории

T

2πlH

,

2π

= ω0,

lH

=

1

.

 

 

 

T

 

 

 

ω0

 

v

v

В этой формуле v – в вырожденном случае скорость электрона на поверхности Ферми (в невырожденном случае эту величину следует заменить на среднюю тепловую скорость). Поскольку время релаксации импульса

l τp v ,

условие ω0τp 1 может быть записано также в виде lH l . Иначе говоря, радиус циклотронной орбиты должен быть много больше длины свободного пробега электрона.

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

Если выполняются неравенства lH l или, что эквивалентно, ω0τp 1 , то магнитные поля называются слабыми.

Если выполняются обратные неравенства, то магнитные поля называются сильными. В этом случае магнитное поле существенно искривляет траекторию движения электронов, но если его влияние можно не учитывать при расчете вероятностей рассеяния, то и в случае сильных магнитных полей даже в условиях ω0τp 1 кинетическое уравнение применимо для описания кинетических явлений в магнитном поле. Естественно, что условие λ l должно оставаться справедливым.

При дальнейшем увеличении магнитного поля нарушается условие (4.25) и магнитное поле становится квантующим.

198 Глава 4. Кинетическое уравнение для электронов

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

§ 3. Определение потоков заряда и тепла. Вычисление кинетических коэффициентов в случае H = 0

Обобщая простейшее выражение для потока заряженных

частиц J = e n v , где n – число частиц, имеющих скоростьv , для плотности потока заряда и тепла получаем выражения

 

J =

 

e

 

 

dp f (p ) v,

(4.27)

 

 

 

 

(2π )

3

 

 

 

 

σ

 

 

 

 

 

 

1

 

 

dp (εp

 

 

 

JQ =

 

 

ζ) f (p ) v.

(4.28)

 

(2π )

3

σ

 

 

 

 

 

 

 

 

Суммирование по спиновому квантовому числу в формулах (4.27), (4.28) дает численный множитель, равный двум, так как мы не учитываем спиновое расщепление уровней. При записи выражения (4.28) мы учли определение (1.12).

Из физических соображений легко понять, что в формулах (4.27), (4.28) отличный от нуля вклад дает только неравновесная поправка к функции распределения f1(p ) , определяемая выражением (4.14) в отсутствие магнитного поля и выражениями (4.15), (4.24) при наличии внешнего магнитного поля.

Рассмотрим сначала кинетические явления в отсутствие внешнего магнитного поля. В этом случае кинетические коэффициенты являются скалярными величинами. Тогда, подставляя выражение (4.14) в формулы, определяющие поток заряда и тепла, имеем

 

2

e

 

(4.29)

 

 

 

 

 

 

 

 

J = e K0ε − T K1 T,

 

 

1

 

(4.30)

 

 

 

 

 

JQ = eK1ε − T K2 T,

§ 3. Вычисление кинетических коэффициентов

199

где интегралы Kl, l = 0, 1, 2 , определены соотношением

Kl =

2 1

 

dp

∂f0

τp v2 (εp − ζ)l.

(4.31)

(2π )3

 

3

∂εp

При выводе формул (4.29) – (4.31) мы учли, что для произвольной функции модуля квазиимпульса электронов Φ(εp) справедливо представление

dp Φ(εp) vi vj = 13 dp Φ(εp) v2δij ,

где δij – символ Кронекера, i, j = x, y z .

Сравнивая формулы (4.29), (4.30) с соответствующими феноменологическими результатами (1.15), (1.35), (1.36), получаем выражение кинетических коэффициентов, определяющих явления теплоэлектропроводности и термоэлектрические явления через введенные выше интегралы Kl :

1

 

 

K1

 

 

ρ =

 

 

, α =

 

 

 

,

e2 K

e T K

0

0

 

 

 

 

κ =

K2 K0 − K12

.

(4.32)

 

 

T K0

Таким образом, для вычисления этих кинетических коэффициентов необходимо вычислить интегралы Kl (4.31). Переходя к интегрированию по энергии в формуле (4.31), после выполнения интегрирования по полярному и азимутальному углам в сферической системе координат имеем

 

2 (2m)1/2

 

 

∂f

0

 

 

 

 

Kl =

 

 

 

p

 

τp εp3/2

(εp

ζ)l. (4.33)

2

3

∂εp

 

3π

 

0

 

 

 

 

Интегралы по энергии, содержащие функцию Ферми или ее

производные, могут быть сведены к так называемым интегра-

лам Ферми Fp

ζ/kБT

индекса p :

 

 

 

 

 

 

 

 

p

 

 

 

 

Fp

 

ζ

 

 

 

x

 

 

dx,

(4.34)

 

=

 

 

 

 

kБT

 

0 exp

x − ζ/kБT

 

+ 1

 

 

200 Глава 4. Кинетическое уравнение для электронов

зависящим от параметра ζ/kБT , для которых хорошо известны различные асимптотические представления [8]. Общие выражения, которые получаются при этом, оказываются достаточно громоздкими, и поэтому мы рассмотрим лишь два предельных случая, для которых легко получаются простые оценки.

Случай сильного вырождения

В этом случае ζ > 0, ζ/kБT 1 и производная по энергии от функции распределения имеет резкий максимум при εp = ζ . На рис. 26 представлена зависимость функции Ферми – Дирака

f0(x) =

1

, x =

εp − ζ

ex + 1

kБT

 

 

и ее первой производной от безразмерного параметра x .

Рис. 26. Графики функция распределения Ферми – Дирака и ее производной:

a – функция распределения Ферми – Дирака f0(x) в зависимости от аргумента x = (εp − ζ)/kБT ; б – производная функции распределения

∂f (x)

в зависимости от того же аргумента

∂x0

Как следует из рис. 26 б, производная от функции распределения отлична от нуля лишь в небольшом интервале энергий εp kБT . Эта особенность производной функции распределения широко используется для построения приближенных формул вычисления интегралов, содержащих в качестве подынтегральной функции произведение гладкой функции Φ(εp ) и производной от функции распределения Ферми – Дирака по