Поступила в редакцию: 20 мая 2002 г.
Выставление онлайн: 19 апреля 2003 г.
При исследовании затухания фотолюминесценци (ФЛ) линий D0,h и D0,x в эпитаксиальных слоях n-GaAs обнаружено, что в затухании присутствует медленная компонента tau~10-5 s [1]. На рис. 1 [1,2] представлена кривая затухания ФЛ линии D0,h. Время быстрого спада tau0 соответствует времени жизни фотовозбужденных носителей tau0=<10-8 s, которое характерно для оптических переходов, формирующих указанные выше линии [2,3]. Проведенные нами эксперименты [4,5] позволили сделать вывод, что медленное затухание интенсивности ФЛ обусловлено существованием метастабильных состояний, захватывающих дырки. Формирование длительного послесвечения линии излучения D0,h происходит в результате безактивационного выброса дырки в валентную зону из метастабильного центра, которая затем рекомбинирует с электроном мелкого донора. Основными характеристиками медленного затухания излучения являются время спада интенсивности taus и величина интеграла по времени от интенсивности излучения под кривой медленного затухания Itau. Время медленного спада taus определяется процессами высвобождения дырок из метастабильных состояний, при этом оно является постоянным для всех исследованных образцов и практически не зависит от внешних воздействий. Интеграл Itau пропорционален количеству дырок, высвободившихся из метастабильных состояний, и, естественно, зависит от концентрации метастабильных центров и их заселенности. Заселение метастабильных центров осуществляется во время действия импульса возбуждения и определяется интенсивностью этого возбуждения (темпом генерации) и эффективностью механизма перехода носителей на эти центры, который связан с величиной сечения захвата sigmaMS. Полное количество захваченных на метастабильные центры дырок после окончания импульса возбуждения пропорционально также концентрации этих центров NMS. Поскольку интегральная интенсивность длительного послесвечения (Itau) определяется концентрацией заселенных в результате возбуждения метастабильных центров, она отличается для разных образцов (рис. 1, b) вследствие того, что отличается NMS. Экспериментально количество захваченных дырок может быть оценено из светосуммы, излученной образцом после окончания импульса возбуждения. Эта светосумма может быть определена как площадь под кривой затухания фотолюминесценции. Она показана на рис. 1, a штриховкой и равна интегралу Itau= I(t)dt. Процессы захвата носителей в метастабильные состояния и их освобождение непосредственно не связаны с процессом излучательной рекомбинации (D0,h; D0,x), следовательно, при анализе экспериментальных данных необходимо разделять влияние внешних воздействий на излучательную рекомбинацию и процессы заселения метастабильных центров. Для того чтобы учесть разнообразные факторы, влияющие на интенсивность линий (D0,h) и (D0,x), мы нормировали этот интеграл на величину интенсивности линий (D0,h) и (D0,x) в момент возбуждения (BMS=Itau/Iex, где Iex - интеграл интенсивности свечения в момент импульса возбуждения). Таким образом, определенная величина BMS не будет зависеть от изменения собственно процесса излучения линии D0,h и, следовательно, BMS может служить мерой, отражающей заселенность метастабильных состояний. Для выяснения особенностей механизмов заселения метастабильных центров исследовано поведение BMS при различных внешних воздействиях. Исследование зависимости величины BMS от интенсивности возбуждения было выполнено ранее [4]. Проведенный в настоящей работе анализ экспериментальных данных совместно с результатами модельного численного расчета позволил описать поведение BMS во всем диапазоне интенсивностей возбуждающего света. Это позволило оценить концентрацию метастабильных центров (NMS~1013 cm-3). Заселение центров определяется (в значительной степени) энергетическим положением уровня захвата, который может находиться как в запрещенной, так и в разрешенной зоне. В последнем случае такие состояния будут действовать как источники резонансного рассеяния [6] свободных носителей. [!tb] a - кривая затухания излучения линии D0,h при T=2 K в одном из образцов n-GaAs. Заштрихованная область обозначает интеграл Itau= I(t)dt, определяемый медленным высвечиванием. На вставке показан типичный вид спектра n-GaAs в экситонной области при T=2 K, где LPB и UPB - нижняя и верхняя ветви поляритона соответственно; линия A0,x - излучение экситона, связанного на нейтральном акцепторе. b - кривые затухания излучения линии D0,h (T=2 K) двух различных образцов n-GaAs. Цель данной работы - выяснение особенностей процесса заселения метастабильных состояний. Для этого изучено поведение величины BMS долговременного затухания (отражающей заселенность метастабильных состояний в данном образце) в зависимости от температуры (Tb), величины одноосного сжатия (Pd) и напряженности магнитного поля (Bf). Эксперимент Исследовались образцы эпитаксиального n-GaAs с Nd-Na=1013-1015 cm-3. Слои выращивались на подложках с ориентаций < 100> различными методами: жидкофазной эпитаксией (LPE), газофазной эпитаксией в хлоридной системе (HVPE) и эпитаксией MOC-гидридным способом. Все образцы, использовавшиеся в экспериментах, имели достаточно высокую интенсивность фотолюминесценции в экситонной области спектра. Люминесценция возбуждалась импульсным YAG-лазером (lambdaex=530 nm, tau=10-7 s) либо полупроводниковым лазером (lambdaex=759 nm, tau=10-6-10-7 s). Спектры и кинетика ФЛ регистрировались с использованием двойного дифракционного спектрометра ДФС-52 методом время-коррелированного счета фотонов. Температурные измерения осуществлялись в диапазоне 2-25 K с использованием криостата УТРЕКС. Зависимость спектров затухания ФЛ (D0,x; D0,h) от величин напряженности внешнего магнитного поля излучалась при температурах 1.8 и 4.2 K, а одноосного сжатия - только при 4.2 K. Регистрация спектров и кинетики в магнитном поле осуществлялась в геометрии Фарадея (магнитное поле перпендикулярно поверхности образца). Величина напряженности магнитного поля менялась в диапазоне от 0 до 4 T. Экспериментальные результаты и обсуждение [!b] Зависимость величины BMS от температуры. [!tb] Изменение BMS при T=4.2 K в зависимости от величины одноосного сжатия ( a) и напряженности магнитного поля ( b). Экспериментальная зависимость BMS от температуры (Tb) приведена на рис. 2. Из рисунка видно, что величина BMS уменьшается с ростом температуры приблизительно на порядок. Исследование зависимости поведения BMS (при T=4.2 K) от величины одноосного сжатия (Pd) и напряженности магнитного поля (Bf) (рис. 3, a и b соответственно) показало, что интеграл BMS также уменьшается (приблизительно в 2 раза) по мере увеличения Pd и Bf. Из рис. 2, 3 видно, что характер изменения BMS в зависимости от Tb,Pd и Bf, хотя и имеет разные функциональные зависимости, качественно идентичен - BMS монотонно уменьшается. Важным обстоятельством является то, что уменьшение величины BMS (отражающей поведение глубокого локализованного центра) происходит при малых значениях (и малом диапазоне изменения) Tb,Pd и Bf. Переходя к обсуждению, заметим, что температурная зависимость сечения захвата на различные центры изучалась подробно для различных полупроводниковых соединений [7-9] и, как показано в этих исследованиях, определяется многими факторами. При этом следует заметить, что характер температурной зависимости сечения sigmaMS=F(Tb) чувствителен к параметрам центра и реализующемуся в конкретных условиях механизму захвата носителей на центр. При термоактивационном механизме сечение захвата (а следовательно, BMS~sigmaMS) растет с увеличением температуры [7], а при каскадном захвате на притягивающий центр [8] сечение уменьшается с ростом температуры, что наблюдалось и в нашем случае. Однако применимость каскадного механима захвата для описания приведенных выше экспериментальных результатов вызывает сомнения, поскольку время высвобождения из метастабильных состояний не зависит от температуры. Кроме того, исследования кинетики в наносекундном диапазоне и ионизации связанного экситона (линия D0,x) в зависимости от температуры [10] показали, что в области температур 6-8 K в поведении BMS наблюдается заметная особенность. Принимая во внимание эти факты, а также формальную независимость сечения захвата в каскадной модели [8] (изменение плотности энергетических уровней при малом изменении Pd и Bf маловероятно), мы использовали для анализа экспериментальных данных по зависимости BMS от Tb,Pd и Bf модель туннельного механизма заселения глубоких центров [11,12] (природа потенциального барьера для приведенного далее феноменологического рассмотрения не столь важна и поэтому здесь не обсуждается). В рамках этой модели заселение метастабильных центров осуществляется в результате туннельного перехода носителя из свободного состояния Ef, принадлежащего зоне (E=hk2/2m*), на уровень захвата метастабильного центра EMS. Вероятность туннелирования определяется в общем виде выражением [13] Wt~exp(-(2m)1/2/h[U(x)-EMS]1/2 dx), где m - масса частицы, U(x) - потенциал барьера, EMS - энергия, для которой рассчитывается вероятность туннельного перехода. Количество свободных носителей, туннелирующих из Ef на уровень метастабильного центра, будет определяться заселенностью состояния (реальное заселение метастабильного центра возможно, если заселено состояние зоны с энергией Et=EMS). [!b] Изменение спектрального положения линии D0,h (T=4.2 K) по мере увеличения одноосного сжатия ( a) и напряженности магнитного поля ( b). [!tb] Уменьшение заселенности N(Ef)=f(E)rho(E) в зоне при значении энергии Et=EMS для четырех различных значений величины одноосного сжатия Pd по мере их увеличения (af1<af2<af3<af4). [!tb] Зависимость величины BMS от напряженности магнитного поля при T=2 K. В этой связи рассмотрим причины изменения заселенности метастабильных состояний (BMS) при изменении внешних воздействий (Tb,Pd,Bf), которые определяются реакцией свободных носителей на эти внешние воздействия и тем самым зависят от функции распределения носителей в разрешенной зоне. Можно полагать, что при использовавшихся в эксперименте сравнительно небольших значениях Pd и Bf состояние локального центра практически не изменяется, а все изменения связаны с функцией распределения свободных носителей. Следует заметить, что увеличение Pd и Bf приводит (несмотря на различную природу воздействия) к одному результату, а именно к "эффективному" увеличению ширины запрещенной зоны в Gamma-точке, на изменение величины которой указывает одинаковый сдвиг всех линий спектра излучения (Fx,D0,x и D0,h) в сторону больших энергий (этот сдвиг показан для линии D0,h на рис. 4). В результате сдвига края зоны изменяется заселенность на энергии Ef=EMS. Заселенность при этом определяется выражением N(E)=rho(E)f(E), где rho(E) - плотность состояний, а f(E) - функция распределения носителей в разрешенной зоне. Для качественной иллюстрации этого процесса функция распределения свободных носителей f(E) принималась нами в виде распределения Максвелла (T=4.2 K). Положение максимума функции f(E) соответствует энергии Ef. Изменение полученной таким образом функции заселенности N(E) при значении энергии Ef=EMS для четырех последовательно возрастающих значений Pd показано на рис. 5. Этот рисунок иллюстрирует уменьшение заселенности в точке зоны с энергией Ef=EMS при увеличении одноосного сжатия и напряженности магнитного поля. Действительно, как видно из рис. 5, по мере увеличения Pd и Bf и, следовательно, максимума функции распределения f(E) заселенность точки N(Ef) будет монотонно уменьшаться (точки af1,af2,af3,af4 на рис. 5). Необходимо отметить два аспекта, следующие из приведенной выше модели туннельного перехода. 1) Заселение центров реально, если энергия туннельного перехода EMS (отвечающая уровню метастабильного состояния) попадает в интервал энергий, равный полуширине функции распределения носителей в зоне. При этом осуществляется резонансное туннелирование. 2) Величина BMS (а следовательно, и заселенность метастабильных состояний) будет уменьшаться при увеличении Pd и Bf только в том случае, если при фиксированной температуре T=2 и 4.2 K и Pd=0, Bf=0 энергия, отвечающая максимуму функции распределения (EF), больше энергии, соответствующей уровню метастабильного состояния (EMS). Резонансный характер туннелирования находит подтверждение в зависимости интеграла BMS от величины напряженности магнитного поля при T=2 K, которая показана на рис. 6. Как видно из рисунка, при увеличении напряженности магнитного поля на фоне монотонного спада интеграла BMS наблюдается заметная особенность при напряженности магнитного поля H=1.5 T. Причиной увеличения заселенности (величины интеграла BMS) при определенном значении напряженности магнитного поля является возникновение уровней Ландау. Действительно, в этом случае выражение для энергии носителей в зоне без учета спиновой составляющей принимает вид E=h2k2z/2m+(n+1/2)homegac, где omegac=eH/m*. Здесь omegac - циклотронная частота, e - заряд электрона, H - величина напряженности магнитного поля, m* - эффективная масса носителя. Из этого выражения следует, что при увеличении внешнего магнитного поля край зоны смещается в сторону больших энергий, вследствие чего и наблюдается, как отмечено выше, монотонный спад заселенности MS при температуре 4.2 K. Появление на фоне монотонного спада локального максимума при B=1.5 T обусловлено модификацией плотности состояний - возникают критические точки для подзон Ландау, которые определяются условием d En(kz)/d kz=0. В результате при совпадении энергии максимума (EF) модифицированной плотности состояний с энергией метастабильного уровня EMS туннелирование возрастает, вследствие чего увеличивается заселенность этого состояния (возрастает интеграл BMS). Величина локального максимума интеграла BMS (и даже само его наличие) зависит от соотношения полуширины функции распределения носителей по энергиям, которая определяется температурой (kBT), и энергетического интервала (homegaB) между уровнями Ландау. Действительно, максимум наблюдается при температуре T=2 K и величине магнитного поля B=1.5 T, однако при температуре 4.2 K, когда величина kT больше homegaC, максимума нет и спад заселенности MS происходит только монотонно. Анализ приведенных выше экспериментальных данных зависимости величины BMS от Pd,Bf и Tb позволяет в первом приближении оценить границы энергии, в пределах которых находится уровень захвата EMS метастабильного центра. Возможность такой оценки основана на следующем: в случае Et<EMS по мере приближения Ef к энергии EMS количество носителей, туннелирующих в состояние с энергией EMS, будет монотонно возрастать и, следовательно, BMS будет увеличиваться, а поскольку с увеличением Pd и Bf наблюдалось монотонное уменьшение интеграла BMS, энергия уровня EMS, на который происходит захват, оказывается меньше энергии, соответствующей максимуму распределения носителей в зоне (EMS<Ef). Только в этом случае при сдвиге функции распределения носителей в сторону высоких энергий (при приложении давления или магнитного поля) количество носителей с энергией, равной энергии уровня захвата на EMS (N(Ef=EMS)), будет уменьшаться (точки af1,af2,af3,af4 на рис. 5). Величина BMS максимальна при af1 и минимальна при af4 (рис. 3). Таким образом, задача определения энергии уровня захвата (EMS) сводится к определению энергетического интервала, в котором изменяется максимум функции распределения и его положение Ef относительно края зоны при (Pd,Bf)=0 и температурах эксперимента T и 4.2 K. Величина изменения энергетического положения (T=const) максимума функции распределения f(E) (af), обусловленная увеличением Eg, определяется из величины сдвига линий излучения (рис. 4) и приблизительно равна 1 meV. Определение реального значения Ef затруднено вследствие того, что функция распределения является неравновесной и эффективная температура носителей (T*) заметно отличается от температуры решетки. Определение вида неравновесной функции распределения является очень трудной задачей. Для грубой оценки мы использовали данные работы [14], в которой приведено значение эффективной температуры для экситонов T*~7 K (0.6 meV). На основании сделанных таким образом оценок можно предположить, что величина EMS находится в пределах нескольких meV от края зоны. Таким образом, совокупность экспериментальных данных, отражающих эволюцию величины BMS при увеличении Pd,Tb и Bf, позволяет предположить существование уровня захвата в разрешенной зоне, который отстоит от края зоны на несколько meV. Можно полагать, что его заселение осуществляется в результате резонансного туннельного обмена носителями между подсистемой свободных носителей и резервуаром метастабильных локализованных состояний. Кроме того, зависимость интеграла BMS от Pd и Bf позволяет качественно определить низкоэнергетический хвост неравновесной функции распределения экситонов.
- А.В. Акимов, А.А. Каплянский, В.В. Криволапчук, Е.С. Москаленко. Письма в ЖЭТФ 46, 35 (1987)
- А.В. Акимов, В.В. Криволапчук, Н.К. Полетаев, В.Г. Шофман. ФТП 27, 310 (1993)
- D. Bimberg, H. Munzel, A. Steckenborn, J. Christen. Phys. Rev. B 12, 7788 (1985)
- А.В. Акимов, Ю.В. Жиляев, В.В. Криволапчук, Н.К. Полетаев, В.Г. Шофман. ФТП 25, 713 (1991)
- В.В. Криволапчук, Н.К. Полетаев, Л.М. Федоров. ФТП 28, 310 (1994)
- Б. Ридли. Квантовые процессы в полупроводниках. Мир, М. (1986). С. 175
- C.H. Henry, D.V. Lang. Phys. Rev. B 15, 989 (1977)
- В.Н. Абакумов, В.И. Перель, И.Н. Яссиевич. ФТП 12, 3 (1978)
- В.Н. Абакумов, И.А. Меркулов, В.И. Перель, И.Н. Яссиевич. ЖЭТФ 89, 1472 (1985)
- В.В. Криволапчук, Н.К. Полетаев. ФТП 32, 307 (1998)
- R. Kubo. Phys. Rev. 86, 929 (1952)
- R. Kubo, Y. Toyozawa. Progr. Theor. Phys. 13, 160 (1957)
- Л.Д. Ландау, Е.М. Лифшиц. Квантовая механика. Наука, М. (1974). С. 215
- В.М. Ботнарюк, Ю.В. Жиляев, В.В. Россин, Т.В. Россина, В.В. Травников. ФТТ 28, 201 (1986)
Подсчитывается количество просмотров абстрактов ("html" на диаграммах) и полных версий статей ("pdf"). Просмотры с одинаковых IP-адресов засчитываются, если происходят с интервалом не менее 2-х часов.
Дата начала обработки статистических данных - 27 января 2016 г.