Батыгин&co
.pdf80 |
Лекция б |
тования по плоским волнам можно записать так:
|
J2 Е UGCt+GaCka. |
(6.22) |
ka- |
kG a |
|
Если сюда включить еще оператор (4.39), описывающий взаимодействие в системе электронного газа, представленный также в необходимой форме, то гамильтониан взаимодействующего электронного газа в периодическом поле решетки имеет вид:
= Е ^с+аска |
+ Е Е uGc++Gi(rck<T+ |
к а |
к G сг |
Ч£ {CU М-+а,аС„аСь, - п) . (6.23)
к к' a-a' G^O
Это полный гамильтониан системы, вторично проквантованный по плоским волнам. Он очень удобен в работе и мы им будем неоднократно пользоваться.
А сейчас рассмотрим, как можно решить задачу о состояниях электронного газа в периодическом поле решетки в формализме операторов вторичного квантования. Мы будем пользоваться гамильтонианом (6.22), считая, что в системе имеется один электрон и большой набор возможных электронных состояний к. Такой прием вызван тем, что нам необходимо найти зависимость одночастичной энергии от волнового вектора. Основное состояние системы задается гамильтонианом ^ е^С^аСка, возмущением служит
ка
периодический потенциал. Обозначим одночастичную волновую функцию состояния г в компактном виде числа частиц в этом состоянии:
Поскольку имеется возмущение, то истинная волновая функция Ф не будет совпадать с одночастичной функцией, а должна быть выражена в форме линейной комбинации этих функций
- «<•••>• |
(6-24) |
Запишем далее уравнение Шредингера с гамильтонианом (6.22), причем,
6.1. Энергетический спектр электрона |
81 |
чтобы показать, что рассматривается одночастичная задача, придадим индексу кзначок «штрих» иопустим спиновый индекс а
к' G'
. • • Щ . . .) =
Умножим левую и правую часть этого уравнения на сопряженную функцию
ек. ( |
... rij ...\C£,Ck> ... |
|
г к' |
* Е |
|
|
Cj,+G,Ck, | щ ) =... ... |
|
г |
к' С |
|
|
|
аг ... |
Используя свойство ортогональности одночастичных волновых функций и раскрывая соответствующие матричные элементы, находим
г |
г |
G' |
|
г |
или, преобразуя, можно записать так: |
|
|||
|
CtiiSjOji — HjOji) + 2_^ai |
/-^tUG'Oj+G',i = U, |
||
i |
|
i |
G' |
|
однако aj(5jj =a^, аг50+о>,г |
=&J+G', тогда |
|||
|
|
|
^ ^ |
= 0 . |
Здесь, как и ранее, индексы нумеруют электронные состояния. Пусть состояние j определено как (к + G), где G — вектор обратной решетки:
ak+G(ek+G |
- |
= 0. |
82 |
Лекция б |
Сумму обратных векторов можно обозначить одним вектором
G + G' = G", |
G' = G" G. |
Таким образом, последнее выражение можно переписать так:
E) + Y^ UG"-Gak+G" = 0. |
(6.25) |
Это есть уже полученная нами ранее система уравнений (6.7) (в других обозначениях) для определения коэффициентов оц в выражении волновой функции возмущенной системы (6.24). Следовательно, тот же результат может быть получен, используя метод операторов вторичного квантования.
Лекция 7
7.1. Приближение Кронига-Пенни
До сих пор мы не делали никаких предположений, касающихся значения периодического потенциала системы ионов U(r). В действительности этот потенциал не представляет собой монотонную функцию, а имеет резкие перевалы вблизи каждого узла решетки. Это значит, что у него имеются фурье-компоненты с очень малой длиной волны, это приводит к плохой сходимости рядов, составленных из фурье-образов UQ- В связи с этим приближение почти свободных электронов в чистом виде не может быть реализовано и становится пригодным благодаря введению приема, связанного с псевдопотенциалом. Тем не менее, все качественные выводы модели почти свободных электронов носят абсолютно всеобщий характер и составляют основу всех последующих приближений. Оставляя рассмотрение указанного приема до следующего параграфа, сделаем сейчас некоторые предположения относительно потенциала U(r). Грубым приближением к реальному распределению его в одномерной решетке является предположение о наличии обрезающего потенциала С/о, позволяющее записать потенциальную энергию электрона в поле решетки в виде:
•
(r-na), (7.1)
здесь S(x) — дельта-функция Дирака. Та-
Рис 4
ким образом, потенциальную энергию электрона можно графически изобразить в форме ступенчатой кривой из одинаковых элемен-
тов-ступенек (рис. 4), а — ширина потенциальной ямы, ао — параметр решетки.
84 |
Лекция 7 |
Запишем одномерное уравнение Шредингера для электрона, движущегося в периодическом поле, аналогичное (5.5)
Мы уже знаем (стр. 63), что это линейное дифференциальное уравнение второго порядка с периодическими коэффициентами. Решениями его, согласно теореме Блоха, являются функции (5.6):
<p = eik-ru(r). |
(7.3) |
Эти решения полностью будут определены, если известна зависимость величины к от коэффициентов уравнения (7.2). Подставим потенциал Крони- га-Пенни (7.1) в уравнение (7.2):
|
2т |
г - па) |
(7.4) |
|
Ч> |
П2 s - |
|||
|
|
Рассмотрим это уравнение в главном интервале ступенчатости, т. е. там где мы выбрали начало координат. Периодичность решетки обеспечивает нам справедливость произвольного выбора начала координат:
|
|
|
|
|
-Ь |
С г С О |
II |
| |
Z//t |
п |
Г» |
^ ^ |
|
ip2 |
Н |
тгЩъ = 0 , |
0 |
$J r ^ |
а. |
|
|
|
h |
|
|
|
|
Подставляем в эти уравнения функцию Блоха (7.3):
u'{{r) + 2iku'1(r) -
и
о,"Ir\ |
_l_ Oi |
hoi* |
(<гЛ I J-2 |
'£m |
I |
I \ |
р, |
2 \ / |
< |
^ ^ 2 |
\ / — |
— — 9~^ |
^2\' |
) |
— *-*• |
Здесь удобно ввести обозначение:
l-(U0-e)=a2,
2т
Н2
(7.5)
(7.6)
(7.7)
(7.8)
7.1. Приближение Кронига - Пенни |
85 |
Тогда последние уравнения можно переписать так:
"(г) + 2г ku[(r) |
- (к2 + а2) иг(г) = 0 |
(7.9) |
м-, |
|
|
«2 (г) + 2iku'2(r) |
- (к2 - (З2) и2(г) = 0. |
(7.10) |
Решения этих уравнений хорошо известны иравны следующим выражениям:
Ul{r) |
=Ае(-гк+а)г |
+Ве-(гк+а)г |
( ? _ п ) |
и2(г) |
= Сег (-к+Р> + D е-1 {к+^г. |
(7.12) |
Решения для других участков потенциальной кривой (7.1) имеют тот же вид, что и (7.11), (7.12), лишь постоянные отличаются нафазовый множитель. Постоянные Л, В, С,D следует выбрать, требуя, чтобы функция и{г) и ее производная и'(г) были непрерывны в точках, соответствующих скачку потенциала U(r), т. е.при г = 0, г = — Ь (г = а):
ui(0)=u 2 (0), wi(0)=u'2(0),
U[(-b)=u'2(a), |
{ " ' |
|
иг(-Ъ) |
=и2{а). |
|
Периодичность решетки позволяет утверждать, чтоусловия непрерывности (7.13) должны выполняться иво всех других точках разрыва потенциала U[г). Присоединяя условия (7.13) крешениям (7.11) и (7.12), находим:
|
|
A+B-C-D |
= 0, |
(ik-a)A |
+ (ik + а)В -i(k-(3)C-i(k +[i)D = 0, |
||
де(гк-а)Ъ |
| |
g e(ik+a)b _ Qg i ( — k+ff)a _ JJg - i (k+/3)a __ g |
|
(ik- a)A e(i k~a)b + (ik + а)В e(i |
k+a)b- |
||
|
-i |
(k - [i)C el (~k+P)a |
-i(k + 0)D e~l (k+^a = 0. |
Запишем определитель этой системы уравнений относительно произволь-
Лекция 7
ных постоянных: |
|
|
|
|
|
(ik — a) |
|
(ik + a) |
|
-i(k-(3) |
|
a-(ik-a)b |
e-(ik+a)b |
|
i(k-f))a |
ег(к+(3)а |
|
~(ik — ct)b |
|
~(ik + a)b |
_pi(-k+0)a |
_p-i(k+0)a = 0. |
|
|
|
|
|
||
(ik-a) |
(ik + a) |
-i(k-P) |
-i(k + (3 |
||
1 |
1 |
- |
1 |
- |
1 |
Разрешим его, последовательно подсчитывая определители третьего порядка, относительно волнового вектора к. Имеем, после перехода к тригонометрическим и параболическим функциям:
cos(k(a +b)) = ch(ab) cos(/3a) 2af3 • sh(o:6) sin(/3a). (7.14)
Это выражение дает важную связь волнового вектора к с параметрами а и /3. Так как, согласно (7.7) и (7.8), имеем
2 2т ТТ |
Д |
2 |
(7.15) |
|
|
|
пг
то, задавая волновому вектору различные значения, можно найти зависимости а(к) иР(к), или е(к).
Используя зависимость (7.14), затем можно было бы определить исами функции и (г). Однако прямое решение уравнения (7.14) невозможно вследствие его трансцендентности. Согласно Кронигу-Пенни, это уравнение может быть значительно упрощено впредельном случае малых толщин потенциальных барьеров. Пусть bстремится кнулю, сдругой стороны можно потребовать, чтобы обрезающий потенциал Щстремился к бесконечности. Сучетом этих условий и условия (7.15) уравнение (7.14) упрощается:
cos(a£;) = cos((3a) + |
^ |
sm((3a) |
(7.16) |
||
~^2p~- |
|||||
|
|
|
|
||
Мы использовали здесь: |
|
|
|
|
|
sh(o;6) « ab, ch(ab) |
1 при |
|
|||
Введем обозначение |
am „ 7 |
|
|
|
|
|
D |
|
(7.17) |
||
|
Uob = Р, |
|
|||
|
п |
|
|
|
7.1. Приближение Кронига - Пенни |
87 |
чтобы величина Р оставалась постоянной, можно потребовать постоянства иф при переменных UQ И Ь.
С учетом введенных обозначений выражение |
(7.16) принимает вид: |
cos(afc) = cos(/3a) + P sm^a>, |
(7.18) |
Это хорошо известное уравнение Кронига-Пенни, определяющее явную связь между собственным значением энергии е и волновым вектором к. Величина Р теперь является приведенным обрезывающим потенциалом.
Очевидно, что если правая часть уравнения (7.18) будет меньше единицы, то волновой вектор является вещественным числом и решения (7.11) и (7.12) имеют смысл, если правая часть больше единицы, то к есть мнимая величина и (7.11) и (7.12) обращаются в бесконечность. Трансцендентное уравнение (7.18) уже можно решить графически для любого значения вектора к. Для этого построим зависимость правой части уравнения от (За; Если Ра = О, то, учитывая, что
lim - ^ ± = 1,
„, ^n (_X
находим
С ростом (Задо +тг эта функция убывает, становясь равной —1 при тг =(За, при /За > тг функция продолжает убывать и, достигая минимума, затем вновь растет, принимая при /За = 2тг значение +1. Далее при /За > 2тг она становится больше единицы, достигает максимума и затем опять убывает и т. д. Аналогичная картина складывается, когда /За изменяется в сторону отрицательных углов. На рис. 5 приведен качественный ход рассматриваемой зависимости
Очевидно, что при |/3a| = 2ттп, где п = ± 1, ±2, ... , ±ЛГ, функция / = +1, при |/За| = (2п + 1)тг имеем / = - 1 .
Значения /За, удовлетворяющие условию (7.18), т. е. корни этого уравнения, получим, проводя прямые параллельные оси /За, и на расстоянии cos ka от нее. Меняя ка от 0 до тг и проводя соответствующие прямые,
Рис. 5
параллельные оси (За, находим точки пересечения прямых с кривой, описываемой функцией /. Эти точки и есть решения трансцендентного уравнения (7.18). При этом видно, что там, где |/| > 1, вещественных корней нет. Это значит, что значения (За,соответствующие |/| > 1, а значит и энергии, не являются в этих интервалах собственными энергиями в уравнении Шредингера (7.2). Следовательно, весь интервал изменения (За, а значит и е, является дискретным, т. е. распадается на зоны разрешенных и запрещенных энергий и значений /За. Разрешенные значения /За на рис. 5 показаны жирной чертой. С возрастанием \(3а\ ширина разрешенных значений \(3а\, а значит и энергий, растет за счет уменьшения запрещенных. Каждому разрешенному значению энергии соответствует два значения(За, отличающиеся знаком. Следовательно, приближение Кронига-Пенни дает нам тот же результат, что и приближение почти свободных электронов, т. е., спектр энергий электрона в периодическом поле решетки состоит из непрерывных полос, разделенных интервалами запрещенных значений энергии. Равноценность выводов обоих приближений позволяет утверждать, что распад энергетического спектра на полосы в приближении Кронига-Пенни не связан с принятыми предельными условиями. Обсуждая уравнение Кро- нига-Пенни, мы не делали никаких заключений о величине приведенного обрезывающего потенциала. Однако энергетический спектр электрона су-
щественно зависит |
от этой величины. Предположим, что Р |
= 0, тогда, |
|
согласно (7.18) |
|
2 |
|
|
cos(fca) = cos((3a), |
s = i^k2 |
(7.19) |
и, следовательно, ±fca + 2тгп = (За. Это значит, что (За может принимать
7.1. Приближение Кронига - Пенни |
89 |
любые значения, т. е. разрешенными являются все значения энергии то нуля до бесконечности. Такая ситуация, как мы знаем, свойственна свободным электронам, когда энергетический спектр непрерывен и отсутствуют интервалы запрещенных энергий.
Предположим другой крайний случай, т. е. Р = со. Мы знаем, что если величина Р растет, то, согласно рис. 5, зоны дозволенных значений (За уменьшаются и когда Р = со эти зоны вырождаются в дискретные уровни. Действительно, если внимательно изучить рис. 5, то легко увидеть, что в этом предельном случае величина (За вообще не зависит от к, а определяется соотношением
(За= дат, где |
п = ±1, ±2, . . . . |
(7.20) |
или, раскрывая значение (3, находим |
|
|
^ 4 |
п = ±1,±2, .... |
(7.21) |
Как известно, эта формула определяет энергетические уровни электрона в изолированном атоме. Следовательно, этот случай соответствует полностью связанному электрону. Сопоставляя оба предельных случая, можно сказать, что величина обрезывающего приведенного потенциала Р характеризует энергию связи электрона, т. е. его свободу или локализацию.
Рассмотрим еще случай, когда величина Р сравнительно велика и электроны сильно связаны. Зоны дозволенных энергий тесно примыкают к значениям /За, равным дат. Пусть ширина этих зон есть Д, тогда разрешенные значения /За можно задать так:
(За = пи + Д, |
где |
Д « 0 . |
(7.22) |
Подставляем это условие сильно связанных электронов в уравнение Кро- нига-Пенни:
cos(afc) = (-1)" + РЩ^А. |
(7.23) |
Действительно, так как
cos(/3a) = cos(n7r + Д) = cos(nn) cos Д + sin(nTr) sin Д = ( — 1)™, sin((3a) = sin(mr + Д) = sin(nTr) cos Д + cos(nTr) sin Д = (—1)" Д, pa = nn + Д ~ П7Г.