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

книги из ГПНТБ / Немошкаленко В.В. Теоретические основы рентгеновской эмиссионной спектроскопии

.pdf
Скачиваний:
14
Добавлен:
25.10.2023
Размер:
13.63 Mб
Скачать

90 Глава 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*

Соседние файлы в папке книги из ГПНТБ