В 1-й лекции мы рассматривали систему невзаимодействующих электронов, взаимодействие которых с атомами/ионами кристалла сводилось к столкновениям со средней частотой
1/ . Электроны не ощущали периодического потенциала решетки. Такую систему мы называли газом свободных электронов.
В этой лекции мы начнем рассмотрение электронов в идеальном кристалле, где атомы расположены в строгом порядке, определяемом решеткой Браве. Следовательно, потенциал взаимодействия также обладает периодичностью решетки Браве:
U(r + R) = U(r),
где R — любой из векторов решетки Браве.
Уравнение Шредингера для электрона в периодической решетке
Н = (-( 2 /2m) 2 + U(r)) = E (r),
где потенциал обладает свойством (2.1).
Одним из основных следствий периодичности потенциала является теорема Блоха.
Собственные волновые функции одноэлектронного гамильтониана
Н = -( 2 /2m) 2 + U(r),
где U(r + R) = U(r), при всех R , принадлежащих решетке Браве, могут быть выбраны в форме плоской волны, умноженной на функцию с периодичностью решетки Браве, т.е.
nk(r) = exp(ikr) unk(r),
unk(r + R) = unk(r)
для всех R, принадлежащих решетке Браве. Здесь n — номер зоны, появление котoрого связано с тем, что для данного k имеется множество решений.
В иной записи теорема Блоха имеет вид
(r + R) = exp(ikR)(r),
что вытекает из (2.3) и (2.4).
Обобщенное граничное условие Борна-Кармана для периодического потенциала.
Вместо макроскопического “ящика” с размерами L для периодических граничных условий Борна-Кармана в теории Зоммерфельда, выбираем “ящик” соразмерный элементарной ячейке соответствующей решетке Браве. Такое граничное условие является естественным обобщением граничных условий Борна-Кармана для периодического потенциала. Запишем его в виде
nk(r + Niai) = nk(r), i = 1, 2, 3
Здесь аi — тройка основных векторов, а все Ni — целые числа достигающие величин порядка N 1/3 , где N = N1N2N3 — полное число элементарных ячеек в кристалле. Как и в теории Зоммерфельда, мы предполагаем, что объемные свойства кристалла не зависят от выбора граничных условий, которые поэтому могут быть выбраны из соображений удобства вычислений.
Применяя к граничным условиям теорему Блоха, находим
nk(r + Niai) = exp(iNikai)nk(r), i = 1, 2, 3
Соотношение (2.7) для произвольных Ni может выполняться только при условии
exp(iNikai) = 1 , где i =1, 2, 3, т.е. Nikai = 2mi,
где mi -целые числа.
Представим k в виде разложения по базису векторов bj обратной решетки (см. Д-2)
k = j=1,2,3 xjbj,
где хj — действительные числа, bjai = 2ij .
Тогда из (2.8) и (2.9) получим xj = mj/Nj и разрешенный блоховский волновой вектор должен иметь следующий общий вид
k = j=1,2,3 (mj /Nj )bj.
Из (2.10) следует, что объем k в k -пространстве, приходящийся на одно разрешенное значение k, равен объему параллелепипеда с ребрами bj/Nj:
k = b1 /N1 [(b2 /N2 ) x (b3 /N3 )] = (1/N) b1 (b2 x b3).
Cледовательно, N = b1 (b2 x b3)/k есть число разрешенных волновых векторов в 1-й элементарной ячейке обратной решетки, т.е. в 1-й зоне Бриллюэна.
Поскольку b1 (b2 x b3) есть объем элементарной ячейки обратной решетки, ф-ла (2.11) означает, что число разрешенных волновых векторов, содержащихся в одной элементарной ячейке обратной решетки, равно числу ячеек в кристалле.
Объем элементарной ячейки обратной решетки равен (2) 3 /v, где v = V/N — объем элементарной ячейки прямой решетки. Следовательно, формулу (2.11) можно записать в виде:
k = (2) 3 /V.
Это совпадает с результатом, полученным в теории Зоммерфельда.
Доказательство теоремы Блоха
Любую функцию, подчиняющуюся граничным условиям Борна -Кармана, можно разложить по базису плоских волн, удовлетворяющим этому условию, и поэтому имеющих волновые векторы типа
q = j=1 3 (mj/Nj)bj:
(r) = qcq exp(iqr)
Поскольку потенциал U(r) имеет периодичность решетки, в его разложении по плоским волнам будут входить только плоские волны с периодичностью решетки, поэтому их волн. векторы являются векторами обратной решетки.
U(r) = G UG exp(iGr)
Cумирование по всем векторам G обратной решетки.
Коэффициенты Фурье UG связаны с U(r) соотношением
UG = 1/v(по ячейке)dr exp(iGr) U(r)
Аддитивную постоянную выберем так, что
U0 = 1/v(по ячейке) dr U(r) = 0
Поскольку U(r) действителен, то U*G = U-G .Пусть кристалл обладает центром инверсии. Тогда при выборе начала координат в центре инверсии
Подставляя в уравнение Шредингера. ( — 2 / 2m 2 + U(r) = E), получаем
q 2 /2mq 2 cq exp(iqr) + (G UG exp(iGr)) (qcq exp(iqr)) = Eqcq exp(iqr)
q exp(iqr) [( 2 /2mq 2 — E)cq + GUGcq-G‘] = 0
<Пояснение: G,qUGcq exp(i(G-q)r)) = G,q‘UGcq‘-G exp(iq‘ r) = (q‘» width=»21″ height=»17″>q) = G,qUGcq-G exp(iq r)>
Поскольку плоские волны, удовлетворяющие граничному условию Борна-Кармана, образуют ортогональный набор, коэффициенты при каждом слагаемом в (2.17) должны быть равны нулю, поэтому для всех разрешенных волновых векторов q получаем
( 2 /2mq 2 — E)cq + GUGcq-G‘ = 0.
Уравнение (2.18) можно переписать в другом виде, представив вектор q в форме q = k — G, где вектор обратной решетки G выбран так, чтобы вектор k лежал в первой зоне Бриллюэна, и заменив G‘ » width=»21″ height=»17″>G‘ — G. А именно,
[ 2 /2m (k-G) 2 — E]ck-G + G‘UG‘-Gck-G‘ = 0.
Уравнения (2.18) и (2.19) — это новая форма записи уравнения Шредингера в импульсном пространстве. Простота записи этих уравнений основывается на том, что, согласно (2.14), коэффициенты Uk отличны от нуля только если k совпадает с одним из векторов обратной решетки.
Для данного разрешенного k из 1-й зоны Бриллюэна имеется N независимых решений ck-G , удовлетворяющих ур-ю (2.19) и отличающихся волновыми векторами на один из вектров обратной решетки. Решив эту систему уравнений, можем теперь, с учетом q = k — G, записать разложение волновой функции (2.13) в виде
(r) = Gck-G exp(i(k-G)r)
k(r) = exp(ikr)Gck-G exp(-iGr).
Т.е. видно, что функция имеет блоховскую форму (2.3), где функция u(r)
u(r) = G ck-G exp(-iGr),
что и требовалось доказать.
Зоны Бриллюэна и энергетические зоны
а) |
Как следует из теоремы Блоха, изменение вектора k на вектор обратной решетки не приводит к изменению состояния, чем он отличается от аналогичной величины в модели свободных электронов. Поэтому вектор k называют не волновым, а квазиволновым вектором. Поскольку вектор k пределен с точностью до векторов обратной решетки, то обычно за “начальный” вектор принимают вектор находящийся ближе других к началу координат k-пространства. Область —bi/2
где m-число вырождений уровня с данным k. Чтобы найти основные поправки к энергиям уровней достаточно решить систему уравнений для m-вырожденных уровней
(E-E 0 k-Gi) ck-Gi = UGj-Gi ck-Gj ,
где i = 1. m
Если же уровни далеки от вырождения, т.е. |Ek — E 0 k-G| >> U, то воздействие потенциала будет пропорционально квадрату потециала |Ek — E 0 k-G|
U 2 , т.е. мало (см. АМ, Гл.9)
Энергетические уровни вблизи одной из брэгговских плоскостей. Энергетическая щель.
Рассмотрим простейший пример, когда два уровня свободных электронов расположены близко друг к другу (по сравнению с U), но далеко от остальных уровней (также по сравнению с U). Т.е. при данном k имеются векторы обратной решетки G1 и G2, для которых |E 0 k-G1 — E 0 k-G2|
U, тогда как |E 0 k-G1,2 — E 0 k-G| >> U, для всех G G1,2. Уравнения (2.29) в этом случае сводятся к 2-м:
(E-E 0 k-G1) ck-G1 = UG2-G1 ck-G2,
(E-E 0 k-G2) ck-G2 = UG1-G2 ck-G1.
q = k — G1 и G = G2 — G1
Тогда (2.30) запишем в виде
(E-E 0 q) cq = UG cq-G
(E-E 0 q-G) cq-G = U-G cq = UG cq
Справедливы следующие соотношения:
E 0 qE 0 q-G, |E 0 q — E 0 q-G‘ | >> U при G‘ G, 0
Энергии E 0 q E 0 q-G лишь при условии |q| = |q — G|. Это означает, что точка q должна лежать на плоскости, которая перпендикулярна отрезку, соединяющему начальную точку в k-пространстве с точкой G обратной решетки, и делит этот отрезок пополам. Эта плоскость называется брэгговской по аналогии с рассеянием рентгеновских лучей. Одновременно, из условия (2.33) |E 0 q — E 0 q-G‘ | >> U при G‘ G, 0 следует, что вектор q должен лежать далеко от места пересечения 2-х и более плоскостей. Следовательно, случай 2-х почти вырожденных уровней соответствует электрону, волновой вектор которого почти точно удовлетворяет условию однократного брэгговского рассеяния. Соответственно, общий случай большого числа вырожденных уровней применим к исследованию такого уровня свободного электрона, волновой вектор которого почти точно удовлетворяет условию одновременного существования многих брэгговских отражений. Поскольку слабый периодический потенциал сильнее всего влияет на почти вырожденные уровни, можно сделать вывод, что главное воздействие слабый периодический потенциал оказывает лишь на те уровни свободного электрона, волновые векторы которых близки к векторам, допускающим брэговское отражение.
Для 2-х уровневой почти вырожденной системы (одна брэгговская плоскость) уравнения (2.32) имеют решение при условии
из которого следует квадратное уравнение
(E-E 0 q) (E-E 0 q-G) = |UG| 2 .
E = (E 0 q + E 0 q-G) /2 + [((E 0 q — E 0 q-G)/2) 2 + |UG| 2 ] 1/2
Рис. 2.2. |
описывают главный результат воздействия периодического потенциала на два уровня E 0 q и E 0 q-G свободных электронов, когда точка q близка к брэгговской плоскости, определяемой вектором G. Графически решения показаны на рис. 2.2. Для точек, лежащих непосредственно на брэгговской плоскости: E 0 q = E 0 q-G, следовательно
E = E 0 q + |UG|.
Т.е. один уровень понижается, а другой повышается на величину |UG|. Исходя из (2.26) можно показать, что при E 0 q = E 0 q-G выполняется соотношение
dE/dq = 2 / m (q — G/2).
Т.е. если точки q лежат на брэгговской плоскости, то градиент Е параллелен этой плоскости. Поскольку градиент перпендикулярен поверхности на которой функция постоянна, изоэнергетические поверхности вблизи брэгговской плоскости перпендикулярны ей. Если точка q принадлежит одной из брэгговских плоскостей, нетрудно также определить вид волновых функций, отвечающих двум решениям E = E 0 q + |UG|. Из уравнения (2.32) следует, что если Е дается выражением (2.37), то два коэффициента cq и сq-G удовлетворяют соотношению
cq = +sign (UG)сq-G.
Поскольку первые два члена доминируют в разложении (2.20) по плоским волнам, то получаем:
Если UG > 0, то:
|(r)| 2
(cos(1/2Gr)) 2 , E = E 0 q + |UG|
|(r)| 2
(sin(1/2Gr)) 2 , E = E 0 q — |UG|,
(sin(1/2Gr)) 2 , E = E 0 q + |UG|
|(r)| 2
(cos(1/2Gr)) 2 , E = E 0 q — |UG| .
Волновую функцию вида | ψ (r)| 2
(cos (1/2Gr)) 2 называют «s-типа» по аналогии с атомными волновыми функциями ( не обращается в ноль в центре ионного остова), а вида | ψ (r)| 2
(sin(1/2Gr)) 2 называют «p-типа». В последнем случае плотность зарядов для малых расстояний пропорциональна квадрату расстояния от иона и обращается в нуль на ионе.
а | д |
б | |
в | е |
г | ж |
Рис. 2.3 а) Параболическая зависимость Е(к), б) Наложение парабол от соседних брэгговских плоскостей при UG = 0 — “пустая” решетка, в) Расщепление в точке G/2 пересечения двух парабол при, г) искажение первоначальной параболы при UG 0 для 2-х брэгговских плоскостей, д) то же что и г) для многих плоскостей. Такой способ представления называют схемой расширенных зон, е) то же что и д) в схеме приведенных зон, ж) то же в схеме повторяющихся зон (Из АМ) |
Энергетическая щель. Если UG = 0, то при пересечении волновым вектором брэгговской плоскости, благодаря вырожденности, электрон плавно переходит с одной ветки на другую (рис. 2.2, 2.3б). Если UG 0, энергия непрерывно меняется только в том случае, если электрон останется на кривой, отвечающей либо нижней, либо верхней ветви, как это видно из (2.37) и рис. 2.2, 2.3в. Чтобы перейти с одной ветви на другую при непрерывном изменении k, энергия должна изменится скачком по крайней мере на величину 2|UG|, которую называют энергетической щелью.
Энергетические зоны в одномерном случае
Рис. 2.4. Построение поверхности Ферми для кристалла Pb для сечений (100) и (110) з.Б., проходящих через точки ГUK а) и ГХW б), соответственно |
На рис.2.1б зонная структура ε (к) приведена в т.н. схеме приведенных зон. В этой схеме зонная структура приведена к 1-й зоне, пользуясь периодичностью зонной структуры. На рис. 2.3 проиллюстрированы способы изображения зонной структуры для 2-кратно вырожденного уровня в одномерном случае в схемах расширенных, приведенных и повторяющихся зон. Если изображать уровни с учетом всех брэгговских плоскостей, то мы получим описание в схеме расширенных зон (рис. 2.3д).
Если задавать уровни с помощью векторов в 1-й зоне Бриллюэна, т.е. сместив уровни на векторы обратной решетки (рис. 2.3 е), то мы получим описание в схеме приведенных зон. Схема повторяющихся зон демонстрирует периодичность описания в к-пространстве, но данная схема избыточна, поскольку каждый уровень повторяется многократно.
Изображение зонной структуры в трехмерном к-пространстве.
Для изображения зонной структуры в 3-мерном пространстве волновых векторов, обычно используют кривые зависимости вдоль определенных направлений k (рис.2.1а,б). Обычно изображение делается в схеме приведенных зон (рис.2.1б), поскольку в общем случае произвольного направления в k-пространстве структура зон не обладает периодичностью. Как видно, даже в приближении “пустой” решетки получаемые зависимости E(k) оказываются достаточно сложными. Большинство кривых вырождено. Из-за конечного, пусть слабого, взаимодействия с решеткой вырождение частично снимается. Описание требует использования аппарата теории групп.
Построение поверхности Ферми. В модели свободных электронов Зоммерфельда энергия Ферми дается (1.26) Е F = (hkF) 2 / 2m, где согласно (1.23)
kF = (3 2 N/ V) 1/3 = (3 2 /0) 1/3 ,
где 0 — объем, приходящийся на один атом. В этом случае поверхность Ферми — сфера. В блоховском кристалле — поверхность Ферми отличается от сферы. Многие свойства кристалла зависят от соотношения размеров Ферми-поверхности и 1-й зоны Бриллюэна. Поверхность Ферми может быть полностью вмещаться в ней не касаясь стенок, касаться и переходить в зоны более высокого порядка.
Рассмотрим построение зонной структуры на примере кристалла Pb c гцк структурой (следуя Анималу, Гл.5, § 2, ААК, Гл.4 § 6). Вычисляем радиус Ферми-сферы по (2.31), содержащей z = 4 электрона на атом в этом случае 0 = a 3 /4 и kF = (12 2 4/a 3 ) 1/3 = 1.241 (2/a). Как видно, полученное значение kF больше радиуса сферы, вписанной в соответствующую данному случаю зону Бриллюэна с радиусом kBril. = (2/a). Это приводит к тому, что поверхность Ферми имеет сложную конфигурацию, проникая во 2-ю и 3-ю зоны Бриллюэна. Построив сферы радиусом kBril. с центром в каждом из соседних узлов обратной решетки, как это сделано на рис. 2.4 для двух сечений (100) (а) (110) (б) з.Б., можно представить заполнение 2-й з.Б. в схеме приведенных зон как перекрытие з.Б. соседней сферой Ферми, 3-й з.Б. двумя соседними сферами Ферми (см. детали в Анималу, Гл.5, § 2 и ААК § 6).
Видео:Урок 455. Уравнение ШрёдингераСкачать
ЗОННАЯ ТЕОРИЯ ТВЕРДЫХ ТЕЛ
2.1. Движение электронов в периодическом поле кристалла.
Уравнение Шредингера для кристалла
В первой главе обсуждалось квантово-механическое описание свободных микрочастиц или частиц, находящихся во внешнем силовом поле. Однако основные успехи квантовой механики связаны с изучением систем взаимодействующих микрочастиц (электронов, ядер, атомов, молекул), из которых состоит вещество. В этой главе мы применим квантовую механику к описанию поведения электронов в твердых кристаллических телах, рассматривая кристалл как систему микрочастиц.
В общем случае эта задача требует решения уравнения Шредингера для системы частиц (электронов и ядер), образующих кристалл. В этом уравнении необходимо учесть кинетическую энергию всех электронов и ядер, потенциальную энергию взаимодействия электронов между собой, ядер между собой, электронов с ядрами. Понятно, что в общем виде решение такого уравнения не представляется возможным, поскольку оно содержит порядка 10 22 переменных. Поэтому задачи, связанные с поведением электронов в кристалле, решаются при некоторых упрощающих допущениях (приближениях), правомерность которых определяется конкретными свойствами кристалла. Рассмотрим основные из этих допущений.
Адиабатическое приближение. В этом приближении предполагается, что электроны движутся в поле неподвижных ядер. Под ядрами здесь подразумевают собственно ядра атомов со всеми электронам, исключая валентные. Правомерность этого допущения определяется тем, что скорости электронов приблизительно на два порядка больше, чем скорости ядер, поэтому для любой, даже неравновесной конфигурации ядер всегда будет успевать устанавливаться соответствующее ей электронное равновесие. В этом представлении исключается обмен энергией между электронной и ядерной системами, поэтому это приближение называется адиабатическим. Естественно, что в адиабатическом приближении нельзя рассматривать такие явления, как диффузия, ионная проводимость и др., связанные с движением атомов или ионов.
Одноэлектронное приближение. В этом приближении вместо взаимодействия данного электрона с остальными электронами и ядрами по отдельности рассматривают его движение в некотором результирующем усредненном поле остальных электронов и ядер. Такое поле называют самосогласованным. В одноэлектронном приближении, таким образом, задача сводится к независимому описанию каждого электрона в среднем внешнем поле с потенциальной энергией U(r). Вид функции U(r) определяется свойствами симметрии кристалла. Основное свойство самосогласованного поля заключается в том, что оно имеет тот же период, что и поле ядер.
Таким образом, адиабатическое и одноэлектронное приближение приводит к задаче движения электрона в некотором периодическом потенциальном поле, имеющем период, равный постоянной решетки кристалла. Уравнение Шредингера в этом случае будет иметь вид
. (2.1)
Следующие два допущения связаны с невозможностью точно определить вид функции U(r). Поэтому обычно при описании свойств электронов в кристалле рассматривают два предельных случая взаимодействия электронов с решеткой.
Приближение слабой связи. В этом приближении электроны в кристалле рассматривают как почти свободные частицы, на движение которых оказывает слабое возмущение поле кристаллической решетки. Данное допущение применимо, когда потенциальная энергия взаимодействия электрона с решеткой много меньше его кинетической энергии. Такой подход, который иногда называют «приближением почти свободных электронов«, позволяет получить решение некоторых задач, связанных с поведением валентных электронов в металлах.
В полупроводниках более приемлемым для анализа их физических свойств является приближение сильной связи. В этом приближении состояние электрона в кристалле мало отличается от его состояния в изолированном атоме. Приближение сильной связи применимо, когда потенциальная энергия электрона значительно больше его кинетической энергии.
Характерным для обоих приближений слабой и сильной связи является то, что оба они приводят к фундаментальному свойству энергетического распределения электронов в кристалле — возникновению разрешенных и запрещенных энергетических зон.
2.2. Энергетические зоны в приближении сильной связи
Несмотря на то, что метод сильной связи применим для электронов глубоких энергетических уровней, он хорошо иллюстрирует общие закономерности образования энергетических зон при сближении изолированных атомов и образования из них кристаллической решетки. Рассмотрим качественно картину возникновения энергетических зон на примере образования кристаллической решетки из изолированных атомов натрия. Электронная структура Na 11 (1s 2 2s 2 2p 6 3s): всего в атоме 11 электронов, по два электрона на 1s и 2s уровнях, 6 электронов на уровне 2р, последний заполненный уровень в атоме натрия — 3s, на котором находится один валентный электрон. Поскольку в приближении сильной связи предполагается, что состояние электрона в кристалле незначительно отличается от его состояния в изолированном атоме, будем в оценке влияния на это состояние кристаллического поля соседних атомов исходить из энергетической структуры изолированного атома. На рис. 2.1,а показаны схематически энергетические уровни и распределение электронов на них для атомов натрия, находящихся на достаточно большом расстоянии друг от друга так, что потенциальные кривые электронов не перекрываются (взаимодействие между атомами пренебрежимо мало). Состояния электронов в этом случае описываются волновыми функциями изолированного атома, разрешенные уровни энергии дискретны и определяются квантовыми числами n, l, m — главным, орбитальным, магнитным соответственно. На каждом невырожденном по энергии уровне могут находиться с учетом спина по два электрона, а на каждом вырожденном по орбитальному квантовому числу уровне 2(2l +1) электронов.
Рис. 2.1. Изменение состояния электронов при сближении атомов
Сблизим теперь эти атомы на расстояние, равное параметру кристаллической решетки натрия (рис. 2.1,б). Взаимодействие с соседними атомами будет оказывать влияние на первоначальные атомные энергетические уровни. В приближении сильной связи предполагается, что потенциальная энергия электрона в кристалле U(r) может быть представлена суммой
, (2.2)
где Ua — потенциальная энергия электрона в изолированном атоме; D U (r) — поправка, учитывающая влияние соседних атомов. Предполагается, что соседние атомы оказывают слабое возмущение на Ua( D U (r) D U (r) приводит к уравнению Шредингера для изолированного атома.
Поскольку в кристалле каждый уровень изолированного атома повторяется N раз, то он становится N-кратно вырожденным. Известно, что электрическое поле снимает вырождение и каждый уровень изолированного атома расщепляется на N близко расположенных (по значениям энергии) энергетических уровней. Здесь имеется аналогия со связанными осцилляторами. Если мы имеем два не связанных между собой каким-либо взаимодействием совершенно одинаковых осциллятора (математические маятники, электрические колебательные контуры и др.), то частоты их собственных колебаний совпадают. Взаимодействие между осцилляторами приводит к расщеплению одной частоты на две близкие частоты (при условии, что энергия взаимодействия между осцилляторами много меньше энергии собственных колебаний). Для N связанных между собой осцилляторов получим полосу из N близко расположенных частот. Аналогичный результат получается для системы взаимодействующих атомов. Число энергетических уровней, на которые расщепляется каждый энергетический уровень изолированного атома, равно числу атомов в кристалле. Величина расщепления тем больше, чем сильнее взаимодействие между атомами, т.е. чем меньше расстояние между ними. На рис. 2.2 показано схематически расщепление двух энергетических уровней атома под воздействием полей соседних атомов. Схема приведена для восьми атомов.
Решение уравнения Шредингера в приближении сильной связи приводит к следующему выражению для энергии электрона в периодическом поле трехмерной кубической решетки
, (2.3)
здесь C — некоторая постоянная величина, которая может принимать положительные и отрицательные значения; А — обменный интеграл, зависящий от перекрытия волновых функций атомов; k x , k y , k z — компоненты волнового вектора электрона; а — параметр решетки кристалла.
Рис.2.2. Расщепление энергетических уровней атома
Экстремальные значения энергии электрона Е имеют место при coskia = ± 1 (i = x, y, z) и определяют ширину энергетической зоны, образованной расщепленным уровнем изолированного атома. Для простой кубической решетки ширина энергетической зоны D E = 12A . Ширина энергетической зоны для более высоких уровней больше, т.к. для этих состояний электронов сильнее перекрываются волновые функции электронов и, следовательно, больше обменный интеграл А. Середина зоны сдвинута относительно положения энергетического уровня изолированного атома на величину С. Направление смещения зависит от знака С. Энергетические зоны в общем случае разделены интервалами энергии D E g , называемыми запрещенными зонами. Иногда энергетические зоны могут перекрываться.
В реальных кристаллах размером приблизительно 1 см 3 содержится
10 22 атомов. Ширина энергетической зоны обычно
1 эВ. В этом случае расстояние между уровнями в зоне составляет
10 -22 эВ. Следовательно, спектр электронов в пределах энергетической зоны можно считать практически непрерывным.
2.3. Общие свойства волновой функции электрона в периодическом потенциале. Теорема Блоха
Для точного решения в одноэлектронном приближении задачи о движении электрона в кристалле необходимо решить уравнение Шредингера вида (2.1), где потенциал U(r) имеет периодичность кристаллической решетки, т.е.
, (2.3)
здесь R — любой вектор прямой кристаллической решетки.
Необходимость решения квантово-механической задачи связана с тем, что длина волны де Бройля электрона по порядку величины совпадает с периодом потенциала U (
10 -8 cм). Можно получить некоторые общие свойства волновой функции электрона в кристалле, используя только свойство периодичности потенциала кристаллического поля, не решая уравнения Шредингера. Мы будем рассматривать здесь идеализированный случай гипотетического кристалла с абсолютно идеальной периодичностью потенциала. Типичный вид потенциала вдоль линии, соединяющей цепочку атомов (одномерный случай) мы получили ранее, анализируя качественно влияние взаимодействия атомов на спектр электронов при сближении изолированных атомов (рис. 2.1,б). Точное определение функции U(r) является очень сложной задачей.
Фундаментальные свойства волновой функции стационарного состояния определяются теоремой Блоха: собственные функции стационарного волнового уравнения с периодическим потенциалом имеют вид произведения функции плоской волны на функцию с периодичностью потенциала:
. (2.4)
Индекс k у волновой функции указывает на то, что эта функция зависит от волнового числа. Появление индекса n связано с тем, что при фиксированных значениях k волновая функция не одинакова для электронов различных зон, образовавшихся из атомных уровней, n часто называют номером зоны. Множитель un,k(r) называют блоховским множителем. Он учитывает влияние кристаллического поля и отражает тот факт, что вероятность нахождения электрона в той или иной области кристалла повторяется от ячейки к ячейке.
Схематическое изображение электронных волновых функций, представленных в теореме Блоха, показано для одномерного случая на рис.2.3. Вверху (рис. 2.3,а) представлен потенциал U(x) вдоль цепочки атомов. Ниже (рис. 2.3,б) приведен пример собственной функции (ее действительной части). Эта функция равна произведению блоховского множителя u(x), имеющего периодичность решетки (рис. 2.3,в) и волновой функции свободного электрона в виде плоской волны (рис. 2.3,г), длина которой определяется волновым числом k. Представление волновой функции в виде (2.4) может быть сделано различными способами. Покажем это для одномерного случая. Одномерная волновая функция по теореме Блоха может быть записана в виде
. (2.5)
Домножим и разделим правую часть равенства (2.5) на функцию , где
а — параметр решетки. Тогда получим
. (2.6)
Рис. 2.3. Схематическое изображение электронных волновых функций в кристалле
В квадратных скобках формулы (2.6) стоит функция, удовлетворяющая требованиям теоремы Блоха: она является периодической с периодом а, т.к. равна произведению двух периодических функций с тем же периодом. Функция описывает плоскую волну, но с другим волновым вектором, отличающимся на величину . Таким образом, одно и то же стационарное состояние электрона в кристаллическом периодическом поле может быть описано как волновой функцией с волновым числом k , так и волновой функцией с волновым числом и другим блоховским множителем. Аналогичные результаты получатся, если k изменить на величину , где n — любое целое число.
Для одномерной цепочки атомов величина совпадает с размером первой зоны Бриллюэна в обратном пространстве. Если ограничиться рассмотрением волновых чисел в пределах первой зоны Бриллюэна, т.е. в интервале от до , то этот набор k исчерпает все физически различные значения волновых чисел в кристалле.
2.4. Модель Кронига-Пенни
Теорема Блоха позволяет аналитически решить задачу об электроне в периодическом поле кристаллической решетки в приближении слабой связи при некоторых упрощающих предположениях. Основная трудность в решении уравнения (2.1) связана с невозможностью точно записать вид функции U(r). Поэтому часто периодическую зависимость функции U(r) заменяют более простой функцией с точно таким же периодом. В модели Кронига-Пенни ограничиваются рассмотрением одномерной задачи, в которой периодический потенциал заменяется цепочкой прямоугольных потенциальных ям (рис. 2.4). Ширина каждой ямы а, они отделены друг от друга прямоугольными потенциальными барьерами высотой U 0 и шириной b. Период повторения ям с = а + b.
Рис.2.4. Изменение потенциальной энергии электрона:
а — в реальном кристалле; б — в модели Кронига-Пенни
Стационарное уравнение Шредингера будет иметь в этом случае вид
. (2.7)
Начало системы координат (точку х = 0) выберем так, чтобы она совпадала с левым краем потенциальной ямы, как это показано на рис. 2.4,б. Tогда потенциальная функция
. (2.8)
В соответствии с теоремой Блоха волновая функция электрона y (x) может быть представлена в виде
. (2.9)
Индексы n и k упущены для простоты записи. Функция u(x) (блоховский множитель) имеет период c
Подставляя (2.9) в уравнение (2.7), получим дифференциальное уравнение для блоховского множителя
(2.10a)
для электронов, находящихся внутри потенциальных ям, и
(2.10б)
для электронов, находящихся вне потенциальных ям. В этих уравнениях Ek — кинетическая энергия электрона
.
Общее решение уравнения (2.10а) для электронов внутри потенциальных ям может быть записано в виде
, (2.11а)
где a — некоторый параметр, который может быть найден подстановкой решения в виде (2.11а) в исходное уравнение (2.10а). Эта подстановка приводит к следующему значению a :
.
В области вне потенциальных ям при условии, что высота потенциального барьера U0 выше полной энергии электрона Е, решение уравнения (2.10б) имеет вид
, (2.11б)
.
Постоянные A, B, C и D в формулах (2.11а) и (2.11б) находятся как обычно из граничных условий. Граничные условия требуют, чтобы функция u(x) и ее первая производная в местах скачков потенциала, т. е. на стенках потенциальных ям, были непрерывны. Эти требования приводят к следующей системе уравнений:
(2.12)
Система уравнений (2.12) после подстановки в нее функций и , согласно равенствам (2.10а) и (2.10б), преобразуется в систему линейных однородных алгебраических уравнений, в которых неизвестными являются коэффициенты A, B, C и D. Определитель этой системы будет равен нулю (только при этом условии система линейных однородных уравнений имеет отличные от нуля решения), если выполняется следующее равенство:
. (2.13)
Выражение (2.13) можно значительно упростить, если допустить, что ширина барьера стремится к нулю , а его высота — к бесконечности , но таким образом, чтобы произведение U0b оставалось постоянным . При этих условиях выражение (2.13) преобразуется к виду:
, (2.14)
.
Поскольку a — параметр, определяемый энергией Е электрона, а k — волновой вектор электрона, то выражение (2.14) представляет зависимость E(k), т. е. дисперсионное соотношение для электрона в кристаллической решетке. Это дисперсионное соотношение можно записать в явном виде, решив уравнение (2.14) относительно a при фиксированном значении параметра p.
2.5. Энергетические зоны в модели Кронига-Пенни
Найдем в явном виде дисперсионное соотношение для электрона в периодическом кристаллическом поле. Исследуя выражение (2.14) находим, что волновое число k может быть вещественным только при условии, что значения левой части этого равенства находятся в интервале от -1 до +1. Зависимость левой части уравнения (2.14) от a для параметра p = 2 приведен на рис. 2.5. Заштрихованные участки соответствуют запрещенным значениям параметра a и, следовательно, энергии электрона в кристалле. Этот результат получен только на основании теоремы Блоха, условием применимости которой является единственное требование периодичности потенциала в стационарном уравнении Шредингера для электрона в кристалле. Таким образом, наличие периодического потенциала приводит к появлению для энергии электрона таких интервалов, для которых нет волнового решения, соответствующего вещественным значениям волнового числа электрона. Результатом этого является чередование разрешенных и запрещенных зон энергии для электрона в кристалле.
Рис. 2.5. Зависимость от параметра a левой части уравнения (2.14)
Рис. 2.6. Зависимость энергии электрона от волнового числа для p =2 и p =0 (штриховая линия)
На рис. 2.6 приведено дисперсионное соотношение для энергии электрона в кристалле. Видно, что зависимость E(k) претерпевает разрывы в точках, где и т. д.
Если параметр p = 0, согласно равенству (2.14) и
.
Последнее равенство соответствует дисперсионному соотношению для свободного электрона. На рис. 2.6 это дисперсионное соотношение изображено штриховой линией.
Поскольку, как подчеркивалось выше, все физически различимые значения волнового числа лежат в пределах первой зоны Бриллюэна, которая в одномерном случае ограничена интервалом значений волнового числа от до , целесообразно перейти от представления расширенных зон Бриллюэна (рис. 2.6) к представлению приведенных зон Бриллюэна (рис. 2.7). Волновые функции, соответствующие вещественным k, могут быть построены только для заштрихованных областей энергии электрона. Эти области представляют собой разрешенные энергетические зоны, которые отделены друг от друга зонами (щелями) запрещенных энергий.
Рис.2.7. Энергия электрона как функция волнового числа в схеме приведенных зон Бриллюэна
Предел P ® ¥ дает дискретный ряд уровней
которые совпадают с полученными в первой главе результатами для частицы в одномерной прямоугольной потенциальной яме (см. уравнение (1.34)). Энергия электронов в периодическом поле кристалла претерпевает разрыв на границах зон Бриллюэна, для которых . Физическая природа разрывов связана с
отражением электронных волн от атомных плоскостей кристаллической решетки. Действительно, с учетом того, что , условие, при котором происходит нарушение непрерывности функции E(k), может быть записано в виде , что совпадает с условием Вульфа-Брэгга при угле падения волн 90 0 .
2.6. Заполнение энергетических зон электронами.
Металлы, диэлектрики и полупроводники
Твердые тела делятся на металлы, диэлектрики и полупроводники прежде всего по величине удельной электропроводности. Для типичных металлов эта величина составляет 10 8 . 10 6 (Ом • м) -1 . В диэлектриках удельная электропроводность ничтожно мала: s -8 (Ом • м) -1 . Для хороших диэлектриков удельная электропроводность достигает величины 10 -11 (Ом • м) -1 . Твердые тела с промежуточной электропроводностью относят к полупроводникам. Оказывается, что столь большие различия в электрических свойствах твердых тел связаны со структурой и степенью заполнения электронами энергетических зон в этих телах.
Несмотря на то, что энергетические зоны квазинепрерывны, они состоят пусть из очень большого, но конечного числа энергетических уровней. Число этих уровней определяется числом атомов N, объединенных в кристалл, и орбитальным квантовым числом l:
(2.15)
В каждой энергетической зоне могут располагаться в соответствии принципом Паули не более 2(2l + 1) электронов — по два с противоположными спинами на каждом уровне. Число электронов в кристалле также конечно и зависит как от числа атомов N, так и от количества электронов в атоме. Поскольку электроны стремятся занять энергетические уровни с наинизшей энергией, то в кристалле нижние энергетические зоны оказываются полностью заполненными, а самые верхние заполнены либо частично, либо совершенно свободны.
Частично заполненная зона образуется, например, у кристалла натрия. Этот элемент имеет полностью заполненные 1s-, 2s- и 2p-уровни, на которых располагается в общей сложности 10 электронов. В кристалле Na соответствующие 1s-, 2s- и 2p-зоны также будут полностью заполнены. Одиннадцатый валентный электрон в атоме Na располагается на 3s-уровне, на котором могут располагаться 2 электрона. Следовательно, 3s-зона кристаллического натрия будет заполнена лишь наполовину. Зонная структура Na приведена на рис. 2.8,a. Заполненные электронами зоны и часть 3s-зоны заштрихованы. Eg — ширина запрещенной зоны.
Часто частично заполненная зона образуется в результате перекрытия полностью заполненной зоны со следующей совершенно свободной. Пример такой зонной структуры приведен на рис. 2.8,б для бериллия, у которого перекрываются заполненная 2s- и свободная 2p-зоны.
Большую группу составляют кристаллы, у которых над целиком заполненным зонами располагаются совершенно пустые зоны, причем ширина запрещенной зоны варьируется у них от нескольких десятков электронвольт до единиц электронвольт. Типичные примеры этой группы кристаллов показаны на рис. 2.8, в, г. Это углерод в модификации алмаза и кремний.
Структура энергетических зон кристалла оказывает решающее влияние на величину его электропроводности. Поскольку электрический ток есть направленное движение зарядов (в металлах — электронов), то возникновение электрического тока связано с увеличением импульса электронов вдоль направления действующей на него силы. Вместе с импульсом электрона меняется его волновой вектор. Поскольку энергия и волновой вектор электрона — две взаимосвязанные величины, связь между которыми определяется дисперсионным соотношением, то увеличение волнового числа должно обязательно сопровождаться увеличением энергии электрона. Нетрудно оценить, каково увеличение энергии электрона за счет его ускорения в электрическом поле, вызывающим электрический ток в проводниках. Если величина напряженности электрического поля равна 10 4 В/м, то на расстоянии, равном средней длине свободного пробега электрона в кристалле, а она обычно составляет
10 -8 м, электрон приобретает энергию приблизительно 10 -4 эВ. Понятно, что эти значения энергии позволяют электрону переходить с уровня на уровень только внутри одной энергетической зоны. Для перехода между зонами необходима энергия больше ширины запрещенной зоны Eg, которая, как указывалось выше, составляет 0.1 . 10 эВ.
Рис.2.8. Заполнение энергетических зон электронами
Эти рассуждения приводят к выводу о том, что для появления у тел высокой проводимости необходимо, чтобы в их энергетическом спектре присутствовали зоны, заполненные частично. На свободные уровни этих зон могут переходить электроны, увеличившие свою энергию под действием внешнего электрического поля (рис. 2.9). Поэтому тела с частично заполненными энергетическими зонами являются проводниками. Частично заполненные зоны имеют все металлы.
Теперь рассмотрим кристаллы, верхняя энергетическая зона которых заполнена электронами полностью (рис. 2.8, в, г). Внешнее электрическое поле не в состоянии изменить характер движения электронов, т. к. оно не в состоянии поднять электроны в вышележащую свободную зону. Внутри же самой полностью заполненной зоны, не содержащей ни одного свободного уровня, оно может вызывать лишь перестановку электронов местами, что не нарушает симметрии их распределения по скоростям. Это не приводит к возникновению электрического тока в таких кристаллах.
Таким образом, твердые тела с полностью заполненными электронами энергетическими зонами являются непроводниками. По ширине запрещенной зоны непроводники делятся на диэлектрики и полупроводники.
К диэлектрикам относят тела, имеющие относительно широкую запрещенную зону. У типичных диэлектриков Eg > 3 эВ. Так, у алмаза Eg = 5,2 эВ; у нитрида бора Eg = 4,6 эВ; у Al2O3 Eg = 7 эВ.
У типичных полупроводников ширина запрещенной зоны менее 3 эВ. Например, у германия Eg = 0,66 эВ; у кремния Eg = 1,12 эВ; у антимонида индия Eg = 0,17 эВ.
Верхняя заполненная зона полупроводников и диэлектриков называется валентной зоной, следующая за ней свободная зона называется зоной проводимости. В металлах частично заполненную зону называют как валентной зоной, так и зоной проводимости.
2.7. Эффективная масса электрона в кристалле и ее физический смысл
Особенности движения электронов в кристалле обусловлены их взаимодействием с кристаллической решеткой. Оказывается, что движение отдельного электрона в кристалле можно описывать тем же уравнением, что и для свободной частицы, т.е. в виде второго закона Ньютона, в котором учитываются только внешние по отношению к кристаллу силы.
Рассмотрим движение электрона в кристалле под действием внешнего электрического поля. Внешнее электрическое поле приводит к увеличению скорости электрона и, следовательно, его энергии. Поскольку электрон в кристалле — это микрочастица, описываемая волновой функцией, то энергия электрона зависит от его волнового вектора. Зависимость между этими двумя характеристиками электрона в кристалле определяется дисперсионным соотношением, которое в свою очередь зависит от строения энергетических зон. Поэтому при расчете движения электрона в кристалле необходимо исходить из закона дисперсии.
Свободный электрон описывается монохроматической волной де Бройля и электрон в этом состоянии нигде не локализован. В кристалле же электрону необходимо сопоставить группу волн де Бройля с различными значениями частот w и волновых векторов k. Центр такой группы волн перемещается в пространстве с групповой скоростью
Эта групповая скорость соответствует скорости перемещения электрона в кристалле.
Задачу о движении электрона будем решать для одномерного случая. Увеличение энергии электрона dE под действием внешней силы F равно элементарной работе dA, которую совершает внешняя сила за бесконечно малый промежуток времени dt:
(2.16)
Учитывая, что для электрона как микрочастицы , имеем следующее выражение для групповой скорости
Подставляя полученное выражение для групповой скорости в формулу (2.16), получим
Распространяя этот результат на трехмерный случай, получим векторное равенство
(2.17)
Как видно из этого равенства, величина ћ k для электрона в кристалле изменяется со временем под действием внешней силы точно так же, как импульс частицы в классической механике Несмотря на это, ћ k нельзя отождествить с импульсом электрона в кристалле, поскольку компоненты вектора k определены с точностью до постоянных слагаемых вида (здесь a — параметр кристаллической решетки, ni=1, 2, 3, . ). Однако в пределах первой зоны Бриллюэна ћ k обладает всемисвойствами импульса. По этой причине величину ћ k называют квазиимпульсом электрона в кристалле.
Вычислим теперь ускорение a, приобретаемое электроном под действием внешней силы F. Ограничимся, как и в предыдущем случае, одномерной задачей. Тогда
При вычислении ускорения учтено, что энергия электрона является функцией времени . Учитывая, что , получим
(2.18)
Сравнивая выражение (2.18) со вторым законом Ньютона, видим, что электрон
в кристалле движется под действием внешней силы так, как двигался бы под действием той же силы свободный электрон, если бы он обладал массой
(2.19)
Величину m* называют эффективной массой электрона в кристалле.
Строго говоря, эффективная масса электрона никакого отношения к массе свободного электрона не имеет. Она является характеристикой системы электронов в кристалле в целом. Вводя понятие эффективной массы, мы реальному электрону в кристалле, связанному взаимодействиями с кристаллической решеткой и другими электронами, сопоставили некую новую свободную “микрочастицу”, обладающую лишь двумя физическими параметрами реального электрона — его зарядом и спином. Все остальные параметры — квазиимпульс, эффективная масса, кинетическая энергия и т.д. — определяются свойствами кристаллической решетки. Такую частицу часто называют квазиэлектроном, электроном-квазичастицей, носителем отрицательного заряда или носителем заряда n-типа, чтобы подчеркнуть ее отличие от реального электрона.
Особенности эффективной массы электрона связаны с видом дисперсионного соотношения электрона в кристалле (рис.2.10). Для электронов, располагающихся у дна энергетической зоны, дисперсионное соотношение можно приблизительно описать параболическим законом
Вторая производная , следовательно, эффективная масса положительная. Такие электроны ведут себя во внешнем электрическом поле подобно свободным электронам: они ускоряются под действием внешнего электрического поля. Отличие таких электронов от свободных состоит в том, что их эффективная масса может существенно отличаться от массы свободного электрона. Для многих металлов, в которых концентрация электронов в частично заполненной зоне мала и они располагаются вблизи ее дна, электроны проводимости ведут себя подобным образом. Если к тому же эти электроны слабо связаны с кристаллом, то их эффективная масса незначительно отличается от массы покоя реального электрона.
Для электронов, находящихся у вершины энергетической зоны (рис.2.10), дисперсионное соотношение можно приблизительно описать параболой вида
и эффективная масса является величиной отрицательной. Такое поведение эффективной массы электрона объясняется тем, что он при своем движении в кристалле обладает не только кинетической энергией поступательного движения Ек, но и потенциальной энергией его взаимодействия с кристаллической решеткой U. Поэтому часть работы A внешней силы может перейти в кинетическую энергию и изменить ее на величину D Eк , другая часть — в потенциальную D U :
Рис. 2.10. Закон дисперсии для электрона в кристалле
Рис. 2.11. Зависимость эффективной массы электрона от волнового числа
Если при движении электрона в потенциальную энергию переходит не только вся работа внешней силы, но и часть кинетической энергии, имевшейся у электрона ( D Eк 0 ), то его скорость будет уменьшаться. В этом случае электрон ведет себя как частица с отрицательной эффективной массой. В случае, когда вся работа внешней силы переходит в потенциальную энергию ( D Eк = 0 ), то приращения кинетической энергии и скорости не происходит. Электрон ведет себя как частица с бесконечно большой эффективной массой. Бесконечно большой эффективной массой обладает электрон в точках перегиба дисперсионной кривой, которые на рис. 2.10 обозначены штриховыми линиями. Схематически зависимость эффективной массы электрона от его волнового числа показана на рис. 2.11.
2.8. Собственные полупроводники. Понятие о дырках
Из структуры энергетических зон полупроводников следует, что при абсолютном нуле они не проводят электрического тока. Нагревание их приводит к тому, что часть электронов валентной зоны приобретает энергию, достаточную для их перехода в зону проводимости, в результате чего появляется заметная электропроводность. С увеличением температуры число электронов в зоне проводимости увеличивается и вместе с этим растет электропроводность полупроводника. Тепловое возбуждение электронов проводимости иллюстрирует рис. 2.12. Ес и Еv обозначают дно зоны проводимости и потолок валентной зоны соответственно. Кроме температуры, возбуждение электронов проводимости может происходить и под действием других факторов, способных сообщить электронам энергию, достаточную для перехода их в зону проводимости. Этими факторами могут быть световое облучение, ионизирующее излучение и др.
Рассмотренный выше механизм возникновения электропроводности полупроводниковых кристаллов, справедлив для абсолютно чистых материалов, не содержащих примесей, влияющих на электропроводность. Такие полупроводники называются собственными, а их электропроводность собственной электропроводностью. К собственным полупроводникам относятся кристаллы чистых химических элементов, таких как германий (Ge), кремний (Si), селен (Se), теллур (Te) и др., а также некоторые химические соединения: арсенид галлия (GaAs), арсенид индия (InAs), антимонид индия (InSb), карбид кремния (SiC) и многие другие.
В разделе 2.8 показано, что электроны, расположенные у по-толка энергетической зоны, об-ладают отрицательной эффектив-ной массой. Именно такие электроны, расположенные у вершины валентной зоны, переходят в зону проводимости и участвуют в электропроводности полупроводника. Каждому электрону, перешедшему в зону проводимости, в валентной зоне соответствует незанятое (вакантное) состояние, которое называют дырочным состоянием. Дырочные состояния изображены на рис. 2.12 светлыми кружками. Наличие вакансий в валентной зоне позволяет электронам этой зоны изменять свое энергетическое состояние под действием внешнего электрического поля. Рассмотрим подробнее этот процесс на примере кристалла, в котором имеется одно вакантное состояние. В отсутствие электрического поля это состояние будет находиться в вершине зоны, т.к. электроны стремятся расположиться на уровнях с наименьшей энергией (рис. 2.13,а). Занятые электронами состояния изображены на рис. 2.13 точками и расположены на дисперсионной кривой, описывающей зависимость энергии электрона от компоненты волнового вектора k x . У вершины энергетической зоны эта кривая приблизительно описывается параболой. Если к полупроводнику приложить внешнее электрическое поле Е (пусть для определенности оно будет направлено вдоль положительного направления оси х, рис. 2.13,б) , то у каждого электрона х-компонента волнового вектора kx одновременно получит отрицательное приращение. Этот вывод следует из уравнения движения, одинакового для каждого электрона:
. (2.20)
Следовательно, электроны валентной зоны будут перемещаться в направлении, указанном стрелкой на рис. 2.13,б. Вакантное состояние в результате этого движения электронов вначале переместится в точку Е, а затем — в точку D и т.д. Таким образом, последовательное перемещение электронов по энергетическим уровням под влиянием электрического поля эквивалентно перемещению вакантного состояния. Квантовое состояние, не занятое электроном в энергетической зоне, называется дыркой. Суммарный волновой вектор электронов в полностью заполненной энергетической зоне равен нулю, поскольку дисперсионная кривая симметрична относительно точки k = 0 и каждому электрону с волновым вектором k всегда найдется электрон с противоположным по знаку волновым вектором — k . Если из состояния с волновым вектором ke удален электрон, то полный волновой вектор системы станет равным — k e . Таким образом, дырке следует приписать волновой вектор
. (2.21)
Учитывая (2.20) и (2.21), уравнение движения дырки будет иметь вид
. (2.22)
Это уравнение движения положительного заряда в электрическом поле. Поскольку дырка перемещается вдоль направления действующей на нее силы, то этой частице следует приписать положительную эффективную массу, равную по абсолютному значению отрицательной эффективной массе электрона, покинувшего вакантное состояние у потолка валентной зоны.
Вычислим ток, создаваемый электронами полностью заполненной энергетической зоны. Вклад в плотность тока от одного электрона, движущегося со скоростью vj равен
.
Ток всех электронов валентной зоны равен сумме токов отдельных электронов:
.
Суммирование производится по всем состояниям, занятым электронами. Поскольку дисперсионные кривые симметричны, каждому электрону с ненулевым значением скорости в положительном направлении всегда найдется электрон с равной по абсолютному значению, но противоположно направленной скоростью. Следовательно, сила тока, создаваемого электронами полностью заполненной зоны, будет равна нулю.
Если в валентной зоне заняты все состояния, кроме одного, характеризующегося волновым вектором ks и скоростью vs (рис. 2.13,г), то суммарную плотность тока в этом случае можно представить в следующем виде:
.
В этой формуле учтено, что первое слагаемое в силу симметричности состояний электронов равно нулю.
Таким образом, движение электронов валентной зоны, в которой есть одно вакантное состояние, эквивалентно движению одной частицы с положительной эффективной массой и положительным электрическим зарядом, помещенной в это состояние.
2.9. Примесные полупроводники
В реальных кристаллах полупроводников всегда присутствуют, пусть и в небольших количествах, дефекты, примеси, некоторые из которых оказывают существенное влияние на их электропроводность. Например, добавление в кремний бора в количестве одного атома на 10 5 атомов кремния увеличивает его электропроводность при комнатной температуре в 1000 раз. Полупроводники, содержащие примеси, существенно влияющие на его электропроводность, называются примесными полупроводниками, а их электропроводность — примесной электропроводностью.
Рассмотрим механизм примесной проводимости на примере полупроводникового кристалла кремния с примесными атомами фосфора. Четыре валентных электрона кремния образуют в химически чистом кристалле парные ковалентные связи с четырьмя своими ближайшими соседями (рис. 2.14,а). Примесный атом фосфора замещает один из атомов кремния в узле кристаллической решетки. У атома фосфора пять валентных электронов, четыре из которых поддерживают связи с соседними атомами кремния, а пятый остается свободным (рис. 2.14,б). Этот избыточный электрон может перейти в зону проводимости кремния и «участвовать» в создании электрического тока. Примеси, поставляющие в зону проводимости дополнительное количество электронов, называются донорными примесями, а полупроводники с такими примесями — донорными полупроводниками или полупроводниками n-типа. Наиболее распространенными донорными примесями в кристаллах кремния и германия являются атомы пятой группы периодической системы элементов Д. И. Менделеева: фосфор (P), мышьяк (As), сурьма (Sb), висмут (Bi). Энергию, которую необходимо затратить, чтобы перевести электрон примесного донорного атома в зону проводимости, называют энергией связи донорной примеси. Оценить энергию связи донорной примеси можно из простой модели, подобной боровской модели атома водорода. Согласно этой модели примесный электрон движется по круговой орбите в кулоновском поле сил иона фосфора подобно электрону в поле ядра атома водорода. Различие заключается в том, что поле примесного иона ослаблено диэлектрическими свойствами кристалла полупроводника. Это влияние учитывается диэлектрической проницаемостью среды, которая для типичных полупроводников составляет 5 . 2000. Необходимо учесть также тот факт, что эффективная масса электрона в кристалле отличается от массы свободного электрона. Для количественных оценок воспользуемся результатами, полученными в теории Бора для атома водорода. Энергия связи электрона в атоме водорода равна . Учитывая диэлектрическую проницаемость полупроводника e и заменяя массу свободного электрона m на его эффективную массу в кристалле m*, получим следующее выражение для энергии ионизации донорной примеси:
Рис. 2.14 Схема проводимости в донорном полупроводнике:
а — ковалентные связи в чистом полупроводнике кремния; б — примесный атом фосфора;
в — зонная структура донорного полупроводника
. (2.23)
Энергия ионизации свободного атома водорода равна 13,6 эВ. В соответствии с формулой (2.23) это значение надо умножить на коэффициент , чтобы получить величину Ed. В кремнии e = 11,7; m*/m » 0,2. В результате получим Ed » 0,02 эВ.
Экспериментальное значение энергии ионизации фосфора в кремнии составляет 0,044 эВ. Другие донорные примеси имеют в кремнии и германии энергию ионизации того же порядка величины (см. таблицу).
Видео:Силаев П. К. - Квантовая теория - Нестационарное уравнение Шредингера (Лекция 12)Скачать
Уравнение Шредингера для кристалла
Возможные значения энергии электрона в кристалле можно найти из решения уравнения Шредингера для стационарных состояний, которое в общем случае имеет вид:
HY = EY, (1.31)
где H – оператор Гамильтона для кристалла, E – полная энергия кристалла, Y — волновая функция кристалла.
Волновая функция Ψ(x,t) обычно определяется из выражения для плотности вероятности нахождения частицы в точке x в момент времени t.
В общем случае волновая функция является комплексной величиной и имеет следующий вид:
(1.32)
Так как вероятность должна быть величиной действительной, то для нахождения плотности вероятности необходимо умножить Ψ(x,t) на комплексно сопряженную с ней функцию Ψ * (x,t). Комплексно-сопряженная волновая функция:
(1.33)
Квадрат амплитуды волновой функции характеризует вероятность (P) нахождения электрона в данной точке пространства, определяемой радиус-вектором .
(1.34)
или более точно:
Квадрат модуля волновой функции, умноженный на элемент объема есть вероятность обнаружить электрон в момент времени t в объеме , образующий бесконечно малую окрестность точки .
Уравнение Шредингера должно быть дополнено ограничениями, накладываемыми на волновую функцию:
— волновая функция должна быть конечной (так как вероятность не может быть больше 1), однозначной (вероятность не может быть неоднозначной величиной) и непрерывной (вероятность не может изменяться скачком).
— производные ¶Y/¶x, ¶Y/¶y, ¶Y/¶z должны быть непрерывны
— функция |Y| 2 должна быть интегрируема. Это условие в простейших случаях сводится к условию нормировки вероятности (т.к. электрон имеется, вероятность его существования где-то в пространстве — 100%-ная): (1.40) | Рис. 1.40. Плотность вероятности нахождения электрона в пространстве. |
Твердое тело представляет собой единую систему легких (электронов) и тяжелых (ядра) частиц. Наиболее полные сведения о свойствах такой системы, в том числе и об ее энергетическом спектре, можно получить решая уравнение Шредингера, соответствующее стационарным состояниям этой системы. Но, из-за того, что в кристалле имеет место большое число взаимодействий, уравнение Шредингера приобретает сложный вид:
(1.41)
В уравнение (2.3) входят следующие виды энергии:
1) Кинетическая энергия электронов: , где -постоянная Дирака ( = h/2 ), m-масса электрона, а — оператор Лапласа для i-го электрона;
2) Кинетическая энергия ядер: , где -масса ядра;
3) Потенциальная энергия попарного взаимодействия электронов между собой: , где -заряд электрона, здесь i ≠ j(!)
4) Потенциальная энергия попарного взаимодействия ядер между собой:
5) Потенциальная энергия взаимодействия электронов с ядрами: ,
где — координаты электронов и ядер соответственно;
6) Полная энергия кристалла E .
— собственная волновая функция кристалла, зависящая от координат всех электронов и атомных ядер.
Уравнение (21) содержит 3N(Z+1) независимых переменных, где N – число атомов в кристалле, Zq0 – заряд ядра. Число атомов в 1 см 3 составляет примерно 5*10 22 и каждый атом содержит большое число электронов. Следовательно, волновая функция зависит от огромного числа (10 24 -10 25 ) независимых переменных (в кремнии: 3N(Z+1) = 3 · 5х10 22 · (14+1) = 2.25х10 24 переменных). Точное решение уравнения Шредингера невозможно даже для отдельных атомов, за исключением атома водорода. Поэтому задача сводится к нахождению приближений в рамках физически оправданных упрощающих предположений. Зонная теория, лежащая в основе современной физики полупроводников, базируется на следующих основных приближениях:
1.3.4. Эффективная масса носителей заряда
На свободный электрон, помещенный в однородное электрическое поле E, действует сила F=-qE, под действием которой электрон приобретает ускорение . Здесь т — масса электрона. Поскольку является единственной силой, определяющей характер движения частицы, вектор ускорения электрона направлен так же, как и вектор внешней силы, т.е. против поля E.
В кристалле, внешнее поле E действует на электрон так же, как на свободный электрон: с силой F=-qE, направленной против поля. Однако, кроме силы -qE, на электрон действуют значительные внутренние силы, создаваемые периодическим полем решетки. Это означает, что ускорение электрона в решетке в общем случае может быть не направлено параллельно внешней силе . Следовательно, движение электрона в кристалле будет более сложным, чем движение свободного электрона. Именно поэтому энергия электрона, движущегося в периодическом поле кристалла, не имеет квадратичной зависимости от волнового вектора E ≠ p 2 /2m. Это также следует из теоретических моделей (например, модели Кронига и Пенни). Однако для практических целей иногда удобно сохранить зависимость энергии электрона от квазиимпульса в классическом виде, а все различия, вызванные влиянием периодического поля, включить в массу электрона. Тогда формулу Е р 2 /2m можно переписать в виде Е=р 2 /2m * , где вместо т появляется некоторая величина т*, которая является функция энергии, и которая называемая эффективной массой. Для определения понятия эффективной массы можно воспользоваться несколькими подходами.
(а) Определение эффективной массы из разложения в ряд Тейлора (формальный, математический метод). В одномерном случае величину т* можно рассчитать из разложения энергии в ряд Тейлора около экстремальных точек (т. к. E(k)- периодическая функция k и E(k)
cos(k) (периодический член))
(1.71)
Так как в точках k=ko энергия имеет максимум или минимум (см. рис. 1.46, 1.47), то первая производная равна нулю. Ограничиваясь вторым приближением, из (1.71) находим:
(1.72)
Или: (1.73)
Если E(k)-E(k0)=E(k1), а k-k0=k1, тогда ,
Следовательно, роль эффективной массы играет величина
(1.74)
В низших точках разрешенных зон Е(k) имеет минимумы, а вторая производная от Е по k больше нуля. Поэтому на дне зоны эффективная масса положительна, а в вершинах зон отрицательна, поскольку d 2 E/ dk 2 удовлетворяют тем же соотношениям, что и соответствующие им классические величины. Так, волновые пакеты, составленные из решений уравнения Шредингера, движутся по траекториям классических частиц. Поэтому уравнению Ньютона должен соответствовать квантовомеханический аналог, то есть квантовомеханическое уравнение движения электрона в кристалле.
Таким образом, можно использовать квазиклассический подход (наполовину классический, наполовину квантовомеханический). В этом подходе квантовомеханическим, является определение электрона как волны. В этом случае движение электрона в кристалле можно описать с помощью волнового пакета, составленного из блоховских функций. Тогда средняя скорость движения электрона равна групповой скорости распространения всего волнового пакета — . Учитывая, что , для групповой скорости получаем:
(1.75)
Аналогичным образом можно ввести усредненное ускорение волнового пакета (электрона) в кристалле:
(1.76)
С учетом того, что время и квазиимпульс независимы, в (1.76) можно поменять местами порядок их дифференцирования:
(1.77)
Классическим,в настоящем подходе,является определение работы внешней силы над рассматриваемым электроном, которая приводит к увеличению энергии (скорости) электрона:
, откуда или (1.78)
Подставляя (1.78) в (1.77) и учитывая, что внешняя сила F не зависит от k, получим:
(1.79)
Перепишем в виде: (1.80)
Т.е. по аналогии с законом Ньютона, величину можно назвать массой электрона:
(1.81)
Величина т* получила название эффективной массы электрона. Эффективная масса отражает влияние периодического потенциала решетки на движение электрона в кристалле под действием внешней силы. В периодическом поле кристаллической решетки электрон движется под действием внешней силы F в среднем так, как двигался бы свободный электрон под действием этой силы, если бы он обладал массой т* Таким образом, если электрону в кристалле вместо массы т приписать эффективную массу т*, то его можно считать свободным и движение этого электрона описывать так, как описывается движение свободного электрона, помещенного во внешнее поле. Разница между т* и m обусловлена взаимодействием электрона с периодическим полем решетки, и, приписывая электрону эффективную массу, мы учитываем это взаимодействие.
Эффективная масса, в отличие от обычной массы, не определяет ни инерционных, ни гравитационных свойств частицы. Она является лишь коэффициентом в уравнении движения и отражает меру взаимодействия электрона с кристаллической решеткой. |
Пользуясь понятием эффективной массы, задачу о движении электрона в периодическом поле решетки V( ) можно свести к задаче о движении свободного электрона с массой т*. Это значит, что вместо уравнения Шредингера с периодическим потенциалом
, нужно решать уравнение
, (1.82)
Метод решения уравнения Шредингера, в котором вид периодического потенциала решетки автоматически учитывается через эффективную массу, называется методом эффективной массы
В общем случае эффективная масса является анизотропной величиной и для разных направлений волнового вектора различна. Она представляет собой тензор второго ранга:
или (1.83)
Таким образом, если зависимость [закон дисперсии] анизотропна, то эффективная масса представляет собой тензор обратных эффективных масс.
Рассмотрим некоторые свойства эффективной массы.
Из формул (1.74), (1.81), (1.85) следует, что эффективная масса определяется видом дисперсионной зависимости E(k). В приближении сильной связи выражение E(k) имеет вид (1.70), как это следует из модели Кронига-Пенни. Периодичность волновой функции (1.51), (1.52) и энергии электрона в кристалле (1.54) позволяет ограничиться рассмотрением зависимости E(k) в первой зоне Бриллюэна. Возьмем первую (нижнюю) ветку дисперсионной кривой E(k) в первой зоне Брилюэна для кристалла кубической сингонии (рис. 1.48а) и проведем дифференцирование. Поскольку E(k) имеет функциональную зависимость, близкую к E(k)
cos (ka) (см. формулу (1.70)), первая производная вблизи точек экстремума будет близкой к зависимости dE/dk
sin (ka). Во всей первой зоне Бриллюэна, зависимость dE/dk есть зависимость скорости электрона от волнового вектора k: (dE/dk=Vg).
Вторая производная и ее обратная зависимость (пропорциональная эффективной массе) приведены на рис. 1.48. Сравнение зависимости E(k) на рис. 1.48(а) и зависимости m*(k) на рис. 1.48(г) позволяет сделать следующие комментарии. Электрон, находящийся внутри идеальной периодической решетки может иметь как положительную, так и отрицательную эффективную массу. Если кривая на диаграмме E-k имеет выпуклость вниз (относительно оси энергий, если принять направление увеличения энергии — вверх), то масса m*>0. Если же кривая имеет выпуклость вверх (около ), тогда m* 2 E/dk 2 (в) и эффективной массы (г) от волнового вектора для кубической решетки.
Понятие дырки поясним следующим примером (рис. 1.49):
Предположим, в исходном состоянии валентная зона полностью заполнена электронами (нет свободных энергетических уровней), а зона проводимости – свободна. Включив внешний приток энергии (нагрев, облучение светом, радиацией и др.), можно инициировать переход электронов из валентной зоны в зону проводимости. Предположим далее, что энергия фотона передаётся электрону в валентной зоне и он переходит в зону проводимости (процесс А на рис. 1.49). С точки зрения химической связи это означает повреждение (разрыв) ковалентной связи и уход электрона в свободное перемещение по кристаллу. С точки зрения зонной теории, электрон уходит в зону проводимости, а на потолке валентной зоны образуется незанятое место (квантовое состояние).
Если внешнее электрическое поле E равно нулю и вследствие того, что электроны стремятся занять самые нижние энергетические состояния, дырка занимает самое верхнее состояние (позицию 1 на рис. 1.49). Под действием электрического поля E на это незанятое состояние перейдет электрон с более низкого энергетического уровня: на рис. 1.49 это обозначено переходом электрона из позиции 2 в позицию 1. Дырка при этом опустится из позиции 1 в позицию 2. Затем этот процесс может повториться переходом (3 à 2) и т.д. по эстафете. При таком перемещении электронов, освобождающееся незанятое место на энергетическом уровне перемещается вниз, вглубь валентной зоны. Этот освобождающийся уровень и отождествляется с некоей фиктивной частицей, которая имеет равный по значению, но противоположный по знаку заряд, по сравнению с электроном. Таким образом, свободное от электрона квантовое состояние вблизи потолка валентной зоны и называют дыркой.
Важно отметить, что при переходах 1 à 2 à 3 … дырка приобретет кинетическую энергию, и полная энергия дырки возрастет (. ), т.е. шкала энергий в валентной зоне возрастает в направлении – противоположном зоне проводимости.
Классический аналог увеличения энергии по зонам: взвешенная капля воды в воздухе (туман) – электрон в зоне проводимости, пузырек воздуха в воде – дырка в валентной зоне.
Таким образом, ток в кристаллах может переноситься не только электронами в зоне проводимости, но и дырками в валентной зоне. Дырочная проводимость наиболее характерна для полупроводников, однако есть и некоторые металлы, которые обладают дырочной проводимостью.
Возвращаясь к рис. 1.48в, отметим, что описывать движение электронов в кристалле, пользуясь понятием эффективной массы, можно только тогда, когда они находятся либо у дна, либо у потолка энергетической зоны. В центре зоны (около значений k = ± π/a) m*à∞, т.е. понятие эффективной массы теряет смысл. На практике почти всегда приходится иметь дело с электронами, располагающимися или у дна, или у потолка зоны. Поэтому использование эффективной массы в этих случаях вполне оправдано.
Ширина разрешенных зон растет, а запрещенных – уменьшается, с увеличением энергии (порядкового номера разрешенной зоны). Т.е. в общем случае зона проводимости, имеет большую ширину, чем валентная зона. Поскольку эффективная масса обратно пропорциональна ширине энергетической зоны, поэтому, как правило, . Следовательно, дырка и электрон проводимости отличаются не только знаком своего заряда, но имеют и разные по величине m*.
Иной результат получается, если в почти заполненной зоне имеются свободные места, то есть не все валентные связи обеспечены электронами. Тогда соседние электроны могут переходить на эти места, а само свободное место как бы перемещается в пространстве. В энергетическом плане это соответствует переходу электронов с низких энергетических уровней на более высокие, а дырок с высоких уровней на более низкие.
Резюмируя полученную информацию:
1. Обозначения носителей зарядов в полупроводнике.
Электроны | n, e (negative, electrons) |
Дырки | p, h (positive, holes) |
Ток частично заполненной зоны может быть представлен как ток положительно заряженных частиц дырок. Заряд дырки положителен и по величине равен заряду электрона. Концентрация дырок обычно обозначается буквой р. |
2. Три представления (определения) дырок:
(a) полноправная положительно заряженная частица, перемещающая в кристалле.
(b) отсутствие электрона в потолке валентной зоны.
(c) физическое отсутствие электрона в том месте, где он должен быть в равновесном состоянии – т.е. в составе ковалентной связи.
3. Направление энергии в зонах: зона проводимости – энергия увеличивается вверх; валентная зона – энергия увеличивается вниз.
4. Величина m* зависит от кривизны зоны (m*
5.Ширина зон увеличивается с E, а m * обратно пропорциональна ширине энергетической зоны, à (a) зона проводимости энергетически шире, чем валентная; (b) (как правило).
💡 Видео
Волновая функция (видео 5) | Квантовая физика | ФизикаСкачать
Структура материи 6: уравнение Шрёдингера. Зачем нужна квантовая механика – Виталий Бейлин | НаучпопСкачать
Лекции 5-6. Уравнение Шредингера и его приближенные решения. Межатомные.Скачать
Классические уравнения | уравнение Шрёдингера (координатное представление) | простейший выводСкачать
96. Уравнение ШредингераСкачать
Урок 32. Уравнение ШрёдингераСкачать
Квантовая механика 41 - Уравнение Шредингера. Гамильтониан.Скачать
10. Уравнение ШрёдингераСкачать
Силаев П. К. - Квантовая теория - Нестационарное уравнение Шредингера (Лекция 11)Скачать
Уравнение ШрёдингераСкачать
Силаев П. К. - Квантовая теория - Нестационарное уравнение Шредингера (Лекция 14)Скачать
7.1 Разделение переменных в уравнении ШрёдингераСкачать
Корректный вывод уравнения Шрёдингера и его физический смысл: Липовка А.А. - Глобальная волнаСкачать
Уравнение Шредингера Стационарные состоянияСкачать
Сущёв И. - Квантовая теория - 4. Уравнение Шредингера (разбор задач)Скачать
5_4. Уравнение ШредингераСкачать
Лекция №4 "Волновая функция. Уравнение Шредингера" (Гавриков А.В.)Скачать
Периодический потенциал 0 (Зоны Бриллюэна)Скачать