Жесткая привязка

редактировать
Модель электронных зонных структур твердых тел

В физике твердого тела, модель сильной связи (или модель TB ) - это подход к расчету электронной зонной структуры с использованием приближенного набора волновых функций на основе при суперпозиции волновых функций для изолированных атомов, расположенных в каждом атомном узле. Метод тесно связан с методом ЛКАО (метод линейной комбинации атомных орбиталей), используемым в химии. Модели с сильным связыванием применяются к широкому спектру твердых тел. Модель дает хорошие качественные результаты во многих случаях и может быть объединена с другими моделями, которые дают лучшие результаты там, где модель сильной привязки терпит неудачу. Хотя модель сильной связи является одноэлектронной, она также обеспечивает основу для более сложных расчетов, таких как вычисление состояний поверхности и применение к различным видам задач многих тел и квазичастиц вычисления.

Содержание
  • 1 Введение
  • 2 Историческая справка
  • 3 Математическая формулировка
    • 3.1 Трансляционная симметрия и нормализация
    • 3.2 Гамильтониан жесткой привязки
    • 3.3 Матричные элементы жесткой привязки
  • 4 Оценка элементов матрицы
  • 5 Связь с функциями Ванье
  • 6 Второе квантование
  • 7 Пример: одномерный s-диапазон
  • 8 Таблица межатомных матричных элементов
  • 9 См. Также
  • 10 Ссылки
  • 11 Дополнительная литература
  • 12 Внешние ссылки
Введение

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

Хотя математическая формулировка одночастичного гамильтониана гамильтониана может показаться сложной на первый взгляд, модель совсем не сложна и может быть легко понята интуитивно. Есть только три вида матричных элементов, которые играют важную роль в теории. Два из этих трех видов элементов должны быть близки к нулю, и ими часто можно пренебречь. Наиболее важными элементами в модели являются элементы межатомной матрицы, которые химик просто назвал бы энергиями связи.

Как правило, в модели участвует ряд атомных энергетических уровней и атомных орбиталей. Это может привести к сложным ленточным структурам, поскольку орбитали принадлежат разным представлениям точечной группы. обратная решетка и зона Бриллюэна часто принадлежат к другой пространственной группе, чем кристалл твердого тела. Точки высокой симметрии в зоне Бриллюэна принадлежат разным представлениям точечных групп. Когда изучаются простые системы, такие как решетки элементов или простые соединения, часто не очень сложно вычислить собственные состояния в точках с высокой симметрией аналитически. Таким образом, модель сильной привязки может служить хорошими примерами для тех, кто хочет узнать больше о теории групп.

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

Историческая справка

К 1928 году идею молекулярной орбитали выдвинул Роберт Малликен, на которого значительное влияние оказали работы Фридриха Хунда. Метод ЛКАО для аппроксимации молекулярных орбиталей был введен в 1928 г. Б. Н. Финклештейном и Г. Е. Горовицем, в то время как метод ЛКАО для твердых тел был разработан Феликсом Блохом в рамках его докторской диссертации в 1928 г. одновременно и независимо от него. подход ЛКАО-МО. Намного более простая схема интерполяции для аппроксимации электронной зонной структуры, особенно для d-зон переходных металлов, - это параметризованный метод сильной связи, разработанный в 1954 году Джоном Кларком Слейтером и Джорджем Фред Костер, иногда упоминаемый как метод плотного связывания SK. При использовании метода сильной связи SK расчеты электронной зонной структуры твердого тела не обязательно должны выполняться с полной строгостью, как в исходной теореме Блоха, а, скорее, расчеты из первых принципов выполняются только при высоких значениях точки симметрии и зонная структура интерполируются по оставшейся части зоны Бриллюэна между этими точками.

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

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

Модель сильной связи обычно используется для расчетов электронной зонной структуры и запрещенной зоны в статическом режиме. Однако в сочетании с другими методами, такими как модель приближения случайной фазы (RPA), также можно изучать динамический отклик систем.

Математическая формулировка

Мы вводим атомные орбитали φ m (r) {\ displaystyle \ varphi _ {m} ({\ boldsymbol {r}})}\ varphi_m (\ boldsymbol { r}) , которые являются собственными функциями гамильтониана H в {\ displaystyle H _ {\ rm {at}}}{\ displaystyle H _ {\ rm {at}}} одиночного изолированного атома. Когда атом помещен в кристалл, эта атомная волновая функция перекрывает соседние атомные узлы, и поэтому не являются истинными собственными функциями гамильтониана кристалла. Перекрытие меньше, когда электроны тесно связаны, что является источником дескриптора «сильная связь». Любые поправки к атомному потенциалу Δ U {\ displaystyle \ Delta U}\ Delta U , необходимые для получения истинного гамильтониана H {\ displaystyle H}H системы, являются предполагается малым:

H (r) = = R n H at (r - R n) + Δ U (r), {\ displaystyle H ({\ boldsymbol {r}}) = \ sum _ {\ boldsymbol { R_ {n}}} H _ {\ mathrm {at}} ({\ boldsymbol {r-R_ {n}}}) + \ Delta U ({\ boldsymbol {r}}) \,}{\ displaystyle H ({\ boldsymbol {r}}) = \ sum _ {\ boldsymbol {R_ {n}}} H _ {\ mathrm {at}} ({\ boldsymbol { r-R_ {n}}}) + \ Delta U ({\ boldsymbol {r}}) \,}

где R n {\ displaystyle {\ boldsymbol {R}} _ {n}}{\ displaystyle {\ boldsymbol {R}} _ {n}} определяет местонахождение атомного узла в кристаллической решетке. Решение ψ m {\ displaystyle \ psi _ {m}}{\ displaystyle \ psi _ {m}} независимого от времени одиночного электрона уравнения Шредингера затем аппроксимируется как линейная комбинация атомных орбитали φ m (r - R n) {\ displaystyle \ varphi _ {m} ({\ boldsymbol {r-R_ {n}}})}{\ displaystyle \ varphi _ {m} ({\ boldsymbol {r-R_ {n}}})} :

ψ m (r) = ∑ R nbm (р N) φ м (г - р n) {\ displaystyle \ psi _ {m} ({\ boldsymbol {r}}) = \ sum _ {\ boldsymbol {R_ {n}}} b_ {m} ( {\ boldsymbol {R_ {n}}}) \ \ varphi _ {m} ({\ boldsymbol {r-R_ {n}}})}{\ displaystyle \ psi _ {m} ({\ boldsymbol {r}}) = \ sum _ {\ boldsymbol {R_ {n}}} b_ {m} ({\ boldsymbol {R_ {n}}}) \ \ varphi _ {m} ({\ boldsymbol {r-R_ {n}}})} ,

где m {\ displaystyle m}mотносится к m-му уровню атомной энергии.

Трансляционная симметрия и нормализация

Теорема Блоха утверждает, что волновая функция в кристалле может изменяться при трансляции только на фазовый множитель:

ψ (r + Р ℓ) знак равно eik ⋅ р ℓ ψ (r), {\ displaystyle \ psi ({\ boldsymbol {r + R _ {\ ell}}}) = e ^ {i {\ boldsymbol {k \ cdot R _ {\ ell} }}} \ psi ({\ boldsymbol {r}}) \,}\ psi (\ boldsymbol {r + R _ {\ ell}}) = e ^ {i \ boldsymbol {k \ cdot R _ {\ ell}}} \ psi (\ boldsymbol {r}) \,

где k {\ displaystyle \ mathbf {k}}\ mathbf {k} - волновой вектор волновой функции. Следовательно, коэффициенты удовлетворяют условию

∑ R n b m (R n) φ m (r - R n + R ℓ) = e i k ⋅ R ℓ ∑ R n b m (R n) φ m (r - R n). {\ displaystyle \ sum _ {\ boldsymbol {R_ {n}}} b_ {m} ({\ boldsymbol {R_ {n}}}) \ \ varphi _ {m} ({\ boldsymbol {r-R_ {n}) + R _ {\ ell}}}) = e ^ {i {\ boldsymbol {k \ cdot R _ {\ ell}}}} \ sum _ {\ boldsymbol {R_ {n}}} b_ {m} ({\ boldsymbol {R_ {n}}}) \ \ varphi _ {m} ({\ boldsymbol {r-R_ {n}}}) \.}{\ displaystyle \ sum _ {\ boldsymbol {R_ {n}}} b_ {m} ({\ boldsymbol {R_ {n}}}) \ \ varphi _ {m} ({\ boldsymbol {r-R_ {n} + R _ {\ ell}}}) = e ^ {i {\ boldsymbol {k \ cdot R _ {\ ell}}}} \ sum _ {\ boldsymbol {R_ {n }}} b_ {m} ({\ boldsymbol {R_ {n}}}) \ \ varphi _ {m} ({\ boldsymbol {r-R_ {n}}}) \.}

Подставляя R p = R n - R ℓ {\ displaystyle {\ boldsymbol {R_ {p}}} = {\ boldsymbol {R_ {n}}} - {\ boldsymbol {R _ {\ ell}}}}\ boldsymbol {R_p} = \ boldsymbol {R_n} - \ boldsymbol {R_ \ ell} , находим

bm (R п + р ℓ) знак равно eik ⋅ р ℓ bm (р п), {\ displaystyle b_ {m} ({\ boldsymbol {R_ {p} + R _ {\ ell}}}) = e ^ {i {\ boldsymbol { k \ cdot R _ {\ ell}}}} b_ {m} ({\ boldsymbol {R_ {p}}}) \,}{\ displaystyle b_ {m} ({\ boldsymbol {R_ {p} + R _ {\ ell}}}) = e ^ {i {\ boldsymbol {k \ cdot R _ {\ ell}}}} b_ {m } ({\ boldsymbol {R_ {p}}}) \,} (где в RHS мы заменили фиктивный индекс R п {\ displaystyle {\ boldsymbol {R_ {n}}}}\boldsymbol{R_n}с R p {\ displaystyle {\ boldsymbol {R_ {p}}}}\ boldsymbol {R_p} )

или

bm ( R l) = eik ⋅ R lbm (0). {\ Displaystyle b_ {m} ({\ boldsymbol {R_ {l}}}) = e ^ {i {\ boldsymbol {k \ cdot R_ {l}}}} b_ {m} ({\ boldsymbol {0}}) \.}{\ displaystyle b_ {m} ({\ boldsymbol {R_ {l}}}) = e ^ {i {\ boldsymbol {k \ cdot R_ {l}}}} b_ {m} ({\ boldsymbol {0}}) \.}

Нормализация волновой функции к единице:

∫ d 3 r ψ m ∗ (r) ψ m (r) = 1 { \ отображает tyle \ int d ^ {3} r \ \ psi _ {m} ^ {*} ({\ boldsymbol {r}}) \ psi _ {m} ({\ boldsymbol {r}}) = 1}{\ displaystyle \ int d ^ {3} r \ \ psi _ {m} ^ {*} ({\ boldsymbol {r}}) \ psi _ {m} ({\ boldsymbol {r}}) = 1}
Знак равно ∑ р nbm * (р n) ∑ р ℓ bm (р ℓ) ∫ d 3 r φ m * (r - R n) φ m (r - R ℓ) {\ displaystyle = \ sum _ {\ boldsymbol {R_ {n}}} b_ {m} ^ {*} ({\ boldsymbol {R_ {n}}}) \ sum _ {\ boldsymbol {R _ {\ ell}}} b_ {m} ({\ boldsymbol {R_ { \ ell}}}) \ int d ^ {3} r \ \ varphi _ {m} ^ {*} ({\ boldsymbol {r-R_ {n}}}) \ varphi _ {m} ({\ boldsymbol { r-R _ {\ ell}}})}{\ displaystyle = \ sum _ {\ boldsymbol {R_ {n}}} b_ {m} ^ {*} ({\ boldsymbol {R_ {n}}}) \ sum _ {\ boldsymbol {R _ {\ ell}}} b_ {m} ({\ boldsymbol {R _ {\ ell}}}) \ int d ^ {3} r \ \ varphi _ {m} ^ {*} ({\ boldsymbol {r-R_ {n}}}) \ varphi _ {m} ({\ boldsymbol {r-R _ {\ ell}}})}
= bm ∗ (0) bm (0) ∑ R ne - ik ⋅ R n ∑ R ℓ eik ⋅ R ℓ ∫ d 3 r φ m ∗ (r - R п) φ м (г - р ℓ) {\ displaystyle = b_ {m} ^ {*} (0) b_ {m} (0) \ sum _ {\ boldsymbol {R_ {n}}} e ^ {- i {\ boldsymbol {k \ cdot R_ {n}}}} \ sum _ {\ boldsymbol {R _ {\ ell}}} e ^ {i {\ boldsymbol {k \ cdot R _ {\ ell}}}} \ \ int d ^ {3} r \ \ varphi _ {m} ^ {*} ({\ boldsymbol {r-R_ {n}}}) \ varphi _ {m} ({\ boldsymbol {r-R _ {\ ell}}) })}{\ displaystyle = b_ {m} ^ {*} (0) b_ {m} (0) \ sum _ {\ boldsymbol {R_ {n}}} e ^ {- i {\ boldsymbol {k \ cdot R_ {n}}}} \ sum _ {\ boldsymbol {R_ {\ ell}}} e ^ {i {\ boldsymbol {k \ cdot R _ {\ ell}}}} \ \ int d ^ {3} r \ \ varphi _ {m} ^ {*} ({\ boldsymbol { r-R_ {n}}}) \ varphi _ {m} ({\ boldsymbol {r-R _ {\ ell}}})}
= N bm ∗ (0) bm (0) ∑ R pe - ik ⋅ R p ∫ d 3 r φ m ∗ (r - R p) φ m (r) {\ displaystyle = Nb_ {m } ^ {*} (0) b_ {m} (0) \ sum _ {\ boldsymbol {R_ {p}}} e ^ {- i {\ boldsymbol {k \ cdot R_ {p}}}} \ \ int d ^ {3} r \ \ varphi _ {m} ^ {*} ({\ bold символ {r-R_ {p}}}) \ varphi _ {m} ({\ boldsymbol {r}}) \}{\ displaystyle = Nb_ {m} ^ {*} (0) b_ {m} (0) \ sum _ {\ boldsymbol {R_ {p}}} e ^ {- i {\ boldsymbol {k \ cdot R_ {p}}}} \ \ int d ^ {3} r \ \ varphi _ {m} ^ {*} ({\ boldsymbol {r-R_ {p}}}) \ varphi _ {m} ( {\ boldsymbol {r}}) \}
= N bm ∗ (0) bm (0) ∑ R peik ⋅ R p ∫ d 3 р φ м * (г) φ м (г - р п), {\ displaystyle = Nb_ {m} ^ {*} (0) b_ {m} (0) \ sum _ {\ boldsymbol {R_ {p} }} e ^ {i {\ boldsymbol {k \ cdot R_ {p}}}} \ \ int d ^ {3} r \ \ varphi _ {m} ^ {*} ({\ boldsymbol {r}}) \ varphi _ {m} ({\ boldsymbol {r-R_ {p}}}) \,}{\ displaystyle = Nb_ {m} ^ {*} (0) b_ {m} (0) \ sum _ {\ boldsymbol {R_ {p}}} e ^ {i {\ boldsymbol { k \ cdot R_ {p}}}} \ \ int d ^ {3} r \ \ varphi _ {m} ^ {*} ({\ boldsymbol {r}}) \ varphi _ {m} ({\ boldsymbol { r-R_ {p}}}) \,}

, поэтому нормализация устанавливает bm (0) {\ displaystyle b_ {m} (0)}{\ displaystyle b_ {m} (0)} как

bm ∗ (0) bm (0) = 1 N ⋅ 1 1 + ∑ R p ≠ 0 eik ⋅ R p α m (R p), {\ displaystyle b_ {m} ^ {*} (0) b_ {m} (0) = {\ frac {1} {N}} \ \ cdot \ {\ frac {1} {1+ \ sum _ {\ boldsymbol {R_ {p} \ neq 0}} e ^ {i {\ boldsymbol {k \ cdot R_ {p}}}} \ alpha _ {m} ({\ boldsymbol {R_ {p}}})}} \,}{\ displaystyle b_ {m} ^ {*} (0) b_ {m} (0) = {\ frac {1} {N} } \ \ cdot \ {\ frac {1} {1+ \ sum _ {\ boldsymbol {R_ {p} \ neq 0}} e ^ {i {\ boldsymbol {k \ cdot R_ {p}}}} \ alpha _ {m} ({\ boldsymbol {R_ {p}}})}} \,}

где α m(Rp) являются атомными интегралами перекрытия, которыми часто пренебрегают, что приводит к

bm (0) ≈ 1 N, {\ displaystyle b_ {m} (0) \ приблизительно {\ frac {1} {\ sqrt {N}}} \,}{\ displaystyle b_ {m} (0) \ приблизительно {\ frac {1} {\ sqrt {N}}} \,}

и

ψ m (r) ≈ 1 N ∑ R neik ⋅ R n φ m (r - R n). {\ displaystyle \ psi _ {m} ({\ boldsymbol {r}}) \ приблизительно {\ frac {1} {\ sqrt {N}}} \ sum _ {\ boldsymbol {R_ {n}}} e ^ { i {\ boldsymbol {k \ cdot R_ {n}}}} \ \ varphi _ {m} ({\ boldsymbol {r-R_ {n}}}) \.}{\ displaystyle \ psi _ {m} ({\ boldsymbol {r}}) \ приблизительно {\ frac {1} {\ sqrt {N}}} \ sum _ {\ boldsymbol {R_ {n}}} e ^ {i {\ boldsymbol {k \ cdot R_ {n}}}} \ \ varphi _ {m} ({\ boldsymbol {r-R_ {n}}}) \.}

Гамильтониан жесткой связи

Используя форму сильной связи для волновой функции и предполагая, что для m-го энергетического диапазона важен только m-й атомный энергетический уровень, блоховские энергии ε m {\ displaystyle \ varepsilon _ {m}}\ varepsilon_m имеют вид

ε m = ∫ d 3 r ψ m ∗ (r) H (r) ψ (r) {\ displaystyle \ varepsilon _ {m} = \ int d ^ {3} r \ \ psi _ {m} ^ {*} ({\ boldsymbol {r}}) H ({\ boldsymbol {r}}) \ psi ({\ boldsymbol {r}})}{\ displaystyle \ varepsilon _ {m} = \ int d ^ {3} r \ \ psi _ {m} ^ {*} ({\ boldsymbol {r}}) H ({\ boldsymbol {r}}) \ psi ({\ boldsymbol {r}}) }
знак равно ∑ р nb * (р n) ∫ d 3 р φ * (r - р n) H (r) ψ (r) {\ displaystyle = \ sum _ {\ boldsymbol {R_ {n}}} b ^ {*} ({\ boldsymbol {R_ {n}}}) \ \ int d ^ {3} r \ \ varphi ^ {*} ({\ boldsymbol {r-R_ {n}}}) H ({\ boldsymbol {r}}) \ psi ({\ boldsymbol {r}}) \}{\ displaystyle = \ sum _ {\ boldsymbol {R_ {n}}} b ^ {*} ( {\ boldsymbol {R_ {n}}}) \ \ int d ^ {3} r \ \ varphi ^ {*} ({\ boldsymbol {r-R_ {n}}}) H ({\ boldsymbol {r}}) \ psi ({\ boldsymbol {r}}) \}
= ∑ R ℓ ∑ R nb ∗ (R n) ∫ d 3 r φ ∗ (r - R n) H at (r - R ℓ) ψ (r) + ∑ R nb ∗ (R n) ∫ d 3 r φ ∗ (r - R n) Δ U (r) ψ ( р). {\ displaystyle = \ sum _ {\ boldsymbol {R _ {\ ell}}} \ \ sum _ {\ boldsymbol {R_ {n}}} b ^ {*} ({\ boldsymbol {R_ {n}}}) \ \ int d ^ {3} r \ \ varphi ^ {*} ({\ boldsymbol {r-R_ {n}}}) H _ {\ mathrm {at}} ({\ boldsymbol {r-R _ {\ ell}} }) \ psi ({\ boldsymbol {r}}) \ + \ sum _ {\ boldsymbol {R_ {n}}} b ^ {*} ({\ boldsymbol {R_ {n}}}) \ \ int d ^ {3} r \ \ varphi ^ {*} ({\ boldsymbol {r-R_ {n}}}) \ Delta U ({\ boldsymbol {r}}) \ psi ({\ boldsymbol {r}}) \. }{\ displaystyle = \ sum _ {\ boldsymbol {R _ {\ ell}}} \ \ sum _ {\ boldsymbol {R_ {n}}} b ^ {*} ({\ boldsymbol {R_ {n}}}) \ \ int d ^ {3} r \ \ varphi ^ {*} ({\ boldsymbol {r-R_ {n}}}) H _ {\ mathrm {at}} ({\ boldsymbol {r-R _ {\ ell}}}) \ psi ({\ boldsymbol {r}}) \ + \ sum _ {\ boldsymbol {R_ {n}}} b ^ {*} ({\ boldsymbol {R_ {n}}}) \ \ int d ^ {3 } r \ \ varphi ^ {*} ({\ boldsymbol {r-R_ {n}}}) \ Delta U ({\ boldsymbol {r}}) \ psi ({\ boldsymbol {r}}) \.}
≈ E m + b ∗ (0) ∑ R ne - ik ⋅ R n ∫ d 3 r φ ∗ (r - R n) Δ U (r) ψ (r). {\ displaystyle \ приблизительно E_ {m} + b ^ {*} (0) \ sum _ {\ boldsymbol {R_ {n}}} e ^ {- i {\ boldsymbol {k \ cdot R_ {n}}}} \ \ int d ^ {3} r \ \ varphi ^ {*} ({\ boldsymbol {r-R_ {n}}}) \ Delta U ({\ boldsymbol {r}}) \ psi ({\ boldsymbol {r }}) \.}\ приблизительно E_m + b ^ * (0) \ sum _ {\ boldsymbol {R_n}} e ^ {- i \ boldsymbol {k \ cdot R_n} } \ \ int d ^ 3 r \ \ varphi ^ * (\ boldsymbol {r-R_n}) \ Delta U (\ boldsymbol {r}) \ psi (\ boldsymbol {r}) \.

Здесь не учитываются члены, включающие атомный гамильтониан в узлах, отличных от его центра. Тогда энергия становится

ε m (k) = E m - N | b (0) | 2 (β м + ∑ р N ≠ 0 ∑ l γ м, l (R n) eik ⋅ R n), {\ displaystyle \ varepsilon _ {m} ({\ boldsymbol {k}}) = E_ {m} - N \ | b (0) | ^ {2} \ left (\ beta _ {m} + \ sum _ {{\ boldsymbol {R_ {n}}} \ neq 0} \ sum _ {l} \ gamma _ { m, l} ({\ boldsymbol {R_ {n}}}) e ^ {i {\ boldsymbol {k}} \ cdot {\ boldsymbol {R_ {n}}}} \ right) \,}\ varepsilon_m (\ boldsymbol {k}) = E_m - N \ | b (0) | ^ 2 \ left (\ beta_m + \ sum _ {\ boldsymbol {R_n} \ neq 0} \ sum_l \ gamma_ {m, l} (\ boldsymbol {R_n}) e ^ {i \ boldsymbol {k} \ cdot \ boldsymbol {R_n}} \ right) \,
= Е м - β м + ∑ р N ≠ 0 ∑ leik ⋅ R n γ m, l (R n) 1 + ∑ R n ≠ 0 ∑ leik ⋅ R n α m, l (R n), {\ displaystyle = E_ {m} - \ {\ frac {\ beta _ {m} + \ sum _ {{\ boldsymbol {R_ {n}}} \ neq 0} \ sum _ {l} e ^ {i {\ boldsymbol {k}) } \ cdot {\ boldsymbol {R_ {n}}}} \ gamma _ {m, l} ({\ boldsymbol {R_ {n}}})} {\ \ 1+ \ sum _ {\ boldsymbol {R_ {n } \ neq 0}} \ sum _ {l} e ^ {i {\ boldsymbol {k \ cdot R_ {n}}}} \ alpha _ {m, l} ({\ boldsymbol {R_ {n}}}) }} \,}= E_m - \ \ frac {\ beta_m + \ sum _ {\ boldsymbol {R_n} \ neq 0} \ sum_l e ^ {i \ boldsymbol {k} \ cdot \ boldsymbol {R_n}} \ gamma_ {m, l} (\ boldsymbo l {R_n})} {\ \ 1 + \ sum _ {\ boldsymbol {R_n \ neq 0}} \ sum_l e ^ {i \ boldsymbol {k \ cdot R_n}} \ alpha_ {m, l} (\ boldsymbol {R_n })} \,

где E m - энергия m-го атомного уровня, а α m, l {\ displaystyle \ alpha _ {m, l}}\ alpha_ {m, l} , β m {\ displaystyle \ beta _ {m}}\ beta_m и γ m, l {\ displaystyle \ gamma _ {m, l}}\ gamma_ {m, l} - матрица жесткой привязки элементы.

Матричные элементы жесткой привязки

Элемент

β m = - ∫ φ m ∗ (r) Δ U (r) φ m (r) d 3 r {\ displaystyle \ beta _ {m} = - \ int \ varphi _ {m} ^ {*} ({\ boldsymbol {r}}) \ Delta U ({\ boldsymbol {r}}) \ varphi _ {m} ({\ boldsymbol {r}}) \, d ^ {3} r \}\ beta_m = - \ int \ varphi_m ^ * (\ boldsymbol {r}) \ Delta U (\ boldsymbol {r}) \ varphi_m (\ boldsymbol {r}) \, d ^ 3 r \ ,

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

Следующий член

γ m, l (R n) = - ∫ φ m ∗ (r) Δ U (r) φ l (r - R n) d 3 r, {\ displaystyle \ гамма _ {m, l} ({\ boldsymbol {R_ {n}}}) = - \ int \ varphi _ {m} ^ {*} ({\ boldsymbol {r}}) \ Delta U ({\ boldsymbol { r}}) \ varphi _ {l} ({\ boldsymbol {r-R_ {n}}}) \, d ^ {3} r \,}\ gamma_ {m, l} (\ boldsymbol {R_n}) = - \ int \ varphi_m ^ * (\ boldsymbol {r}) \ Delta U (\ boldsymbol {r}) \ varphi_l (\ boldsymbol {r - R_n}) \, d ^ 3 r \,

- это элемент межатомной матрицы между атомные орбитали m и l на соседних атомах. Его также называют энергией связи или интегралом двух центров, и это наиболее важный элемент в модели сильной связи.

Последние члены

α m, l (R n) = ∫ φ m ∗ (r) φ l (r - R n) d 3 r {\ displaystyle \ alpha _ {m, l} ({\ boldsymbol {R_ {n}}}) = \ int \ varphi _ {m} ^ {*} ({\ boldsymbol {r}}) \ varphi _ {l} ({\ boldsymbol {r-R_ {n }}}) \, d ^ {3} r \}\ alpha_ {m, l} ( \ boldsymbol {R_n}) = \ int \ varphi_m ^ * (\ boldsymbol {r}) \ varphi_l (\ boldsymbol {r - R_n}) \, d ^ 3 r \ ,

обозначают интегралы перекрытия между атомными орбиталями m и l на соседних атомах.

Оценка элементов матрицы

Как упоминалось ранее, значения элементов матрицы β m {\ displaystyle \ beta _ {m}}\ beta_m не настолько велика по сравнению с энергией ионизации, потому что потенциалы соседних атомов на центральном атоме ограничены. Если β m {\ displaystyle \ beta _ {m}}\ beta_m не относительно мало, это означает, что потенциал соседнего атома на центральном атоме тоже не мал. В этом случае это признак того, что модель сильного связывания по какой-то причине не очень хорошая модель для описания структуры ленты. Например, межатомные расстояния могут быть слишком малы или заряды на атомах или ионах в решетке неправильные.

Межатомные матричные элементы γ m, l {\ displaystyle \ gamma _ {m, l}}\ gamma_ {m, l} могут быть вычислены напрямую, если атомные волновые функции и потенциалы известны в деталь. Чаще всего это не так. Есть множество способов получить параметры для этих матричных элементов. Параметры могут быть получены из данных по энергии химической связи. Энергии и собственные состояния в некоторых точках высокой симметрии в зоне Бриллюэна могут быть оценены, а интегралы значений в матричных элементах могут быть сопоставлены с данными зонной структуры из других источников.

Элементы матрицы межатомного перекрытия α m, l {\ displaystyle \ alpha _ {m, l}}\ alpha_ {m, l} должны быть довольно маленькими или пренебрежимо малыми. Если они большие, это снова указывает на то, что модель жесткой привязки имеет ограниченную ценность для некоторых целей. Например, большое перекрытие является признаком слишком короткого межатомного расстояния. В металлах и переходных металлах широкая s-зона или sp-зона может быть лучше приспособлена к существующему расчету зонной структуры путем введения матричных элементов следующего ближайшего соседа и интегралов перекрытия, но подобные подходы не дают очень полезной модели. для электронной волновой функции металла. Широкие зоны в плотных материалах лучше описываются моделью почти свободных электронов.

Модель сильной связи особенно хорошо работает в случаях, когда ширина зоны мала, а электроны сильно локализованы, как в случае d-зон и f-диапазоны. Модель также дает хорошие результаты в случае открытых кристаллических структур, таких как алмаз или кремний, где количество соседей невелико. Модель может быть легко объединена с моделью почти свободных электронов в гибридной модели NFE-TB.

Связь с функциями Ванье

Функции Блоха описывают электронные состояния в периодической кристаллической решетке. Функции Блоха могут быть представлены как ряд Фурье

ψ m (k, r) = 1 N ∑ nam (R n, r) eik ⋅ R n, {\ displaystyle \ psi _ {m} \ mathbf { (k, r)} = {\ frac {1} {\ sqrt {N}}} \ sum _ {n} {a_ {m} \ mathbf {(R_ {n}, r)}} e ^ {\ mathbf {ik \ cdot R_ {n}}} \,}\ psi_m \ mathbf {(k, r)} = \ frac { 1} {\ sqrt {N}} \ sum_ {n} {a_m \ mathbf {(R_n, r)}} e ^ {\ mathbf {ik \ cdot R_n}} \,

где Rnобозначает атомную позицию в периодической кристаллической решетке, k - волновой вектор блоховского По теореме r - положение электрона, m - индекс зоны, а сумма - по всем N атомным позициям. Теорема Блоха представляет собой точное собственное решение для волновой функции электрона в периодическом кристаллическом потенциале, соответствующем энергии E m(k), и распространяется по всему объему кристалла.

Используя анализ преобразования Фурье, пространственно локализованную волновую функцию для m-го энергетического диапазона можно построить из нескольких теорем Блоха:

am (R n, r) = 1 N ∑ ke - ik R n ψ m (k, r) = 1 N ∑ keik ⋅ (r - R n) um (k, r). {\ displaystyle a_ {m} \ mathbf {(R_ {n}, r)} = {\ frac {1} {\ sqrt {N}}} \ sum _ {\ mathbf {k}} {e ^ {\ mathbf {-ik \ cdot R_ {n}}} \ psi _ {m} \ mathbf {(k, r)}} = {\ frac {1} {\ sqrt {N}}} \ sum _ {\ mathbf {k }} {e ^ {\ mathbf {ik \ cdot (r-R_ {n})}} u_ {m} \ mathbf {(k, r)}}.}a_m \ mathbf {(R_n, r)} = \ frac {1} {\ sqrt {N}} \ sum _ {\ mathbf {k}} {e ^ {\ mathbf {-ik \ cdot R_n}} \ psi_m \ mathbf {(k, r)}} = \ frac {1} {\ sqrt {N}} \ sum _ {\ mathbf {k}} {e ^ {\ mathbf {ik \ cdot (r-R_n)} } u_m \ mathbf {(k, r)}}.

Эти волновые функции реального пространства am (R n, r) {\ displaystyle {a_ {m} \ mathbf {(R_ {n}, r)}}}{a_m \ mathbf {(R_n, r)}} называются функциями Ванье и довольно близко локализованы на атомный сайт Rn. Конечно, если у нас есть точные функции Ванье, точные функции Блоха могут быть получены с использованием обратного преобразования Фурье.

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

Второе квантование

Современные объяснения электронной структуры, такие как модель t-J и модель Хаббарда, основаны на модели сильной связи. Тесную привязку можно понять, работая в рамках формализма второго квантования.

Используя атомную орбиталь в качестве базового состояния, оператор гамильтониана второго квантования в структуре сильной связи можно записать как:

H = - t ∑ ⟨i, j⟩, σ (ci, σ † cj, σ + час. c.) {\ displaystyle H = -t \ sum _ {\ langle i, j \ rangle, \ sigma} (c_ {i, \ sigma} ^ {\ dagger} c_ {j, \ sigma } ^ {} + hc)}{ \ displaystyle H = -t \ sum _ {\ langle i, j \ rangle, \ sigma} (c_ {i, \ sigma} ^ {\ dagger} c_ {j, \ sigma} ^ {} + hc)} ,
ci σ †, cj σ {\ displaystyle c_ {i \ sigma} ^ {\ dagger}, c_ {j \ sigma}}c ^ \ dagger_ {i \ sigma}, c_ {j \ sigma} - операторы создания и уничтожения
σ {\ displaystyle \ displaystyle \ sigma}\ displaystyle \ sigma - спиновая поляризация
t {\ displaystyle \ displaystyle t}\ displaystyle t - интеграл перескока
⟨i, j⟩ { \ displaystyle \ displaystyle \ langle i, j \ rangle}\ displaystyle \ langle i, j \ rangle - индекс ближайшего соседа
h. c. {\ displaystyle \ displaystyle hc}{\ displaystyle \ displaystyle hc} - эрмитовское сопряжение другого члена (ов)

Здесь интеграл перескока t {\ displaystyle \ displaystyle t}\ displaystyle t соответствует интеграл переноса γ {\ displaystyle \ displaystyle \ gamma}\displaystyle\gammaв модели с сильной привязкой. Рассматривая крайние случаи t → 0 {\ displaystyle t \ rightarrow 0}t \ стрелка вправо 0 , электрон не может прыгать в соседние узлы. Этот случай - изолированная атомная система. Если параметр скачкообразного изменения включен (t>0 {\ displaystyle \ displaystyle t>0}\displaystyle t>0 ) электроны могут оставаться в обоих узлах, понижая их кинетическую энергию.

В сильно коррелированной электронной системе необходимо учитывать электрон-электронное взаимодействие. Этот член можно записать в виде

H ee = 1 2 ∑ n, m, σ ⟨n 1 m 1, n 2 m 2 | e 2 | r 1 - r 2 | | n 3 m 3, N 4 м 4⟩ сп 1 м 1 σ 1 † сп 2 м 2 σ 2 † сп 4 м 4 σ 2 сп 3 м 3 σ 1 {\ displaystyle \ displaystyle H_ {ee} = {\ frac {1} { 2}} \ sum _ {n, m, \ sigma} \ langle n_ {1} m_ {1}, n_ {2} m_ {2} | {\ frac {e ^ {2}} {| r_ {1} -r_ {2} |}} | n_ {3} m_ {3}, n_ {4} m_ {4} \ rangle c_ {n_ {1} m_ {1} \ sigma _ {1}} ^ {\ dagger} c_ {n_ {2} m_ {2} \ sigma _ {2}} ^ {\ dagger} c_ {n_ {4} m_ {4} \ sigma _ {2}} c_ {n_ {3} m_ {3} \ sigma _ {1}}}\ displaystyle H_ {ee} = \ frac {1} {2} \ sum_ {n, m, \ sigma} \ langle n_1 m_1, n_2 m_2 | \ frac {e ^ 2} {| r_1-r_2 |} | n_3 m_3, n_4 m_4 \ rangle c ^ \ dagger_ {n_ 1 m_1 \ sigma_1} c ^ \ dagger_ {n_2 m_2 \ sigma_2} c_ {n_4 m_4 \ sigma_2} c_ {n_3 m_3 \ sigma_1}

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

Пример: одномерный s-диапазон

Здесь модель сильной связи проиллюстрирована с помощью модели s-диапазона для цепочки атомов с одной s-орбиталью по прямой с интервалом а и σ-связями между атомными узлами.

Чтобы найти приблизительные собственные состояния гамильтониана, мы можем использовать линейную комбинацию атомных орбиталей

| k⟩ = 1 N ∑ n = 1 N e i n k a | n⟩ {\ displaystyle | k \ rangle = {\ frac {1} {\ sqrt {N}}} \ sum _ {n = 1} ^ {N} e ^ {inka} | n \ rangle}| k \ rangle = \ frac {1} {\ sqrt {N}} \ sum_ {n = 1 } ^ N e ^ {inka} | n \ rangle

где N = общее количество сайтов, а k {\ displaystyle k}k - реальный параметр с - π a ≦ k ≦ π a {\ displaystyle - {\ frac {\ pi} { a}} \ leqq k \ leqq {\ frac {\ pi} {a}}}- \ frac {\ pi} {a} \ leqq k \ leqq \ frac {\ pi} {a} . (Эта волновая функция нормирована к единице главным множителем 1 / √N при условии, что перекрытие атомных волновых функций игнорируется.) Если предположить, что перекрываются только ближайшие соседи, единственные ненулевые матричные элементы гамильтониана могут быть выражены как

⟨ п | H | п⟩ знак равно E 0 = E я - U. {\ displaystyle \ langle n | H | n \ rangle = E_ {0} = E_ {i} -U \.}\ langle n | H | n \ rangle = E_0 = E_i - U \.
⟨n ± 1 | H | n⟩ знак равно - Δ {\ displaystyle \ langle n \ pm 1 | H | n \ rangle = - \ Delta \}\ langle n \ pm 1 | H | n \ rangle = - \ Delta \
⟨n | n⟩ = 1; {\ displaystyle \ langle n | n \ rangle = 1 \;}\ langle n | n \ rangle = 1 \; ⟨n ± 1 | п⟩ = S. {\ displaystyle \ langle n \ pm 1 | n \ rangle = S \.}\ langle n \ pm 1 | n \ rangle = S \.

Энергия E i - это энергия ионизации, соответствующая выбранной атомной орбитали, а U - энергетический сдвиг орбитали. в результате потенциала соседних атомов. ⟨n ± 1 | H | n⟩ = - Δ {\ displaystyle \ langle n \ pm 1 | H | n \ rangle = - \ Delta}\ langle n \ pm 1 | H | n \ rangle = - \ Delta элементы, которые являются элементами межатомной матрицы Слейтера и Костера, являются энергии связи E i, j {\ displaystyle E_ {i, j}}E_{i,j}. В этой одномерной модели s-диапазона у нас есть только σ {\ displaystyle \ sigma}\ sigma -связи между s-орбиталями с энергией связи E s, s = V ss σ {\ стиль отображения E_ {s, s} = V_ {ss \ sigma}}E_ {s, s} = V_ {ss \ sigma} . Перекрытие состояний на соседних атомах равно S. Мы можем получить энергию состояния | k⟩ {\ displaystyle | k \ rangle}| k \ rangle используя приведенное выше уравнение:

H | k⟩ = 1 N ∑ n e i n k a H | n⟩ {\ displaystyle H | k \ rangle = {\ frac {1} {\ sqrt {N}}} \ sum _ {n} e ^ {inka} H | n \ rangle}H | k \ rangle = \ frac {1} {\ sqrt {N}} \ sum_n e ^ {inka} H | n \ rangle
⟨k | H | k⟩ знак равно 1 N ∑ n, m e i (n - m) k a m ​​| H | п⟩ {\ displaystyle \ langle k | H | k \ rangle = {\ frac {1} {N}} \ sum _ {n, \ m} e ^ {i (nm) ka} \ langle m | H | n \ rangle}\ langle k | H | k \ rangle = \ frac {1} {N} \ sum_ {n, \ m} e ^ {i (nm) ka} \ langle m | H | n \ rangle = 1 N ∑ n ⟨n | H | n⟩ + 1 N ∑ n ⟨n - 1 | H | п⟩ е + я к а + 1 н ∑ п ⟨п + 1 | H | п⟩ е - ika {\ displaystyle = {\ frac {1} {N}} \ sum _ {n} \ langle n | H | n \ rangle + {\ frac {1} {N}} \ sum _ {n } \ langle n-1 | H | n \ rangle e ^ {+ ika} + {\ frac {1} {N}} \ sum _ {n} \ langle n + 1 | H | n \ rangle e ^ {- ika}}= \ frac {1} {N} \ sum_n \ langle n | H | n \ rangle + \ frac {1} {N} \ sum_n \ langle n-1 | H | n \ rangle e ^ {+ ika} + \ frac {1} {N} \ sum_n \ langle n + 1 | H | n \ rangle e ^ {- ika} = E 0–2 Δ cos ⁡ (ka), {\ displaystyle = E_ {0} -2 \ Delta \, \ cos (ka) \,}= E_0 -2 \ Delta \, \ cos (ka) \,

где, например,

1 N ∑ n ⟨n | H | n⟩ знак равно E 0 1 N ∑ N 1 знак равно E 0, {\ displaystyle {\ frac {1} {N}} \ sum _ {n} \ langle n | H | n \ rangle = E_ {0} {\ frac {1} {N}} \ sum _ {n} 1 = E_ {0} \,}\ frac {1} {N} \ sum_n \ langle n | H | n \ rangle = E_0 \ frac {1} {N} \ sum_n 1 = E_0 \,

и

1 N ∑ n ⟨n - 1 | H | N⟩ е + я К а знак равно - Δ е я К а 1 N ∑ N 1 знак равно - Δ е я К а. {\ displaystyle {\ frac {1} {N}} \ sum _ {n} \ langle n-1 | H | n \ rangle e ^ {+ ika} = - \ Delta e ^ {ika} {\ frac {1 } {N}} \ sum _ {n} 1 = - \ Delta e ^ {ika} \.}\ frac {1} {N} \ sum_n \ langle n-1 | H | n \ rangle e ^ {+ ika} = - \ Delta e ^ {ika} \ frac {1} {N} \ sum_n 1 = - \ Delta e ^ {ika} \.
1 N ∑ n ⟨n - 1 | N⟩ E + I K A = S e i K a 1 N ∑ N 1 = S E i K a. {\ displaystyle {\ frac {1} {N}} \ sum _ {n} \ langle n-1 | n \ rangle e ^ {+ ika} = Se ^ {ika} {\ frac {1} {N}} \ sum _ {n} 1 = Se ^ {ika} \.}\ frac {1} {N} \ sum_n \ langle n-1 | n \ rangle e ^ {+ ika} = S e ^ {ika} \ frac {1} {N} \ sum_n 1 = S e ^ {ika} \.

Таким образом, энергия этого состояния | k⟩ {\ displaystyle | k \ rangle}| k \ rangle можно представить в знакомой форме энергетической дисперсии:

E (k) = E 0 - 2 Δ cos ⁡ (ka) 1 + 2 S соз ⁡ (ка) {\ Displaystyle E (к) = {\ гидроразрыва {E_ {0} -2 \ Delta \, \ cos (ka)} {1 + 2S \, \ cos (ka)}}}E (k) = \ гидроразрыв {E_0-2 \ Delta \, \ cos (ka)} {1 + 2 S \, \ cos (ka)} .
  • Для k = 0 {\ displaystyle k = 0}k = 0 энергия равна E = (E 0–2 Δ) / (1 + 2 S) {\ displaystyle E = (E_ { 0} -2 \ Delta) / (1 + 2S)}E = (E_0 - 2 \ Delta) / (1 + 2 S) и состояние состоит из суммы всех атомных орбиталей. Это состояние можно рассматривать как цепочку связывающих орбиталей.
  • . Для k = π / (2 a) {\ displaystyle k = \ pi / (2a)}k = \ pi / (2 a) энергия равна E = E 0 {\ displaystyle E = E_ {0}}E = E_0 и состояние состоит из суммы атомных орбиталей, которые являются множителем ei π / 2 {\ displaystyle e ^ { i \ pi / 2}}e ^ {i \ pi / 2} не в фазе. Это состояние можно рассматривать как цепочку несвязывающих орбиталей.
  • Наконец для k = π / a {\ displaystyle k = \ pi / a}k = \ pi / a энергия E = (E 0 + 2 Δ) / (1-2 S) {\ displaystyle E = (E_ {0} +2 \ Delta) / (1-2S)}E = (E_0 + 2 \ Delta) / (1-2 S) и состояние состоит из переменная сумма атомных орбиталей. Это состояние можно рассматривать как цепочку антисвязывающих орбиталей.

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

Таблица межатомных матричных элементов

В 1954 году Дж. Слейтер и Г.Ф. Костер опубликовал, в основном для расчета d-полос переходного металла, таблицу межатомных матричных элементов

E i, j (r → n, n ') = ⟨n, i | H | n ′, j⟩ {\ displaystyle E_ {i, j} ({\ vec {\ mathbf {r}}} _ {n, n '}) = \ langle n, i | H | n', j \ rangle}{\displaystyle E_{i,j}({\vec {\mathbf {r} }}_{n,n'})=\langle n,i|H|n',j\rangle }

, который также может быть получен из кубических гармонических орбиталей напрямую. В таблице представлены матричные элементы как функции от LCAO двухцентровых интегралов связей между двумя кубическими гармоническими орбиталями, i и j, на соседних атомах. Интегралами связи являются, например, V ss σ {\ displaystyle V_ {ss \ sigma}}V_{ss\sigma}, V pp π {\ displaystyle V_ {pp \ pi}}V_ {pp \ pi} и V dd δ {\ displaystyle V_ {dd \ delta}}V_{dd\delta}для связей сигма, pi и дельта (обратите внимание, что эти интегралы также должны зависеть от расстояния между атомами, т.е. являются функцией от (l, m, n) {\ displaystyle (l, m, n)}{\ displaystyle (l, m, n)} , хотя это не указывается явно каждый раз.).

Межатомный вектор выражается как

r → n, n '= (rx, ry, rz) = d (l, m, n) {\ displaystyle {\ vec {\ mathbf {r}) }} _ {n, n '} = (r_ {x}, r_ {y}, r_ {z}) = d (l, m, n)}{\displaystyle {\vec {\mathbf {r} }}_{n,n'}=(r_{x},r_{y},r_{z})=d(l,m,n)}

где d - расстояние между атомами и l, m и n - направляющие косинусы к соседнему атому.

E s, s = V ss σ {\ displaystyle E_ {s, s} = V_ {ss \ sigma}}E_ {s, s} = V_ {ss \ sigma}
E s, x = l V sp σ {\ displaystyle E_ {s, x} = lV_ {sp \ sigma}}E_ {s, x} = l V_ {sp \ sigma}
E x, x = l 2 V pp σ + (1 - l 2) V pp π {\ displaystyle E_ {x, x} = l ^ {2} V_ {pp \ sigma } + (1-l ^ {2}) V_ {pp \ pi}}E_ {x, x} = l ^ 2 V_ {pp \ sigma} + (1 - l ^ 2) V_ {pp \ pi}
E x, y = lm V pp σ - lm V pp π {\ displaystyle E_ {x, y} = lmV_ {pp \ sigma } -lmV_ {pp \ pi}}E_ {x, y} = lm V_ {pp \ sigma} - lm V_ {pp \ pi}
E x, z = ln V pp σ - ln V pp π {\ displaystyle E_ {x, z} = lnV_ {pp \ sigma} -lnV_ {pp \ pi}}E_ {x, z} = ln V_ {pp \ sigma } - ln V_ {pp \ pi}
E s, xy = 3 лм V sd σ {\ displaystyle E_ {s, xy} = {\ sqrt {3}} lmV_ {sd \ sigma}}E_ {s, xy} = \ sqrt {3} lm V_ {sd \ sigma}
E s, x 2 - y 2 = 3 2 (l 2 - m 2) V sd σ {\ displaystyle E_ {s, x ^ {2} -y ^ {2}} = {\ frac {\ sqrt {3}} {2}} (l ^ {2 } -m ^ {2}) V_ {sd \ sigma}}E_ {s, x ^ 2-y ^ 2 } = \ frac {\ sqrt {3}} {2} (l ^ 2 - m ^ 2) V_ {sd \ sigma}
E s, 3 z 2 - r 2 = [n 2 - (l 2 + m 2) / 2] V sd σ {\ displaystyle E_ { s, 3z ^ {2} -r ^ {2}} = [n ^ {2} - (l ^ {2} + m ^ {2}) / 2] V_ {sd \ sigma}}E_ {s, 3z ^ 2-r ^ 2} = [n ^ 2 - (l ^ 2 + m ^ 2) / 2] V_ {sd \ sigma}
E x, xy = 3 l 2 m V pd σ + m ( 1 − 2 l 2) V pd π {\displaystyle E_{x,xy}={\sqrt {3}}l^{2}mV_{pd\sigma } +m(1-2l^{2})V_{pd\pi }}E_ {x, xy} = \ sqrt {3} l ^ 2 m V_ {pd \ sigma} + m (1-2 l ^ 2) V_ {pd \ pi}
E x, yz = 3 lmn V pd σ − 2 lmn V pd π {\displaystyle E_{x,yz}={\sqrt {3}}lmnV _{pd\sigma }-2lmnV_{pd\pi }}E_ {x, yz} = \ sqrt {3 } lmn V_ {pd \ sigma} - 2 lmn V_ {pd \ pi}
E x, z x = 3 l 2 n V p d σ + n ( 1 − 2 l 2) V p d π {\displaystyle E_{x,zx}={\sqrt {3}}l^{2}nV_{pd\sigma }+n(1-2l^{2})V_{pd\pi }}E_ {x, zx} = \ sqrt {3} l ^ 2 n V_ {pd \ sigma} + n (1-2 l ^ 2) V_ {pd \ pi}
E x, x 2 − y 2 = 3 2 l ( l 2 − m 2) V p d σ + l ( 1 − l 2 + m 2) V p d π {\displaystyle E_{x,x^{2}-y^{2}}={\frac {\sqrt {3}}{2}}l(l^{2}-m^{2})V_{pd\sigma }+l(1-l^{2}+m^{2})V_{pd\pi }}E_ {x, x ^ 2-y ^ 2} = \ frac {\ sqrt {3}} {2} l (l ^ 2 - m ^ 2) V_ {pd \ сигма} + l (1 - l ^ 2 + m ^ 2) V_ {pd \ pi}
E y, x 2 − y 2 = 3 2 m ( l 2 − m 2) V p d σ − m ( 1 + l 2 − m 2) V p d π {\displaystyle E_{y,x^{2}-y^{2}}={\frac {\sqrt {3}}{2}}m(l^{2}-m^{2})V_{pd\sigma }-m(1+l^{2}-m^{2})V_{pd\pi }}E_ {y, x ^ 2-y ^ 2} = \ frac {\ sqrt {3}} {2} m (l ^ 2 - m ^ 2) V_ {pd \ sigma} - m (1 + l ^ 2 - m ^ 2) V_ {pd \ pi}
E z, x 2 − y 2 = 3 2 n ( l 2 − m 2) V p d σ − n ( l 2 − m 2) V p d π {\displaystyle E_{z,x^{2}-y^{2}}={\frac {\sqrt {3}}{2}}n(l^{2}-m^{2})V_{pd\sigma }-n(l^{2}-m^{2})V_{pd\pi }}E_ {z, x ^ 2-y ^ 2} = \ frac {\ sqrt {3}} {2} n (l ^ 2 - m ^ 2) V_ {pd \ sigma} - n (l ^ 2 - m ^ 2) V_ {pd \ pi}
E x, 3 z 2 − r 2 = l [ n 2 − ( l 2 + m 2) / 2 ] V p d σ − 3 l n 2 V p d π {\displaystyle E_{x,3z^{2}-r^{2}}=l[n^{2}-(l^{2}+m^{2})/2]V_{pd\sigma }-{\sqrt {3}}ln^{2}V_{pd\pi }}E_ {x, 3z ^ 2-r ^ 2} = l [n ^ 2 - (l ^ 2 + m ^ 2) / 2] V_ {pd \ sigma} - \ sqrt {3} ln ^ 2 V_ {pd \ pi}
E y, 3 z 2 − r 2 = m [ n 2 − ( l 2 + m 2) / 2 ] V p d σ − 3 m n 2 V p d π { \displaystyle E_{y,3z^{2}-r^{2}}=m[n^{2}-(l^{2}+m^{2})/2]V_{pd\sigma }-{\sqrt {3}}mn^{2}V_{pd\pi }}E_ {y, 3z ^ 2-r ^ 2} = m [n ^ 2 - (l ^ 2 + m ^ 2) / 2 ] V_ {pd \ sigma} - \ sqrt {3} mn ^ 2 V_ {pd \ pi}
E z, 3 z 2 − r 2 = n [ n 2 − ( l 2 + m 2) / 2 ] V p d σ + 3 n ( l 2 + m 2) V p d π {\displaystyle E_{z,3z^{2}-r^{2}}=n[n^{2}-(l^{2}+m^{2})/2]V_{pd\sigma }+{\sqrt {3}}n(l^{2}+m^{2})V_{pd\pi }}E_ {z, 3z ^ 2- r ^ 2} = n [n ^ 2 - (l ^ 2 + m ^ 2) / 2] V_ {pd \ sigma} + \ sqrt {3} n (l ^ 2 + m ^ 2) V_ {pd \ pi }
E x y, x y = 3 l 2 m 2 V d d σ + ( l 2 + m 2 − 4 l 2 m 2) V d d π + ( n 2 + l 2 m 2) V d d δ {\displaystyle E_{xy,xy}=3l^{2}m^{2}V_{dd\sigma }+(l^{2}+m^{2}-4l^{2}m^{2})V_{dd\pi }+(n^{2}+l^{2}m^{2})V_{dd\delta }}E_ {xy, xy} = 3 l ^ 2 m ^ 2 V_ { dd \ sigma} + (l ^ 2 + m ^ 2 - 4 l ^ 2 m ^ 2) V_ {dd \ pi} + (n ^ 2 + l ^ 2 m ^ 2) V_ {dd \ delta}
E x y, y z = 3 l m 2 n V d d σ + l n ( 1 − 4 m 2) V d d π + l n ( m 2 − 1) V d d δ {\displaystyle E_{xy,yz}=3lm^{2}nV_{dd\sigma }+ln(1-4m^{2})V_{dd\pi }+ln(m^{2}-1)V_{dd\delta }}E_ {xy, yz} = 3 lm ^ 2 nV_ {dd \ sigma} + ln (1 - 4 m ^ 2) V_ {dd \ pi} + ln (m ^ 2 - 1) V_ {dd \ delta}
E x y, z x = 3 l 2 m n V d d σ + m n ( 1 − 4 l 2) V d d π + m n ( l 2 − 1) V d d δ {\displaystyle E_{xy,zx}=3l^{2}mnV_{dd\sigma }+mn(1-4l^{2})V_{dd\pi }+mn(l^{2}-1)V_{dd\delta }}E_ {xy, zx} = 3 l ^ 2 mn V_ {dd \ sigma} + mn (1 - 4 l ^ 2) V_ {dd \ pi} + mn (l ^ 2 - 1) V_ {dd \ delta}
E x y, x 2 − y 2 = 3 2 l m ( l 2 − m 2) V d d σ + 2 l m ( m 2 − l 2) V d d π + [ l m ( l 2 − m 2) / 2 ] V d d δ {\displaystyle E_{xy,x^{2}-y^{2}}={\frac {3}{2}}lm(l^{2}-m^{2})V_{dd\sigma }+2lm(m^{2}-l^{2})V_{dd\pi }+[lm(l^{2}-m^{2})/2]V_{dd\delta }}E_ {xy, x ^ 2-y ^ 2} = \ frac {3} {2} lm (l ^ 2 - m ^ 2) V_ {dd \ sigma} + 2 l m (m ^ 2 - l ^ 2) V_ {dd \ pi} + [l m (l ^ 2 - m ^ 2) / 2] V_ {dd \ delta}
E y z, x 2 − y 2 = 3 2 m n ( l 2 − m 2) V d d σ − m n [ 1 + 2 ( l 2 − m 2) ] V d d π + m n [ 1 + ( l 2 − m 2) / 2 ] V d d δ {\displaystyle E_{yz,x^{2}-y^{2}}={\frac {3}{2}}mn(l^{2}-m^{2})V_{dd\sigma }-mn[1+2(l^{2}-m^{2})]V_{dd\pi }+mn[1+(l^{2}-m^{2})/2]V_{dd\delta }}E_ {yz, x ^ 2-y ^ 2 } = \ frac {3} {2} mn (l ^ 2 - m ^ 2) V_ {dd \ sigma} - mn [1 + 2 (l ^ 2 - m ^ 2)] V_ {dd \ pi} + mn [1 + (l ^ 2 - m ^ 2) / 2] V_ {dd \ delta}
E z x, x 2 − y 2 = 3 2 n l ( l 2 − m 2) V d d σ + n l [ 1 − 2 ( l 2 − m 2) ] V d d π − n l [ 1 − ( l 2 − m 2) / 2 ] V d d δ {\displaystyle E_{zx,x^{2}-y^{2}}={\frac {3}{2}}nl(l^{2}-m^{2})V_{dd\sigma }+nl[1-2(l^{2}-m^{2})]V_{dd\pi }-nl[1-(l^{2}-m^{2})/2]V_{dd\delta }}E _ {zx, x ^ 2-y ^ 2} = \ frac {3} {2} nl (l ^ 2 - m ^ 2) V_ {dd \ sigma} + nl [1 - 2 (l ^ 2 - m ^ 2)] V_ {dd \ pi} - nl [1 - (l ^ 2 - m ^ 2) / 2] V_ {dd \ delta}
E x y, 3 z 2 − r 2 = 3 [ l m ( n 2 − ( l 2 + m 2) / 2) ] V d d σ − 2 l m n 2 V d d π + [ l m ( 1 + n 2) / 2 ] V d d δ ] {\displaystyle E_{xy,3z^{2}-r^{2}}={\sqrt {3}}\left[lm(n^{2}-(l^{2}+m^{2})/2)]V_{dd\sigma }-2lmn^{2}V_{dd\pi }+[lm(1+n^{2})/2]V_{dd\delta }\right]}E_ {xy, 3z ^ 2-r ^ 2} = \ sqrt { 3} \ left [lm (n ^ 2 - (l ^ 2 + m ^ 2) / 2)] V_ {dd \ sigma} - 2 lmn ^ 2 V_ {dd \ pi} + [lm (1 + n ^ 2) / 2] V_ {dd \ delta} \ right]
E y z, 3 z 2 − r 2 = 3 [ m n ( n 2 − ( l 2 + m 2) / 2) V d d σ + m n ( l 2 + m 2 − n 2) V d d π − [ m n ( l 2 + m 2) / 2 ] V d d δ ] {\displaystyle E_{yz,3z^{2}-r^{2}}={\sqrt {3}}\left[mn(n^{2}-(l^{2}+m^{2})/2)V_{dd\sigma }+mn(l^{2}+m^{2}-n^{2})V_{dd\pi }-[mn(l^{2}+m^{2})/2]V_{dd\delta }\right]}E_ {yz, 3z ^ 2-r ^ 2} = \ sqrt {3} \ left [mn (n ^ 2 - (l ^ 2 + m ^ 2) / 2) V_ {dd \ sigma} + mn (l ^ 2 + m ^ 2 - n ^ 2) V_ {dd \ pi} - [mn (l ^ 2 + m ^ 2) / 2] V_ {dd \ delta } \ right]
E z x, 3 z 2 − r 2 = 3 [ l n ( n 2 − ( l 2 + m 2) / 2) V d d σ + l n ( l 2 + m 2 − n 2) V d d π − [ l n ( l 2 + m 2) / 2 ] V d d δ ] {\displaystyle E_{zx,3z^{2}-r^{2}}={\sqrt {3}}\left[ln(n^{2}-(l^{2}+m^{2})/2)V_{dd\sigma }+ln(l^{2}+m^{2}-n^{2})V_{dd\pi }-[ln(l^{2}+m^{2})/2]V_{dd\delta }\right]}E_ {zx, 3z ^ 2-r ^ 2} = \ sqrt {3} \ left [ln (n ^ 2 - (l ^ 2 + m ^ 2) / 2) V_ {dd \ sigma} + ln (l ^ 2 + m ^ 2 - n ^ 2) V_ {dd \ pi} - [ln (l ^ 2 + m ^ 2) / 2] V_ {dd \ delta} \ right]
E x 2 − y 2, x 2 − y 2 = 3 4 ( l 2 − m 2) 2 V d d σ + [ l 2 + m 2 − ( l 2 − m 2) 2 ] V d d π + [ n 2 + ( l 2 − m 2) 2 / 4 ] V d d δ {\displaystyle E_{x^{2}-y^{2},x^{2}-y^{2}}={\frac {3}{4}}(l^{2}-m^{2})^{2}V_{dd\sigma }+[l^{2}+m^{2}-(l^{2}-m^{2})^{2}]V_{dd\pi }+[n^{2}+(l^{2}-m^{2})^{2}/4]V_{dd\delta }}E_ {x ^ 2-y ^ 2, x ^ 2-y ^ 2} = \ frac {3} {4} (l ^ 2 - m ^ 2) ^ 2 V_ {dd \ sigma} + [l ^ 2 + m ^ 2 - (l ^ 2 - m ^ 2) ^ 2] V_ {dd \ pi} + [n ^ 2 + (l ^ 2 - m ^ 2) ^ 2/4] V_ {dd \ delta}
E x 2 − y 2, 3 z 2 − r 2 = 3 [ ( l 2 − m 2) [ n 2 − ( l 2 + m 2) / 2 ] V d d σ / 2 + n 2 ( m 2 − l 2) V d d π + [ ( 1 + n 2) ( l 2 − m 2) / 4 ] V d d δ ] {\displaystyle E_{x^{2}-y^{2},3z^{2}-r^{2}}={\sqrt {3}}\left[(l^{2}-m^{2})[n^{2}-(l^{2}+m^{2})/2]V_{dd\sigma }/2+n^{2}(m^{2}-l^{2})V_{dd\pi }+[(1+n^{2})(l^{2}-m^{2})/4]V_{dd\delta }\right]}E_ {x ^ 2-y ^ 2,3z ^ 2-r ^ 2} = \ sqrt {3} \ left [(l ^ 2 - m ^ 2) [n ^ 2 - (l ^ 2 + m ^ 2) / 2] V_ {dd \ sigma} / 2 + n ^ 2 (m ^ 2 - l ^ 2) V_ {dd \ pi} + [(1 + n ^ 2) (l ^ 2 - m ^ 2) / 4] V_ {dd \ delta } \ right]
E 3 z 2 − r 2, 3 z 2 − r 2 = [ n 2 − ( l 2 + m 2) / 2] 2 V dd σ + 3 n 2 (l 2 + m 2) V dd π + 3 4 (l 2 + m 2) 2 V dd δ {\ displaystyle E_ {3z ^ {2 } -r ^ {2}, 3z ^ {2} -r ^ {2}} = [n ^ {2} - (l ^ {2} + m ^ {2}) / 2] ^ {2} V_ { dd \ sigma} + 3n ^ {2} (l ^ {2} + m ^ {2}) V_ {dd \ pi} + {\ frac {3} {4}} (l ^ {2} + m ^ { 2}) ^ {2} V_ {dd \ delta}}E_ {3z ^ 2-r ^ 2,3z ^ 2-r ^ 2} = [n ^ 2 - (l ^ 2 + m ^ 2) / 2] ^ 2 V_ {dd \ sigma} + 3 n ^ 2 (l ^ 2 + m ^ 2) V_ {dd \ pi} + \ гидроразрыв {3} {4} (l ^ 2 + m ^ 2) ^ 2 V_ {dd \ delta}

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

См. Также
Ссылки
На Викискладе есть материалы, относящиеся к Дисперсионные отношения электронов.
  • N. У. Эшкрофт и Н. Д. Мермин, Физика твердого тела (Thomson Learning, Торонто, 1976).
  • Стивен Бланделл Магнетизм в конденсированных средах (Оксфорд, 2001).
  • С.Маекава и др. Физика оксидов переходных металлов (Springer-Verlag Berlin Heidelberg, 2004).
  • Теория зон Джона Синглтона и электронные свойства твердых тел (Оксфорд, 2001).
Дополнительная литература
  • Уолтер Эшли Харрисон (1989). Электронная структура и свойства твердых тел. Dover Publications. ISBN 0-486-66021-4.
  • N. У. Эшкрофт и Н. Д. Мермин (1976). Физика твердого тела. Торонто: Thomson Learning.
  • Дэвис, Джон Х. (1998). Физика низкоразмерных полупроводников: Введение. Кембридж, Соединенное Королевство: Издательство Кембриджского университета. ISBN 0-521-48491-X.
  • Goringe, C M; Bowler, D.R; Эрнандес, Э (1997). «Плотно связывающее моделирование материалов». Отчеты о достижениях физики. 60 (12): 1447–1512. Bibcode : 1997RPPh... 60.1447G. doi : 10,1088 / 0034-4885 / 60/12/001.
  • Slater, J.C.; Костер, Г. Ф. (1954). «Упрощенный метод ЛКАО для периодической потенциальной проблемы». Физический обзор. 94 (6): 1498–1524. Полномочный код : 1954PhRv... 94.1498S. doi : 10.1103 / PhysRev.94.1498.
Внешние ссылки
Последняя правка сделана 2021-06-11 11:59:14
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте