Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лекции / Lections_nano_3.ppt
Скачиваний:
76
Добавлен:
13.02.2015
Размер:
623.1 Кб
Скачать

Российский химико-технологический университет им. Д.И. Менделеева Кафедра квантовой химии

Компьютерное моделирование процессов нанотехнологий.

Лекция 3. Неэмпирические методы расчета строения и свойств молекул и кластеров. Свойства электронной волновой функции. Приближение Борна-Оппенгеймера. Методы Хартри-Фока и Кона-Шэма.

Цирельсон В.Г., Бобров М.Ф. «Многоэлектронный атом». Цирельсон В.Г. , Бобров М.Ф. «Квантовая химия молекул».

Москва 2007 г.

Методы вычислительной химия наноразмерных систем

молекулярная

механика

полуэмпирическая квантовая химия

квантовая

молекулярная динамика

статистическая

и метод Монте-Карло

механика

 

2

Приближение независимых частиц. Одноэлектронная модель

Оператор потенциальной энергии электронов:

N

N

e

2

 

 

VЭЭ (r)

 

 

,

 

4 r

i

j

 

0

ij

 

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

Как можно найти волновые функции и уровни энергии неподвижного N- электронного атома, максимально близкие к точным?

3

Гамильтониан атома с зарядом ядра Z|е| (начало координат на ядре):

 

2 N

N

Ze2

N N

e2

H Tэ (r) Vэя (r, R) Vээ (r)

 

i2

 

 

 

 

2m

4π r

4π r

 

 

i

i

0 i

i j

0 ij

Исключим Vээ

hiχi (ri ) εiχi (ri )

i 1,2,3,...,N

 

hi

2

2

Ze

2

одноэлектронный

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

2m

0ri

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

Электрон i описывается волновой функцией χi (ri ) и имеет энергию εi

Иначе: поведение каждого электрона не зависит от поведения остальных электронов и описывается некоторой волновой функцией подобно единственному электрону в атоме водорода.

χi (rВi )этомназываютсясостоит сутьодноэлектроннымиприближения независимыхволновыми функциямичастиц.или орбиталями.

4

Гамильтониан атома с зарядом ядра Z|е| (начало координат на ядре):

 

2 N

N

Ze2

N N

e2

H Tэ (r) Vэя (r, R) Vээ (r)

 

i2

 

 

 

 

2m

4π r

4π r

 

 

i

i

0 i

i j

0 ij

Исключим Vээ

hiχi (ri ) εiχi (ri )

i 1,2,3,...,N

 

hi

2

2

Ze

2

одноэлектронный

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

2m

0ri

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

Электрон i описывается волновой функцией χi (ri ) и имеет энергию εi

Иначе: поведение каждого электрона не зависит от поведения остальных электронов и описывается некоторой волновой функцией подобно единственному электрону в атоме водорода.

χi (rВi )этомназываютсясостоит сутьодноэлектроннымиприближения независимыхволновыми функциямичастиц.или орбиталями.

5

Полный гамильтониан атома в приближении независимых частиц

 

N

 

H hi

 

1

 

 

 

 

Собственная функция Н

χ1(r1) χ1(r2) χ3(r3) χ N (rN )

 

 

 

 

 

(волновая функция Хартри)

 

 

 

 

 

Энергия атома

Ψ

 

H

 

Ψ ε1 ε2 ε3 ... εN E

 

 

6

Метод самосогласованного поля

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

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

 

 

 

 

 

 

 

 

 

 

 

 

Одноэлектронный гамильтониан:

N

 

 

2

 

 

 

ССП

 

 

2

2

 

Ze

2

e

 

 

 

 

 

 

 

 

 

 

 

 

 

hi

 

 

i

 

 

 

 

 

 

 

 

 

, j 1,2,3,...,(N -1)

2m

4

 

r

4

0

r

 

 

 

 

 

 

 

0 i

j

 

 

 

ij

j

Оператор описывает отталкивание между электронами i и j, усредненное по

всем положениям электрона j.

 

 

n

 

 

 

 

 

 

 

 

 

 

Собственные функции

 

 

Н hССПi

-

гамильтониана:

 

 

 

i

 

 

 

Орбитальные произведения:

χ1 (r1 ) χ 2 (r2 ) χ 3 (r3 ) χ N (rN )

 

n

 

n

| h

ÑÑÏ

 

 

Собственные значения

E' ε

 

 

|

 

i

i

 

Н:

 

 

 

 

 

i

 

i

 

 

 

 

i есть сумма кинетической энергии i-го электрона, потенциальной энергии его

притяжения к ядру и средняя потенциальная энергия его отталкивания от

 

остальных электронов.

7

 

Отталкивание между электронами i и j учтено дважды: как среднее по j и среднее по i hвÑÑÏj .

С учетом этого, полная энергия атома равна:

 

 

 

1

n

n

 

 

 

 

e2

 

 

 

 

E E'

 

 

 

 

 

 

 

 

 

 

i j

2

4 r

 

 

 

i

j

 

 

 

j

 

 

 

 

 

 

 

 

 

 

0

ij

i,

Гамильтониан атома должен иметь

 

 

 

 

 

 

 

 

 

 

вид:

 

 

 

 

 

 

e2

 

 

 

 

 

 

 

 

1

 

n

n

 

 

 

 

 

 

 

H H -

 

 

 

 

 

 

 

 

 

 

 

,

i j

2

4

 

r

 

 

 

 

 

i

j

 

0

 

 

i,

j

 

 

 

 

 

 

 

ij

 

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

 

2

 

 

.

следует прежде рассчитать усредненныеe

величины

 

 

 

 

 

 

 

4

0

r

 

 

 

 

ij j

 

8

Как это сделать?

Вероятность того, что электрон j с

волновой функцией j(rj) находится в

χ 2j dv j

.

бесконечно малом объеме dvj, равна

Электрон i в точке, последовательно “пробегающей” все положения в пространстве

Отталкивание электрона i, усредненное по всем положениям электрона j, равно:

 

e

2

 

 

 

 

χ

j

2

(r

)

 

 

 

 

 

e

2

 

 

 

 

j

 

 

dv j

 

4

r

 

 

 

4

r

 

 

 

 

 

0 ij j

 

 

 

 

 

 

 

0 ij

 

9

Однако, чтобы вычислить этот интеграл,

волновые функции χ

j

(r

) должны уже быть известны!

 

j

 

Задаются некоторым набором

 

N одноэлектронных функций, максимально

близких к правильным

χ 0j (rj )

и строят оператор 0 ССП

(hi )

 

2

 

 

 

χ

j

2 (r

)

 

. С их помощью вычисляют eинтеграл

2

 

 

 

j

 

 

 

 

 

 

 

e

 

 

 

 

 

dv j

 

4

r

 

 

 

4

r

 

 

 

 

0 ij j

 

 

 

 

 

 

0 ij

 

. Затем решают набор одноэлектронных уравнений

Хартри, возникающий из условия минимума среднего значения гамильтониана,

вычисляемого с волновой функцией Хартри:

 

 

 

2

Ze2

 

N

 

0

 

2

 

 

 

 

 

 

 

 

 

j

(rj)]

 

 

 

 

 

 

{

 

i2

 

 

 

e2

 

 

 

 

dv j1i

(ri) εiχ1i (ri),

i 1,2,3,..., N

 

 

2m

4

r

4

r

 

 

 

 

 

 

o i

 

 

 

 

 

o ij

 

 

 

 

 

 

 

 

 

 

 

 

 

j i

 

 

 

 

 

 

 

Полученные решения1

(ri )

используют, чтобы построить "исправленный"

оператор

 

χ i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

ССП

 

1

)

ССП

 

 

 

 

 

 

 

 

 

 

 

 

(hi

)

 

(hi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, вновь решают систему уравнений, но теперь – с и т.д. Этот процесс называется самосогласованием, а результирующее поле, создающее усредненный потенциал в (40), называется самосогласованным10

Соседние файлы в папке Лекции