Файл: Хабердитцл, В. Строение материи и химическая связь.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 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*