Электронные свойства твердых тел

advertisement
Кафедра твердотельной электроники и радиофизики МФТИ
Курс лекций для студентов IV курса ФФКЭ
"Электронные свойства твердых тел"
профессор С.Н.Артеменко
Текст лекций был набран в 2000 — 2001 гг. усилиями студентов ФФКЭ:
1.
2.
3.
4.
5.
6.
7.
8.
Деменцова Александра
Иванова Михаила
Чугурова Анатолия
Труфанова Андрея
Дронова Александра
Файзулина Романа
Коломейченко Игоря
Малышева Дмитрия
2007 г.
Глава 1
Основные методы и приближения для
описания электронных состояний в твердом
теле.
1.1
Адиабатическое приближение.
Твердое тело — сложная многочастичная система, состоящая из ядер и электронов. Его можно представить
себе как совокупность сблизившихся атомов так что электронные волновые функции стали перекрываться.
При этом электроны внешних оболочек перестают быть локализованными около своего атома.
Многочастичное уравнение движения нельзя точно решить как в квантовой механике, так и в классической. Поэтому применяют приближения.
Первое приближение, которое мы изучим, позволяет отделить движение электронов от движения ядер.
Это приближение называется адиабатическим. Оно основано на том, что масса электронов m ¿ M , массы
ионов (m/M ∼ 10−4 ÷ 10−5 ).
Мы не будем учитывать релятивистские и спиновые эффекты (парамагнитный и ядерный спиновыe
резонансы и т.д.). Тогда надо решить уравнение Шредингера:
b Ψ̃ = E Ψ̃;
H
Ψ̃ ({R}, {r})
(1.1)
2 X
X ~2 O2
1 X e2
1 X ZJ ZK e2 X ZJ e2
J
b =−~
O2i −
+
+
−
,
H
2m
2MJ
2
rik
2
RJK
riJ
i
r ik = |r i − r k |
J
i6=k
RJK = |RJ − RK |
J6=K
r iJ = |r i − RJ |,
J,i
(1.2)
где m — масса электрона, M - масса иона, Z - заряд иона. Адиабатическое приближение состоит в том,
что мы рассматриваем сначала движение электронов, считая ядра в силу m ¿ M неподвижными. Тогда
b ):
в (1.2) можно пренебречь кинетической энергией ядер (2 слагаемое в гамильтониане H
"
#
~2 X 2
−
Oi + V (r, R) Ψ̃ = εΨ̃,
(1.3)
2m
i
где R является параметром, который входит в Ψ̃ и в ε →, Ψ̃ ({R}, {r}) , ε(R). Полную волновую функцию
ищем в виде
Ψ̃(r, R) = Φ(R)Ψ(r,R)
(1.4)
1
Глава 1. Основные методы и приближения для описания электронных состояний в твердом теле.
Подставив (1.4) в (1.1), используя (1.2) и (1.3), получим
−
¤
~2 X 1 £ 2
ΨOJ Φ + 2(OJ ΨOJ Φ) + ΦO2J Ψ + εΦΨ = EΨΦ.
2
MJ
(1.5)
J
R
R
R
Используя условие нормировки d3 r |Ψ(r, R)|2 = 1 и OJ d3 r |Ψ(r, R)|2 = 2 d3 rΨOJ Ψ = 0 (Ψ всегда
можно выбрать вещественной), умножим (1.5) на Ψ и проинтегрируем по всем r. Получим


Z
2
2
~ X 1 2
~ X 1


−
OJ Φ + ε(R) −
d3 rΨO2J Ψ Φ = EΦ
(1.6)
2
MJ
2
MJ
J
| J
{z
}
Второе слагаемое в [ ] мало, и в адиабатическом приближениии отбрасывается. Оценим его в двух противоположных случаях, полагая, что его величина ограничена полученными при оценке предельными
значениями. Если электроны почти свободны, их волновая функция почти не зависит от R и O2J Ψ → 0.
В худшем случае, когда электроны жестко связаны со своими ионами Ψ(ri , RJ ) = Ψ(ri − RJ )
Z
Z
mi
O2
mi
1
3
2
d rΨOJ Ψ =
d3 rΨ i Ψ =
εel .
(1.7)
2MJ
MJ
2m
MJ
|
{z i }
Т.е. слагаемое неадиабатичности можно отбросить по параметру m/M . Итак, из (1.7) получим
"
#
~2 X O2J
−
+ ε(R) Φ = EΦ
2
MJ
(1.8)
J
Таким образом, получим уравнение для волновых функций ядер, движущихся в потенциале ε(R), созданном только ядрами. Т.е. приближение по m/M оказалось замкнутым. ε включает в себя энергию электронной системы при заданных координатах ядер R и кулоновскую энергию ядер. Таким образом, в
адиабатическом приближении в силу m/M ¿ 1 движение электронов определяется мгновенным положением
ионов, тогда как на медленное движение ионов влияет усредненное распределение электронов и поэтому
задача разбивается на 2 этапа: 1) решение уравнений (1.3) для электронов; 2) затем – уравнения (1.8) для
ядер.
1.2
Приближение самосогласованного поля, метод Хартри-Фока.
Задача после адиабатического приближения все равно остается многочастичной, а электронные переменные
в уравнении Шредингера не разделяются. Приближение самосогласованного поля состоит в предположении,
что каждый электрон "чувствует"некоторое среднее поле, созданное всеми остальными электронами, то
есть в замене многоэлектронного взаимодействия эффективным самосогласованным полем Uef f . Это поле
следует вычислить при помощи самосогласованной процедуры.
Таким образом система взаимодействующих электронов сводится к системе невзаимодействующих
частиц, находящихся в эффективном поле Uef f , а гамильтониан системы распадается на сумму одноэлектронных гамильтонианов. В результате уравнение Шредингера имеет решения в виде произведения
одноэлектронных функций. Но для вычисления этих одноэлектронных функций надо знать Uef f , который
сам зависит от электронных волновых функций ψi (n), где для краткости буквой n обозначена совокупность
координат rn и спинов электрона номер n. Потенциал Uef f находят с помощью самосогласованной процедуры.
Ĥ =
N
X
i=1
Ĥi ,
Ĥi = −
~2 2
O + V (ri ) + Uef f (ri ),
2mi i
2
(1.9)
Глава 1. Основные методы и приближения для описания электронных состояний в твердом теле.
где V (ri ) — потенциальная энергия в поле ядер.
Согласно принципу Паули, волновую функцию нельзя брать в виде произведения функций ψi (n), так
как такое произведение не обладает свойством менять знак на противоположный при перестановке двух
частиц. В качестве решения следует брать антисимметризованную линейную комбинацию произведений
в виде определителя Слетера
¯
¯
¯ ψ1 (1) ψ1 (2) . . . ψ1 (N ) ¯
¯
¯
1 ¯¯ ψ2 (1) ψ2 (2) . . . ψ2 (N ) ¯¯
Ψ(r1 , . . . , rN ) = √ ¯
(1.10)
¯,
...
... ...
N! ¯ . . .
¯
¯ ψN (1) ψN (2) . . . ψN (N ) ¯
где снова 1 = {r1 , s1 }, ψ — одночастичная функция.
Мы будем искать ψi (n) и Uef f вариационным методом, согласно которому решение уравнения Шредингера можно найти с помощью минимизации выражения для энергии системы. Вычислим энергию
системы с помощью (1.10) и исходного гамильтониана
ε = hΨ|Ĥ|Ψi =
N XZ
X
·
ψi∗ (1)
i=1 s1
¸
Z
N
~2 2
1 X X
e2
3
−
Or1 + V (r1 ) ψi (1)d r1 +
d3 r1 d3 r2 |ψi (1)|2
|ψj (2)|2 −
2m1
2
r
12
s ,s
i,j=1
1
2
Z
N
e2
1 X X
d3 r1 d3 r2 ψi∗ (1)ψj (1) ψi (2)ψj∗ (2).
−
2
r12
s ,s
i,j=1
1
(1.11)
2
Здесь использовано условие ортогональности функций {ψi }. Во второе слагаемое дадут вклад только
e2
электроны с параллельными спинами. Это видно из того, что r не зависит от спина, а ψi = ψn ν(s)
12
¡
¢
P ∗
и при суммировании по спину
νi (s)νj (s) = δij ν1 ( 21 ) = 1, ν1 (− 21 ) = 0; ν2 ( 12 ) = 0, ν2 (− 21 ) = 1 . В этом
s↑,↓
можно
µ убедиться
¶
µ и
¶ другим способом, перемножив волновые функции, представленные в виде спиноров
1
0
ψn
, ψn
.
0
1
Первый интеграл в (1.11) — кулоновский — описывает кулоновскую энергию взаимодействия распределений зарядов с плотностями |ψi |2 и |ψj |2 . Второй — обменный — не имеет аналогии в классическом
пределе, он возникает из-за принципа Паули. НаилучшимиR приближениями для ψi являются те, что
3
дают минимум ε при выполнения условий ортогональности
R i (r)ψj (r)d r 3= 0. Для того, чтобы найти
P δψ
условный экстремум проварьируем по ψi функционал ε + εi δψi (r)ψj (r)d r, где εi — неопределенные
i
множители Лагранжа (то есть сделаем замену ψi → ψi + δψi и линеаризуем функционал по δψi ). Получим

−
~2
2m
O2r
+ V (r) +
X 0 Z e2 ψj∗ (r 0 )ψj (r 0 )
0
j
|r − r |

0
d r  ψi (r)−
3
X0
Z
ψj (r)
j
0
0
e2 ψj∗ (r )ψi (r )
0
|r − r |
0
d3 r = εi ψi (r)
(1.12)
(1.12) и есть уравнение самосогласованного поля с
Uef f
X 0 Z e2 |ψj (r 0 )|2
X 0 ψj (r) Z e2 ψj∗ (r 0 )ψi (r 0 )
0
3 0
=
d r −
d3 r .
0
0
ψ
(r)
|r − r |
|r − r |
i
j
j
(1.13)
Видно, что εi имеет смысл одноэлектронной энергии. Причем полная энергия не является суммой энергий
εi , соответствующих одночастичным состояниям. Видно также, что Uef f само зависит от ψj (r), поэтому
3
Глава 1. Основные методы и приближения для описания электронных состояний в твердом теле.
найти его — непростая задача. Обычно пользуются приближениями, чаще всего теорией возмущений.
Надо помнить, что самосогласованное поле — это приближение, существуют также вклады от корреляций
(корреляционная энергия). Приближение самосогласованного поля, называемое также приближением
Хартри-Фока, хорошо работает в приближении плотного газа, когда кулоновская энергия не доминирует
над кинетической.
Можно показать, что кулоновская энергия ECoul ∼ r1s Ry, где rs = ar00 — безразмерное расстояние
между электронами, измеренное в единицах a0 (a0 — боровский радиус a0 = ~2 /me2 = 5.29 × 10−9 cm, r0
— среднее расстояние между электронами, которое находится из условия 1/n = 43 πr03 , n — концентрация
4
' 13.5 эВ — Ридберг), а кинетическая энергия Ekin ∼ r12 Ry . Это видно из
электронов, Ry = me
2~2
s
структуры соответствующих членов, так как кулоновская энергия содержит в знаменателе расстояние
между частицами, а кинетическая энергия — вторую производную по координатам.(Еще имеется энергия
1
нулевых колебаний 3/2 Ry ).
rs
При малых rs (плотный электронный газ) кинетическая энергия больше, а при достаточно больших
rs (разреженный газ) кинетическая энергия становится пренебрежимо малой. Вигнер показал, что в этом
случае электроны располагаются упорядоченным образом — в Вигнеровский кристалл (такой кристалл
был впервые обнаружен в разреженном двумерном электронном газе на поверхности жидкого гелия).
В реальных металлах типичное значение rs составляет несколько единиц, то есть приближение самосогласованного поля применимо качественно, но не количественно, то есть корреляционные поправки не
малы.
Какой характер носит движение электрона в приближении Хартри-Фока? Электрон окружен облаком
других электронов, при этом электроны с тем же спином из-за принципа Паули будут стремиться уйти с
его пути (образуется "обменная дырка"вокруг каждого электрона). Такие "модифицированные"электроны
называются квазичастицами.
Нормальными возбуждениями в системе взаимодействующих частиц являются квазичастицы, волновые
функции которых представляют собой суперпозицию волновых функций многих электронов. Свойства
квазичастиц описываются с помощью теории ферми-жидкости Ландау. Энергия квазичастиц может иметь
мнимую часть (это означает конечное время жизни возбуждений). Как следует из теории ферми-жидкости
Ландау, квазичастицы хорошо определены (затухание мало) в металле только вблизи энергии Ферми.
Квазичастицы обладают многими свойствами обычных частиц.
1.3
Волновая функция электрона в периодическом поле.
Обсудим общие свойства одноэлектронной волновой функции, вытекающие из того, что потенциал в
уравнении Шредингера является периодической функцией координат.
ĤΨ(r) = εΨ(r),
Ĥ = −
~2 2
O + V (r).
2m
(1.14)
Гамильтониан Ĥ инвариантен относительно сдвига решетки на период решетки: V (r) = V (r + an )), an
— любой вектор решетки . Заменим в (1.14) r на r + an , получим, что ψ(r + an ) удовлетворяет тому же
уравнению, что и ψ(r). Если ψ(r) — невырожденное состояние, то
ψ(r + an ) = c(an )ψ(r)
(1.15)
R
В силу условия нормировки |ψ(r)|2 d3 r и (1.15), ψ(r) может отличаться от периодической функции
только фазовым множителем
ψ(r) = eiϕ(r ) u(r);
c(an ) = ei[ϕ(r +an )−ϕ(r )]
4
(1.16)
Глава 1. Основные методы и приближения для описания электронных состояний в твердом теле.
Так как ψ(r + an + am ) = c(an + am )ψ(r), но с другой стороны
ψ(r + an + am ) = c(an ) (c(an )ψ(r)) ,
(1.17)
из (1.16) и (1.17) следует
ei[ϕ(r +an +am )−ϕ(r )] = c(an + am ) = c(an )c(am ) = ei[ϕ(r +an )+ϕ(r +am )−2ϕ(r )] .
Отсюда видно, что ϕ(r) — линейная функция: ϕ = ~kr. Величина ~k является квантовым числом, которое
называется квазиимпульс. Таким образом волновая функция электрона в периодическом поле имеет вид
ψk (r) = eikr uk (r).
(1.18)
Соотношение (1.18) называется теоремой Блоха. Теорема Блоха доказывается и для вырожденных состояний,
при этом вместо (1.15-1.16) следует иметь дело с линейными суперпозициями волновых функций для
состояний с одинаковой энергией.
5
Глава 2
Спектр электронов в твердом теле, зонная
структура.
2.1
Спектр электронов в твердом теле.
Итак, мы показали, что волновая функция одноэлектронного состояния имеет вид
ψmk (r) = eikr umk (r);
k=
p
,
~
(2.1)
где p — квазиимпульс, m — другие квантовые числа (номер зоны), umk (r) — периодическая функция
координат r с периодом кристаллической решетки. Из (2.1) следует, что электрон в кристалле не обладает
определенным импульсом. Действительно, состояние с заданным импульсом описывается плоской волной
exp(ikr), а (2.1) может быть разложено в ряд Фурье
ψmk (r) = eikr
X
uGn eiGn r =
X
uGn ei(k + Gn )r ,
(2.2)
n
n
где uGn — коэффициенты разложения unk (r) в трехмерный ряд Фурье. Формулу (2.2) можно проиллюстрировать на примере одномерного пространства
X
ψnk (x) = eikx
uGn eiGn x ,
n
где Gn = 2π
a n, a — период вдоль оси x для простого кубического кристалла. В трехмерном пространстве
nx ny nz
Gn = 2π ( a , a , a ) , ai — периоды по осям x, y, z. Векторы Gn называются векторами обратной решетки.
x
y
z
Их множество образует множество точек, называемых узлами обратной решетки.
Итак, согласно (2.2), блоховский электрон является суперпозицией состояний с импульсами, отличающимися на величину векторов обратной решетки Gn .
Для произвольного кристалла с набором основных векторов прямой решетки (a1 , a2 , a3 ), обратную
[a ,a ]
решетку образуют вектора b1 = 2π (a [a2 ,a3 ]) ; остальные вектора получаются циклической перестановкой
1 2 3
[a
,a
]
bi = 2π (a [ai+1 ,ai−1 ]) . Построение обратной решетки рассматривалось в прошлогоднем курсе по физике
i i+1 i−1
твердого тела. В частности, было показано, что для объемноцентрированной кубической решетки обратная
решетка являлась кубической гранецентрированной.
Из (2.2) также видно, что значения k, отличающиеся на один из векторов обратной решетки Gn ,
описывают ту же самую волновую функцию. Поэтому вектор k достаточно выбрать в пределах одной
6
Глава 2. Спектр электронов в твердом теле, зонная структура.
Рис. 2.1: Простой пример ячейки Вигнера-Зейца
элементарной ячейки обратной решетки. Чаще всего выбирают элементарную ячейку в виде ячейки
Вигнера-Зейтца. Ячейка Вигнера-Зейтца (ЯВЗ) — это область пространства, лежащая к данной точке
обратной решетки ближе, чем к какой-либо другой точке решетки.
Чтобы построить ЯВЗ, необходимо соединить один узел со всеми соседними и провести перпендикулярные прямые через середины прямых, соединяющих узлы. Область, ограниченная перпендикулярами
и есть ЯВЗ. ЯВЗ обладает симметрией решетки. ЯВЗ для обратной решетки называют первой зоной
Бриллюэна (из чего следует, что аналогичным образом можно определить вторую, третью и следующие
зоны Бриллюэна). В задачах квантовой механики часто использовались периодические граничные условия
Борна-Кармана для волновых функций свободного электрона — плоских волн. Эти периодические условия
часто используются и для функций Блоха. Если число элементарных ячеек в кристалле N = Nx Ny Nz
(где Ni — число ячеек по оси i), то из условий Борна-Кармана следует
ψnk (r + Ni ai ) = eiNi ki ai ψnk (r).
Таким образом eiNi ki ai = 1 и допустимыми значениями k являются
¾
½
mi 2π
k=
.
Ni ai
Так как разница между соседними значениями k мала, порядка 1/Ni , суммирование по k можно
заменить интегрированием
Z
X
d3 k
→V
.
(2π)3
k
2.2
Модель почти свободных электронов.
Изучим, какой вид имеет энергетический спектр электронов и его основные свойства на примере простейших
моделей.
Рассмотрим электроны в слабом периодическом потенциале V (r), который можно учитывать по теории
возмущений.
В отсутствии взаимодействия, когда электрон свободен, его энергия и волновая функция имеют вид
ε0 (k) =
~2 k 2
,
2m
ψ(r) = eikr .
(2.3)
7
Глава 2. Спектр электронов в твердом теле, зонная структура.
Отметим, что ψ(r) в (2.3) может рассматриваться как удовлетворяющая теореме Блоха с u(r) = 1).
По теории возмущений, по потенциалу V (r) спектр описывается в виде
0
ε(k) = ε (k) + hk|V |ki +
X 0 |hk|V |k0 i|2
k
0
ε0 − ε0 0
k
k
= ε0 (k) + V0 +
X
|VG |2
,
ε0 − ε0k+G
G6=0 k
(2.4)
0
0
RP
0
V (r)ei(k − k )r d3 r =
VG ei(G + k − k )r d3 r 6= 0 только, если k − k
G
равен одному из векторов обратной решетки.
Из (2.4) видно, что даже при самом малом потенциале V поправка второго порядка расходится при
ε0k = ε0k − G , т.е. при |k| = |k − G|. Поэтому надо использовать теорию возмущений для вырожденных
состояний. Ищем волновую функцию в виде
0
где мы учли, что hk |V |ki =
R
Ψk = ak eikr + ak − G ei(k − G)r .
(2.5)
Подставим (2.5) в уравнение Шредингера
· 2 2
¸
~ O
−
+ V (r) Ψk = ε(k)Ψk ,
2m
умножим уравнение Шредингера сначала на e−ikr , а потом на e−i(k − G)r и проинтегрируем. Получим
два уравнения на коэффициенты a:
( £
¤
ε0 (k) − ε(k) ak + VG ak − G = 0
£ 0
¤
(2.6)
V−G ak + ε (k − G) − ε(k) ak − G = 0
Решение существует, когда определитель (2.6) равен нулю.
¤ 1q
1£ 0
0
ε(k) =
ε (k) + ε (k − G) ±
[ε0 (k) + ε0 (k − G)]2 + 4|VG |2
2
2
(2.7)
На рис. 2.2. представлен спектр для 1D случая. При малых VG вдали от границы зоны Бриллюэна ε(k) ≈
ε0 (k), а возле границы, вследствие брэгговского отражения, возникает запрещенная зона (энергетическая
щель) шириной 2VG . Энергетическая щель возникает на границе всех зон Бриллюэна. Можно ограничиться
описанием спектра только в пределах первой зоны Бриллюэна, сместив куски, отвечающие следующим
зонам, на соответствующий вектор обратной решетки (см. рис.2.2).
2.3
Приближение сильной связи.
Приближение сильной связи(tight-binding approximation) применимо к случаям, когда перекрытие атомных
волновых функций в кристалле достаточно мало, то есть когда вблизи атома волновая функция электрона
в кристалле мало отличается от атомной (применимо для многих диэлектриков и частично заполненных
d- оболочек переходных металлов).
Итак, пусть связанные электронные состояния атома в начале координат описываются уравнением
Шредингера:
Ĥat Ψn (r) = En Ψn (r),
где Ĥat = −
1 2 X
∇ +
Ûat (r).
2m
n
8
(2.8)
Глава 2. Спектр электронов в твердом теле, зонная структура.
Рис. 2.2: Электронный спектр в приближении слабой связи для 1D случая
Полный гамильтониан уравнения Шредингера, описывающего электронные состояния кристалла, отличается
от атомного
1 2 X
Ûat (r − R),
ĤΨ = εΨ,
(2.9)
Ĥ = Ĥat (r)+ M U (r) ≈ −
∇ +
2m
R
где M U (r) является разницей между потенциалом атома и периодическим потенциалом кристалла, а знак
примерного равенства используется потому, что в M U (r) кроме потенциалов других атомов входят также
поправки, связанные с изменением потенциала при сближении атомов, например, вследствие деформации
электронных оболочек и других причин. Величина M U (r) мала в области локализации Ψn из (2.8), а там,
где потенциал M U (r) велик, малы Ψn , поэтому его матричные элементы на Ψn , локализованных на
разных атомах, малы. Используя эти обстоятельства, можно построить теорию возмущений.
В пренебрежении перекрытиями волновых функций спектр состоит из N -кратно вырожденных уровней
с волновыми функциями Ψn (r − R), где N — число атомов в кристалле.
Действуя в духе теории возмущений, для вырожденных состояний функции нулевого приближения
следует искать в виде линейных комбинаций атомных функций:
X
(0)
Ψnk (r) =
eikR Ψn (r − R)
(2.10)
R
Функций (2.10) всего N штук (по числу допустимых значений k), причем функция вида (2.10) удовлетворяет
условию Блоха:
Ψnk (R + r) = eikR Ψnk (r),
поскольку
Ψnk (r + R) =
X
0
eikR Ψn (r + R − R) =
0
R0
X
R̃
9
eikR eikR̃ Ψn (r − R̃) = eikR Ψnk (r),
Глава 2. Спектр электронов в твердом теле, зонная структура.
0
где R̃ = R − R Вообще говоря, любую блоховскую функцию можно представить в виде, похожем на
(2.10), а именно
Ψnk (r) =
X
eikR ϕn (r − R)
(2.11)
R
Эта формула соответствует разложению функции от R в трехмерный ряд Фурье по k, где ϕn (r − R)
являются коэффициентами разложения, которые можно найти по обычной формуле обращения:
Z
Z
3
1
1
d3 k
−ikR
ik(r − R) u (r) d k ,
=
(2.12)
e
Ψnk (r)
e
ϕn (r − R) =
nk
V0
(2π)3
V0
(2π)3
где V0 - объем первой зоны Бриллюэна в k-пространстве. Функции ϕn (r) называются функциями Ваннье.
Непосредственным интегрированием можно показать, что функции Ваннье ортонормированы (и, следовательно,
образуют полную систему функций)
Z
0
ϕ∗n (r − R)ϕm (r − R )d3 r = δnm δRR0 .
Функции Ваннье локализованы. Это видно из (2.12), так как при больших |r − R| последний интеграл
содержит быстро осциллирующую функцию от k и поэтому мал. Из формулы (2.10) следует, что в нулевом
приближении по взаимодействию между атомами функция Ваннье совпадает с атомной функцией.
Приступим теперь к решению уравнения Шредингера. Естественно, его решение надо искать в форме,
которая сохраняет общий вид решения (2.11)
Ψ(r) =
X
eikR φ(r − R).
(2.13)
R
При этом нам надо будет найти функцию φn , которую мы будем искать в виде ряда по атомным функциям:
X
bm Ψm (r)
(2.14)
φ(r) =
m
что
Умножим теперь уравнение (2.9) на Ψ∗n (r), проинтегрируем по всем r и используем то обстоятельство,
Z
Ψ∗n (r)Hˆat Ψ(r)d3 r
Z h
Z
i∗
3
ˆ
=
Hat Ψn (r) Ψ(r)d r = En Ψ∗n (r)Ψ(r)d3 r.
Получим
Z
(ε(k) − En )
Z
Ψ∗n (r)Ψ(r)d3 r =
Ψ∗n (r) M U (r)Ψ(r)d3 r
(2.15)
Подставив (2.13) и (2.14) в (2.15) и используя ортонормированность атомных волновых функций
hm|ni = δmn ,
приходим к задаче на собственные значения

(ε(k) − En ) bn = − (ε(k) − En )
X
m

XZ

Ψ∗n (r)Ψm (r − R)eikR d3 r  bm +
R6=0
10
Глава 2. Спектр электронов в твердом теле, зонная структура.
+
X µZ
m
+
X


m
Ψ∗n (r)
XZ
¶
ikR
3
M U (r)Ψm (r)e
d r bm +

Ψ∗m (r) M U (r)Ψm (r − R)eikR d3 r  bm
(2.16)
R6=0
Видно, что первый и третий интегралы в правой части (2.16) малы, т.к. содержат интегралы, содержащие
волновые функции, относящиеся к разным атомам. Второй интеграл мал, т.к. содержит потенциал M
U (r), который не мал только в той части пространства, где малы Ψn (r). Таким образом левая часть
(ε(k) − En ) bn всегда мала. То есть мало либо ε(k) − En , либо bn . Иными словами, если энергия εk близка
к энергиям атомного уровня En , то bn может быть не мало. Если ε(k) − En велико, то соответствующее
bn можно выбросить за малостью. Итак, выбросим bm с m 6= n. Это соответствует тому, что в методе
сильной связи функция Ваннье заменяется атомной функцией.
В случае, если имеется несколько вырожденных уровней или уровней с близкой энергии, наша задача
свелась бы к системе уравнений. Мы для простоты ограничимся случаем, когда имется единственный
невырожденный атомный уровень. В этом случае получим довольно простую формулу для энергии
P
εn (k) = En −
R6=0
1+
γ(R)eikR + β
P
(2.17)
α(R)eikR
R6=0
Z
β=−
Z
3
2
d r M U (r)|Ψn (r)| ;
d3 rΨ∗n (r) M U (r)Ψn (r − R);
γ(R) =
Z
α(R) =
d3 rΨ∗n (r)Ψn (r − R).
Взяв в (2.17) R, относящиеся к ближайшим соседям, получим результат в приближении сильной связи.
Например, пренебрегая малыми поправками с α, получим
X
ε(k) = En − β −
γ(R) cos kR
(2.18)
ближ.соседи
Спектр ε(k) в (2.18) зависит от структуры кристалла. Например, для простой кубической решетки, сумма
в (2.18) включает 6 атомов
ε(k) = En − β − 2γ (cosx a + cosy a + cosz a)
(2.19)
В объемно-центрированной кубической решетке, суммируя вклады от всех ближайших соседей, мы получим
из (2.18)
ε(k) = En − β − 8γ cos
ky a
kx a
kz a
cos
cos
2
2
2
(2.20)
Обратим внимание на то, что спектр ε(k) обладает той же симметрией, что и сам кристалл. Например,
если кристалл обладает центром симметрии, то выполняется соотношение εs (−k) = εs (k). Как видно из
(2.18), а также из приближения слабой связи (2.7), закон дисперсии имеет нулевой наклон ∂ε/∂k вблизи
границы зоны Бриллюэна. Это является общим правилом, не зависящим от приближения, и следует из
того, что состояния, отличающиеся по импульсу на вектор обратной решетки, эквивалентны.
Вывод из двух противоположных подходов к вычислению спектра одинаков:
11
Глава 2. Спектр электронов в твердом теле, зонная структура.
— в приближении сильной связи, при сближении атомов, атомные уровни расщепляются в разрешенные
энергетические зоны;
— в приближении слабой связи появляются запрещенные зоны энергии, разделяющие разрешенные
зоны.
Таким образом спектр имеет зонный характер: чередуются запрещенные зоны энергии и разрешенные
зоны энергии с законом дисперсии ε(k). Отличия спектров в двух рассмотренных предельных случаях
являются не качественными, а всего лишь количественными: при слабой связи разрешенные зоны оказываются больше запрещенных, а в пределе сильной связи, наоборот, разрешенные зоны меньше запрещенных.
12
Глава 3
Свойства блоховских электронов
3.1
Металлы, диэлектрики, полупроводники.
Итак, при сближении атомов, уровни энергии расщепляются и образуются зоны. Становится ясным откуда
возникают понятия 2s-зона, 3p-зона и т.д., это указания на атомные термы, из которых произошла данная
зона.
Различные зоны могут перекрываться или оставаться разделенными зонами запрещенных энергий.
Пусть зоны не перекрываются. Тогда из полностью заполненных (полностью пустых, частично заполненных) атомных термов образуются полностью заполненные (соответственно, полностью пустые или
частично заполненные) зоны. Если же зоны перекрываются (гибридизация зон), то из атомного терма,
занятого электроном, и терма с незаполненным состоянием может образоваться частично заполненная
электронами полоса энергии. Согласно принципу Паули, при T = 0 в зоне будут заняты ZN/2 энергетических состояний с наименьшей энергией, где N — число атомов, Z — число электронов на соответствующих
уровнях в атоме, 2 возникла из-за спина. Всего же в одной зоне Бриллюэна содержится N состояний с
различными значениями k. Таким образом, по заряду иона Z можно судить о характере заполнения
зоны. Например, при нечетном Z заведомо появятся частично заполненные зоны. И действительно, такая
ситуация имеет место, например, в щелочных металлах, где на верхнем заполненном уровне имеется один
электрон (Z = 1).
В зависимости от степени заполнения зоны кристалл можно поделить на металлы и неметаллы. Дело в
том, что на основе вида зонной структуры можно судить о характере отклика кристалла на электрическое
поле. Характерная величина электрического поля в атоме eEa ∼ Vchar /a, где a — межатомное расстояние
∼ 10−8 ÷ 10−7 см, а Vchar — характерные атомные энергии ∼ 1 ÷ 10 эВ. Таким образом, атомное поле
— огромное, Ea ∼ 107 ÷ 108 В/см. Если приложенное внешнее поле E ¿ Ea (что бывает почти всегда),
его влияние можно учитывать по теории возмущений, согласно которой, вероятность перехода между
0
состояниями с энергией εn (k) и εm (k ) равна
W =
2π
0
0
|hn, k|V |m, k i|2 δ(εn (k) − εm (k ) − ~ω)
~
где ω — частота поля. В слабых V (r), плавно изменяющихся в пространстве (то есть мало изменяющихся
на масштабе a) и имеющих частоты ω ¿ εG /~ (εG — запрещенные зоны), согласно (3.1), переходы
возможны только между состояниями одной зоны и близкими k. Кроме того, согласно принципу Паули,
в одном квантовом состоянии не может находиться более одного электрона, поэтому переходы возможны
только из заполненного состояния в пустое. Таким образом, если все зоны либо полностью заполнены,
либо полностью пустые, переходов между состояниями, вызванных электрическим полем, не будет. Такой
кристалл является диэлектриком.Если же зона заполнена частично, в ней возможны переходы между
состояниями с разными p. Это — металл.
13
Глава 3. Свойства блоховских электронов
Если εG не очень велика и возможны перебросы из заполненной зоны в пустую вследствие термической
активации (тепловых забросов), то появляется возможность для небольшого числа переходов — это полупроводник (как бы частный случай диэлектрика). При слабом перекрытии пустой и заполненной зон
образуется металл с малым количеством заполненных состояний в верхней зоне и пустых в нижней. Это
— полуметалл.
3.2
Динамика блоховского электрона.
Выясним теперь, как связана скорость электрона hṙi со спектром E(k).
Вычислим сначала вид оператора r̂ в импульсном представлении, для чегоPразложим
произвольную
R
csk ψsk d3 k.
функцию ψ, на которую действует оператор r̂, по блоховским функциям, ψ =
s
rψ =
XZ
csk rψsk d3 k =
s
XZ
s
µ
¶
∂ψsk
ikr ∂usk
csk −i
+ ie
d3 k,
∂k
∂k
(3.1)
Затем первый член проинтегрируем по частям, а во втором разложим периодическую функцию ∂usk /∂k
по usk с тем же k :
X
∂usk
0
= −i
hsk|Ω̂|s kius0 k .
∂k
0
s
Тогда
rψ =
XZ
s

X
∂c
0
i sk +
hsk|Ω̂|s kics0 k  ψsk d3 k.
∂k
0

(3.2)
s
Так как по определению
XZ
(rcsk )ψsk d3 k,
rψ =
s
то
r̂ = i
∂
+ Ω̂,
∂k
(3.3)
0
где Ω̂ задается матрицей hsk|Ω̂|s ki. Оператор
Ω̂ диагонален
по k поэтому коммутирует с k̂.
h
i
i
Вычислим оператор скорости: v̂ = ~ Ĥr̂ − r̂Ĥ . Вычисляя коммутатор, получим:
1
v̂ = −
~
µ
¶
¤
∂
∂
i£
0
εs (k)
−
εs (k) +
εs (k) − εs0 (k) hsk|Ω̂|s ki.
∂k ∂k
~
(3.4)
0
Первое слагаемое дает ∂εs (k)/~∂k, а второе обращается в нуль при s = s . Значит средняя скорость в
зоне s, которая описывается диагональным матричным элементом, имеет вид
v̄s =
∂εs (k)
∂εs (k)
=
,
~∂k
∂p
(3.5)
как в классической физике. Более того, в квазиклассическом пределе, в (3.4) ∂/∂k велико, так как
характерные изменения волнового вектора M k малы в случае плавных и медленных возмущений, не
14
Глава 3. Свойства блоховских электронов
вызывающих межатомные переходы. То есть Ω̂ можно отбросить по сравнению с ∂/∂k. В результате
в квазиклассическом пределе получаем для блоховского электрона такую же связь квазиимпульса с
координатой, как для импульса с координатой для обычного электрона. Отсюда, в частности, следуют
правила коммутации как в обычной квантовой механике и возможность рассматривать εn (k) как оператор.
В r - пространстве кинетическая энергия имеет вид εn (−iO), а уравнение Шредингера:
[εn (−iO) + Uext (r)] ψ = Eψ.
(3.6)
Конечно, формула (3.6) справедлива тогда и только тогда, когда Uext (r) — плавная функция, мало
меняющаяся на межатомных расстояниях.
Можно учесть и влияние слабого магнитного поля (такого, что радиус циклотронной орбиты À a).
Известно, что при градиентных преобразованиях векторного потенциала A → A + Of , меняется волновая
ie
функция ψ → ψ exp( f ). В слабом поле H медленно меняется в пространстве A(r). В пределе H →
~c
0 A = const ≡ A0 . Переход от A = 0 к A = A0 соответствует преобразованию с f = A0 r, поэтому
e
e
из функций ψsk = usk eikr получается usk exp(i(k + A0 )r). K = k + A0 можно назвать обобщенным
~c
~c
квазиимпульсом. Таким образом волновые функции имеют вид ψs,K−eA0 /~c exp(i(Kr)), а соответствующие
им значения энергии εs (k) = εs (K − eA/~c). При непостоянном, но плавном A, такие ψ являются
волновыми функциями нулевого приближения по H, причем usk из-за A(r) уже не периодичны. Энергии
εs (K − eA/~c) можно рассматривать как оператор, образующий гамильтониан в K — представлении
r̂ = i∂/∂K
Ĥs = εs (K −
eA
) + U (r)
~c
r̂ = i
∂
∂K
или в r - представлении:
Ĥs = εs (−i~O −
eA
) + U (r)
~c
~ω ¿ εG
(3.7)
Операторы r̂ = i∂/∂K и K̂ ≡ Kимеют те же правила коммутации, что и координаты, и обобщенные
импульсы свободных частиц. Поэтому сохраняются все обычные правила квантовой механики, в частности,
ˆ = − ∂ Ĥ ,
~K̇
∂r
~ˆṙ =
∂ Ĥ
∂K
(3.8)
и в квазиклассическом пределе имеет смысл понятие траектории частицы, а также импульса и энергии,
зависящих от координат. В этом случае справедливы обычные уравнения движения:
~K̇ = ṗ = F,
v=
∂ε
.
∂p
(3.9)
Можно также говорить, например, что дно зоны проводимости и потолок валентной зоны зависят от
координат и рисовать диаграммы, в которых уровни энергии зависят от координат (не путать с законом
дисперсии E(p)!) Важно помнить, что такой подход справедлив лишь в квазиклассическом пределе, когда
отбрасывается недиагональный по номеру зоны оператор Ω̂!
3.3
Эффективная масса.
В металлах, где зона заполнена на величину порядка половины зоны, играют роль электроны вблизи
энергии, соответствующей границе между заполненными и незаполненными состояниями в зоне — энергия
Ферми. Характерная величина εF — порядка атомных энергий ∼ 1 ÷ 10 эВ → 105 K. Электронное
15
Глава 3. Свойства блоховских электронов
Рис. 3.1: Функция распределения Ферми
распределение Ферми f = 1/(exp(
ε − εF
)+1) — ступенька, размытая на величину порядка T (T < 103 K ¿
kT
εF /k) (см. рис. 3.1).
Поэтому все интересные процессы, в том числе проводимость, происходят за счет электронов с энергиями вблизи εF и можно забыть о наличии других зон. Именно вблизи εF хорошо определены квазичастицы
и хорошо работает теория ферми-жидкости Ландау. В этом причина того, что очень часто качественно
хорошо работает и модель свободных электронов Зоммерфельда, в которой металл рассматривается
как состоящий из свободных, невзаимодействующих электронов. Отклонения от модели Зоммерфельда,
связанные с взаимодействиями квазичастиц, не включенными в приближение самосогласованного поля,
приводят к таким эффектам, как сверхпроводимость.
В полупроводниках имеется небольшое число занятых и (или) пустых максимумов вблизи экстремума
зоны Бриллюэна. Поэтому важную роль в процессах переноса играет именно эта область энергий. Но
вблизи экстремума, в невырожденной зоне, E(p) можно разложить в ряд:
¯
1 X −1
1 ∂ 2 En (p) ¯¯
−1
mαβ (p − p0 )α (p − p0 )β , mαβ = 2
En (p) = En (p0 ) +
.
(3.10)
2
~ ∂pα ∂pβ ¯p=p0
α,β
Коэффициенты разложения m−1
αβ образуют тензор эффективной массы. Вблизи экстремума в главных
осях эта формула приобретает совсем простой вид:
En (p) − En (p0 ) =
X (p − p0 )2
α
α
2mα
,
а скорость электрона определяется как vα = ∂E/∂pα = pα /mα .
Эффективную массу можно вычислить и связать с шириной запрещенной зоны. Это делается с помощью
подхода, называемого kp - методом.
Запишем уравнение Шредингера для блоховской волновой функции в виде уравнения для периодической
функции upn (r) :
1
(−i~O + p)2 upn (r) + V (r)upn (r) = En (p)upn (r)
2m
(3.11)
Уравнение для волновой функции u0n (r), относящейся к экстремуму энергетической зоны имеет вид
¸
·
1
(−i~O)2 + V (r) u0n (r) = En (0)u0n (r).
(3.12)
2m
16
Глава 3. Свойства блоховских электронов
Формула (3.11) отличается от (3.12) членами −i~pO/m + p2 /2m, которые вблизи экстремумов малы и их
можно вычислить по теории возмущений.
Во втором порядке теории возмущений из уравнения (3.11) получим
En (p) = En (0) +
p2
1 X 0 |phum0 | − i~O|un0 i|2
+ 2
2m m m
En (0) − Em (0)
(3.13)
Так как по определению эффективной массы
1
En (p) = En (0) + m−1
pα pβ ,
2 αβ
(3.14)
из уравнения (3.13) находим
m−1
αβ =
δαβ
2 X 0 hun0 | − i~Oα |ul0 ihul0 | − i~Oβ |un0 i
+ 2
m
m
En (0) − El (0)
(3.15)
l
Матричные элементы по порядку величины равны ~/a, поэтому для эффективной массы получаем оценку
¶¸
·
µ
2~2
∗−1
−1
,
m
=m
1+O
ma2 ME
При ~2 /(ma2 ) ∼ 10 eV, и M E ∼ 0.2 eV получим m∗ ∼ 10−2 m0 . Таким образом, в веществах с малой EG
мала эффективная масса.
3.4
Зонная структура типичных полупроводников.
1) Полупроводники InSb, GaSb, GaAs имеют кубическую структуру и в зоне проводимости p0 = 0, а
mx = my = mz ;
2) У полупроводников Ge, Si в зоне проводимости p0 6= 0. В этом случае существует несколько экстремумов, расположенных симметрично. Изоэнергетическая поверхность состоит из нескольких эллипсоидов.
В Ge дно зоны проводимости находится на объёмных диагоналях куба [111] (всего 6 штук), а в Si на
оси [100] (8 штук) и ей эквивалентных. Изоэнергетические поверхности — эллипсоиды вращения с осью
вращения вдоль [111] у Ge и [100] в Si.
E(p) =
p21 + p22
p2
+ 3
2m1
2m3
В Ge m1 = m2 = 0.082m0 , m3 = 1.58m0 Во многих полупроводниках валентная зона имеет довольно
сложный характер, связанный с вырождением, снятым спин - орбитальным взаимодействием. Так и в Ge,
и в Si вблизи p = 0 в валентной зоне
o
1 n 2 q 2 4
1
E1,2 (p) = −
Ap ± B p + C 2 (p2x p2y + p2y p2z + p2z p2x ) ; E3 (p) = −Es0 −
Ap2
2m0
2m0
Ge : Es0 = 0.28 eV A = 13, B = 8.7, C = 11.4
Si : Es0 = 0.035 eV A = 4.1, B = 1.4, C = 3.7
Изоэнергетические поверхности гофрированы (эффективная масса зависит от направления).
17
Глава 3. Свойства блоховских электронов
Рис. 3.2: Плотность состояний в приближении сильной связи для кубической решетки
3.5
Плотность состояний.
Плотность состояний N (E) определяется таким образом: N (E)dE — это число квантовых состояний в
единице объема в маленьком интервале энергий от E до E + dE.
Мы знаем, как найти число состояний в элементе объема в фазовом пространстве, поэтому можем
вычислить, сколько состояний находится в интервале энергий dE:
Z
X
d3 k
≡ N (E)dE
(3.16)
(2π)3
n
E<En (k)<E+dE
Если выбрать элемент объема в k - пространстве d3 k = dsdk⊥ , где ds — элемент изоэнергетической
поверхности, а dk⊥ — приращение k по нормали к этой поверхности, то
Z
ds
1 X
(3.17)
N (E) =
3
(2π) n
|Ok En (k)|
En (k)=E
В простейшем случае сферической изоэнергетической поверхности N (E) вычисляется путем простой
замены переменных k → E) в (3.16):
E = ~2 k 2 /2m, d3 k = 4πk 2 dk, dE = ~2 kdk/m,
√ ∗3/2
√
2m
N (E) =
E.
2
3
π ~
В простой кубической решетке (в приближении сильно связанных электронов):
Z
ds
1
q
,
N (E) =
16π 3 γa
2
2
2
sin
ak
+
sin
ak
+
sin
ak
x
y
z
S
где интегрирование ведется по поверхности
E = −2γ(cos akx + cos aky + cos akz )
Точки излома на зависимости N (E) называются сингулярностями Ван - Хова (см.рисунок).
18
Глава 4
Приближение эффективной массы.
4.1
Электроны и дырки.
Ток, создаваемый одним электроном в состоянии k есть:
−ev(k).
Сумма всех токов, создаваемых электронами в полностью заполненной зоне,
XZ
X
dk
v(k) = −eV
v(k) 3 = 0,
j = −e
2π
σ
(4.1)
k,σ
где интеграл обращается в ноль из-за того, что он представляет собой интеграл от производной от
периодической функции по периоду,
1 dE(k)
.
v(k) =
~ dk
Если вероятность заполнения данного состояния k, σ электроном (называемая функцией распределения
электронов) есть ne (k, σ), то выражение для тока
X
X
X
X
nh (k)v(k),
(4.2)
nh (k)v(k) = e
v(k) + e
ne (k)v(k) = −e
j = −e
k,σ
k,σ
k,σ
k,σ
где nh (k) = 1 − ne (k) называется функцией распределения дырок, то есть состояний, не занятых электронами. Видно, что ток можно описать, как результат движения положительных зарядов — дырок.
Аналогично можно вычислить поток энергии W при наличии электрического потенциала ϕ(x)
X
X
ne (k, σ)[E(k) + eϕ]v(k) =
[1 − nh (k, σ)][E(k) + eϕ]v(k) =
W=
k,σ
k,σ
X
X
[E(k) + eϕ]v(k) −
nh (k, σ)[E(k) + eϕ]v(k),
где
P
k,σ
(4.3)
k,σ
[E(k) + eϕ]v(k) = 0, поскольку E(k)v(k) = ∇E 2 (k)/(2~), так что суммирование дает ноль по тем
k,σ
же причинам, что и в (4.1). Последнее слагаемое в (4.3) может быть интерпретировано как поток энергии
дырок — частиц с энергией −E(k) и зарядом |e|.
P
Вблизи максимума Ek = Ek0 +
α
~2 kα2
2|mα | , где |mα | — эффективная масса дырки. Мы раньше пришли
к выводу , что в классическом приближении
e
~k̇ = ṗ = F = [v, H] + eE,
c
19
Глава 4. Приближение эффективной массы.
∂E(k)
где v = ~∂k . Это уравнение вблизи максимума зоны также можно интерпретировать как уравнение
∂2E
движения частицы с положительной массой −1/
, положительным зарядом и положительной энергией,
∂pα ∂pβ
отсчитанной вниз от максимума. Вычисление эффективной массы вблизи экстремума можно провести с
помощью kp метода(см.3.11 - 3.15).
4.2
Уравнение приближения эффективной массы.
Как было показано, в случае плавно меняющихся на межатомных расстояниях полей приближенно справедливы уравнения, аналогичные квантовой механике свободного электрона, в которых квазиимпульс
играет роль импульса:
r̂ = i
∂
∂
, и K̂ = −i ,
∂k
∂r
а гамильтониан для электронов в зоне n в пренебрежении межзонными переходами (слагаемыми типа Ω)
имеет вид (см. (3.6),(3.7) )
[εn (−i~Or ) + Uext (r)]ψ = Eψ,
(4.4a)
в r-представлении или
[εn (p) + Uext (i~Op )]ψ = Eψ,
(4.4b)
в p-представлении.
Вблизи экстремума εn (p) можно разложить в ряд
εn (p) = εn (p0 ) +
X (p − p0 )α (p − p0 )β
(4.5)
2mαβ
αβ
Обозначим E = εn (p) + ε , где ε — энергия, отсчитанная от экстремума. Подставив (4.5) в (4.4b), получим


X (p − p0 )α (p − p0 )β
∂

+ Uext (i~
) ψ(p) = εψ(p).
(4.6)
2mαβ
∂(p − p0 )
αβ
То есть для ψ(p) получено почти обычное уравнение Шредингера, которое отличается от обычного
заменой массы электрона на эффективную. Это уравнение записано в квазиимпульсном представлении.
Можно вычислить ψ(r), для чего разложим ее по блоховским функциям Ψnp (r):
X
X
ψ(r) =
Φn (p)Ψnp (r) =
Φn (p)eipr unp (r).
(4.7)
np
n,p
Так как мы всюду пренебрегали межзонными переходами, то вклад других зон надо опустить и в (4.7).
Получим, что
X
ψ(r) '
Φn (p)eipr unp (r).
(4.8)
p
Так как в сумму по p входят только малые p−p0 , то в unp (r) можно в главном приближении не учитывать
также и зависимость от p и вынести unp0 (r) за знак суммы
X
ψ(r) ' unp0 (r)eip0 r
Φn (p)ei(p−p0 )r ≡ Φn (r)Ψnp0 (r),
(4.9)
p
20
Глава 4. Приближение эффективной массы.
где Φn (r) и Φn (p) связаны друг с другом преобразованием Фурье.
Вычислим теперь, какой вид имеет функция ψ(r) в квазиимпульсном представлении, используя ее
приближенный вид (4.9).
Z
Z
∗
ψ(p) ≡ drψ(r)Ψnp (r) = drΦn (r)Ψnp0 (r)Ψ∗np (r) =
Z
Z
i(p0 −p)r
∗
= drΦn (r)e
unp0 (r)unp (r) ≈ drΦn (r)ei(p0 −p)r unp0 (r)u∗np0 (r) ≈
Z
X
X
i(p0 −p)rn
≈
Φn (rn )e
drunp0 (r)u∗np0 (r) =
Φn (rν )ei(p0 −p)rν ≈ Φn (p),
(4.10)
rn
rν
где мы вынесли часть подынтегрального выражения, медленно меняющуюся в пределах одной элементарной
ячейки, в точках rν , относящихся к элементарной ячейке номер ν и учли нормировку функций Ψnp0 (r).
Так как ψn (p) удовлетворяет уравнению Шредингера (4.6) в импульсном представлении, то подставив
(4.10) в (4.6)и сделав обратное преобразование Фурье над (4.6), мы получим уравнение Шредингера вида


2
X
Oα Oβ
− ~
+ Uext (r) Φn (r) = εΦn (r).
(4.11)
2
mαβ
αβ
Таким образом функция Φn (r), которая называется огибающей функцией, удовлетворяет уравнению Шредингера с эффективной массой в координатном пространстве, причем это уравнение не содержит периодического потенциала. Учет периодического потенциала свелся к появлению эффективной массы
электрона вместо массы свободного электрона. Полная волновая функция, согласно (4.9), вычисляется
как произведение огибающей функции на блоховскую функцию, относящуюся к экстремуму зоны.
Уравнение (4.11) называется уравнением Шредингера в приближении эффективной массы. Следует
помнить, что приближение эффективной массы применимо только в том случае, если Uext (r) и, соответственно, Φn (r) — плавные функции, которые мало меняются на расстояниях порядка межатомных, и
когда характерные энергии ∼ p2 /2m∗ малы по сравнению с расстояниями до других зон. Только в этом
случае можно пренебречь межзонными матричными элементами типа hm|Ω|ni.
В дальнейшем мы будем опускать индекс n у огибающей, чтобы не было путаницы с другими квантовыми числами, которые характеризуют функцию Φ и энергию ε.
4.3
Примесные атомы.
Различные точечные дефекты (вакансии, лишние атомы в междоузлиях, примеси замещения и т.д.)
создают отклонение от периодического потенциала — помимо периодического потенциала появляется
внешний потенциал U (r) в уравнении Шредингера. Если этот потенциал плавный, то изменения в спектре
можно найти по методу эффективной массы. Важным примером такого случая являются примесные
атомы, создающие кулоновский потенциал. Например, элементы III (B,Al,Ga,...) и V (P,As,Sb,...) групп,
добавленные к Ge и Si, образуют твердый раствор замещения, замещая атом вещества-хозяина. В результате
у элемента V группы (донор) остается "лишний"электрон, не участвующий в ковалентных связях. У
элемента III группы, наоборот, будет не хватать электрона для заполнения связей. При удалении лишнего
электрона (в случае III группы — дырки) ионизированный примесный атом будет создавать кулоновское
поле. Причем на расстояниях больше межатомных кулоновское поле будет ослабляться поляризацией
кристалла, что учитывается введением макроскопической диэлектрической проницаемости κ
U (r) = −
e2
,
κr
(4.12)
21
Глава 4. Приближение эффективной массы.
где потенциал U (r) (4.12) надо вставить в уравнение приближения эффективной массы. Для лишнего
электрона (донорный атом) уравнение Шредингера в приближении эффективной массы записывается
как


2
2
X
Oα Oβ
e
− ~
+  Ψ(r) = εΨ(r).
(4.13)
2
mαβ
κr
αβ
Решение (4.13) легко найти для случая изотропной эффективной массы mαβ = m∗ δαβ . Тогда (4.13)
√
сводится к уравнению Шредингера для атома водорода с массой m∗ и зарядом e∗ = e/ κ и можно
получить решение (4.13), подставив m∗ и e∗ в известные из квантовой механики решения. Для спектра
получим
µ
¶
m∗ (e∗ )4
m∗ 1 me4
13.5 m∗ 1
εn = −
=−
=− 2
eV,
(4.14)
2~2 n2
m κ2 2~2 n2
n
m κ2
где n — главное квантовое число. Так как энергия отсчитывается от дна зоны проводимости, образуются
связанные состояния в запрещенной зоне. Типичные значения κ велики, а m∗ /m малы. Например, в Ge κ =
16, m∗ ≈ 0.2m0 . Поэтому энергии связи примесного атома на 3-4 порядка меньше ширины запрещенной
зоны. Такие примеси называют мелкими водородоподобными примесями.
Волновая функция в основном состоянии (n = 1) имеет вид:
−
1
f0 = √
e
πaB
|r − Rc |
aB ,
(4.15)
где R — узел, в котором сидит примесь, а
µ 2 ¶³
³m ´
~2
~
m´
−8
aB = ∗ ∗ 2 =
κ
=
0.53
×
10
κ cm,
m (e )
me2
m∗
m∗
где aB —¡ боровский
радиус, a = ~2 /me2 — боровский радиус в атоме M . Видно, что aB À a за счет
¢
фактора mm∗ κ примерно на 2 порядка. А так как волновая функция (4.15) экспоненциально спадает на
расстояниях больше aB , то получится, что она локализована на расстояниях много больше постоянной
решетки. Это означает, что приближение эффективной массы действительно применимо, поскольку возмущение является плавным по сравнению с межатомным расстоянием, а энергия мала по сравнению
с шириной запрещенной зоны. Это является также оправданным для использования макроскопической
диэлектрической проницаемости κ.
Аналогично мы могли бы рассуждать о спектре дырки у акцепторного атома, надо только отсчитывать
энергию в противоположном направлении и связанные состояния появятся в запрещенной зоне выше
максимума валентной зоны.
Мелкие примеси играют огромную роль в полупроводниках. Так как ширина запрещенной зоны EG
обычно велика (∼ 1 eV) по сравнению с kT (так как 300 К ∼ 0.025 eV), термические забросы электронов
из валентной зоны в зону проводимости из-за малой вероятности не приведут к заметной плотности
свободных электронов и дырок. А энергии ионизации примесных атомов (4.14) сравнимы с kT , поэтому
электроны или дырки связанных состояний могут с большой вероятностью переходить в непрерывный
спектр (то есть в зону проводимости или валентную зону), определяя проводимость кристалла.
На самом деле всё не так просто. Потенциал примеси не является чисто кулоновским и отличается
у разных атомов. Это приводит к так называемому потенциалу центральной ячейки и химическому
сдвигу локальных уровней. Химический сдвиг важен особенно для основного состояния, так как оно более
локализовано и волновая функция в возбужденных состояниях обращается в ноль в начале координат
22
Глава 4. Приближение эффективной массы.
(и не чувствует потенциал центра ячейки). Примеры энергии основного состояния трех доноров и трех
акцепторов (в электрон-Вольтах):
Примесь
Ge
Si
P
0,0120 0,044
As
0,0127 0,049
Sb
0,0096 0,039 .
B
0,0104 0,045
Al
0,0102 0,057
Ga
0,0108 0,065
До сих пор мы считали примесь изолированной. Если же имеется много примесей с концентрацией N ,
то среднее расстояние между ними r̄ ∼ 1/N 1/3 . Ясно, что примеси можно считать изолированными при
r̄ À aB , или при N a3B ¿ 1. В противном случае волновые функции начинают перекрываться и уровни
расщепляются, образуя внутри запрещенной зоны разрешённую примесную зону. Так как она образована
атомами, расположенными не периодически, её нельзя описывать блоховскими функциями и пользоваться
квазиимпульсом. В ней много локализованных состояний.
Ещё одно усложняющее обстоятельство — это анизотропия эффективной массы. В этом случае уже
нельзя свести решение к задаче об атоме водорода и изменяется даже систематика состояний, например,
квантовое число `, определяющее момент перестаёт быть "хорошим"квантовым числом. А если зона
вырождена в экстремуме, как это часто бывает в полупроводниках AIII B V , то надо учитывать сразу
две зоны и решать матричное уравнение Шредингера.
Ещё один вид отклонений от простого случая водородоподобного атома — это случай, когда число
валентных электронов примесного атома отличается больше чем на единицу. Тогда возникает несколько
различных уровней энергии и возникают глубокие уровни энергии, которые, строго говоря, не описываются
в методе эффективной массы. Например в Ge атом меди может существовать в четырех зарядовых
состояниях: Cu0 , Cu− , Cu2− , Cu3− .
4.4
Экситоны Ваннье-Мотта.
В рамках метода эффективной массы можно выйти за пределы приближения самосогласованного поля
и исследовать взаимодействие электрона и дырки, концентрация которых предполагается маленькой.
Электрон с дыркой могут образовывать связанное состояние, которое называется экситоном. Если размер
экситона велик по сравнению с межатомными расстояниями, он называется экситоном Ваннье - Мотта
и описывается в приближении эффективной массы. Противоположный случай — экситон Френкеля
— мы рассматривать не будем. Пусть rn и rp — радиус-векторы электрона и дырки. Тогда уравнение
Шредингера для них
·
¸
~2 2
~2 2
e2
O −
O −
Ψ(rn rp ) = EΨ(rn rp ).
(4.16)
−
2mn n 2mp p κ|rn − rp |
Введем координату центра тяжести и расстояния между электроном и дыркой
R=
mn r n + mp r p
,
mn + mp
ρ = rn − rp .
В новых координатах (4.16) принимает вид
¸
·
~2 2
e2
~2 2
O −
O −
Ψ(R, ρ) = EΨ(R, ρ),
−
2M R 2µ ρ κ|ρ|
−1
где M = mn + mp , µ−1 = m−1
n + mp .
23
(4.17)
Глава 4. Приближение эффективной массы.
В уравнении (4.17) разделяются переменные
Ψ(R, ρ) = χ(R)ϕ(ρ)
E = W + ε.
~2 2
O χ(R) = W χ(R),
−
2M R
· 2
¸
~ 2
e2
− Oρ −
ϕ(ρ) = εϕ(ρ),
2µ
κ|ρ|
(4.18)
(4.19)
Очевидно, что (4.18) описывает свободное трансляционное движение центра тяжести (то есть экситона
как целого) частицы с массой M и энергией W .(Решение — плоская волна, W = ~2 k 2 /2M ) Уравнение
(4.19) описывает водородоподобную структуру (структуру экситона), в которой роль массы играет приведёная
1 ³µ´ 1
1 ³ µ ´ 1 me4
≈− 2
·
масса µ. Собственные значения внутренней энергии экситона равны εn = − 2
n m κ2 2~2
n m κ2
2
m ~
m
13.5 эВ, а радиус равен κ 2 = κaB .
µ me
µ
Численные оценки для Ge дают εn ∼ 0.001 эВ, а размер экситона ∼ 42 Å, то есть много больше
межатомных расстояний. Таким образом, приближение эффективной массы применимо. В типичных
полупроводниках экситоны Ваннье-Мотта существуют, но только при очень низких температурах, при
температурах больше энергии экситонов последние ионизованы, то есть распадаются на свободные электрон
и дырку. В диэлектриках, где диэлектрическая проницаемость невелика (почему — будет видно из следующих
лекций), энергия экситона увеличивается, а радиус уменьшается и мы приходим к экситону Френкеля, к
которому неприменимо приближение эффективной массы.
24
Глава 5
Статистика носителей заряда в металлах и
полупроводниках.
5.1
Распределение Ферми-Дирака.
В приближении самосогласованного поля носителями заряда являются невзаимодействующие одноэлектроннные возбуждения — квазичастицы: электроны и дырки. Их часто называют электронным газом. Так
как спин электронов равен 1/2, эти квазичастицы описываются в равновесном случае (т.е. в отсутствие
диссипации, в частности, в отсутствие электрического тока) статистикой Ферми-Дирака. Пусть данное
квантовое состояние характеризуется номером зоны j, спином σ и квазиимпульсом p = (px , py , pz ). Тогда
среднее число электронов в этом состоянии определяется функцией распределения :
f0 (εjk ) =
1
εjk
−µ
exp(
)+1
kT
,
(5.1)
зависящей только от энергии состояния, здесь µ — энергия Ферми (синонимы: химический потенциал,
уровень Ферми).
Функция (5.1) экспоненциально приближается к 1 при ε − µ < kT и к 0 при ε − µ < kT . Т.е. все
состояния, находящиеся на расстояниях много больше kT от энергии Ферми, либо практически полностью
пусты, либо практически полностью заняты. Поэтому все самые существенные и интересные эффекты
связаны с состояниями, находящимися в области |ε − µ| . kT .
Формула (5.1) позволяет найти полное число электронов в объёме кристалла:
X
N=
f0 (εjk,σ ).
(5.2)
j,k,σ
В кристалле макроскопических размеров ( много больше межатомных размеров) спектр энергии квазинепрерывен и суммирование в (5.2) можно заменить интегрированием, воспользовавшись тем, что
V d3 k
плотность состояний в объёме фазового пространства d3 k есть
, где V — объём кристалла. Поэтому
(2π)3
для концентрации электронов n = N/V получим:
n=
XZ
j,σ
d3 k
f0 (εjk,σ ) =
(2π)3
Z∞
dεN (ε)f0 (ε),
0
25
(5.3)
Глава 5. Статистика носителей заряда в металлах и полупроводниках.
где N (ε) — плотность
в интервале dε. Для случая, когда ε = p2 /2m∗ = ~2 k 2 /2m∗ из (5.3)
√ ∗3/2состояний
√
2m
ε
следует N (ε) =
и
2
3
2π ~
√
√
Z∞
2(m∗ kT )3/2
dx x
n=
,
~3
1 + ex−z
(5.4)
0
µ
ε
и в интеграле сделана подстановка x =
.
kT
kT
Интеграл (5.4) может быть вычислен в элементарных функциях в предельных случаях z À 1 и z ¿ 0.
где z =
5.2
Вырожденный электронный газ. Металл.
Вырожденным называется газ, в котором µ À kT,
(z À 1) . В этом случае вид функции Ферми похож
1
→ θ(z − x) и
на ступеньку, а её производной на δ - функцию. В первом приближении в (5.4)
1 + ex−z
√
√
Zz
2(m∗ kT )3/2 √
2 2 (m∗ µ)3/2
n=
xdx =
,
~3
3
~3
(5.5a)
0
или
µ(T = 0) =
~2
(3π 2 n)2/3 .
2m∗
(5.5b)
Следующей поправка квадратична по 1/z к (5.5b) и с ее учетом химический потенциал приобретает вид
"
µ ¶ #
π 2 kT 2
µ = µ0 1 −
,
(5.6)
12 µ
где µ0 = µ(0).
Условие вырождения можно записать в виде
~2 (3π 2 n)2/3
µ0
=
À 1.
kT
2m∗
kT
(5.7)
Видно, что вырождению способствует большая концентрация электронов. В металлах, где n ∼ 1022 см−3 ,
для m∗ ∼ 10−27 г получим µ0 ∼ 104 ÷ 105 К (или µ0 ∼ 1 ÷ 10 эВ). В первом случае энергия измерена
в градусах, т.е. фактически приведена не энергия µ0 , а µ0 /k. Такая большая величина µ0 означает,
что в металлах электроны вырождены при всех температурах, вплоть до температуры плавления. В
полупроводниках электроны тоже могут быть вырождены при низких температурах, но электронов в
полупроводнике значительно меньше, чем в металлах, поэтому для того, чтобы получить вырожденный
полупроводник, требуется большая концентрация мелких примесей в полупроводнике.
5.3
Невырожденный электронный газ.
Этот случай соответствует статистике Больцмана |z| À 1и z < 0.
¶
µ
µ−ε
.
f0 ' exp
kT
26
(µ < 0).
Глава 5. Статистика носителей заряда в металлах и полупроводниках.
Интетрируя (5.4) в этом приближении, получим:
√
Z∞
√
2 (m∗ kT )3/2 µ
n= 2
dx xe−x
e kT
3
π
~
0
или
³ µ ´
³ µ ´
(2πm∗ kT )3/2
exp
≡
N
exp
,
C
4π 3 ~3
kT
kT
·
¸
4π 3 ~3 n
n
µ = kT ln
≡ kT ln
.
3/2
∗
NC
(2πm kT )
n=
(5.8)
(5.9)
Критерий применимости классической статистики:
(2πm∗ kT )3/2
À 1.
4π 3 ~3 n
Видно, что невырожденности газа способствует малая концентрация n, высокая температура и большая
эффективная масса m∗ . Для n = 1017 см−3 , m∗ = 10−27 г и T ≈ 300 К получим exp(−µ/kT ) = 300 À 1,
газ не вырожден.
Оценим величину NC , которая называется эффективной плотностью состояний в зоне проводимости
(не путать с плотностью электронных состояний):
µ
NC = (2.5)
m∗
m
¶3/2 µ
T
300K
¶3/2
· 1019 cm−3 .
(5.10)
Формула (5.8) не даёт в явном виде концентрацию электронов в полупроводнике, она лишь выражает
n через неизвестную величину µ. Для явного вычисления n, µ и их температурных зависимостей надо
воспользоваться условием нейтральности, включив в рассмотрение электроны и дырки как в зонах, так
и на примесных уровнях. При этом снова заметим, что дырки можно описывать как ферми-частицы,
энергия которых отсчитывается в отрицательном направлении. Вероятность того, что состояние не занято
электроном (т.е. занято дыркой) есть
f00 (ε) = 1 − f0 (ε) =
1
1
=
,
0
0
µ−ε
ε −µ
exp
+1
exp
+1
kT
kT
(5.11)
0
где µ0 = εG − µ, ε = −ε − εG — химический потенциал и энергия, соответственно, отсчитанные от потолка
валентной зоны в отрицательном направлении ( εG — ширина запрещённой зоны). Видно, что в таком
0
0
виде функция распределения дырок f0 (ε ) описывается распределением Ферми. Формула (5.11) позволяет
получать все результаты для дырок заменой обозначений в формулах, полученных для электронов.
Условие нейтральности запишем в виде:
+
,
n + NA− = p + ND
(5.12)
где n и p — концентрации электронов и дырок в зонах, NA− — концентрация ионизованных акцепторов
+
(т.е. захвативших электрон), ND
— концентрация ионизованных доноров (отдавших электрон в зону
проводимости).
Концентрации n и p задаются формулой (5.8):
µ
¶
³ µ ´
−εG − µ
;
p = NV exp
.
(5.13)
n = NC exp
kT
kT
27
Глава 5. Статистика носителей заряда в металлах и полупроводниках.
Найдём теперь плотность донорных состояний, с которых ушли электроны: Nd+ = Nd (1−nd ). Можно было
εD − µ
бы предположить, что nd описывается функцией распределения Ферми nd = 1/(1 + exp
), однако
kT
это было бы неверно , так как распределение Ферми применимо для невзаимодействующих электронов, а
примесные состояния нельзя рассматривать без учёта кулоновского взаимодействия. Всего с примесным
уровнем связано 4 состояния: 0 электронов, 1 электрон со спином вверх, 1 электрон со спином вниз и 2
электрона с противоположными спинами. Однако последнее состояние имело бы огромную энергию из-за
кулоновского
отталкивания
электронов, помещённых на один атом, поэтому вероятность его заполнения
µ
¶
−Ej − µNj
∝ exp
¿ 1 (здесь Nj — число частиц в состоянии j, а Ej — энергия состояния).
kT
Согласно принципам статистической физики, среднее число частиц определяется как
P
−(Ej − µNj )
−(εD − µ)
2 exp
kT
j
kT
nD =
=
P
−(Ej − µNj )
−(εD − µ)
1 + 2 exp
exp
kT
kT
j
Nj exp
Следовательно,
ND
µ − εD
1 + 2 exp
kT
Аналогично можно получить плотность ионизованных акцепторов
+
ND
=
NA+ =
NA
εA − µ
1 + 2 exp
kT
(5.14)
(5.15)
Подставляя (5.13 - 5.15) в (5.12) получим уравнение для величины x = exp(µ/(kT )):
µ
¶
Na
−εG 1
N
µ D¶
.
NC x +
= NV exp
+
³ε ´ 1
|εD |
kT
x
A
2 exp
+1
2 exp
x+1
kT x
kT
(5.16)
Уравнение (5.16) — уравнение 4-ой степени относительно x. В общем случае его решение слишком громоздко,
поэтому мы рассмотрим несколько предельных случаев.
1. Собственный полупроводник. Этот случай соответствует отсутствию примесей или достаточно
высоким температурам, когда свободных носителей заряда так много, что примесями можно пренебречь в
силу условия NC , NV , À NA , ND , выполняющегося в достаточно чистых полупроводниках при достаточно
|ND − NA |
высоких температурах. Тогда при x À √
вторыми слагаемыми в обеих частях (5.16) можно
NC NV
пренебречь. В результате получим:
r
³ ε ´
³ ε ´ µ m ¶3/4
NV
p
G
G
x=
=
exp −
;
exp −
NC
2kT
mn
2kT
µ
¶
p
mp
2εG
εG 3
p = n = NC NV exp −
; µ=−
+ kT ln
(5.17)
kT
2
4
mn
2. Рассмотрим теперь, что происходит при понижении температуры. Предположим, что ND À NA
(нескомпенсированный полупроводник донорного типа). В этом случае можно пренебречь влиянием акцепторов и дырок. Тогда (5.16) сводится к квадратному уравнению
µ
¶
|εD |
NC x 2xe kT + 1 = ND ,
28
Глава 5. Статистика носителей заряда в металлах и полупроводниках.
или
1 |εD |
ND − |εD |
e kT = 0
x2 + e− kT x −
2
2NC
s
1 − |εD |
1 − 2|εD |
ND − |εD |
x = − e kT +
e kT +
e kT .
(5.18)
4
16
2NC
µ
¶
Nc
|εD |
2а. При
exp −
À 1 (что возможно при не слишком низких температурах), первый член
ND
kT
под корнем в (5.18) доминирует и
µ
¶
1 − |εD | 1 − |εD |
ND − |εD |
ND
kT
kT
kT
x≈− e
+ e
1+4
e
=
4
4
NC
NC
n = xNC = ND ,
µ = kT ln
ND
.
NC
(5.19)
Согласно (5.19) практически все доноры при таких относительно высоких температурах ионизованны
(n = ND ) и химпотенциал
µ
¶ находится ниже донорного уровня.
|εD |
NC
exp −
¿ 1, что возможно при более низких температурах. Под корнем в (5.18) доминирует
2б.
ND
kT
2-ое слагаемое.
r
r
µ
¶
µ
¶
|εD |
|εD |
ND
NC ND
exp −
;
n = xNC =
exp −
(5.20)
x≈
2NC
2kT
2
2kT
µ=−
|εD | kT
NC ND
+
ln √ .
2
2
2
Таким образом, концентрация электронов уменьшается при понижении температуры, а уровень Ферми
лежит посередине между донорным уровнем и зоной проводимости.
2в. Отметим, однако, что в реальных материалах всегда имеется какое-то количество компенсирующих
примесей (то есть NA 6= 0 ), которые расположены очень низко (при ε ∼ −εG ) и, поэтому должны быть
заполнены электронами. Это означает, что часть электронов перейдёт на них с доноров. т.е. часть доноров
окажется ионизованной. Т.е. при T → 0 уровень Ферми будет проходить по границе пустых и заполненных
состояний, находящейся на донорных уровнях. В результате при самых низких температурах
¶
µ
|εD |
.
(5.21)
n ∝ exp −
kT
µ
¶
Ea
Зависимости вида n ∝ exp −
, как в (5.17), (5.20) и (5.21) называются термоактивационными (или
kT
законом Аррениуса), при этом Ea называется энергией активации. Так как активационные зависимости
соответствуют изменениям n на многие порядки величины, зависимости n(T ) обычно изображают в
координатах (ln n, 1/T ). Найденные зависимости n(T ) в таких координатах показаны на рисю 5.1. Из
рисунка видно, что концентрация электронов уменьшается при понижении температуры. Из-за этого
растёт сопротивление полупроводников при понижении температуры. Ситуацию в случае легирования
акцепторами можно понять, заменив nD , εD и mn в полученных формулах на nA , εA и mp . Отметим, что
согласно (5.8) и (5.17) в невырожденном полупроводнике всегда
µ
¶
√
(2π mp mp kT )3
−εG
np = n2i =
exp
,
(5.22)
(4π 3 ~3 )2
kT
29
Глава 5. Статистика носителей заряда в металлах и полупроводниках.
Рис. 5.1: Температурная зависимость концентрации электронов в полупроводнике. 1) собственная область,
энергия активации εG ; 2a) концентрация не зависит от T ; 2б) энергия активации εG /2; 2в) энергия
активации εG .
то есть произведение концентраций электронов и дырок для данного материала не зависит от наличия
примесей. При вычислении n и p в (5.8) мы предполагали изотропный параболический закон дисперсии
электронов. В других случаях выражения для плотностей NC и NV , а значит и предэкспонента в (5.22)
изменятся. Если, например, имеется несколько эквивалентных минимумов энергии νc и эффективная
масса анизотропна, то мы получили бы для такое же выражение как в (5.8), если под m∗ понимать
2/3
"эффективную массу плотности состояний"mef f = νc (m1 m2 m3 )1/3 . Аналогично для NV .
30
Глава 6
Диэлектрическая проницаемость твердого
тела.Формула Линдхарда.
6.1
Пространственная и временная дисперсия.
При приложении к твердому телу внешнего электрического поля появляется электрический ток j или
1 ∂(D − E)
вектор смещения D. Эти величины связаны между собой соотношением j =
, следующим
4π
∂t
из уравнений Максвелла, поэтому одни и те же рассуждения можно применять и к току, и к смещению.
В линейном приближении по полю связь между током (или вектором смещения) и приложенным электрическим полем называется линейным откликом. В случае однородной и стационарной среды линейный
отклик в общем виде может быть записан как
Z
Dα (r, t) =
Zt
3
0
0
d r
0
0
dt εαβ (r − r , t − t )Eβ (r, t),
(6.1)
−∞
то есть D ( или j) определяется не только локальным в пространстве и времени полем, но и полем в
предыдущие моменты времени и полем в соседних точках. Такая нелокальная в пространстве и/или во
времени зависимость называется, соответственно, пространственной и временной дисперсией.
Соотношение (6.1) приобретает более простой и привычный вид, если к нему применить преобразование
Фурье по координатам и по времени:
Z
d3 k dω i(kr−ωt)
e
Dα (k, ω),
(6.2)
Dα (k, ω) = εαβ (k, ω)Eβ (k, ω), Dα (r, t) =
(2π)3 2π
где Eβ (k, ω) и εαβ (k, ω) также связаны с εαβ (r, t) и Eβ (r, t) формулой обращения Фурье. Зависимость ε от
k и ω отражает наличие, соответственно, пространственной и временной дисперсии. Задание комплексного
ε эквивалентно заданию комплексной проводимости σ:
εαβ = δαβ + i
4π
σαβ .
ω
(6.3)
Так что можно описывать линейный отклик твердого тела на электрическое поле либо на языке комплексной
диэлектрической проницаемости, либо на эквивалентном языке комплексной проводимости.
31
Глава 6. Диэлектрическая проницаемость твердого тела.Формула Линдхарда.
6.2
Вычисление диэлектрической проницаемости с помощью теории возмущений.
Ниже мы найдем отклик системы электронов в твердом теле на электрическое поле, то есть, какие
возникают возмущения плотности заряда δρ и потенциала U при приложении стороннего возмущающего
потенциала V . Внешний потенциал вызовет перераспределение электронов с плотностью δρ, которая, в
свою очередь, создаст дополнительный электрический потенциал Φ. При этом электроны будут чувствовать полный потенциал
U = Φ + V.
(6.4)
При вычислении линейного отклика, то есть отклика в линейном приближении по электрическому
полю, можно воспользоваться преобразованием Фурье. В результате мы вычислим отклик на гармоническое электрическое поле. А отклик на возмущение с произвольной зависимостью от r и t может быть
вычислен как сумма откликов на гармоники exp(ikr − iωt), на которые он разложен по Фурье. Итак,
ищем отклик на внешний потенциал Vq exp(iqr + iωt + 0t), где i0 – бесконечно малая добавка к частоте,
учитывающая, что система “теряет память” о событиях произошедших бесконечно давно, и позволяющая
в последующих расчетах правильно выполнять интегрирование в формально расходящихся интегралах.
Пусть в отсутствие внешнего потенциала V электроны описываются полным набором блоховских
волновых функций Ψnk ≡ |ki в состояниях с энергией Ek (в схеме расширенных зон можно не писать
номер зоны, поскольку номера различных зон можно задавать различными векторами обратной решетки).
Из теории возмущений для потенциала, периодически зависящего от времени поправка к блоховской
волновой функции |ki под действием возмущения принимает вид:
Ψ(r, t) = |ki + eUq
X
k
0
0
hk |eiqr |ki
0
|k i,
Ek − Ek0 − ~(ω − i0)
(6.5)
где используется тот факт, что электроны чувствуют полный потенциал U .
Возмущение плотности заряда δρ = ρ − ρ(V = 0)
X
f (k)[|Ψ(r, t)|2 − 1].
δρ = −e
k
В линейном приближении по Uq для фурье-компоненты, соответствующей q, получим:


¯
¯2
¯ 0
¯
f (k) ¯hk |eiqr |ki¯
X

2
δρq = −e
− 1 + o(Uq ) .
|hk|ki| + 2Uq e
Ek − Ek0 − ~(ω − i0)
0
(6.6)
k,k
Обратите внимание, что вследствие вычисления фурье-компоненты в (6.6) появился квадрат матричного
элемента.
0
0
В матричном элементе при суммировании по k сохранятся лишь компоненты с k = k + q + G
¯
¯
X f (k) ¯hk + q + G|eiqr |ki¯2
2
.
(6.7)
δρq = −2e Uq
Ek − Ek+q+G − ~(ω − i0)
k,G
Воспользуемся теперь свойством преобразования Фурье для действительной величины ρ(q, ω) = ρ∗ (−q, −ω).
Получим:
¯
¯2
¯
¯
¯
¯
X f (k) ¯hk − q + G|eiqr |ki¯2
X f (k̃ + q + G) ¯hk̃|eiqr |k̃ − q + Gi¯
2
2
δρq = −2e Uq
= −2e Uq
,
Ek − Ek−q+G + ~(ω − i0)
Ek̃+q+G − Ek̃ + ~(ω − i0)
k,G
k̃,G
32
Глава 6. Диэлектрическая проницаемость твердого тела.Формула Линдхарда.
(6.8)
где во втором равенстве (6.8) мы сделали замену переменных k − q + G = k̃; G → −G Возьмем в качестве
δρq полусумму уравнений (6.7) и (6.8):
2
δρq = −e Uq
¯
¯
X [f (k) − f (k + q + G)] ¯hk|eiqr |k + q + Gi¯2
Ek − Ek+q+G − ~(ω − i0)
k,G
.
(6.9)
Подставим (6.9) в уравнение Пуассона −O2 Φ = 4πδρ (которое после преобразования Фурье имеет вид
q = 4πδρq ) и подставим из (6.4) Φq в виде Φq = Uq − Vq . В результате получим
q2Φ
¯2
¯
4πe2 X [f (k) − f (k + q + G)] ¯hk|eiqr |k + q + Gi¯
Uq = Vq + Uq 2
.
q
Ek − Ek+q+G − ~(ω − i0)
(6.10)
k,G
По определению диэлектрической проницаемости она описывает ослабление потенциала средой Uq =
Vq /ε(q, ω). Tаким образом из (6.10) следует
¯2
¯
4πe2 X [f (k) − f (k + q + G)] ¯hk|eiqr |k + q + Gi¯
.
ε(q, ω) = 1 + 2
q
Ek+q+G − Ek + ~(ω − i0)
(6.11)
k,G
Это и есть формула Линдхарда для диэлектрической проницаемости. Видно, что вклад в диэлектрическую
проницаемость дают виртуальные переходы электронов между состояниями |ki и |k + q + Gi.
Используем теперь формулу Линдхарда для вывода нескольких фундаментальных свойств металлов
и полупроводников.
6.3
Экранирование статического (ω = 0) поля в проводниках.
Рассмотрим сначала диэлектрическую проницаемость ε(0, q) в длинноволновом пределе q → 0, т.е. отклик
на плавный потенциал. В пределе q → 0 можно упростить формулу Линдхарда (6.11). Эта формула
содержит вклады, относящиеся к разным зонам (G 6= 0) и от состояний в одной зоне G = 0). Причем,
вклад от состояний (|ki и |k + qi, относящихся к одной зоне, возможен только, если f (k) 6= f (k + q), то
есть если в зоне имеются как пустые, так и занятые состояния. Рассмотрим вклад в ε от состояний одной
зоны, иными словами, вклад электронов проводимости. Обозначим сумму слагаемых с G (а она в пределе
q → 0, как мы сможем убедиться позднее, стремится к константе) κ − 1. Разложим теперь все величины
в (6.11) в ряд по малым q с точностью до первых неисчезающих членов:
E(k + q) − E(k) ≈ q∇k E(k); f (k) − f (k + q) ≈ −q
∂f
∇k E(k); hk|eiqr |k + qi ≈ hk|ki;
∂ε
Получим:
4πe2
ε(0, q) w κ + 2
q
¶
¶
µ
¶
Z µ
Z µ
4πe2
λ2
∂f
∂f
3
d k=κ+ 2
N (E)dE ≡ κ 1 + 2 ,
−
−
∂E
q
∂E
q
(6.12)
где введено обозначение λ2 для интеграла с сомножителем 4πe2 . Подобный интеграл мы вычисляли
∂f
в предыдущей лекции. В металле, где f близка к ступеньке (вырожденный газ)
∼ δ(E − EF ) и
∂E
33
Глава 6. Диэлектрическая проницаемость твердого тела.Формула Линдхарда.
4πe2
~2 k 2
=
N (EF ), что при спектре E =
дает плотность состояний N (E) =
κ
2m
4
2 · ( 3 πkF3 )
что n =
, мы получим
(2π)3
λ2
λ2 =
4πe2 3 n
4me2 n1/3
≡
.
κ 2 EF
~2 κ
(6.13a)
Вычисление интеграла в (6.12) в общем случае дает
¯
√
√
¯
Z∞ √
x−z
∗3/2
¯
4 2 2m
kB T
xdxe
2
¯
λ =
.
3
x−z
2
π
κ~
(1 + e ) ¯¯
EF
0
√
√
2(m∗ )3/2 E
, учитывая,
2π 2 ~3
(6.13b)
z= k
BT
Для невырожденного электронного газа (полупроводники), когда в (6.13b) единица отбрасывается по
сравнению с экспонентой,
λ2 =
4πe2 n
.
κkB T
(6.13c)
Из (6.12) видно, что чем более плавное поле (то есть содержащее гармоники с меньшими q), тем больше ε
и, соответственно, сильнее ослабление поля. Если в качестве внешнего поля выступает кулоновское поле
e2
4πe2
точечного заряда V (r) = , Vq = 2 , полное поле
r
q
Uq =
Vq
=
ε(0, q)
4πe2
e2 −λr
4πe2
µ
¶ =
⇒
U
(r)
=
e ,
2
κ(q 2 + λ2 )
κr
λ
κ 1 + 2 q2
q
(6.14)
т.е. кулоновский потенциал превращается в быстро, экспоненциально, спадающий потенциал Юкавы.
Длина спадания потенциала r0 = 1/λ называется радиусом экранирования. Такое экранирование в невырожденном газе называют дебаевским экранированием, в вырожденном — экранированием Томаса Ферми. В типичных металлах r0 порядка межатомного расстояния.
Экранирование (6.13 - 6.14) можно получить и из более наглядных соображений: Решим уравнение
Пуассона:
div(κE) = 4π(ρ + ρext ),
(6.15)
где ρext – сторонний заряд, а ρ — плотность заряда, созданного из-за перераспределения электронов
проводимости, возникшего из-за появления электрического потенциала Φ в проводнике
Z
ρ = −e[n(r) − n0 ] = −e [f (E + eΦ) − f (E)]N (E)dE,
(6.16)
здесь n(r) — плотность электронов, а n0 — невозмущенная плотность электронов (при Φ = 0). Поскольку
речь идет о равновесной системе (нет тока), функция распределения остается фермиевской, но энергия
изменяется на величину потенциальной энергии eΦ. Разложим (6.16) в ряд по Φ. Получим:
Z
∂f
λ2 κ
2
ρ = −e Φ
N (E)dE ≡ −
Φ.
(6.17)
∂E
4π
Подставив (6.17) в (6.15), получим уравнение для Φ (E = −OΦ):
O2 Φ − λ2 Φ = −
4π
ρext .
κ
(6.18)
34
Глава 6. Диэлектрическая проницаемость твердого тела.Формула Линдхарда.
Если в качестве ρext выступает точечный заряд ρext = eδ(r), решение легко найти, сделав преобразование
Фурье уравнения (6.18) (учитывая, что преобразование Фурье от δ-функции — единица):
−q 2 Φ − λ2 Φ = −
4πe
4πe
Φ=
.
2
κ
κ(q + λ2 )
Отсюда в координатном представлении получаем
Φ(r) =
e2 −λr
e ,
κr
что совпадает с (6.14).
Рассмотрим теперь случай экранирования постоянного внешнего поля E0 , перпендикулярного поверхности проводника (например, проводник помещен в конденсатор). Пусть проводник занимает полупространство x > 0. Тогда нам надо решить (6.18) с граничными условиями E(x → ∞) = 0 и E(x = 0) = E0 .
Решение имеет вид
E = E0 e−λx ,
то есть поле проникает вглубь проводника только на длину экранирования. Так как поле проникает в
металл только на расстояния в несколько ангстрем, то говорят, что оно не проникает в металл, а только
создает поверхностный заряд.
Эти рассуждения наглядно демонстрируют, что мы учитываем кулоновское взаимодействие электронов
в приближении самосогласованного поля eΦ. А что, если не ограничиваться малыми q, то есть учесть
изменение U (r) на малых расстояниях? Честное вычисление ε(0, q) по формуле Линдхарда (6.11) для газа
свободных электронов при T = 0 (модель Зоммерфельда для металла) дает
¯
¯¾
½
1 4kF2 − q 2 ¯¯ 2kF + q ¯¯
4πe2 n
+
ln ¯
ε(q) = 1 + 2 2
(6.19)
q 3 εF 2
8kF q
2kF − q ¯
Зависимость ε(q) имеет особенность при |q| = 2kF , которая называется коновской аномалией. Эта особенность приводит к сингулярности в фононных частотах ωph (q) (эффект Кона). Особенность возникает по
следующим причинам: ионы при колебаниях создают потенциал Vq , который экранируется электронами.
Поэтому ионы взаимодействуют между собой посредством экранированного потенциала Uq = V (q)/ε(q),
имеющего особенность при |q| = 2kF . В результате собственные частоты колебаний ионов (фононные
частоты) зависят от ε и на них появляется коновская аномалия.
Наличие особенности в ε при |q| = 2kF приводит к изменению характера спадания поля U . Дополнительно к плавной зависимости вида (6.14) появляется добавка, осциллирующая с волновым числом 2kF ,
т.е. индуцированный потенциал и индуцированная плотность заряда на больших расстояниях спадают
R iqr 4πe2 d3 q
cos 2kF r
как
.
В
этом
можно
убедиться,
проинтегрировав
U
(r)
=
e
с ε(q) в виде (6.19).
r3
q 2 ε(q) (2π)3
6.4
Низкочастотная диэлектрическая проницаемость диэлектриков.
Рассмотрим вклад в ε от состояний из разных зон при q → 0. Такой вклад возможен, если f (k) 6=
f (k + q + G), т.е. если одна зона заполнена, а вторая пустая. В этом случае в знаменателе (6.11) стоит
величина ε(k + q + G) − ε(k) & εG .
Нам надо еще оценить числитель. При q → 0 и G 6= 0 матричный элемент в числителе формулы
Линдхарда (6.11) обращается в 0, что приведет к сокращению q 2 в знаменателе, так что вклад в диэлектрическую проницаемость за счет виртуальных переходов между различными зонами стремится при
35
Глава 6. Диэлектрическая проницаемость твердого тела.Формула Линдхарда.
q → 0 к константе. Для того, чтобы убедиться в этом более строгим образом, воспользуемся обобщенным
"правилом сумм для сил осцилляторов справедливым для любого полного набора волновых функций |ni
~2 O2
оператора Ĥ = −
+ U (r):
2m
X
¯
¯2 ~2 q 2
(En − Es ) ¯hn|eiqr |si¯ =
.
2m
n
(6.20)
(Здесь m — масса свободного электрона !) Правила такого рода используются очень часто; мы использовали
его разновидность, когда получали связь эффективной массы m∗ в полупроводниках с m. Докажем (6.20).
Для этого вычислим двумя способами двойной коммутатор:
h
i
[Ĥ, eiqr ], e−iqr ≡ 2Ĥ − eiqr Ĥe−iqr − e−iqr Ĥeiqr .
(6.21)
1) С одной стороны, вычислим действие коммутатора на некую ψ, непосредственно учитывая, что
O2 [e±iqr ψ(r)] = [−q 2 ± 2iqO + O2 ]ψ,
а гамильтониан имеет вид Ĥ = −
~2 O2
+ U (r):
2m
h
i
~2 q 2
2q 2 ~2
+ 2iqO − 2iqO] = −
.
(6.22)
[Ĥ, eiqr ], e−iqr = [2Ĥ − 2Ĥ −
2m
m
h
i
2)С другой стороны, вычислим матричный элемент (6.21) hs| [Ĥ, eiqr ], e−iqr |si, используя представление
P
единичного оператора в виде Ê = |nihn| ≡ 1. (Последнее тождество легко доказывается1 ). Итак,
n
hs|2Ĥ − eiqr Ĥe−iqr − e−iqr Ĥeiqr |si =
= 2Es −
i
Xh
hs|eiqr |nihn|Ĥ|mihm|e−iqr |si + hs|e−iqr |nihn|Ĥ|mihm|eiqr |si =
m,n
= 2(Es −
X
|hs|eiqr |ni|2 En ) = 2
n
P
X
(Es − Em )|hs|eiqr |ni|2 .
n
P
hs|eiqr |nihn|e−iqr |si = hs|eiqr e−iqr |si = hs|si = 1.
¿ ¯ 2 2¯ À
¯~ q ¯
¯ s (6.22), мы
Сравнивая теперь (6.23) и матричный элемент того же коммутатора в виде − s ¯¯
m ¯
получим правило сумм (6.20).
Домножив числитель и знаменатель (6.11) на E(k + q) − E(k) − ~ω и используя опять симметрию
выражения под суммой относительно комплексного сопряжения и взаимозамены переменных суммирования k → k + q + G (как это уже использовали в (6.8 - 6.10)), можно привести (6.11) к виду
¯
¯2
4πe2 X [f (k) − f (k + q + G)][E(k) − E(k + q + G)] ¯hk|eiqr |k + q + Gi¯
ε(q, ω) = 1 + 2
. (6.24)
q
~2 ω 2 − [E(k) − E(k + q + G)]2
В последнем равенстве учтено, что
|hs|eiqr |ni|2 ≡
(6.23)
n
n
k,G
Положив теперь f (k) = 1, f (k + q + G) = 0 (заполненная и пустая зоны), сделаем грубую оценку для ω →
0. А именно: заменим в знаменателе E(k + q + G) − E(k) ≈ EG (что является неплохим приближением
1
действительно, для любой f (r), f (r) =
P
fm |mi, Êf (r) =
P
m
n
f (r), где hn|mi = δnm — в силу ортонормированности
36
|nihn|f (r) =
P
n,m
|nihn|mifm =
P
n,m
fm |niδnm =
P
m
fm |mi =
Глава 6. Диэлектрическая проницаемость твердого тела.Формула Линдхарда.
при больших запрещенных и маленьких разрешенных зонах). Применим теперь (6.20) к полному набору
волновых функций |k + Gi. Тогда
X
¯
¯2 ~2 q 2
[E(k) − E(k + q + G)] ¯hk|eiqr |k + q + Gi¯ ≈
.
2m
G
Выполняя теперь в (6.24) при ω = 0 суммирование по k, которое сводится к умножению на число состояний
в зоне, равное числу электронов n, получим
µ
¶
~ωp 2
4πe2 ~2 q 2
ε(0, q) ≈ 1 + 2n 2
=1+
,
(6.25)
q 2mEG
EG
µ
¶1/2
4πne2
где введено обозначение ωp =
— плазменная частота, смысл которой выяснится ниже. Из
m
оценочной формулы (6.25) следует, что в полупроводниках с узкой запрещенной зоной следует ожидать
большую диэлектрическую проницаемость, а в изоляторах с широкой EG диэлектрическая проницаемость
должна быть меньше. Именно эти ε из (6.25) фигурируют в качестве κ в (6.12) и (6.15).
6.5
Экранирование на высоких частотах.(q → 0, ω - велика)
Во-первых, отметим, что ε(ω, q) имеет сингулярность и мнимую часть
1
1
=P
+iπδ[E(k + q + G)−E(k)+~ω]
E(k + q + G) − E(k) + ~ω − i0
E(k + q + G) − E(k) + ~ω
— значит на частотах ω = ([E(k + q + G) − E(k)])/~ имеется поглощение (диссипация), ведь мнимая
часть ε связана с действительной частью проводимости σ, описывающей поглощение. Такое поглощение
возможно за счет переходов между разными зонами. А так как скорость света очень велика, то на частотах
~ω ∼ EG , волновой вектор света q ¿ G, поэтому говорят о вертикальных междузонных переходах из E(k)
в E(k + G). (Напомним, что мы работаем в схеме расширенных зон, где k + G обозначает состояние с
тем же квазиимпульсом k, но в другой зоне.)
Рассмотрим теперь случай более высоких частот, больше энергетической разности в (6.11) или (6.24).
При этом в качестве знаменателя можно взять просто ω 2 , и снова произвести суммирование по G с
помощью (6.20) как при получении (6.25), полагая в числителе f (k) = 1 и f (k + G + q) = 0. Выполняя
суммирование по k, как при выводе (6.25), получим
ε(q, ω)|q→0 = 1 −
ωp2
.
ω2
(6.26)
Такую же формулу можно получить для электронов, находящихся в одной зоне, когда справедлива модель
Зоммерфельда для газа свободных электронов. В этом случае матричный элемент равен 1 и, разлагая
числитель в ряд по q и выполняя интегрирование по k по частям, получим
µ
¶
Z
ωp2
4πe2
4πne2
d3 k df (k) dE 2 q 2
ε(q, ω) = 1 + 2
=
1
−
=
1
−
.
(6.27)
q
(2π)3 dE
dk
(~ω)2
mω 2
ω2
Из (6.26) следует, что ε = 0 при ω = ωp /. Так как divD = div(εE) = 0 в этом случае можно сделать
вывод, что даже в отсутствие внешнего поля на частоте ωp возможно ненулевое электрическое поле, так
как оно удовлетворяет уравнениям Максвелла без внешних зарядов и полей. То есть на частоте ω = ωp
возможны собственные колебания продольного (потенциального, rotE = 0,) электрического поля. Такие
37
Глава 6. Диэлектрическая проницаемость твердого тела.Формула Линдхарда.
колебания называются плазменными. Этот результат можно проиллюстрировать с помощью простой
модели. При смещении электронов возникает дипольный момент и поляризация P = nex, а значит и
электрическое поле E = −4πP. Из уравнений движения для электрона
mẍ = eE = −4πne2 x
следует, что x = x0 e−iωp t . Плазменная частота в металлах составляет огромную величину ∼ 10 eV À kB T .
Плазменные колебания возбуждаются при прохождении быстрых электронов через пленки. Частота в 5 −
10 эВ превышает типичные EG для полупроводников, поэтому и справедливы рассуждения, приведенные
к (6.24), при этом в качестве n в ωp входит концентрация всех электронов валентной зоны.
А в полупроводниках возможны и плазменные колебания с малой частотой, в которых участвуют
только электроны проводимости. Для них плазменные колебания получаются путем рассуждений, приведших к (6.27). При этом в качестве n войдет концентрация электронов в зоне проводимости, а вместо
m — эффективная масса электронов проводимости. Такие же рассуждения можно провести и для дырок.
Плазменная частота может быть определена и в оптических экспериментах. Дело в том, что при
√
ω < ωp , ε < 0 и коэффициент преломления n = ε — мнимый, это соответствует полному отражению
света при ω > ωp , ε > 0 и n — действительно, т.е. проводник становится прозрачным. И действительно
металлы становятся прозрачными в ультрафиолетовом свете.
Плазменные колебания можно рассматривать как квазичастицы - плазмоны, они ведь такие же слабозатухающие собственные возбуждения (нормальные колебания, нормальные моды ), как фононы, электроны и дырки. Закон дисперсии плазмонов мы вычисляли в пределе q → 0, а если бы мы учли поправки
следующего порядка по q, то получили бы для спектра плазмонов в вырожденном электронном газе
q2v2
ωp (q) = ωp (0)[1+a 2F +. . .], (vF = pF /m), где a – постоянная порядка единицы. В случае параболического
ωp
закона дисперсии электронов a = 3/10.
6.6
Переход Мотта-Хаббарда.
Когда мы обсуждали образование энергетической зонной структуры в рамках модели сильно связанных
электронов (уширение уровней при сближении атомов), мы совсем не учитывали в наших рассуждениях
короткодействующее кулоновское взаимодействие. Учет кулоновского взаимодействия приводит к возможности перехода металл - диэлектрик. Действительно, пусть для помещения второго электрона на атом
требуется преодолеть энергию кулоновского отталкивания U . Если эта энергия отталкивания U À t, где t
— ширина разрешенной зоны, вызванная перекрытием волновых функций, то электрон не сможет перейти
на следующий атом (и участвовать в проводимости). То есть получится диэлектрик. Если же t À U , то
справедливы наши прежние рассуждения, что при нечетном числе электронов на атоме образуется металл.
Хаббард показал, что переход между состояниями металла и диэлектрика происходит при U = 1.15t.
При больших значениях U зона расширяется на 2 (пустую и заполненную), разделенные энергетической
щелью. Так как величина t увеличивается при уменьшении межатомного расстояния a, следует ожидать,
что при малых a должно быть металлическое состояние, а при больших - диэлектрическое (моттовский
диэлектрик). Такой переход (переход Мотта, иногда говорят Мотта-Хаббарда) наблюдался в некоторых
окислах при повышении температуры (и уменьшении a).
Получим критерий Мотта, показывающий, когда будет металл, а когда диэлектрик. В металлическом
состоянии свободные электроны стремятся заэкранировать кулоновский потенциал на длине r0 (r0−2 ≡
4me2 n1/3
λ2 =
), (см. формулу(6.13a)). С другой стороны, неэкранированный кулоновский потенциал
~2
~2
. Если r0 < aB , то это связанное
образует связанное состояние с радиусом локализации aB , aB =
me2
состояние не сможет образоваться, и получится металл, если r0 > aB , то электроны окажутся связанными
38
Глава 6. Диэлектрическая проницаемость твердого тела.Формула Линдхарда.
на своих ионах и получится изолятор. Таким образом, критерий Мотта для перехода металла в диэлектрическое состояние есть: r0 = aB , или подставляя aB и λ:
4aB n1/3 < 1.
То есть при малых концентрациях электронов (большом расстоянии между атомами) получится металл,
а при больших концентрациях (малых межатомных расстояниях) — диэлектрик. Переход Мотта может
происходить также в примесной зоне полупроводника, надо только под a понимать расстояние между
примесями.
39
Глава 7
Явление переноса в твердых телах.
Кинетическое уравнение.
7.1
Кинетическое уравнение Больцмана.
Когда в предыдущих лекциях мы говорили о функции распределении электронов по импульсам и координатам, мы пользовались функцией распределения Ферми fF (или Больцмана в пределе невырожденного газа), которая описывает распределение электронов в состоянии термодинамического равновесия
и не зависит от характера взаимодействия и характера установления равновесия в системе. Сейчас мы
рассмотрим поведение электронов в неравновесном случае, когда система выведена из равновесия внешним
воздействием и в системе происходят процессы диссипации и, следовательно, функция распределения
отличается от равновесной. Оказывается, что в этом случае существенную роль играют процессы взаимодействия электронов между собой, с другими квазичастицами и примесями.
Рассмотрим сначала электрон в идеальном кристалле. Если к кристаллу приложено постоянное и
dp
однородное электрическое поле E, то согласно квазиклассическому уравнению движения
= eE кваdt
зиимпульс электрона линейно растет со временем. Если электрон находится у дна зоны, то p = m∗ v
и скорость растет как функция времени, а значит растет и ток. При дальнейшем увеличении импульса
электрон удалится от минимума энергии и уже нельзя сказать, что скорость электрона постоянно растет.
∂E(p)
будет изменяться периодически во времени, так как E(p) — периодическая функция
Скорость v =
∂p
квазиимпульса. Соответственно, будет периодически осциллировать и ток j = evn, возникающий в однородном поле. Эта картина, конечно, противоречит опытным данным, в частности, закону Ома. На самом
деле, в реальном кристалле, электрон, взаимодействуя с различными рассеивателями ( то есть с отклонениями потенциала от периодичности, возникающими из-за наличия примесей, колебаний решетки, других
электронов и т. д.) будет испытывать силу трения F f riction = −ηv. В результате уравнение движения
примет вид:
dv
m∗
= eE − ηv,
(7.1)
dt
решение которого
µ
¶
eEτ
eEτ
− τt
(7.2)
v=
+
v
−
e
0
m∗
m∗
показывает, что через время t À τ = m∗ /η скорость электрона становится постоянной независимо от
начальной скорости v0 . То есть в кристалле установится постоянный ток и будет выполняться закон Ома
j=
e2 nτ
E = σE.
m∗
(7.3)
40
Глава 7. Явление переноса в твердых телах. Кинетическое уравнение.
Время τ , за которое "забывается"первоначальная скорость (и импульс) электрона, называется временем
рассеяния импульса.
Приведенные рассуждения слишком примитивны. Ведь на самом деле электроны имеют различные
импульсы и требуется менее грубый подход, чтобы описать процессы, связанные с протеканием тока, более
строго. (Явления, связанные с протеканием тока, называются явлениями переноса или кинетическими
явлениями, или транспортными явлениями).
7.2
Кинетическое уравнение Больцмана.
Рассмотрим явления переноса в наиболее распространенном случае плавного изменения внешнего поля,
когда оно мало меняется на длине волны электрона. В этом случае применимо квазиклассическое приближение и электроны можно рассматривать как волновые пакеты, одновременно задавая их импульс
и координату. Применимость такого подхода ограничено соотношением неопределенности M p M r ∼ ~.
Так как неопределенность координат M r должна быть меньше всех характерных длин, возникающих
в явлениях переноса, в том числе длины свободного пробега ` = v̄τ ( v̄ -средняя скорость электрона),
то есть должно выполняться соотношение ` = v̄τ À ~/p, где pv̄ ∼ ε есть характерная энергия задачи
(в невырожденном газе ε ∼ kT , в вырожденном – ε ∼ EF ). Таким образом мы приходим к условию
применимости квазиклассического подхода kT τ À ~ (для невырожденного газа) и EF τ À ~ (для вырожденного газа).
Введем функцию распределения электронов в фазовом пространстве f (p, r, t). Плотность числа частиц
в области (p, p + d p), (r, r + d r) задается как
dn(p, r, t) =
2f (p, r, t) 3 3
d pd r.
(2π~)3
(7.4)
От точки к точке фазового пространства f изменяется за счет движения электронов в координатном
пространстве (v = dr/dt 6= 0) и в импульсном пространстве (dr/dt = eE 6= 0). Полная скорость изменения
f за счет таких процессов составляет
df
∂f
∂f dr ∂f dp
=
+
+
.
dt
∂t
∂r dt
∂p dt
(7.5)
Кроме того, происходит полное изменение f из-за рассеяния, когда импульс электрона изменяется изза столкновения, например, с примесью, или из-за испускания (поглощения) фонона (процессы такого
рода не учитываются квазиклассическим приближением). Поэтому полное изменение f в соответствие с
законом сохранения числа частиц (теоремой Лиувилля) записывается как
¯
∂f
∂f
∂f ¯¯
+ vOf + F
+
= 0,
(7.6)
∂t
∂p
∂t ¯st
где мы представили
dp
= F,
dt
F = eE +
e
[v, H] ,
c
(7.7)
F —сила, действующая
на электрон, E и H — электрическое и магнитное поля.
¯
∂f ¯¯
— интеграла столкновений зависит от характера столкновений. Рассмотрим сначала столВид
∂t ¯st
кновения с неподвижными дефектами (примеси, дислокации и т.д.). Примеси можно считать неподвижными,
так как их масса много больше электронной. Поэтому рассеяние электрона на примесях является упругим
(происходящим без изменения энергии).
41
Глава 7. Явление переноса в твердых телах. Кинетическое уравнение.
Обозначим вероятность (вообще говоря, квантовомеханическую) перехода электрона из состояния p
2π
в p0 как W (p, p0 ). По правилам квантовой механики W (p, p0 ) =
|hp|V̂ |p0 i|2 δ(εp − εp0 ), где V̂ (r) —
~
потенциал примеси (или другого дефекта кристалла).
Число переходов из состояния p в p0 запишем с учетом принципа Паули, который учитывает вероятность того, что начальное состояние заполнено электроном, (∼ f ) и того, что в конечном состоянии есть
незаполненные состояния (∼ (1 − f ))
W (p, p0 )f (p, r, t)[1 − f (p0 , r, t)]
2d3 p0
.
(2π~)3
(7.8)
Скорость уменьшения числа частиц в состоянии p дается интегралом по всем состояниям p0 , в которые
происходит уход электронов из-за рассеяния:
¯
Z
∂f ¯¯
2d3 p0
= − W (p, p0 )f (p, r, t)[1 − f (p0 , r, t)]
.
(7.9)
¯
∂t p→p0
(2π~)3
Число частиц в состоянии p будет также изменяться из-за переходов в это состояние электронов из других
состояний p0 . Рассуждая аналогично случаю ухода, получим
¯
Z
∂f ¯¯
2d3 p0
.
(7.10)
=
−
W (p0 , p)f (p0 , r, t)[1 − f (p, r, t)]
¯
∂t p0 →p
(2π~)3
Как известно, вероятности переходов p → p0 и p0 → p в квантовой механике одинаковы (принцип
детального равновесия. Учитывая, что W (p, p0 ) = W (p0 , p) из (7.7) и (7.10) получим интеграл столкновений
со статическими дефектами:
¯
¯
¯
Z
∂f ¯¯
∂f ¯¯
∂f ¯¯
2d3 p0
Ist {f } ≡
=
−
W (p, p0 )[f (p0 ) − f (p)].
(7.11)
=
∂t ¯st
∂t ¯p0 →p
∂t ¯p→p0
(2π~)3
Отметим, что результат (7.11) оказался таким же, каким он был бы, если бы мы не учитывали принцип
Паули (члены f (p)f (p0 ) в членах прихода и ухода взаимно уничтожили друг друга ).
7.3
Время релаксации импульса.
Выражение (7.11) можно упростить в случае изотропной энергетической поверхности (ε = ε(|p|)). Пусть
отклонение f (p) от равновесной функции распределения f0 (p) связано с присутствием слабой силы F ,
действующей на электрон. Тогда
f (p) = f0 (p) + f1 (p),
(7.12)
где f1 (p) должна быть скаляром и при этом пропорциональна силе F . Единственно возможный вид
f1 (p) = pF g(|p|) ≡ pG(ε), где g(|p|) – некоторая функция, зависящая только от модуля p (то есть от
энергии, но не от направления импульса). Подставляя f (p) в виде (7.12) в формулу (7.11) и учитывая,
что ε(p) = ε(p0 ) (упругое рассеяние), получим.
Z
2d3 p0
Ist {f } = −
W (p, p0 )G(ε)(p − p0 ) =
(2π~)3
Z
2d3 p0
f − f0
1
= −pG(ε)
W (p, p0 )(1 − cos θpp0 ) ≡ −f1 (p)
=−
,
(7.13a)
(2π~)3
τ (p)
τ (p)
42
Глава 7. Явление переноса в твердых телах. Кинетическое уравнение.
θpp0 — угол между p0 и p0 , где введено обозначение
1
=
τ (p)
Z
2d3 p0
W (p, p0 )(1 − cos θpp0 ).
(2π~)3
(7.13b)
τ (p) называется временем релаксации импульса. Форма интеграла столкновений (7.13) выведена нами для
изотропного энергетического спектра и упругого рассеяния. Но ей пользуются часто и в общем случае,
рассматривая ее как приближение. Такое приближение называют τ - приближением или приближением
времени релаксации.
Проиллюстрируем физический смысл τ . Запишем кинетическое уравнение в пространственно-однородном электрическом поле (F = eE):
∂f
∂f
f − f0
+E
=−
.
∂t
∂p
τ (p)
(7.14)
В линейном приближении по E во втором слагаемом f можно заменить равновесной функцией f0 , тогда,
∂f0
запишем опять f в виде f = f0 + f1 , и учитывая
= 0, из (7.14) получим:
∂t
∂f0
∂f1 f1
+
= −eE
.
∂t
τ
∂p
(7.15)
Решение (7.15) имеет вид:
µ
¶
t
∂f0
∂f0
f1 (p, t) = f1 (p, 0) + eEvτ
e− τ − eEvτ
,
∂ε
∂ε
(7.16)
∂f0
∂f0 ∂ε
∂f0
=
=
v.
∂p
∂ε ∂p
∂ε
Из (7.16) видно, что за время порядка τ устанавливается стационарная функция распределения (последнее слагаемое), а всякие отклонения от стационарного значения затухают (релаксируют) за время
τ.
Вычисляя ток с помощью (7.16), мы получим
Z
Z
3
2d3 p
e2 nτ
2
2 ∂f0 2d p
j = −2 evf (p)
=
2e
v
τ
E
=
E = σE,
(7.17)
(2π~)3
∂ε (2π~)3
m∗
где мы использовали, что f0 зависит от ε и где
e2 nτ
где в σ =
, как и в примитивной модели (см.(7.3)), входит время релаксации импульса τ . Таким
m∗
образом, закон Ома получается благодаря наличию рассеяния. Время τ называют также временем свободного пробега. При τ → ∞ (редкие столкновения) проводимость большая, а при частых столкновениях
(τ → 0) проводимость мала.
Отметим ряд упрощений и предположений, сделанных нами при выводе кинетического уравнения. Вопервых, вероятности рассеяния W (p, p0 ) могут зависеть от поля E, но чаще всего эта зависимость может
не учитываться, поскольку обычно E мало по сравнению с атомными масштабами электрических полей.
Во-вторых, может быть несколько видов носителей заряда, например, относящихся к разным зонам. В
этом случае надо решать несколько кинетических уравнений для функций распределения электронов в
каждой из этих зон. В случае рассеяния электрона из разных зон друг на друге эти уравнения окажутся
связанными.
43
Глава 7. Явление переноса в твердых телах. Кинетическое уравнение.
7.4
Вид интеграла столкновений при рассеянии на фононах.
Выражения (7.7) и (7.10) для скорости ухода и прихода электронов в состоянии с импульсом p можно
использовать и для случая рассеяния электронов на потенциале, созданном отклонением атомов решетки
от положения равновесия, т.е. для рассеяния на фононах. Для этого надо вычислить вероятность перехода
W (p, p0 ) с учетом зависимости потенциала, созданного колебаниями решетки от времени. Из общих
соображений ясно, что величина такого потенциала при малых смещениях пропорциональна величине
смещений атомов из положения равновесия. Как известно, вектор смещения в точке r, u(r, t) обычно
раскладывают по нормальным колебаниям (называемым фононами), имеющим собственную частоту ω(k)
i
Xh
−i(kr−ωj (k)t)
u=
uk,j ei(kr−ωj (k)t) + u+
e
,
(7.18)
k,j
k,j
а затем в соответствии с правилами квантовой механики начинают рассматривать u как оператор. При
этом очень удобным оказывается представление
чисел заполнения ( или иначе - вторичного квантования),
r
2ρωj (k)
в котором используются операторы ak,j =
uk,j (где ρ – плотность кристалла), и a+
k,j , называемые,
~
соответственно операторами уничтожения и рождения фонона ( т.е. колебания с волновым вектором k
и частотой ωj (k)). В курсе квантовой механики показано, что действие оператора ak,j (a+
k,j ) уменьшает
(увеличивает) число частиц в состоянии (k, j):
q
q
ak,j |Nj (k)i = Nj (k) |Nj (k) − 1i , a+
Nj (k) + 1 |Nj (k) + 1i ,
|N
(k)i
=
j
k,j
где |Nj (k)i— волновая функция состояния с Nj (k) числом частиц типа j (в нашем случае – фононов ветви
j) с импульсом k.
Так как оператор смещения (7.17) линеен по полевым операторам a и a+ , то и оператор потенциальной
энергии, созданной смещениями решетки будет линеен по операторам a и a+ . В самом общем виде с учетом
эрмитовости оператор электрон-фононного взаимодействия запишется как
i
Xh
0∗
0
−i(kr−ωj (k,t))
He−ph =
Hk,j
ak,j ei(kr−ωj (k,t)) + Hk,j
a+
e
,
(7.19)
k,j
k,j
где коэффициенты Hk,i зависят от конкретного механизма взаимодействия электрона со смещениями
решетки, они имеют разную зависимость от k в случаях акустических, оптических фононов и вида полей,
создаваемых колебаниями решетки (пьезоэлектрическое поле, потенциал деформации и т.д.). Так как
энергия взаимодействия должна быть пропорциональна плотности электронов, оператор He−ph должен
быть билинеен по электронным операторам.
По правилам квантовой механики вероятность перехода из начального состояния |ii в конечное |f i под
действием потенциала зависящего от времени (7.19) выражается через матричные элементы гамильтониана
возмущения:
i
¯2
¯2
¯
2π X h¯¯
0∗
0
−ikr
¯ δ(εp − εp0 + ~ω) . (7.20)
hi|Hk,j
ak,j eikr |f i ¯ δ(εp − εp0 − ~ω)+ ¯ hi|Hk,j
a+
e
|f
i
W (i → f ) =
k,j
~
k,j
Волновая функция невзаимодействующих электронов и фононов записывается в виде произведения электронной Ψ и фононной ΦiN функций
|ii = Ψp (r)ΦiN .
(7.21)
При вычислении матричных элементов в (7.20) мы учтем, что операторы a и a+ имеют ненулевые
матричные элементы только между состояниями с числом фононов N , отличающимися на единицу,
q
q
f
hΦiNj |ak,j |ΦfNj +1 i = Nj (k) + 1;
hΦiNj |a+
Nj (k),
|Φ
i
=
k,j Nj −1
44
Глава 7. Явление переноса в твердых телах. Кинетическое уравнение.
а вычисленные на блоховских функциях матричные элементы, Hk,j hΨp (r)|e±ikr |Ψp0 (r)i ∼ δp−p0 ,±k+G ,
не равны нулю только при выполнении закона сохранения квазиимпульса с точностью до одного из
векторов обратной решетки. Таким образом, вероятность перехода между состояниями p и p0 состоит из
двух частей, соответствующих первому и второму членам (7.19): W (p, p0 ) = W − + W + , и описывающих
поглощение и испускание фонона:
W − (p, p0 ) =
¯2
2π X ¯¯
Hp−p0 ,j ¯ Nj (p − p0 )δ(εp − εp0 + ~ω(p − p0 ))
~
j
W + (p, p0 ) =
¯2
2π X ¯¯
Hp−p0 ,j ¯ (Nj (p − p0 ) + 1)δ(εp − εp0 − ~ω(p + p0 )).
~
(7.22)
j
В результате получим, что при рассеянии на фононах сохраняется как энергия, так и квазиимпульс.
Подставляя (7.11) в (7.7) и (7.10) мы можем получить интеграл столкновений электронов с фононами:
¯
¯
Z
¯2
∂f ¯¯
∂f ¯¯
2π
2d3 p X ¯¯
0 ,j ¯
−
=
H
Ist {f } =
p−p
¯
¯
∂t p0 →p
∂t p→p0
~
(2π~)3
j
©
Nj (p0 − p)f (p0 )[1 − f (p)]δ(εp0 − εp + ~ω(p0 − p))+
+[1 + Nj (p0 − p)]f (p0 )[1 − f (p)]δ(εp0 − εp − ~ω(p0 − p))−
−Nj (p − p0 )f (p)[1 − f (p0 )]δ(εp − εp0 + ~ω(p − p0 ))−
(7.23)
ª
−[1 + Nj (p − p0 )]f (p)[1 − f (p0 )]δ(εp − εp0 − ~ω(p + p0 )) .
Все слагаемые в { } имеют простой физический смысл. Так, первое слагаемое описывает изменение
функции распределения электронов с импульсом p за счет перехода из состояния p0 с поглощением фонона
с импульсом p0 −p, второе — переход электрона p0 → p с испусканием фонона с импульсом p0 −p , третье —
переход электронов p → p0 с поглощением фонона с импульсом p − p0 , четвертое — p → p0 с испусканием
фонона.
Как и должно быть, вероятности поглощения фонона пропорциональны функции распределения фононов N , а вероятность испускания пропорциональны 1 + N , то есть имеет место вынужденное рождение
фононов, как и положено для бозонов. Вынужденное испускание хорошо известно из теории лазеров,
принцип действия которых основан на вынужденном рождении фотонов.
Отметим, что подстановка в (7.23) в качестве N и f равновесных функций распределений Бозе и
Ферми, соответственно, обращает интеграл столкновений в ноль (это можно проверить непосредственным
вычислением). При отсутствии внешнего поля обращается в ноль и полевая (левая) часть кинетического
уравнения. Таким образом, в равновесии решением кинетического уравнения является функция распределения Ферми.
7.5
Интеграл межэлектронных столкновений.
Рассуждая аналогично случаям рассеяния электронов на примесях и фононах, можно написать и интеграл
столкновений электронов между собой:
X ©
Ie−e {f } =
W (p01 , p02 ; p, p2 )f (p01 )f (p2 )[1 − f (p)][1 − f (p2 )]−
p01 ,p02 ,p2
ª
−W (p, p2 ; p01 , p02 )f (p)f (p2 )[1 − f (p01 )][1 − f (p02 )] ,
45
(7.24)
Глава 7. Явление переноса в твердых телах. Кинетическое уравнение.
где первый член описывает приход электронов за счет столкновений электронов p01 и p02 , которые после
рассеяния друг на друге оказываются в состояниях с импульсом p и p2 , а второй - уход электронов за
счет процесса p + p2 → p01 + p02 . Анализ вероятностей переходов W показывает, что в столкновениях
сохраняется импульсы и энергии электронов, а подстановка в качестве f функции распределения Ферми
обращает (7.24) в ноль.
7.6
Время рассеяния импульса на фононах.
Вернемся к электрон-фононному рассеянию и рассмотрим его подробнее. Механизмов взаимодействия с
колебаниями решетки много. Рассмотрим некоторые из них. Колебания решетки сопровождаются деформацией решетки, а как мы знаем, изменение расстояния между атомами изменяет зонный энергетический
спектр. Изменение спектра может быть записано как He−ph = εik uik , где uik - тензор деформации, а εik тензор деформационного потенциала. В простейшем случае изотропной зоны и акустических фононов эта
энергия имеет вид He−ph = εdivu, где ε - константа деформационного потенциала, u - смещение решетки
в волне колебаний (divu — описывает сжатие кристалла). Сделав преобразование Фурье, получим, что
i
X h
−i(kr−ωj (kt)
He−ph =
ε ikuk,j ei(kr−ωj (kt) − iku+
e
.
(7.25)
k,j
k,j
Сравнивая (7.25) с (7.17) и (7.19) видим, что, действительно, гамильтониан взаимодействия пропорционален операторам смещения, а значит и полевым операторам, причем коэффициент пропорциональности
H 0 (7.19) в имеет вид:
s
~
0
HDef
= ikE
.
(7.26)
2ρωj (N )
Потенциал деформации присутствует во всех кристаллах, но в некоторых случаях имеются и другие,
более сильные взаимодействия. Например, в ряде кристаллов механическая деформация сопровождается
их электрической поляризацией и возникновением электрического поля. Такие кристаллы называются
пьезоэлектриками (примеры: CdS, CdSe - AII B V I , GaAs, InSb - AIII B V ). Удельная поляризация, возникающая при деформации, равна Pi = βijk ujk , где βijk - тензор пьезоэлектрических модулей. Вектор
индукции равен Di = εEi +4πβijk ujk . Из уравнения divD = 0, E = −Oϕ следует уравнение для потенциала
ε M ϕ = 4πOi βijk ujk ,
из которого легко находится фурье-компонента потенциала ϕ(k)() и гамильтониана взаимодействия:
Hpiezo = eϕ(k) =
4πeiki βijk ujk (k)
.
εk 2
(7.27)
Видно, что µ
опять оператор
¶ электрон - фононного взаимодействия линеен по смещению, а учитывая, что
1 ∂uj
∂uk
1
ujk (r) =
−
; ujk (k) = (kk uj − kj uk ), ясно, что Hpiezo не зависит от модуля |k|, хотя и
2 ∂xk
∂xj
2
зависит от углов в импульсном пространстве.
В ионных кристаллах атомы обладают ненулевыми зарядами, поэтому колебания решетки сопровождаются колебаниями дипольного момента, т.е. колебаниями вектора поляризации среды P (r). Такие
колебания называются поляризационными колебаниями. Т.к. они создают электрический потенциал, они
приведут к дополнительному взаимодействию электронов с фононами. Этот вид взаимодействия возникает
с продольными оптическими фононами, т.к. именно они сдвигают ионы в элементарной ячейке навстречу
друг другу. Вычисление гамильтониана взаимодействия основано на решении уравнения Пуассона
O2 ϕ = −4πρ,
ρ = −divP ,
46
Глава 7. Явление переноса в твердых телах. Кинетическое уравнение.
где ρ — поляризационный заряд, а P = d/r0 — вектор поляризации, V0 - объем элементарной ячейки,
d = Ze [u(r, 1) − u(r, 2)], Ze — заряд атомов, u(r, 1) и u(r, 2) — смещение атомов номер 1 и номер 2.
0 ∼ 1/k.
Расчет показывает, что для полярный оптических фононов Hpol
Далее для вычисления электрон-фононных столкновений надо вычислять матричные элементы и подставлять в (7.23). Оказывается, что при взаимодействии с акустическими фононами рассеяние электронов
является почти упругим. Покажем это. Согласно (7.23), при рассеянии сохраняются энергия и импульс.
Допустим, что в результате рассеяния электрон сильно изменил направление движения, например повернул в противоположную сторону. Тогда по порядку величины | M p| = |p − p0 | ≡ k ∼ p. Изменение
энергии электронов равно как раз энергии фонона
M ε = ~ωk = sk ∼ sp.
Относительное изменение энергии электрона
Mε
sp
s
∼
= ,
ε
vp
v
q
где s — скорость звука, v — характерная энергия электронов. В невырожденном электронном v ∼
kT
m∗
— тепловая скорость, она больше скорости звука s ∼ 105 см/с при T > 1 K. В вырожденном газе v =
pF
∼ 107 ÷ 108 см/с. Таким образом, из-за того, что скорость электрона больше скорости звука, рассеяние
m
на акустических фононах является почти упругим, M ε/ε ¿ 1, т.е. при столкновениях импульс меняется
значительно, а энергия незначительно. А вследствие квазиупругости можно корректно определить (с
точностью до s/v) время релаксации электронов при рассеянии на акустических фононах, как это было
сделано для упругого рассеяния на примесях:
Z
1
2d3 p0
=
W (p, p0 )(1 − cos θpp0 ).
τ (p)
(2π~)3
Мы не будем производить подстановку соответствующих гамильтонианов в матричные элементы и вычислять интегралы. Подробный расчет приведен в учебниках, мы приведем лишь результаты для зависимостей τ от температуры и энергии электронов:
1. Для рассеяния электронов на акустических фононах через:
а)потенциал деформации — τ ∼ T −1 ε−1/2 (полупроводники), τ ∼ T −1 ε3/2 (металлы)
б) пьезопотенциал — τ ∼ T −1 ε1/2 .
Разная степень в зависимости от энергии связана с разной зависимостью от k при разных взаимодействиях.
2. Расчет времени релаксации импульса на примесях дает
а) на заряженных примесях: τ ∼ ε1/2 , а при рассеянии, б) на нейтральных примесях τ от ε не зависит.
3. Для рассеяния электронов на оптических фононах время релаксации импульса, очевидно, можно
ввести при kT À ~ωopt , когда характерные энергии электронов велики по сравнению с фононными. В
результате получится
а)τ ∼ T −1 ε−1/2 для неполярных оптических фононов (потенциал деформации),
б) τ ∼ T −1 ε1/2 для полярных оптических фононов
При kT . ~ωopt рассеяние не является упругими и, строго говоря, рассмотренная схема введения
времени релаксации не годится. Но в этом случае большинство электронов не может испустить оптических
фонон, так как энергия электрона меньше энергии оптических фононов. Мала также вероятность поглощения оптических фононов, так как при kT ¿ ~ωopt мало таких фононов. Рассеяние при низких
температурах на оптических фононах может оказаться существенным для малого числа наиболее "горячих"электронов. Такое рассеяние приводит к интересным эффектам в сильных электрических полях,
когда наблюдается разогрев электронного газа.
47
Глава 7. Явление переноса в твердых телах. Кинетическое уравнение.
При наличии нескольких независимых друг от друга механизмов рассеяния полная вероятность
P столWi . В
кновений W равна сумме вероятностей столкновений за счет различных механизмов W =
приближении времени релаксации отсюда сразу следует
1 X1
=
,
τ
τi
i
(7.28)
i
где τi — времена релаксации за счет различных механизмов рассеяния. Формула (7.28) называется правилом Матиссена. Отметим, что при рассеянии на фононах мы предполагали, что функция распределения
фононов равновесная. Это справедливо, если существуют какие-то эффективные механизмы передачи
импульса, энергии от фононов в термостат. В противном случае функция распределения фононов сама
становится неравновесной и для нее надо решать свое кинетическое уравнение для фононов. Подобная
ситуация наблюдается, например, в случае увлечения фононов электронами.
48
Глава 8
Кинетические явления. Решение уравнения
Больцмана. Проводимость.
Термоэлектрические эффекты.
8.1
Решение кинетического уравнения в τ -приближении. Отклик на однородное поле E
Ограничимся рассмотрением случая изотропного закона дисперсии электронов. Все качественные результаты, полученные в этом случае, сохраняются и в более общих случаях.
Для упругого рассеяния интеграл столкновений может быть записан в τ -приближении :
µ
¶
∂f
∂f
1
∂f
f − f0
+v
− e E + [v, H]
=−
(8.1)
∂t
∂r
c
∂p
τ
(здесь e — абсолютная величина заряда электрона.)
В общем случае зависимость f от направления импульса может быть произвольной. Функцию f можно
разложить по плоской системе функций, зависящих от угла в импульсном пространстве. Наиболее удобно
разлагать по полиномам Лежандра.
X
f (l) (|p|, r, t)Pl (cos θ).
(8.2)
f (p, r, t) =
l
Оказывается, что в линейном приближении по внешним полям можно ограничиться в разложении (8.2)
двумя первыми членами. Имея в виду отклик на слабое поле, мы так и поступим. Вычислив линейный
отклик системы на внешнее поле, мы обсудим, при каких условиях можно отбрасывать квадратичные по
полю E слагаемые, описывающие, в частности, эффект разогрева электронного газа. Таким образом
f = f (0) + f (1) cos θ.
(8.3)
Рассмотрим сначала отклик на однородное внешнее поле E, зависящее от времени. Разложим поле E
и функцию f в ряд Фурье по времени (В силу линейности задачи отклик на произвольное E(t) является
суммой по частотам откликов на поле Re {E ω exp(−iωt)}). В линейном приближении по E из (8.1) получим
∂f1
∂f0
f1
− eE
=− ,
∂t
∂p
τ
(8.4)
49
Глава 8. Кинетические явления. Решение уравнения Больцмана. Проводимость. Термоэлектрические эффекты.
где f = f0 + f1 , f0 — равновесная функция распределения. f1 — добавка. Из структуры уравнения (8.4)
видно, что добавка f1 должна быть пропорциональна углу между p и E. Поэтому решение (8.4) будет
иметь вид (8.3), где f0 = f (0) , f1 = f (1) cos θ. Обычно принято f1 записывать в виде
f1 = −
∂f0
∂f0
χ(ε)p = −
|χ(ε)||p| cos θ,
∂ε
∂ε
(8.5)
где χ(ε(|p|)) — пока ещё неизвестный вектор, вдоль которого выбирается некоторая ось (ось z для
сферической системы координат) в импульсном пространстве.
Вообще говоря, направление χ определяется направлением внешней силы, т.е. комбинацией внешних
полей E, H и т.д. В случае, когда есть только электрическое поле χ k E.
+∞
R dω
Итак, разложив (8.4) по Фурье E(t) =
E ω e−iωt , получим:
−∞ 2π
−iωf1 − eE ω
f1
∂f0
=− ,
∂p
τp
(8.6)
или
f1 =
eE ω vτp ∂f0
eE ω ∂f0
=
,
1
∂p
1 − iωτp ∂ε
− iω
τp
(8.7)
∂ε
eE ω
. Видно, что результат (8.7) действительно имеет вид (8.5) с χ =
. С помощью (8.7)
∂p
ω(1 − iωτp )
легко вычислить плотность тока j
À
¿
Z
Z
τp
2d3 p
2d3 p
ne2
Eω ,
(8.8)
j = −e
vf = −e
vf1 (p) = ∗
(2π~)3
(2π~)3
m
1 − iωτp
где v =
где
(2m∗ )3/2
hg(ε)i = −
3π 2 ~3 n
Z∞
dεε3/2
0
∂f0
g(ε),
∂ε
(8.9)
n — концентрация электронов, hg(ε)i — кинетическое усреднение по энергии. При этом мы воспользовались
p2
тем, что ε =
, а при усреднении по углам в импульсном пространстве
2m
Z
Z
Z
Z
1
2
2
dϕ sin θdθv(v, E) = dϕ sin θdθv cos θE = dϕ sin θ cos2 θdθv 2 E = 2π v 2 E.
3
Таким образом, из (8.8) получаем j(ω) = σ(ω)E(ω), где проводимость σ может быть вычислена, если
известен механизм рассеяния импульса и, следовательно зависимость τp (ε). Если считать, что τp не зависит
ne2
τp
от |p| ( только от ε), то σ(ω) = m∗ 1−iωτ . Видно, что частотная дисперсия играет роль при ω & 1/τp .
p
ne2
При ω ¿ 1/τp (мы в дальнейшeм ограничимся этим случаем) σ =
hτp i определяется усредненным
m∗
∂f0
по энергиям временем релаксации. В вырожденном электронном газе, где
= δ(ε − εF )
hτp (ε)i =
∂ε
τp (εF ). В невырожденном газе усреднение (8.9) надо выполнить честно, причем результат зависит от вида
зависимости τp .
50
Глава 8. Кинетические явления. Решение уравнения Больцмана. Проводимость. Термоэлектрические эффекты.
Формулу j = σE часто записывают в виде
j = nev dr ,
v dr = µn E,
j = neµn E,
σ = neµn ,
(8.10)
где введена подвижность µn как коэффициент между E и (средней) дрейфовой
p скоростью v dr (не путать
со скоростью v = p/m отдельного электрона и тепловой скоростью vth ∼ T /m!!!). Согласно (8.10),
µn = ehτp i/m∗ и зависимость подвижности µ от температуры определяется температурной зависимостью
hτp (ε)i. Часто зависимость τ (ε) = τ0 εr , где τ0 может зависеть от температуры. Например, для рассеяния
электронов на акустических фононах r = −1/2, τ0 ∼ 1/T , что после подстановки в (8.9) даёт µ(T ) ∼ T −3/2 ,
(такой закон удовлетворительно описывает µ(T ), например, в n-Ge.
В полупроводниках, вообще говоря, надо учитывать вклад электронов и дырок
j = e(nµn + pµp )E,
σ = e(nµn + pµp ),
где n, µn и p, µp — концентрации и подвижности электронов и дырок, соответственно.
8.2
Стационарное решение кинетического уравнения при наличии электрического и магнитного полей и градиента температуры.
∂f0
∂f
χ(ε)p и, подставив f (p) в (8.1) при
= 0, будем искать
∂ε
∂t
∂f
χ(ε). Если величина поля или температура зависят от координат, то
6= 0. Мы рассмотрим случай
∂r
слабой пространственной зависимости, когда внешние поля слабо меняются на длине свободного пробега
∂f1
f1
` = vT τ ( vF τ — для вырожденного случая). В этом случае можно пренебречь v
по сравнению с .
∂r
τ
∂f
При подстановке f0 в
получим:
∂r
½
¾
∂f0 ∂ε ∂f0 ∂T
∂f0 µ − ε
∂f0
=
+
=
OT − Oµ ,
(8.11)
∂r
∂ε ∂r
∂T ∂r
∂ε
T
Запишем опять f (p) = f0 (ε) + f1 (p), f1 = −
где опять введено обозначение µ для химического потенциала.
В полевое слагаемое с E в линейном приближении достаточно подставить f0 в
с H этого недостаточно, так как результат подстановки f0 обращается в 0:
µ
¶
∂f0
∂f0 ∂ε
∂f0
[v, H]
= [v, H]
=
([v, H], v) ≡ 0.
∂p
∂ε ∂p
∂ε
∂f
, а вот в слагаемом
∂p
∂f
надо сохранить и
∂p
½
¾
µ
½
¾¶
∂f1
∂ ∂f0
∂f0
∂ ∂f0
=−
χ(ε)p =
χ(ε) − p
χ(ε)
,
∂p
∂p ∂ε
∂ε
∂p ∂ε
Поэтому для учета влияния магнитного поля в
причём при подстановке в слагаемое с силой Лоренца второе слагаемое выпадает, поскольку ([v, H], v) ≡
0.
В интеграле столкновений, естественно, останется только f1 . В результате для χn (ε) (для определённости
рассмотрим электроны)
½
¾
1
µ−ε
e
χn
OT + O(eϕ − µ) + [H, χn ] =
.
(8.12)
mn
T
c
τn
51
Глава 8. Кинетические явления. Решение уравнения Больцмана. Проводимость. Термоэлектрические эффекты.
Это линейное уравнение имеет решение


τn
χn = −
mn




(
)
µ
¶

eτn
eτn 2

[H, Pn ] + H
(HPn ) ,
µ
¶  Pn +
cmn
cmn
eτn H 2 
1+
mn c
1
(8.13)
где введена "обобщенная сила, действующая на электрон"
Pn = O(µ − eϕ) +
ε−µ
OT.
T
(8.14)
Как и раньше, решение кинетического уравнения для дырок можно написать, не решая, а заменив
индексы ”n” на ”p”, поменяв знак заряда и заменив химпотенциал µ в (8.14) на −εG − µ (так как энергия
дырок отсчитывается вниз, а начало отсчета энергии сдвинуто на εG от вершины валентной зоны).
8.3
Ток в неоднородном проводнике и градиент электрохимического потенциала.
Отметим очень важное обстоятельство. При OT = 0 и H = 0 из кинетического уравнения мы опять
получим закон Ома, но в качестве внешнего поля будет стоять не электрическое поле E = −Oϕ, как в
однородном случае (8.8), а градиент электрохимического потенциала:
´
³µ
−ϕ .
(8.15)
j = σO
e
Отсюда, во-первых, следует известный из статистической физики критерий равновесия: электрохимический
потенциал в равновесных условиях постоянен. Во-вторых, из (8.15) следует, что наличие электрического
поля в проводнике не обязательно приводит к электрическому току. С такой ситуацией E 6= 0, j = 0
мы уже сталкивались, когда изучали экранирование. Физический смысл (8.15) становится понятнее, если
воспользоваться выражением для плотности числа частиц (5.4) через химпотенциал µ. Мы для простоты
−εG − µ
n
µ
, p = Nv exp
. Тогда On =
Oµ.
рассмотрим невырожденный случай, когда n = Nc exp
kT
kT
kT
Подставляя Oµ в (8.15), получим:
j = σE + e (Dn On − Dp Op) ,
(8.16)
µ kT
µ kT
где мы воспользовались определением σ = neµn + peµp , и ввели Dn = ne , Dn = ne — коэффициенты
диффузии электронов и дырок. Формула (8.16) справедлива и для вырожденного случая, отличие лишь
в том, что коэффициент диффузии по другому связан с подвижностью (не выполняется соотношение
Эйнштейна, справедливое для невырожденного газа).
8.4
Термоэлектрические эффекты.
Вычислим ток при наличии градиента температуры (и H = 0). Для этого подставим (8.13) в (8.5) и в
выражение для тока. Получим для тока электронов
i
en
ετn
en h
µ
j n = ∗ hτn Pn i = ∗ O(µ − eϕ)hτn i − hτn iOT +
OT =
mn
mn
T
T
µ
¶
¸
·
³
hετn i
OT
µ´
+
−µ
= σ Ẽ + ηn OT,
(8.17)
= nµn −eO ϕ −
e
hτn i
T
52
Глава 8. Кинетические явления. Решение уравнения Больцмана. Проводимость. Термоэлектрические эффекты.
где время рассеяния импульса для электронов обозначено как τn и введен градиент электрохимического
µ
потенциала Ẽ = E + O . (Не путать подвижность µn с химическим потенциалом, который как правило,
e
обозначается той же буквой µ!).
Допустим, что вклад в ток дают только электроны. Из (8.17) следует, что OT вызывает плотность
тока, которая называется термотоком. В разомкнутом образце полный ток отсутствует, значит омический
ток σ Ẽ должен компенсировать термоток, т.е. возникает градиент электрохимического потенциала и
напряжение
Z1
U12 =
Ẽdx = αn (T2 − T1 ),
(8.18)
2
где αn = ηn /σn называется коэффициентом Зеебека или дифференциальной термо-э.д.с. (Здесь T2 и T1
температуры точек 1 и 2 образца). Отметим, что возникает именно э.д.с. (разность электрохимических
потенциалов, которая вызывает ток и может совершать работу), а не просто электрическое поле. Такая
э.д.с. называется термо-э.д.с. (иначе — эффект Зеебека).
В случае, если есть несколько носителей заряда в проводнике
P
P
ηi
αi σi
P
α=
≡ P
.
(8.19)
σi
σi
Формула (8.19) получается тривиальным образом суммированием вкладов выражений (8.17) для различных типов носителей заряда. Обратим внимание на одно важное обстоятельство. Так как формула тока
(8.17) для дырок получается заменой знака перед e и µ на −µ − εG (химический потенциал отсчитывается
от экстремума зоны), то коэффициенты ηn и ηp (и, следовательно, αn и αp ) имеют противоположные
знаки, в то время как проводимости σn и σp — одинаковые. Значит, знак термо-э.д.с. определяется знаком
доминирующих носителей заряда, при этом возможна взаимная компенсация различных вкладов.
Оценим величину α. Из (8.18) и (8.19) имеем:
¸
µ ¶µ ¶·
1
hτ εi
k
−µ .
(8.20)
α=
e
kT
hτ i
В вырожденном случае (металл, вырожденный полупроводник) усреднение по энергиям h i (см.(8.9))
сводится к интегрированию τ (ε) δ - образной функций, так f0 (ε) представляет собой ступеньку, размытую
на величину ∼ kT . В главном приближении по kT /µ
hτ εi = hτ iµ и α = 0. Для получения ненулевого
вклада надо учесть зависимость τ (ε). При τ (ε) = τ0 εr вычисление интегралов hτ εi/hτ i (см. (8.9)) дает
µ ¶
µ ¶ 3
k π
3
kT
(r + )
.
(8.21)
α=
e 3
2
µ
Для невырожденного полупроводника вычисление h i сводится к интегрированию с экспонентой и
даёт
µ
¶
k
5
µ
α=
r+ −
.
(8.22)
e
2 kT
Согласно (8.22), характерный масштаб α в полупроводниках значительно больше (на несколько порядков),
чем в металлах , так как (8.21) содержит малый фактор kT /µ (в металлах). Типичные величины дифференциальной термо-э.д.с. α в металлах — мкВ/К (микровольты на градус), а в полупроводниках - мВ/К
(милливольты на градус).
53
Глава 8. Кинетические явления. Решение уравнения Больцмана. Проводимость. Термоэлектрические эффекты.
Отметим также, что из (8.18) следует, что в замкнутом однородном полупроводнике суммарная термоэ.д.с. равна нулю. Зато в неоднородном случае, когда цепь состоит из 2-х материалов A и B (см. рис.)
(термопара), наводится ненулевая термо-э.д.с.
V = αA (T1 − T2 ) − αB (T1 − T2 ) = (αA − αB )(T1 − T2 ).
Рассмотрим теперь другие тепловые эффекты. С помощью найденной нами функции распределения
можно вычислить поток энергии электронов
Z
2d3 p
(n)
Wn =
(ε − eϕ)vf1 (p).
(8.23)
(2π~)3
Полный поток энергии вычисляется суммированием по всем видам квазичастиц. Подставив в (8.23) функцию распределения и проинтегрировав, как это было сделано для плотности тока, мы получим в потоке
тепла (как в токе) слагаемые, пропорциональные полю Ẽ и OT . Исключив слагаемое с Ẽ помощью
выражения для тока, после длинных, но тривиальных вычислений получим для плотности потока энергии:
³
µ´
W = ϕ−
j − Πj − κOT,
(8.24)
e
где Π называется коэффициентом Пельтье, и κ — коэффициентом электронной теплопроводности. При
этом для Π будет выполняться равенство Π = αT , называемое соотношением Томсона, а для κ получится
выражение
κ=σ
1 hε2 τ ihτ i − hετ i2
.
e2 T
hτ i2
(8.25)
В вырожденном случае из (8.25) получится
κ=σ
π 2 kT
3e2
(8.26)
κ
Согласно соотношению (8.26), называемому законом Видемана-Франца, отношение
должно быть
σT
κ
одинаковым для всех металлов независимо от механизма рассеяния. (L =
называется числом Лоренца).
σT
Отклонения от закона Видемана - Франца могут возникнуть, если имеются отклонения от приближения
времени релаксации, а также в случаях, когда электронный вклад в теплопроводность не является доминирующим. Например в полуметаллах и полупроводниках поток тепла, переносимого фононами, может
быть сравним или быть больше электронного вклада. Для невырожденного полупроводника из (8.25) при
τn ∼ εrn , τp ∼ εrp (τn и τp времена рассеяния импульса электронов и дырок)
½·µ
¶
µ
¶ ¸ h
¾
σn σp
5
5
εG i2
k2 T
+ rn σn +
+ rp σp + (5 + rn + rp ) +
.
(8.27)
κ= 2
e
2
2
kT
σn + σp
Если доминирует один тип носителей, то из (8.27) опять получим закон Видемана-Франца, но с другим
числом Лоренца, зависящим теперь от механизма рассеяния.
Рассмотрим теперь эффекты Пельтье и Томсона . Учтем, что в стационарном случае из уравнений
∂ρ
∂ρε
непрерывности для тока и для энергии следует divj = 0 и divW = 0 (так как
=0и
= 0). Возьмём
∂t
∂t
div от потока тепла (8.24). Получим:
´i
h
³µ
−ϕ ,
(8.28)
div(−κOT ) = j OΠ + O
e
54
Глава 8. Кинетические явления. Решение уравнения Больцмана. Проводимость. Термоэлектрические эффекты.
dΠ
где OΠ = dT OT . Подставим O
³µ
´
− ϕ в (8.28) из выражения для тока, получим:
e
µ
¶
Π
j2
j2
div(−κOT ) =
+ j OΠ − OT =
+ τT (j, OT ),
(8.29)
σ
T
σ
µ ¶
dΠ Π
d Π
dα
где τT =
− =T
=T
называется коэффициентом Томсона.
dT
T
dT T
dT
Из (8.29) видно, что выделение тепла, точнее нагрев в смысле изменения температуры. (т.е. div от
потока тепла, стоящая в левой части) состоит из двух членов: тепла Джоуля (1-ый член справа) и тепла
Томсона (2-ой член),знак которого, в отличие от Джоулева тепла, изменяется при изменении направления
тока или OT .
Применим теперь (8.28) к границе 2-х проводников. Пусть при x < 0 Π = Π1 , а при x > 0 Π = Π2 .
Проинтегрировав (8.29) по x от x − δx до x + δx, получим при δx → 0:
¯
¯
∂T ¯¯
∂T ¯¯
+ κ2
= (Π1 − Π2 )j,
(8.30)
−κ1
∂x ¯1
∂x ¯2
где мы учли, что ток j непрерывен на контакте.
Слева в (8.30) стоит сумма потоков тепла к границе проводников. Согласно (8.30), в контакте выделяется
(или поглощается, в зависимости от знака правой части) тепло, пропорциональное плотности тока. Этот
эффект называется эффектом Пельтье.
55
Глава 9
Гальваномагнитные явления
9.1
Эффект Холла
Рассмотрим явления переноса тока в случае, когда к образцу приложены электрическое и магнитное поля.
Как известно из квантовой механики, в магнитном поле спектр электрона расщепляется на подзоны,
которые называются уровнями Ландау. Энергия электрона на разных уровнях отличается на величину
eH
~ωc , где циклотронная частота ωc =
. Этот результат справедлив и для магнитных полей, в которых
mc
радиус циклотронной орбиты велик по сравнению с межатомным расстоянием (что выполняется для всех
доступных полей H), и для блоховских электронов в твердом теле, если под m понимать эффективную
массу. Если ~ωc ¿ kT , то в кинетических процессах участвуют электроны со многих уровней Ландау,
и можно пользоваться квазиклассическим приближением . Такие магнитные поля называются классическими. В противоположном случае очень больших H, при которых ~ωc & kT , начинают играть роль
отдельные уровни Ландау и классическое описание не годится. Такие поля называются квантующими
магнитными полями.
Будем считать, что OT = 0 и, следовательно, опустим термомагнитные эффекты.
∂f0
χ(ε)p,
Воспользуемся общим решением для функции распределения (8.13) при OT = 0, f = f0 −
∂ε
где для электронов
χn = −
eτn
mn
E+
eτn
e2 τn2
[H, E] +
H(H, E)
mn
mn c2
,
µ
¶2
eτn
1+
H
mn c
(9.1)
а для дырок
χp =
eτp
mp
E−
e2 τp2
eτp
[H, E] +
H(H, E)
mp
mp c2
.
µ
¶2
eτp
1+
H
mp c
(9.2)
Если HkE, то вторые члены в числителе выпадают, после чего числитель и знаменатель сокращаются.
Таким образом, получаем ,что магнитное поле, приложенное вдоль электрического, никак не сказывается.
Это утверждение верно для сферических энергетических зон. В случае анизотропного электронного
спектра магнитное поле в числителе и знаменателе не сократилось бы, что означало бы влияние магнитного
поля на сопротивление (продольное магнетосопротивление). Таким образом, можно сделать вывод, что
56
Глава 9. Гальваномагнитные явления
продольное магнетосопротивление отсутствует в проводниках с изотропным энергетическим спектром
электронов.
Рассмотрим теперь случай H ⊥ E. Тогда равен нулю последний член в числителях (9.1) и (9.2).
Направим теперь ось z вдоль H. Тогда поле E может иметь компоненты Ex и Ey . Вычислим с помощью
(9.1) и (9.2) ток. Получим для тока электронов
(n)
(n)
jx(n) = σxx
Ex + σxy
Ey ;
(n)
(n)
jy(n) = σyx
Ex + σyy
Ey ,
где
(n)
σxx
=
(n)
σyy
À
τn
,
2 τ2
1 + ωcn
n
¿
À
ne2
ωcn τn2
=−
2 τ2
mn 1 + ωcn
n
ne2
=
mn
(n)
(n)
σxy
= −σyx
¿
(9.3)
(9.4)
и введена циклотронная частота для электронов ωcn =
eH
. h i означает кинетическое усреднение:
mn c
(2m∗ )3/2 R∞
∂f0
dεε3/2
g(ε).
3π 2 ~3 n 0
∂ε
Для дырок имеют место соотношения, аналогичные (9.3),(9.4), которые получаются заменой в (9.3),(9.4)
"n"на "p"и изменением знаков σxy и σyx (см. (9.1),(9.2)), которые пропорциональны нечетной степени
заряда e.
Рассмотрим теперь образец, в котором течет ток вдоль оси x, а ток jy = 0. Так как под действием
силы Лоренца электроны приобретают компоненту скорости вдоль оси yk[j, H], наряду с приложенным
полем Ex возникает и поле Ey , обеспечивающее выполнение условия jy = 0. Появление напряженности в
направлении, перпендикулярном току и магнитному полю называется эффектом Холла. Эффект Холла
характеризуют постоянной Холла R, которая определяется как коэффициент в соотношении Ey = RHjx .
Вычислим постоянную Холла , воспользовавшись соотношением j = j(n) +j(p) , где j(n) и j(p) определяются
формулами (9.3) и (9.4). Получим
hg(ε)i = −
(p)
(n)
1
σxy + σxy
R=
.
(n)
(p) 2
(n)
(p)
H (σxx
+ σxx ) + (σxy + σxy )2
(9.5)
Отметим, что слагаемые в числителе (9.5) имеют разные знаки:
*
+
¿
À
ωcp τp2
pe2
ne2
ωcn τn2
(p)
(n)
σxy =
,
σ
=
−
.
xy
2 τ2
2 τ2
mp 1 + ωcp
mn 1 + ωcn
p
n
Поэтому при полной электрон - дырочной симметрии вклады электронов и дырок в эффект Холла
полностью компенсируют друг друга и R = 0.
Рассмотрим важные частные случаи:
а) Слабое магнитное поле, ωc τ ¿ 1. При этом время свободного пробега много меньше периода
обращения электрона или дырки по циклотронной орбите. Пренебрегая σxy в знаменателе (9.5) и (σc τ )2
по сравнению с 1 в σ, получим
p 2
n
hτp i − 2 hτn2 i
2
mn
1 mp
R=
µ
¶2 .
ec n
p
hτn i +
hτp i
mn
mp
(9.6)
57
Глава 9. Гальваномагнитные явления
В случае доминирования носителей одного знака
R=
1 hτp2 i
1 hτn2 i
(p
−
тип),
R
=
−
(n − тип).
cep hτp i2
cen hτn i2
(9.7)
Формулы (9.7) часто используются для экспериментального определения концентрации и знака носителей
заряда. Если τ ∼ εr , то

1
для вырожденного газа,




2r
5
hτp2 i
hε i
√
=
= 3 π Γ(2r + 2 )
b=

hτp i2
hεi2
для невырожденного газа.


 4 Γ2 (r + 5 )
2
1
3π
Например, при рассеянии на акустических фононах (r = − ), b =
≈ 1.18, а при рассеянии на
2
8
3
315π
ионизованных примесях (r = − ), b =
≈ 1.93,.
2
512
Угол между j и E называется углом Холла
tgθ =
Ey
= RσH.
Ex
Для случая носителей одного знака
|tgθ| =
1 hτn2 i ne2 hτn i
µn H hτn2 i
=
¿ 1.
nec hτn i2 mn
c hτn i
б) Классически сильное поле, ωc τ À 1. (Но классический случай, когда квантование Ландау не учитывается, то есть kT À ~ωc . В этом случае принято говорить о классически сильном магнитном поле.).
Отбросим в (9.5) слагаемые порядка 1/ωc
µ
R≈
H
1
ne2
pe2
1
1
−
mn ωc
mp ωc
¶=
1
.
(p − n)ec
(9.8)
Отметим, что (9.7) не работает при p → n, так как при её выводе мы считали, что второе слагаемое
nc
в знаменателе больше первого, т.е. |n − p| À
. Как видно из (9.5) при точном равенстве
(µn + µp )H
n = p, R = 0 (а не ∞).
9.2
Поперечное магнетосопротивление.
Поперечным магнетосопротивлением называется зависимость сопротивления, определенного как
ρ(H) =
Ex
(E, j)
=
jx
j2
от магнитного поля H ⊥ E. Определим
Mρ
ρ(H) − ρ(0)
Ex
=
= σxx (H = 0)
−1
ρ
ρ(0)
jx
(9.9)
58
Глава 9. Гальваномагнитные явления
и подставим Ex /jx в (9.9) из (9.3) и (9.4) при условии jy = 0. Получим для случая носителей одного типа:
¿
À
τ
hτ i
Mρ
1 + ωc2 τ 2
=¿
(9.10)
À2
¿
À2 − 1.
ρ
τ2
τ
2
+ ωc
1 + ωc2 τ 2
1 + ωc2 τ 2
Из (9.3) видно ,что магнетосопротивление четно по H.
а) Слабое магнитное поле, ωc τ ¿ 1.
Из (9.10)получим
Mρ
hτ ihτ 3 i − hτ 2 i2
.
= ωc2
ρ
hτ i2
(9.11)
В пренебрежении зависимостью τ (ε), M ρ = 0. В вырожденном случае ненулевой результат получается во
втором порядке по kT
π2 2
Mρ
=
r [ωc τ (ε = µ)]2
ρ
3
µ
kT
µ
¶2
.
В невырожденном полупроводнике и τ ∼ εr

1

 9π (1 − π ),
r = − (акустические фононы)
Mρ
16
4
2
= (ωc τ )2 15π
= O((ωc τ )2 ).
6615π
3

ρ

(1 −
), r = (ионизованные примеси)
8
32768
2
б) Классически сильное поле (ωc τ ) À 1.
Mρ
1
При (ωc τ ) À 1(9.10) преобразуется к виду
= hτ ih i − 1. Опять M ρ = 0, если τ не зависит от энергии,
ρ
τ
и опять в вырожденном газе магнетосопротивление мало по сравнению с невырожденным
π2 2
Mρ
=
r
ρ
3
µ
kT
µ
¶2
− вырожденный газ

1
 32
M ρ  9π − 1, r = − 2 (акустические фононы)
= 32
− невырожденный газ.
3

ρ

− 1, r = (ионизованные примеси)
3π
2
Таким образом, в пределе сильных магнитных полей магнетосопротивление выходит на насыщение и не
зависит от H.
59
Глава 10
Разогрев электронного газа.
10.1
Время рассеяния энергии.
Как уже отмечалось, при приложении тока функция распределения начинает зависеть от углов, и удобно
представить её в виде ряда по полиномам Лежандра. В линейном приближении по электрическому полю
в функции распределения присутствовали две гармоники: гармоника с ` = 0, которая описывалась
равновесной функцией распределения, и гармоника с ` = 1, которая описывала электрический ток, все
остальные гармоники с ` > 1 — малые следующих порядков по полю E. Когда поле становится сильным,
то появляются высшие гармоники, а что важнее всего — изменяется распределение электронов по энергии,
то есть происходит разогрев электронного газа — изменяется нулевая гармоника функции распределения.
Электронный газ получает от поля мощность jE в единице обьема кристалла. Электроны передают эту
мощность фононам. Скорость отдачи энергии электронным газом характеризуется временем релаксации
энергии τe . Обсудим величину τe .
Как уже отмечалось, столкновения с примесями являются упругими, а столкновения с акустическими
фононами — квазиупругими. Столкновения с оптическими фононами при низких температурах kT <
~Ωopt вообще не играют роли, так как электронам не хватает энергии, чтобы испустить оптический фонон,
а поглотить оптический фонон невозможно, так как оптических фононов при низких температурах нет.
Рассмотрим столкновение электрона с акустическим фононом. Пусть электрон испустил фонон с
импульсом порядка его собственного (то есть фиктивно рассеял импульс и сильно изменил направление
движения) |M p| ∼ |p| . Из закона сохранения импульса следует,что M p = q, где q — импульс фонона.
Значит p ∼ q. Изменение энергии электрона M ε = ~Ωq = sq ∼ sp. Относительное изменение энергии
электрона
sp
s
Mε
∼
= .
ε
vp
v
5
Как мы
v в невырожденном полупроводнике v ∼
p уже отмечали s ∼ 10 см/c , а характерная скорость
∗
vT ∼ T /m превышает скорость звука уже при T > 10 K, в металле v ∼ vF ∼ 108 см/с À s. Таким
образом, в каждом столкновении сильно рассеивается импульс и слабо — энергия.
Но есть и вторая причина, по которой энергия рассеивается медленнее, чем импульс. Для рассеяния
импульса все равно, поглотился фонон или испустился, а скорость рассеяния энергии определяется разностью двух скоростей — испускания энергии при испускании фононов и поглощения энергии при поглощении
фононов. Вероятность поглощения фонона пропорциональна функции распределения фононов
µ
¶
~Ωq
1
1
Nq =
=
cth
−
1
,
~Ω
2
2kT
exp( kTq ) − 1
60
Глава 10. Разогрев электронного газа.
а вероятность испускания фонона пропорциональна 1 + Nq . Относительная доля превышения испускания
фононов над их поглощением равна
(NQ + 1) − NQ
~Ωq
1
=
= th
.
(NQ + 1) + NQ
2NQ + 1
2kT
Поэтому отношение обратного времени рассеяния энергии к обратному времени рассеяния импульса,
примерно совпадающему с частотой столкновений с фононами, может быть оценено как:
τe−1
s
1
s ~Ωq
−1 ≈ v 2N + 1 = v th 2kT .
τp
Q
(10.1)
Мы будем рассматривать разогрев (10.1) электронного газа в полупроводниках, где эффекты разогрева
существеннее, чем в металлах. Поэтому, учитывая, что в полупроводниках электроны находятся вблизи
экстремумов энергетических зон (энергия электрона ε ∼ vp ∼ T, Ωq ∼ sq ∼ sp ¿ T ) можно разложить
~Ωq
th
в (10.1) в ряд:
2kT
³ s ´2 ms2
τe−1
s ~Ωq
∼
∼
¿ 1. (при T & 1o K)
∼
v kT
v
kT
τp−1
Таким образом, время рассеяния энергии на фононах больше времени рассеяния импульса на тех же
фононах на фактор порядка kT /ms2 À 1.
Есть еще одна причина, по которой τe À τp . При низких температурах, когда фононы "вымораживаются"(то есть их вклад в рассеяние становится малым), рассеяние импульса часто определяется рассеянием
на примесях. Это значит, что результирующее время рассеяния импульса (то есть время рассеяния и на
примесях и на фононах) становится гораздо меньше времени рассеяния импульса на фононах, которое, в
(ph)
свою очередь, меньше времени рассеяния энергии электронами на фононах: τe À τp À τp .
Рассеяние энергий электрона может происходить из-за межэлектронных столкновений, причем при
очень низких температурах соответствующее время релаксации становится меньше времени рассеяния
энергии на фононах. Однако межэлектронное рассеяние приводит лишь к перераспределению энергии
между электронами, не приводя к отдаче энергии от электронного газа в термостат.
10.2
Горячие электроны, электронная температура.
Рассмотрим эффекты второго порядка по электрическому полю. Напишем кинетическое уравнение:
∂f
∂f
∂f
(ph)
(ee)
+v
− eE
= Ist {f } + Ist {f },
∂t
∂r
∂p
(ph)
(10.2)
(ee)
где Ist — интеграл столкновений с фононами и примесями, Ist — интеграл межэлектронных столкновений. Представим функцию распределения в виде суммы симметричной части (при малой неупругости
столкновений она с большой точностью равна функции распределения по энергиям, но не совпадает с
равновесной) и антисимметричной части, которая описывает ток. f = f0 (ε, r, t) + f1 (p, r, t). Разобьем все
слагаемые в (10.2) на симметричные и антисимметричные:
∂f0
∂f1
∂f1
(ph)
(ee)
+v
− eE
= Ist {f0 } + Ist {f0 },
∂t
∂r
∂p
(10.3)
∂f1
∂f0
∂f0
(ph)
(ee)
+v
− eE
= Ist {f1 } + Ist {f1 }
∂t
∂r
∂p
(10.4)
61
Глава 10. Разогрев электронного газа.
Рассмотрим второе уравнение, которое в линейном приближении мы уже решали. В приближении
(ph)
(ee)
времени релаксации Ist = −f1 /τp , а Ist практически не влияет на рассеяние импульса, так как в
столкновениях полный импульс электронов не меняется. Таким образом, (10.4) позволяет выразить антисимметричную часть f1 через симметричную и вычислить ток по известной f0 .
·
¸
1
∂f0
∂f0
f1 =
eE
−v
.
(10.5)
−iω + 1/τp
∂p
∂r
(ee)
В уравнении (10.3) для f0 интеграл Ist {f0 } описывает рассеяние энергии на фононах и может быть
оценен как (f0 − f0F )/τe (где f0F — фермиевская функция). Интеграл межэлектронных столкновений
обращается в 0 при подстановке равновесной функции распределения f0F (ε, T ∗ ) с любой температурой
T ∗ и имеет порядок величины (f0 −f0F (ε, T ∗ ))/τee . Общий путь решения (10.3) - (10.4) состоит в том, чтобы
подставить (10.5) в (10.3) и получить уравнение для f0 , содержащее члены ∼ E 2 , и решить его. Найденное
решение f0 позволит вычислить и среднюю энергию электронов и (с помощью (10.5)) проводимость,
которая теперь будет зависеть от электрического поля. Мы ограничимся частным, но физически простым
(ph)
случаям, когда τee ¿ τe . При очень маленьких τee последний член в (10.3) становится доминирующим, а
(ee)
значит приближенное решение (10.3) при малых τee сводится к решению уравнения Ist {f0 } = 0, которое
удовлетворяется функцией f0F (ε, T ∗ ) с любой температурой T ∗ . Таким образом, в первом приближении
f0 = f0F (ε, T ∗ ).
Значение температуры электронов T ∗ находят из уравнения баланса энергии, которое получается умножением всех слагаемых в (10.3) на энергию электронов и последующим интегрированием по импульсам:
R
Здесь u =
∂nu
+ div w − jE = −P (T ),
∂t
(10.6)
2d3 p
(2π~)3 εf0 — средняя энергия одного электрона,w =
R
2d3 p
(2π~)3 εvf1 — плотность потока энергии
электронов, j — плотность электрического тока (при получении третьего слагаемого было использовано
∂n
∂
∂ε
интегрирование по частям: ε ∂p = ∂p (εn) − n ∂p ). P (T ) =
тронами фононам. Ее можно представить в виде
P (T ) = nCe
T∗ − T
,
τe (T )
R
2d3 p
ph
(2π~)3 εIst — мощность, передаваемая элек-
(10.7)
3
где τe — время релаксации энергии (точнее электронной температуры), Ce = k — теплоемкость электрон2
ного газа, приходящаяся на один электрон. Интеграл с последнего слагаемого в (10.3) обращается в нуль,
так как он описывает уход энергии из электронной системы в результате межэлектронных столкновений.
В стационарном и однородном случае (10.6) с учетом (10.7) дает простое уравнение для электронной
температуры:
σE 2 = n
3 k(T ∗ − T )
,
2
τe
σ=
ne2 τp
,
m
то есть
T∗ = T +
3 e2 E 2
τp τe
2 mk
(10.8)
По порядку величины (10.8) можно применять для средней энергии электронного газа в поле E и для
(ph)
случая, когда τee не является доминирующим над τe , то есть представление об электронной температуре
не является строгим.
62
Глава 10. Разогрев электронного газа.
Уравнение (10.8) также дает ответ на вопрос, когда можно пренебрегать разогревными эффектами,
ограничиваясь линейными по E членами в кинетическом уравнении. Это возможно, если T ∗ − T ¿ T , то
есть, согласно (10.8), при
µ
¶
e2 E 2 τp kT
e2 E 2 `2e
e2 E 2
τp τe =
=
,
(10.9)
kT À
m
kT
m
kT
√
где мы ввели длину энергетической релаксации `e = Dτe , D — коэффициент диффузии, согласно
соотношению Эйнштейна, равный µn kT /e, где µn = eτp /m.
√
Длина энергетической релаксации `e ∼ v τp τe много больше длины свободного пробега ` = vτp .
Длина ` имеет смысл расстояния, проходимого электроном за время τe (при этом электрон движется
диффузионным образом, как и положено броуновской частице, на которую действуют
случайные силы
√
со стороны рассеивателя; путь, проходимый при диффузионном движении s ∼ Dt). Условие малости
разогревных эффектов (10.9) в результате принимает вид:
eE`e ¿ kT,
то есть, если работа, которую совершает поле над электроном на длине релаксации энергии, меньше его
тепловой энергии, то разогревом электронного газа можно пренебречь.
63
Глава 11
Контактная разность потенциалов.
Неоднородные электронные системы
11.1
Работа выхода
Для того, чтобы удалить электрон через поверхность проводника наружу, требуется совершить работу
W . Эта работа называется работой выхода, она зависит не только от объёмных свойств проводника, но
и от вида поверхности (кристаллического направления, вида обработки). Рассмотрим проводник (x < 0),
граничащий с вакуумом (x > 0). Для электрического потенциала имеем уравнение Пуассона:
d2 ϕ
= −4πρ.
d2 x
(11.1)
Интегрируем (11.1) два раза, и при втором интегрировании интегрируем по частям:
dϕ
= −4π
dx
Zx
ρdx
−∞
Zx
ϕ(x) − ϕ(−∞) = −4πx
Zx
ρdx + 4π
−∞
xρdx.
(11.2)
−∞
Ввиду электронейтральности всего проводника первый интеграл в (11.2) стремится к 0 при x → ∞.
Поэтому
+∞
Z
W = ϕ(+∞) − ϕ(−∞) = 4π
xρdx.
(11.3)
−∞
Таким образом, разность потенциалов, равная работе выхода, определяется дипольным моментом зарядов,
распределенных на поверхности. То есть работа выхода электрона связана с образованием на поверхности
двойного слоя, в котором возникают противоположные заряды на поверхности.
11.2
Контакт металл-полупроводник.
Рассмотрим контакт двух проводников, например, металла и полупроводника. В равновесии разность
электрохимических потенциалов между ними должна отсутствовать. В то же время, до контакта, каждый
64
Глава 11. Контактная разность потенциалов. Неоднородные электронные системы
Рис. 11.1: Контакт металл-полупроводник
из них имеет свой уровень химического потенциала и, вообще говоря, µS 6= µM (рис. 11.1. Это значит,
что в первый момент после контакта потечет диффузионный ток и будет течь до тех пор, пока не
наберется разность зарядов, создающая электрическое поле, препятствующее дальнейшему протеканию
тока. Равновесие установится,
когда возникнет постоянный электрохимический потенциал (уровень Ферми).
³
µ´
При этом ток j = −σO ϕ −
= 0. Или, что то же самое, j = enµE + eDOn = 0, то есть дрейфовый и
e
диффузионный токи компенсируют друг друга.
При этом в приконтактной области возникнет контактная разность потенциалов:
+∞
+∞
Z
Z
³µ´
ϕk =
Oϕdx =
= µS − µM .
O
e
−∞
(11.4)
−∞
Отметим, что в зависимости от вида материалов может быть как ϕk < 0, так и ϕk > 0. Используя
потенциальность электрического поля, можно утверждать, что разность потенциалов должна быть той
же самой, если мы переместим электрон не прямо из одного проводника в другой, (т.е. от x = −∞ в
x = +∞), а сначала переместив электрон в вакуум (достаточно далеко от проводников), а затем в другой
проводник. При этом мы получим, что ϕk = W1 − W2 , то есть контактная разность потенциалов равна
разности работ выхода материалов.
Работу выхода часто определяют как разность между уровнем химического потенциала и уровнем
потенциала вакуума.
Длина, на которой спадает приконтактное электрическое поле, очевидно, определяется длиной экранирования. В металле эта длина очень мала, порядка межатомных расстояний, а в полупроводнике
— велика, так как концентрация носителей заряда мала по сравнению с металлом. Поэтому, с большой
степенью точности, в контакте металл-полупроводник вся контактная разность потенциалов падает в
полупроводнике. Координатную зависимость ϕ(x) легко найти из условия, что в равновесии,
Z
2d3 p
µ = eϕ,
n=
fF (ε + eϕ)
(11.5)
(2π~)3
и уравнение Пуассона:
d2 ϕ
4πe
=
(n − n0 ),
2
dx
κ
(11.6)
65
Глава 11. Контактная разность потенциалов. Неоднородные электронные системы
Рис. 11.2: p-n переход
где κ — диэлектрическая проницаемость, n0 — концентрация электронов в глубине проводника. Если есть
дырка, то (11.6) надо дополнить зарядом дырок ∝ e(p − p0 ). В результате, выполняя интегрирование в
(11.5), как мы это делали в разделах 5 и 8, подставляя (11.5) в (11.6), получим нелинейное дифференциальное
уравнение второго порядка для ϕ, которое нужно решать с граничными условиями ϕ(∞) = 0, ϕ(0) = ϕk .
Если контактная разность потенциалов мала: eϕk ¿ kT , то задача сводится к линейному экранированию.
eϕ(x)
),
Например, для невырожденного полупроводника n = n0 exp(
kT
µ
¶
eϕ(x)
d2 ϕ
4πe
n0 exp
−1 .
(11.7)
=−
dx2
κ
kT
В предельном случае линейного экранирования экспоненту можно разложить в ряд:
d2 ϕ
1
= 2 ϕ,
2
dx
LD
1
4πe2 n0
=
κkT
L2D
(11.8)
ϕ(x) = ϕk exp (−
x
).
LD
(11.9)
LD — длина экранирования Дебая (смотрите разделы 6 и 8).
Аналогично контакту металл - полупроводник можно рассмотреть и контакт двух полупроводников,
например, p-n переход. В этом случае,конечно, поле E будет существенным образом проникать в оба
материала, так как длины экранирования оказываются сравнивыми (рисунок 11.3).
11.3
Низкоразмерные электронные системы.
Рассмотрим ещё один интересный тип контактов — гетеропереход — контакт полупроводников с разной
εG . На практике это обычно GaAs/Alx Ga1−x As (типичные значения x ≈ 0.3, при этом M εG ≈ 0.4
эВ). [рисунок 11.3] Современная техника позволяет получать резкую границу гетероперехода размером
порядка межатомных расстояний. В области изгиба зон в GaAs возникает ситуация, когда электроны
в зоне проводимости оказываются в потенциальной яме, по форме близкой к треугольной. Характерная
66
Глава 11. Контактная разность потенциалов. Неоднородные электронные системы
Рис. 11.3: Гетеропереход
ϕk
∼ 105 В/см. Спектр электронов в яме
LD
можно найти с помощью уравнения Шредингера, справедливого в приближении эффективной массы,
применимом, если область локализации волновой функции много больше межатомных расстояний. В
полупроводниках это, как мы уже отметили, возможно за счёт малой эффективной массы m∗ и большой
диэлектрической проницмости κ.
Спектр в яме квантуется и имеет вид
Ã
!
ky2
~2 kx2
0
εi (kx , ky ) = εi +
+
,
(11.10)
2 mx my
величина электрического поля может быть оценена как E ∼
где x − y — плоскость, параллельная границе раздела. Если kT ¿Mε0i , то все электроны сосредоточены
в одной энергетической подзоне и ведут себя как двумерные . Таким образом, в результате размерного
квантования возникает двумерный электронный газ.
Ещё один способ создания 2D-газа — приложить напряжение к металлическому электроду, отделённому
от полупроводника тонким слоем диэлектрика. Такая структура называется МДП - структура. Чаще
всего она создаётся на основе p − Si, а роль диэлектрика играет тонкий слой окисла SiO2 (∼ 10−5
см.). Приложенное к металлическому электроду ( затвору ) напряжение создаёт электрическое поле,
которое экранируется в полупроводнике, в результате чего возникает изгиб зон и потенциальная яма,
величина которой определяется приложенным к затвору напряжением VG (G — от gate — затвор). Если
напряжение таково, что дно ямы опускается ниже химического потенциала, возникает инверсионный слой
(с носителями заряда противоположного знака по сравнению с объёмом полупроводника). При низких
температурах в инверсионных слоях наблюдается электронный 2D - газ.
Плотность состояний N (ε) в этом случае (для изотропной массы):
d2 k
2 · 2π2kdk
2
2m∗
=
=
kdk
=
dε ≡ N (ε)dε
(2π)2
(2π)2
π
π~2
µ
¶
k 2 ~2
~2 k
ε=
; dε = ∗ dk; k = (kx , ky ) .
2m∗
m
2
(11.11)
Видно, что N (ε) не зависит от ε.
Аналогично двумерному электронному газу можно рассмотреть одномерный электронный газ в квантовых
проволоках, когда размерное квантование происходит по двум направлениям и энергия зависит от одной
67
Глава 11. Контактная разность потенциалов. Неоднородные электронные системы
Рис. 11.4: МДП - структура
Рис. 11.5: плотность состояний в двумерном электронном газе
68
Глава 11. Контактная разность потенциалов. Неоднородные электронные системы
q
компоненты импульса. В этом случае в пределах одной подзоны N (ε) ∝ 1/ ε − ε0i . В 0-мерной квантовой
яме (квантовой точке) N (ε) ∝ δ(ε − ε0i ).
Существуют также естественные низкоразмерные электронные системы. Например, графен - монослой
графита - это чисто двухмерная система. Дело в том, что графит имеет слоистую структуру, он состоит
из довольно слабо связанных между собой проводящих слоев углеродных атомов, образующих внутри
слоев структуру, напоминающую пчелиные соты. А углеродные нанотрубки, которые можно представить
себе, как свернутый в трубочку лист графена, являются одномерной электронной системой.
Эффекты, связанные с пониженной размерностью могут проявляться также в сильно анизотропных
трехмерных кристаллах. Если кристалл состоит из слабо связанных проводящих слоев, он называется
квазидвумерным проводником. Примером подбного материала является упомянутый выше графит, а
также многие высокотемпературные сверхпроводники. Соответственно, квазиодномерными называются
проводники, состоящие из проводящих нитей, электронные переходы между которыми затруднены. Примером квазиодномерных проводников являются многие органические проводники и проводящие полимеры,
а также ряда неорганических материалов.
Исследования свойств электронов в низкоразмерных систем являются сейчас одной из самых горячих
областей физики твердого тела. Интерес к ним вызван тем, что в одномерном и двумерном случаях в
электронных системах наблюдаются качественно новые физические явления, отсутствующие в обычных
трехмерных системах электронов, кроме того, ожидается, что низкоразмерные системы окажутся перспективными для практического использования в приборах электронной техники.
Ярким примером явления, отсутствующего в трехмерных системах, но наблюдающегося в двухмерных
проводниках, является квантовый эффект Холла. Одномерные электронные системы изучены меньше,
чем двухмерные, причем это связано не только с технологическими трудностями изготовления таких
систем. Дело в том, что к одномерным системам неприменима теория ферми-жидкости Ландау, хорошо
описывающая электроны в трехмерных проводниках. Согласно теории ферми-жидкости низкоэнергетическими элементарными возбуждениями в системе взаимодействующих электронов являются свободные
электроны и дырки. А в одномерных проводниках такие одноэлектронные возбуждения оказываются
неустойчивыми при сколь угодно слабом межэлектронном взаимодействии, так вместо ферми-жидкости
образуется состояние, называемое жидкостью Латтинджера, в котором элементарными возбуждениями
являются коллективные колебания - волны плотности заряда и плотности спина. С подобной ситуацией мы
уже встречались при описании возбуждений кристаллической решетки, где элементарными возбуждениями
оказывались не движения отдельных атомов, а коллективные фононные моды.
Следует также отметить, что при понижении размерности системы увеличивается роль тепловых и
квантовых флуктуаций, а также влияние примесей и дефектов. В одномерных системах это приводит к
очень болшим отличиям от привычного поведения обычных электронных систем.
69
Глава 12
Сверхпроводимость.
12.1
Эффективное электрон-электронное взаимодействие через фононы.
Сверхпроводимость возникает из-за межэлектронного притяжения. Откуда может возникнуть притяжение,
если известно, что электроны отталкиваются по закону Кулона? Притяжение возникает из-за взаимодействия электронов с кристаллической решёткой (с фононами). Дело в том, что электроны поляризуют
решетку, и возникающая поляризация чувствуется другими электронами, в результате чего и возникает
эффективное притяжение. Отметим, что сказанное относится к обычным, не высокотемпературным,
сверхпроводникам. Механизм высокотемпературной сверхпроводимости сложнее и до конца еще не выяснен.
Строгий расчет межэлектронного взаимодействия в присутствии фононов в реальном веществе —
очень сложная задача. Но основные качественные черты, полностью отражающие физическую природу
взаимодейстия, можно изучить на примере простейшей модели, в которой ионы описываются как набор
частиц, взаимодействующих по закону Кулона, как между собой, так и с электронами. Такое приближение
неплохо описывает продольные акустические волны в металле.
Мы до сих пор изучали экранирующее действие электронов. В частности когда выводили формулу
Линдхарда для диэлектрической проницаемости. Сейчас мы учтем экранировку взаимодействия всеми
заряженными частицами в металле — электронами и ионами, для чего вычислим полную диэлектрическую
проницаемость (εt ) с учетом вкладов электронов и ионов ("t— oт total).
Вспомним определение диэлектрической проницаемости:
εt ϕt = ϕext ,
(12.1)
где ϕt — полный потенциал, возникающий в результате реакции среды с диэлектрической проницаемостью
εt на внешний потенциал ϕext . Если бы мы рассматривали как среду только электроны, то потенциал ионов
ϕion играл бы вместе с ϕext роль внешнего потенциала, и мы написали бы
εel ϕt = ϕext + ϕion ,
(12.2)
где εel — диэлектрическая проницаемость электронов в отсутствии ионов, которую мы вычисляли по
формуле Линдхарда.
С другой стороны, мы можем считать средой "голые"(по-английски "bare") ионы, а электроны отнести
к внешнему потенциалу:
εion ϕt = ϕext + ϕel .
(12.3)
Таким образом, мы имеем три уравнения, позволяющие вместе с очевидным равенством
ϕt = ϕel + ϕion + ϕext
(12.4)
70
Глава 12. Сверхпроводимость.
найти εt через εel и εion . Сложив (12.2) и (12.3) и, вычтя из них (12.1), получим
(εel + εion − εt )ϕt = ϕel + ϕion + ϕext .
(12.5)
С учетом (12.4) из (12.5) следует
εt = εel + εion − 1,
(κt = κel + κion )
(12.6)
Отметим, что нам понадобится εel в пределе частот много меньше плазменной, поэтому воспользуемся
статическим пределом формулы Линдхарда, описывающим экранирование Томаса-Ферми:
εel = 1 +
λ2
;
q2
λ2 =
6πne2
,
εF
(12.7)
где 1/λ — радиус экранирования. Вычислим теперь диэлектрическую проницаемость голых ионов. Запишем уравнение движения голых (в отсутствии электронов) ионов в пределе высоких частот, когда ионы
ведут себя как свободные
M Ẍi = ZeE.
(12.8)
Это уравнение не описывает реальное движение ионов, так как не учитывает экранировку электронами.
Мы решаем (12.8) для вычисления εion голых ионов, а не для вычисления спектра фононов. Из (12.8)
получим:
−ω 2 M X = ZeE.
(12.9)
Для вектора поляризации P получим
µ
¶
4πZe2 N
D = E + 4πP = E 1 −
.
M ω2
N e2 Z
P = eXN = −
,
M ω2
(12.10)
Т.к. D = εE, из (12.10) получим
εion = 1 −
Ω2pl
ω2
,
(12.11)
4πZe2 N
где Ω2pl =
отличается от электронной плазменной частоты ωpl тем, что стоит заряд иона Ze и
M
масса иона M . Поэтому Ωpl ¿ ωpl .
Отметим, что было бы неправильно пытаться искать спектр колебаний ионов из условия εion = 0
(Обычное условие предельных колебаний ε = 0 следует из условия divD = divεE = 0, которое выполняется
либо при E = 0, когда нет колебаний, либо при ε = 0, E 6= 0, когда колебания есть). Правильный
способ вычислять спектр состоит в рассматривании "одетых"("dressed") ионов. Ведь в адиабатическом
приближении электроны успевают подстраиваться к мгновенному распределению ионов и, следовательно,
поля, созданные смещением ионов, ослабляются в εel раз. Это также приводит к уменьшению частоты
голых ионов Ω2pl в εel раз. В простейшем виде это легко понять из следующих рассуждений: Для голых
ионов M Ẍ = −kX/, где справа — возвращающая сила, k/M = Ω2pl . Для одетых ионов сила ослаблена
экранировкой ионов: kXεel , соответственно, фононные частоты:
Ω2q =
Ω2pl
k
=
=
εel M
εel
Ω2pl
λ2
1+ 2
q
=
Ω2pl q 2
λ2 + q 2
.
(12.12)
71
Глава 12. Сверхпроводимость.
Согласно (12.12) спектр фононов имеет в пределе малых q привычный звуковой вид: Ωq = sq, где s =
Ωpl /λ — скорость звука. Итак, подставляя (12.7) и (12.11) в (12.6), получим полную диэлектрическую
проницаемость
εt = 1 +
2
λ2 Ωpl
−
.
q2
ω2
(12.13)
С помощью (12.12), выразим εt в (12.13) через фононные частоты. Получим стандартное выражение
1
ω2
q2
.
= 2
εt (ω, q)
(q + λ2 ) ω 2 − Ω2q
(12.14)
Формула (12.14) позволяет вычислить полную энергию взаимодействия электронов. Для этого надо кулоновское взаимодействие поделить на εt (ω, q):
Ã
!
Ω2q
4πe2
4πe2 ω 2
4πe2
Vcoulomb
= 2
= 2
= 2
(12.15)
vef f (ω, q) =
1+ 2
εt
q εt (ω, q)
(q + λ2 )(ω 2 − Ω2q )
q + λ2
ω − Ω2q
4πe2
есть экранированное Кулоновское взаимодействие, а второй член — добавка
q 2 + λ2
от взаимодействия с решеткой, т.е. электрон-фононное взаимодейстие, отвечающее при малых частотах
ω < Ωq притяжению. Матричный элемент взаимодействия электронов с энергиями ε и ε + ω и импульсами
p и p + q есть
Z
hε, p|vef f (r, t)|p + q, ε + ωi = d3 rdte−iωt+iqr vef f (r, t) ≡ vef f (ω, q).
Первый член в (12.15),
Таким образом, (12.15) и определяет матричный элемент взаимодействия двух электронов, разность
энергий которых есть ω, а разность импульсов — q.
Если учесть, что плотность состояний фононов при малых q и Ω растет с Ωq :
2π 2
1 Ω2q
q
dq
=
dΩq ,
(2π)2
(2π) s3
N (Ωq ) ∝ Ω2q ∝ q 2 , то ясно, что основной вклад во взаимодействие (после суммирования по всем фононам)
дадут коротковолновые (высокочастотные) фононы с q ∼ qD ∼ π/a и Ω ∼ ωD (ωD — Дебаевская фононная
частота ωD ∼ Ωpl ). А значит притяжение будет при ω < ωD . То есть притягиваются электроны с энергиями
ε в узком слое шириной ∼ ωD около поверхности Ферми. (Напомним, что в металлах εF À ωD /).
Отметим, что кулоновское отталкивание и фононное притяжение в (12.15) при малых ω почти компенсируют друг друга, но все-таки фононное притяжение побеждает. Если выйти за рамки нашей простой
модели, то окажется, что и кулоновское отталкивание, и фононное притяжение несколько отличаются от
(12.15) в зависимости от параметров вещества может доминировать либо то, либо другое. С этим связан
тот факт, что примерно половина металлов — сверхпроводники, вторая половина — нет.
Проанализируем еще фононное притяжение, забыв о кулоновском отталкивании:
vef f (ω, q) = Vq2
Ω2q
.
ω 2 − Ω2q
(12.16)
При малых ω и q (в металлах расстояния 1/λ ∼ межатомного расстояния) из (12.16) получаем отрицательную
константу, то есть в координатном представлении мы получили бы точечное взаимодействие vef f (x) ∼
−|V0 |δ(x). Именно такой модельный потенциал в виде притяжения в одной точке соответствует модели
БКШ (Бардина-Купера-Шриффера). Дальше мы будем пользоваться моделью БКШ.
72
Глава 12. Сверхпроводимость.
12.2
Куперовские пары.
Итак, электроны притягиваются и можно было бы ожидать, что могут появиться связанные состояния в
виде пары электронов. Какие пары могли бы образоваться?
Так как взаимодействие не зависит от угла вектора q, то есть оно изотропно, оно соответствует S
- состоянию. Следовательно, волновая функция относительного движения оказалась бы зависящей от
|r1 − r2 | и была бы симметрична относительно перестановки частиц. Но поскольку полная волновая
функция должна быть антисимметрична, то значит спиновая функция не должна быть симметричной
относительно перестановки частиц, то есть у электронов не должны быть одинаковые спины. То есть в
пару могут связаться только электроны с противоположными спинами.
Иначе к этому выводу можно было придти, исходя из того, что принцип Паули запрещает двум
одинаковым фермионам находиться в одной точке, чтобы образовать связанное состояние, они должны
иметь противоположные спины.
Решение в виде связанных пар было найдено Купером. Следующий шаг мог бы быть рассмотрение
Бозе-конденсации пар, которая к моменту работы Купера уже была объяснена для Бозе-частиц. Однако
куперовские пары нельзя рассматривать как независимые частицы. Действительно, притяжение имеет
место в интервале энергий ~ωD , или импульсов электронов ~ωD /v (v — скорость на поверхности Ферми).
Согласно соотношению неопределенности это соответствует размеру потенциальной ямы v/ωD ∼ (v/ωD )a,
где a — межатомное расстояние. Значит размер пары окажется не меньше размеров ямы (на самом деле,
если посчитать, то даже намного больше). Так как среднее расстояние между электронами в металле
∼ a, то получается, что у Куперовских пар размер много больше расстояния между парами, значит
представление о сверхпроводнике, как о Бозе-газе пар — неверное.
12.3
Фазовый переход и спонтанное нарушение симметрии.
Между экспериментальным открытием (1911 г.) и теоретическим объяснением сверхпроводимости прошло
46 лет. Наряду с трудностью в объяснении природы притяжения между электронами существовала еще
одна причина, по которой потребовалось так много лет для объяснения сверхпроводимости. Дело в том,
что ее нельзя было (как и всякий фазовый переход) получить по теории возмущений.
Всякий фазовый переход второго рода связан с понижением симметрии системы. Например, в ферромагнитном состоянии, существует выделенное направление магнитного момента, то есть состояние не
является изотропным. (Симметрия восстанавливается после усреднения по всем возможным направлениям
магнитного момента). В то же время выше температуры Кюри система изотропна, то есть ее симметрия
совпадает с симметрией шара.
При образовании кристалла из газа происходит нарушение трансляционной инвариантности. В газе
система трансляционно инвариантна, а в кристалле — нет, инвариантность имеется лишь по отношению
к сдвигу на период решетки, а не на любое расстояние.
Если бы мы пытались по теории возмущений получить менее симметричное состояние из симметричного,
ничего бы не вышло, так как каждый член в ряду теории возмущений обладал бы симметрией исходного
состояния. Поэтому нужно с самого начала предположить, что в системе имеется ненулевая величина,
нарушающая симметрию, которая отсутствует в симметричной фазе. Такая величина называется параметром порядка. Например, в случае ферромагнетика параметром порядка может служить магнитный
момент, а в случае кристалла — Фурье-компонента плотности частиц с волновым числом, равным 2π/период.
Еще одно важное обстоятельство — с фазовым переходом связано нарушение какого-нибудь закона
сохранения в системе. Например, при образовании кристалла теряется закон сохранения импульса (остается
более слабый закон сохранения квазиимпульса), а в случае с ферромагнетиком теряется закон сохранения
момента импульса, связанный с изотропностью пространства. Во всех случаях идет речь о сопряженных
73
Глава 12. Сверхпроводимость.
величинах, которые связаны между собой соотношением неопределенности в квантовой механике. Теряется
инвариантность относительно смещения по х — нарушается закон сохранения импульса в примере о
кристалле. А в примере о ферромагнетике речь идет о паре угловых моментов импульса.
В сверхпроводнике в роли такой пары выступают фаза волновой функции и числом частиц, которые
в квантовой механике связаны соотношением неопределенности. Возникает определенная фаза волновой
функции и теряется закон сохранения числа частиц.
12.4
Метод самосогласованного поля в теории сверхпроводимости.
Из квантовой механики известно представление чисел заполнения (или иначе — метод вторичного квантования), в котором вводятся операторы рождения и уничтожения частицы в состоянии n: â+
n и ân .
+
Причем в качестве n можно выбирать разные переменные. Например, импульс p и спин α (âpα , âpα ). А
можно и координату. Операторы рождения и уничтожения частицы в точке r со спином α (α = 12 или 12 )
X
X
ψ̂ + (r, α) =
e−ikr â+
ψ̂(r, α) =
eikr âkα
(12.17)
kα ,
k
k
0
0
Операторы ψ̂(r, α) антикоммутируют с ψ̂(r , α) (а ψ̂ + (r, α) с ψ̂ + (r , α)), а для ψ̂ и ψ̂ + выполняется соотношение
0
0
0
ψ̂ + (r, α)ψ̂(r , α) + ψ̂(r , α)ψ̂ + (r, α) = δαβ δ(r − r ),
(12.18)
которое следует из соотношения
+
â+
pα âp0 α + âp0 α âpα = δαβ δpp0 .
(12.19)
Как показано в учебниках квантовой механики, гамильтониан в представлении операторов ψ̂ и ψˆ+ записывается как выражение для энергии системы в состоянии с волновой функцией ψ(r, α), где вместо волновой
функции стоит ψ̂(r, α). Гамильтониан, таким образом, имеет вид
Ĥ = Ĥ0 + Ĥint ,
#
"¡
¢2
XZ
p − ec A
3
+
+ U0 (r) ψ̂(r, α)
Ĥ0 =
d rψ̂ (r, α)
2m
α
Z
1X
0
0
0
0
Ĥint =
d3 rd3 r ψ̂ + (r, α)ψ̂ + (r , β)vint (r − r , α, β)ψ̂(r , β)ψ̂(r, α),
2
(12.20)
(12.21)
(12.22)
αβ
где A — векторный потенциал, U0 (r) — потенциал, в котором находятся электроны, vint — потенциал
0
0
взаимодействия. В приближении БКШ vint (r − r , α, β) = −V δ(r − r )δ↑↓ . Поэтому
Z
1
BSC
Ĥint = − V
d3 rψ̂ + (r, α)ψ̂ + (r, β)ψ̂(r, β)ψ̂(r, α)
(12.23)
2
Переход от взаимодействия к приближению самосогласованного поля состоит в замене четверки операторов на пару, умноженную на функцию, описывающую самосаглосованный потенциал, так что слагаемое
типа (12.22) или (12.23) заменяется на слагаемое типа второго члена в (12.21) с самосогласованным
потенциалом U (r), который затем определяется из минимума энергии системы.
Боголюбов предложил заменить (12.23) наR эффективный гамильтониан самосогласованного поля, в
котором введено два самосогласованных поля d3 rU (r)ψ̂ + (r, α)ψ̂(r, α) — что соответствует приближению
Хартри-Фока,
Z
h
i
ef f
Ĥint = d3 r M (r)ψ̂ + (r, ↑)ψ̂ + (r, ↓)+ M∗ (r)ψ̂(r, ↓)ψ̂(r, ↑) , где
(12.24)
74
Глава 12. Сверхпроводимость.
M (r) = V hN |ψ̂(r, ↑)ψ̂(r, ↓)|N + 2i
(12.25)
(|N i — подразумевает усреднение по состоянию системы N частиц, а |N + 2i — соответственно, из
N + 2 частиц). То есть введено предположение о существовании ненулевого среднего поля вида (12.25),
которое в обычных (несверхпроводящих) системах равно нулю, так как оно не сохраняет число частиц,
первое слагаемое в (12.24)увеличивает число частиц на 2, а второе (добавленное так, чтобы (12.24)
было эрмитовым) — уменьшает на два. Несохранение числа частиц для большой системы не является
"криминалом".
Итак, рассмотрим эффективный гамильтониан:
ef f
Ĥef f = Ĥ0 + Ĥint
(12.26)
и опустим слагаемые, содержащие внешнее поле U0 (и эффективное поле Хартри-Фока, не имеющее
отношения к сверхпроводимости, а лишь перенормирующая спектр электронов в нормальном состоянии).
Для нахождения собственных состояний гамильтониана (12.26) надо привести его к диагональному
виду
X
+
εn γ̂nα
γ̂nα ,
(12.27)
Ĥef f = E0 +
n,α
+ и γ̂
где εn — энергия квазичастицы в состоянии n, E0 — энергия основного состояния, γ̂nα
nα — полевые
операторы, соответствующие квазичастицам и связанные с исходными электронными полевыми операторами унитарным преобразованием:
X
+ ∗
[γ̂n↑ un (r) − γ̂n↓
ψ̂(r, ↑) =
vn (r)],
(12.28a)
n
X
+ ∗
ψ̂(r, ↓) =
[γ̂n↓ un (r) + γ̂n↑
vn (r)],
(12.28b)
n
где коэффициенты un и vn надо выбрать так, чтобы при подстановке (12.28) в (12.26) гамильтониан
привелся к виду (12.27). Новые операторы γ̂ и γ̂ + должны удовлетворять фермиевским соотношениям,
чтобы иметь смысл операторов квазичастиц:
+
+
γ̂nα
γ̂mβ + γ̂mβ γ̂nα
= δnm δαβ
(12.29)
γ̂nα γ̂mβ + γ̂mβ γ̂nα = 0
Чтобы найти уравнение для u(x) и v(x) вычислим коммутатор [Ĥef f , ψ̂], получим, используя (12.18)
[ψ̂(r ↑), Ĥef f ] = Ĥ0 ψ̂(r ↑)+ M (r)ψ̂ + (r ↓)
(12.30)
[ψ̂(r ↓), Ĥef f ] = Ĥ0 ψ̂(r ↓)− M (r)ψ̂ + (r ↑)
Подставим теперь в (12.30) ψ̂ и ψ̂ + через γ̂ и γ̂ + с помощью (12.28), подставив в левую часть коммутаторы
для γ̂ и γ̂ + , которые легко найти с помощью (12.27)
[Ĥef f , γ̂nα ] = −εn γ̂nα ,
(12.31)
+
+
[Ĥef f , γ̂nα
] = εn γ̂nα
.
Сравнивая затем коэффициенты перед полевыми операторами слева и справа, в получившихся равенствах,
мы получим уравнения Боголюбова
½
Ĥ0 u(r)+ M (r)v(r) = εu(r),
(12.32)
−Ĥ0∗ v(r)+ M∗ (r)u(r) = εv(r),
75
Глава 12. Сверхпроводимость.
или в матричном виде
Ĥϕ = εϕ,
где
µ
¶
Ĥ0
M
Ĥ =
;
M∗ −Ĥ0
µ
¶
u(r)
ϕ=
v(r)
(p̂ − ec A)2
p2
− F (энергия отсчитывается от уровня Ферми). В отсутствии внешних полей и в
2m
2m
пространственно однородной системе M не зависит от координат. Тогда уравнение (12.32) легко решить,
пытаясь найти u(r) и v(r) в виде плоских волн.
Ĥ0 =
u(r) = uk exp(ikr),
v(r) = vk exp(ikr).
Подставляя (12.33) в (12.32) получим:
½
(ε − ξk )uk − M vk = 0,
− M∗ uk + (ε + ξk )vk = 0,
(12.33)
(12.34)
p2
p2 − p2F
~2 k 2
(p + pF )
где ξk =
− F ≡
= (p − pF )
≈ (p − pF )vF ( при p не далеко от pF ). Система (13.32)
2m
2m
2m
2m
имеет решение, если определитель равен нулю, то есть при
p
(12.35)
ε = ± | M |2 + ξ 2 .
Итак, мы нашли энергетический спектр системы и видим, что при |ε| <M имеется энергетическая щель
(запрещенная зона шириной 2 M вокруг поверхности Ферми).
uk и vk можно найти, воспользовавшись условием нормировки
|uk |2 + |vk |2 = 1
(12.36)
Получим
u2k
1
=
2
µ
¶
ξ
1+
,
ε
vk2
1
=
2
µ
¶
ξ
1−
.
ε
(12.37)
Если сравнить (12.28) с (12.16), можно найти выражение для связи полевых операторов электронов âk и
+
â+
k с полевыми операторами нормальных возбуждений γ̂k и γ̂k :
+
âk↑ = uk γ̂k↑ − vk γ̂−k↓
,
+
âk↓ = uk γ̂k↓ − vk γ̂−k↑ .
(12.38)
В пределе |(p − pF )v| ÀM для p > pF и p < pF , (uk , vk ) переходят в (1, 0) и (0, 1) соответственно. Иными
словами, при p > pF оператор γ̂k↑ переходит в âk↑ , то есть оператор γ̂ переходит в электронный при
p > pF . При p < pF (и |p − pF |v ÀM) γ̂ переходит в a+
−k , то есть уничтожение квазичастицы означает
рождение дырки. А операторы γ̂ и γ̂ + при произвольных ξ являются линейной комбинацией операторов
рождения и уничтожения электронов, значит квазичастицы являются линейной комбинацией электрона
и дырки.
До сих пор мы не обсуждали величину M. Она должна определяется через u и v с помощью (12.25). При
нулевой температуре h i означает вычисление квантовомеханического среднего (матричного элемента),
а при T > 0, согласно принципам статистической физики
P
hn|Â|nie−En /T
n
P −E /T
hÂi =
.
(12.39)
e n
n
76
Глава 12. Сверхпроводимость.
Если вычислить M в (12.25) с помощью (12.39), подставив в явном виде ψ̂ и ψ̂ + в виде (12.28), получится
условие самосогласования для M:
X
ε
un (r)vn∗ (r)th
M (r) = V
.
(12.40)
2kT
n
Для пространственно
однородного случая u и v в виде (12.33) и (12.37) из (12.40) получится ( после
P R
перехода от
к )
~ω
ZD
M= λ
0
p
| M |2 + ξ 2
p
th
,
2kT
| M |2 + ξ 2
M dξ
(12.41)
где λ = νN (εF ), (N (εF ) — плотность состояний на поверхности Ферми). Интеграл обрезан на ~ωD ,
поскольку взаимодействуют электроны в полосе шириной ωD у поверхности Ферми.
При T = 0 из (12.41) следует
M (0) ≈ ~ωD e−1/λ ¿ ~ωD ,
при
λ < 1.
Уравнение (12.41) имеет правильное (M6= 0) решение лишь при T < Tc , где критическая температура,
при которой M обращается в ноль равна:
kT = 1.14~ωD e−1/λ .
[рисунок 13.1]
График решения уравнения (12.41) для M (T ) приведен на рисунке. M (T ) часто называют параметром
порядка сверхпроводника.
12.5
Незатухающий ток в сверхпроводнике.
Отметим сначала, что комплексный параметр M обладает свойствами, благодаря которым его часто
называют волновой функцией сверхпроводящих пар. Это будет видно из следующих соображений.
Рассмотрим сверхпроводник во внешнем поле H. Для учета H надо заменить в H0 в (12.32) p̂ = −i~O
на p̂ − ec A, где A — векторный потенциал.
Сделаем теперь унитарное преобразование функций u, v и M:
u, v, M−→ ueiχ(r)/2 , ve−iχ(r)/2 , M eiχ(r) .
(12.42)
Для новых функций u, v и M мы получим уравнения (как можно убедиться непосредственной подстановкой),
~2 ³
e ´2
~2 ³
e ´2
~2
в которых члены Ĥ0 =
p̂ − A и Ĥ0∗ =
−p̂ − A заменяются на Ĥ0 =
(p̂ − p̂s )2 и
2m
c
2m
c
2m
~
e
~2
(p̂ + p̂s )2 , соответственно, где мы ввели обозначение p̂s = Oχ − A
Ĥ0∗ =
2m
2
c

~2


(p̂ − p̂s )2 u+ M v = εu;
2m
(12.43)
2

 − ~ (p̂ + p̂s )2 v+ M∗ u = εv.
2m
Oχ в (12.43) появилась из-за дифференцирования
χ
¶ χ
´ iχ µ
³ Oχ
i
i(Oχ)
u + Ou e 2 ≡ p̂s u +
u e 2.
i~O(e 2 u) = i~ i
2
2
i
77
Глава 12. Сверхпроводимость.
Комбинация (12.43) является калибровочно-инвариантной. Действительно, из квантовой механики известно, что при изменении фазы волновой функции системы одновременно нужно добавлять слагаемое,
пропорциональное градиенту фазы к векторному потенциалу. С другой стороны из (12.43) следует, что
если мы изменим калибровку A, добавив к нему градиент произвольной функции Of , то мы получим
2e
опять уравнение (12.43) для нового набора u, v, M, в котором сделано преобразование (12.42) с χ =
f.
c~
Иными словам, при калибровочных преобразованиях M изменяется как волновая функция частиц с
зарядом 2e. Поэтому M часто называют волновой функцией куперовской пары.
Величина ps называется сверхпроводящим импульсом, а vs = ps /m — сверхпроводящей скоростью.
Причина этому следующая. Рассмотрим квантовомеханический оператор плотности тока:
³
´
ie~ ³ + ´
e2
ĵ(r) =
[ Oψ̂ (r) ψ̂(r) − ψ̂ + (r) Oψ̂(r) ] −
Aψ̂ + (r)ψ̂(r).
(12.44)
2m
mc
Он получается по обычному правилу вторичного квантования — подставить операторы ψ̂ и ψ̂ + вместо
волновых функций ψ и ψ ∗ в диагональный матричный элемент тока.
Вычислить с помощью (12.44) среднее значение оператора тока, выделив предварительно фазу у операторов ψ̂ и ψ̂ + (согласно(12.28) изменение фазы u и v изменяет фазу ψ̂ и ψ̂ + ). Произведя дифференцирование
χ
χ
´
i
i ³ Oχ
ψ̂ + Oψ̂ , получим
по правилу O(e 2 ψ̂) = e 2 i
2
µ
¶
e ~
e
e
js ≡ hĵi =
Oχ − A hψ̂ + ψ̂i ≡ ps Ns ,
(12.45)
m 2
c
m
где мы ввели обозначение Ns = hψ̂ + ψ̂i, где Ns называется плотностью сверхпроводящих электронов в
соответствии с физическим смыслом (12.45)
js = eNs vs .
(12.46)
Вычисление Ns с помощью подстановки ψ̂ и ψ̂ + из (12.38) и усреднения (12.39) приводит к результату
µ
¶
M 2
Ns → n при T → 0 и Ns →
N при T → TC (при этом M→ 0). (Отметим, что из (12.46) и
kTC
e
4π
j, rotps = − H, (так как H = rotA), rotrotH =
(12.45) немедленно следует эффект Мейснера: rotH =
c
c
4πe2
4πe2 Ns
−x/λ
−2
−
Ns H; решение у границы H = H0 e
, где λ =
.
mc2
mc2
Уравнение (12.45) является еще одной иллюстрацией к тому, что M имеет свойства волновой функции.
Если выбирать M в виде | M |e2ips r (соответствующем сверхтекучему импульсу ps на один электрон), то
параметр порядка M имеет вид волновой функции частицы, движущейся с импульсом 2ps . Это еще одно
основание для наглядного (но не вполне строгого) представления о сверхпроводнике, как о конденсате
куперовских пар.
Таким образом, мы получили, что в сверхпроводящем состоянии ток может быть отличен от нуля
без электрического поля. При этом ток создается "конденсатом"электронов, а квазичастицы с энергией
(12.35) могут вообще отсутствовать, например при T = 0 их вообще нет, так как функция распределения
(
1 при ε < − M,
fF (ε) =
0 при ε >M,
то есть нет ни электронов, ни дырок, а ток есть.
Возникает вопрос, чему равно сопротивление сверхпроводника. Представим себе сверхпроводник в
виде двух жидкостей: сверхпроводящей и нормальной (электроны и дырки — квазичастицы со спектром
(12.35)), каждая из которых вносит вклад в ток:
j = js + jN ;
js = eNs vs ;
jN = σE,
(12.47)
78
Глава 12. Сверхпроводимость.
где σ — проводимость электронов и дырок. В электрическом поле следует ожидать ускорения "пар":
dvs
= eE; −iωmvs = eE;
dt
µ 2
¶
e Ns
e
vs = i
E; j = i
+ σ E.
mω
ωm
m
(12.48)
Таким образом, комплексная проводимость сверхпроводника состоит из обычной проводимости нормальne2
ных квазичастиц σ, приводящей к джоулевым потерям (при ω 6= 0, σ =
), а также из
m(1/τ − iω)
индуктивного тока сверхпроводящего конденсата. Отметим, что сверхпроводник — это не идеальный
проводник (идеальный проводник – это проводник без дейектов, в котором время свободного пробега
τ → ∞ в σ).
Рассмотренная двухжидкостная модель не является строгой, но дает в некоторых случаях правильные
выводы, например, при T → TTC и в вихревом электрическом поле E.
Еще один хорошо известный вид проводников с коллективным механизмом проводимости, отличный от
сверхпроводников – это проводники с электронным кристаллом. Примеры: вигнеровский кристал, волны
зарядовой и спиновой плотности в квазиодномерных проводниках.
79
Related documents
Download