- •Российский химико-технологический университет им. Д.И. Менделеева Кафедра квантовой химии
- •Методы вычислительной химия наноразмерных систем
- •Приближение независимых частиц. Одноэлектронная модель
- •Гамильтониан атома с зарядом ядра Z|е| (начало координат на ядре):
- •Гамильтониан атома с зарядом ядра Z|е| (начало координат на ядре):
- •Полный гамильтониан атома в приближении независимых частиц
- •Метод самосогласованного поля
- •Отталкивание между электронами i и j учтено дважды: как среднее по j и
- •Как это сделать?
- •Однако, чтобы вычислить этот интеграл,
- •Приближение центрального поля
- •Атомные орбитали и их характеристики
- •Радиальные нормированные функции водородоподобных атомов Rnl
- •Радиальные составляющие 1s (а), 2s (б), 3s (в) орбиталей атома водорода
- •Радиальные составляющие
- •Основные свойства радиальных функций.
- •Угловые функции Ylm ( , ) собственные функции оператора квадрата углового момента L2
- •Угловые части волновой функции атома, обладающего центральным полем
- •Угловые части волновой функции атома, обладающего центральным полем
- •Обозначения орбиталей с различными угловыми зависимостями
Российский химико-технологический университет им. Д.И. Менделеева Кафедра квантовой химии
Компьютерное моделирование процессов нанотехнологий.
Лекция 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 |
4π 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 |
4π 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 j}χ1i |
(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