Φ(1, 2,...2n) =
(2n!)− |
1 |
det{Ψ1 (1) |
|
1 (2)Ψ 2 (3) |
|
2 (4)...Ψ n (2n − 1) |
|
n (2n)} |
(III.64) |
2 |
|
|
|
|
|||||
Ψ |
Ψ |
Ψ |
|
Черта над спин-орбиталью означает, что ей отвечает противоположный спин. Например, если ψ (i) = ϕ (qi ) α (i) , то
ψ (i) = ϕ (qi ) β (i) , где qi – набор пространственных коорди-
нат электронов. Уравнения Хартри-Фока для такой системы будут:
|
ϕk (q1 ) + |
|
|
|
|
|
||
H (1) |
dv2 |
ϕk |
(q1 ) − ∫ϕi |
(q2 ) ϕk (q2 ) |
dv2 ϕi (q1 ) = |
|||
∑ 2∫ | ϕi |
(q2 ) | |
|||||||
n |
|
|
2 |
|
|
* |
|
|
|
|
|
|
|
|
|
r12 |
|
i=1 |
|
|
r12 |
|
|
|
|
|
ε k ϕk (q1 ) |
|
|
|
|
|
(III.65) |
||
Причины суммирования до n , а не до 2n, и наличие двойки в кулоновском члене, как и ее отсутствие в обменном члене, прояснится дальше при рассмотрении выражения для полной энергии такой системы.
Как мы уже говорили εк являются орбитальными энергиями. Выясним более полно физический смысл εк. Энергия такой системы, как мы видели раньше,
|
|
2n |
|
* |
|
|
|
1 |
2n 2n |
|ψ i (1) |2 |
|ψ j |
(2) |2 |
|
||
E0 |
= ∑∫ψ i |
(1)H (1)ψ i (1)dτ1 |
+ |
|
∑∑∫ |
|
|
|
dτ1dτ 2 − |
||||||
|
|
r12 |
|
||||||||||||
|
|
i=1 |
|
|
|
|
|
|
2 i=1 j=1 |
|
|
|
|||
1 |
2n |
2n |
|
ψ i* (1) ψ *j |
(2) ψ j (1) ψ i (2) |
|
|
|
|
|
|||||
∑∑ |
∫ |
|
|
|
|
|
|
|
dτ1dτ 2 |
|
|
|
|||
|
|
|
r12 |
|
|
|
|
|
|
||||||
2 i=1 |
j=1 |
|
|
|
|
|
|
|
|
|
|
|
|
||
(III.66)
106
Проинтегрируем выражение (III.66) по спинам в предположении, что спин-орбитали являются произведением пространственной и спиновой функций. Получим следующее выражение для энергии:
E0 =
n |
∫ |
|
|
|
|
|
|
|
|
|
n |
∫ |
|
|
|
|
|
|
|
|
||||
|
ϕi*(q1)H(q1)ϕi (q1)dv1 |
|
|
|
|
ϕi*(q1)H(q1)ϕi (q1)dv1 β | β + |
|
|||||||||||||||||
∑ |
|
α |α +∑ |
|
|||||||||||||||||||||
i=1 |
|
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
α |α α |α |
+ |
|
|
|
|
|
|
|
|
|
|
|
ϕi*(q1) ϕi (q1) ϕ*j (q2) ϕj (q2) |
|
|
|
|
|
|
|
|
|
|||||||
|
1 |
|
n |
|
n |
|
|
|
|
β | β β | β + |
|
|
|
|||||||||||
|
|
|
∑∑ |
|
|
|
|
|
|
|
|
1 2 |
|
|
|
|
|
|
|
|||||
+ |
2 i 1 |
j 1 ∫ |
|
r |
|
|
|
|
|
|
dv dv |
|
|
|
α |α β | β |
+ |
− |
|
||||||
|
|
|
|
= |
|
= |
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
β | β α |α |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
n |
|
n |
ϕi*(q1) ϕ*j (q2) ϕj (q1) ϕi (q2) |
|
|
α |α β | β + |
|
|
|
||||||||||||
− |
|
|
∑∑ |
|
|
|
|
|
|
|
|
1 2 |
|
|
= |
|
||||||||
|
2 i 1 |
j 1 ∫ |
|
r |
|
|
|
|
|
|
|
|
|
β | β α |α |
|
|
|
|
||||||
|
|
|
|
= |
|
= |
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
* |
|
|
|
|
1 |
n n |
|
|
ϕi*(q1) ϕi |
(q1) ϕ*j (q2) ϕj (q2) |
|
||||||||
2∑ |
|
ϕi (q1)H(q1)ϕi (q1)dv1 |
+ |
|
|
∑∑4 |
|
|
|
|
|
|
|
dv1dv2 |
− |
|||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
i=1 |
∫ |
|
|
|
|
|
|
2 i=1 j=1 |
∫ |
|
|
|
|
r12 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
1 |
n |
|
n |
|
|
* |
* |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
∑∑2∫ |
ϕi (q1) |
ϕj (q2) ϕj (q1) ϕi (q2) |
dv1dv2 = |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||
2 i=1 |
j=1 |
|
r12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
nn n
=2∑Hii +∑∑(2Jij −Kij ).................................................................(III.67)
i=1 i=1 j=1
где
|
|
* |
|
ˆ |
|
|
|
|
|
|
|
|
|
|
|
|
Hii = ∫ϕi |
|
(q1 )H (q1 )ϕi (q1 )dv1 |
|
|
|
|
|
|
||||||||
Jij = ∫ |
ϕ* |
(q ) ϕ |
(q ) ϕ* (q |
) ϕ |
(q |
) |
|
|
|
|||||||
|
i |
|
1 |
i |
1 |
|
j |
2 |
j |
|
2 |
|
|
dv1dv2 |
(III.68) |
|
|
|
|
|
|
|
r12 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Kij = ∫ |
ϕ* (q ) ϕ* (q |
) ϕ |
(q ) ϕ |
(q |
) |
|
|
|||||||||
i |
|
1 |
|
j |
2 |
j |
|
1 |
i |
2 |
|
|
dv1dv2 |
|
||
|
|
|
|
|
|
r12 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
107 |
|
|
|
|
|
Прокомментируем наличие четырех комбинаций спиновых функций в кулоновском интеграле и только двух – в обменном. В кулоновском интеграле каждый электрон находится на одной орбитали и поэтому в этом случае возможны четыре различные комбинации: оба электрона имеют одинаковые спины (либо оба имеют спин α, либо оба – β), либо они имеют разные спины ( первый имеет спин α, а второй имеет спин β, либо наоборот). В обменном интеграле электроны распределены по обеим орбиталям и поэтому возможны только две комбинации. Например, пусть электрон с условным номером 1 на орбитали φ*i имеет спин α. Тогда на орбитали φj этот электрон тоже имеет спин α , поскольку это один и тот же электрон. А электрон с номером два на двух других орбиталях может по принципу Паули иметь только спин β. Вторая возможная комбинация соответствует спинам, противоположным первой комбинации. Никакая третья комбинация спинов здесь невозможна без нарушения принципа Паули. В чем можно легко убедиться, следуя приведенным логическим рассуждениям. Эти рассуждения проясняют и вид выражения (III.65). Итак, энергия системы с замкнутыми оболочками имеет вид:
n |
n |
n |
|
E0 = 2∑Hii |
+ ∑∑(2Jij − Kij ) |
(III.69) |
|
i=1 |
i=1 |
j=1 |
|
Удалим теперь из системы один электрон в состоянии
ψk с каким-нибудь определенным спином.
Φ( k ) (1, 2,..., 2n − 1) =
|
− |
1 |
|
|
|
|
α |
|
|
|
|
|
|
|
[(2n − 1)!] |
|
det{ϕ1 ϕ1 ...ϕ k −1 ϕ k −1 |
ϕ k +1 |
ϕ k +1 ...ϕ 2 n −1 ϕ 2 n −1} , (III.70) |
||||||||||
2 |
ϕ k |
|
||||||||||||
|
|
|
|
|
|
|
|
β |
|
|
|
|
|
|
α
где β означает α или β .
108
Энергия такой системы является разностью между энергией системы, содержащей 2n электронов и вкладом в эту энергию электрона в состоянии ψ k . Нетрудно показать,
что разность между энергиями системы с 2n электронами и 2n-1 электроном равна
n |
|
E0 − Ek = H kk + ∑(2J ki − Kki ) |
(III.71) |
i=1 |
|
Эта разность представляет собой (с обратным знаком) величину энергии ионизации электрона из состояния ψ k .
С другой стороны, если уравнение (III.65) умножить на ψ k и проинтегрировать по всему пространству, то полу-
чим выражение для орбитальной энергии:
n |
|
ε k = H kk + ∑(2J ki − Kki ) |
(III.72) |
i=1 |
|
Сравнивая (III.71) и (III.72), видим, что ε k |
= E0 − Ek , |
то есть ε k – есть энергия ионизации из системы электрона, находящегося в состоянии ψ k . В этом и заключается теорема
Купменса. Чаще ее произносят так:
Орбитальная энергия равна потенциалу ионизации электрона с этой орбитали, взятому с противоположным знаком
Заметим, что
n |
n |
n n |
|
∑ε k = ∑ H kk + ∑∑(2J ki − Kki ) |
(III.73) |
||
k =1 |
k =1 |
k =1 i=1 |
|
Следовательно, сумма орбитальных энергий εк не равна полной энергии системы, как видно из сравнения выражений (III.69) и (III.73). Это очень важное замечание, показы-
109
вающее, что простота всех приведенных рассуждений является скорее кажущейся. Дальше нам предстоит познакомиться с практическим методом решения уравнений ХартриФока, которым является линейный вариационный метод.
Линейный вариационный метод, уравнения Рутаана
В квантово-химической литературе он известен больше как метод МО ЛКАО, что означает “Молекулярные ор-
битали в виде линейных комбинаций атомных орбита-
лей”. Такой способ решения уравнений Хартри-Фока был предложен Рутааном. И по существу его правильно называть методом Хартри-Фока-Рутаана. Оба эти названия (МО ЛКАО и Хартри-Фока-Рутаана) до сих пор широко используются в литературе.
Но, в связи с использованием в современных программных комплексах обозначений метода Хартри-Фока- Рутаана или МО ЛКАО как UHF или RHF (ROHF), названия МО ЛКАО и Хартри-Фока-Рутана постепенно заменяются этими новыми названиями. Это связанно с тем, что в современной квантовой химии исследователи используют готовые программные комплексы и, естественно, мало задумываются над сутью этих методов, их физическим смыслом.
Когда мы рассматривали теорему Купменса, мы уже фактически освободились от спиновых функций.
H (1) ϕk (q1 ) + |
|
|
|
|
|
|
|
|
|
|
|
|
|||
n |
| ϕ |
(q |
) |2 |
|
ϕ* |
(q |
) |
ϕ |
|
(q |
) |
|
|
|
|
∑ 2∫ |
i |
2 |
|
dτ 2ϕk (q1 ) − ∫ |
i |
2 |
|
|
k |
2 |
|
dτ 2 |
ϕ1 |
(q1 ) |
= |
|
r12 |
|
r12 |
|
|||||||||||
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
ε k ϕk (q1 ) |
|
|
|
|
|
|
|
|
|
|
|
|
(III.74) |
||
|
|
|
|
|
110 |
|
|
|
|
|
|
|
|
|
|