Пространственно-временные и спектральные характеристики нестационарных волновых процессов в неоднородных конденсированных средах

Штыгашев, Александр Анатольевич

Пространственно-временные и спектральные характеристики нестационарных волновых процессов в неоднородных конденсированных средах тема автореферата и диссертации по физике, 01.04.07 ВАК РФ

Штыгашев, Александр Анатольевич АВТОР

доктора физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Новосибирск МЕСТО ЗАЩИТЫ

2012 ГОД ЗАЩИТЫ

01.04.07 КОД ВАК РФ

Диссертация по физике на тему «Пространственно-временные и спектральные характеристики нестационарных волновых процессов в неоднородных конденсированных средах»

Автореферат диссертации на тему "Пространственно-временные и спектральные характеристики нестационарных волновых процессов в неоднородных конденсированных средах"

НОВОСИБИРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

На правах рукописи

Штыгашев Александр Анатольевич

Пространственно-временные и спектральные

характеристики нестационарных волновых процессов в неоднородных конденсированных

средах

01.04.07 - физика конденсированного состояния

АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора физико-математических наук

г 4 ЯНВ 2013 005048622

Новосибирск - 2012

005048622

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Новосибирский государственный технический университет».

Научный консультант: доктор физико-математических наук, доцент, Пейсахович Юрий Григорьевич

Официальные оппоненты:

Вальков Валерий Владимирович, доктор физико-математических наук, профессор, Федеральное государственное бюджетное учреждение науки Институт физики им. Л. В. Киренского Сибирского отделения Российской академии наук, г. Красноярск, зав. лабораторией

Слабко Виталий Васильевич, доктор физико-математических наук, профессор, ФГАОУ ВПО «Сибирский федеральный университет», г. Красноярск, профессор

Сатанин Аркадий Михайлович, доктор физико-математических наук, профессор, ФГБОУ ВПО «Нижегородский государственный университет», г. Нижний Новгород, профессор

Ведущая организация: ФГБОУ ВПО «Новосибирский национальный исследовательский государственный университет» (Новосибирский государственный университет), г. Новосибирск

Защита состоится « 5~ » апреля 2013 г. в ^^ часов на заседании диссертационного совета Д 003.055.02 в Федеральном государственном бюджетном учреждении науки Институте физики им. Л.В. Киренского Сибирского отделения Российской академии наук, расположенном по адресу: 660036, Красноярск, Академгородок 50, стр. 38, в конференц-зале главного корпуса.

С диссертацией можно ознакомиться в библиотеке Федерального государственного бюджетного учреждения науки Институт физики им. Л.В. Киренского Сибирского отделения Российской академии наук

Автореферат разослан « -У » января 2013 г.

Ученый секретарь диссертационного совета, доктор физико-математических наук

Общая характеристика работы

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

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

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

Основное содержание диссертации составляет исследование волновых процессов в твердотельных пространственно-неоднородных структурах. Такие структуры, как правило, имеют выделенное направление, поэтому при описании волновых процессов используется приближенное или точное разделение переменных двух- или трехмерных волновых уравнений и сведение их к одномерным моделям. Одномерные волновые уравнения являются основными моделями в диссертационном исследовании.

Целью диссертационной работы является исследование стационарных и нестационарных квантовых и классических линейных волновых процессов, происходящих в пространственно-неоднородных структурах, посредством разработки математических моделей и алгоритмов расчета пространственно-временных и спектральных характеристик волновых полей и процессов переноса.

Научная новизна результатов диссертационной работы состоит в следующем:

1. Изучен процесс заселения и распада квазистационарных состояний в многобарьерных структурах при прохождении квантового волнового пакета. Найдена зависимость от длины решетки для времени жизни квазистационарных состояний и времени прохождения пакетов. Эта зависимость

различна в случаях, когда резонансное состояние и спектральный центр пакета лежат в середине, либо на краю полосы пропускания.

2. Исследован коллапс резонансов полной прозрачности в двойной периодической зеркально-симметричной сверхрешетке.

3. Предложено объяснение двух режимов временной зависимости вероятности нахождения электронов в системе «молекула сенсибилизатора - слой наночастиц» фотоэлектрохимического элемента, исследовано влияние буферного слоя наночастиц на процесс переноса электрона и среды на электронные переходы в молекуле.

4. Построена теория нестационарной фотоэмиссии через поверхностную трехбарьерную гетероструктуру, основанная на методе матрицы плотности и учитывающая неупругие процессы рассеяния электронов. Рассчитаны параметры осцилляции фотоэмиссионного тока.

5. Для монохроматической циркулярно-поляризованной ультразвуковой волны, падающей на многослойную среду, состоящую из чередующихся ферромагнитных и немагнитных материалов, получены спектры коэффициентов прохождения и отражения в области магнитоакустического резонанса (МАР) без учета и с учетом затухания упругих и спиновых возбуждений в системе. Показано, что взаимодействие упругой и спиновой волн может приводить к полному отражению ультразвука. Изучены характеристики магнитоупругих квазистационарных состояний ферромагнитной пластины при изменении её толщины и внешнего магнитного поля. Исследована динамика прохождения ультразвукового импульса через ферромагнитный слой в области МАР.

6. Получено спектральное уравнение и проведены численные расчеты спектра андреевских состояний в несимметричном ¿¡^-переходе при разных радиусах поверхностей Ферми в слоях структуры.

Научная и практическая значимость работы. Разработанный в диссертации подход позволяет расширить круг решаемых задач при моделировании волновых процессов в многослойных пространственно-неоднородных структурах и увеличить точность их описания.

Алгоритмы и результаты расчета состояний электронов в решетках конечной длины во внешних полях могут быть применены при описании баллистических и диссипативных токов через контакты и гетероструктуры, при разработке многослойных фотодиодов, фотоприемников электромагнитных излучений и квантовых каскадных лазеров.

Численное моделирование процесса переноса фотоэлектрона через буферный слой наночастиц в коллектор фотоэлектрохимической ячейки позволяет продемонстрировать и оценить влияние на перенос заряда размеров и дисперсии размеров наночастиц. Результаты моделирования могут учитываться при разработке эффективных солнечных элементов на органических красителях.

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

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

Научные положения и результаты, выносимые на защиту.

1. Результаты моделирования распространения волновых пакетов в многослойных структурах в отсутствие и при наличии внешних силовых полей. Диаграммы локальной устойчивости стационарных функций элек-

трона в решетке в сильном электрическом поле. Асимптотические и численные оценки характеристик распространения пакетов, формирования и распада квазистационарных состояний в зависимости от длины решетки, спектрального состава пакета и его положения по отношению к полосам пропускания. Показано, что время жизни пороговых квазистационарных состояний пропорционально третьей степени длины решетки, а для квазистационарных состояний средней части полосы пропускания время жизни пропорционально длине решетки. Распадный ток из многобарьерных структур может осциллировать с частотами, определяемыми энергиями смешиваемых квазистационарных состояний.

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

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

4. Основанная на методе матрицы плотности теория нестационарной импульсной фотоэмиссии и стационарной фотоэмиссии из фотокатода с поверхностной гетероструктурой. Из этой теории следует, что фотоэмиссионный ток через поверхностную трехбарьерную гетероструктуру может осциллировать.

5. Решение системы связанных волновых уравнений для упругих и спиновых волн в многослойной структуре. Результаты исследования процессов рассеяния стационарной ультразвуковой волны и ультразвукового импульса на ферромагнитной пластине и магнитной сверхрешетке в области магнитоакустического резонанса (МАР) при произвольных условиях

закрепления спинов на поверхности. Показано, что при резонансных частотах вблизи МАР происходит полное отражение ультразвука (антирезо-нансы Фано), вызываемое взаимодействием упругой волны со спиновыми волнами, локализованными в магнитных слоях. Изменением параметров исходного волнового пакета можно целенаправленно возбуждать в магнитных слоях квазистационарные состояния и вторичные импульсы, обусловленные интерференцией и перекачкой энергии между упругой и спиновой подсистемами.

6. Численно-аналитическое исследование спектра андреевских состояний в несимметричном 5^52 переходе при разных поперечных импульсах Ферми в слоях. Показано, что при увеличении разности поперечных импульсов Ферми фазовая зависимость спектра энергий стационарных андреевских состояний ослабевает.

Апробация результатов работы. Основные результаты исследований докладывались и обсуждались на Всесоюзном симпозиуме с международным участием по новым методам лазерной спектроскопии в низкотемпературных средах (Таллин, 1987); Всесоюзном совещании по люминесценции молекул и кристаллов (Таллин, 1987); III Российско-корейском симпозиуме по науке и технике (Новосибирск, 1999); 32 Всероссийском совещании по физике низких температур (Казань, 2000); международных конференциях по актуальным проблемам электронного приборостроения (Новосибирск, 2004, 2006, 2008, 2010); совещании «Кремний-2004»(Иркутск, 2004); IX Российской конференции по физике полупроводников (Новосибирск-Томск, 2009); международной конференции с элементами научной школы для молодежи «Нанофизика и наноэлектроника» «Мезоскопические структуры в фундаментальных и прикладных исследованиях» (Новосибирск, 2010); IV Байкальской международной конференции «Магнитные материалы. Новые технологии» (Иркутск, 2010); X Всероссийской научной конференции «Краевые задачи и математическое моделирование» (Новокузнецк, 2010); Российской конференции и школе по актуальным проблемам полупроводниковой нанофотоэлектроники «Фотоника-2011» (Новосибирск, 2011); между-

народной зимней школе физиков-теоретиков «Коуровка-34» (Новоуральск, 2012); научных сессиях и семинарах профессорско-преподавательского состава НГТУ и научных учреждений СО РАН.

Публикации. По теме диссертации опубликовано 28 научных работ, из них - 19 статей из списка ВАК, 9 работ — в материалах и трудах российских и международных научных и научно-технических конференций и совещаний, в сборниках научных трудов и монография.

Личный вклад автора диссертации в получении выносимых на её защиту положений заключался в постановке физических проблем для их численного моделирования и получении ряда аналитических решений указанных задач; в разработке вычислительных алгоритмов и проведении вычислительных экспериментов; в анализе физического содержания представленных результатов численного моделирования; формулировке окончательных выводов исследований.

Структура и объем диссертации. Диссертация состоит из введения, шести глав, заключения, приложения и списка литературы. Общий объем работы составляет 317 страниц машинописного текста, включая 86 рисунков, 2 таблицы и список литературы из 254 наименований.

Содержание работы

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

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

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

ных структурах во внешних силовых полях. Эти методы далее применяются при расчетах пространственно-временных и спектральных характеристик стационарных и нестационарных волновых полей в открытых и замкнутых системах.

Исходную систему связанных линейных стационарных волновых уравнений при частоте и или энергии Е = Ни можно представить в виде неавтономной системы дифференциальных уравнений первого порядка

^Ъ{х) = £)(ы,х)Ъ(х), (1)

где Ф(ж) = (^х, гр2, ■ ■ ■ , ~ вектор-столбец состояния, а Е)(и>,х) - динамическая матрица размерности N х N. Решение Ф(ж) в точке х связано с решением Ф(а) в точке а соотношением

Щх) = М(х,а)Ща), (2)

где М(х, а) - матрица перехода (трансфер-матрица) из а в х.

Расчетная область разбивается на сегменты [¡£¿,£¿+1], внутри которых элементы динамической матрицы Е)(ш,х) непрерывны. Вычисление трансфер-матрицы на сегментах производится либо на основе фундаментальной системы решений (ФСР) волновых уравнений, входящих в (1), либо на основе конечно-разностной (КР) схемы интегрирования системы уравнений (1). В зависимости от применяемой квадратурной схемы используются одношаговые и многошаговые трансфер-матрицы, обеспечивающие необходимую численную устойчивость. Матрица перехода из а = х^ в х = = X] находится последовательным перемножением трансфер-матриц через сегменты

М(х, а) = М^-хМ^-г . . . М;+2,г+1^+1,г- (3)

Для задач, рассмотренных в главе 3, методы ФСР и КР, применяемые совместно, образуют основу гибридного алгоритма расчета волновых полей.

В третьей главе рассмотрены задачи об эволюции квантовых волновых пакетов в многослойных структурах, с учетом влияния внешних си-

ловых полей, при условии, что длина фазовой когерентности волн де Бройля больше толщины слоев.

В разделе 3.1 рассмотрена задача одномерного рассеяния электронного волнового пакета, падающего слева на периодическую решетку конечной длины, содержащую N одинаковых ячеек с периодом элементарной трансляции (1. Данная модель описывает композиционные сверхрешетки, являющиеся важными элементами приборов фотоники, опто- и наноэлектроники.

Эволюцию волнового пакета описывает решение нестационарного уравнения Шредингера в форме интеграла

Ф(ж, 0 =

с(Е)ф(Е,х)ехр{-іЕі/ГЇ)(ІЕ, (4)

где ф(Е,х) - решения стационарного уравнения Шредингера, Е - энергия стационарных состояний, с(Е) - спектральная функция

с(Е) =

■ф\Е,х)Ъ(х,0)<1х, (5)

где 0) - волновая функция пакета в момент времени í = 0. При помощи метода трансфер-матрицы получены базисные решения ф(Е,х) уравнения Шредингера, описывающие стационарное рассеяние падающей слева монохроматической волны. Парциальные амплитуды падающей А0, отраженной В0 и прошедшей 1 волн связаны через эффективную трансфер-матрицу Ме/. Амплитуды и коэффициенты отражения г(Е) = В0/А0, И(Е) и прохождения Ь(Е) = Лдг+х/Ло, Т{Е) выражены через элементы эффективной трансфер-матрицы М^. Найдены нули амплитуды отражения г(Е) и полюсы амплитуд рассеяния г(Е) и Ь{Е) в плоскости комплексной переменной Е вблизи порогов и в середине полос (зон) пропускания. Показано, что с увеличением N полюсы приближаются к вещественной оси энергий. Мнимая часть полюсов пороговых резонансных состояний много меньше мнимой части полюсов резонансных состояний середины полосы, соответственно, пороговые квазистационарные состояния являются более долгоживущими.

Рис. 1. Полоса пропускания Т(к), к = Ы1\/2тЕ, полюсы кл, нули к°я амплитуды отражения, А^ - седловая точка, АВСБ - перевальный контур

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

в (4) от полюсов амплитуд рассеяния и фазовой седловой точки, объясняющие наличие порогов и пространственно-временных осцилляций ¿). При заданных длине и других параметрах решетки по положению полюсов проведены аналитические и численные оценки энергии Ед, время жизни Тд квазистационарных состояний вблизи порогов и в середине й-й полосы пропускания. Получены также аналитическое выражение и оценки групповой скорости ьд волнового пакета в решетке конечной длины. В частности, для решетки ¿-барьеров в приближении сильной связи (ш/йя-) » 1 результаты представлены в таблице для двух случаев: (1) = вблизи потолка

или Рх = (И-р) «Ж вблизи дна полос и (2) V « N/2, « - Н/2\ с N в середине полос: {и - номер квазистационарного состояния, отсчитанный от потолка полосы)

Таблица 1. Оценки Ед, тд, уд и

Ец Тя Д*

1 Ч-жН-ч-Нф)) тПЧ^3 втг5/^2 2тг3Пз2и1 тШ2 N тШ3И2 2тг Ч2з2иг

2 Ч'-шН-Ч'т?)) шП'^АГ 8тг 2тт2Пз2 т£1 А2 тШлЫ 2тг Ч'з1

Здесь Еео = (Я27г2й2) / (2тсР) -пороговая энергия, Г2 = Ьг22тЩЬ - мощность ¿-барьера, Щ - высота, а Ь- толщина моделируемого потенциального барьера, с? - период решетки.

Проведено численное моделирование динамики процесса рассеяния на решетке дельта-барьеров гауссовых волновых пакетов с разными спектральными ширинами и положениями относительно полосы пропускания решетки. Если спектрально узкий пакет возбуждает пороговое долгожи-вущее квазистационарное состояние, то распад возбуждения в решетке проходит по экспоненциальному временному закону, причем время жизни, оцененное по мнимой части полюсов амплитуды рассеяния с погрешностью менее 1% близко к значению 7д, полученному из расчета вероятности Р(Ь) нахождения электрона в решетке. При возбуждении пакетом двух или более пороговых квазистационарных состояний в решетке формируется сложная структура возбуждения с несколькими максимумами, между которыми происходит перекачка интенсивности, причем при достижении границ происходит периодическая эмиссия вторичных пакетов из решетки, частичное внутреннее отражение и осцилляции распадного тока. Построены пространственно-временные диаграммы движения максимумов пакетов вне и внутри решетки и по ним определены групповые скорости пакетов. Величины скоростей движения максимумов внутри решетки и вне

решетки с хорошей степенью точности (в пределах нескольких процентов) совпадают с величинами, рассчитанными по соответствующим аналитическим формулам.

В разделе 3.2 рассмотрена задача о распаде возбуждения, образованного суперпозицией квазистационарных состояний. Закон распада определяется наиболее долгоживущими квазистационарными состояниями. Изучен распад возбужденного состояния в открытой двухбарьерной структуре. Расчет полюсов амплитуды пропускания дает значение времени жизни возбужденного состояния с погрешностью менее 2%, совпадающее со значением, определенным из зависимости Р(£). Показано, что при возбуждении двух квазистационарных состояний распадный ток осциллирует с частотой,

определяемой разностью энергий этих состояний.

В разделе 3.3 изучено прохождение стационарных электронных волн и волновых пакетов через двойные зеркально-симметричные сверхрешетки. Подробно изучены две модели: в одной потенциальная энергия электрона аппроксимирована кусочно-постоянной функцией (с параметрами сверхрешетки Si/CaF?), в другой - равноотстоящими ¿-функциями. Проведен анализ перестройки полос пропускания двойной сверхрешетки при изменении параметров, не нарушающем её зеркальную симметричность. Показано, что в каждой полосе имеется пара смежных резонансных пиков единичной прозрачности, которые могут сближаться, сливаться, а величина результирующего пика может уменьшаться вплоть до исчезновения (коллапс резонансов прозрачности). Показано, что увеличение отношения а мощности барьера второй внутренней подрешетки к мощности барьера первой подрешетки приводит к смещению пиков прозрачности, положения которых определяются нулями амплитуды отражения г. Для двойной симметричной сверхрешетки дельта-барьеров с произвольным N выведено уравнение, определяющее энергии резонансов и обобщающее результат, полученный для трехбарьерной структуры1. При рассеянии волнового пакета на симметричной двойной сверхрешетке в условиях коллапса время жизни разных квазистационарных состояний в полосах пропускания имеет разную зависимость от а. Если при прохождении волнового пакета возбуждаются квазистационарные состояния, соответствующие неколлапсирующим резонансам, то с увеличением а время жизни этих состояний увеличивается. Если же возбуждаются квазистационарные состояния в области коллапсирующих пиков, то с увеличением а время жизни таких состояний уменьшается.

В разделе 3.4 рассмотрена задача о рассеянии волнового пакета периодической решеткой конечной длины, находящейся в электрическом поле,

1 Горбацевич A.A., Журавлев М.Н., Капаев В.В. Коллапс резонансов в полупроводниковых гетеро-структурах как переход с нарушением симметрии в открытой квантовой системе // ЖЭТФ. 2008. Т.134. N° 2(8). С.338-353.

при учете возбуждения квазистационарных состояний в ней. Предварительно решена задача Штурма-Лиувилля при нулевых граничных условиях и проведен анализ перестройки стационарных состояний электрона в одномерной решетке потенциальных барьеров во внешнем однородном поле У(х) — Рх. При кусочно-постоянной аппроксимации У(х) — Fd[x/d] и применении приближения сильной связи исследованы решетки для (И » 1), вычислена энергетическая плотность состояний, построены качественные диаграммы локальной устойчивости решений, получены аналитические оценки границ пространственной локализации электрона в решетке для произвольных стационарных состояний, изучена перестройка состояний при изменении внешнего электрического поля. Применение гибридного алгоритма позволило исследовать периодическую решетку в электрическом поле без использования приближения сильной связи и кусочно-постоянной аппроксимации потенциальной энергии электрона. Структура энергетической плотности состояний электронов в решетке д{Е)=йп/йЕ сильно изменяется в зависимости от приложенного к решетке напряжения V: при \q\f \ < 2АЕв0, (где 9 - элементарный заряд, АЕз0 - полуширина й-й энергетической зоны в отсутствие внешнего поля) состояния можно подразделить на локализованные (краевые) и делокализованные (зонные); при \цУ\ > 2АЕ8о все состояния локализованные, но вблизи порогов полосы остаются краевые состояния, а в середине полосы образуются лестничные Ванье-Штарковские состояния, плотность состояний которых не зависит от энергии. Разработанный алгоритм применен для расчета тока через короткопериодную сверхрешетку 5г/Са^ в электрическом поле. Решетка является интерференционным фильтром и перераспределяет энергию распространяющегося волнового возбуждения, усиливая одни спектральные компоненты и подавляя другие в прошедших, отраженных и захваченных решеткой импульсах. Последний тип импульсов может существовать достаточно долго после ухода основных отраженного и прошедшего импульсов. Исследованы полюсы амплитуды рассеяния в комплексной энергетической плоскости и рассеяние волновых пакетов разного спектрального состава и

положения в зависимости от приложенного внешнего поля. Показано, что при большом значении внешнего поля, такого, что \qV\ > 2AEso туннельная прозрачность решетки резко уменьшается, и пакет практически полностью отражается от решетки, несмотря на сильное тянущее поле. Поведение импульса внутри решетки зависит от положения спектрального максимума Ес пакета относительно полосы пропускания, в частности, пакет, спектрально расположенный вблизи верхнего порога полосы, проникает достаточно глубоко в решетку (до классической точки поворота), образуя внутри решетки долгоживущее возбуждение. Время жизни этого возбуждения определяется соответствующими полюсами, близкими к вещественной оси в сильном поле. Отметим, что подавляется когерентный механизм переноса пакета, при наличии неупругого рассеяния в решетке будет происходить последовательный перенос возбуждения по спектральным ступенькам лестницы Ванье-Штарка.

Четвёртая глава посвящена изучению электронных переходов в фотоэлектрохимической ячейке (ФЭХЯ) нового поколения солнечных элементов на органических красителях, предложенной М. Гретцелем2. Исследование процессов переноса заряда в таких устройствах позволяет указать пути повышения эффективности солнечных элементов.

В разделе 4.1 дано описание ФЭХЯ, принцип действия которого заключается в разделении заряда на границе между полупроводниковой на-ночастицей и молекулой сенсибилизатора (хромофором). Описаны основные процессы электронного переноса после оптического возбуждения хромофора: обратный излучательный переход электрона (флуоресценция), ин-жекция фотовозбужденного электрона в наночастицу, перенос электрона через слой наночастиц в коллектор.

В разделе 4.2 на основе подхода3 рассмотрено температурное ушире-

2 O'Regan В., Gratze! М. A low-cost, high-efficiency solar cell based on dye-sensitized colloidal ТЮ2 films // Nature. 1991. V.353. N.6346. P.737-739.

3 Осадько И.С. Теория формы однородной оптической полосы примесного центра в аморфных и кристаллических средах // ЖЭТФ. 1986. Т.90. №4. С.1453-1465. Осадько И.С. Селективная спектроскопия одиночных молекул. М.: Физматлит, 2000. 320 с.

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

В разделе 4.3 предложена модель бездиссипативного прохождения электронного волнового пакета от хромофора через систему потенциальных барьеров, разделяющих наночастицы, до коллектора. Для объяснения экспериментально наблюдаемых двух стадий распада электронного возбуждения - сначала быстрого (десятки фемтосекунд), а затем более медленного (десятки, сотни пикосекунд) - предложенная модель позволяет выделить основные факторы, влияющие на перенос электрона и оценить значения временных характеристик этого процесса. Изучены вероятности нахождения электрона в хромофоре (Р\(1)) и в системе «хромофор-слой наноча-стиц» {Р1у(£)). Вероятность характеризуется начальным быстрым распадом, переходящим затем в медленное убывание, которое сопровождается колебаниями, обусловленными интерференцией спектральных компонент волнового пакета при их отражении от потенциальных барьеров, ограничивающих наночастицы (обратный электронный перенос). Вероятность РдК^) имеет начальный горизонтальный участок, соответствующий прохождению волновым пакетом слоя наночастиц, который затем сменяется относительно быстрым спадающим участком, отражающим уход основного тела пакета в коллектор, и относительно медленным участком, формируемым долгоживу-щими и распространенными на всю цепочку наночастиц квазистационарными состояниями. Показано, что большое влияние на динамику волнового пакета оказывают размеры и разброс размеров наночастиц, а также их количество в транспортной цепочке от молекулы фотосенсибилизатора до коллектора. Для системы одинаковых по размеру наночастиц существует оптимальный размер, при котором реализуется наиболее быстрый уход

электрона в коллектор.

В пятой главе рассмотрены задачи расчета нестационарного тока из квантово-размерных многослойных гетероструктур.

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

В разделе 5.2 проводится дальнейшее развитие идеи получения гармонической составляющей нестационарного тока из гетероструктуры. Построена теория нестационарной фотоэмиссии, основанная на применении метода матрицы плотности, используемого при расчетах быстро протекающих квантовых кинетических процессов, обобщающая трехшаговую теорию фотоэмиссии4. Элементы матрицы плотности находятся из решения квантового кинетического уравнения. Рассчитаны недиагональные элементы матрицы плотности, зависящие от времени и учитывающие различные неупругие процессы, ответственные за электронную релаксацию. Плотность тока j(r, t) фотоэмиссии выражается через элементы матрицы плотности pf^p,(t) во втором порядке теории возмущений по электромагнитному полю и через матричные элементы плотности стационарного тока jp^r) как

j(r,i) = 2^p®(i)jp>(r). (6)

р,р'

Фотоэмиссионный ток равен сумме (6) по стационарным состояниям р и р' непрерывного спектра задачи рассеяния, принадлежащим регистрируемой

4 Solid-State Photoemission and Related Methods. Theory and Experiment. / Eds. Schattke W., Van Hove M.A. Wiley-VCH. 2003. 515 p.

1.5

1.6

1.7

1.8

1.9 Кс(£"„), эВ

в)

Рис. 2. а) Схема фотоэмиттера с поверхностной гетероструктурой, б) спектр коэффициента прохождения Тр через гетероструктуру; в) надвакуумные низкоэнергетические дублеты полюсов амплитуд рассеяния в комплексной энергетической плоскости Ец 1)Я2 и Ецз^щ

полосе энергии фотоэлектронов выше вакуумного уровня (рис.2,а). Элементы матрицы плотности рнесут информацию о фотовозбуждении электронов из глубоко лежащих стационарных состояний, а также о процессах неупругого рассеяния. Численный расчет тока фотоэмиссии проведен для одномерной модели (рис.2,а). Амплитуда прохождения через гетероструктуру и матричные элементы тока ^р>р имеют полюсные особенности в нижней полуплоскости комплексной энергии электрона при значениях энергии Ея = Т1еЕ{1 + йтЕц. С этими полюсами связано положение узких пиков коэффициента прохождения через гетероструктуру Тр (рис.2,б,в). Проведено численное моделирование тока фотоэмиссии при импульсном освещении светом длительностью ¿о- Показано, что при этом могут возбуждаться надвакуумные близколежащие квазистационарные состояния, характеристики которых определяются полюсами амплитуд рассения. Времена жизни квазистационарных состояний дублетов тт « К/\\тЕт\, тЕ2 и Н/\1тЕН2\. Разность энергий дублетов Д-Ед12 = йеЕдг — Ке£д1 эВ определяет разностную частоту 1^12 = АЕЯ12/2тт}1 и период Т12 = 1/^12 осцилляций фототока. Осцилляции тока будут эффективно наблюдаться при выполнении неравенств тльТн2,тр 3> Т12, где тр = Н/чР - время релаксации электронов, определяемое неупругим рассеянием. При 0 < < < ¿о значение тока возрастает и осциллирует. Определенные из (6) частота и период осцилляций близки к

значениям Щ2 и Тг2 с погрешностью менее 4% (рис.3,а). При увеличении длительности импульса накачки ^ ток асимптотически стремится к стационарному значению объемного фототока. При й > ¿о» после прекращения

бГ

2 4 6 8 10 Г, пс 12 14 16 18 пс

Рис. 3. Производная плотности тока по времени а) после включения импульса накачки и б) после выключения импульса накачки

действия светового импульса Е{£)= 0, фотоэмиссионный ток определяется смешиванием и распадом заселившихся при накачке квазистационарных состояний в гетероструктуре. Он убывает до нуля, испытывая осцилляции почти той же разностной частоты, что и при накачке (рис.3,6). Показано, что для генерации разностной гармонической составляющей переменного фототока длительность импульсов накачки должна быть порядка времени релаксации, а интервалы между ними - порядка времени жизни квазистационарных состояний. Характеристики тока определяются параметрами гетероструктуры. Для толщин слоев 1 — 10 нм и высоты барьеров 0.5 — — 2.5 эВ можно обеспечить времена жизни квазистационарных состояний 10~2 —(3 х 102) пс и генерируемые разностные частоты для них 1011 —1014 Гц. Варьируя параметры поверхностной гетероструктуры, можно изменять времена жизни тд1,тд2 и разностные частоты дублетных квазистационарных состояний, что меняет форму кривых зависимости фототока от времени, анализ этих кривых может дать информацию о значениях времен релаксации тр = возбужденных электронов в фотокатоде.

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

Ранее теоретически и экспериментально изучалось распростанение маг-

нитоупругих волн в неограниченных объемных5 и планарных средах 6, в сверхрешетках7. В разделах 6.1-6.3 излагается развитая нами теория распространения упругих волн в структурах, состоящих из чередующихся магнитных (М) и немагнитных (D) слоев конечной длины в области MAP, получены спектры прохождения и отражения ультразвука, изучено рассеяние ультразвукового волнового пакета М-слоем.

В разделе 6.1 представлены результаты исследования дифракции монохроматической с частотой uj правополяризованной ультразвуковой волны, нормально падающей на ферромагнитный слой в области магнитоакустиче-ского резонанса (MAP). Ферромагнетик намагничен до насыщения перпендикулярно плоскости слоя. Связанные магнитоупругие волны описываются системой уравнений, состоящей из уравнения Ландау-Лифшица для намагниченности m и волнового уравнения для упругой волны и в линейном приближении. Из решения этой системы получены четыре магнитоупру-ГИХ ветви &1,2,3,4- В отсутствие МЭГНИТОупруГОГО взаимодействия (£?2 = 0) пересечение спинволновой кю и упруговолновой &2о ветвей определяет координату ко, шо нижней точки MAP.

В М-слое магнитоупругие волны представлены суперпозицией четырех плосковолновых компонент, трансфер-матрица через слой четырехмерна, в немагнитных средах распространяется только двухкомпонентная упругая волна. Для построения редуцированной двухмерной матрицы перехода упругой волны через Af-слой учтены условия закрепления спинов на границах раздела слоев, позволяющие получить необходимые соотношения между парциальными амплитудами магнитоупругих волн. Получены аналитические выражения для амплитуд рассеяния г{ш) и t(uj) упругой волны при различных условиях закрепления спинов на границах раздела

5 Ахиезер А.И., Барьяхтар В.Г., Пелетминский C.B. Спиновые волны. М.: Наука, 1967. Гуревич А.Г., Мелков Г.А. Магнитные колебания и волны. М.: Физматлит, 1994.

6 Гуляев Ю.В., Зильберман П.Е. Магнитоупругие волны в пластинах и пленках ферромагнетиков // Известия вузов. Физика. 1988. Т. 31. С. 6-23.

7 Ignatchenko V.A., Laletin O.N. Magnetoelastic ground state and waves in ferromagnet-nonmagnetic dielectric multilayer structures // Phys. Rev. B. 2007. V.76. N.10. P. 104419

слоев. Ширина области MAP характеризуется параметром относительной когерентной связи I = 1{ш) между упругими и магнитными колебаниями, который в резонансе по величине близок к единице и быстро спадает вне области резонанса. Показано, что при прохождениии ультразвуковой волны через М-слой в области MAP возбуждается локализованная спиновая волна, при этом реализуются условия антирезонанса Фано8. В спектре Т этот эффект проявляется в наличии узких линий полного отражения соседствующих с более широкими максимумами (рис. 4, г). Аналогичная структура спектра может возникать при отражения электронной волны в электронном волноводе с примесью конечных размеров9. В модельных расчетах использованы физические параметры железо-иттриевого граната (ЖИГ). Исследованы случаи тонкого и толстого слоя толщины а по сравнению с длиной ультразвуковой волны Л. В тонком слое а ~ А узкие линии отражения в области MAP становятся более широкими по отношению к таким же линиям вне резонанса. Учет затухания для тонких слоев ЖИГ практически не влияет на коэффициенты отражения и пропускания. В толстом слое а»А в области MAP, как правило, формируется полоса сильного отражения. Ширина и количество линий отражения определяются величиной магнитоупругой связи, толщиной М-слоя, параметрами затухания волн в структуре и расстоянием по частоте от точки MAP. Интерференция волн в многослойных структурах приводит к эффективному усилению влияния магнитоупругого взаимодействия на рассеяние упругих волн, как ранее было отмечено в работе10.

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

8 Miroshnichenko А.Е., Flach S., Kivchar Y.S. Fano resonances in nanoscale structure // Rev. Mod. Phys. 2010. V. 82. N.3 P.2257-2298.

9 Ким Ч.С., Сатанин A.M., Джо Ю.С., Косби P.M. Коллапс резонансов в квазиодномерных квантовых каналах // ЖЭТФ. 1999. Т.116. №1. С.263-275.

10 Игнатченко В.А., Кузьмин Е.В. Спектр связанных магнитоупругих колебаний в тонкой магнитной пленке // ЖЭТФ. 1964. Т.47. №11. С.1814-1821.

МАР в случае, когда М-слои намагничены до насыщения вдоль направления оси сверхрешетки. Получены рекуррентные формулы расчета трансфер-матрицы М(хм+1,х\) через структуру, амплитуды прохождения и отражения г (и/) , а также выражения для связанных волновых полей в М-слоях и упругих полей в £>-слоях структуры при различных условиях закрепления спинов на границах раздела слоев. Показано, что при прохождении ультразвука с частотой, близкой к частоте МАР, усиливаются резонансные эффекты, определяемые магнитоупругой связью в М-слоях и рассеянием волн на границах раздела магнитных и немагнитных сред. Увеличение толщины каждого М-слоя приводит к увеличению густоты и смещению влево по частоте узких линий полного отражения в спектре коэффициента пропускания ультразвука в любом частотном интервале. Показано, что затухание приводит к уменьшению ультразвуковой прозрачности структуры, тогда как на отражение ультразвука затухание практически не оказывает влияние. Изучено влияние величины внешнего магнитного поля на коэффициент прохождения ультразвуковой волны Т в области МАР при фиксированных остальных параметрах многослойной магнитной сверхрешетки. Отмечено, что согласованный подбор числа и толщины М- и .О-слоев в многослойных структурах может позволить создать устройства магнитной электроники, пропускающие или отражающие ультразвуковые волны в заданных спектральных полосах вблизи МАР, положение которых зависит от значения внешнего магнитного поля.

В разделе 6.3 рассматривается задача о рассеянии упругого правопо-ляризованного волнового импульса, нормально падающего на ферромагнитный слой, намагниченный перпендикулярно его плоскости. Волновое поле ультразвукового импульса в момент времени Ь можно представить в виде

и(х, t) =

c(Lü)e-mtu(io,x)du>, (7)

где c{üj) - спектральная функция, центрированная на частоте Í20 со спектральной шириной равной v/a, v - скорость максимума пакета вне М-слоя, а - пространственная ширина пакета. Вычисления стационарных волно-

х 10

l4l

ваш»

x 10"4 ^

¡шшш ¡Щ111

co"/coH -5 B)

1.0365

1.0370 1.0375 C0/00H

Рис. 4. Поверхность |í| в области MAP в нижней комплексной полуплоскости частот для а) тонкого М-слоя, б) слоя промежуточной толщины (а > А), в) толстого М-слоя, г) фрагмент спектра пропускания Т(и) ультразвуковой волны через тонкий М-слой

вых полей и(со,х) и m(uj,x) во всей пространственной области проведены с помощью метода трансфер-матрицы (раздел 6.1). Рассмотрены волновые пакеты со спектральным максимумом вблизи нижней точки МАР. При £>2 = 0 амплитуда t(uj) имеет полюсы, которые являются решениями уравнения т~1 = 0, где r~1=cosfc„oa — §(^noA+ k/kno)smknoa, п= 1,2 и лежат в нижней полуплоскости комплексной частоты При Lüja = Lú'ja + ш"а. При увеличении толщины слоя а эти полюсы, сгущаясь, смещаются влево и к действительной оси. Полюсы ujja лежат точно под нулями функции sinA;na на акустической ветви, которые определяют частоты полного резонансного

пропускания упругих волн. Уравнение sinfcna = 0 на спинволновой ветви (ниже точки MAP c¿o ей соответствует къ, а выше к\) есть условие образования стоячих спиновых волн (полуволновой слой). При £?2 ф 0 у амплитуд А„, г(ш) и t(uS) имеется два типа полюсов шrv, совпадающих с нулями знаменателя. Полюсы первого типа происходят от параметров г"1 = 0 и связаны с эволюцией полюсов Wja квазиакустической ветви (квазиупругие резонансы). Полюсы второго типа появляются за счет их отделения от нулей sinfc„a на спинволновой ветви и схода в нижнюю полуплоскость (квазиспинволновые резонансы), при удалении от точки бифуркации нули t{tü) продолжают оставаться на реальной оси, а отделившиеся полюсы с ростом £?2 отходят вниз. Показано, что спектральный вклад полюса ujrv пропорционален значению групповой скорости на соответствующей ветви. В области нижнего MAP групповая скорость упругих волн ~ vt много больше групповой скорости спиновых волн ~ 2r/LJo/vt, поэтому квазиакустические полюсы значительно мощнее квазиспинволновых полюсов. При изменении толщины слоя а и напряженности магнитного поля Я полюсы обоих типов движутся в гиперпространстве {uj',uj",a, Я} по своим траекториям. С ростом Я область MAP перемещается к более высоким частотам и коротким длинам волн. Показано, что изменяя толщину а, напряженность поля Я, а также параметры пакета fio и ст, можно наблюдать возбуждение и эволюцию вторичных волновых структур вне и внутри пластины, связанных с возбуждением квазистационарных состояний в пластине. В зависимости от толщины М-слоя определены три режима.

а) Тонкий М-слой толщины а~2л/к0. Проведен расчет квазиспинволновых полюсов в области MAP (рис. 4,а). При рассеянии на пластине происходит отражение и прохождение импульса, при этом возможно возбуждение квазистационарного состояния, распадающееся по экспоненциальному закону. Время жизни этого состояния с погрешностью менее 3% совпадает с величиной, определенной по мнимой части полюса tr = 2-k/lo" .

б) Слой промежуточной толщины a ~ 27r/(fci(wo) —^г(^о)). В этом случае в области MAP квазиакустические полюсы, локализованные вдоль

-0.4 -0.2 0 0.2 -0.4 -0.2 0 0.2

X, СМ X, СМ

Рис. 5. Процесс рассеяния волнового пакета толстым М-слоем (П0/<^н = 1 0075, а = 0.075 см, х0 = -0.25 см). Границы пластины обозначены вертикальными пунктирами.

дуги (рис. 4,6), приближаются к реальной оси настолько, что начинают определять рассеивающие свойства системы и структуру квазистационарных состояний. Квазиспинволновые полюсы расположены как около квазиакустических полюсов, так и недалеко от реальной оси вблизи нулей Исследовано поведение полюсов при изменении магнитного поля и рассеяние волнового пакета при его разном спектральном положении относительно центра МАР. Показано, что при возбуждении двух спектрально близких резонансов могут возникнуть осцилляции, приводящие к пульсациям ультразвукового сигнала на разностной частоте и'т —и>'т.

в) Толстый М-слой а » 2ж/ (^(шо) — По сравнению с преды-

дущим случаем квазиакустические полюсы еще сильнее сгущаются, об-

разуя дугообразный хребет, приближающийся к реальной оси в области MAP (рис. 4,в). В этой области групповые скорости на обеих ветвях меньше скорости упругих волн v, поэтому квазиакустические полюсы уменьшаются по мощности и смешиваются с квазиспинволновыми полюсами, нули t(uj) становятся реже, а участки вокруг них шире, что может привести к уширению линий полного отражения. При таком сгущении полюсов спектрально узким пакетом перекрывается много полюсов и нулей t(uj). Изучено прохождение волнового пакета (7) через М-слой в области MAP. Внутри слоя пакет при каждом отражении от границы раздела разбивается на два импульса, движущихся со скоростями, близкими к групповым скоростям Vgl и vg2 на двух спектральных ветвях при частоте спектрального центра пакета, которые значительно меньше фазовой скорости упругой волны в М-слое vt и вне М-слоя V. Показано, что скорости максимумов хорошо разрешенных импульсов внутри слоя практически равны указанным групповым скоростям (вычислительная погрешность меньше 1%). Расчет относительной интенсивности упругого поля внутри М-слоя показывает, что при внутренних отражениях может существенно увеличиться энергия упругих колебаний в пленке. Увеличение происходит за счет перекачки энергии из магнитной подсистемы в упругую подсистему.

В разделе 6.4 на основе метода четырехмерных трансфер-матриц изучаются пространственно-неоднородные сверхпроводящие системы. Волновые поля и и v элементарных возбуждений в сверхпроводящем слое описываются системой уравнений Боголюбова - де Жена11. Для однополостной поверхности Ферми радиуса hkFX, где X = N,S\,S2, обладающей центром инверсии, характеристическое уравнение системы уравнений Боголюбова -де Жена имеет четыре корня ¿1,2,3,4. а общее решение для и и v полей есть суперпозиция четырех компонент вида Ajexp(ikjx). Вектор-столбец состояния Ф(х) = (u,u',v,v')T в произвольной точке х связан с решением Ф(а) в некоторой точке а соотношением (2). Для стационарной задачи андреев-

11 Свидзинский A.B. Пространственно-неоднородные задачи теории сверхпроводимости. М.: Наука, 1982. 312 с.

ского отражения в несимметричном 2 переходе, с учетом различных поперечных импульсов Ферми в слоях ЛГ, 5г, построена 4-мерная трансфер-матрица М(а, 0) перехода через Л^слой, толщиной й, получено спектральное уравнение для энергии квазичастиц е = е(Дь Дг, &Р52)> где ДьДг - потенциалы спаривания в 5х и 5г слоях. Показано, что в андреевском приближении крх = крБ! = крБ2 спектральное уравнение для е переходит в уравнение, полученное в работе12, а для симметричного перехода в уравнение, полученное в работе13. Проведен численный расчет е в зависимости от величины параметра асимметрии импульсов Ферми 5кр = = {кры/крз + йгвАглг) /2 — 1. В зависимости от величины параметра асимметрии 5кр и фазы изучены асимптотики спектра е и показано, что с увеличением расстройки дкр ослабевает фазовая зависимость е.

Основные результаты и выводы диссертации

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

1. Для волновых уравнений с переменными коэффициентами построен рекуррентный алгоритм вычисления одношаговых трансфер-матриц. Изучены спектральные свойства таких матриц, установлена их численная устойчивость, определены локальные и глобальные погрешности расчета трансфер-матрицы. Построенный алгоритм является частью гибридного алгоритма расчета волновых полей, энергетических спектров, полос пропускания и других характеристик дифракции волн в многослойных системах.

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

12 Куплевахский С.В., Фалько И.И. Андреевские состояния индуцированные током в сверхпроводящих контактах // ФНТ. 1991. Т.17. N.8. С.961-970.

13 Кулик И.О.Пространственное квантование и эффект близости в Б-Г^-Б-контактах // ЖЭТФ 1969. Т.57. №.11. С.1745-1758.

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

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

4. Для периодической решетки конечной длины, помещенной в электрическое поле, проведен численно-аналитический расчет волновых функций стационарных состояний электрона при различных значениях энергии электрона и напряженности внешнего поля. Построены диаграммы локальной устойчивости волновых функций этих состояний. Изучена перестройка однозонной плотности состояний в зависимости от величины внешнего электрического поля. Показано, что при усилении поля увеличивается число краевых локализованных состояний зоны, а при превышении определенного критического значения напряженности поля в середине зоны формируется горизонтальное плато плотности состояний, соответствующее образованию спектральной лестницы Ванье-Штарка в этой области энергий.

5. Исследованы два основных электронных процесса в фотоэлектрохимической ячейке Гретцеля - излучательный переход из возбужденного

электронного состояния в основное состояние и инжекция электрона в систему полупроводниковых наночастиц. Получены асимптотические выражения для полуширины бесфононной линии флуоресценции при учете взаимодействия с низкоэнергетическими возбуждениями в туннельных двухуровневых системах аморфной среды в низкотемпературной области. Обнаружено существенное влияние упорядоченности нанослоя и дисперсии размеров наночастиц на скорость процесса переноса. Показано, что после оптического возбуждения вероятность Р{Ь) нахождения электрона в молекуле красителя сначала быстро уменьшается в диапазоне десятков фемтосекунд, а затем медленно убывает в диапазоне десятков и сотен пи-косекунд, последний участок обусловлен распадом долгоживущих квазистационарных состояний в нанослое и эмиссией в коллектор.

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

7. Изучено поведение коэффициента прохождения монохроматической правополяризованной ультразвуковой волны через один М-слой и через сверхрешетку «ферромагнетик/немагнетик» конечной длины в области МАР. Показано, что при прохождении ультразвука с частотой, близкой к частоте МАР М-слоев, в сверхрешетке усиливаются резонансные эффекты, определяемые магнитоупругой связью в магнитных слоях и рассеянием волн на границах раздела магнитных и немагнитных сред. Увеличение толщины каждого М-слоя приводит к увеличению густоты узких линий полного отражения в спектре коэффициента пропускания ультразвука. Кроме того, вблизи МАР возникают антирезонансы Фано, вызываемые взаимодействи-

ем рассеиваемой упругой волны со спиновыми волнами, локализованными в магнитных слоях. Затухание волн приводит к уменьшению ультразвуковой прозрачности структуры, а на отражение ультразвука оказывает меньшее влияние. Согласованный подбор числа и толщины магнитных и немагнитных слоев в многослойных структурах может позволить создать устройства, пропускающие или отражающие ультразвуковые волны в заданных спектральных полосах вблизи MAP, положение которых зависит от значения внешнего магнитного поля. Проведен анализ резонансного взаимодействия правополяризованных ультразвуковых импульсов со спиновыми возбуждениями ферромагнитной пластины в области MAP. Выделены два типа полюсов. Резонансные свойства тонких пластин определяются, в основном, полюсами спинволнового происхождения, а свойства толстых пластин смешиванием акустических и спинволновых полюсов. Проведенные численные расчеты распространения ультразвуковых правополяризованных импульсов через ферромагнитные пластины показали, что в пластинах при прохождении пакетов могут возбуждаться квазистационарные состояния и вторичные импульсы, обусловленные интерференцией и перекачкой энергии между упругой и спиновой подсистемами.

8. Проведен расчет спектра локализованных андреевских состояний в несимметричном S1NS2 переходе с разными поперечными импульсами Ферми в слоях. При увеличении разности поперечных импульсов Ферми в слоях, характеризуемых параметром асимметрии ÖIcf, ослабевает фазовая зависимость энергии стационарных состояний элементарных возбуждений и при SkF 1 эта зависимость исчезает.

Публикации по теме диссертации

1. Osad'ko I.S., Shtygashev A.A. The tunnelon theory of optical dephasing and impurity homogeneous optical band in amorphous solids: finding the two-level system parameters // J.of Luminescence 1987. V.36. P.315-325.

2. Osad'ko I.S., Shtygashev A.A. Theoretical aspects of optical dephasing in polymers and glasses. // J.of Luminescence 1987. V.37. P.255-267.

3. Осадько И.С., Штыгашев A.A. Извлечение информации о туннельных системах стекол и полимеров из данных по низкотемпературному уширению оптических линий // ФТТ. 1987. Т.29, №5. С.1550-1554.

4. Осадько И.С., Штыгашев A.A. Влияние туннельных степений свободы полимеров и стекол на форму однородной оптической полосы // Оптика и спектроскопия. 1989. Т.66. №1. С.86-91.

5. Пейсахович Ю.Г., Штыгашев A.A. Спектр андреевских состояний в несимметричном SNS-переходе // Физика низких температур. 1999. Т.25. №5. С.455-458.

6. Пейсахович Ю.Г., Штыгашев A.A. Полное отражение ультразвука от ферромагнитной пластины при закреплении спинов на поверхности // ЖЭТФ. 2000. Т.118. №1(7). С.213-222.

7. Shtygashev A.A., Ovchinnikov Yu.E., Shklover V. Simple quantum models of electron injection from a sensitizer molecule to semiconductor nanocrysta // Solar Energy Materials and Solar Cells. 2003. V.76. P.75-84.

8. Величко A.A., Илюшин В.А., Филимонова Н.И., Штыгашев A.A. Вольт-амперная характеристика сверхрешеток на основе структуры CaF^/Si // «Актуальные проблемы электронного приборостроения» / Материалы VII международной конференции. 2004. N.2. С. 16-19.

9. Пейсахович Ю.Г., Штыгашев A.A. Распад низкоэнергетического квазистационарного состояния в квантовой яме // «Актуальные проблемы электронного приборостроения» / Материалы VIII международной конференции. Новосибирск, 2006. Т.2. С.47-52

10. Peisakhovich Y.G., Shtygashev A.A. The Gauss wave packets scattering on the regular lattice of finite length // «Актуальные проблемы электронного приборостроения» / Материалы VIII международной конференции. Новосибирск, 2006 Т.1. С.19-21.

11. Пейсахович Ю.Г., Штыгашев A.A. Одномерная квантовая механика. Новосибирск: Издательство НГТУ, 2007. 476с.

12. Peisakhovich Yu.G., Shtygashev A.A. Formation of the quasistationary state by scattering of wave packets on a finite lattice // Phys. Rev. B.

2008. V.77. P.075326 (lip.).

13. Peisakhovich Yu.G., Shtygashev A.A. Formation of quasi-stationary state by Gaussian wave packet scattering on the lattice of N identical deltapotentials // Phys. Rev. B. 2008. V.77, P. 075327 (9p.).

14. Штыгашев А.А. Применение неравномерной сетки в численном анализе спектральной задачи Штурма-Лиувилля. // «Актуальные проблемы электронного приборостроения» / Материалы IX международной конференции. Новосибирск, 2008. Т.6. С.117-121.

15. Штыгашев А.А. Распад квазистационарного состояния в решетке дельта-барьеров // Математическое моделирование. 2009. Т.21. №5. С. 67-76.

16. Штыгашев А.А. Осцилляции тока из двойной квантовой ямы при распаде квазистационарного состояния. // Известия вузов. Физика.

2009. Т.52. №9. С.46-50.

17. Штыгашев А.А. Коллапс резонансов в сверхрешетках. // Известия вузов. Физика. 2010. Т.53. №8. С.17-23.

18. Штыгашев А.А., Пейсахович Ю.Г. Рассеяние ультразвуковой волны ферромагнитной пластиной // «Актуальные проблемы электронного приборостроения» / Материалы X международной конференции. Новосибирск, 2010 Т.6. С.29-33.

19. Штыгашев А.А. Коллапс резонансов в симметричной двойной сверхрешетке (CaF2/Si)n/CaF2 конечной длины // Международная конференция с элементами научной школы для молодежи «Нанофизика и наноэлектроника» «Мезоскопические структуры в фундаментальных и прикладных исследованиях». Новосибирск, 2010. С.26-30.

20. Штыгашев А.А. Гибридный алгоритм в методе трансфер-матрицы. // Актуальные проблемы электронного приборостроения АПЭП-2010 / Материалы X международной конференции. Новосибирск, 2010. Т.6. С.161-165.

21. Штыгашев A.A. Математическое моделирование стационарных состояний частицы в многослойной структуре во внешних силовых полях // X Всероссийская научная конференция «Краевые задачи и математическое моделирование». Новокузнецк, 2010. С.115-122.

22. Peisakhovich Y.G., Shtygashev A.A. Formation of quasistationary states in a ferromagnetic plate by scattering of ultrasound impulse in the vicinity of magnetoacoustic resonance // J. Appl. Phys. 2011. V.110. N.5. P.053904

23. Штыгашев A.A. Прохождение ультразвуковой волны через магнитную сверхрешетку конечной длины // ФТТ. 2011. Т.53. №5. С.964-970.

24. Штыгашев A.A., Овчинников Ю.Э. Моделирование процесса переноса заряда в системе «фотосенсибилизатор-наночастицы-коллектор» с помощью волнового пакета // Научный вестник НГТУ. 2011. №2(43). С.107-114.

25. Штыгашев A.A. Применение гибридного алгоритма в спектральной задаче Штурма-Лиувилля для многослойных структур во внешних силовых полях // Научный вестник НГТУ 2011. №1(42). С.77-88.

26. Штыгашев A.A. Рассеяние волнового пакета решеткой дельта-барьеров в электрическом поле // ДАН ВШ РФ 2011. №1(16). С.64-74.

27. Пейсахович Ю.Г., Штыгашев A.A. Осцилляции тока при фотоэмиссии через трехбарьерную структуру // ДАН ВШ РФ. 2011. №1(16). С.26-40.

28. Пейсахович Ю.Г., Штыгашев A.A. Взаимодействие ультразвукового импульса с ферромагнитной пластиной в условиях магнитоакустиче-ского резонанса // Журнал сибирского федерального университета. Математика и физика. 2011. №3. С.371-381.

Отпечатано в типографии Новосибирского государственного Технического университета 630092, г. Новосибирск, пр. К. Маркса, 20, тел./факс: (383) 346-08-57 формат 60x84 1\16, объем 2.25 пл., тираж 120 заказ № 109 подписано в печать 26.12.12 г.

Текст научной работы диссертации и автореферата по физике, доктора физико-математических наук, Штыгашев, Александр Анатольевич, Новосибирск

НОВОСИБИРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ

УНИВЕРСИТЕТ

На правах рукописи

05201350758 Штыгашев Александр Анатольевич

Пространственно-временные и спектральные

характеристики нестационарных волновых процессов в неоднородных конденсированных

средах

01.04.07 - физика конденсированного состояния

ДИССЕРТАЦИЯ на соискание ученой степени доктора физико-математических наук

Научный консультант д. ф.-м. н., доцент Пейсахович Юрий Григорьевич

Новосибирск - 2012

Содержание

ВВЕДЕНИЕ.................................. 8

1. Стационарные и нестационарные волновые процессы в многослойных пространственно-неоднородных структурах....... 13

1.1. Метод трансфер-матрицы..................... 13

1.2. Динамика квантовых волновых пакетов в многослойных пространственно-неоднородных структурах............. 17

1.2.1. Рассеяние волнового пакета периодической решеткой . 19

1.2.2. Возбуждение и распад квазистационарных состояний . 21

1.2.3. Зеркально-симметричные гетероструктуры. Явление коллапса резонансных пиков прозрачности......... 24

1.2.4. Распространение волнового пакета в решетке при наличии внешнего силового поля ............. 29

1.3. Перенос заряда в фотоэлектрохимической ячейке Гретцеля. Температурное уширение оптических линий примесных центров в аморфных средах...................... 34

1.4. Осцилляции электронного тока из многобарьерных гетерострук-тур, связанные с квазистационарными состояниями...... 42

1.4.1. Распад квазистационарного состояния из поверхностных многобарьерных гетероструктур.......... 42

1.4.2. Осцилляции тока фотоэмиссии через поверхностную трехбарьерную гетероструктуру............. 45

1.5. Дифракция связанных волн в многослойных структурах .... 47 1.5.1. Рассеяние ультразвуковой монохроматической волны

ферромагнитным слоем вблизи магнитоакустического резонанса.......................... 48

1.5.2. Рассеяние ультразвуковой монохроматической волны магнитной сверхрешеткой вблизи магнитоакустическо-

го резонанса........................50

1.5.3. Рассеяние ультразвукового пакета ферромагнитной пластиной вблизи магнитоакустического резонанса.....51

1.5.4. Волновые функции и спектр энергии возбуждений бо-голюбовских квазичастиц в переходе........53

Основные результаты и выводы главы 1.................55

2. Метод трансфер-матрицы и гибридный алгоритм для многослойных систем...............................57

2.1. Метод трансфер-матрицы для многокомпонентных волновых

полей в многослойных структурах................57

2.1.1. Метод трансфер-матрицы для случая N — 2.......60

2.2. Конечно-разностный метод трансфер-матрицы..........67

2.2.1. Трансфер-матрица, построенная с помощью явной схемы Эйлера.......................68

2.2.2. Трансфер-матрица, построенная с помощью неявной схемы Адамса для случая N = 2..............69

2.2.3. Трансфер-матрица, построенная с помощью неявной схемы Адамса для случая N > 2..............71

2.3. Гибридный алгоритм в методе трансфер-матрицы.........71

Основные результаты и выводы главы 2..................72

3. Динамика квантовых волновых пакетов вмногослойных пространственно-неоднородных структурах.................73

3.1. Рассеяние волнового пакета на периодической решетке.....73

3.1.1. Численно-аналитическое описание волнового пакета в решетке конечной длины..................74

3.1.2. Рассеяние гауссова волнового пакета решеткой из N одинаковых дельта-барьеров..................90

3.2. Формирование и эволюция квазистационарных состояний . . . 104

3.2.1. Образование и распад квазистационарного состояния при рассеянии волнового пакета на решетке из N одинаковых дельта-барьеров ..................106

3.2.2. Положение нулей и полюсов амплитуды отражения . . 112

3.2.3. Распад квазистационарных состояний из решетки N дельта-барьеров....................118

3.3. Рассеяние волнового пакета в условиях коллапса резонансов в симметричной двойной сверхрешетке................125

3.3.1. Симметричная двойная сверхрешетка...........126

3.3.2. Коллапс резонансов в сверхрешетке 5-барьеров .... 128

3.3.3. Коллапс резонансов в сверхрешетке (CaF2/Si)n/CaF2 . 135

3.3.4. Рассеяние волнового пакета в условиях коллапса резонансов .........................140

3.4. Распространение волнового пакета в решетке при наличии внешнего силового поля ..........................141

3.4.1. Стационарные состояния электрона в периодических

структурах конечной длины, помещенных в электрическое поле.........................142

3.4.1.1. Решетка Кронига-Пенни в электрическом поле ........................143

3.4.1.2. Решетка ö-барьеров в приближении кусочно-постоянной аппроксимации потенциала внешнего электрического поля............147

3.4.1.3. Решетка б-барьеров конечной длины в электрическом поле. Расчет стационарных состояний при помощи гибридного алгоритма 159

3.4.1.4. Перестройка резонансных состояний полосы пропускания при увеличении напряженности внешнего электрического поля .... 163

3.4.2. Рассеяние волнового пакета на решетке в электрическом поле..........................167

Основные результаты и выводы главы 3..................171

4. Электронные переходы в примесных центрах аморфной среды и туннельный перенос заряда через мезопористый слой наночастиц 175

4.1. Фотоэлектрохимическая ячейка Гретцеля.............175

4.2. Влияние туннельных переходов на форму оптических полос примесных центов в полимерах и стеклах..............178

4.2.1. Модельный гамильтониан, полуширина и сдвиг бесфо-нонной линии примесной молекулы............178

4.2.2. Однородная полуширина БФЛ примесного центра в твердом растворе......................182

4.3. Перенос заряда в цепочке «молекула сенсибилизатора — слой наночастиц — коллектор»......................189

4.3.1. Математическая модель инжекции и переноса электрона189

4.3.2. Влияние буферного слоя наночастиц на перенос фотовозбужденного электрона.................190

Основные результаты и выводы главы 4.................195

5. Осцилляции электронного тока из многобарьерных гетерострук-тур, связанные с квазистационарными состояниями.......197

5.1. Распад квазистационарного состояния в двойной квантовой яме 197

5.2. Фотоэмиссия через поверхностную гетероструктуру.......203

5.2.1. Эмиссионный ток.....................205

5.2.2. Матрица плотности в теории нестационарной фотоэмиссии ............................208

5.2.3. Стационарная фотоэмиссия..................213

5.2.4. Осцилляции электронного тока при фотоэмиссии через поверхностную трехбарьерную гетероструктуру.....216

Основные результаты и выводы главы 5..................220

6. Связанные волны в многослойных структурах. Магнитоупругие волны в магнитных слоях и элементарные возбуждения в сверхпроводящих слоях пространственно-неоднородных структур .... 222

6.1. Рассеяние ультразвуковой монохроматической волны, проходящей через магнитный слой вблизи магнитоупругого резонанса 223

6.1.1. Описание модели.....................223

6.1.2. Закон дисперсии магнитоупругих волн..........227

6.1.3. Эффективная трансфер-матрица...............228

6.1.4. Нули и полюсы амплитуд рассеяния............231

6.1.5. Спектры коэффициентов рассеяния...........235

6.2. Рассеяние ультразвуковой монохроматической волны, проходящей через магнитную сверхрешетку вблизи магнитоупругого резонанса...........................237

6.2.1. Постановка задачи.......................238

6.2.2. Волновые поля в й- и М-слоях...............238

6.2.3. Метод трансфер-матрицы.................239

6.2.4. Спектр коэффициента прохождения...........240

6.3. Рассеяние ультразвукового пакета ферромагнитной пластиной вблизи магнитоупругого резонанса................244

6.3.1. Постановка задачи..................245

6.3.2. Асимптотические оценки спектральных интегралов . . 247

6.3.3. Эволюция и распад магнитоупругих возбуждений в магнитном слое..................250

6.4. Моделирование спектра андреевских состояний в SN структурах262

6.4.1. Андреевские состояния в SN кольцах........264

6.4.2. Андреевские состояния в Б,N8., структурах......265

Основные результаты и выводы главы 6...............270

ЗАКЛЮЧЕНИЕ...........................273

ПРИЛОЖЕНИЯ............................277

ПРИЛОЖЕНИЕ А. Гибридный алгоритм для решения задач Коши и Штурма-Лиувилля (к главе 2)....................277

ПРИЛОЖЕНИЕ В. Полиномы Чебышева второго рода (к главе 3) 281

ПРИЛОЖЕНИЕС. Алгоритм вычисления энергийиволновых функций стационарных состояний частицы в решетке б-потенциалов и решение нестационарного уравнения Шредингера (к главе 3) . . 283

ПРИЛОЖЕНИЕ Б. Вычисление базисных функций {ф^ (к главе 4) 287

Литература

289

ВВЕДЕНИЕ

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

структурах «магнетик/немагнетик».

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

Научные положения и результаты, выносимые на защиту.

1. Результаты моделирования распространения волновых пакетов в многослойных структурах в отсутствие и при наличии внешних силовых полей.

Диаграммы локальной устойчивости стационарных функций электрона в решетке в сильном электрическом поле. Асимптотические и численные оценки характеристик распространения пакетов, формирования и распада квазистационарных состояний в зависимости от длины решетки, спектрального состава пакета и его положения по отношению к полосам пропускания. Показано, что время жизни пороговых квазистационарных состояний пропорционально третьей степени длины решетки, а для квазистационарных состояний средней части полосы пропускания время жизни пропорционально длине решетки. Рас-падный ток из многобарьерных структур может осциллировать с частотами, определяемыми энергиями смешиваемых квазистационарных состояний.

ферромагнитной пластине и магнитной сверхрешетке в области магнитоаку-стического резонанса (MAP) при произвольных условиях закрепления спинов на поверхности. Показано, что при резонансных частотах вблизи MAP происходит полное отражение ультразвука (антирезонансы Фано), вызываемое взаимодействием упругой волны со спиновыми волнами, локализованными в магнитных слоях. Изменением параметров исходного волнового пакета можно целенаправленно возбуждать в магнитных слоях квазистационарные состояния и вторичные импульсы, обусловленные интерференцией и перекачкой энергии между упругой и спиновой подсистемами.

6. Численно-аналитическое исследование спектра андреевских состояний в несимметричном S \NS2 переходе при разных поперечных импульсах Ферми в слоях. Показано, что при увеличении разности поперечных импульсов Ферми фазовая зависимость спектра энергий стационарных андреевских состояний ослабевает.

Структура диссертации

Диссертация состоит из введения, шести глав, заключения, четырех приложений и списка литературы.

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

да квазистационарных состояний в решетке. В четвертой главе рассмотрены процессы температурного уширения оптических линий излучения примесных молекул в аморфных средах, инжекции фотовозбужденных электронов в полупроводниковые наночастицы и переноса электронов в мезопористом слое наночастиц. Пятая глава посвящена описанию процессов распада возбужденных электронных состояний из структур, ограниченных потенциальными барьерами, а также описанию нестационарной фотоэмиссии при накачке и распаде возбуждений через трехбарьерную поверхностную гетероструктуру. В шестой главе проведено исследование связанных волн (N>2). В первых трех разделах рассмотрены процессы стационарного и нестационарного рассеяния ультразвуковых волн на магнитном слое и магнон-фононной сверхрешетке вблизи MAP без и с учетом затухания упругих и спиновых волн. Исследованы полюсы и нули амплитуд рассеяния упругих волн с учетом взаимодействия со спиновыми волнами в магнитных слоях. Изучена нестационарная задача рассеяния упругих волновых пакетов на магнитном слое в области MAP при разных толщинах магнитного слоя и при жестком закреплении спинов на поверхности раздела сред. В четвертом разделе решается стационарная задача рассеяния связанных боголюбовск