Твердое тело можно рассматривать как очень большую молекулу, в которой соблюдается либо дальний (идеальный), либо ближний порядок расположения атомов или молекул. Первый случай отвечает идеально регулярным кристаллам, симметрия которых описывается 230 федоровскими пространственными группами. Второй, более общий случай, охватывает реальные кристаллы с дефектами кристаллической решетки, аморфные тела, полимеры и т.д.: для них характерна лишь локальная симметрия определенных фрагментов. Для каждого из упомянутых случаев существуют свои способы описания электронной структуры. Идеально регулярные кристаллы считаются бесконечными (наличие поверхности игнорируется) и их волновые функции определяются с явным учетом трансляционной симметрии кристалла. В кристаллах с ближним порядком часто достаточно использовать модифицированные молекулярные модели, выделив некоторый атомный или молекулярный фрагмент – кластер. Рассмотрим оба этих подхода подробнее.
Одноэлектронные волновые функции в бесконечных периодических кристаллах и методы их расчета
Симметрийные свойства идеальных
кристаллов состоят в следующем. В идеальном
кристалле всегда можно ввести три вектора трансляций: ,
и
так, что физические свойства кристалла в некоторой произвольно выбранной
точке
точно воспроизводятся
в любой другой точке
удовлетворяющей
условию:
Параллелепипед, имеющий в
качестве ребер вектора ,
и
,
называется элементарной ячейкой (рис.5.1).
Перемещение всего кристалла как целого параллельно самому себе, описываемое
вектором
,называется
трансляцией. Вектор трансляции кристаллической решетки связывает
любые две точки ре-шетки.
(5.2)
Посредством
соответствующих операций трансляций элементарной ячейкой можно заполнить
все пространство кристаллической структуры. Такое свойство кристалла названо
трансляционной
симметрией. На рис.5.2 представлены
структуры некоторых атомных и молекулярных кристаллов.
Рис.5.2. Кристаллические структуры NaCl(а), алмаза (б), нафталина и I2 (г).
Трансляционная симметрия предполагает бесконечную протяженность кристалла. Конечно, регулярные структуры не являются бесконечными, а при отсутствии бесконечности теряется важное свойство трансляционной симметрии.
Один из способов сохранения трансляционной симметрии конечных систем - наложение циклических граничных условий (условия Борна-фон-Кармана). Суть их в том, что эквивалентные группы атомов отождествляются друг с другом так, чтобы граница отсутствовала.
Важнейшие электронные свойства периодических кристаллов хорошо видны уже в случае одномерных систем. Для цепочки атомов это равносильно изгибу цепочки в окружность большого радиуса и соединение между собой концевых атомов. Рассмотрим цепочку, состоящую из атомов, разделенных расстоянием а и имеющих один валентный электрон (Li, Na, K и т.д.), описываемый атомной орбиталью c (х). Каждое ядро в этих атомах хорошо экранировано внутренними электронами и валентный электрон лишь слабо связан со “своим” остовом. Влияние соседних атомом в цепочке приводит к тому, что валентные электроны могут передвигаться от атома к атому, как если бы они были свободными (модель почти свободных электронов Зоммерфельда). Волновая функция свободного электрона есть бегущая волна exp(ikx), а его кинетическая энергия равна
u0(x)=N [c ( x )+c ( x + a)+ c (x + 2a)+ c (x + 3a) +…], (5.6)
а при b=1 ®k= ± p /a и
up/a (x)=N [c ( x )-c ( x + a)+ c (x + 2a)- c (x + 3a) +…] (5.7)
(см. рис5.3). В первом случае АО на разных атомах интерферируют в фазе и результирующая функция Блоха является связывающей и отвечает более низкой энергии. Второй случай отвечает деструктивной интерференции АО и,
a) k = 0 u0 = N [ c (x) + c (x + a) + c ....]
б) k = + p/a up/a = N [ c (a) - c (x + a) + c .....]
Рис.5.3.
соответственно, - разрыхляющей функции Блоха. Таким образом, волновые функции электрона в периодическом потенциале можно ассоциировать со стоячими волнами вида u1=sinpx/a и u2 =cospx/a, каждая из которых есть сумма бегущих волн (рис.5.4). Эти стоячие волны имеют узлы и пучности в разных областях пространства по отношению к создающим внешний потенциал атомным остовам. Это значит, что они отвечают концентрации электронов (которая пропорциональна u2) в местах, отвечающих различным значениям потенциальной энергии. Так, u2 описывает концентрации электронов вблизи остовов, где потенциальная энергия минимальна, а u1 - между остовами.
Число
возможных значений k определяется
числом узлов кристаллической решетки. Область в k
–пространстве, лежащая между ±p /a,
называется
первой зоной Бриллюэна. Поведение электронов в более
высоких зонах Бриллюэна
Рис.5.4.
получено из информации о первой зоне (рис.5.5). В пределах первой зоны Бриллюэна энергия по-прежнему имеет вид, близкий к (2). Однако на границах зон Бриллюэна среди возможных (“разрешенных”) значений энергии появляются разрывы из-за того, что средняя потенциальная энергия электрона, описываемого функцией u2 меньше, чем для случая u1. Разность энергий Eg называется энергетической щелью.
Кристаллические орбитали – аналоги МО – строятся из функций Блоха (4), которые таким образом играют в кристалле роль базисных функций:
Распределение энергетических уровней в зоне характеризует плотность состояний:
Разность между высшим и низшим энергетическими уровнями называют шириной зоны (полосы), она также характеризует дисперсию энергии. Ширина зоны определяется величиной взаимодействия между атомами в соседних ячейках и зависит от соответствующих интегралов перекрывания АО.
Промежуток
между верхней энергетической зоной, заполненной элект-ронами (валентной),
и нижней, не заполненной электронами (зоной проводи-мости) называется запрещенной
зоной.
атомов
с g валентными электронами (g=1 для s оболочки; g=3 для p оболочки и т.д.),
то зона состоит из Ng энергетических уровней, каждый из которых может быть
занят 1 или 2 электронами (рис5.6). Если валентная зона полностью запол-нена
электронами, а ширина запрещенной зоны Еg?2
эВ, то вещество называется диэлектриком. Если валентная зона полностью
заполнена электронами, а ширина запрещенной зоны Еg?2
эВ, то вещество называется полупроводником. Если валентная зона заполнена
электронами частично, либо перекрывается с зоной проводимости, то вещество
называется металлом.
здесь VN –потенциал ядер, Vxc – обменно-корреляционный потенциал, для которого используются различные приближения – функционалы ЭП. Базисные функции различны: это могут быть как комбинации атомных орбиталей, так и наборы плоских волн. Недостаток: недооценка ширины запрещенной зоны.
Кластерные методы
Часто (например, при расчетах нестехиометрических кристаллов) целесообразно рассчитать лишь некоторый атомный или молекулярный фрагмент кристалла – кластер. Основной метод, используемый для это цели– метод Хартри-Фока. Одна из распространенных его кластерных модификаций – так называемый Xa -метод. Корреляция электронов здесь по-прежнему игнорируется, а обменный потенциал представляется в виде: