
книги из ГПНТБ / Немошкаленко В.В. Теоретические основы рентгеновской эмиссионной спектроскопии
.pdf90 Глава 1. Элементы еонпой теории
функциями вектора к, а величина VE (к) непрерывна в зонах или частях зон, принадлежащих одному и тому же неприводимому пред ставлению группы G (к), и после пересечения эта непрерывность сохраняется.
Причины касания энергетических зон подробно исследованы Херрингом [11].
Особый интерес представляет касание двух зон, соответствую щих одному и тому же неприводимому представлению. Две зоны,
Е |
соответствующие одному и тому же |
непри |
|||
водимому представлению, |
могут касаться, |
||||
|
если а) в кристалле нет центра инверсии, но |
||||
|
есть ось симметрии второго порядка (каса |
||||
|
ние происходит в изолированных точках, |
||||
|
лежащих в плоскостях |
симметрии или в |
|||
|
плоскости, перпендикулярной оси симмет |
||||
|
рии второго порядка); б) в |
кристалле есть |
|||
|
центр инверсии |
(контакты |
энергетических |
||
|
зон существуют |
вдоль |
некоторых |
линий |
|
Рис. 11. Пересечение зон |
в к-пространстве— линий вырождения). |
||||
Рассмотрим понятие |
о совместности. В |
||||
с разной симметрией. |
ряде случаев при помощи этого понятия мож |
||||
|
но высказать обоснованные предположения об изменении энергетических зон в кристалле. Допустим, что известно неприводимое представление, по которому преобразуются волновые функции в точке к = 0. Необходимо установить, на какие представ ления оно может быть разложено при выходе из точки к = 0 вдоль какого-либо направления высокой симметрии. Вдоль направления высокой симметрии группа симметрии является подгруппой группы симметрии при к = 0. Поэтому, если использовать в качестве матриц представления матрицы большей группы симметрии, то в меньшей группе они являются приводимым представлением, которое можно разложить обычным способом на неприводимые представления, ха рактеризующие меньшую группу симметрии. Поскольку к точке к = 0 или какой-либо другой точке подходит несколько направле ний симметрии, указанную процедуру можно выполнить для каждо го из них. Обычно эти результаты записываются в виде таблиц сов местности.
Заметим, что можно определить совместность неприводимых представлений, соответствующих оси симметрии, с неприводимыми представлениями плоскостей, содержащих эту ось. Но так как это теоретико-групповые соображения, то они не могут дать ответ на вопрос о том, какому неприводимому представлению соответствует то или иное значение энергии.
Свойства симметрии энергетических зон |
91 |
Некоторые особенности симметрии энергетических зон
Остановимся кратко на вопросах, связанных с обращением в нуль величины VE (к) и формой пересечения поверхностей постоянной энергии с осями симметрии.
При отражении в плоскостях симметрии, лежащих в зоне Бриллюэна, каждому значению энергии с данным к и некоторым неприво-
димым представлением группы G (к) со |
к' |
|
||
ответствует |
энергетический уровень с |
|
||
|
|
|||
таким же значением энергии и таким |
г |
■ |
||
же неприводимым представлением. Это |
Рис. 12. Симметрия энерге |
|||
справедливо и для плоскостей симмет |
тических зон |
по отношению |
||
рии, лежащих на поверхности |
зоны |
к плоскостям отражения. |
||
Бриллюэна. |
На рис. 12 точка к |
нахо |
|
|
дится непосредственно у плоскости симметрии, на границе зоны Бриллюэна. Параллельно этой плоскости проходит плоскость через центр зоны Г, относительно которой Е (к) имеет симметрию отраже
ния. |
Точка к' является отражением точки к в этой плоскости. Точка |
||||||
|
|
|
к" получается при трансляции |
||||
|
|
|
точки к' на вектор обратной |
ре |
|||
|
|
|
шетки вдоль линии, соединяю |
||||
|
|
|
щей точки к ' |
и к, следовательно, |
|||
|
|
|
к" |
является |
отражением |
точ |
|
|
|
|
ки к |
в плоскости симметрии, |
ле |
||
|
|
|
жащей на границе зоны Брил |
||||
|
|
|
люэна. Как отмечалось выше, |
||||
|
|
|
энергетический уровень в точ |
||||
|
|
|
ке к' имеет то же самое значение |
||||
|
|
|
и соответствует тому же непри |
||||
|
|
|
водимому представлению, что и |
||||
|
|
|
энергия в точке к. |
|
|||
Рис. |
13. Обращение величины v £ |
(к) |
Величину V Е (к) можно |
ис |
|||
следовать вблизи плоскости сим |
|||||||
|
в нуль: |
|
метрии в направлении, перпен |
||||
а — в |
невырожденной зоне, б — при |
дву |
|||||
кратном вырождении, в — прикасании зон в |
дикулярном этой плоскости. |
Ес |
|||||
одной точке (зоны соответствуют разным не |
ли зона (п) |
не |
касается другой |
||||
приводимым представлениям группы О (к)), |
|||||||
г — при касании зон в одной точке (зоны |
зоны в точках, |
расположенных |
|||||
соответствуют одному н тому же представ |
|||||||
лению группы О (к)). |
|
на плоскости симметрии, то в си |
|||||
|
|
|
лу симметрии, связанной с отра |
||||
жением и непрерывностью ^ (k ), |
проекция величины V Е (к) на на |
правление, перпендикулярное плоскости симметрии, обращается в нуль (рис. 13, а). Этот результат справедлив также в том случае, если зоны (п) и (/г -f 1) соответствуют двукратному вырождению неприводи мого представления группы G (к) при всех к , расположенных вдоль нормали к рассматриваемой плоскости (рис. 13, б). Величина VE (к)
02 Глава 1. Элементы зонной теории
обращается в нуль, если зоны (л) и (п + 1) касаются только в точке, лежащей на плоскости симметрии, но соответствуют разным непри водимым представлениям группы G (к) при к, изменяющемся вдоль направления нормали (рис. 13, в). Нормальная компонента вектора V Е (к) исчезает в обеих зонах, так как Е (к) является непрерывной функцией к в каждом неприводимом представлении и, следователь но, в каждой зоне. Если зоны (п) и (л + 1) касаются в точке, лежа щей на плоскости симметрии, и соответствуют одному и тому же неприводимому представлению группы G (к) при к, изменяющемся вдоль нормали (рис. 13, г), нормальная компонента вектора VE (к), как правило, не обращается в нуль, поскольку V Е (к) не является непрерывной функцией к в каждой зоне.
Пересечение контуров постоянной энергии, лежащих в плоско стях симметрии, с осями симметрии, лежащими в этих плоскостях, подробно исследовано в работах [12— 14].. Отметим, что проведенное Еыше рассмотрение тонких особенностей изменения Е (к) применимо также к несимморфным пространственным группам. Несколько от личается исследование свойств совместности несимморфных про
странственных |
групп. |
В этом |
случае |
необходимо рассматривать |
|
, |
G (k) |
|
|
С(к,) |
на соответствующие пред- |
фактор-группу |
У (К) |
и разложение т . |
|||
|
|
при .к |
У (Ks) |
|
|
ставления группы |
G(к) |
.ks. |
|
СИММЕТРИЯ ПО ОТНОШЕНИЮ К ИЗМЕНЕНИЮ ЗНАКА ВРЕМЕНИ
Симметрия по отношению к изменению знака времени является дополнением пространственной симметрии. Гамильтониан Н веще ствен и все его собственные значения вещественны. Если выполнить комплексное сопряжение по отношению к уравнению Шредингера, то
Нщ (г) = Е (к) фк (г).
Однако (г) — блоховская функция, соответствующая волновому
Еектору —к, поэтому Е (к) = Е (—к) при всех к. Следовательно, Е (к) обладает симметрией по отношению к инверсии к, даже если группа F не содержит инверсию /.
Понятие «симметрия по отношению к изменению знака времени» названо так потому, что если выполнить комплексное сопряжение по отношению к зависящему от времени уравнению Шредингера, то получим
ch|)*
H\f = £ д { - 0 '
В новом уравнении знак у времени изменен.
Симметрия по отношению к изменению знака времени |
93 |
Равенство Е (к) = Е (— к) представляет собой вырождение, |
до |
полнительное к вызванному симметрией по отношению к враще ниям, если точечная группа F не содержит инверсию I. Однако сим метрия по отношению к изменению знака времени может привести к еще большему вырождению.
Вопросы, связанные с изменением знака времени, рассмотрены Вигнером [7]. Пусть Г — матрицы, образующие представление группы симметрии уравнения Шредингера. Возможны следующие случаи:
а) представление Г эквивалентно представлению с вещественны ми матрицами;
б) представление Г не эквивалентно комплексно-сопряженному представлению Г*;
в) |
представление Г |
эквивалентно Г *, но не эквивалентно ника |
кому |
представлению, |
построенному полностью из вещественных |
матриц.
В случае (а) нет дополнительного вырождения. В случаях (б) и (в) есть дополнительное вырождение, собственные значения энер гии, соответствующие Г и Г *, равны, собственные функции, соответ ствующие Г*, комплексно сопряжены с функциями, соответствую щими Г, и обе системы функций линейно независимы.
Существует простой критерий [15], позволяющий устанавливать,
к какому случаю относится представление |
Г: |
|
g — случай |
(а), |
|
О — случай |
(б), |
(47) |
— g — случай |
(в) |
|
(суммирование проводится по всем преобразованиям Т группы (?)• Для электронов в кристалле группой уравнения Шредингера является группа G. Каждое неприводимое представление группы G описывается волновым вектором к и неприводимым представлением Г точечной группы G (к). Критерий (47) этого представления группы
G сводится [16] к сравнительно простой форме
|
g (к) — случай |
(а), |
|
2 ^ Р({0| Q}2) — ' |
0 — случай |
(б), |
(48) |
. — g (к) — случай |
(в) |
|
|
(суммирование проводится по всем элементам [0[ |
Q) |
точечной груп |
|
пы G (к), которая преобразует |
вектор к в вектор, |
эквивалентный |
|
— к). Заметим, что хотя {0 1Q} |
может и не быть элементом группы |
G (к), величина [0| Q)2 всегда является элементом этой группы. Рассмотрим некоторые частные случаи применения критерия (48).
1.Вектор — к не входит в звезду вектора к и не эквивалентен к.
Вэтом случае нет дополнительного вырождения Е (к). Вырож
дение, соответствующее Е (к) = Е (— к), означает совпадение
94 |
Глава 1. Элементы войной теории |
собственных значений, соответствующих различным волновым векто рам. Это возможно лишь тогда, когда группа F не содержит инвер сию I. Если cpks (г) являются базисными функциями неприводимо го представления Гр группы G (к), то базисные функции, соответ ствующие неприводимому представлению Гкр пространственной
группы G, являются функциями Р ({0 1г,-)) фь (г), вращения {0 j л£} генерируют звезду вектора к. Базисными функциями представления
(Гкр)* тогда являются функции Р ({0| г,-)) ср^ (г)*, - и ни одна из этих функций не содержит к как их волновой вектор, если вектор
— к не содержится в звезде вектора к или не эквивалентен |
к. |
2. Вектор — к содержится в звезде вектора к. В этом |
случае, |
как следует из соображений, развитых в первом пункте, существует дополнительное вырождение Е (к) вследствие изменения знака вре
мени. Так как |
представление Гкр* |
соответствует теперь двум век |
торам к и — к, |
то оно эквивалентно некоторому представлению Гкс|. |
|
В случае (б) р |
случае (в) р = |
q. Дополнительное вырождение |
в точке к есть не что иное, как вырождение между уровнями Е (к), соответствующими неприводимым представлениям Гр и Г17 группы G (к). Поэтому говорят, что Г р и Г 7 «склеиваются» в результате симметрии относительно изменения знака времени.
3. Вектор— к эквивалентен к . В этом случае существует допол нительное вырождение Е (к), вызванное симметрией по отношению к изменению знака времени, и если Гк? — представление, эквивалент ное Гкр*. то Гр и Г 7 «склеиваются».
Часто в случаях 2 и 3 существует только два неприводимых представления группы G(k) типа (б) и очевидно, что именно они явля ются «склеивающимися» представлениями. Если существует больше двух неприводимых представлений группы G (к) типа (б), необхо дим дальнейший анализ.
В случае 2, так как — к находится в звезде вектора к, существует вращение { 0 !/у} из группы G такое, что /ук = — к. По правилу построения звезды только один элемент г, из тех, что позволяют построить звезду, обладает таким свойством. Обозначим его г_. Если Гр и Г17 — два «склеивающихся» неприводимых представле ния группы G (к) типа (б), а %р и V — их характеры, то
Хр ({0 И ) |
= Г ({0| л Е ,гг1.}) |
(49) |
для каждого {0 1г} из группы |
G (к). Поэтому определить, |
какие Г 7 |
типа (б) «склеиваются» с Гр, можно по таблице характеров группы G (к). Справедливость соотношения (49) можно доказать следующим образом. Если Г кр* эквивалентно Гк<7, то функции cp{L (г)* явля ются линейными комбинациями функций Р ({01г_}) ср£, (г). Соот
Симметрия по отношению к изменению знака времени |
95 |
ветствующим преобразованием подобия функции cpw (г) могут быть
упорядочены так, что |
|
Ч& (гу* = Р ( { 0 1г_}) cph (г) |
(50) |
для каждой строки s. Для любого вращения {0 1г} из группы G (к)
|
^({0|г))сРкр;(г ) |
= |
2 г р* ((о | л })^ ;(г ). |
|
|||
Тогда |
|
|
|
|
|
|
|
р ({01Г}) [Р ({0 1г_}) <pL (г)] = |
Р ({0 1г_}) Р ((01г='гг_}) cpL (г) = |
||||||
|
= Ъ Г? ({0-| л=‘/г_,))й [Р ({01г-}) ФЬ (г)], |
|
|||||
|
t |
|
|
|
|
|
|
так |
как {0|rl|rr_i} принадлежит группе G (к). Следовательно, |
||||||
(50) |
выполняется, |
если |
|
|
|
|
|
|
|
ГР* ({01r))/s = Г 7 ((01гГ’гг_)) |
|
||||
для всех t и s, для которых справедлива формула (48). |
|
||||||
В случае 3 (— к эквивалентно к) Гкр |
эквивалентно |
Гкр, если |
|||||
функции фь (г)* |
являются |
линейными |
комбинациями |
функций |
|||
Фк/ |
(г), следовательно, |
условию (49) соответствует условие |
|||||
|
|
Хр ,({0|г}) = Хр ({0| л}) |
|
||||
для каждого {0 1/-} |
из группы G (к). |
|
|
||||
Выведем уравнение (48) из уравнения (47). Так как |
|
||||||
то, |
используя соотношение |
(31), получаем |
|
||||
2 % kp( { t » 2) = |
2 |
Xp({0|A-VV/})exp { — r7k - (rt„ + |
t„)}. (51) |
||||
|
|
R,Rptn |
|
|
|
|
Суммирование в левой части равенства проводится по всем преобра
зованиям ( t j r ) |
пространственной группы, в |
правой части — по |
|
всем гиг/ таким, что {0|г/'1г2г,-} является элементом группы G (к). |
|||
Однако |
|
|
|
2 ехр {— гг/к •(rt„ + |
t j j = 2 exp {— i (г-1 + 1) r,k •У . |
||
*n |
|
*/i |
|
Это выражение |
равно N, |
если вектор г-1г(к эквивалентен — г,к, и |
|
нулю, если вектор r- V/k |
не эквивалентен — г,к. В первом случае |
||
г7‘г2г/ является |
элементом группы G (к). |
|
|
Условие, при котором выражение (51) не обращается в нуль, |
|||
приобретает теперь более ясный смысл. Пусть |
Q — вращение, при |
||
чем такое, что |
Qk эквивалентно — к. Если |
г = r/Qr/-1, а г_1г;-к |
96 |
Глава 1. Элементы зонной теории |
|
|
эквивалентно — г/к, |
то только такие вращения сохраняются в (51). |
||
Так как r~l r2rj = |
Q2, то (51) преобразуется к виду |
|
|
^ X kp({t„|r}2) = / V M (k )SxP ({О |Q)2), |
|
||
|
|
Q |
|
откуда, учитывая, |
что порядок группы G равен g = |
Ng0 = |
|
= NM (k) g0 (k), |
получаем соотношение (48). |
|
|
Для несимморфных пространственных групп критерий, |
соответ |
ствующий критерию (48) для симморфных пространственных групп, может быть записан в виде
|
go |
случай |
(а), |
|
|
М (к) |
|
||
|
|
|
|
|
2 l p ({t |Q}2) = |
О — случай |
(б), |
(52) |
|
|
go |
•случай (в), |
|
|
|
М (к) |
|
||
|
|
|
|
где суммирование проводится по представителям (t| Q) группы
в которых Qk эквивалентно вектору — k, g0 — порядок группы F, М (к) — число векторов в звезде вектора k, %р ({Q Q}2) — характер
элемента {t |Q)2 группы G (к) в неприводимом представлении Гр. В двойных симморфных пространственных группах комплексно сопряженная спинорная функция cpft (г)* не является, вообще го воря, собственной функцией гамильтониана, содержащего члены
спин-орбитального взаимодействия. Можно показать, что
cr2 'Н (г) ст2 = Н* (г),
следовательно, если
Н (г) Фк (г) = Е (к) срк (г),
то, выполнив комплексное сопряжение, получаем
Н (г) {ст2фк (г)) = Е (к) {<т2фк (г)).
Поэтому а2фк (г) — собственная функция, соответствующая тому же собственному значению энергии, что и срк (г).
В двойных группах также возможны три случая [7]:
а) Г эквивалентно представлению, построенному из веществен
ных матриц; |
|
||
б) |
Г |
не эквивалентно |
Г *; |
в) |
Г |
эквивалентно Г *, |
но не эквивалентно представлению, по |
строенному из вещественных матриц.
Как видим, эти случаи аналогичны рассмотренным выше, однако соответствующие им вырождения в случае спин-орбитальной зави симости гамильтониана носят другой характер. В случае (а) суще ствует дополнительное вырождение. В случае (б) также существу ет дополнительное вырождение, поскольку собственные значения,
Применение |
теории групп |
в зонных вычислениях |
97 |
соответствующие представлениям Г |
и Г *, совпадают. В случае (в) |
||
нет дополнительного |
вырождения. |
|
|
При помощи неприводимых представлений двойных симморфных пространственных групп критерий, аналогичный (52), можно
записать в виде |
|
|
[ |
(к) — случай |
(а), |
2 %Р ([0 |Q]2) = |
0 — случай |
(б), |
1 — So (к) — случай |
(в), |
где суммирование ведется по всем вращениям Q точечной группы F, преобразующим к в вектор, эквивалентный— k, g0 (к) — порядок
группы G (к), %р ([01Q]2)— характер элемента [0|Q]2 в неприводи мом представлении Гр группы G (к). Элемент [01Q] может не при
надлежать группе G (к), однако элемент ЮIQ]2 всегда принадлежит этой группе.
ПРИМЕНЕНИЕ ТЕОРИИ ГРУПП В ЗОННЫХ
ВЫЧИСЛЕНИЯХ
При вычислениях энергетических зон уравнение Шредингера заменяют алгебраическими уравнениями. Для этого волновую функ цию представляют в виде разложения по базисным функциям qy (г):
Ф(г) = 2 с,-ф£(г). |
(53) |
i |
|
Подставляя это разложение в уравнение (19), умножая его левую часть последовательно на все функции qy (г) и интегрируя по объему кристалла, получают систему уравнений
|
2 с, (Hij — ESij) = 0, |
(54) |
|
|
|
где |
|
|
н а = \ |
ФС(г) H(Pi (г) dv, Sij = j фГ (г) ф/ (г) dV. |
|
V |
V |
|
Разложение (53) содержит конечное число членов, и поэтому систе ма линейных однородных уравнений (54) конечна. Для того чтобы эта система имела решение, отличное от нулевого, необходимо и до статочно, чтобы ее определитель был равен нулю, т. е.
Нп £ S n , |
Н12 ESVi . . . |
Det|//tf- E S t/| H21 — ES21, |
H22 — ES22 . . . = 0. |
Детерминант обращается в нуль только при определенных значе ниях Е, которые и являются собственными значениями энергии первоначального дифференциального уравнения. Подставляя эти
4 3 - 2 0 2 3
98 Глава 1. Элементы зонной теории
значения Et в систему (54), можно определить коэффициенты сс, вхо дящие в разложение (53).
Функции срг (г) обычно выбираются из соображений сходимости разложения (53). Однако вычислительные возможности ограничива ют число членов в этом разложении. Преодолеть возникающие труд ности можно при помощи методов теории групп. В качестве функ ций ср,- (г) в (53) следует использовать функции, преобразующиеся по какой-либо строке неприводимого представления пространст венной группы. Действительно, любую функцию можно представить
в виде ф = 2 |
2 |
2 с?тф£п (г), где индексы /и, р означают, что функ- |
( |
р |
m |
ции ф?т (г) преобразуются по m-й строке неприводимого представ
ления Гр группы G, a i — набор функции |
одного типа симметрии. |
|
В этом случае секулярное уравнение (54) |
имеет вид |
|
Det |(фjn, Яф?ш) — Е (ф/„, фГт) |= 0. |
(55) |
Столбцы и строки определителя можно поменять местами так, что все члены, соответствующие данной строке данного неприводимого представления, будут сгруппированы. При этом, согласно теореме о матричных элементах гамильтониана и ортогональности функций, преобразующихся по разному типу симметрии, величины (фfn, ф?ш) и (ф/п, Яф/'ш), входящие в (55), обратятся в нуль и определитель примет вид
D ( 1, 1) |
0 |
0 |
. |
|
0 |
£>(1.2) |
0 . |
|
|
0 |
0 |
£>( i./i) |
0 |
0 . . |
0 |
0 |
0 |
D (2, 1) |
0 . . |
0 |
0 |
0 |
0 |
D (2, 2) . |
где D [р, т)ц = (q>U Я Ф?т) - |
Е (Ф/Рт , у Ц . |
Таким образом, |
||
Det |D (р, т) |= |
0 |
|
(5б) |
|
Для каждого р и т. Заметим, |
что D (р,т ) для данного |
р |
не зависит |
|
от т. Поэтому необходимо решить только одно уравнение из урав |
||||
нений (56) для каждого неприводимого представления |
р, |
и каждое |
полученное значение энергии будет /„-кратно вырожденным. Преимущество уравнений типа (56) заключается в том, что поря
док их меньше порядка системы (54), а точность определения соб ственного значения такая же. Разумеется, каждому неприводимому представлению соответствует определенная система уравнений, но решение их проще решения исходной системы (54).
Г Л А В А 2
РАСЧЕТ ЭНЕРГЕТИЧЕСКИХ ЗОН В ТВЕРДЫ Х ТЕЛАХ
ПОСТРОЕНИЕ УРАВНЕНИЙ КРИСТАЛЛИЧЕСКОГО ПОТЕНЦИАЛА
Уравнения Хартрп — Фока
Любое твердое тело представляет собой совокупность огромного числа атомных ядер и электронов. Предположим, что ядра находятся в покое, в тех положениях, которые они занимают при абсолютном нуле температуры, и попытаемся решить уравнение Шредингера для электронов. В гамильтониане
N |
V2 |
|
, |
|
(=1 |
------г |
+ J/(',-) + 2 |
7L , |
(57) |
|
/>£ |
4 |
|
где N — число электронов, первое слагаемое представляет собой оператор кинетической энергии электронов, второе — их потенци альную энергию в поле всех ядер, третье — электростатическое взаимодействие электронов. Одним из эффективных методов нахож дения решения уравнения Шредингера Я ¥ = £ ¥ с гамильтониа ном (57) является метод Хартри — Фока [17, 18]. Волновая функ ция аппроксимируется детерминантом
¥ (хъ х2, |
. . . , xN) = - у = |
<Pl (*i) |
• • ■ |
ФN ( * д |
(58) |
|
^ t ( x N ) |
. . . |
cpN ( x N ) ’ |
||||
|
|
|
состоящим из ортоиормированных одноэлектронных функций — двухкомпонентных спиноров, которые будем считать заданными в виде
Ф,- (*/) = Ф; (г/) Че (Ф)
(пространственные части ср, (г/) для обоих направлений спина сов падают). Функция ¥ (xlt ..., xn) не является, вообще говоря, соб ственной функцией операторов S 2 (S — оператор спина) и проекции спина на ось Sz. Исключение составляют случаи, когда рассматри ваются состояния с максимальным значением спина и синглетные состояния. Во втором случае электроны кристалла делятся на две группы, в каждую из которых входят электроны с одинаковым зна чением спина. Матрица плотности для каждой группы электронов
4*