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

книги / Принципы динамической теории решетки

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
19.18 Mб
Скачать

142

 

Глава 2

 

а

6

в

г

Лошш/щиная

Локализованная

 

м ода

10(1)

.

 

О

05

1,0 О

0,5

1,0 О

015

1,0 0

0,5

1,0

Состав

Рис. 2.16. Типы смешанного поведения кристаллов, a разрешены ка к щеле­ вая, так и локализованная моды (двухмодовое поведение); б — разрешена

только локализованная мода (смешанно-модовое поведение); в — разреше­ на только щелевая мода (смешанно-модовое поведение); г - н е разрешены ни щелевая, ни локализованная моды (одномодовое поведение).

в следующем разделе. Вообще говоря, основные особенности поведе­ ния мод в настоящий момент поняты довольно хорошо, по крайней ме­ ре количественно. Однако существуют некие тонкие детали в спектрах смешанных кристаллов, которые все еще представляют интересные объекты исследований.

Вернемся теперь к вкладу в однофононные инфракрасный и ком­

бинационный спектры смешанных кристаллов от плотности состояний, определяемой примесями. Он возникает вследствие устранения пра­ вил отбора по импульсу = 0) кристалла совершенно аналогично слу­ чаю одиночного дефекта (см. разд. 2.2.2). Совсем не просто идентифи­ цировать экспериментально такой вклад, потому что он может быть замаскирован вкладом зон с к « 0, зон резонансных мод, а также пе­

рекрытием однофононного и многофононного спектров. Для иллюстра­ ции на рис. 2.17 изображена низкочастотная область спектра комбина­

ционного рассеяния в смешанном кристалле Gel , где зона ниже

120 см "1 приписывается рассеянию ТА-фононов Ge благодаря наличию беспорядка, а ниже 170 см "1 - рассеянию ТА-фононов Si. Отметим, что остальные пики на этом рисунке идентифицируются как двухфонон-

Динамика решетки кристаллов с беспорядком замещения

143

Рис. 2.17. Спектры первого и второго порядков комбинационно! о рассеяния

в кристаллическом плаве Ge1 _ <,Sic при Т - 330 К для С О^340 с м '1. (Соглас­

но [231] .)

ные. Вклад плотности состояний в спектры будет также обсуждаться в последующих разделах.

2 .3 .2 . Приближение виртуального кристалла и

модель равных смещений

Как обсуждалось в предыдущем разделе, смешанные кристаллы можно классифицировать по отношению к зависимости от концентрации длин­ новолновых оптических фононов. Сейчас мы рассмотрим два простых

144 Глава 2

приближения для описания такой зависимости от концентрации: при­ ближение когерентного потенциала (ПКП) и модель равных смещений (PC). В ПКП смешанный кристалл заменяется на упорядоченный крис­ талл с массами и силовыми постоянными, усредненными по конфигу­

рационному пространству, в то время как в модели РСРС смешанный

кристалл заменяется на упорядоченный кристалл с псевдоэлементарной ячейкой, также усредненной по конфигурационному пространству,, Ока­ зывается, что ПКП удобен для описания одномодового поведения, а

PC — двухмодового.

Рассмотрим смешанный кристалл типа АВ1_ с Сс и ограничимся

для простоты одномерными системами. Последующая процедура мо­ жет быть легко применена к системам высшей размерности. Соглас­ но (2.1.27) и (1.3.52), функция Грина такой системы подчиняется урав­ нению

co2M(lx) G(lx, Vx\ со) =

 

 

= <5/А *' + Е

I"*") [<?(*"*", IV, со) -

G(lx, Vx', со)], (2.3.1)

/"*"(=Ик)

 

 

 

где полагается, что подрешетка

1 полностью занята атомами А, а

подрешетка х = 2 случайно занята атомами В, с

концентрацией 1 - с

и атомами С с концентрацией с .

 

 

Используя числа заполнения, определяемые (2.1.66), можно выра­ зить массы и силовые постоянные, входящие в (2.3.1), в следующем

виде:

 

(2.3.2)

M ( l x ) = £ v{lx)Mi(x),

i

 

 

ф(1х, iv ) =

Е rU fird W q*, V*').

(2.3.3)

 

W

 

Здесь индекс ;

пробегает значение от 1 до 3, и мы определяем

М’ (1) = МА, М2(1) = М3(1) = М1 (2) = 0, М2(2) = МВ, М3(2) = Мс , где

МА, и Мс обозначают соответственно массы атомов А, В и С.

Уравнение (2.3.3) основывается на предположении, что константа взаи­ модействия определяется только двумя атомами, которые оно связы­

вает. Разумно предположить, что Ф” '(1х, Г х ') = Ф,] У ~

)• Мы

будем рассматривать взаимодействие между атомами как сумму

ко­

роткодействующего отталкивания и дальнодействующей электростати­ ческой части. В двух- и трехмерных системах последняя часть приво­ дит к расцеплению LO- и ТО-колебаний. (см, разд. 1.5). Учет этого эф­

Динамика решетки кристаллов с беспорядком замещения

145

фекта важен для интерпретации оптических спектров таких систем (см., например, [228]).

Усредняя (2.3.1), получаем следующие функции Грина:

(2.3.4)

(2.3.5)

(2.3.6)

Здесь символ <• ••> означает усреднение по конфигурационному про­ странству. Функции Грина (2.3.5) и (2.3.6) называются условно усред­ ненными функциями Грина, потому что некоторые атомы при усредне­ нии остаются фиксированными. Подставляя теперь (2.3.2) и (2.3.3) в (2.3.1) и усредняя результат, получаем

со2Г 1Щ х) ctfklx, Гх', со) =

.7

*

(1х, 1'х', <о)1;

=

V'x") с’ с'» Ы !"(Гх", 1'х', со) -

 

 

(2.3.7)

здесь мы использовали (2.3.5), (2.3.6) и (2.1.96). Для рассматривае­

мого смешанного кристалла мы имеем с 1 = 1, с 2 = с3 = с 1 = 0.

«

, г

1

1

1

2

С2 ~ ^

~ с )> С2~ С '

 

 

 

 

Мы видим, что уравнение (2.3.7) связывает две различные функ­

ции Грина. Следовательно, это уравнение нельзя решить без какоголибо удобного приближения. Заметим, что иерархию связанных уравне­ ний для условно усредненных функций Грина (когда остаются фикси­ рованными один, два, . . . атома) можно получить из уравнения (2.3.1), если его умножить на ц и произведения величин т\и усреднить резуль­ тирующие уравнения. Для решения нам нужно будет получить из этой иерархии конечный замкнутый набор уравнений. С этой целью обычно условно усредненную функцию Грина с п + 1 фиксированными атомами заменяют на функцию с п фиксированными атомами. Самым грубым таким приближением является ПКП, которое состоит в полном пренеб­ режении различием между функциями Грина (2.3.4) - (2.3.6).

В ПКП уравнение (2.3.7) сводится к виду

со2М(х) $(1х, 1'х', со) = = <М**'+ £ Ф(1хЛпх")[$(Г,х",1'х\о>)--Ц1х,1,х\о>)\, (2.3.8)

10-297

146 Глава 2

где _

М (х) = E c j M

4 p ) ,

(2.3.9)

}

 

(2.3.10)

Ф(1х, 1'х') = £

JV) сМ ’

Сравнивая (2.3.8) и (2.3.1), мы видим, что (2.3.8) описывает кристалл с усредненными массами (2.3.9) и усредненными силовыми постоян­ ными (2.3.10).

Разложим теперь § в ряд Фурье

S(lx, Vx\ ш) = -i- Е У(**'. к><°)

,

(2.3.11)

™ к

где N означает число атомов в подрешетке, а х(1« ) описывает поло­ жение атома (/*)• Подставляя (2.3.11) в (2.3.8). мы получаем

Г {[со2Ж(*) +

Д (1 - <5,,"'5».'") Щ*, I"V")]

-

- £ (1 -

дН"дхх") Щ х , l" x " )

Ш) =

 

 

(2.3.12)

Полюса § или соответственно нули det| § ~ 1 1 определяют соб-

ственные частоты системы. Так как выражение в фигурных скобках в (2.3.12) - действительная матрица 2 x 2 , мы получаем из равенст­ ва det |§ “ 1|==0 две фононные ветви: акустическую и оптическую.

Итак, ПКП приводит к возбуждениям типа плоской волны, которые не подавляются беспорядком. Но это противоречит тому факту, что такие возбуждения не могут быть собственными модами неупорядо­ ченной системы, поскольку периодичность решетки отсутствует.

Исследуем сейчас длинноволновые колебания в ПКП. Для этого

положим в (2.3.12) k ■+0 и предположим для простоты, что взаимодей­

ствие существует только между ближайшими соседями.

Если опреде­

лить силовые постоянные FAB и рАС посредством

 

ФА*(Ц, 12) =

#*A(Z2,11) = ФА*(И, (I - 1) 2)

=

 

=

Ф* ((1 - 1) 2 ,11) г - F A*,

f = в, С,

(2.3.13)

то секулярное уравнение для длинноволновых колебаний примет вид

|й)2ЛГл — (1 — с) F AB cF AC (1 -

гЛ

+ c F AC

|

|(1 - с) FAB + cF AC

G>2[(1 -

с) M в +

сМ с] -

(1 - с) F** - c F AC

= 0.

 

 

 

(2.3.14)

Динамика решетки кристаллов с беспорядком замещения

147

Это уравнение имеет два корня:

сох = 0 (акустическое колебание),

- -[<“ -°> +■*">(i+

(2.3.15)

 

(оптическое колебание).

 

Следовательно, ПКП для любой концентрации дает одну длинноволно­ вую оптическую моду. Поэтому это приближение может быть исполь­ зовано только для объяснения (но не для предсказания) одномодового поведения. Оно представляет собой наиболее грубую интерполяцион­ ную схему во всей области концентраций.

Возвратимся к модели РССЭ. Снова начнем с уравнения (2.3.1). Заменяя в этом уравнении массы и силовые постоянные с помощью

(2.3.2) и (2.3.3), умножая получившееся уравнение на л; . и усредняя

{М)

его, мы имеем

(о*М>(х) s!(lx,l'x',w) = м » . ' +

Е Ф”"(1х,1"х")

х

*

Гх"(ф/х)

 

х

fSi’r(l"x", l'x', a>) -

sif (lx, l'x', со)].

(2.3.16)

<%,)>

L «*"

J

 

Здесь иерархическая система уравнений, упомянутая выше, разбивает^ ся с помощью следующего приближения в (2.3.16):

сfit”ДЬе, l’x', со) ян $i(lx, Vx', со).

(2.3.17)

Это представляет собой одноузельное приближение, так как атомы, оставленные фиксированными при усреднении, не "чувствуют”, явля­ ется ли другой атом также фиксированным или нет. То есть мы пре­ небрегаем соответствующими корреляциями между двумя или более атомами. С помощью (2.3.17) уравнение (2.3.16) приобретает вид

аРМЧх) $i(lx, l'x', at) = бц'дхк' 4- 27 ф”'(lx, l"x") с{ х

х

I'V'HWx)

 

X

[s ’r(l"x", l'x', со) - Si(lx, l'x', со)],

(2.3.18)

где также использовано (2.1.96).

Заметим, что в (2.3.18) входит только один элемент функции Гри­ на &iH(lx, I V , со) с / ~ = / и х " = х. Но при выводе приближений ПУТ и ПКП (см. разд. 2.3.4) помимо таких диагональных элемен­

148 Глава 2

тов мы должны также учитывать недиагональные элементы (см. [ 61]).

Величину $1 (1х, 1'х', со), можно разложить в ряд Фурье анало­ гично §(/*, Гх', со) в (2.3.11). Уравнение (вытекающее из (2.3.18)) для соответствующих коэффициентов Фурье имеет вид

27 {[*W M (1 - д ц - д т ~ ) < % № Г ( Ы , Г '*'")] M r -

-

2; (1 - <5/Г<$*к")

Гх")

J.

х

X # > " * ', А;, си) =

дкн>.

 

 

(2.3.19)

 

.. ^

с помощью равенства

 

 

Вводя функцию Грина gJJ

 

 

27 {. •J /'> > "*', *, (О) = M

r >

 

(2.3.20)

*"Г

 

 

 

 

 

где I .. .1

означает выражение в фигурных скобках в (2.3.19), полу­

чаем.

 

 

 

 

 

$(xx't k,co)=2J cj&(xxr9kt o))=2J cjg№(xx'9к, со) .

(2.3.21)

 

j

 

W

 

 

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

ны из уравнения det |{•••11 = 0, потому что полюса функций § и g”

совпадают. Так как {••• J - действительная матрица 3 x 3 , уравнение

det | |= 0 дает три ветви незатухающих фононов: одну акусти­ ческую и две оптические.

Секулярное уравнение, вытекающее из (2.3.19), совпадает с урав­

нением для собственных значений для кристалла с псевдоэлементар-

ной ячейкой, составленной из атомов А, (1 - с)В и с С (рис. 2.18).

Атомы В и С расположены в тех же самых узлах псевдоэлементарной ячейки. В уравнение движения для такой системы массы и силовые постоянные входят с весом, равным концентрации рассматриваемых атомов. Для более детального ознакомления с формулировкой такой модели псевдоэлементарной ячейки и строгим вычислением с ее по­ мощью дисперсионных кривых отошлем читателя к работе [ 401].

Рассмотрим теперь длинноволновый предел уравнения (2.3.19).

Предположим, что межатомные силы простираются вплоть до вторых соседей. С помощью

ФиЩ (I + 1) 2) = ФЩ1 +

1) 2,12) = Ф>?(12, (I -

1) 2) =

= ФЩ1 -

1) 2, i2) = -Fn'12,

j, j' = В, С, (2.3.22)

Динамика решетки кристаллов с беспорядком замещения

149

( 1-с)В

(1-с)В

 

2 (n -J)a

(2п-1)а

Зла

(2л+1)а

Рис. 2.18. Псевдоэлементарная ячейка смешанной двухатомной цепочки

АВ 1 _ с Сс. (Согласно [401 ] .)

 

 

и (2.3.13) мы получаем из (2.3.19) секулярное уравнение

ш2М к - (1 — с) F** -

cFAC (1 -

с) FAB

cFAC

FAB

со2Мъ — F** cFBC сFBC

FAC

(1 - c ) F BC

co2M c - F AC- ( l - c ) F BC

= 0.

 

 

(2.3.23)

Это уравнение идентично секулярному уравнению, полученному в ра­ боте [91] для модели PC, которая представляет собой описанную выше модель псевдоэлементарной ячейки для к = 0. Интересно отметить, что (2.3.23) не содержат силовых постоянных F ^ , FBB и Fc c для вторых соседей. Вклад этих силовых постоянных, возникающий из первого и второго члена в фигурных скобках (2.3.19), при k = 0 сокра­ щается. Это отражает тот факт, что в рассматриваемом приближении при к = 0 все атомы данного типа колеблются в фазе, как и в идеаль­ ном кристалле. Эти равные смещения и являются основой предположе­ ния при выводе (2.3.23) с помощью PC модели. Здесь этот результат вытекает из разбиения (2.3.17).

Уравнение (2.3.23) имеет три корня, причем один из них равен нулю (акустическая мода). Вблизи предельного случая кристалла одна из остающихся двух мод представляет собой колебание примеси, в то время как вторая является оптическим колебанием. Для системы, об­ наруживающей двухмодовое поведение, колебание примеси представ­ ляет собой локальное колебание,, когда состав смеси близок к одному

150

Глава 2

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

В первом случае, например, полагают, что колебанием примеси явля­ ется резонансная мода, а во втором случае, по-видимому, более под­ ходит комбинация резонансной и щелевой или локализованной мод.

Преимущества модели PC таковы: 1) двухмодовое поведение может быть объяснено, 2) может быть учтено изменение силовых постоянных, 3) возможна интерполяция на всю область концентраций. Недостатки же этой модели следующие: 1) все колебания описываются незатухаю - щими возбуждениями типа плоских волн независимо от того, являет­ ся ли данная мода примесной или зонной, 2) сравнение с эксперимен­

том дает для силовой постоянной F BC для взаимодействия атомов, следующих за ближайшими соседями,величину, вообще говоря, того же самого порядка, что и силовые постоянные для ближайших сосе­ дей. Для критического изучения модели PC отошлем читателя, напри­ мер, к работам [192, 249]; см. также примечание на с. 000.

2.3.3.'Колебания, 'активированные беспорядком,в инфракрасном и комбинационном спектрах

Как уже отмечалось в разд. 2.3.1, однофононные спектры инфракрас­ ного поглощения и комбинационного рассеяния в неупорядоченных системах иногда недвусмысленно показывают наличие мод, ативированных беспорядком. В частности, такие спектры могут отражать не­ которые особенности плотности колебательных состояний. Прежде чем рассматривать детально данное явление, приведем некоторые основ­ ные формулы для комбинационного рассеяния и инфракрасного погло­ щения.

Используя (П3.56), (П3.57), (2.1.45) и (2.1.22), мы получаем сле­ дующее выражение для сечения (нерезонансного) однофононного ком­ бинационного рассеяния:

d2g _ соа

da>s dQ

a)j

Е e«se .!/« .'!s 'K -< U i)e ;se.1.,

(2.3.24)

\771CJ ао 'аа'

 

Динамика решетки кристаллов с беспорядком замещения

151

ГД6

h

 

Z

Р *Ы *) X

 

 

= —— п(а>) £

 

 

71

м

т

 

 

 

X Im Gpjj(lx, 1к, со + is),

е -> +0,

(2.3.25)

И

[eAttt^ - I]"*1.

 

 

 

(2.3.26)

п(со) =

 

 

 

Здесь п(со) обозначает среднее число фононов с частотой со (см. 1.4.50));

еа и еа “* это а "я декартова компонента векторов поляризации соот­

ветственно падающего и рассеянного света; со. и со

- частоты падаю­

щего и рассеянного света, a P^V р (/к) - коэффициент разложения

первого порядка (по компоненте смещения W^ (/K)), оса' -

компоненты

тензора поляризуемости. Таким образом, для со = со

-со .

> 0 величи-

S

1

 

на /(со) из (2.3.25) пропорциональна числу фононов п(со) и описывает стоксову компоненту спектра. Если со отрицательна, то, вследствие

того что ImG(—|со |)=—ImG(| со]) (см. (2.1.22) и (2.1.45)), сечение рас­ сеяния пропорционально - п (- |со |) = п(| со |) + 1 и описывает антистоксову компоненту.

Коэффициент однофононного инфракрасного поглощения дается

выражением (см. (П3.20), (П3.21) или (2.2.10) соответственно)

« И

=

с

i

Е м^(Ы )

х

 

 

 

"

1x0

 

 

 

 

 

 

1к0

 

 

 

 

Xlm Gpp(lK, 1к, со + ге),

е -»• + 0 .

(2.3.27)

Здесь

 

* р (1к) -

коэффициент разложения (относительно

(/*)) пер­

вого порядка а -й компоненты электрического дипольного момента (т.е. эффективный заряд).

Таким образом, мы видим, что однофононное комбинационное рас­ сеяние и инфракрасное поглощение описываются выражением типа

О(со) сс - 27 Qa(lx) Qa'il'x') Im Gaa’(lK, Vx', со),

(2.3.28)

Ixa 1'н’а

где величина Q представляет собой коэффициент разложения первого порядка либо величины поляризуемости, либо дипольного момента. Согласно известной концепции усреднения по конфигурационному про­ странству, мы должны усреднить (2.3.28):

(О(со)) ос — 27 {Qa(lK) Qail'x') Im £««'(&> /V , (О));

(2.3.29)

1ха

Соседние файлы в папке книги