Эффективная масса (физика твердого тела)

Чтобы узнать об отрицательной массе в теоретической физике , см. Отрицательную массу .

В физике твердого тела частицы эффективная масса (часто обозначаемая ${\textstyle m^{*}}$) — это масса , которую он имеет при реакции на силы, или масса, которую он имеет при взаимодействии с другими идентичными частицами в тепловом распределении . Одним из результатов зонной теории твердого тела является то, что движение частиц в периодическом потенциале на большие расстояния, превышающие период решетки, может сильно отличаться от их движения в вакууме. Эффективная масса — это величина, которая используется для упрощения зонных структур путем моделирования поведения свободной частицы с этой массой. Для некоторых целей и некоторых материалов эффективную массу можно считать простой константой материала. Однако в целом значение эффективной массы зависит от цели, для которой она используется, и может варьироваться в зависимости от ряда факторов.

Для электронов или электронных дырок в твердом теле эффективная масса обычно указывается как коэффициент, умножающий массу покоя электрона , m _e (9,11 × 10 ⁻³¹ кг). Этот коэффициент обычно находится в диапазоне от 0,01 до 10, но может быть ниже или выше — например, достигая 1000 в экзотических материалах с тяжелыми фермионами или от нуля до бесконечности (в зависимости от определения) в графене . Поскольку это упрощает более общую зонную теорию, электронную эффективную массу можно рассматривать как важный базовый параметр, который влияет на измеримые свойства твердого тела, включая все, от эффективности солнечного элемента до скорости интегральной схемы.

Простой случай: параболическое

При самых высоких энергиях валентной зоны во многих полупроводниках (Ge, Si, GaAs,...) и самых низких энергиях зоны проводимости в некоторых полупроводниках (GaAs,...) зонная структура E ( k ) может быть локально аппроксимировано как

$E(\mathbf {k} )=E_{0}+{\frac {\hbar ^{2}\mathbf {k} ^{2}}{2m^{*}}}$

где E ( k ) — энергия электрона на волновом векторе k в этой зоне, E ₀ — константа, задающая край энергии этой зоны, и m ^* – константа (эффективная масса).

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

Одним из замечательных свойств является то, что эффективная масса может стать отрицательной , когда полоса изгибается вниз от максимума. В результате отрицательной массы электроны реагируют на электрические и магнитные силы, набирая скорость в направлении, противоположном нормальному; даже несмотря на то, что эти электроны имеют отрицательный заряд, они движутся по траекториям, как если бы они имели положительный заряд (и положительную массу). Это объясняет существование дырок валентной зоны с положительным зарядом и положительной массой, — квазичастиц которые можно обнаружить в полупроводниках.

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

Промежуточный случай: параболическое

Эллипсоиды постоянной энергии в кремнии вблизи шести минимумов зоны проводимости. Для каждой долины (минимум полосы) эффективные массы составляют m _ℓ = 0,92 м _e («продольное»; вдоль одной оси) и m _t = 0,19 м _e («поперечное»; вдоль двух осей).

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

$E\left(\mathbf {k} \right)=E_{0}+{\frac {\hbar ^{2}}{2m_{x}^{*}}}\left(k_{x}-k_{0,x}\right)^{2}+{\frac {\hbar ^{2}}{2m_{y}^{*}}}\left(k_{y}-k_{0,y}\right)^{2}+{\frac {\hbar ^{2}}{2m_{z}^{*}}}\left(k_{z}-k_{0,z}\right)^{2}$

где оси x , y и z совпадают с главными осями эллипсоидов, а m ^*
_х , м ^*
_й и м ^*
_z — инерционные эффективные массы вдоль этих разных осей. Смещения k _{0, x} , k _{0, y} и k _{0, z} отражают то, что минимум зоны проводимости больше не центрируется на нулевом волновом векторе. (Эти эффективные массы соответствуют главным компонентам тензора инерционной эффективной массы, описанного позже. )

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

• Для расчета полной плотности состояний и полной плотности носителей через среднее геометрическое в сочетании с коэффициентом вырождения g , который подсчитывает количество впадин (в кремнии g = 6 ):

  $m_{\text{density}}^{*}={\sqrt[{3}]{g^{2}m_{x}m_{y}m_{z}}}$

  (Эта эффективная масса соответствует эффективной массе плотности состояний, описанной позже.)

  Для плотности состояний на долину и плотности носителей на долину фактор вырождения не учитывается.
• Для расчета проводимости, как в модели Друде, через среднее гармоническое

  $m_{\text{conductivity}}^{*}=3\left[{\frac {1}{m_{x}^{*}}}+{\frac {1}{m_{y}^{*}}}+{\frac {1}{m_{z}^{*}}}\right]^{-1}$

  Поскольку закон Друде также зависит от времени рассеяния, которое сильно варьируется, эта эффективная масса используется редко; Вместо этого проводимость обычно выражается через плотность носителей и эмпирически измеряемый параметр — подвижность носителей .

Общий случай

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

массы

Классическая частица под действием силы ускоряется по второму закону Ньютона , a = m ⁻¹ F или, альтернативно, импульс изменяется в соответствии с d / d т п знак равно F . Этот интуитивный принцип тождественно проявляется в квазиклассических приближениях, основанных на зонной структуре, когда межзонными переходами можно пренебречь при достаточно слабых внешних полях. Сила дает скорость изменения импульса кристалла p _{кристалла} :

$\mathbf {F} ={\frac {\operatorname {d} \mathbf {p} _{\text{crystal}}}{\operatorname {d} t}}=\hbar {\frac {\operatorname {d} \mathbf {k} }{\operatorname {d} t}},$

где ħ = h /2π — приведенная постоянная Планка .

Ускорение волнообразной частицы становится скоростью изменения групповой скорости :

$\mathbf {a} ={\frac {\operatorname {d} }{\operatorname {d} t}}\,\mathbf {v} _{\text{g}}={\frac {\operatorname {d} }{\operatorname {d} t}}\left(\nabla _{k}\,\omega \left(\mathbf {k} \right)\right)=\nabla _{k}{\frac {\operatorname {d} \omega \left(\mathbf {k} \right)}{\operatorname {d} t}}=\nabla _{k}\left({\frac {\operatorname {d} \mathbf {k} }{\operatorname {d} t}}\cdot \nabla _{k}\,\omega (\mathbf {k} )\right),$

где ∇ _k — оператор del в обратном пространстве . Последний шаг следует из использования цепного правила для полной производной для величины с косвенными зависимостями , поскольку прямым результатом действия силы является изменение k ( t ), приведенное выше, что косвенно приводит к изменению E ( k )= ħ ω( k ) . Объединение этих двух уравнений дает

$\mathbf {a} =\nabla _{k}\left({\frac {\mathbf {F} }{\hbar }}\cdot \nabla _{k}\,{\frac {E(\mathbf {k} )}{\hbar }}\right)={\frac {1}{\hbar ^{2}}}\left(\nabla _{k}\left(\nabla _{k}\,E(\mathbf {k} )\right)\right)\cdot \mathbf {F} =M_{\text{inert}}^{-1}\cdot \mathbf {F}$

используя правило скалярного произведения с равномерной силой ( ∇ _k F =0 ). $\nabla _{k}\left(\nabla _{k}\,E(\mathbf {k} )\right)$ — матрица Гессе E k ( в ) обратном пространстве. Мы видим, что эквивалент ньютоновской обратной инерционной массы для свободной частицы, определяемой формулой a = m ⁻¹ F стала тензорной величиной

$M_{\text{inert}}^{-1}={\frac {1}{\hbar ^{2}}}\nabla _{k}\left(\nabla _{k}\,E(\mathbf {k} )\right).$

чьи элементы

$\left[M_{\text{inert}}^{-1}\right]_{ij}={\frac {1}{\hbar ^{2}}}\left[\nabla _{k}\left(\nabla _{k}\,E(\mathbf {k} )\right)\right]_{ij}={\frac {1}{\hbar ^{2}}}{\frac {\partial ^{2}E}{\partial k_{i}\partial k_{j}}}\,.$

Этот тензор позволяет ускорению и силе быть в разных направлениях, а величина ускорения зависит от направления силы.

• Для параболических полос недиагональные элементы M _inert ⁻¹ равны нулю, а диагональные элементы являются константами
• Для изотропных полос все диагональные элементы должны быть равны, а все недиагональные элементы должны быть равны.
• Для параболических изотропных зон M _inert ⁻¹ = 1 / м ^* Я , где м ^* — скалярная эффективная масса, а I — единица.
• В целом элементы М _{инертны} ⁻¹ являются функциями k .
• Обратное: M _inert = ( M _inert ⁻¹ ) ⁻¹ , известен как тензор эффективной массы . не всегда возможно. Обратите внимание, что инвертировать M _inert ⁻¹

Для полос с линейной дисперсией $E\propto k$ например, с фотонами или электронами в графене , групповая скорость фиксирована, т.е. электроны, движущиеся параллельно k направлению силы F, не могут быть ускорены, а диагональные элементы M _{инертны.} ⁻¹ очевидно, равны нулю. Однако электроны, движущиеся с компонентой, перпендикулярной силе, могут ускоряться в направлении силы, а недиагональные элементы M _{инертны.} ⁻¹ ненулевые. Фактически, недиагональные элементы масштабируются обратно пропорционально k , т.е. они расходятся (становятся бесконечными) при малых k . Вот почему иногда говорят, что электроны в графене имеют бесконечную массу (из-за нулей на диагонали M _{инертного} ⁻¹ ) и иногда называют безмассовым (из-за расходимости на недиагоналях).

циклотрона

Классически заряженная частица в магнитном поле движется по спирали вдоль оси магнитного поля. Период Т его движения зависит от его массы т и заряда е ,

$T=\left\vert {\frac {2\pi m}{eB}}\right\vert$

где B — плотность магнитного потока .

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

Квазиклассическое движение частицы можно описать замкнутым контуром в k-пространстве. На протяжении всего этого цикла частица сохраняет постоянную энергию, а также постоянный импульс вдоль оси магнитного поля. Определив A как площадь k -пространства , заключенную в эту петлю (эта площадь зависит от энергии E , направления магнитного поля и волнового вектора на оси k _B ), можно показать, что эффективная масса циклотрона зависит от зонной структуры через производную этой площади по энергии:

$m^{*}\left(E,{\hat {B}},k_{\hat {B}}\right)={\frac {\hbar ^{2}}{2\pi }}\cdot {\frac {\partial }{\partial E}}A\left(E,{\hat {B}},k_{\hat {B}}\right)$

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

В двумерных электронных газах эффективная масса циклотрона определена только для одного направления магнитного поля (перпендикулярного), и внеплоскостной волновой вектор выпадает. Таким образом, эффективная масса циклотрона является только функцией энергии и оказывается точно связанной с плотностью состояний при этой энергии соотношением $\scriptstyle g(E)\;=\;{\frac {g_{v}m^{*}}{\pi \hbar ^{2}}}$, где g _v — вырождение долины. Такая простая зависимость неприменима к трехмерным материалам.

Плотность состояний эффективных масс (слаболегированные полупроводники)

Плотность состояний эффективной массы в различных полупроводниках

Группа	Материал	Электрон	Дыра
IV	Si (4  K)	1.06	0.59
	Си (300   К)	1.09	1.15
	Ге	0.55	0.37
III–V	GaAs	0.067	0.45
	InSb	0.013	0.6
II–VI	ZnO	0.29	1.21
	ZnSe	0.17	1.44

В полупроводниках с низким уровнем легирования концентрация электронов в зоне проводимости обычно определяется выражением

$n_{\text{e}}=N_{\text{C}}\exp \left(-{\frac {E_{\text{C}}-E_{\text{F}}}{kT}}\right)$

где E _F — уровень Ферми , EC _— минимальная энергия зоны проводимости, а NC _— коэффициент концентрации, зависящий от температуры. приведенное выше соотношение для n _e применимо для любой формы зоны проводимости (включая непараболические асимметричные зоны) при условии, что легирование слабое ( EC _{Можно показать, что} - E _F ≫ kT ); это следствие статистики Ферми – Дирака ограничения в сторону статистики Максвелла – Больцмана .

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

$\quad N_{\text{C}}=2\left({\frac {2\pi m_{\text{e}}^{*}kT}{h^{2}}}\right)^{\frac {3}{2}}$

В полупроводниках с непростой зонной структурой это соотношение используется для определения эффективной массы, известной как эффективная масса электронов в плотности состояний . масса плотности состояний» используется, поскольку приведенное выше выражение для NC _{Название « эффективная} получено через плотность состояний для параболической зоны.

На практике эффективная масса, экстрагированная таким способом, не совсем постоянна по температуре ( не _NC меняется точно при T ^3/2 ). В кремнии, например, эта эффективная масса варьируется на несколько процентов от абсолютного нуля до комнатной температуры, поскольку сама зонная структура слегка меняет форму. Эти искажения зонной структуры являются результатом изменения энергий электрон-фононного взаимодействия, при этом тепловое расширение решетки играет незначительную роль.

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

$n_{\text{h}}=N_{\text{V}}\exp \left(-{\frac {E_{\text{F}}-E_{\text{V}}}{kT}}\right),\quad N_{\text{V}}=2\left({\frac {2\pi m_{\text{h}}^{*}kT}{h^{2}}}\right)^{\frac {3}{2}}$

где E _V — максимальная энергия валентной зоны. На практике эта эффективная масса имеет тенденцию сильно различаться между абсолютным нулем и комнатной температурой во многих материалах (например, в два раза в кремнии), поскольку существует множество валентных зон с отчетливым и существенно непараболическим характером, все из которых достигают максимума вблизи одной и той же энергии. .

Определение

Экспериментальный

Традиционно эффективные массы измерялись с помощью циклотронного резонанса — метода, при котором микроволновое поглощение полупроводника, погруженного в магнитное поле, проходит через острый пик, когда микроволновая частота равна циклотронной частоте. $\scriptstyle f_{c}\;=\;{\frac {eB}{2\pi m^{*}}}$. В последние годы эффективные массы чаще определялись путем измерения зонных структур с использованием таких методов, как фотоэмиссионная спектроскопия с угловым разрешением ( ARPES ) или, что наиболее непосредственно, эффект де Хааса-Ван Альфена . Эффективные массы также можно оценить, используя коэффициент γ линейного члена низкотемпературной электронной удельной теплоемкости при постоянном объеме. $\scriptstyle c_{v}$. Удельная теплоемкость зависит от эффективной массы через плотность состояний на уровне Ферми и, как таковая, является мерой вырождения, а также кривизны зоны. Очень большие оценки массы носителя на основе измерений удельной теплоемкости привели к появлению концепции тяжелых фермионных материалов. носителей Поскольку подвижность зависит от соотношения времен жизни столкновений носителей $\tau$ до эффективной массы, массы в принципе могут быть определены на основе транспортных измерений, но этот метод непрактичен, поскольку вероятности столкновения носителей обычно не известны априори. Оптический эффект Холла — это новый метод измерения плотности свободных носителей заряда, эффективной массы и параметров подвижности в полупроводниках. Оптический эффект Холла является аналогом квазистатического электрического эффекта Холла, индуцированного электрическим полем, на оптических частотах в проводящих и сложных слоистых материалах. Оптический эффект Холла позволяет также охарактеризовать анизотропию (тензорный характер) параметров эффективной массы и подвижности.

Теоретический

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

Значение

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

Некоторые соединения групп III – V , такие как арсенид галлия (GaAs) и антимонид индия (InSb), имеют гораздо меньшие эффективные массы, чем тетраэдрические группы IV, материалы такие как кремний и германий . В простейшей картине электронного транспорта Друде максимально достижимая скорость носителя заряда обратно пропорциональна эффективной массе: ${\textstyle {\vec {v}}\;=\;\left\Vert \mu \right\Vert \cdot {\vec {E}}}$, где ${\textstyle \left\Vert \mu \right\Vert \;=\;{e\tau }/{\left\Vert m^{*}\right\Vert }}$ с ${\textstyle e}$ электронный заряд . Предельная скорость интегральных схем зависит от скорости несущей, поэтому низкая эффективная масса является основной причиной того, что GaAs и его производные используются вместо Si в с высокой пропускной способностью приложениях , таких как сотовая телефония .

В апреле 2017 года исследователи из Университета штата Вашингтон заявили, что создали жидкость с отрицательной эффективной массой внутри конденсата Бозе-Эйнштейна , разработав соотношение дисперсии .

См. также

Модели твердых тел и кристаллов:

• Модель жесткой привязки
• Модель свободных электронов
• Модель почти свободных электронов

Ссылки

• Пастори Парравичини, Г. (1975). Электронные состояния и оптические переходы в твердых телах . Пергамон Пресс . ISBN  978-0-08-016846-3 . Эта книга содержит исчерпывающее, но доступное обсуждение темы с обширным сравнением расчетов и эксперимента.
• Пекар С., Метод эффективной массы электронов в кристаллах, Журн. Эксп. Теор. Физ. 16 , 933 (1946).

Внешние ссылки

• Архив НСМ