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

f_5y40we-arphclqbl0c

.pdf
Скачиваний:
18
Добавлен:
15.02.2016
Размер:
1.46 Mб
Скачать

101

Аналогично (8.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

3

q D

 

k 2dk

 

 

ν

(ω) =

 

 

 

 

δ[ω s p k ].

(9.4)

 

2π

2

 

 

 

 

p =1 0

 

 

 

 

 

Совершая замену переменных ωp = s p k , получаем

 

 

1

3

 

s p q D

 

 

 

ν(ω) =

s p3

 

 

δ[ω ωp ]ωp2dωp .

(9.5)

 

2π 2

 

 

 

 

p =1

 

0

 

 

 

 

Для ω <<θD аргумент δ-функции обращается в ноль на интервале интегрирования.

Согласно свойству δ-функции,

ω2

 

f (a), если ω1

< a < ω2,

.

f (ω)δ(a ω)dω =

 

 

или a < ω1

ω1

 

0, если a > ω2

 

Поэтому для ω <<θD

 

 

 

 

 

 

 

 

1

 

3

3

 

 

 

ν(ω) =

ω2

s p3 =

 

ω2 .

(9.6)

 

2π 2s 3

2π 2

 

p =1

 

 

Вид плотности фононных состояний для реального закона дисперсии в трехмерном кристалле изображен на рис.9.1.

ν

ω

Рис.9.1. Плотность фононных состояний в трехмерном кристалле

На нем видны изломы - особые точки. Они называются осо-

102

бенностями Ван-Хова и возникают при значениях ω, соответствующих особым точкам закона дисперсии ωp (k ) - таким точкам,

в которых обращается в нуль групповая скорость фононов или, другими словами, одновременно обращаются в нуль частные

производные

ωp (k )

,

ωp (k )

,

ωp (k )

.

k x

k y

k z

 

 

 

 

Пусть точка k 0 - особая точка для закона дисперсии ωp (k ) . Разложим функцию ωp (k ) вблизи k 0 с точностью до квадратич-

ных по k k 0

членов.

При этом разность ωp (k ) ωp (k 0 ) будет

представлять

собой

квадратичную

 

форму по

переменным

kx k0 x ,

 

ky k0 y ,

 

kz k0 z . Выбором ориентации ортогональной

системы координат приведем ее к диагональному виду. Тогда

 

 

 

 

 

 

ω

(k ) =ω

(k

) +α (k

x

k

0x

)2

+

 

 

 

 

 

 

 

 

p

p

0

1

 

 

 

,

(9.7)

 

 

 

 

 

 

+α2 (ky k0 y )2 +α3 (kz k0 z )2

 

 

 

 

 

 

 

 

 

 

где a = 1

2ωp (k )

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

2

 

ki2

 

k =k0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проклассифицируем эти особые точки в зависимости от знаков коэффициентов αi .

Случай положительных αi отвечает минимуму закона дисперсии для оптических ветвей.

Случай отрицательных αi соответствует максимуму за-

кона дисперсии.

Если два коэффициента αi из трех отрицательны, а один

положителен, то это седловая точка первого типа, а если наоборот: один отрицателен, а два положительны - то седловая точка второго типа.

Рассмотрим поведение функции ν(ω) вблизи максимума закона дисперсии одной из ветвей. Подставляя (9.7) в (9.3), получаем

 

103

 

ν p (ω) = ∫

dqxdqydqz

δ[ω ωp (k0 ) α1qx2

α2q2y α3qz2 ], (9.8)

 

 

(2π)3

 

где q = k k0 .

Если ω > ωp (k 0 ) , то в силу неравенства αi <0 аргумент δ- функции в ноль не обращается и ν p (ω) = 0.

Если же ω < ωp (k 0 ), то сделав замену переменных

~

=

 

αi

 

1/ 2

qi ,

(9.9)

 

 

qi

 

 

 

находим

 

 

 

 

 

 

 

~

 

 

~

~

 

 

 

 

 

 

 

 

 

~

2

 

 

 

 

 

 

 

dqx dq y dqz

 

 

 

 

 

 

 

 

ν p (ω) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ[ω ωp (k

0 ) +q

] =

 

(2π)3

 

α α α

3

 

1/ 2

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ ~

2

 

 

 

 

 

 

~2

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

qdq

δ[ω ωp (k 0 ) +q

] =

(2π)2

 

α α α

3

 

 

1/ 2

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

[ωp (k 0 ) ω]1/ 2

.

 

 

 

 

 

 

 

 

 

 

 

 

(2π)2

 

α α α

3

 

1/ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

Таким образом, вблизи максимума

 

 

 

 

 

 

 

0, если ω > ωp (k 0 ),

 

[ωp (k 0 ) ω]1/ 2

 

ν p (ω) =

, еслиω < ωp (k 0 ).

 

 

 

 

 

 

(2π)2

 

α α α

 

1/ 2

 

 

3

 

 

 

 

1 2

 

 

(9.10)

(9.11)

Функция ν p (ω) имеет корневую особенность (см. рис.9.1). Производная функции ν p (ω) слева при ω ωp (k0 ) 0 стремится к -

, а справа при ω ωp (k0 ) + 0 равна нулю.

В точке минимума оптического закона дисперсии ситуация полностью аналогична рассмотренной, и для получения ν p (ω) надо

изменитьзнаки перед ω и ωp (k0 ) в правойчасти выражения (9.11).

104

Исследуем теперь седловую точку первого типа. Пусть для определенности α3 >0, а α1, α2 <0.

После замены (9.9) приведем выражение (9.8) к виду

 

 

 

~

 

 

 

~

~

 

 

 

 

 

 

 

 

 

 

 

~2

~2

~2

 

 

 

 

 

 

 

 

dqxdqy dqz

 

 

 

 

 

 

 

 

 

 

 

 

ν p (ω) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ

[ω ωp (k0 ) + qx + qy qz ].

(9.12)

(2π)3

 

α α α

3

 

1/ 2

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

~2

~2

1/ 2

и выражая d

3~

Вводя новую переменную q = (qx

+ qy )

q

в цилиндрических координатах, получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

~

 

 

~

 

 

 

 

 

 

~2

~2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q dq dqz

 

 

 

 

 

 

 

ν p

(ω) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ[ω ωp (k0 )

+ q

qz

] =

(9.13)

 

(2π)2

 

α α α

3

 

1/ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

qгр

 

~

qгр

 

~2

 

 

 

~2

~2

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dqz

d(q )δ[ω ωp

(k0 )

+ q

qz

]

 

8π 2

 

α α α

 

1/ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

гр

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где qгр – величина порядка qD.

Если ω < ωp (k0 ) , то аргумент δ-функции обращается в ноль

~2

для любых

 

~

 

 

<< qгр. При этом

 

 

 

 

при интегрировании по q

 

qz

 

 

 

 

1

 

 

 

 

 

qгр ~

 

ν p =

 

 

 

 

 

 

 

 

 

 

 

dqz = const .

(9.14)

(2π)2

 

α α α

 

1/ 2

 

 

 

 

3

 

0

 

 

 

 

 

1

2

 

 

 

 

 

 

 

Более точное рассмотрение показывает, что формула (9.14) является оценкой по порядку величины, так как при q~z ~ qгр уже не

применимо разложение (9.7). Однако ясно, что при ω < ωp (k0 )

νp (ω) не имеет особенностей и слабо зависит от ω при малых

ωωp (k0 ) :

ν p (ω) = const +O[ω ωp (k0 )],

(9.15)

где O(γ) – совокупность членов первой и более высоких степеней по параметру γ.

Если ω > ωp (k0 ), аргумент дельта-функции в (9.13) обраща-

105

ется

 

~2

, только если

в нуль при интегрировании по q

 

 

~

 

1/ 2

. Следовательно,

 

 

 

 

 

 

qz

> [ω ωp (k0 )]

 

 

 

 

1

 

 

 

qгр ~

 

 

 

[ω ωp (k0 )]1/ 2

 

ν p =

 

 

 

 

 

 

 

 

dqz

= const

 

 

 

 

 

 

, (9.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

1/ 2

 

 

(2π)2

 

α1α2α3

 

1/ 2 [ωω p (k0 )]1/ 2

(2π)2

 

α1α2α3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то есть ν p (ω) имеет корневую особенность при ω > ωp (k0 ):

 

ν p (ω)

 

 

 

 

 

 

 

→ −∞;

 

ν p (ω)

 

 

 

= const .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

ω

 

 

 

 

 

 

 

 

ωωp (k0 )+0

 

 

 

ωωp (k0 )0

 

Седловая точка второго типа отличается от рассмотренной выше заменой знаков перед ω и ωp (k0 ) в правой части выраже-

ния (9.13).

Для нее

 

 

 

 

1/2

 

 

 

const

[ωp (k

0 ) ω]

, если

 

 

1/2

 

 

 

 

(2π)2

α1α2α3

 

 

 

 

 

 

νp (ω) =

 

 

 

 

 

 

если ω

 

 

 

 

 

 

 

const +O[ωp (k0 ) ω],

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

ν p (ω)

 

 

 

 

= const ;

ν p (ω)

 

ω

 

ωωp (k0 )+0

 

ω

ω <ωp (k0 );

>ωp (k0 );

(9.17)

 

→ +∞

ωωp (k0 )0

Графически зависимость ν p (ω) в седловых точках обоих

типов изображена на рис.9.1.

Согласно теореме Ван-Хова, в трехмерном случае закон дисперсии обладает, как минимум, одной седловой точкой каждого типа.

9.2. Локальные колебания

До сих пор мы изучали колебания идеальной решетки. Но реальные кристаллы содержат примеси, дислокации, плоские де-

106

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

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

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

Локальные колебания возникают также, если силовые постоянные взаимодействия примеси с атомами матрицы превосходят силовые постоянные матрицы, а ее масса примерно равна или меньше массы атомов матрицы.

Если же частота колебаний примесного атома лежит в области частот колебаний матрицы, то возникает так называемое квазилокальное колебание. Амплитуда колебаний на такой моде также спадает по мере удаления от примесного атома, но выходит вдали от примеси на постоянное значение (рис. 9.26).

A A

x

x

а б

Рис.9.2. Зависимость амплитуды локальных (а) и квазилокальных (б) колебаний от расстояния до примесного атома

9.3. Методы исследования фононных спектров

Основным методом исследования фононных спектров явля-

107

ется неупругая дифракция нейтронов на кристалле. Для проведения эксперимента выбираются "тепловые" нейтроны с дебройлевской длиной волны порядка межатомного расстояния d. Энергия таких нейтронов порядка 0,1 эВ. При неупругом взаимодействии с кристаллической решеткой (в отличие от упругой Брэгговской дифракции) нейтрон поглощает или испускает фонон. При этом его энергия изменяется на энергию фонона (примерно 0,01 эВ), а импульс - на величину импульса фонона (с точностью до вектора обратной решетки).

Измеряя направление и энергию вылетевших из кристалла нейтронов и сравнивая их с первоначальными значениями, можно определить из законов сохранения энергию и импульс поглощенного (или испущенного) фонона и найти закон дисперсии фононов (рис.9.3).

Рис.9.3. Законы дисперсии фононов в алмазе

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

108

9.4. Оценка величины параметра ulj,s / d

Теперь настала пора ликвидировать старую задолженность и оценить величину параметра разложения потенциальной энергии

кристалла в ряд по смещениям, то есть параметра ulj,s / d (пара-

метра ангармонизма).

 

Характерная частота колебаний атомов ω

порядка

ω ~ (G / M )1/ 2 , где М - масса иона, а силовая постоянная G, как

уже отмечалось, имеет порядок величины G ~ Eат / d 2 .

 

Следовательно

 

ω ~ d 1(Eат / M )1/ 2 .

(9.18)

С другой стороны, характерная энергия колебаний составляет величину Gu2/2, где u - амплитуда смещения атома в процессе колебаний.

В области низких температур (Т<<θD) это энергия нулевых колебаний, то есть

Gu2 ~ ω ~ d 1(Eат / M )1/ 2 .

(9.19)

Отсюда

 

u ~ d( 2 / d 2MEат)1/ 4 .

(9.20)

Отношение /d по порядку величины совпадает с импульсом pB на границе зоны Бриллюэна. Как будет показано позднее, характерная энергия электронов c импульсом pB также имеет

атомный масштаб: pB2 / me ~ Eат (me - масса электрона). Подстав-

ляя в (9.20)

2 / d 2 ~ p2

~ m E

ат

, получаем

 

 

 

B

e

 

 

 

 

 

u ~ d(m / M )1/ 4

,

(9.21)

 

 

 

 

e

 

 

или окончательно

u / d ~ (me / M )1/ 4 <<1.

 

 

 

(9.22)

109

Характерное значение параметра u/d изменяется от нескольких сотых для тяжелых ионов до одной десятой для гелия.

В области высоких температур (Т>>θD) характерная энергия колебаний атома становится порядка температуры.

Теперь

Gu2 ~ T ~ ω(T / ω)

(9.23)

Таким образом, величина u/d увеличивается, по сравнению с (9.22) в (T / ω)1/ 2 раз. Учитывая, что характерное значение ω ~ θD , получаем

(m u / d = e(me

/ M )1/

4 (T /θD )1/ 2 , если T >>θD ;

/ M )1/

4 , еслиT

(9.24)

θD.

Сомножитель (T /θD )1/ 2 увеличивает значение параметра

u/d в области высоких температур, но вплоть до точки плавления кристалла u/d остается много меньшим единицы.

Таким образом, разложение (6.7) справедливо во всей области существования кристаллического состояния.

10. КИНЕТИЧЕСКОЕ УРАВНЕНИЕ БОЛЬЦМАНА

10.1. Вводные замечания

До сих пор мы рассматривали равновесные свойства фононной подсистемы, используя функцию распределения (8.1). Такое рассмотрение справедливо далеко не всегда. Под действием внешних сил или вследствие изменения условий окружающей среды (например, температуры) система может быть выведена из равновесия. Для описания неравновесной ситуации и процессов релаксации к равновесию нам необходимо уравнение, описывающее эволюцию неравновесной функции распределения при изменяющихся внешних условиях и (или) под действием внешних сил, приложенных к системе. Таким уравнением является кинетическое уравнение Больцмана.

110

Мы получим его для классической идеальной системы частиц (не важно каких), а затем рассмотрим условия его применимости к реальным квантовым системам (рамки квазиклассического приближения).

Исследуем систему частиц, выведенных из равновесия посредством внешнего воздействия. Пусть частицы системы слабо взаимодействуют между собой, так что потенциальной энергией этого взаимодействия можно пренебречь по сравнению с их кинетической энергией (смотри §7.2). Для того, чтобы такая система релаксировала, необходимы процессы обмена энергией и импульсом. Таковыми являются столкновения частиц. Будем предполагать, что столкновения происходят практически мгновенно, то есть характеристики частиц изменяются скачком от старых значений к новым. Это справедливо, если время столкновения намного меньше времени свободного пробега между столкновениями, то есть при изучении разреженных систем. Именно такие системы и являются слабонеидеальными.

Ограничимся рассмотрением точечных объектов, не имеющих внутренней структуры. Состояние такой частицы во внешнем поле полностью характеризуется координатами и импульсом и его можно задать точкой в шестимерном пространстве координат (3 измерения) и импульсов (3 измерения), называемом фазовым пространством. Точку в фазовом пространстве, соответствующую данной частице, называют изобразительной точкой.

Если бы мы учли взаимодействие частиц друг с другом, то для определения состояния данной частицы мы должны были бы задать координаты и импульсы окружающих частиц, с которыми она взаимодействует. То есть шестью переменными дело не ограничилось бы. Поскольку все частицы оказались бы связанными друг с другом посредством взаимодействия, то потребовалось бы 6N-мерное фазовое пространство, где N - число частиц в системе.

Функцией распределения называют функцию семи переменных: времени, координат и импульсов или, другими словами, - времени и координат фазового пространства. Она задает плотность распределения изобразительных точек, отвечающих не-

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]