Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
246
Добавлен:
22.08.2013
Размер:
6.91 Mб
Скачать

3.9.2. Метод функционала плотности

Необходимо решить

Hψ = Eψ

уравнение Шредингера:

Ψ - волновая функция, зависит от

 

 

координат и спина всех электронов

Плотность электронов

Поиск многочастичной функции –

проводимости в металле ~1023см-3

дело безнадежное.

 

К счастью

Знание такой ψ совершенно необязательно

Можно упростить задачу

Перейти к описанию системы

на языке плотности заряда.

 

Гамильтониан

T −оператор кинетической энергии

H = T + U + V U −оператор электрон− электронного взаимодействия

V −оператор взаимодействия электронов с внешним полем

в т.ч. с положительным фоном

Считаем основное состояние не вырождено

Распределение электронов в пространстве однозначно определяется внешним потенциалом v(r)

Справедливо и

v(r) - однозначный функционал плотности

 

 

обратное

 

 

v[n(r)]

 

с точностью до аддитивной постоянной

 

 

 

 

утверждение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ(r) также функционал n(r)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ[n(r)]

 

)

ψ = T[n(r)]

 

Средние, вычисленные на функциях ψ,

 

 

 

 

T =

ψ

T

 

 

 

 

 

 

 

 

 

 

E = E[n(r)]

 

 

однозначные функционалы плотности

 

 

 

 

 

 

Выделим из Е части, связанные с дальнодействующим взаимодействием

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r r

1

 

n(r)n(r1)

 

1

 

r

 

 

 

 

 

 

 

E[n(r )] = òv(r )n(r )dτ +

2

òò

 

r r1

 

 

dτdτ

 

+ G[n(r )]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r − r

 

 

 

 

 

 

 

 

 

 

 

 

 

G[n(r)] - некоторый функционал

Хоэнберг и Коэн доказали следующую теорему:

Существует универсальный функционал G[n(r)] такой, что плотность, c

оответствующая любому внешнему потенциалу, должна минимизировать функционал при дополнительном условии:

Физический смысл G[n(r)]

G[n(r )] = T[n(r )] + E

 

[n(r )] + E

[n(r )]

 

r

r

x

r

c

r

T[n(r)] - функционал кинетической энергии

Ex[n(r)] функционал обменной энергии.

Электрон окружен пространством, в котором имеется недостаток (если забыть о нашем электроне) частиц с параллельным нашему электрону спином.

Принцип Паули

Обменная дырка

EС[n(r)] функционал корреляционной

энергии, Электроны - заряженные частицы

ЕС - разность между точной величиной энергии и ее Хартри-Фоковским значением

Расталкиваются

Движение электронов не является независимым

кулоновскими силами

В среднем находятся дальше друг от друга, чем это

 

было бы при хаотическом движении.

Корреляционная дырка Обменно-корреляционная дырка

Трудность при использовании МФП

Удалось в двух предельных случаях

Электронный газ с почти постоянной плотностью

~

~

n(r ) = n + n(r )

 

òn(r)dτ = 0

n~ << n

Для электронного газа с медленно изменяющейся плотностью

Ñn(r) << kF n(r)

kF - волновое число на поверхности Ферми

n =

объем в простр.импульсов

4

3

πp3

× 2 = 2 ×

 

F

 

объем, занимаемый одним сост.

 

 

h3

Конкретный вид функционала G[n(r)].

=8π pF3

3 h3