Рамановское рассеяние на интерфейсных модах

Для наблюдения ин­терфейсных мод с помощью рамановской спектроскопии волновой вектор рассеянно­го света должен иметь компоненту в плоскости структуры, так чтобы выполнялось условие θ >0. Однако, впервые интерфейсные моды наблюдались в геометрии рассеяния назад вдоль направления z. При этом необходи­мый волновой вектор в плоскости структуры возникал вследствие шероховатостей на индивидуальных интерфейсах (благодаря которым в плоскости нарушалась трансля­ционная инвариантность). На рис. 23 приведены спектры комбинационного рассеяния для AlAs-подобных мод в сверхрешетке GaAs/AlAs. Широкая полоса, расположенная посередине между частотами LO и ТО фононов, характерна для интерфейсных мод в симмет­ричных структурах (т.е. при d _A= d _B). Наблюдаемыйпик соответствует соотношению ε_A (ω_IF)= –ε_B (ω_IF). В асимметричном случае для получения длинноволновых решений при θ = π /2 следует использовать выражения (26) и (27). В отличие от симметричного случая, при d _A≠ d _B по­являются две различных моды. Им отвечают две полосы в спектрах, приведенных в нижней части рис. 23. Доминирует полоса с большим или с меньшим рамановским сдвигом в зависимости от условия d _A> d _B или d _A< d _B.

 

 

Рис. 23. Рамановские спектры в области частот оптических фононов объемного AlAs, полученные в резонансных условиях при ёь || es для СР GaAs/AlAs (А/В) с различными отношениями между толщиной слоев: 1) йд = с/в; 2) ds — Зс/д; 3) d\ = 3ds при Т = 10 К

 

Спектры на рис. 23 можно качественно объяснить как проявление простых электростатических интерфейсных мод типа, предсказанного уравнением (18), без учета взаимодействия с квантованными модами с нечетными значениями m.

 

 

Рис. 24. Сравнение расчетных (сплошные кривые) и измеренных (пунктир) рамановских спектров, обусловленных GaAs-подобными оптическими модами в МКЯ GaAs/AlAs с раз­личной толщиной слоев, в условиях выходного резонанса с переходами el-lhl при 10 К и параллельными поляризациями падающего и рассеянного света. Спектры нормализованы к одинаковой высоте и для наглядности сдвинуты вертикально друг относительно друга.

 

Проявление полярных колебаний в спектрах комбинационного рассеяния связано в первую очередь с электон-фононым взаимодействием, описываемым фрелиховским потенциалом. Гамильтониан фрёлиховского взаимодействия может быть выражен в виде скалярного электростатического потенциала, умноженного на заряд электрона. Теоретически этот потенциал можно вычислить, используя микроскопические ди­намические модели решетки с подходящими граничными условиями. Однако, на практике подобные расчеты занимают очень много времени, а получаемые результаты будут справедливы только для заданных в расчете ширин ямы и барьера. Поэтому желательно найти простые, хотя и приближенные, выражения для скалярного потенциала, которые можно было бы использовать для образцов с различной шириной ям. Было предложено несколько таких моделей. Они называются «макроскопическими», поскольку на начальном этапе в них предполагается, что образец представляет собой континуум. Основное различие между этими моделями заключается в подходе к граничным условиям, которые накладываются на оптические фононы на интерфейсах в квантовых ямах или сверхрешетках. В зависимости от применяемых граничных условий одни из них получили название «механических моделей», а другие — «диэлектрических континуальных моделей».

Грубо говоря, в «механических моделях» выдвигается требование, чтобы смещения квантованных LO фононов обращались в нуль на интерфейсе, даже если это приводит к нарушению на нем уравнений Максвелла. Пример картины смещений и электростатического потенциала квантованных LO-фононов, удовлетворяющей механическим граничным условиям, был приведен на рис. 17. Из рисунка видно, что электростатический потенциал ф не исчезает на интерфейсе. Для фонона, квантованного в среде А и обладающего не равным нулю волновым вектром q_x, можно представить скалярный потенциал как

 

ф (х, z)= ф _о ехр (iq_xx) cos (к_mz)для четного т,

ф (х, z)= ф _о ехр (iq_xx) sin (к_mz)для нечетного т.

 

 

Компонента электрического поля, параллельная интерфейсу {Е_х), имеет вид

 

E_x=–dф/dx=– (iq_x) ф (х, z)

 

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

В «диэлектрических континуальных моделях», как уже говорилось, в качестве отправной точки исполь­зуются уравнения Максвелла. На их основании получаются интерфейсные моды как часть решений уравнения Лапласа Ñ ²ф(r)= 0. Хотя такие модели и нарушают механические граничные условия (требующие, чтобы атомные смещения квантованных фононов были равны нулю на интерфейсе), однако для интерфейсных мод они являются довольно хорошим приближением к результатам микроскопических вычислений, поскольку в действительности смещения атомов становятся равными нулю только в непосредственной близости от интерфейса.

Макроскопическая модель, в которой предпринята попытка воспроизвести результаты микроскопических вычислений, была предложена Хуаном и Джу. Сделав модельные микроскопические вычисления динамики решетки для определения смещений атомов и электростатического потенциала, они обратили внимание на то, что диэлектрическая континуальная модель дает довольно хорошее приближение к результатам микроскопического расчета, за исключением нарушения механических граничных условий. Интерфейсные моды особенно хорошо описывались этой моделью. Для одновременного выполнения механических и максвелловских граничных условий необходимо, чтобы как ф, так и его производная dф/dz обращались в нуль на интерфейсе. Хуан и Джу заметили, что этого можно достичь, вычитая подходящую константу из ф с четной симметрией (по отношению к отражению в плоскости, проходящей через центр слоя) или подходящий член, линейный по z, из ф с нечетной симметрией.

 

 

Рис. 25. Сравнение атомных смещений электростатических потенциалов, связанных с самыми низкими квантованными и интерфейсными фононами в GaAs/AlAs, рассчитанных с помощью «макроскопических» моделей: a) диэлектрического континуума; б) механических граничных условий; в) Хуана-Джу; и г) в рамках микроскопической теории.

 

На рис. 25 потенциалы, связанные с квантованными LO фононами самого низкого порядка и с интерфейсными модами в сверхрешетке GaAs/AlAs, полученные с помощью трех обсуждавшихся выше макроскопических моделей, сравниваются с результатами микроскопической модели. Мы видим, что модель Хуана-Джу является наилучшей аппроксимацией микроскопической модели, за ней следует диэлектрическая континуальная модель. Хотя и существуют различия в вероятностях рассеяния, вычисленных с помощью этих моделей, однако различия между моделью Хуана-Джу и диэлектрической континуальной моделью исчезают в случае ям с малой шириной. Для таких ям основной вклад в вероятность рассеяния вносят интерфейсные моды, а они в этих двух моделях почти идентичны. Полученные в настоящее время экспериментальные резуль­таты лучше всего согласуются с предсказаниями модели Хуана-Джу.

 

 

Модель диэлектрического континуума для предельных фононов с k =0

 

Сравнение с экспериментом требует, как известно, рассмотрения исключительно центрозонных фононов, что соответствует синфазным колебаниям атомов в каждой элементарной ячейке. Поэтому важным является рассмотрение состояния сверхрешетки с однородной поляризацией в каждом слое. Именно такие состояния соответствуют атомным смещениям вдоль собственных векторов центрозонных фононов, и именно такие состояния имеет смысл рассматривать, когда речь идет о возбуждениях, проявляющихся в эксперименте (например, в комбинационном рассеянии). Соответствующую задачу можно рассмотреть в приближении диэлектрического континуума. В рамках такого подхода пренебрегают микроскопическими деталями кристаллической структуры и рассматривают вещество, из которого состоят слои сверхрешетки, как однородную упругую среду с заданной диэлектрической проницаемостью. Это модель, где использованы формулы, выведенные много раньше С. М. Рытовым [ЖЭТФ, 29, N5, 605 (1955)] при решении задачи о распространении радиоволн в слоистой среде. Рассмотрим полярные колебания упругого диэлектрического континуума в бесконечной среде, состоящей из плоских слоев толщины d _A и d _B с диэлектрической проницаемостью ε _A и ε _B периодически повторяющихся в направлении z, перпендикулярном границам раздела. Решая совместно уравнения движения и уравнения Максвелла, и усредняя полученные зависимости по объему, можно получить следующие выражения для эффективной диэлектрической постоянной в плоскости решетки ε_x,y и перпендикулярно к ней ε_z:

 

, (24)

, (25)

 

где .

Полученный результат можно пояснить следующим образом. Электрическое поле в слоях описывается векторами E _A и E _B и векторами индукции D _A= ε _A E _A и D _B= ε _B E _B. В случае поля, направленного параллельно оси x, из уравнения rot E =0 следует напряженности поля на границах раздела

.

И, следовательно, среднее значение индукции будет определяться выражением:

 

,

 

из которого непосредственно следует выражение (24).

В случае поля, направленного параллельно оси z, из уравнения div D =0 следует непрерывность индукции на границах раздела

 

.

 

И, следовательно, среднее значение напряженности будет определяться выражением:

 

,

 

из которого непосредственно следует выражение (25).

Выражение (24) соответствует случаю длинноволновых колебаний с поляризацией, перпендикулярно оси z, т.е. колебаниям типа E, а выражение (25) соответствует случаю длинноволновых колебаний с поляризацией, параллельной оси z, т.е. колебаниям типа А. Решение уравнения ε (ω)=0 соответствует LO модам, а частоты, при которых ε_z (ω) обращается в бесконечность, соответствуют TO модам.

Из формулы (25) непосредственно следует вывод:

 

ε_z (ω) → ε_A(ω)=0 или ε _B(ω)=0

Это означает, что частоты колебаний A(LО) в СР, совпадают с частотами колебаний A(LО), локализованных в слое A или B. Иными словами, спектр колебаний A (LO) в сверхрешетке состоит из двух мод, частоты которых совпадают с частотами A (LО) мод в объемной структуре.

С другой стороны, из формулы (24) следует,

 

Ε_xy (ω) → ε_A(ω)=∞ или ε _B(ω)=∞.

 

Это означает, что частота колебаний E (TO) в сверхрешетке совпадает с частотами таких колебаний, локализованных в каждом из слоев. Иными словами, в сверхрешетке спектр колебаний Е (ТО) состоит из двух мод, частоты которых совпадают с частотами колебаний Е (ТО) в объемном материале.

Таким образом, из уравнений (24 и 25), описывающих волны поляризации в слоистой среде, следует важный вывод: в сверхрешетках существуют длинноволновые оптически активные колебания, частоты которых совпадают с частотами A (LO) и Е (ТО) чистых кристаллов. Собственные векторы этих колебаний включают атомные смещения, локализованные в соответствующих слоях. Иными словами, в сверхрешетках моды A (LO) и Е (ТО) локализованы в отдельных слоях (confined modes) и не зависят от конкретной структуры сверхрешетки (относительной толщины слоев, например).

Иной характер имеют решения уравнений ε_xy (ω)=0 и ε_z (ω)=∞ которые соответствуют модам А (ТО) и Е (LO). Для колебаний Е (LO) из соотношения (24) непосредственно следует:

 

, (26)

 

а для колебаний A (TO) из соотношения (2) непосредственно следует:

. (27)

Зависимость ε_z (ω) для каждой ОС описывается зависимостью

(28),

причем, компоненты ε_xy (ω) определяются вкладом Е -мод, а ε_z (ω) определяются вкладом A -мод. Пользуясь этими соотношениями и решая уравнения (26 и 27) относительно ε (ω), можно вывести выражение для частот A (TO) и E (LO) в сверхрешетке через параметры спектра объемных структур и значения толщины слоев. Вместо вывода этих громоздких выражений, рассмотрим качественно характер зависимости этих решений от отношения d _A/ d _B. Прежде всего заметим, что отношение ε_A (ω)/ ε_B (ω) (см. рис. 19) принимает отрицательные значения в интервалах ТО1< ω <ТО2 и LO1< ω <LO2. Поскольку в каждом из этих интервалов отношение ε_A (ω)/ ε_B (ω) монотонно изменяется от нуля до бесконечности, то нетрудно понять, что при любом значении d _A/ d _B каждое из уравнений (26) и (27) имеет по одному решению в этих интервалах. В каждой паре таких решений есть «нормальные» компоненты и «аномальные». Поэтому положение интерфейсных мод болжно зависить от отношения периодов слоев А и В. Это хорошо видно из рис. 26.

По виду кривой, представленной на рис. 26 можно качественно определить характер зависимости частот нормальных и аномальных мод от соотношения толщины слоев в сверхрешетке. Рассмотрим, к примеру, моды A(TO), частоты которых определяются уравнением (4). Нормальное решение этого уравнения из интервала (TO1, TO2) при d ₁®0 стремится к TO2, а при d ₂®0 к TO1, т. е. частота этой моды изменяется пропорционально составу сверхрешеки. Аномальное же решение этого уравнения из интервала (LO1, LO2) при d ₁®0 стремится к LO1, а при d ₂®0 к LO2, т. е. частота этой моды изменяется обратно пропорционально составу СР. Аналогичные заключения можно вывести и относительно решений E(LO)-I и E(LO)-II.

Используя выражения (24 и 25) и зависимости (28), можно определить положение нулей и полюсов функции ε (ω)при произвольном направлении волнового вектора и, таким образом, исследовать угловую дисперсию нормальных колебаний сверхрешетке.

Рис. 26. Частотная зависимость отношения диэлектрических проницаемостей кристаллов GaN и AlN, из которой можно получить частоты интерфейсных мод сверхрешетки GaN/AlN.

 

 

Рис. 26 а. Рассчитанные КР-спектры СР (AlN)n/(GaN)m, симметрии E(TO) (а) и A(LO) (б) при разных значениях отношениях n/m, указанных на рисунках.

 

 

 

Рис. 26 б. Рассчитанные КР-спектры СР (AlN)n/(GaN)m, симметрии E(LO) (а) и A(TO) (б) при разных значениях отношения n/m, указанных на рисунках.

 

 

 

 

 

Рис. 26 с. Амплитуды смещений атомов в модах A(LO) (а) и E(LO) (б) при различных соотношениях толщины слоев m/n, указанных на рисунке. Значения рассчитанных частот указаны мелкими цифрами.

 

 

Фононы в нанокристаллах

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

Говоря о квантовых точках или нанокристаллах (нульмерные структуры), имеют в виду объекты, размеры которых составляют от 20 до 100Ǻ. Развитая технология получения полупроводниковых нанообъектов позволяет направленно варьировать размер таких нанокристаллов, а, следовательно, и энергетический спектр и их оптические свойства. В настоящее время это уже находит применение в приборах наноэлектроники. Все это вызывает повышенный интерес к нанообъектам как теоретиков, так и экспериментаторов.

Полупроводниковые нанокристаллы обычно исследуются в стеклянной матрице, в которой концентрация полупроводников группы А₂В₆, не превышает 0,1–1,5 %. Температура синтеза фосфатного стекла составляет около 1100С. Отжиг осуществляется при температуре стеклования Т_g в течение времени от 5 до 60 минут, затем образцы охлаждаются до комнатной температуры. Стекла до отжига всегда полностью бесцветные и микрокристаллы в них ни оптическими, ни методами рентгеноструктурного анализа не обнаруживаются. После отжига исследуемые образцы приобретают цвет от соломенно-желтого до темно-красного, что позволяет говорить о появлении нанокристаллов полупроводника. Отжиг стекол приводит к переконденсационному росту квантовых точек, причем средний размер нанокристаллов увеличивается при увеличении времени отжига. Исследование структуры таких образований производилось неоднократно.

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

 

 

Рис. 27. Электронно-микроскопическое изображение нанокристалла CdS в стеклянной матрице.

 

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

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

 

 

⇐ Предыдущая123456789Следующая ⇒

Поделиться:

Воспользуйтесь поиском по сайту: