- •Санкт-петербургский государственный университет физический факультет
- •С.В.Карпов фононы в кристаллах и гетероструктурах
- •Санкт-Петербургский государственный университет
- •Фононы в кристаллах и гетероструктурах
- •1. Симметрия кристаллов
- •1.1. Кристаллическая решетка
- •1.2. Элементы симметрии кристалла
- •Типы плоскостей скольжения
- •1.3. Сингонии и кристаллические классы
- •Кристаллические системы – сингонии
- •1 Тип решетки Браве
- •1 Тип кристаллического класса
- •1 Тип выбора частичной трансляции r для каждой операции группы r
- •Распределение кристаллических классов по сингониям
- •1.4. Классификация возбуждений в кристаллах
- •Неприводимые представления группы трансляций
- •1.5. Классификация возбуждений для фактор-группы
- •2. Ристаллический периодический потенциал
- •2.1. Общая модель твердого тела. Гамильтониан
- •2.2. Адиабатическое приближение
- •3. Зонные состояния периодических систем
- •3.1. Линейная моноатомная цепочка
- •Постановка решения в виде функции Блоха
- •3.2. Дисперсионные соотношения (закон дисперсии)
- •Двухпроводная электрическая линия
- •2. Акустические колебания в системе резонаторов
- •3. Связанные маятники
- •Электромагнитные волны в атмосфере
- •5. Многоатомная линейная цепочка
- •Однородный упругий стержень и стержень с периодической плотностью
- •Волны де-Бройля
- •3.3. Уравнение Матье и зонная структура
- •3.4. Фазовая и групповая скорость волн в диспергирующей среде
- •4. Фононы в идеальных кристаллах
- •4.1. Линейная двухатомная цепочка
- •4.2. Колебания трехмерной решетки
- •4.3. Обратная решетка и зона Бриллюэна
- •4.4. Ход ветвей колебаний в зоне
- •4.5. Расчеты колебаний кристаллов
- •Как известно, коэффициенты Lkl являются элементами матрицы, для которой выполнено:
- •4.6. Функция распределения плотности частот
- •Особенности функции g(), обусловленные различными критическими точками
- •5. Полярные колебания в кристаллах
- •5.1. Продольные и поперечные акустические колебания
- •Поэтому:
- •5.2. Поперечные и продольные оптические колебания
- •5.3. Соотношения Лиддейна-Сакса-Теллера
- •Отсюда следует, что
- •5.4. Реальные состояния. Эффект "запаздывания". Поляритон
- •Первые два уравнения, как известно, дают
- •6. Квантовомеханическое представление колебаний
- •6.1. Нормальные колебания.
- •6.2. Фононы
- •6.3. Гармонический осциллятор
- •Решение стационарного уравнения Шредингера
- •6.4. Операторы рождения и уничтожения фононов
- •6.5. Ангармонический осциллятор и кристалл
- •6.6. Фонон-фононные взаимодействия
- •7.1. Низкоразмерные 3d, 2d, 1d, 0d системы
- •7.2. Фононы в объемных и ограниченных структурах
- •7.3. Размерно-ограниченные кристаллические среды.
- •7.4. Приближение упругого континуума.
- •7.5. Рамановское рассеяние на сложенных акустических фононах (folding phonons)
- •7.6. Приближение механического континуума.
- •7.7. Рамановское рассеяние на квантованных оптических фононах
- •7.8. Приближение диэлектрического континуума
- •7.9. Рамановское рассеяние на интерфейсных модах
- •8.1. Модель упругого континуума. Лэмбовская мода
- •8.2. Модель механического континуума
- •8.3. Модель диэлектрического континуума
- •8.4. Расчеты колебательных спектров нанокристаллов
- •Оглавление
- •I. Симметрия и структура кристаллов
- •II. Кристаллический периодический потенциал
- •III. Зонные состояния периодических систем
Особенности функции g(), обусловленные различными критическими точками
-
Тип точки
Обозна-чение
поведение g() вблизи c
c
c
Минимум
Po
0
(c2-2)1/2
Седловая точка
P1
(c2-2)1/2
const
Седловая точка
P2
const
(c2-2)1/2
Максимум
P3
(c2-2)1/2
0
Вид функции плотности частот вблизи этих точек показан на рис. 42a. Существует важная теорема Ван-Хова, утверждающая, что в трехмерном случае спектр частот колебаний каждой ветви должен содержать, по крайней мере, точки минимума (P0) и максимума (P3) и по трем критическим седловым точкам каждого типа (P1 и P2). Важно, что производная на высокочастотном конце спектра должна стремиться к – (критическая точка P3).
Существование особенностей функции g() является следствием дискретности кристаллической решетки. В трехмерном случае изолированная критическая точка приводит к разрыву лишь производной g(), а не самой функции. Кроме аналитических критических точек существует другие критические точки, которые могут давать более сильные особенности на функции распределения плотности частот. Это относится к тем случаям, когда в критической точке вторая производная по волновому вектору равна нулю. В частности, критические точки объемоцентрированной кубической решетки, для которой дисперсионная формула имеет вид:
,
где i=1/2aoki, и a – постоянная решетки. Уравнения, определяющие критические точки, имеют вид:
cos1 cos2 sin3 = 0
cos1 sin2 cos3 = 0
sin1 cos2 cos3 = 0.
Решение этой системы показывает, что точки (0,0,0), (,,), (/2,/2,/2) являются критическими, причем первые две из них аналитические, последняя – неаналитическая. Эта точка (/2,/2,/2) является точкой пересечения трех взаимно-перпендикулярных плоскостей постоянной частоты 1=/2, 2=/2, 3=/2 и обуславливает появление у функции g() особенности вида ln2|2–1/2o2|.
Рис. 35. Появление особенностей Ван-Хова в простой кубической решетке. а) зона Бриллюэна простой кубической решетки и симметричные точки зоны Г, X, M и R. б) дисперсионные зависимости одной ветви, нарисованные для случая направлений (100), (110) и (111), и имеющие приблизительный функциональный вид 1–cos(kiai) для каждого направления ki. Справа показана получающаяся для этого случая функция плотности состояний g().
На рис. 36 показан спектр собственных частот кристалла меди отдельно для трех ветвей. Он получен путем решения векового уравнения для 5600 значений волнового вектора k, равномерно распределенных по неприводимой 1/48 части зоны Бриллюэна для различных отношений силовых постоянных взаимодействий между атомами. Очевидно, что такое большое количество вычисленных частот дает спектр, в котором отчетливо видны особенности Ван-Хова. Однако в ряде случаев возможна компенсация особенностей функции плотности частот. В трехмерном случае, например, возможна компенсация особенностей, обусловленных минимумом и седловой точкой P2, а также максимумом и седловой точкой P1 в различных ветвях.
Существование особенностей функции распределения собственных частот было показано в расчетных спектрах для определенны моделей кристаллов. Однако, окончательно было не ясно, присуще ли наличие критических точек только рассмотренным частным моделям или их наличие носит общий характер. Ответ на этот вопрос дал Ван-Хов, который показал, что существование критических точек в семействе плоскостей постоянной частоты в k-пространстве, а, следовательно, и существование особенностей функции распределения частот является необходимым следствием периодичности решетки в обычном пространстве.
Рис. 36. а) Спектр трех ветвей собственных частот кристалла меди, вычисленный для модели с центральным взаимодействием между ближайшими и следующими за ними соседями. б) Спектр кристалла натрия, рассчитанный для 24.576.000 значений волнового вектора в первой зоне Бриллюэна.
Топологическое обоснование существования аналитических критических точек в двумерном случае состоит в следующем. На рис. 42б показано несколько элементарных ячеек обратного k-пространства. Поскольку дисперсионная зависимость (k1,k2) непрерывна и периодична в пространстве волновых векторов k, то в каждой ячейке она должна иметь, по крайней мере, один максимум и один минимум. Пусть расположение максимума в каждой ячейке отмечено знаком , а положение минимума функции (k) – черными кружками. Если максимумы A и B в соседних ячейках соединить произвольной кривой, то на ней будет один минимум, т.е. точка, в которой значение функции (k) принимает меньшее значение, чем в соседних точках. Аналогичные точки есть и на любых других кривых, соединяющих рассматриваемые максимумы A и B. Геометрическое место всех таких точек образует непрерывную кривую, проходящую через абсолютные минимумы C и D в элементарно ячейке обратного пространства. Ясно, что на этой кривой минимумов будет точка, в которой функция (k) принимает наибольшее значение. Эта точка должна быть седловой. Действительно, если двигаться вдоль кривой 4 от C к D, то она будет соответствовать относительному максимуму, а если двигаться по кривой 2 от A к B – относительному минимуму. В двумерной случае в элементарной ячейке обратного пространства имеется две таких седловых точки. Аналогично эвристическое рассмотрение можно провести и в случае трехмерной решетки, что приводит к результатам, установленным Ван-Ховом. В трехмерном случае седловые точки будут двух типов, а дисперсионная функция будет иметь по три точки каждого типа.