ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 146
Скачиваний: 0
188 Часть II. Введение в квантовохимические расчеты
б) Ортогонализация. Из тех же соображений, что и выше, мы используем длязаписи интегралов обозначения Дирака:
J |
= |
<Фа I Фб > = <“ !&>; |
< ху | jc* - у* > = ~ |
< (2) - |
(-2 ) | (2) + (-2 ) > = |
= -j - [ <2 |2 >- <- 2 |2 > + <2 |- 2 >- <- 2 |- 2 >] =
= 4 " [1 - 0 + 0 -1 ] =
=0 .
2.2.Радиальная составляющая
Теперь запишем дифференциальное уравнение, ре шением которого является радиальная функция; мы получим его, если подставим в уравнение (1) <р = R - Y f и используем решение (3):
d*R |
2 |
dR |
, |
f 8 я 2т |
f D , |
е* \ |
1(1+ 1) |
/?= о. |
dr* |
г |
dr |
[ |
h* |
\ С + |
г ) |
г* |
Это дифференциальное уравнение решается при по мощи следующих функций (которые при г -> оо должны превращаться в нуль):
Rn (р) = рЧ-n+i (р) ё~ .
1—I |
2meV |
= |
2аг , |
„ |
При этом р = |
|
—-— (а было введено |
в разд. 5.2 |
части I), a Ln+i1— так называемый присоединенный по лином Лагерра, который является решением дифферен циального уравнения
+ ( 2 / + 2 - р ) 4 ^ + |
1) L = 0 |
при условии, что п — I — 1 равно 0 или целому положи тельному числу; п = 1, 2, 3 ... (где п ^ / + 1), причем п называется главным квантовым числом.
Обобщим эти результаты на «водородоподобные» ато мы и ионы с одним валентным электроном в рамках одно
2. Одноэлектронный атом: функции водорода |
189 |
электронной модели. Введем эффективный заряд ядра Z'
и атомарную единицу длины |
|
ct0 = |
^ |
— боровский |
||||||
радиус (см. разд. 3.2); |
тогда р оказывается равным 2Z'r . |
|||||||||
|
|
|
|
|
|
|
|
|
|
Таблица 3 |
Нормированные одноэлектронные атомные |
функции |
|||||||||
|
__р_ |
|
|
|
|
|
|
|
|
|
Is = |
е |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
_ _ L |
|
|
|
_ J L |
|
|
||
2s = |
(32) |
2 (2 — р) е |
2 |
|
|
|||||
|
|
|
|
___р _ |
|
|
|
|||
2рг = |
(32) |
2 ре |
2 cos % |
|
|
|||||
|
|
|
1 |
|
р |
|
|
|
|
|
2рх = |
(32) Y ре |
2 sin ft cos <p |
|
|||||||
|
|
_ i _ |
_ |
JL |
|
|
|
|
|
|
2py = |
(32) |
2 |
pe |
2 |
sin ft sin Ф |
|
||||
|
|
_ _1_ |
|
|
|
|
|
|
||
3s = |
(972) |
2 (6 — 6p + |
p2)e |
2 |
|
|||||
|
|
|
|
|
|
|
|
_ jl |
|
|
3рг --= (648) |
2 (4p — p2)e |
2 |
cos ft |
|
||||||
|
|
_ i _ |
|
|
|
_ JL |
|
|
||
3px = (648) |
2 |
(4p — p2)e |
2 sin ft cos Ф |
|||||||
|
|
_ J_ |
|
|
|
|
|
|
||
3pu = (648) |
2 (4p — p2)e |
2 sin ft sin Ф |
||||||||
|
|
|
_ J_ |
_p_ |
|
|
|
|||
3d.%= |
(2592) |
|
2 p2e |
|
|
|
(3 cos2 ^ — 1) |
|||
|
|
|
___ L |
_ |
|
£_ |
|
|
|
|
3dxz = |
(2592) |
2 p2e |
|
2 sin 2ft cos Ф |
|
|||||
|
|
|
_ |
_L |
_ |
£ . |
|
|
|
|
3dyz = |
|
(2592) |
2 p2e |
|
2 |
sin 2ft sin Ф |
|
|||
|
|
|
_ _ i_ |
_ |
|
|
|
|
||
3dx2 _ lj2 = (2592) |
2 p2e |
|
2 |
sin2 8 cos 2ф |
|
|||||
|
|
|
___1_ |
_ |
_P |
|
|
|
3dxy = (2592) 2 p2e 2 sin2 ft sin 2ф
190Часть II. Введение в квантовохимические расчеты
Втабл. 3 приведены нормированные атомные функ ции для одноэлектронного атома [общий коэффициент
Расчет собственных значений по уравнению (2) здесь не дается, так как при квантовохимических расчетах используют прежде всего собственные функции. Собст венные значения бывают нужны только для сравнения, а для полуэмпирических оценок чаще всего применяют экспериментальные значения. В разд. 5.3 (часть I) мы
вывели формулу Еп — - ^ Е г, которая для одноэлектрон
ного атома с эффективным зарядом ядра Z' преобразовы вается в выражение
1Г*(?
па 2а0
Точный расчет многоэлектронных атомов (методХартри — Фока) выходит за рамки данного курса (здесь будут приведены только некоторые предварительные за мечания).
2.3. Переход к многоэлектронным атомам
Мы уже пытались рассматривать водородоподобные атомы с одним валентным электроном путем введения эффективного заряда ядра. Подробное изложение этого вопроса не является целью данного введения. Однако здесь будут коротко описаны два метода расчетов. Пер вый метод используется в модели самосогласованного поля (ССП) и в вариационной технике. Второй метод исполь зует подход (гл. 4), основанный на теории возмущений. Запишем уравнение Шредингера в форме
т = Е У ,
где Y зависит от координат нескольких электронов, а оператор Гамильтона Н имеет вид
3. Многоэлектронные атомы: модель CCI7 |
191 |
Складываем и вычитаем член «экранирования» |
U0 — |
= — 'EiUJfi), который в большинстве случаев должен i
быть по возможности равен —U. Тогда можно U = £/„ + + U1 расчленить на аддитивный электронный член U0 и небольшой возмущающий член U^. Теперь в первом при ближении можно записать Y как антисимметричное про изведение одноэлектронных функций. Эта исходная идея возмущений лежит в основе разработанной Слейтером теории атомов и ионов, которая особенно важна для ио нов тяжелых металлов и их комплексов.
3. М Н О Г О Э Л Е К Т Р О Н Н Ы Е А Т О М Ы : М О Д Е Л Ь СС П
Если удалить электрон i из многоэлектронного атома, его потенциальная энергия £/г выражается как
1Ф1
где rt — расстояние между ядром |
и t-м электроном, а |
гij — расстояние между i-u и /-м |
электронами. Второй |
член, представляющий взаимодействие электронов, при водит к тому, что в данном случае в уравнении Шредингера нельзя разделить переменные. Выход из этой ситуа ции, чрезвычайно важный для квантовой химии, состоит в том, что переходят к уравнению Шредингера с «настоя щим» оператором Гамильтона, который включает ука занное выше выражение Uu но, несмотря на это, исполь зуют произведения одноэлектронных функций (модель независимых частиц). Тогда в виде решений получают функции, которые снова можно интерпретировать как атомные орбитали. В упрощенной форме эти орбитали содержат наряду с зарядом ядра константу экранирова ния, которая рассчитывается для разных атомов по про стым правилам (слейтеровские орбитали) (см. ниже).
Ниже коротко рассмотрен ход такого расчета. Мы отрываем один электрон из «-электронной системы и для этого электрона решаем уравнение Шредингера со сред ним потенциалом, создаваемым другими электронами. Мы не знаем его точно, поэтому сначала применяем для
192 Часть П. Введение в квантовохимические расчеты
его расчета одноэлектронные функции состояния. При этом для отдельного электрона получается функция со стояния в первом приближении. Эту функцию для элект рона 1 вместе с другими уже использованными п — 2
Рис. 2. Радиальная плотность в зависимости от г дл я i c i c e Li+, Na+, К+.
функциями применяют для расчета нового, «улучшен ного» потенциального поля. С этим последним далее решают уравнение Шредингера для другого отдельного электрона 2. Этот процесс продолжается до тех пор, пока не будет получен полный набор из п «улучшенных» функ
3. Многоэлектронные атомы: модель ССП |
193 |
ций. Затем снова начинают всю процедуру сначала, пока не будет получен аналогичный набор. И это повторяется до тех пор, пока новый набор функций уже не будет более отличаться от предыдущего на сколько-нибудь заметную величину. Этот случай называют самосогласованием, а итерационный метод — методом самосогласо ванного поля (методом ССП). В этом методе мы можем сохранить классификацию орбиталей на s, р и d и т. д. в предположении, что потенциал ССП сделан сферически
симметричным путем усреднения. Естественно, что такой подход также является приближением, требующим, смот ря по обстоятельствам, критической оценки. Например, он ближе к истине в случае тяжелых атомов, так как вследствие больших зарядов ядер на поле отдельного электрона при усреднении приходится меньший вес. На рис. 2 показаны радиальные функции плотности, полученные методом ССП по Хартри.
Выведем качественно выражение для потенциала в уравнении Шредингера для одного электрона в методе
ССП. Предположим, что |
электрон 1 |
временно находится |
в точке Р на расстоянии |
г1 от ядра |
(рис. 3). Рассмотрим |
электронное облако ф* электрона 2. Оно создает в не большом элементе объема dxz на расстоянии г12 от точки Р заряд величиной
p2-flfx = —еф2^т;
потенциальная энергия по отношению к заряду —е в точке Р равна
Г12
1 3 — 2 0 8
194 Часть If. Введение в квантовохимические расчеты
Вклад электронного облака ф2 в потенциальную энергию получают интегрированием по всему объему
Далее суммируют все п — 1 вклады этого рода:
складывают с потенциальной энергией я д р а ---- — и
усредняют по всем направлениям. Таким образом, полу чают исходное выражение для потенциала в методе ССП.
Метод ССП сопряжен с трудоемкими вычислениями, которые чрезмерно возрастают с увеличением п, однако эти трудности можно будет преодолеть с совершенствова нием вычислительной техники (см. разд. 5.4).
3.1. Слейтеровские функции
Слейтеру удалось показать, что построенные по про стым рецептам и приведенные ниже волновые функции очень хорошо аппроксимируют численные значения, по лученные для одного атома методом ССП. В табл. 4 с обозначает эффективный заряд ядра, называемый также орбитальной экспонентой с = Z' = Z9(M) = Z — о; о — постоянная экранирования, а N — постоянная норми ровки.
Неизвестные а или с определяют по правилам Слей тера.
а. Подразделение электронов на следующие группы: Is
2s, 2р 3s, 3р
Ы
4s, 4р
4d
4/
5s, 5р
3. Многоэлектронные |
атомы: модель |
ССП |
195 |
||
|
|
|
|
Таблица 4 |
|
Слейтеровские |
функции <р( n l ) |
и постоянные |
нормировки N |
i |
|
Ф (Is) = |
N^e-*' |
СГ |
|
|
|
|
|
|
|
|
|
ф (2s) = |
N^re |
2 |
|
|
|
|
|
СГ |
|
|
|
Ф (3s) = |
N ^ e |
3 |
|
|
|
|
|
СГ |
|
|
|
Ф (2px) = Nwxe |
2 |
|
|
|
|
|
|
СГ |
|
|
|
Ф (3px) = |
N3Pxre |
3 |
|
|
|
|
—СГ |
|
|
|
|
'l5(3dxtJ) = |
N3(i2xye |
2 |
СГ |
|
|
|
|
|
|
|
|
Ф (3^2 _ fr2) = |
JV3d(JC3—y1)e |
3 |
|
|
3 Z2 _ r2 —i l Ф(3ф“ ^ у = ^ e 3
б. Для данного электрона а состоит из следующих вкладов:
1)электроны внешних оболочек, расположенных даль ше от ядра, чем оболочка рассматриваемого электрона, не дают вклада;
2)для каждого электрона в той же группе, что й рас сматриваемый электрон, устанавливается вклад 0,35 (за
исключением ls-группы, где он равен 0,3);
3)ближайшая внутренняя оболочка: каждый элект рон дает вклад 0,85 (если рассматриваемый электрон принадлежит к s- или р-группе) или 1,00 (если рассмат риваемый электрон принадлежит к d- или /-группе);
4)следующие внутренние оболочки: все электроны дают вклад 1,00.
Слейтеровские функции, построенные таким образом, часто используют при квантовохимических расчетах; в последнее время применяют также (см. разд. 5.5.4) функции Гаусса, имеющие форму e~ari, так как из их
13*