книги из ГПНТБ / Татевский В.М. Классическая теория строения молекул и квантовая механика
.pdf§ 4. Общие замечания о молекулярных орбиталях
Как было изложено выше, на варьируемые функции <р в рас смотренных методах накладываются только следующие предвари тельные условия — однозначность, непрерывность и нормированность (в методе Фока — обычно ортонормированность). Опти мальные функции ф должны удовлетворять уравнениям Хюккеля. (XXVIII,9), Хартри (XXVIII,23) или Фока (XXVIII, 35) в зависи мости от выбранного метода.
Оптимальные функции ф, вообще говоря, могут быть (в рамках указанных условий) весьма разнообразными по своей форме и распределению их значений, отличных от нуля, в отдельных обла стях пространства вокруг ядер. Некоторые из этих функций могут
оказаться приближенно центрированными на |
отдельных |
ядрах |
(т. е. обладающими значениями, отличными |
от нуля |
только |
в сравнительно небольшой области пространства, включающей одно определенное ядро). Другие функции ф могут оказаться суще ственно отличными от нуля в области пространства, включающей два ядра, три ядра и, наконец, все ядра молекулы.
Никаких ограничений такого рода на функции ф заранее не накладывается. Конкретная форма их определяется только урав нениями {XXVIII, 9) или (XXVIII, 35) в зависимости от выбран ного варианта метода МО и заранее (без конкретного исследования данной задачи) форму функции ф определить детальнее нельзя.
Следует - отметить далее, что если приближенная многоэлек тронная функция Ф конструируется в виде простого произведения
одноэлектронных функций ф (как в вариантах Хюккеля и |
Хар |
три), то она заведомо не обладает свойством антисимметрии, |
тре |
буемым квантовой механикой, и не может правильно описать со стояний систем в отношении вероятности распределения электронов в пространстве вокруг ядер. Поэтому использовать входящие в та кую функцию Ф одноэлектронные функции ф для описания «состоя ний» отдельных электронов неправомерно. Результаты такого опи сания будут противоречить общим теоремам квантовой механики.
Если |
же |
приближенная |
функция |
конструируется |
в |
виде |
||
определителя |
Фока и обладает требуемым квантовой |
механикой |
||||||
свойством |
антисимметрии по |
отношению к перестановке простран- • |
||||||
ственных |
и- спиновых координат любой |
пары |
электронов, |
то |
и |
|||
в этом случае |
конкретный вид одноэлектронных |
функций ф, из |
ко |
|||||
торых строится приближенная многоэлектронная функция Ф для системы, не представляет существенного интереса и не может дать никаких дополнительных данных о состоянии и свойствах системы или ее отдельных электронов по сравнению с самой многоэлектрон ной функцией Ф для системы.
|
Действительно,- если определена |
некоторая |
система |
функций |
Ф,, |
ф д , / 2 как оптимальная для |
построения |
функции |
Фока Ф |
(XXVIII, 2.7), описывающей основное состояние системы, то любая система N/2 линейно независимых ортонормированных функций
Фі, . . . » |
ф^/2 > |
являющихся |
линейными комбинациями |
ИСХОДНЫХ |
функции |
ф р |
ф д , / 2 , |
дает определитель Фока |
Ф' вида |
(XXVIII,27), описывающий то же состояние системы, что и ис
ходный определитель |
Ф. |
|
|
|
|
При таком переходе от исходного |
набора |
одноэлектронных |
|||
функций |
ф р . . . » q>N/2 |
и соответствующей этому |
набору |
прибли |
|
женной функции Ф к другому набору |
одноэлектронных |
функций |
|||
ФР |
ф ^ и соответствующей этому |
набору функции Ф' описа |
|||
ние состояния системы не меняется. Не |
меняются |
и значения всех |
|||
физических величин, |
которые могут быть вычислены методами кван |
||||
товой механики для системы при вычислении их с помощью функ ции Ф или Ф' и сопоставлены с соответствующими эксперименталь
ными значениями. Следовательно, описания состояния |
системы |
||
наборами одноэлектронных функций |
ф , |
Ф ^ и ф р |
ф ^ / 2 |
совершенно эквивалентны, в то время как свойства |
отдельных |
||
функций этих разных наборов могут |
быть существенно |
различны. |
|
Например, функции ф первого набора могут быть более или менее центрированы на отдельных ядрах или парах ядер, а функции ф '
второго набора «размазаны» по |
всему пространству |
вокруг ядер |
и т. п. |
|
функций ф |
Таким образом, конкретный |
вид одноэлектронных |
не является существенным и не может дать никаких более деталь ных сведений о состоянии системы или отдельных ее частей, чем те, что даются многоэлектронной волновой функцией Ф системы.
§ .6. Вариант ЛКАО метода молекулярных орбиталей (прямой вариационный метод Ритца решения уравнений Хюккеля или Фока)
Мы изложили выше три варианта метода молекулярных орби талей (варианты Хюккеля, Хартри и Фока). Собственно метод молекулярных орбиталей для основного состояния и состоит в из
ложенных выше |
методах |
нахождения |
приближенной функции Ф |
и приближенного |
значения |
энергии Е. |
Однако для практического |
осуществления этого метода нужно уметь решать уравнения Хюк келя (XXVIII, 9), Хартри (XXVIII, 23) или Фока (XXVIII, 35). Эти уравнения много проще исходного электронного уравнения, так как искомые функции ф в этих уравнениях зависят только от трех де картовых координат х, у, г, в то время как функция W в электрон ном уравнении (или приближенная функция Ф) зависит от 3N де картовых координат. Однако решение и этих уравнений, т.е. опре деление из них функций ф и собственных значений є, представляет
собой даже для |
простейших молекул столь сложную задачу, что |
|||||
для ее решения |
снова нужно находить специальные приближенные |
|||||
приемы. Одним |
из |
таких |
специальных |
чисто математических |
при |
|
емов является так |
называемый прямой |
вариационный |
метод |
Рит |
||
ца. Этот метод |
применим |
к решению |
математических |
уравнений |
||
разного характера, а не только уравнений типа Хюккеля, Хартри или Фока. Сущность этого метода (в его линейном варианте) состоит в том, что какая-либо неизвестная функция представляется в виде линейной комбинации каких-либо известных более или менее про извольно выбранных функций % тех же переменных, что и искомая функция. Коэффициенты этой линейной комбинации варьируются затем так, чтобы линейная комбинация заранее выбранных, извест ных, более или менее произвольных функций х была в определен ной области значений аргументов или во всей области значений аргументов наиболее близка к искомой функции.
Применение прямого вариационного метода Ритца к решению уравнений Хюккеля, Хартри или Фока и представляет собой сущ ность так называемого варианта ЛКАО («линейная комбинация атомных орбиталей») метода молекулярных орбиталей. Вариант ЛКАО мы поясним на примерах его приложения к решению урав нений Хюккеля. Приложение его к решению уравнений Хартри и уравнений Фока, в принципе, аналогично, но сложнее. Мы не будем рассматривать применение этого метода для уравнений Хартри, применение его к решению уравнений Фока иллюстрируем только конспективно, поясняя главным образом конечные результаты.
§ 6. Приложение прямого вариационного метода Ритца (вариант ЛКАО) к решению уравнения Хюккеля
Уравнение Хюккеля (XXVIII, 9) имеет вид:
Н0 (х, у, г) ф (х, у, г) = е<р (х, у, z) |
( X X V I I I , 41) |
Будем искать приближенное решение этого уравнения в форме
т |
|
Ф (*> У, г) = ^ Ср%р (*- У> 2) |
(XXVIII, 42) |
где %р (Х, у, г) — заранее заданные функции, непрерывные, одно значные и нормированные, т. е. удовлетворяющие условию
{ Хр (х, у, z) %р (х, у, z)dx=l |
(XXVIII, 43) |
ав остальном произвольные.
Вчастности, функции %р могут быть центрированы в любых точках пространства вокруг ядер *. Число т функций %р должно быть не меньше числа искомых функций ф, чтобы последние были линейно независимы, но может быть больше. Функции ф предпола-
* Напомним, что рассматривается задача приближенного решения электрон ного уравнения и соответствующего уравнения Хюккеля (XXVIII, 9) для задан ной ядерной конфигурации, т. е. при фиксированных координатах ядер. Это важно отметить потому, что при разных ядерных конфигурациях можно выби рать разные наборы функций %р , и никаких асимптотических условий, например условий, соответствующих диссоциации молекулы на свободные атомы, в задаче, сформулированной таким образом, не накладывается,
В этом уравнении все вариации 6СР и ЬСЧ можно рассматривать как независимые, поэтому коэффициент при каждой вариации дол жен обращаться в нуль и, следовательно:
2 Шм - eSpq) Сч = 0 |
( X X V I I I , 53а) |
я |
|
И |
|
|
|
|
|
р=\, |
2, ..., т |
|
|
|
|
|
|||
|
|
|
|
2 ( ^ p ? - e s p ? ) c ; = o |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
р |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q— |
1, |
2 |
|
т |
|
|
|
( X X V I I I , 536) |
|
Можно показать, |
что достаточно |
рассматривать одну |
из этих си |
||||||||||||
стем, гак как они эквивалентны. |
|
|
|
|
|
|
|
|
|||||||
|
Система |
(XXVIII, 53а) в развернутом виде будет |
|
|
|||||||||||
|
{Нп |
— е) С] + |
(Я 12 — eS1 2 ) С 2 |
+ |
. . . +{Hlm-zSm)Cm |
|
= 0 |
||||||||
|
|
|
(Нт |
- |
e S m l ) С, + . . . + |
(Нтт |
- |
є) Ст = |
О |
( X X V I I I , 54) |
|||||
|
Эта система однородных уравнений имеет решения |
Си |
. . . , С т , |
||||||||||||
отличные от нуля, только если определитель ее равен |
нулю: |
||||||||||||||
|
|
|
#п — |
е # ] 2 |
— eSij ... Hlm |
— |
eSlm |
0 |
( X X V I I I , 55) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
= |
|||
|
|
|
Нml |
в$ті |
|
|
|
Нтт |
е |
|
|
|
|
||
Это уравнение в общем случае имеет т решений |
|
|
|
||||||||||||
|
|
|
|
|
|
ё ь |
ё 2 , ...,lm |
|
|
|
|
(XXVIII, 56) |
|||
которые и будут приближенными |
собственными значениями є урав |
||||||||||||||
нения Хюккеля |
(XXVIII, 41). |
Подставив |
одно из них, например |
||||||||||||
Ik, |
в уравнение |
(XXVIII, 54), |
определим |
набор |
коэффициентов |
||||||||||
Сій, |
C2h, |
• • • , Cmh, |
ОТВЄЧаЮЩИХ |
ЭТОМУ решению eft. |
|
|
Хюккеля |
||||||||
|
Тогда |
одно из приближенных |
решений уравнения |
||||||||||||
(XXVIII, 41) будет |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
Ф ^ |
= |
2 |
С |
Л |
|
|
|
( X X V I I I , 57) |
|
|
|
|
|
|
|
|
р |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k = 1, 2 |
|
т |
|
|
|
|
|
||
|
Выбирая из найденных приближенных |
функций |
ФЙЛ72 функции, |
||||||||||||
отвечающие |
наименьшим собственным значениям |
(пусть |
это будет |
||||||||||||
її, |
8 2 , |
|
ёлг/г), |
получим |
систему |
N/2 |
|
искомых приближенных |
|||||||
решений |
уравнения |
Хюккеля ф,, |
|
vpNj2 |
и N/2 |
соответствующих |
|||||||||
им приближенных собственных значений. Если искомые прибли
женные функции ер,, |
cpN/2 найдены, |
т . е . определены коэффи |
|
циенты CPk в выражениях этих функций |
(XXVIII,57) |
через функ |
|
ции %р, Т О , подставив выражения ф й вида |
(XXVIII, 57) |
в выражение |
|
Тогда |
выражение |
для энергии |
(XXVIII, 30) |
может |
быть при |
||||||||||||||
ведено к виду |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
о, Э |
1 |
|
р, q |
|
V |
ft |
/ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+2 2 G P « r S |
ffi |
С;кС,Л ф С ; , С , Д |
( X X V I I I , 70) |
||||||||||||
где |
|
|
|
|
р. Я г, s |
|
|
\ |
k |
|
] |
\ l |
|
|
j |
|
|
|
|
|
|
|
|
Gpqrs — 2Jpqrs |
|
— Kpqrs |
|
|
|
|
( X X V I I I , 71) |
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
Если |
выразить |
HVQ |
в виде |
(XXVIII, 64), |
то выражение |
для Е° |
|||||||||||||
будет |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
о, р |
р |
|
р, ? |
|
\ |
ft |
|
I |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
*Ь 2 2 ^ P ? Y (SC p f t C 9 f e ) + |
S 2 GPI?''S (S c pftc ?s J ( 2 C r ' C * ' |
j |
|||||||||||||||||
|
P,q |
У |
\ k |
|
J |
|
p, q r, s |
|
\ |
ft |
/ |
\ |
J |
|
|
/ |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( X X V I I I , 72) |
||
§ 8 . Замечания о выборе функций % в прямом |
|
|
|
|
|
||||||||||||||
вариационном методе Ритца (вариант ЛКАО метода МО) |
|
||||||||||||||||||
В прямом |
вариационном |
методе |
Ритца по его существу |
|
не мо |
||||||||||||||
жет быть никаких однозначных рецептов для выбора набора |
зара |
||||||||||||||||||
нее заданных |
функций |
%, через |
которые представляется |
прибли |
|||||||||||||||
женно |
некоторая |
искомая |
функция. |
При |
решении |
уравнений |
|||||||||||||
Хюккеля или уравнений Фока выбор функций % ограничен |
|
только |
|||||||||||||||||
тем, что они должны |
|
иметь |
интегрируемый |
квадрат |
модуля, т. е. |
||||||||||||||
могут |
быть выбраны |
|
нормированными |
и должны |
удовлетворять |
||||||||||||||
общим |
требованиям |
однозначности |
и |
непрерывности, |
поскольку |
||||||||||||||
этим требованиям |
должны удовлетворять функции ср. |
|
|
|
|
||||||||||||||
В остальном конкретный вид функций |
%, выбираемых для |
||||||||||||||||||
апроксимаций |
функций |
ср в |
методе |
Ритца, в общем |
|
случае |
одно |
||||||||||||
значно определен быть не может. В зависимости от особенностей
функций |
ср, которые можно |
предполагать, |
исходя |
из уравнений, |
||
которым |
они должны удовлетворять |
(в рассматривавшихся слу |
||||
чаях это уравнение Хюккеля |
или уравнение Фока), каждый, ре |
|||||
шающий подобную задачу, может пытаться |
выбрать набор функций |
|||||
X таким |
образом, чтобы |
их линейная |
комбинация |
могла передать |
||
особенности неизвестных |
функций ф более |
или менее удовлетвори |
||||
тельно. Ясно, что разные наборы одинакового числа функций % будут приводить к разным приближениям и что дополнение неко торого набора из щ функций % еще одной или несколькими функ циями может улучшить приближение, но не может его ухудшить.
Если никакие специальные аналитические свойства функций ф неизвестны заранее и не могут быть сформулированы из общих физических соображений, то выбор функций х может основываться
только |
на математической |
интуиции. |
Если некоторые (хотя бы |
|||
качественные) соображения |
о свойствах функции |
ф |
могут |
быть |
||
сделаны на основании исследования уравнений, |
которым |
они |
||||
должны |
удовлетворять, или |
из общих |
физических |
соображений, |
||
то они могут быть учтены при выборе функций х- |
|
|
|
|||
Мы остановимся только |
на одном вопросе — на |
вопросе о том, |
||||
в каких областях пространства вокруг ядер функции % должны иметь наибольшие по модулю значения или вообще достаточно большие значения, чтобы по возможности лучше были апроксимированы функции ф методом Ритца. Иными словами, как целе сообразно центрировать функции %.
Решить этот вопрос строго и однозначно математически на осно вании исследования свойств уравнений, определяющих функции ф, очень трудно. Поэтому, вообще говоря, для конкретной молекулы можно пробовать использовать набор функций х> центрированных в любых точках пространства вблизи от ядер (порядка атомных размеров), например вдоль линий, соединяющих пары ядер атомов молекулы, связанных химическими связями (с точки зрения клас сической теории строения), если предполагается определенная фор мула химического строения молекулы.
Однако из общих физических соображений можно предполагать, что полная функция W, описывающая какое-либо электронное со
стояние молекулы, должна быть такова, что |
квадрат ее модуля |
||
1 F*4f должен иметь большие значения |
для |
таких |
конфигураций |
электронов, когда они распределены в |
непосредственной близости |
||
от ядер (не обязательно точно в месте |
нахождения |
ядер). Иными |
|
словами, при значениях координат электронов, близких к коорди
натам ядер, можно предположить, что |
должно |
иметь относи |
тельно большие значения. |
|
|
Квадрат модуля Ч*, т.е. ЧГ *ЧГ , дает |
вероятность |
определенного |
распределения электронов в пространстве, а поскольку электроны
находятся |
в |
поле ядер, |
вероятность их |
нахождения вблизи ядер |
|
должна быть больше, |
чем в областях пространства, более далеких |
||||
от ядер. Предполагая |
это для функции |
можно предполагать это |
|||
также и |
для |
функций |
ф, из которых |
Ч? конструируется. Отсюда |
|
представляется целесообразным из указанных соображений выби рать набор функций х т а к > чтобы каждая функция набора была центрирована на каком-либо ядре. Тогда в набор функций % будут
входить группы функций таких, что каждая функция |
определенной |
||||
группы, центрирована на одном определенном |
ядре *. |
|
|
|
|
* Естественно, можно вводить в |
рассмотрение |
функции |
% |
более |
общего |
вида, например центрированные сразу |
на нескольких |
ядрах или |
на всех |
ядрах |
|
(т. е. имеющие отличные от нуля относительно большие значения только в об ластях пространства, примыкающих к соответствующим ядрам). Указанный выше способ центрирования функций % является поэтому частным способом, соответ ствующим изложенному рассуждению, которое его в какой-то мере объясняет качественно и обосновывает возможность указанного выбора вида функций ф из общих физических соображений.
§ 9. Распределение отрицательного |
электрического заряда |
|
в пространстве вокруг ядер в варианте ЛАКО |
|
|
метода молекулярных орбиталей |
|
|
Выше было выведено выражение |
для |
pe(x,y,z)—плотности |
отрицательного электрического заряда в пространстве вокруг ядер, создаваемого электронами. Это выражение имеет одинаковую ма тематическую форму как в методе Хюккеля *, так и в методе Фока. Именно для четного N
|
|
2 |
|
|
|
|
|
|
|
Ре{х, у, |
z) = |
- 2 2 « Р І ( * . У . z)q>k(x,y,z) |
|
(XXVIII,73) |
|||||
Отличия в функции ре, полученной в этих двух методах, будут |
|||||||||
состоять в том, что функции фй(х, у, z), |
полученные |
из уравнения |
|||||||
Хюккеля и из уравнения Фока, будут разными. |
|
|
|
||||||
В варианте ЛКАО |
функции |
щ(х,у„г) |
представляются в виде |
||||||
линейных комбинаций |
некоторых** |
заданных |
функций |
%P(x,y,z) |
|||||
|
|
|
т |
|
|
|
|
|
|
|
|
Ч>* - |
2 |
с |
Р к Ь |
|
|
|
|
|
|
|
P=I |
|
|
|
|
|
|
Подставляя ф/, в этой |
форме |
в выражение для ре |
(XXVIII, 73), |
||||||
получим |
|
|
|
|
|
|
|
|
|
Ре (*> У. *) = |
2 |
Ь (*• У> z) |
X, (*• У> 2 ) ( 2 |
СІкСдЛ |
( X X V I I I , 74) |
||||
|
P. q |
|
|
|
|
\ f t = l |
|
/ |
|
Это выражение для ре(х,у,г) мы используем в одной из следую щих глав при рассмотрении вопроса о распределении отрицатель ного электрического заряда в пространстве вокруг ядер для разных фрагментов, относящихся к одному и тому же типу и виду, встре чающемуся в разных молекулах или в одной молекуле.
Выражение для ріг в приближении Фока—Рузана. Из выраже ния для ріг в приближении Фока, приведенном в § 3 гл. XXVIII, и выражения для щ (XXVIII, 60) следует, что в приближении Фока— Рузана ріг будет выражаться формулой
Ріг (*і. Уи Zu х 2 , уъ z2 ) =
|
= 2 2 |
WP (U X, (1) *; (2) X, (2) - |
Хр ( I ) X, (2) х; (2) X, О)] |
PpqP„ |
||||
|
p. q г, s |
|
|
|
|
|
|
|
ГДЄ |
|
|
, |
|
|
. |
' |
|
|
|
Ppq ~ |
2 CpkCqk' |
Prs = |
2 |
^rfisl |
|
|
|
|
|
к |
|
I |
|
|
|
* Поскольку |
р в (х, у, г) |
зависит |
только |
от |
вида приближенной |
волновой |
||
функции |
Ф, а этот |
вид в методах Хюккеля и Хартри одинаков, выражение для |
||||||
р е (х, у, г) |
в методе |
Хартри |
имеет ту ж е математическую форму, что и в |
методе |
||||
Хюккеля. Отличие состоит в конкретном виде функций фА (х, у, г) в методах Хюк келя, Хартри и Фока.
** Естественно, что при использовании одного и того ж е набора функций %Р коэффициенты СРк, апроксимирующие функции Ф А , полученные как решения - уравнения Хюккеля, будут отличаться от коэффициентов С» А, апроксимирующих функции ф*, полученные как решения уравнений Фока.
