Материал: Боженко Основы квантовой химии

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Φ(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

=2Hii +∑∑(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 = 2Hii

+ ∑∑(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