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

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

Φ(x) = ciΨi (x)

(III.16)

i

 

Здесь под “ x ” подразумевается вся совокупность координат, характеризующих систему. Поскольку Φ описывает некоторое состояние системы, то средняя величина энергии в этом состоянии:

Φ

 

 

 

 

 

Φ ( x )

=

 

( x ) | H |

 

 

 

 

 

 

 

 

c jψ j ( x )

 

c iψ i ( x ) | H |

=

i

 

 

 

 

 

 

i

 

∑ ∑ c i c *j

 

 

 

ψ i ( x ) | H | ψ j ( x )

=

i

j

 

 

 

 

(III.17)

∑ ∑ c i c *j ψ i ( x ) | ψ j ( x ) E j

=

i

j

 

 

 

 

 

∑ ∑ c i c *j δ i j E j =

 

 

i

j

 

 

 

 

 

=

 

c i

 

2 E i

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

Заменим в (III.17) все значения E наименьшим собственным значением E0 . В результате это выражение или

уменьшится или останется неизменным. Таким образом:

 

E0 = E0 (| ci |2 = E0

 

Φ | H | Φ ≥ | ci |2

(III.18)

i

i

 

Очевидно знак равенства в (III.18) имеет место, когда Φ совпадает с ψ 0 , то есть является решением уравнения Шре-

дингера, отвечающего минимуму полной энергии E0 . Так что теорема доказана.

86

Вариационный метод

Этот метод очень важен для химических задач, и основан на вариационном принципе.

Допустим теперь, что мы не знаем ни собственных значений, ни собственных функций уравнения (III.15), а только

ˆ

гамильтониан H . Возьмем некоторую произвольную нормированную функцию Φ и запишем равенство:

 

Согласно условию (III.18) E E0 . Возьмем

E = Φ | H | Φ

другую нормированную функцию Φ ' и опять запишем равенство:

 

 

, где E ' E0 .

 

(III.19)

E ' =

Φ ' | H | Φ '

 

Если

окажется,

что

E ' < E ,

то это

значит, что

E0 E ' E , то есть величина E ' более близка к искомому

значению E0 . Следовательно, функция Φ ' выбрана более

удачно чем Φ .

 

 

 

 

Иными словами мы

должны

подобрать пробную

функцию Φ так, чтобы интеграл I

 

был равен

Φ | H | Φ

его наименьшему возможному значению для данной функции. Для этого надо, чтобы вариация интеграла δI по функции Φ обращалась в ноль, то есть, должно быть

δ Φ Φ = .

| H | 0

Варьируя пробную функцию Φ , надо сохранять условие её нормировки, то есть, чтобы выполнялось требование

Φ | Φ = 1. Это можно учесть в соответствии с методом Ла-

гранжа, введя неопределенный множитель Лагранжа ε , и записав условие варьирования:

87

 

 

 

= 0

(III.20)

δ

Φ | H | Φ − ε ( Φ | Φ −1)

 

 

 

 

 

 

Здесь стоит напомнить, что E = Φ | H | Φ зависит от

вида функции Φ , поэтому E является функцией от функции, т.е. функционалом. Поскольку мы хотим найти наименьшее значение энергии E E0 , эта задача состоит в

поиске экстремума E . А исследование экстремальных значений функционалов проводится методами вариационного исчисления. Понятие «вариация» является обобщением понятия «дифференциал». Обозначается символом δ и операции с вариациями можно проводить также как с обычными дифференциалами. Практически для того, чтобы удовлетворить условию (III.20), из каких-либо физических соображений подбирают пробную функцию Φ как функцию некоторых параметров так, что её варьирование осуществляется с помощью варьирования этих параметров. Чем удачнее выбран вид Φ и чем больше параметров, тем глубже получится минимум интеграла. Итак, уравнение (III.20) служит для нахождения приближения Φ → Φ0 к наименьшему состоя-

нию

ψ 0 . Интеграл Φ 0 | H | Φ 0 = ε 0 определяет энергию

ε 0 ,

которая при должном выборе пробной функции Φ

имеет наиболее близкое к E0 значение. В этом заключается суть этого метода.

Одноэлектронное приближение, уравнения Хартри

Запишем гамильтониан системы электронов и ядер – то есть атома, молекулы ( i, j – номера электронов, α, β – номера ядер ):

88

ˆ

 

 

2

 

 

 

 

2

 

2

 

 

 

α

2

 

 

 

 

 

 

H = −

 

 

 

i

 

 

 

 

 

 

 

 

 

2m

2

 

Mα

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

(III.21)

1

∑∑

Z

 

e2

 

1

∑∑

Zα Zβ

 

 

1

∑∑

e2

 

+

 

+

при i j

 

 

α

 

 

 

 

 

 

 

 

 

2

 

riα

2

 

Rαβ

 

 

2

rij

i α

 

 

 

α β

 

 

 

 

i j

 

Здесь первый член – это оператор кинетической энергии электронов, второй член – оператор кинетической энергии ядер, третий член – оператор потенциальной энергии взаимодействия ядер и электронов между собой, четвертый член – оператор потенциальной энергии взаимодействия ядер между собой, пятый член – оператор потенциальной энергии взаимодействия электронов между собой. В силу приближения Борна-Оппенгеймера, оператор кинетической энергии ядер равен нулю. Оператор потенциальной энергии взаимодействия ядер между собой можно положить также равным нулю, поскольку его вклад в полную электронную волновую функцию при фиксированных положениях ядер в пространстве постоянен и не зависит от состояния системы.

Тогда гамильтониан может быть записан следующим образом:

2

i

2

1

∑∑

Zα e

2

 

1

∑∑e

2

 

H = −

 

 

+

 

,

(III.22)

2m

2

riα

 

2

 

 

 

 

i α

 

 

i j rij

 

 

при i≠j

Введем единицы Хартри:

= 2hπ = 1 m = 1

e = 1

1 а.е. = 0,529177 Å

89

Тогда:

ˆ

1

 

2

 

1

 

 

Zα

 

 

1

 

1

 

 

 

H = −

2 i

 

∑∑i α

 

 

+

 

∑∑i j

 

где

i j

(III.23)

2

r

 

2

r

 

 

 

 

 

 

 

 

 

 

iα

 

 

 

 

ij

 

 

 

и

 

1

 

 

 

 

Zα

 

 

1

 

 

 

 

 

 

 

ˆ

 

2

+

 

 

 

 

 

 

 

Hi = −

 

i

 

 

 

 

 

 

 

 

где

i j

(III.24)

2

 

r

r

 

 

 

 

 

 

 

α

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iα

 

 

 

ij

 

 

 

 

 

 

Запишем формулы, выражающие основные положения одноэлектронного приближения:

H = H(qi )

(III.25)

i

Φ(q1,q2 ,q3,..,qi ,...) =φ1(q1) φ2 (q2 ) φ3 (q3 ) φi (qi )

 

причем qi – набор пространственных координат частиц –

электронов. Можно сказать, что одноэлектронное прибли-

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

Причем каждый одноэлектронный гамильтониан действует только на одноэлектронную функцию того же самого электрона:

ˆ

) ϕi

(qi ) = Ei

ϕi (qi )

H (qi

Посмотрим, как связана Ei

i

(III.26)

в (III.26) с полной энер-

гией системы E . Запишем стационарное уравнение Шредингера в одноэлектронном приближении с учетом того, что

90