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