Материал: Эффект фотонной лавины в кристаллах и наноструктурах. Монография (Перлин), 2007, c.120

лавины, описанной в § 1.1, то n-фотонные переходы v → c играют роль слабых нерезонансных переходов 1 → 2, l-фотонные переходы c → с1 соответствуют быстрым резонансным переходам 2 → 3, а переходы c1 + (n –

l)hω → ccv играют ту же роль, что кросс-релаксация.

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

Рассматриваемая модель, естественно, включает также релаксационный процесс c1 → c с участием фононов и процессы рекомбинации нерав-

новесных фотовозбуждённых электронов и дырок.

Очевидно, что предложенная модель может работать лишь в области очень сильных электромагнитных полей (Fω t 107 В/см). При hω = 1.17 эВ,

принимает значение γ ≈ 6.4. Большая величина параметра γ позволяет заключить, что в области интенсивностей меньших 1013 Вт/см2 в рассматриваемой ситуации адекватной является картина многофотонных межзонных переходов, а не межзонного туннелирования в поле сильной электромаг-

нитной волны. При этом неравенство 1(4γ )2 <<1 позволяет утверждать,

что между заданными зонами фактически будут идти переходы с наименьшим числом фотонов, допускаемым законами сохранения3 . Для оценки вероятностей многофотонных межзонных переходов далее будут использоваться результаты работ [86], [87], где получены формулы для ве-

роятностей n-фотонных межзонных переходов Wvc(n) ≡σvc(n) jn при произ-

вольных n. В области малых j для n = 3 эти формулы совпадают с полученными в рамках стандартной теории возмущений. При n >> 1 формулы из [86, 87] дают асимптотику, близкую к той, что получена Л.В. Келдышем в [84] и Ю.А. Бычковым и А.М. Дыхне [85].

Как уже отмечалось выше, ключевую роль для запуска процесса многофотонной лавины играют переходы оже-типа c1 + shω → ccv ,

s = n −l . Вычислению вероятностей этих переходов Wa(s) ≡ γ (s) jsnc1 посвящен следующий параграф.

§ 5.3. Расчет вероятностей многофотонных переходов оже-типа

Представим гамильтониан электрон-фотонной системы в виде

3 Такая ситуация является типичной, но бывают исключения. Так в работе [110] показано, что в непрямозонном материале AgBr благодаря особенностям электронной зонной структуры и проявлениям оптического эффекта Штарка в условиях двойного межзонного резонанса вероятности n-фотонных межзонных переходов с n = 3 и n = 4 оказываются в при j t 1010 Вт/см2 (λ = 560 нм) бóльшими, чем вероятности переходов с n = 2.

73

– гамильтонианы невзаимодействующих подсистем свободных электронов и фотонов,

квадратичная по полю часть электрон-фотонного взаимодействия, отличная от нуля лишь при учете малого волнового вектора фотона κ,

между зонами. В формулах (5.4-5.10) использованы следующие обозначения: m – масса свободного электрона; V0 – объем элементарной ячейки;

εi(k) – блоховская энергия электрона в i-ой зоне; ξi+k , ξik – операторы рождения и уничтожения для состояний с волновым вектором k в i-й зоне, которые заменяются на электронные операторы ai+k , aik или дырочные операторы bik , bi+k , когда индекс i пробегает значения, соответствующие зонам

проводимости либо валентным зонам; κ, ωκ, λ и eλ – волновой вектор, частота, индекс поляризации и единичный вектор поляризации фотонов;

cκ+λ ,cκλ – операторы рождения и уничтожения фотонов; εL(ω,q) и εT(ω,q) – продольная и поперечная диэлектрические проницаемости, зависящие от частоты и волнового вектора, переданная при межэлектронных столкнове-

ниях частота Ω% дается выражением Ω% = h−1[εi (k0 ) −εi′(k1)].

74

Оценка матричных элементов операторов He′-phot и He′-e с использова-

нием k p - теории возмущений приводит к следующему результату: отношение внутризонного матричного элемента оператора электрон-фотонного взаимодействия He′-phot к межзонному матричному элементу пропорцио-

нально величине βkii′ <<1, а в случае кулоновского взаимодействия He′-e отношение межзонного матричного элемента к внутризонному пропорционально величине βqii′ [см. формулы (5.6-5.10)].

Выражение для вероятности многофотонного перехода из зоны v в зону c с участием свободных носителей из зоны c1 запишем в виде:

чает усреднение по состояниям фотонной подсистемы. Для усреднения используется т.н. диагональное представление Глаубера [102, 103] с оператором плотности электромагнитного поля

где ξκ и ξκ собственные значения и собственные функции оператора cκ.

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

где nκ – среднее число фотонов в моде.

В формуле (5.11) учтено, что из законов сохранения энергии и импульса следует, что конечные состояния электронов в зоне c удалены от дна зоны и могут считаться незаполненными.

Матричные элементы перехода будем вычислять в (l + 1)-м порядке теории возмущения – l порядков по взаимодействию He′-phot и один порядок

по взаимодействию He′-e . Вводя обозначение

(5.14)

запишем критерий применимости борновского приближения по межэлектронному взаимодействию He′-e для рассматриваемой задачи в виде

75

где pcv – межзонный матричный элемент оператора импульса, Fω – амплитуда электрического поля световой волны, ωR(s) – s-фотонная частота Раби (s = 2 или 3) для переходов между зонами c и c1, τp – время релаксации импульса электронов, ωR(s) js2 , где j – интенсивность излучения. Первое из

условий (5.16) выполняется для всех актуальных интенсивностей, тогда как второе при s = 2 выполняются лишь при j d 1012 Вт/см2. При более высоких интенсивностях, вообще говоря, требуется учет эффектов переизлучения фотонов.

Рис. 5.2. Диаграммы, представляющие прямой вклад в процесс в случае модели А: сплошные линии со стрелками вверх (вниз) – электроны (дырки), волнистые – фотоны, пунктирные – кулоновское взаимодействие

Из множества диаграмм Фейнмана выберем те, которые вносят основной вклад в процесс. В случае модели А это диаграммы, в которых присутствуют только межзонные матричные элементы оператора электронфотонного взаимодействия и внутризонный матричный элемент оператора кулоновского взаимодействия. На рис. 5.2 представлены диаграммы, дающие прямой вклад в составной матричный элемент. Для модели В при выборе диаграмм считаем, что переходы из зоны v в зону c1 запрещены. Как показывает анализ, в этом случае удается выделить две диаграммы, вклад

которых в (βqii′)−2 раз больше вклада остальных диаграмм (рис. 5.3).

76

Рис. 5.3. Диаграммы, представляющие прямой вклад в случае модели В (обозначения такие же, как на рис. 5.2)

Для получения формул в виде, удобном для численных расчетов, используем следующие упрощающие предположения: будем считать, что зоны v, c, c1 являются параболическими с экстремумами в центре зоны Бриллюэна в точке k = 0, а электроны в зоне c1 сосредоточены вблизи экстремума зоны, так что можно положить k0 = 0 .

Прямой вклад составного матричного элемента запишем в виде:

элементы оператора импульса, A0 – амплитуда вектор-потенциала свето-

вой волны, с – скорость света в вакууме, ζ = mc mv , ζ ′ = mc mc1 , mi – эффективная масса частиц в i-й зоне. Выражение для обменного вклада составного матричного элемента можно получить из формул (5.17), (5.18) путем замены k1 ↔ k2 .

77