Особенности различных вариантов метода Хартри-Фока.
Неограниченный и ограниченный методы ХартриФока.
Вспомним, прежде всего, сами уравнения Хартри - Фока. Это базовые уравнения современной квантовой химии. Последовательным усложнением или упрощением их получаются все остальные вычислительные методы. Эти уравнения и все, основанные на них, подходы называются также методом молекулярных орбиталей.
|
|
N |
| (2) |2 |
N |
* (2) |
|
(2) |
|
H(1) |
k (1) |
|
i |
d 2 k (1) |
i |
|
k |
|
r |
|
r |
|
|
||||
|
|
i 1 |
12 |
i 1 |
|
12 |
|
|
k 1,2,.........N
N |
|
|
d 2 i (1) ki |
i |
(1) |
i 1 |
|
(1) |
|
|
N |
| (2) |2 |
N |
* (2) |
|
(2) |
|
|
|
|
|
H (1) |
k (1) |
|
i |
d 2 k (1) |
i |
|
k |
|
d 2 |
i (1) |
k |
k (1) |
r |
|
r |
|
|
||||||||
|
|
i 1 |
12 |
i 1 |
|
12 |
|
|
|
|
|
|
(2)
при k 1,2,.........N .
При рассмотрении систем с замкнутыми оболочками, после интегрирования по спиновым переменным получаются уравнения:
|
n |
|
| |
(q ) |2 |
|
* (q ) |
|
(q ) |
|
|
|
H (1) k (q1 ) 2 |
i |
2 |
dv2 k (q1 ) |
i 2 |
k |
2 |
dv2 |
i |
(q1 ) |
||
|
r12 |
r12 |
|
|
|||||||
|
i 1 |
|
|
|
|
|
|
|
|
||
k k (q1 ) |
|
|
|
|
|
|
|
|
|
|
(3) |
Далее вводятся кулоновский и обменный операторы
6
ˆ |
k (1) |
| i |
(2) |2 |
d 2 k (1) |
|
|
||||
Ji |
|
|
|
|
(4а) |
|||||
|
|
|
|
|||||||
|
|
|
|
r12 |
|
|
|
|
|
|
ˆ |
k (1) |
i* (2) k |
(2) |
d 2 |
i |
(1) |
(4б) |
|||
Ki |
|
r12 |
|
|||||||
|
|
|
|
|
|
|
|
|
||
ˆ
Оператор Ji имеет довольно простую интерпретацию: это
кулоновский потенциал, создаваемый в области нахождения первого электрона зарядом второго электрона, находящегося на другой орбитали. Причем, плотность распределения этого электрона задается квадратом модуля спин – орбитали
i (2) i (2) 2 . По этой причине этот оператор называют
кулоновским.
ˆ
Оператор Ki по сравнению с кулоновским оператором имеет более сложный характер - при действии на
функцию k (1) |
он переводит еѐ в функцию i (1) . Этот |
оператор не имеет классического аналога, в отличие от кулоновского оператора, хотя также определяет потенциал кулоновского типа, создаваемый в области нахождения первого электрона зарядом второго электрона. Но здесь второй электрон находится не на одной орбитали, а распределен по двум орбиталям. Или, как говорят
ˆ
делокализован по двум орбиталям. Оператор Ki выполняет
перестановку электронов, то есть их обмен. Поэтому этот оператор называют обменным.
И сами уравнения Хартри – Фока (2) для систем не обязательно имеющих замкнутые оболочки будут такими:
эф |
k |
k k , |
|
|
|
(5) |
|
H |
|
|
|
||||
|
ˆ ýô |
|
ˆ |
N |
ˆ |
ˆ |
|
где |
|
|
) |
||||
H |
H (1) |
(Ji |
Ki |
||||
|
|
|
|
i 1 |
|
|
|
Оператор |
ˆ ýô |
часто записывают так: |
|||||
H |
|||||||
7
ˆ ˆ |
ˆ |
ˆ |
(6) |
F H (1) |
J |
K |
И его называют оператором Фока или фокианом. Уравнения
(5) при этом записываются так:
|
|
|
F k |
k k |
(7) |
Ограниченный метод Хартри –Фока для замкнутых оболочек.
Уравнения Хартри - Фока в форме (2) записаны для спин-орбиталей в виде произведения пространственной и спиновой частей ψi = φiα или ψi = φiβ. Поскольку заполнение спин - орбиталей мы задаем сами, то и число электронов с разными проекциями спина известно заранее. То есть в результате решения уравнений Хартри - Фока надо найти пространственные части – молекулярные орбитали φi.
В системах с замкнутыми оболочками число электронов со спином α равно числу электронов со спином β. Считается, что каждая спин – орбиталь занята двумя электронами, один из которых имеет спин α, а второй – спин β. При этом пространственная орбиталь для обоих электронов одна и та же. Иначе говоря, на каждой пространственной орбитали находится два электрона. В этом
изаключается суть ограниченного метода Хартри –Фока для систем с замкнутыми оболочками. Сокращенно - метод ОХФ или restricted Hartree - Fock, что в англоязычной литературе
иквантово-химических программах обозначается как RHF. Эти рассуждения справедливы только для атомных и молекулярных систем с замкнутыми оболочками. Отметим, что, хотя для каждой орбитали имеется свой собственный оператор Фока, вид этих операторов для разных электронов одинаков, то есть можно говорить об одном операторе Фока. В результате процесса самосогласования электроны располагаются по два на орбиталях, соответствующих n
8
низшим энергиям εi, где n = N/2, а N – общее число электронов в системе. При этом возникает вопрос – почему данный метод называется ограниченным методом Хартри – Фока? Прежде всего, потому, что в процессе получения самосогласованных решений обеспечивается требование сохранения правильной симметрии волновой функции. Она должна быть собственной функцией оператора квадрата полного спина системы. Именно в этом и заключается его ограниченность. И одноэлектронные функции также должны отвечать определенным требованиям, то есть принадлежать неприводимым представлениям пространственной группы симметрии данной системы. Однако, при этом возникает одна сложность принципиального характера. Ее смысл заключается в том, что при получении наилучшей энергии системы теряются свойства симметрии ее волновой функции, а при получении волновой функции правильной симметрии - повышается энергия данной системы. Это соображение известно, как дилемма симметрии Левдина в методе Хартри
– Фока. Левдин это крупнейший теоретик химии.
Неограниченный метод Хартри –Фока.
Вернемся к уравнениям (2) и после интегрирования по спину запишем их отдельно для электронов со спином α и спином β:
|
|
|
|
|
|
|
|
|
|
F |
|
|
|
|
|||||
i |
(1) |
i |
i |
(1) |
|
(8) |
|||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|||
F i |
(1) |
i i |
(1) |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Операторы Фока F |
|
и F |
определяются соотношениями: |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
H (1) J |
K |
|
||||||
F |
(9) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
H (1) J |
K |
|
||||||
F |
|
||||||||
9
Эти операторы Фока различаются обменными частями, поскольку в обменном взаимодействии участвуют только электроны с одинаковыми спинами. В результате для α и β спинов получаются разные пространственные функции. Эти уравнения не являются независимыми - они связаны между собой через кулоновские взаимодействия между электронами с различными спинами. Их решениями будут два набора орбиталей - φiα и φiβ. В ходе процедуры самосогласования электроны со спином α размещают по одному на орбиталях с низшими энергиями εiα , а электроны со спином β – также по одному на орбиталях с низшими энергиями εiβ. Такой метод называется неограниченным методом Хартри – Фока. Здесь не накладываются условия сохранения правильной симметрии волновой функции, и он называется неограниченным (НХФ или Unrestricted Hartee –Fock – UHF).
То есть волновая функция системы не является собственной функцией оператора квадрата полного спина системы, что является ее недостатком. Суть этого метода может быть выражена так – разные орбитали для электронов с разными спинами.
Схема размещения электронов по уровням энергии молекулярных орбиталей выглядит так:
Рис.1. Схема размещения электронов по уровням энергии МО в неограниченном методе Хартри - Фока.
10