- •Санкт-петербургский государственный университет физический факультет
- •С.В.Карпов фононы в кристаллах и гетероструктурах
- •Санкт-Петербургский государственный университет
- •Фононы в кристаллах и гетероструктурах
- •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. Зонные состояния периодических систем
7.5. Рамановское рассеяние на сложенных акустических фононах (folding phonons)
Акустические фононы в сверхрешетках получаются из акустических фононов в объемных материалах в результате усреднения акустических дисперсионных соотношений для каждого q и далее путем сложения ветвей столько раз, сколько потребуется в пределах мини-ЗБ. В случае Ge2Si2 (см. рис. 57) появляется только одна щель на границе ЗБ. Отметим, что в общем случае эта граница соответствует продольному вектору k с величиной π/d/, в то время как в объемном кристалле для направления [001] соответствующий вектор имеет величину 2π/a0, т.е. вдвое больший, чем на рис. 57.
В объемном случае обычно предполагается, что волновой вектор фононов, активных в рамановском рассеянии, очень мал, т.е. только фононы, относящиеся к центру зоны, могут быть раман-активными. Обычное объяснение этому утверждению в пределах дипольного приближения довольно очевидно: длина волны света много больше, чем характеристические длины в материале (т.е. радиус экситона и/или постоянная решетки). Ограничение k ~ 0 остается справедливым для короткопериодных сверхрешеток, но должно быть изменено при достаточно большой величине d. Определим приведенную величину волнового вектора рассеяния k для рассеяния назад как
k = [4πn/λ]/[π/λ]=4πn/λ .
Для типичных значений показателя преломления п ~ 3, 5 и длины волны лазера λ=500 нм, k=1 для d=36 нм, что соответствует приблизительно 130 монослоям GaAs. В этом случае рассеяние происходит на фононах, относящихся к границе мини-ЗБ. Таким образом, меняя d или λ, можно охватить весь диапазон изменений приведенного волнового вектора и даже достичь границы мини-ЗБ. Наблюдаемые частоты в рамановском рассеянии для геометрии назад, когда волновой вектор фонона равен 4πn/λ, легко получить из анализа рис. 59, где показана сложенная акустическая ветвь и значение k.
Экспериментальные спектры комбинационного рассеяния приведены на рис. 60, из которого видно, что наблюдаются низкочастотные линии, позволяющие с высокой точностью определять период сверхрешетки.
Рис. 60. Соответствие дисперсионных сложенных ветвей и спектра комбинационного рассеяния сверхрешетки GaAs/AlAs.
7.6. Приближение механического континуума.
Квантованные неполярные оптические моды
Квантование неполярных оптических мод математически выражается наложением граничного условия, требующего, чтобы амплитуды колебаний обращались в нуль в непосредственной близости от границ между слоями А-В и В-А. При этих условиях мы получаем А-подобные и В-подобные квантованные оптические моды, зависимость частот которых от дискретных векторов k можно найти из соответствующих дисперсионных соотношений для объемного случая:
Это уравнение можно использовать для LO и для ТО мод. Соответствующие смещения атомов меняют знак при переходе от одного атомного слоя к соседнему, как это требуется для оптических мод, с величинами, определяемыми огибающими функциями
um(z) = cos kmz, m=1, 3, 5.. ; um(z) = sin kmz, m=2, 4, 6,... ,
где z указывает на положение атомной плоскости вдоль оси сверхрешетки, а значение z=0 находится в середине слоя.
Рис. 61. Картина смещений для AlAs-подобных квантованных LO мод в сверхрешетке (GaAs)5/(AlAs)5; kх = 0, kг→∞. а) – слой GaAs, б) – слой AlAs. Крупным кружками (голубыми) показана величина смещения в модах LO1, LO2, LO3 для атомов As, а маленькими (зелеными) кружками — для атомов Ga и А1.
Обычно написанное выражение является хорошим приближением для эффективного волнового вектора. Однако, в частном случае короткопериодных сверхрешеток (когда период состоит только из нескольких атомных слоев) следует понять, будет ли приближение, в котором на границе каждого слоя материала волновая функция в точности равна нулю, наилучшим. Расчет показывает (см. рис.61), что разумнее считать огибающую функцию равной нулю на первых атомах «другого сорта» (т.е. на атомах А для В-подобных колебаний и наоборот). Колебательные амплитуды ведут себя в соответствии с упомянутым выражением, но эффективные векторы kт нужно сделать слегка меньше, чем это указано: амплитуды исчезают не на номинальной границе слоя AlAs, а около первого слоя Ga, вне границы.
Это связано с тем, что граница образована слоями As, которые должны все еще колебаться на частотах AlAs, в то время как Ga, имеющий гораздо большую массу, чем Аl, не будет колебаться. Последний факт можно учесть, заменив толщину слоя dA, в выражении на dA,B+ао /4, где ао — объемная постоянная решетки, которая соответствует четырем атомным слоям для материалов типа алмаза и цинковой обманки.