Эффект Ганна - Левинштейн М.Е.
Скачать (прямая ссылка):
(2.13), устанавливаем различные случайные моменты рассеяния электрона ta.
Теперь необходимо определить, какой именно процесс рассеяния произошел в
момент времени t0. Для того чтобы ответить на этот вопрос, выкинем еще
одно случайное число d (принадлежащее к случайным числам, распределенным
с одинаковой плотностью вероятности между 0 и Г). После этого проверим
выполнение серии неравенств по следующей схеме. Если с?<Ао, то будем
считать, что произошел фиктивный процесс рассеяния, если то произошел
пер-
вый процесс рассеяния, если ta + Ai^d^Ao + Ai-j-Aa, второй и т. д. Из
этой схемы понятно, что выгодно выбирать величину Г возможно меньшей (но
все же достаточно большой, чтобы обеспечить положительность Ао) для того,
чтобы фиктивный процесс рассеяния происходил возможно реже.
После того, как определено, какой именно процесс рассеяния произошел,
необходимо определить конечное состояние электрона, т. е. значение его
энергии <§ и импульса k' после рассеяния.
Для фиктивного процесса рассеяния конечное состояние, разумеется,
известно - оно совпадает с начальным. Для физических же процессов
рассеяния энергии конечное состояние <§ (k') определяется законом
сохранения энергии. Значение волнового вектора k' в конечном состоянии
зависит от распределения вероятности волновых векторов конечного
состояния по углам. Оно в свою очередь определяется механизмом рассеяния.
Для того чтобы определить значение k', также используют случайные числа.
Рассмотрим в качестве примера простейший случай: рассеяние электрона на
акустических фононах в параболической и сферически симметричной долине
зоны проводимости. Для модели зонной структуры GaAs, используемой при
расчетах методом Монте-Карло, функция распределения электронов
симметрична относительно направления электрического поля. Поэтому
достаточно характеризовать волновой вектор k двумя составляющими - kz,
направленной вдоль поля, и перпендикулярной полю составляющей ?р. Для
акустического рассеяния в параболической долине зоны проводимости все
значения волновых векторов k' равновероятны. Поэтому вероятность того,
что вектор k' будет составлять угол 0 с направлением поля Z,
пропорциональна длине окружности радиусом \k'\ sin6. Отсюда вероятность
того, что волновой вектор к' будет составлять с осью Z угол, лежащий
между значениями 0 и Q + dQ, равна
где константа 7'2 определена из условия нормировки jP(0)c?0 = l. Из
о
(2.14) видно, что можно определить угол 0 и, следовательно, конечное
состояние k', генерируя случайные числа Ь, распределенные с равной
плотностью в интервале между 0 и 1 при использовании соотношения
P(Q)dQ = ilz sin QdQ,
(2.14)
b = j" P (6') db = */* 0 - cos e)-
0
(2.15)
25
В более сложных случаях, например при рассеянии электрона на полярных
оптических фононах, необходимо учитывать угловую зависимость плотности
распределения вероятности конечных состояний k',. что усложняет описанную
выше процедуру [5]. Основная идея расчета остается, однако, той же:
определить конечное состояние электрона, используя случайные числа. После
того, как волновой вектор электрона в конечном состоянии определен, мы
вновь выкидываем случайное число г, определяя случайный момент следующего
рассеяния. Затем выбрасываем случайное число d, определяя, какой именно
процесс рассеяния произошел. Наконец, выкидывая одно (для акустического
рассеяния) или несколько (в более сложных случаях) случайных чисел,
определяем конечное состояние электрона. Затем весь процесс повторяется и
т. д.
Таким образом, осуществляется моделирование движения электрона в
импульсном пространстве. При этом подсчитывается время, которое проводит
электрон в каждом малом объеме импульсного пространства. Это время будет
пропорционально функции распределения электронов f(k). Зная функцию
распределения электронов, легко подсчитать все физические наблюдаемые
величины, такие, как средняя дрейфовая скорость электронов
р= / vhf(k)dk, (2.16)
где Vh - скорость электрона в состоянии k, и другие*).
Можно показать, что описанный выше метод Монте-Карло эквивалентен точному
решению уравнения Больцмана для невырожденного-электронного газа при
заданных скоростях рассеяния Sq(k, k') и заданной зависимости энергии
электронов от импульса <§ {к) [5].
2.2-2. Результаты расчетов эффектов переноса в GaAs, выполненных методом
Монте-Карло
Наиболее подробные расчеты характеристик эффектов переноса, выполненные
методом Монте-Карло, содержатся в работах {3-5]. Во всех этих работах в
основу расчета положена многодолинная модель зоны проводимости GaAs,
описанная в приложении I. Поверхности равной энергии в нижней долине
считались сферическими, однако з работе [5] была учтена непараболичность
нижней долины [7].
Далее, принималось, что верхние <100> долины, переход электронов в
которые обусловливает отрицательную проводимость в сильном поле,
относятся к точке Xj зоны Бриллюэна. Это предположение согласуется с
данными по изучению фотоэмиссии GaAs [8, 9], с результатами по измерению
