- •Санкт-петербургский государственный университет физический факультет
- •С.В.Карпов фононы в кристаллах и гетероструктурах
- •Санкт-Петербургский государственный университет
- •Фононы в кристаллах и гетероструктурах
- •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. Зонные состояния периодических систем
6.6. Фонон-фононные взаимодействия
Учет ангармонизма при рассмотрении отдельного осциллятора приводит только к перенормировке собственных значений системы (и, естественно, к изменению правил отбора). При рассмотрении же набора ангармонических осцилляторов, которые имею место в кристалле, задача становится более сложной. Набор ангармонических осцилляторов никаким выбором новых координат невозможно свести к совокупности независимых (невзаимодействующих) мод. Перекрестные члены в гамильтониане, которые были исключены в случае набора гармонических осцилляторов выбором нормальных осцилляторов и тем самым вводят взаимодействие между отдельными осцилляторами (фононами). Вид матричных элементов в обоих случаях выглядит следующим образом:
Поскольку в ангармоническом случае точно диагонализировать гамильтониан невозможно, используют нормальные координаты, полученные в гармоническом приближении. Это позволяет выделить основные несвязанные члены, входящие в матричные элементы. Добавки H′nn на диагонали фактически указывают на перенормировку собственных значений (т.е. энергий фононов). Однако, в случае набора ангармонических осцилляторов эта поправка комплексна. Действительная ее часть характеризует изменение энергии реального физического фонона, а мнимая часть – указывает на конечное время жизни. Недиагональные члены матричных элементов характеризуют взаимодействие с другими фононами.
Гамильтониан набора ангармонических осцилляторов имеет вид H′= H0+ H′, где H0 – гармоническая часть, а H′ представляет собой кубический член и член четвертого порядка при разложении потенциальной энергии в ряд по смещениям:
Если использовать при учете ангармонизма нормальные (т.е. гармонические) координаты Qj(k), то гамильтониан H0 будет иметь квадратичную форму, а добавка H′ вид:
Ясно, что в этом случае кубический ангармонический член связывает три нормальные координаты и, следовательно, характеризует процесс взаимодействия трех фононов Qj1(k1), Qj2(k2), Qj3(k3) с волновыми векторами k1, k2, k3 , принадлежащие ветвям j1, j2, j3 соответственно. Член четвертой степени описывает взаимодействие четырех фононов Qj1(k1), Qj2(k2), Qj3(k3), Qj4(k4).
Поскольку общий вид нормального колебания представляет собой функцию Блоха с множителем exp[i(krn)], то в произведении Qj1(k1)*Qj2(k2)*Qj3(k3) для трехфононного процесса будет входить множитель
exp[i(k1rn)]*exp[i(k2rn)]*exp[i(k3rn)],
который из-за условий периодичности кристалла должен быть инвариантен относительно добавлению к вектору rn произвольного вектора ri:
exp[i(k1+k2+k3,ai)]=1 и, значит, k 1+k2+k3=K m ,
где Km=b1m1+b2m2+b3m3 – целочисленный вектор обратной решетки.
Аналогично сказанному, для четырехфононных процессов будет выполнено:
k1+k2+k3+k4=Km .
Эти выражения представляют собой закон сохранения импульса k фонона: при взаимодействии сумма импульсов взаимодействующих фононов сохраняется с точностью до целочисленного вектора обратной решетки Km. Отличие этого закона сохранения от классического связано с наличием трансляционной симметрии кристалла. Именно поэтому импульс фонона (как и других частиц) называется квазиимпульсом, а сами частицы носят название квазичастиц. В частном случае, когда Km =0, процесс носит название N-процесса (normal), в случае Km0 – U-процесса ( umklap).
Члены, описывающие взаимодействия трех и четырех фононов, можно записать через оператора рождения и уничтожения фононов и использовать формализм, развитый для отдельного ангармонического осциллятора. Для представления взаимодействия фононов часто используются диаграммы, указывающие временное развитие процесса взаимодействия. Помимо выполнения закона сохранения квазиимпульса должен быть выполнен закон сохранения энергии:
ћ(k1)=ћ (k2)+ћ (k3) .
Первый процесс на рис. 43 представляет собой процесс ангармонического распада фона с энергией ћ(k) и квазиимпульсом ћk на два других фонона с импульсами ћk1 и ћk2. Возможен также и обратный процесс аннигиляции двух фононов с образованием третьего. Эти процессы определяются кубическим членом гамильтониана H′3 и описываются первым порядком в теории возмущения - член <n′|H′3|n>. Возможно рассмотреть вклад в фонон-фононное взаимодействие и члена второго порядка теории возмущения, включающего все виртуальные процессы, в результате которых при распаде одного фонона появляются два (рис. 43). Очевидно, имеется бесконечное число таких процессов. На рис. 44 показаны некоторые из возможных четырехфононных процессов. Они классифицируются как по порядку взаимодействия (число фононов), так и по величине вкладов. Кроме того, необходимо указание на относительную величину параметра , появляющегося в разложении ангармонического потенциала H′=H′3+2H′4+3H′5. На рисунке указаны также некоторые многофононные процессы высоких порядков и указан порядок вклада этих процессов по возмущению и по параметру . При учете взаимодействия во втором порядке теории возмущения необходимо правильно суммировать по всем промежуточным (виртуальным) состояниям, приводящим к конечному результату. Для учета членов 2-го порядка в 4-х фононном процессе (рис. 44) нужно просуммировать следующие матричные элементы:
где k1 - промежуточное виртуальное состояние, для которого не требуется выполнения закона сохранения энергии. Выполнение же закона сохранения квазиимпульса требуется в каждом взаимодействии (т.е. в каждой точке диаграммы).
Рис. 43. Трехфононные процессы, которые могут давать вклад в изменение числа фононов n(k,j). Km – целочисленный вектор обратной решетки. Процессы 1б и 2б не отличаются от процессов 1a и 2a. Процессы 1 и 4 увеличивают число n(k,j) фононов, а процессы 2 и 3 – уменьшают его.
Рис. 44. Некоторые четырехфононные процессы, которые могут давать вклад в изменение числа фононов n(k,j). Km – целочисленный вектор обратной решетки.
YII. ФОНОНЫ В ГЕТЕРОСТРУКТУРАХ И СВЕРХРЕШЕТКАХ