Методичка Зубкова
.pdf
распределения электронов по скоростям |
|
||||
симметрична |
относительно оси |
орди- |
f |
||
|
|
|
|
|
|
нат: dn vx dn vx , и никакой ток в |
F0 |
||||
0 |
|||||
f |
|||||
этих условиях не может существовать. |
F |
||||
|
|||||
В условиях |
внешнего |
возмущения |
Ff11 |
||
(например, приложенной разности по- |
|||||
тенциалов) все электроны получают до- |
|
||||
полнительную скорость в направлении |
p,, k |
||||
vx , при этом возникает дрейфовый ток, |
EE |
||||
внеш |
|||||
и функция распределения смещается в |
Рис. 4.3. Равновесная и неравновесная |
||||
функции распределения электронов |
|||||
соответствующем |
направлении. |
Если |
по скоростям для невырожденного |
||
возмущение носит стационарный ха- |
электронного газа |
||||
рактер, то устанавливается динамическое равновесие, когда действия хаоти- |
|||||
ческих соударений и внешнего поля компенсируют друг друга: |
|||||
df |
|
f |
|
|
|
f |
|
|
0 . |
|
|
||||||||
dt |
|
t |
|
пол |
|
t |
|
столк |
|
|
|
|
|
|
Как правило, дрейфовая добавка к скорости электрона в поле мала по сравнению с тепловой скоростью, поэтому неравновесная и равновесная функции распределения отличаются на малое приращение f1 :
ff0 f1.
Всилу аддитивного характера кинетического уравнения Больцмана все кинетические коэффициенты могут быть рассчитаны через f1 .
Из-за столкновения носителей зарядов с атомами кристаллической решетки их движение в направлении действия электрического поля будет прерывистым и может быть охарактеризовано средней скоростью. Средняя скорость движения носителей заряда в поле единичной напряженности называется подвижностью носителей заряда. В приближении времени релаксации подвижность выражается формулой [5]*
q m * . |
(4.13) |
Подвижность электронов и дырок в полупроводнике зависит от зонной структуры материала и от механизма их рассеяния в кристаллической решетке. Характер зависимости подвижности от температуры, концентрации примесей, а также смысл эффективной массы m* в (4.13) рассматриваются в следующих главах.
*См. также: Ю П., Кардона М. Основы физики полупроводников. М.: Физмалит, 2002. 61
vxdt |
Покажем типовое решение кинетиче- |
|
|
ской задачи с помощью неравновесной |
|
|
функции распределения f ( p,r, t) . Для этого |
|
|
вычислим плотность тока в точке r образца |
|
x |
в момент t [4]. Выделим единичную площад- |
|
ку, перпендикулярную оси x (рис. 4.4). Чис- |
||
Рис. 4.4. К расчету плотности тока |
||
ло электронов внутри цилиндра длиной vxdt , |
построенного на этой площадке, равно f ( p,r, t)vxdtdvxdvydvz . Все эти элек-
троны за время dt пересекут рассматриваемую площадку. Полное число электронов всех скоростей, пересекающих площадку за время dt , равно
dt f ( p,r, t)vxdvxdvydvz ,
при этом учитываются электроны, пересекающие площадку как слева направо, так и справа налево. Так как каждый электрон имеет заряд –q, то плотность тока в направлении x выражается интегралом
jx q f ( p,r, t)vxdvxdvydvz .
Так как f f0 f1, а f0 (E ) – четная функция vx , то интеграл по dvx от f0 (E )vx равен нулю и окончательно получим:
jx q f1( p,r,t)vxdvxdvydvz .
Таким образом, при рассмотрении кинетических явлений основной задачей является определение неравновесной функции распределения f ( p,r,t)
(или приращения f1( p, r,t) ). Эта задача решается при условии, что известен
интеграл столкновений (либо время релаксации), следовательно, необходимо знание механизма рассеяния носителей заряда.
4.3.Контрольные вопросы
1.Объясните физический смысл кинетического уравнения Больцмана.
2.Что является решением кинетического уравнения Больцмана?
3.Раскройте понятие фазового пространства.
4.Что характеризует неравновесная функция распределения?
5.Что такое интеграл столкновений?
6.В чем заключается приближение времени релаксации?
7.Что такое подвижность носителей заряда?
62
5.МЕХАНИЗМЫ РАССЕЯНИЯ НОСИТЕЛЕЙ ЗАРЯДА
ВПОЛУПРОВОДНИКЕ
5.1.Основные механизмы рассеяния. Общие положения
Движущиеся в кристалле электроны проводимости характеризуются
средней кинетической энергией E , которая имеет величину порядка 32 kBT , и
соответствующей этой энергии длиной волны |
h |
. Если электрон- |
* |
||
|
3kBTmn |
|
ная волна встречает на своем пути дефект, протяженность которого соизмерима с
4 , то она эффективно рассеивается, т. е. частично изменяет или
полностью утрачивает направленную составляющую скорости. Рассеяние носителей заряда в твердых телах – процесс взаимодействия электрона (дырки) с нарушениями идеальной периодичности кристалла, приводящий к потере направленной составляющей скорости. В результате акта рассеяния электрон переходит из состояния с импульсом p в состояние с импульсом p , кроме
того, также может измениться и его энергия.
Рассеяние называется упругим, если энергии электрона в начальном и
конечном состояниях равны: |
|
, или неупругим, если |
|
E p E p |
|||
|
. |
|
|
E p E p |
|
|
|
Есть триосновных источника, приводящих крассеянию носителей заряда:
1.Статические дефекты решетки. К ним относятся точечные дефекты: примесные атомы замещения, вакансии, междоузельные атомы. Кроме того,
кним относятся одномерные дефекты – дислокации, а также двумерные дефекты – межкристаллитные границы в поликристалле, интерфейсы гетероструктур и поверхность кристалла. Статические дефекты рассмотрены в гл. 2.
2.Динамические дефекты решетки. Эти дефекты связаны с движением атомов кристалла (фононы) или их спинов (магноны).
3.Кулоновское взаимодействие квазисвободных носителей заряда друг с другом [6].
Пусть свободные электроны, которые характеризуются концентрацией
n, движутся со средней скоростью vт в данном направлении. Тогда nvт есть плотность потока электронов, т. е. количество электронов, проходящих в единицу времени через единичную площадку, перпендикулярную направле-
63
нию их скорости. Положим, что на пути потока существуют центры рассеяния с концентрацией Nдеф . Для количественной оценки вероятности рассея-
ния электрона на дефекте вводится эффективное сечение рассеяния (σ) – отношение числа электронов, удаленных из пучка в результате рассеяния на одном центре в единицу времени, к плотности падающего потока частиц. Для точечных дефектов величину сечения рассеяния можно оценить как площадь квадрата со стороной в несколько единиц или десятков постоянных решетки
(10–12…10–15 см2).
Количество рассеянных в единицу эффективным сечением, концентрацией стью падающего потока электронов nvт ,
времени электронов определяется рассеивающих центров и плотно- т. е.
n1 Nдефnvт . |
(5.1) |
Если W – вероятность рассеяния одной частицы в единицу времени, то количество рассеянных в единицу времени электронов есть
n1 Wn . |
(5.2) |
Из сопоставления (5.1) и (5.2) получим выражения для эффективного сечения рассеяния частицы:
W
Nдефvт
идля вероятности рассеяния (выражается в секундах в минус первой степени) одной частицы:
W N дефvт . |
(5.3) |
||||
Величина, обратная W, называется временем релаксации носителей заряда: |
|||||
1 |
1 |
|
|
||
|
|
|
|
. |
(5.4) |
W |
N v |
||||
|
|
|
деф т |
|
|
Рассеяние носителей заряда является причиной того, что любое неравновесное по энергии или импульсу распределение электронов, созданное внешним возмущением, релаксирует с характерной постоянной времени τ к равновесному фермиевскому распределению f0 E , соответствующему тем-
пературе кристалла Т.
Далее рассмотрим различные механизмы рассеяния носителей заряда.
64
5.2. Рассеяние носителей заряда на динамических дефектах
Основные сведения о тепловых колебаниях решетки. При темпера-
туре выше абсолютного нуля атомы в кристалле колеблются около своих положений устойчивого равновесия. Так как все атомы твердого тела связаны друг с другом упругими силами, то колебания любого из них передаются соседнему, и, таким образом, по всему кристаллу во всевозможных направле-
ниях распространяются упругие волны. Каж- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
Е |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
дую волну можно характеризовать квазивол- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
новым вектором q и частотой ω, зависящей |
|
|
ЕС |
|
|
|
|
|
|
|
|
|
|
|
|||||||
от q : j q , где индекс j = 1, 2, ..., 3r (r – |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
ЕF |
|
|
|
|
|||||||
число атомов в элементарной ячейке кри- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
сталла) нумерует ветвь колебания. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|||||
Аналогично тому, |
как для электронов |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
ω |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
зависимость энергии от квазиволнового век- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
LO |
|
|
|
|
|
||||||||
тора |
рассматривают в |
приведенной зоне |
|
|
|
|
|
|
|
|
|
|
|
|
|
TO2 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Бриллюэна, дисперсионные зависимости ча- |
|
|
|
|
|
|
|
|
|
|
|
|
|
TO1 |
|||||||
|
|
|
|
|
|
|
LA |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
стоты |
колебаний решетки также принято |
|
|
|
|
|
|
|
|
|
|
|
|
|
TA2 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
TA1 |
|
|||||||
строить в пределах первой зоны Бриллюэна. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
В элементарных ячейках полупроводников |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
L |
|
|
Г |
|
|
|
|
|
|
X |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
q |
|
||||||||||
типа алмаза и цинковой обманки содержатся |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Рис. 5.1. Дисперсионные кривые |
|||||||||||||||||||||
два атома, поэтому в них имеется 6 различ- |
|||||||||||||||||||||
ных типов колебаний решетки (рис. 5.1). |
фононного спектра для трехмерной |
||||||||||||||||||||
решетки с базисом (cверху – |
|||||||||||||||||||||
Три нижние ветви называют акустиче- |
|||||||||||||||||||||
дисперсионная зависимость |
|||||||||||||||||||||
скими, остальные – оптическими. Оптиче- |
|
электрона в прямозонном |
|||||||||||||||||||
ские и акустические ветви разделены поло- |
|
|
полупроводнике) |
||||||||||||||||||
сой запрещенных частот, так что во всем интервале изменений волновых чисел частота оптических колебаний больше
частоты акустических. Для каждого типа (акустических или оптических) колебаний в трехмерном кристалле возможны три поляризации колебаний: одна продольная (LA, LO) и две поперечные (TA, TO). Для акустических мод соседние атомы в элементарной ячейке движутся в одном направлении (в фазе). Для оптических мод колебаний соседние атомы движутся в противофазе.
Энергия колебаний решетки квантована. Квант энергии тепловых колебаний решетки (фонон) обладает энергией Eф и квазиимпульсом
Pф
.
65
Как известно из квантовой механики, электроны являются фермионами и подчиняются статистике Ферми–Дирака (3.3). В отличие от этого фононы являются бозонами, т. е. в кристаллической решетке может быть возбуждено одновременно сколь угодно одинаковых фононов. Число фононов в кристалле является переменной величиной и определяется из условий теплового равновесия. В состоянии термодинамического равновесия среднее число фононов данного типа с частотой j и волновым вектором q (т. е. заполненность
фононной моды) выражается распределением Бозе–Эйнштейна
n jq |
|
1 |
|
|
, |
(5.5) |
|
jq |
|
||||
|
1 |
|
||||
|
exp |
|
|
|
||
|
k T |
|
||||
|
|
B |
|
|
|
|
где j – номер фононной моды (LA, TA и т. д., см. рис. 5.1).
На каждое колебание приходится средняя энергия jq n jq . Строго говоря, к этой энергии надо прибавить энергию нулевых колебаний, равную
12 jq . Но энергией нулевых колебаний кристалл обладает всегда, и ее
обычно принимают за начало отсчета.
Можно показать, что максимальная частота фонона для кристалла, элементарная ячейка которого состоит из двух атомов с массами m1 и m2 , вы-
числяется как max |
m1 m2 |
, где |
– упругая постоянная кристалла, а |
|
|||
|
m1m2 |
|
|
максимальная энергия max kB D , где D – температура Дебая. Фактиче-
ски температура Дебая разделяет две температурные области. В области низких температур на энергию и теплоемкость решетки сильное влияние оказывают квантовые эффекты («вымерзание» высокочастотных колебаний). В области высоких температур эти эффекты несущественны и теплоемкость может быть вычислена в классическом приближении. Для большинства кристаллов D лежит в интервале от 100 до 300 К.
Чтобы вычислить среднюю энергию всех колебаний кристаллической решетки, нужно просуммировать среднюю энергию всех типов колебаний (всех состояний фононов):
E |
jq n jq |
|
jq |
|
|
. |
|
|
|
jq |
|
||||
j,q |
j,q exp |
|
|
|
1 |
||
|
k T |
||||||
|
|
|
|
B |
|
|
|
|
66 |
|
|
|
|
|
|
Проще всего данную сумму вычислить при высоких температурах, когда для частот всех колебаний выполняется неравенство jq kBT (классиче-
ский предел). Тогда средняя энергия, приходящаяся на каждое колебание, равна kBT , всего колебаний 3rN (N – число примитивных ячеек кристалла в единице объема) и для полной энергии получаем E 3rNkBT .
Механизм взаимодействия носителей заряда с фононами. При взаимо-
действии с колебаниями решетки энергия электрона изменяется порциями, крат-
ными энергии соответствующего кванта колеба- |
|
|
|||
ний решетки. В рассеянии частицы могут участ- |
kx |
|
|||
вовать несколько фононов, однако наиболее ве- |
q |
kz |
|||
роятен однофононный процесс. При этом долж- |
|||||
ны быть выполнены законы сохранения энергии и |
k2 |
1 |
|||
|
|||||
импульса(рис. 5.2): |
|
|
|
|
|
|
|
|
|
ky |
|
E k E k |
|
; |
|
||
2 |
|
|
|||
|
1 |
|
|
|
|
|
k1 q k2. |
|
|
|
|
Поглощение или испускание i фононов приводит к изменению числа заполнения фононной моды n jq i n jq .
Рис. 5.2. Изменение волнового вектора электрона в процессе упругого рассеяния на фононе
Как было показано в предыдущих главах, невырожденный электронный газ заполняет узкий диапазон энергетических состояний ( ~kBT ) и квазивол-
новых векторов вблизи абсолютного экстремума зоны проводимости (см. рис. 5.1). Поэтому при рассмотрении различных механизмов рассеяния следует иметь в виду, что наиболее эффективно в конкретном процессе будут участвовать те акустические либо оптические фононы, которые обладают соответствующими значениями энергии и квазиволнового вектора.
5.3. Рассеяние на акустическом деформационном потенциале
Продольные акустические волны приводят к локальным изменениям постоянной решетки, создавая в кристалле участки сжатия и разрежения, что в соответствии с рис. 1.3 вызывает локальные колебания ширины запрещенной зоны (рис. 5.3). Эти колебания проявляют себя как наличие дополнительного деформационного потенциала, на котором происходит рассеяние носителей заряда. Данный процесс рассеяния носит фундаментальный характер, поскольку не связан с наличием разного рода несовершенств, а существует даже в идеальном кристалле.
67
Деформационный потенциал Edef есть, по определению, коэффициент
пропорциональности между величиной сдвига энергии соответствующей зоны EC и относительным изменением первоначального объема кристалла
EC Edef V
V . Анализируя рассеяние электронов на деформационном по-
тенциале решетки, Бардин и Шокли показали [7], что для простого изотропного случая зоны проводимости основную роль играет рассеяние электронов на длинноволновых колебаниях решетки (вблизи точки Г), длина волны которых существенно больше межатомного расстояния. Поэтому в качестве волновой
функции электрона можно использовать плоскую волну k (r) Ceikr , а блоховскую периодическую составляющую рассматривать как константу.
отр 2
1
EC
EV
x = 0
а |
б |
Рис. 5.3. Продольные акустические колебания решетки кристалла: a – вариации энергии дна зоны проводимости и потолка валентной зоны при наличии продольной акустической волны; б – рассеяние электрона на деформационном потенциале
Среднюю энергию движущихся электронов в невырожденном полупроводнике можно оценить как E 32 kBT . При не очень низких температурах
она намного превышает энергию акустических фононов jq в центре зоны
Бриллюэна. Поэтому рассеяние электронов на акустических колебаниях решетки, вызывающих изменение объема, можно рассматривать как упругое рассеяние, при котором электрон сохраняет свою энергию, но изменяет квазиимпульс, передавая или принимая его дискретными порциями от решетки.
Для оценки эффекта рассеяния на деформационном потенциале приведем рассуждения, рассматривая одномерный случай отражения плоской волны от прямоугольного барьера. Пусть частица, описываемая плоской волной
1 x eik1x , налетает на прямоугольный потенциальный барьер. Тогда су-
ществует вероятность, что она может либо отразиться, либо пройти через него, что выражается соответствующими волновыми функциями (рис. 5.3, б):
отр x re ik1x , 2 x teik2 x .
68
В силу общих условий, накладываемых на волновую функцию и ее производную, они должны оставаться непрерывными в точке разрыва потенциальной энергии. Поэтому при x 0 должны иметь место равенства
|
0 |
|
отр |
0 |
2 |
0 ; |
|
|
|
|
|
||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
0 |
|
|
|
0 0 , |
|
|
|
|
|
|||||||||||||
|
1 |
|
|
отр |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
из которых получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 r t; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
k rk |
|
tk |
2 |
. |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Решение системы имеет вид r |
k2 |
k1 |
, |
|
t |
|
|
2k1 |
|
. Для рассматривае- |
|||||||||||||
k |
|
|
k |
|
|
k |
|
|
|
||||||||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
k |
|
|
||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
1 |
|
2k . |
|||||
мого случая положим 2k k |
2 |
k и k k |
2 |
k |
|
, тогда r k |
|||||||||||||||||
|
|
1 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
По определению, коэффициент отражения есть отношение интенсивно-
стей отраженной и падающей волн, т. е. R r 2 k22 . Выражение для k
4k
можно получить, продифференцировав по k дисперсионную зависимость
|
2 |
2 |
|
m 2 E 2 |
|||
E |
|
k . Отсюда |
R |
|
|
|
. |
|
4 |
k |
4 |
||||
|
2m |
|
4 |
|
|
||
Подставив выражение EC Edef V
V для изменения энергии, вы-
званное деформационным потенциалом, получим R m*2 Edef 2 V / V 2 .
4 4k4
Можно показать, что при тепловом расширении V
V 2 kВT
V , где –
коэффициент пропорциональности (объемный модуль упругости). Анализируемый механизм отражения от барьера повторяется в кристал-
ле с периодичностью
4 . Поэтому в качестве рассматриваемого объема естественно взять V
4 3 . Тогда
|
m 2 |
E |
2 k T |
|
m 2 E |
2 k T 42 |
|
|
|
|
R |
|
|
def |
В |
|
def |
В |
|
. |
(5.6) |
4 |
3 4 4 2 4 |
h4 |
|
|||||||
|
|
|
|
|
||||||
Как видно, коэффициент отражения от барьера R пропорционален температуре решетки T и не зависит от энергии электрона.
69
Вероятность рассеяния электрона на барьере равна коэффициенту отражения: W R . Если бы электрон с вероятностью W 1 отразился от единственного деформационного барьера, тогда его длина свободного пробега равнялась точно
4 . Таким образом, вероятность рассеяния можно выразить
как отношение периода деформационного потенциала к длине свободного
пробега электрона: W |
|
. Тогда из (5.6) получим, что длина свободно- |
|||
4l |
|||||
|
|
|
|
||
|
|
|
св.пр |
|
|
го пробега |
l |
T 1 |
и, как и коэффициент отражения R , не зависит от |
||
|
св.пр |
|
|
|
|
энергии электрона.
Длина свободного пробега электрона есть произведение средней скорости его теплового движения vт 3kВT
m * и времени свободного пробега: lсв.пр vт . Тогда окончательно получим для постоянной времени рассея-
ния |
электронов на деформационном |
акустическом потенциале решетки: |
||
ак |
|
lсв.пр |
T 3/2 и, соответственно, |
ак T 3/2 . |
|
||||
|
|
vт |
|
|
5.4. Проявление зонной структуры полупроводника и типа химической связи в решеточном рассеянии
Междолинное рассеяние электронов. Этот динамический механизм
рассеяния с участием фононов может наблюдаться в непрямозонных полупроводниках, зона проводимости которых содержит несколько эквивалент-
ных долин. В равновесных условиях электроны в долинах имеют одинаковую энергию и поэтому распределены поровну. Однако всякое воздействие (поле, давление, градиент температуры) приводит к перезаселению долин за счет
междолинного рассеяния.
В кремнии и других алмазоподобных полупроводниках n-типа различа-
ют два типа междолинного рассеяния:
– g-рассеяние, происходящее между данной долиной и долиной на противоположной стороне той же оси, например между долинами в направлени-
ях 100 и 100 ;
– f-рассеяние, когда электрон переходит из данной долины в одну из остальных(перпендикулярных) долин, напримериздолины 100 вдолину 010 .
Носители заряда в долинах отличаются направлением волнового вектора, и переход электрона из одной долины в k-пространстве в другую сопровожда-
70
