В физике конденсированного состояния плотность состояний ( DOS ) системы описывает количество разрешенных режимов или состояний на единицу энергетического диапазона. Плотность состояний определяется как , где – число состояний в системе объёма , энергии которых лежат в диапазоне от до . Математически оно представляется как распределение с помощью функции плотности вероятности и обычно представляет собой среднее значение по пространственным и временным областям различных состояний, занимаемых системой. Плотность состояний напрямую связана с дисперсионными соотношениями свойств системы. Высокий DOS на определенном энергетическом уровне означает, что многие состояния доступны для оккупации.
В общем случае плотность состояний материи непрерывна. Однако в изолированных системах , таких как атомы или молекулы в газовой фазе, распределение плотности дискретно , как спектральная плотность . Локальные вариации, чаще всего возникающие из-за искажений исходной системы, часто называют локальными плотностями состояний (LDOS).
В квантово-механических системах волны или волноподобные частицы могут занимать моды или состояния с длинами волн и направлениями распространения, определяемыми системой. Например, в некоторых системах межатомное расстояние и атомный заряд материала могут позволять существовать только электронам определенной длины волны. В других системах кристаллическая структура материала может позволять волнам распространяться в одном направлении, подавляя распространение волн в другом направлении. Часто разрешены только определенные состояния. Таким образом, может случиться так, что на определенном энергетическом уровне доступно множество состояний, в то время как на других энергетических уровнях нет доступных состояний.
Если посмотреть на плотность состояний электронов на краю зоны между валентной зоной и зоной проводимости в полупроводнике, то для электрона в зоне проводимости увеличение энергии электрона делает больше состояний доступными для заполнения. Альтернативно, плотность состояний является прерывистой для определенного интервала энергии, что означает, что электронам не доступны состояния, которые они могли бы занять в запрещенной зоне материала. Это условие также означает, что электрон на краю зоны проводимости должен потерять как минимум энергию запрещенной зоны материала, чтобы перейти в другое состояние в валентной зоне.
Это определяет, является ли материал изолятором или металлом в измерении распространения. Результат количества состояний в зоне также полезен для прогнозирования свойств проводимости. Например, в одномерной кристаллической структуре нечетное число электронов на атом приводит к наполовину заполненной верхней зоне; на уровне Ферми находятся свободные электроны, в результате чего получается металл. С другой стороны, четное число электронов заполняет ровно целое число зон, оставляя остальные пустыми. Если тогда уровень Ферми находится в заполненной запрещенной зоне между самым высоким занятым состоянием и самым низким пустым состоянием, материал будет изолятором или полупроводником .
В зависимости от квантово-механической системы плотность состояний может быть рассчитана для электронов , фотонов или фононов и может быть задана как функция либо энергии, либо волнового вектора k . Чтобы преобразовать DOS как функцию энергии и DOS как функцию волнового вектора, необходимо знать специфическое для системы соотношение дисперсии энергии между E и k .
В общем, топологические свойства системы, такие как зонная структура, оказывают большое влияние на свойства плотности состояний. Наиболее известные системы, такие как нейтронное вещество в нейтронных звездах и свободные электронные газы в металлах (примеры вырожденного вещества и ферми-газа ), имеют трехмерную евклидову топологию . Менее знакомые системы, такие как двумерные электронные газы (2DEG) в графитовых слоях и система квантового эффекта Холла в устройствах типа MOSFET , имеют двумерную евклидову топологию. Еще менее известны углеродные нанотрубки , квантовая проволока и жидкость Латтинжера с их одномерной топологией. Системы с 1D и 2D топологиями, вероятно, станут более распространенными, если предположить, что развитие нанотехнологий и материаловедения продолжится.
Плотность состояний, связанных с объемом V и N счетными уровнями энергии, определяется как:
Для изотропных одномерных систем с параболической дисперсией энергии плотность состояний равна . В двух измерениях плотность состояний постоянна , а в трех измерениях она становится .
Эквивалентно, плотность состояний также можно понимать как производную микроканонической статистической суммы (то есть общего числа состояний с энергией меньше ) по энергии:
Число состояний с энергией (степенью вырождения) определяется выражением:
Существует большое разнообразие систем и типов состояний, для которых можно выполнить расчеты DOS.
Некоторые системы конденсированного состояния обладают структурной симметрией в микроскопическом масштабе, которую можно использовать для упрощения расчета их плотностей состояний. В сферически-симметричных системах интегралы от функций одномерны, поскольку все переменные в расчете зависят только от радиального параметра дисперсионного уравнения. Жидкости , стекла и аморфные твердые тела являются примерами симметричной системы, дисперсионные соотношения которой обладают вращательной симметрией.
Измерения на порошках или поликристаллических образцах требуют оценки и расчета функций и интегралов по всей области , чаще всего зоне Бриллюэна , дисперсионных соотношений интересующей системы. Иногда симметрия системы высока, что приводит к многократному появлению формы функций, описывающих дисперсионные соотношения системы, во всей области действия дисперсионного уравнения. В таких случаях усилия по расчету DOS могут быть значительно уменьшены, если расчет ограничен уменьшенной зоной или фундаментальной областью . [1] Зона Бриллюэна гранецентрированной кубической решетки (ГЦК) на рисунке справа имеет 48-кратную симметрию точечной группы Oh с полной октаэдрической симметрией . Такая конфигурация означает, что интегрирование по всей области зоны Бриллюэна может быть сведено к 48-й части всей зоны Бриллюэна. Как показывает периодическая таблица кристаллической структуры , существует множество элементов с кристаллической структурой FCC, таких как алмаз , кремний и платина , а их зоны Бриллюэна и дисперсионные соотношения имеют 48-кратную симметрию. Двумя другими известными кристаллическими структурами являются объемно-центрированная кубическая решетка (BCC) и гексагональные закрытые упакованные структуры (HCP) с кубической и гексагональной решетками соответственно. Структура ОЦК имеет 24-кратную пиритоэдрическую симметрию точечной группы Th . Структура ГПУ имеет 12-кратную призматическую диэдральную симметрию точечной группы D 3h . Полный список свойств симметрии точечной группы можно найти в таблицах символов точечных групп .
В целом DOS легче рассчитать, когда симметрия системы выше, а количество топологических размерностей дисперсионного уравнения меньше. DOS дисперсионных соотношений с вращательной симметрией часто можно рассчитать аналитически. Этот результат удачен, поскольку многие материалы, представляющие практический интерес, такие как сталь и кремний, обладают высокой симметрией.
В анизотропных конденсированных системах, таких как монокристалл соединения, плотность состояний может быть разной в одном кристаллографическом направлении, чем в другом. Из-за этого анизотропную плотность состояний становится сложнее визуализировать, и могут потребоваться такие методы, как расчет DOS только для определенных точек или направлений или расчет проецируемой плотности состояний (PDOS) для определенной ориентации кристалла.
Плотность состояний зависит от размерных ограничений самого объекта. В системе, описываемой тремя ортогональными параметрами (3-мерное измерение), единицами DOS являются [Энергия] −1 [Объем] −1 , в двумерной системе единицами DOS являются [Энергия] −1 [Площадь] −1. , в одномерной системе единицами DOS являются [Энергия] −1 [Длина] −1 . Указанный объем — это объем k -пространства; пространство, ограниченное поверхностью постоянной энергии системы, полученной посредством дисперсионного соотношения , которое связывает E с k . Пример трехмерного k -пространства приведен на рис. 1. Видно, что размерность системы ограничивает импульс частиц внутри системы.
Расчет для DOS начинается с подсчета N разрешенных состояний при определенном k , которые содержатся в пределах [ k , k + d k ] внутри объема системы. Эта процедура выполняется путем дифференцирования всего объема k-пространства в n-мерностях при произвольном k относительно k . Объем, площадь или длина в 3-, 2- или 1-мерном сферическом k- пространстве выражаются выражением
для n -мерного k -пространства с топологически определенными константами
Согласно этой схеме плотность состояний волнового вектора N путем дифференцирования по k выражается выражением
1-, 2- и 3-мерная плотность состояний волнового вектора для линии, диска или сферы явно записывается как
Одно состояние достаточно велико, чтобы содержать частицы с длиной волны λ. Длина волны связана с k посредством соотношения.
В квантовой системе длина λ будет зависеть от характерного расстояния системы L, удерживающей частицы. Наконец, плотность состояний N умножается на коэффициент , где s — постоянный коэффициент вырождения, учитывающий внутренние степени свободы, обусловленные такими физическими явлениями, как спин или поляризация. Если такого явления нет, то . V k — это объем в k-пространстве, волновые векторы которого меньше наименьших возможных волновых векторов, определяемых характерным расстоянием системы.
Чтобы завершить расчет для DOS, найдите количество состояний на единицу объема выборки при энергии внутри интервала . Общий вид DOS системы задается как
Закон дисперсии электронов в твердом теле задается электронной зонной структурой .
Кинетическая энергия частицы зависит от величины и направления волнового вектора k , свойств частицы и среды, в которой частица движется. Например, кинетическая энергия электрона в ферми -газе определяется выражением
где m — масса электрона . Дисперсионное соотношение представляет собой сферически-симметричную параболу и постоянно возрастает, поэтому DOS можно легко рассчитать.
Для продольных фононов в цепочке атомов закон дисперсии кинетической энергии в одномерном k- пространстве, как показано на рисунке 2, определяется выражением
Когда энергия
С помощью преобразования и малого это отношение можно преобразовать в
Два упомянутых здесь примера можно выразить так:
Это выражение является своего рода дисперсионным соотношением , поскольку оно связывает два волновых свойства и является изотропным, поскольку в выражении фигурирует только длина, а не направление волнового вектора. Величина волнового вектора связана с энергией следующим образом:
Соответственно, объем n-мерного k -пространства, содержащего волновые векторы, меньшие k , равен:
Замена изотропного энергетического соотношения дает объем занятых состояний
Дифференцирование этого объема по энергии дает выражение для ПЭС изотропного дисперсионного уравнения
В случае параболического закона дисперсии ( p = 2), который применяется к свободным электронам в ферми-газе, результирующая плотность состояний для электронов в n-мерных системах равна
для , с для .
В одномерных системах DOS расходится в нижней части полосы и падает до . В двумерных системах DOS оказывается независимой от . Наконец, для трехмерных систем DOS возрастает как квадратный корень из энергии. [4]
С учетом префактора выражение для 3D DOS имеет вид
где – полный объем, включая 2-кратное спиновое вырождение.
В случае линейной зависимости ( p = 1), например, применимой к фотонам , акустическим фононам или некоторым особым видам электронных зон в твердом теле, DOS в 1, 2 и 3-мерных системах связана с энергией как :
Плотность состояний играет важную роль в кинетической теории твердого тела . Произведение плотности состояний и функции распределения вероятностей представляет собой количество занятых состояний в единице объема при заданной энергии для системы, находящейся в тепловом равновесии. Эта величина широко используется для исследования различных физических свойств материи. Ниже приведены примеры с использованием двух общих функций распределения того, как применение функции распределения к плотности состояний может привести к появлению физических свойств.
Статистика Ферми-Дирака : Функция распределения вероятностей Ферми-Дирака, рис. 4, используется для определения вероятности того, что фермион займет определенное квантовое состояние в системе, находящейся в тепловом равновесии. Фермионы — это частицы, подчиняющиеся принципу запрета Паули (например, электроны, протоны, нейтроны). Функцию распределения можно записать как
– химический потенциал (также обозначаемый как EF и называемый уровнем Ферми , когда T = 0), – постоянная Больцмана и – температура. На рис. 4 показано, как произведение функции распределения Ферми-Дирака и трехмерной плотности состояний полупроводника может дать представление о физических свойствах, таких как концентрация носителей заряда и ширина запрещенной зоны.
Статистика Бозе-Эйнштейна : функция распределения вероятностей Бозе-Эйнштейна используется для определения вероятности того, что бозон займет определенное квантовое состояние в системе, находящейся в тепловом равновесии. Бозоны — это частицы, которые не подчиняются принципу Паули (например, фононы и фотоны). Функцию распределения можно записать как
Из этих двух распределений можно рассчитать такие свойства, как внутренняя энергия единицы объема , число частиц , удельная теплоемкость и теплопроводность . Соотношения между этими свойствами и произведением плотности состояний и распределения вероятностей, обозначающего плотность состояний через вместо , имеют вид
— размерность, — скорость звука и — средняя длина свободного пробега .
Плотность состояний появляется во многих областях физики и помогает объяснить ряд квантово-механических явлений.
Расчет плотности состояний для малых структур показывает, что распределение электронов меняется с уменьшением размерности. Для квантовых проводов DOS для определенных энергий фактически становится выше, чем DOS для объемных полупроводников, а для квантовых точек электроны квантоваются до определенных энергий.
Плотностью состояний фотонов можно манипулировать, используя периодические структуры с масштабами длины порядка длины волны света. Некоторые структуры могут полностью подавлять распространение света определенных цветов (энергий), создавая фотонную запрещенную зону: DOS равна нулю для этих энергий фотонов. Другие структуры могут препятствовать распространению света только в определенных направлениях, создавая зеркала, волноводы и полости. Такие периодические структуры известны как фотонные кристаллы . [5] [6] [7] [8] В наноструктурированных средах концепция локальной плотности состояний (LDOS) часто более актуальна, чем концепция DOS, поскольку DOS значительно варьируется от точки к точке.
Интересные системы вообще сложны, например соединения, биомолекулы, полимеры и т. д. Из-за сложности этих систем аналитический расчет плотности состояний в большинстве случаев невозможен. Компьютерное моделирование предлагает набор алгоритмов для оценки плотности состояний с высокой точностью. Один из этих алгоритмов называется алгоритмом Ванга и Ландау . [9]
В рамках схемы Ванга и Ландау требуются любые предварительные знания о плотности состояний. Поступают следующим образом: функция стоимости (например, энергии) системы дискретизируется. Каждый раз, когда достигается интервал i , гистограмма плотности состояний обновляется на
Алгоритм Ванга и Ландау имеет некоторые преимущества перед другими распространенными алгоритмами, такими как мультиканоническое моделирование и параллельное регулирование . Например, плотность состояний получается как основной продукт моделирования. Кроме того, моделирование Ванга и Ландау полностью не зависит от температуры. Эта функция позволяет вычислять плотность состояний систем с очень грубым энергетическим ландшафтом, таких как белки. [10]
Математически плотность состояний выражается в виде башни покрывающих карт. [11]
Важной особенностью определения DOS является то, что его можно распространить на любую систему. Одним из его свойств является трансляционная инвариантность, означающая, что плотность состояний однородна и одинакова в каждой точке системы. Но это всего лишь частный случай, и LDOS дает более широкое описание с неоднородной плотностью состояний в системе.
Локальная плотность состояний (LDOS) описывает плотность состояний с пространственным разрешением. В материаловедении, например, этот термин полезен при интерпретации данных сканирующего туннельного микроскопа (СТМ), поскольку этот метод способен отображать электронные плотности состояний с атомным разрешением. В соответствии с кристаллической структурой эту величину можно предсказать с помощью вычислительных методов, например, с помощью теории функционала плотности .
В локальной плотности состояний вклад каждого состояния взвешивается плотностью его волновой функции в точке. становится
Фактор означает, что каждый штат вносит больший вклад в регионы с высокой плотностью населения. Усреднение этого выражения восстановит обычную формулу для DOS. LDOS полезен в неоднородных системах, где содержится больше информации, чем в одиночку.
Для одномерной системы со стенкой синусоидальные волны дают
где .
В трехмерной системе с выражением
Фактически, мы можем обобщить локальную плотность состояний дальше:
это называется спектральной функцией , и это функция, в которой каждая волновая функция находится отдельно в своей переменной. В более продвинутой теории это связано с функциями Грина и обеспечивает компактное представление некоторых результатов, таких как оптическое поглощение .
LDOS можно использовать для получения прибыли от твердотельного устройства. Например, рисунок справа иллюстрирует LDOS транзистора при его включении и выключении в баллистическом моделировании. LDOS имеет четкую границу в истоке и стоке, что соответствует положению края зоны. В канале DOS увеличивается по мере увеличения напряжения на затворе и снижения потенциального барьера.
В оптике и фотонике понятие локальной плотности состояний относится к состояниям, которые может занять фотон. Для света его обычно измеряют методами флуоресценции, методами сканирования ближнего поля или методами катодолюминесценции. Различные фотонные структуры имеют разное поведение LDOS с разными последствиями для спонтанного излучения. В фотонных кристаллах ожидается околонулевая LDOS, подавляющая спонтанное излучение. [12] Аналогичное усиление LDOS ожидается и в плазмонной полости. [13] Однако в неупорядоченных фотонных наноструктурах LDOS ведут себя иначе. Они колеблются в пространстве в зависимости от своей статистики и пропорциональны силе рассеяния структур. [14] Кроме того, связь со средней длиной свободного пробега рассеяния тривиальна, поскольку на LDOS все еще могут сильно влиять короткие детали сильных нарушений в форме сильного Парселловского усиления излучения. [15] и, наконец, для плазмонного беспорядка этот эффект гораздо сильнее для флуктуаций LDOS, поскольку его можно наблюдать как сильную ближнепольную локализацию. [16]