конспект лекций - Кафедра квантовой химии

advertisement
1
МНОГОЭЛЕКТРОННЫЙ АТОМ (конспект лекций)
1.
Принципы квантовой механики
Квантовая химия атомов, молекул, молекулярных ансамблей и кристаллов опирается на следующие постулаты квантовой механики:
1. Каждое состояние системы n частиц полностью описывается функцией координат частиц xi и времени t
Ψ(x1, x2, …, xn, t) ≡ Ψ({x},t), называемой волновой функцией. Волновая функция существует во всем интервале
изменения переменных; она непрерывна, конечна и однозначна. Выражение Ψ*({x},t)Ψ ({x},t)dx имеет смысл
вероятности того, что в момент времени t i-я частица находится в интервале координат от xi до xi+dxi. Если
справедливо выражение
∞
∫ Ψ ({x}, t)Ψ ({x}, t)dx dx ,...,dx
*
1
2
n
= 1,
(1)
-∞
то говорят, что волновые функции нормированы на единицу. Поскольку физический смысл имеет лишь плотность вероятности Ψ*Ψ, то волновая функция определена с точностью до произвольного фазового множителя
типа eiα. Ψ* - функция комплексно сопряженная с Ψ.
2. Каждой доступной измерению величине А в любом из возможных состояний соответствует линейный
эрмитов оператор А. Оператор обозначается символом, определяющим математическую операцию, с помощью которой из одной функции получается друга. В квантовой химии используются операторы, которым отвечает уравнение типа
Аf = af,
(2)
где а − в общем случае комплексное число, называемое собственным значением оператора А; f называется
собственной функцией оператора А.
Оператор, обладающий свойством
∫f
*
1
(x) Af 2 ( x) dx = ∫ f 2 ( x) A* f1* ( x)dx
(3)
называется эрмитовым; собственные значения эрмитовых операторов −
действительные числа, а их собственные функции образуют полный ортонормированный набор, т.е.
∫f
*
i
(x)f j (x)dx = δ ij
=
1, если i = j
0, если i ≠ j.
(4)
А * – оператор, комплексно сопряженный с А; δij – символ Кронекера.
Действуя на волновую функцию, оператор превращает ее в другую волновую функцию. Иными словами,
действие оператора переводит систему в другое состояние; в частном случае система может остаться в том же
состоянии.
Таблица 1
Операторы основных физических величин
Переменная
Обозначение
Обозначение
Производимая
переменной
оператора
операция
2
Координата
r
r
Умножение на r
Момент
p
p
⎛ ∂
∂
∂ ⎞
− i=⎜⎜
+
+ ⎟⎟
⎝ ∂x ∂y ∂z ⎠
T
T
∂2
∂2
=2 ⎛ ∂2
⎜
−
+
+
2m ⎜⎝ ∂x 2 ∂y 2 ∂z 2
V(r)
V(r)
Умножение на V(r)
H
∂2
∂2 ⎞
=2 ⎛ ∂2
⎜
⎟ + V(r)
−
+
+
2m ⎜⎝ ∂x 2 ∂y 2 ∂z 2 ⎟⎠
Кинетическая
энергия
Потенциальная
энергия
Полная
E
энергия
⎞
⎟
⎟
⎠
Между операторами квантовой механики сохраняются те же соответствия, какие имеют место в классической
механике между динамическими переменными (координата, импульс, энергия и т.д.): это утверждение составляет принцип соответствия.
3. Независящая от времени волновая функция удовлетворяет стационарному уравнению Шредингера:
HΨ = ЕΨ.
(5)
Эрмитов оператор полной энергии системы (гамильтониан) H=T+V есть сумма оператора кинетической энергии всех частиц системы Т и оператора их потенциальной энергии V; Е − полная энергия системы.
Атомы, молекулы и кристаллы состоят из положительных ядер и отрицательных электронов, потенциальная
энергия которых определяется кулоновским взаимодействием. Операторы кинетической энергии системы, содержащей М ядер,
TЯ ( R ) и N электронов TЭ (r ) выглядят следующим образом:
TЯ (R ) = −
TЭ (r ) = −
где
Ma
=2
2
M
1
∑M
a
∇ a2 ,
(6)
a
=2 N 2
∑ ∇i ,
2m i
(7)
− масса ядра a; m − масса электрона; ∇ 2 = ∂
2
∂x 2
+
∂2
∂y 2
+
∂ 2 − оператор Лапласа (лапласиан), ħ =
∂z 2
h/2π = 1,0545·10-34 Дж·с, где h – постоянная Планка.
Дифференцирование в уравнении (6) ведется по координатам ядер
R a , а в (7) − по координатам элек-
тронов ri.
Вид операторов потенциальной энергии (в системе СИ) следующий:
M
M
a
b
VЯЯ (R) = ∑∑
3
2
Za Z b e
, a ≠ b,
4πε 0 R ab
(8)
Z ae2
,
VЭЯ (R, r) = −∑∑
a
i 4πε 0 rai
M
N
N
VЭЭ (r) = ∑∑
i
где
j
N
e2
4πε 0rij
,
i≠j
(9)
,
(10)
Za и Z b − атомный номер элемента, e − заряд электрона, R ab − расстояние между ядрами, rai
− рас-
стояние между ядром а и электроном i, rij − расстояние между электронами i и j, ε0 − диэлектрическая проницаемость вакуума. Оператор (8) описывает отталкивание ядер, (9) − притяжение электронов к ядрам, (10) − отталкивание электронов.
Строго говоря, в гамильтониане следует учесть релятивистские эффекты, обусловленные близкой к световой скоростью движения электронов в низкоэнергетических состояниях (электронов остова), спинорбитальным взаимодействием и малыми флуктуациями движущихся электронов относительно средних позиций. Однако в хорошем приближении их можно игнорировать при рассмотрении многих квантовохимических
задач.
Зависящая от времени волновая функция удовлетворяет нестационарному уравнению Шредингера
HΨ = i =
∂Ψ
∂t
,
(11)
4. Все одинаковые частицы тождественны. Именно поэтому можно говорить о неразличимости электронов:
замена одного из них другим не может быть обнаружена экспериментально.
5. Значения величины А, которые могут быть измерены, являются собственными значениями аi уравнения
на собственные значения
АΨi = аiΨi ,
(12)
где собственные функции Ψi есть волновые функции, описывающие возможные состояния системы, в которых
проводятся измерения. Это означает, что решение уравнения Шредингера (5) есть не что иное, как решение задачи на собственные значения для оператора полной энергии системы Н. Набор (спектр) собственных значений
Еi и набор собственных функций Ψi гамильтониана полностью характеризуют систему
HΨ i = E i Ψ i ,
E 0 ≤ E1 ≤ E 2 ≤ ... ≤ E n .
(13)
6. Среднее значение величины А для системы, находящейся в состоянии i, определяется выражением
a i = ∫ Ψ *i (x)AΨ i (x)dx
(предполагается, что волновые функции
Ψ *i (x)
(14)
ортонормированы). Если же система находится в нескольких
состояниях, то
a = ∑ Wi ∫ Ψ*i (x)AΨ i (x)dx ,
i
(15)
4
где
Wi
− вероятность пребывания системы в состоянии i, причем
∑ Wi = 1 .
i
Соотношение (15), выражающее принцип суперпозиции состояний. Определение измеряемых характеристик системы с помощью волновых функций производится с использованием выражений (14) и (15).
7. Координату и импульс частицы в любом состоянии одновременно определить точно невозможно
(принцип неопределенности). Количественно этот принцип записывается следующим образом:
ΔpxΔx≥ ħ,
(16)
где px – проекция импульса частицы на ось х, Δ обозначает неопределенность значений величины.
2.
Вариационный принцип. Решение уравнения Шредингера
Согласно принципам квантовой механики для изучения электронных свойств химических систем необходимо решить уравнение Шредингера для всех возможных электронных состояний. Любая система стремится
занять состояние с минимальной энергией. Поэтому решения ищут с помощью вариационного принципа, минимизируя энергию системы и определяя функции, максимально близкие к собственным функциям оператора Н.
Вариационный принцип утверждает, что среднее значение энергии Еi любого из возможных i состояний системы, вычисленное с приближенной волновой функцией, не может быть меньше нижнего собственного значения
Е0 оператора Н. Покажем это.
Среднее значение оператора Н для некоторой приближенной волновой функции Ψ, нормированной на 1,
равно
E = ∫ Ψ * (x)HΨ (x)dx ≡ Ψ | H | Ψ
.
(17)
которое сокращенно записывается как E = Ψ | H | Ψ . Представим Ψ в виде разложения по собственным
функциям оператора Н:
n
Ψ = ∑ ci Ψ i .
(18)
i
Это всегда возможно, поскольку функции Ψi составляют полную ортонормированную систему. В силу
этого
n
*
*
Ψ
(
x
)
Ψ
(
x
)
d
x
=
c
c
Ψ
∑
i
i
i
∫
∫ Ψi dx = 1,
*
(19)
i
откуда следует:
n
∑
i
c*ic i =
n
∑ ci
2
=1 .
(20)
i
Подставляя (20) в (19), имеем:
n
E = ∫ Ψ (x)HΨ (x)dx = ∑ ci
*
i
2
∫ Ψ ( x ) H Ψ ( x ) dx ,
*
i
i
(21)
5
где Еi − энергия i-го состояния.
С другой стороны,
n
2
E = ∑ ci E i ,
(22)
i
откуда
n
E ≥ E 0 ∑ ci
2
= E0 ,
(23)
i
где Е0 − нижнее (наименьшее) собственное значение оператора Н, так как энергии возбужденных состояний
выше, чем энергия самого низкого по энергии состояния, называемого основным. Что и требовалось доказать.
Итак, чтобы решить уравнение Шредингера, нужно минимизировать выражение для энергии (17), т.е.
подобрать такие волновые функции, для которых энергия будет минимальна. Это можно сделать, если волновая
функция включает некоторые переменные параметры, изменяя которые можно обеспечить минимум. Например, принцип суперпозиции позволяет записать волновую функцию в виде линейной комбинации
n
Ψ = ∑ c i ϕi ,
(24)
i
где ϕi − n выбранных из некоторых соображений функций, часто называемых базисными (например, для атомов
− это могут быть атомные орбитали, которые мы рассмотрим далее), ci − переменные комплексные параметры.
Чтобы найти минимум энергии (17), нужно, чтобы все первые производные энергии по параметрам ci
равнялись нулю, а вторые производные были положительны. Однако обычно ограничиваются равенством нулю
первых производных, что отвечает условию стационарности энергии:
∂Ε ∂Ε
∂Ε
=
= ... =
=0,
∂c1 ∂c2
∂cn
(25)
∂Ε ∂Ε
∂Ε
=
=
...
=
=0
∂c1* ∂c2*
∂cn*
На языке вариационного исчисления, обычно используемого в квантовой химии, условие (25) эквивалентно требованию обращения в нуль первой вариации (вариация – бесконечно малое изменение функции):
δ ∫ Ψ* (x)HΨ (x)dx = 0 .
(26)
Из условия стационарности (25) должен следовать набор уравнений, решая которые совместно, можно
найти параметры ci . При этом следует принять во внимание, что параметры ci не независимы, а связаны в силу
ортономированности функции Ψ соотношением:
∫Ψ
*
n n
(x)Ψ(x)dx = ∑∑
i
j
c*ic j ϕ*i (x)ϕ j (x)dx =
∫
n n
∑∑ c*i c jSij = 1 .
i
j
(27)
(27)
6
Интеграл
Sij = ∫ ϕ*i (x)ϕ j (x)dx называется интегралом перекрывания функций ϕi и ϕj.
При минимизации с учетом ограничений, каковым является условие (25), в математике используется метод неопределенных множителей Лагранжа. Вводя такой множитель Е, представим уравнение для определения параметров ci в виде:
[ ∫ Ψ* (x)HΨ (x)dx − E(∑∑ ∫ ϕi* (x)ϕ j (x)dx − 1) ] =
i
j
δ
=δ [
n
n
i
j
∑∑ c c ∫ ϕ
*
i j
*
i
(x)Hϕ j (x)dx − E(Sij − 1)
(28)
]=0.
Теперь все параметры ci можно считать независимыми. Производя варьирование по ci и cj, имеем в компактной записи:
n
n
∑ дс ∑ с [H
*
i
j
i
ij
− ESij ] = 0,
n
n
j
i
∑ дс j ∑ с*i[Hij − ESij ] = 0,
где
i = 1,2,3,..., n
(29)
j
H ij = ∫ ϕ*i (x)Hϕi (x)dx
i = 1,2,3,..., n
− матричные элементы оператора Н в базисе функций
элементы матрицы интегралов перекрывания, вычисленной с тем же набором функций
⎡ H11
⎢H
H = ⎢ 21
⎢ .....
⎢
⎣H n1
H12
H 22
.....
H n2
..... H1n ⎤
⎡S11 S12
⎢S
S 22
..... H 2n ⎥
⎥ S = ⎢ 21
⎢..... .....
..... ..... ⎥ ,
⎢
⎥
..... H nn ⎦
⎣S n1 S n2
ϕi (x) , Sij −
ϕi ( x ) :
..... S1n ⎤
..... S 2n ⎥
⎥
..... ..... ⎥
⎥
..... S nn ⎦
(30)
Поскольку все вариации δсi в (29) независимы, матричные уравнения справедливы лишь, если коэффициенты при вариациях равны нулю. Каждое матричное уравнение в (29) получается из другого операцией комплексного сопряжения, поэтому достаточно рассматривать только одно из них:
n
∑ с [H
j
j
ij
− ESij ] = 0,
i = 1,2,3,..., n
(31)
7
Полученная система однородных линейных уравнений позволяет найти параметры ci, обеспечивающие
минимум функционала (17). Чтобы ее решить, необходимо приравнять нулю определитель (детерминант) из
коэффициентов при ci:
H 11 − ES11
H 21 − ES 21
.....
H 12 − ES12
H 22 − ES 22
.....
..... H 1n − ES1n
..... H 2 n − ES 2 n
=0
.....
.....
H n1 − ES n1
H n 2 − ES n 2
..... H nn − ES nn
(32)
или в других обозначениях
H ij − ES ij = 0
(33)
Уравнение (33) впервые возникло в задачах классической механики о вековых (секулярных) возмущениях планетарных орбит. Отсюда его название − вековое или секулярное уравнение. При разложении определителя
получается многочлен n-й степени по Е, значит, вековое уравнение имеет n корней − n различных значений Е.
Подставляя их в (31), можно найти набор параметров ci . Величины Еi играют роль энергий состояний системы.
Чтобы найти волновую функцию основного состояния, нужно взять наименьшее из полученных значений энергии и, подставив его в (31), найти коэффициенты сi, а затем и волновую функцию
Ψ
по формуле (24).
Волновые функции возбужденных состояний ищут аналогичным образом, учитывая, что они должны быть ортогональны друг другу.
Точное решение уравнения Шредингера возможно только для атома водорода и водородоподобных ионов. Для многоэлектронных атомов необходимо применять приближения, рассмотренные ниже.
3.
Приближение независимых частиц. Одноэлектронная модель
Рассмотрим, как можно найти волновые функции и уровни энергии неподвижного N-электронного атома, максимально близкие к точным. Выберем начало координат на его ядре с зарядом
Ze
. Гамильтониан в
этом случае имеет вид:
(34)
H = Tэ (r ) + Vэя (r, R ) + Vээ (r ) = −
2
2
2
1 N N e
= N 2 N Ze
∇i − ∑
+ ∑∑
∑
2m i
2 i j 4πε 0 rij
i 4πε 0 ri
Вначале просто исключим из (34) оператор энергии межэлектронного взаимодействия Vээ. Многоэлектронное уравнение Шредингера в этом случае распадается на систему из N одноэлектронных уравнений
hi χ i (ri ) = ε i χ i (ri )
i = 1,2,3,...,N
(35)
c одноэлектронными гамильтонианами
hi = −
=2 2
Ze2
∇i −
2m
4πε 0 ri
.
(36)
8
Согласно (35) электрон i описывается волновой функцией
χ i (ri )
и имеет энергию
ε i . Это означает,
что поведение каждого электрона не зависит от поведения остальных электронов и описывается некоторой волновой функцией подобно единственному электрону в атоме водорода. В этом состоит суть приближения независимых частиц. Решения одноэлектронных уравнений (35)
χi (ri )
называются одноэлектронными волновыми
функциями или орбиталями (в атоме − атомными орбиталями, в молекуле − молекулярными, в кристалле −
кристаллическими).
Полный гамильтониан атома в этом приближении есть просто сумма одноэлектронных гамильтонианов:
N
H = ∑ hi
1
.
(37)
Собственные функции Н представляют собой произведение атомных орбиталей, заселенных электронами:
Ψ = χ1(r1 ) ⋅ χ1(r2 ) ⋅ χ3 (r3 ) ⋅⋅⋅ χN (rN ) ,
(38)
а энергия атома является суммой индивидуальных орбитальных энергий:
Ψ H Ψ = ε1 + ε 2 + ε 3 + ... + ε N = E.
(39)
Приближенная волновая функция вида (38) называется волновой функцией Хартри.
Гамильтониан (36) является чрезмерно упрощенным: энергия электрон-электронного отталкивания не
мала и пренебрегать ею нельзя.
4.
Метод самосогласованного поля
Более близкие к истинным решения получают с помощью метода самосогласованного поля (ССП), предложенного Хартри. В методе ССП межэлектронным отталкиванием не пренебрегают, но действие на данный
электрон всех остальных электронов заменяют средним полем, приближенно воспроизводящим их суммарное
действие; последнее зависит от координат только рассматриваемого электрона. Это дает возможность разделить в сферической системе координат переменные в уравнении Шредингера. С формальной точки зрения это
достигается следующим образом. Одноэлектронный гамильтониан записывают в виде:
h ССП
=−
i
N
Ze 2
=2 2
∇i −
+∑
2m
4πε0 ri
j
⎛ e2 ⎞
⎜
⎟ , j = 1,2,3,..., (N - 1)
⎜ 4πε0 rij ⎟
⎝
⎠ j
(40)
Последнее слагаемое описывает отталкивание между электронами i и j, усредненное по всем положениям электрона j и, следовательно, зависящее только от координат электрона i . Последствия этого состоят в следующем.
Рассмотрим гамильтониан
Н=
n
.
∑ h ССП
i
(41)
i
Его собственные функции (функции Хартри) имеют вид орбитальных произведений:
Ψ = χ1(r1 ) ⋅ χ1 (r2 ) ⋅ χ3 (r3 ) ⋅⋅⋅ χN (rN ) .
(42)
Собственные значения Н представляются суммой собственных значений
h
9
ССП
:
i
.
n
n
i
i
E' = ∑ е i = ∑ Ψ | h ССП
|Ψ
i
(43)
Энергия εi есть сумма кинетической энергии i-го электрона, потенциальной энергии его притяжения к
ядру и средней потенциальной энергии его отталкивания от остальных электронов. Следовательно, Е′ есть сумма кинетической энергии всех электронов, потенциальной энергии их притяжения к ядру и удвоенной потенциальной энергии их усредненного отталкивания от остальных электронов. Удвоение возникло потому, что отталкивание между электронами i и j учтено дважды: как среднее по j в
h ССП
i
(40) и среднее по i в
hССП
. С учеj
том этого, полная энергия атома равна:
⎛ e2 ⎞
⎜
⎟
⎜ 4πε r ⎟
0 ij ⎠
⎝
1 n n
∑∑
2 i j
E = E'−
i ≠ j.
(44)
i, j
Соответственно, гамильтониан атома должен иметь вид:
H=H−
1 n n
∑∑
2 i j
⎛ e2 ⎞
⎜
⎟
⎜ 4πε r ⎟
0
ij
⎝
⎠
i ≠ j.
,
(45)
i, j
Таким образом, необходимо решить систему одноэлектронных уравнений с гамильтонианом (45), включающим усредненное межэлектронное взаимодействие – систему уравнений Хартри. Для этого нужно построить набор операторов
h ССП
, для чего следует прежде рассчитать усредненные величины
i
⎛ e2 ⎞
⎟ . Как
⎜
⎜ 4πε r ⎟
0 ij ⎠
⎝
j
это сделать? Вероятность того, что электрон j с волновой функцией χj(rj) находится в бесконечно малом объеме
dvj, равна
χ 2j dv j
(рис. 1).
Рис. 1. Иллюстрация к вычислению усредненного межэлектронного кулоновского взаимодействия. Электрон i в
точке, последовательно “пробегающей” все положения в пространстве
Значит, отталкивание электрона i, усредненное по всем положениям
электрона j, равно:
⎛ e2
⎜⎜
⎝ 4πε 0 rij
⎞
⎟⎟
⎠
⎛ χ j 2 (rj ) ⎞
dv j .
=e ∫⎜
⎜ 4πε r ⎟⎟
0 ij ⎠
−∞ ⎝
∞
2
j
Однако, чтобы вычислить этот интеграл, волновые функции
(46)
χ j (rj )
должны уже быть известны! Это
противоречие преодолевается следующим образом. Сначала задаются некоторым набором N одноэлектронных
функций, максимально близких к правильным
10
χ (rj ) (позже мы увидим, что сделать это легко). С их по0
j
мощью вычисляют интеграл (46) и строят оператор
(h 0i ) ССП . Затем решают набор одноэлектронных урав-
нений Хартри, возникающий из условия минимума среднего значения гамильтониана (40), вычисляемого с волновой функцией Хартри (42):
N ∞ ⎛ [ χ 0 (r )]2 ⎞
⎧⎪ = 2 2
⎫⎪ 1
Ze 2
j
j
2
1
e
dv
−
∇
−
+
⎜
⎟
⎨
⎬ χ i (ri ) = ε i χ i (ri ), i = 1,2,3,...,N .
∑
i
j
∫
⎜
⎟
4πε o ri
j ≠ i −∞ ⎝ 4πε o rij ⎠
⎩⎪ 2m
⎭⎪
Полученные решения
χ i1 (ri )
(47)
используют, чтобы построить "исправленный" оператор
(h 1i ) ССП
решают ту же систему уравнений, но теперь – с
(h 1i ) ССП , вновь
и т.д. – до тех пор, пока получаемые собственные
значения уравнений Хартри (т.е. энергии состояний) будут отличаться от полученных на предыдущей итерации
лишь на очень незначительную величину (~ 10-6 а.е.). Этот процесс называется самосогласованием, а результирующее поле, создающее усредненный потенциал в (40), называется самосогласованным полем – отсюда и название метода. Отметим, что сходимость метода не гарантируется теорией, но, как правило, достигается на
практике. Существуют довольно хорошо разработанные методы, которые позволяют обойти встречающиеся
здесь иногда затруднения.
Одноэлектронное приближение и метод ССП на первый взгляд кажутся довольно грубыми, однако это не
так. Дело в том, что быстро движущийся электрон “чувствует” скорее среднее эффективное поле остальных
частиц, чем реагирует на мгновенные изменения их позиций. Принципиально важно, что самосогласованные
решения удовлетворяют вариационному принципу, т.е. приводят к средним значениям энергии состояний, которые не ниже, чем точные энергии.
5. Приближение центрального поля
N
Потенциал
∑
j
⎛ e2 ⎞
⎜
⎟
⎜ 4πε r ⎟
0 ij ⎠
⎝
в (40) только в частных случаях (положительные одноэлектронные иоj
ны, атомы инертных газов, атомы N, P и т.д.) является сферически симметричным, т.е. не зависит от углов θ и ϕ
в сферической системе координат. Опыт, однако, показывает, что учет асферичности электронного распределения не улучшает заметно результат расчета. Поэтому обычно используют дополнительно усредненный по всем
направлениям потенциал в (40), интегрируя его по углам θ и ϕ:
N
∑
j
⎛ e2
⎜⎜
⎝ 4πε 0 rij
⎞
⎟⎟
⎠
j ,θ ,ϕ
e2
=
4π
⎛ χ j 2 (rj ) ⎞ 2
∫0 ⎜⎜ 4πε 0 rij ⎟⎟ rj sinθ jdθ jdϕ jdrj .
⎝
⎠
π 2π
∑∫
i
0
(48)
Вводимое таким образом приближение центрального поля имеет весьма важные последствия: оно позволяет рассматривать ССП-решения для любого атома как модифицированные решения для одноэлектронного
11
Ze
водородоподобного атома с потенциалом
. В этом приближении потенциальная энергия зависит
4πε 0 r
только от расстояния до ядра, т.е. сила притяжения к ядру носит центральный характер. Тогда переменные в
уравнении Шредингера в сферических координатах разделяются, и волновые функции (атомные орбитали)
имеют вид:
χ(r) = N(n,l) Rnl (r)Ylm (θ, ϕ),
(49)
где N(n,l) – нормировочный множитель, Rnl (r) и Ylm (θ, ϕ) – радиальная и угловая части волновой функции, n, l
и m – главное, орбитальное и магнитное квантовые числа, соответственно.
Рассмотрим свойства атомных орбиталей (АО) подробнее.
6. Атомные орбитали и их характеристики
Точное значение нормированной радиальной функции Rn,l для водородоподобного атома дается выражением:
1
Zr
⎛ 2 Zr ⎞ − na 0 2l +1 ⎛ 2 Zr ⎞
⎟⎟ ,
⎟⎟ e L n +l ⎜⎜
R nl (r ) = − N Znl ⎜⎜
na
n
a
⎝ 0⎠
⎝ 0⎠
где
(50)
N Znl – нормировочный множитель, зависящий от Z, n и l; L2nl++l1 ⎛⎜⎜ 2Zr ⎞⎟⎟
⎝ na 0 ⎠
Лягерра, a 0 =
– присоединенные полиномы
4πε 0 = 2
= 0.529·10–10 м – радиус Бора, l = 0, 1, 2, 3,…, n - 1.
2
em
Выражение (50) есть решение радиального уравнения Шредингера
∂ 2 R 2 ∂R 2Z l (l + 1)
2mE
,
+
+
−
=
−
r ∂r ra 0
∂r 2
r2
=2
(51)
конкретный вид которого возник после разделения переменных в сферических координатах. Несколько нормированных радиальных волновых функций, описывающих основное (n = 1) и первые возбужденные (n = 2, n = 3)
состояния, приведены в табл. 2, а их зависимость от расстояния изображена на рис. 2 и 3.
12
Таблица 2
Радиальные нормированные функции водородоподобных атомов Rnl (r)
n
1
3
2
1
0
1
3
2
Zr
⎛ Z ⎞ − a0
2⎜⎜ ⎟⎟ e
⎝ a0 ⎠
0
0
2
Rnl (r)
l
1 ⎛Z⎞
⎜ ⎟
2 ⎜⎝ a 0 ⎟⎠
3
2
Zr
⎛
Zr ⎞ − 2a 0
⎜⎜1 −
⎟⎟e
2a
0 ⎠
⎝
5
2
Zr
1 ⎛ Z ⎞ − 2a 0
⎜⎜ ⎟⎟ re
2 6 ⎝ a0 ⎠
2 ⎛Z⎞
⎜⎜ ⎟⎟
3 3 ⎝ a0 ⎠
3
2
Zr
⎛ 2 Zr 2 Z 2 r 2 ⎞ − 3a 0
⎟e
⎜⎜1 −
+
2 ⎟
3
a
27a 0 ⎠
0
⎝
8 ⎛Z⎞
⎜ ⎟
27 6 ⎜⎝ a 0 ⎟⎠
3
2
Zr
⎛ Zr Z 2 r 2 ⎞ − 3a 0
⎟e
⎜ −
2 ⎟
⎜a
6
a
0 ⎠
⎝ 0
7
2
Zr
4 ⎛ Z ⎞ 2 − 3a 0
⎜⎜ ⎟⎟ r e
81 30 ⎝ a 0 ⎠
Свойства радиальных функций.
1) Как следствие свойств полиномов Лягерра радиальные функции с различными n и l ортогональны.
2) Имеются точки (поверхности), где функции Rnl (r) обращаются в нуль; они называются узловыми точками (поверхностями) или просто узлами. Вероятность найти электрон в узле равна нулю. Радиальные функции с
(n=1, l=0), (n=2, l=1), (n=3, l=2) и т.д. не имеют узловых точек; функции с (n=2, l=0), (n=3, l=1) и т.д. имеют одну
узловую точку; функция с (n=3, l=0) – две узловые точки. Таким образом, число узлов радиальной функции
равно n-l-1.
13
Узел
Узел
r
a0
r
a0
а
б
r
a0
в
Рис. 2. Радиальные составляющие 1s (а), 2s (б), 3s (в) орбиталей атома водорода
3) Вероятность нахождения электрона в пространственном слое между значениями r и r+dr равна:
2π π
Pnl (r) dr = R nl (r) r dr ∫
2
2
∫
2
2
Ylm ( θ , ϕ ) sin θ d θ d ϕ = R nl (r) r 2 dr (52)
0 0
(интеграл по углам θ и ϕ равен единице). Функция Pnl (r), определяющая плотность вероятности нахождения
электрона в слое dr на расстоянии r от ядра, называется радиальной функцией распределения (рис. 4). Приравнивая нулю производную Pnl по r, можно найти наиболее вероятное положение электрона на соответствующей
орбитали. Для основного состояния атома водорода это расстояние равно радиусу Бора
a0 .
4) Вблизи ядра электрон-ядерный потенциал Vэя (9) становится неопределенным из-за стремления знаменателя к нулю. Чтобы волновая функция на ядре была конечна (как это имеет место для функций s-типа), необходимо, чтобы ее радиальная часть удовлетворяла асимптотическому условию
⎛ Z⎞
⎛ ∂R ⎞
= −⎜⎜ ⎟⎟R r → 0 .
⎟
⎜
⎝ ∂r ⎠ r → 0
⎝ a0 ⎠
(53)
5) На больших расстояниях от ядра атомная орбиталь зависит от расстояния как
R(r) ~ exp [–
где I1 – первый потенциал ионизации.
2I 1
r],
(54)
14
Рис. 3. Радиальные составляющие
Рис. 4. Радиальные функции распределения атома во-
2р и 3р орбиталей атома водорода
дорода
Угловые функции Ylm (θ, ϕ) − собственные функции оператора квадрата углового момента L2 − описывают в сферических координатах (θ, ϕ) угловую зависимость вероятности нахождения электронов в центральном поле атома. Они представляют собой сферические гармоники:
m+ m
2
Ylm(θ,ϕ)=(−1)
⎡ (2l + 1) (1 − m )!⎤
⎢
⎥
⎣⎢ 4π (1 + m )!⎥⎦
где l = 0, 1, 2,..; m = − l,...0,…+ l;
1
2
Pl
m
(cos θ) exp (imϕ) , (55)
m
Pl (cos и) − присоединенные полиномы Лежандра.
Это комплексные ортонормированные функции, из которых легко построить действительные комбинации, оставляющие АО собственными функциями того же одноэлектронного уравнения:
ylm+ = (
1
)[(−1)m Ylm + Yl-m],
2
ylm- = − (
(56)
i
)[(−1)m Ylm − Yl-m], l = 0, 1,2, ...; m = ±1, ±2,…
2
В табл. 3 и на рис. 5-7 представлены действительные угловые функции. Они имеют простую интерпретацию в декартовых координатах. Для них, также как и для радиальных функций, характерно наличие узлов и узловых плоскостей, число которых равно l. Узлы полной атомной орбитали χ(r) (49) определяются узлами её радиальной и угловой составляющих (рис. 6).
При классификации электронных состояний водородоподобного атома придерживаются следующих
представлений. Главное квантовое число n характеризует энергию орбитали. Орбитальное квантовое число l
характеризует угловую зависимость орбитали (орбитальный момент). Для каждого l приняты свои обозначения
15
(табл. 4). Важно понимать, что как результат приближения центрального поля угловая зависимость АО c
одним и тем же l для всех атомов одинакова.
Таблица 3
Угловые части волновой функции атома, обладающего центральным полем
l
Линейная
yl m
m±
Обозначение
Z
комбинация
1
4π
Y
0
0
-
s
1
0
6
cos θ
2 2π
-
pz
1
−
6
sin θ sin ϕ
2 2π
1
+
6
sin θ cos ϕ
2 2π
2
0
15
(3 cos 2 θ − 1)
4 π
-
dz2
2
+
15
sin 2θ cos ϕ
4 π
1
(Y2 −1 + Y21 )
2
dxz
2
−
15
sin 2θ sin ϕ
4 π
2
+
15
sin 2 θ cos 2ϕ
4 π
2
−
15
sin 2 θ sin 2ϕ
4 π
i
(Y11 + Y1−1 )
2
−
−
1
(Y11 − Y1−1 )
2
i
(Y21 − Y2 −1 )
2
1
(Y22 + Y2− 2 )
2
−
i
(Y22 − Y2 − 2 )
2
X
py
px
dyz
dx 2 −y2
dxy
Таблица 4
Обозначения орбиталей с различными угловыми зависимостями
l
0
1
2
3
4
5
Символ
s
р
d
f
g
h
16
a
б
Рис. 5. Графическое представление угловых частей 2р орбиталей атома водорода: а) − с помощью поверхности, ограничивающей 90% электронной плотности (2pZ); б) − с помощью изолиний (2рX)
17
py
pz
d yz
d xz
px
а
б
dz2
d xy
dx 2 −y2
в
fx(x 2 −y2
f xyz
fz(x 2 − y2 )
f y(x 2 − y2 )
f yz 2
fz3
f xz 2
г
Рис. 6. Угловые части орбиталей атома водорода: а) 1s, б) 2p, в) 3d, г) 4f
а
б
18
Рис. 7. Контурные карты электронной плотности орбиталей: а) − 2pZ орбитали атома углерода; б) − 3pZ орбитали атома хлора. Следы радиальных и угловых узловых поверхностей показаны штриховыми линиями
7. Антисимметричность многоэлектронной волновой функции
До этого момента мы игнорировали тот факт, что электроны обладают собственным моментом количества движения (спином), который в единицах ħ равен ±
1
. В квантово-химической литературе два возможных
2
значения спина электрона α и β обозначаются символами ↑ и ↓ . Часто говорят "спин вверх" или "спин вниз",
указывая разные значения спина электрона. Чтобы учесть наличие спина, в аргумент одноэлектронных функций вводят спиновую переменную s, а функцию χi (xi) называют спин-орбиталью (xi = ri , si).
Из-за наличия спина имеются специальные ограничения на электронную волновую функцию, известные
как принцип Паули:
1.
Электронная волновая функция системы из нескольких электронов должна быть антисимметричной (ме-
нять знак) относительно обмена положениями и спинами пары из любых двух электронов i и j:
Ψ (x 1 ,...x i ,...x j ,...) = − Ψ (x 1 ,...x j ,...x i ,...) .
2.
(57)
Никакие два электрона не могут быть в одном и том же квантовом состоянии (принцип исключения).
19
Чтобы понять, как нужно строить волновую функцию с учетом принципа Паули, рассмотрим двухэлектронный атом гелия. Две эквивалентных волновых функции Хартри для этой системы имеют вид:
Ψ 1 = χ 1 (x1 ) ⋅ χ 2 (x 2 )
Ψ 2 = χ 1 (x 2 ) ⋅ χ 2 (x1 )
(58)
.
Очевидно, что ни одна из них не является антисимметричной. Однако легко заметить, что связанная с ними
функция
Ψ=
антисимметрична (
1
2
1
[χ1 (x1 ) ⋅ χ 2 (x 2 )-χ1 (x 2 ) ⋅ χ 2 (x1 )]
2
(59)
– нормировочный коэффициент).
Если мы попытаемся поместить электроны 1 и 2 на одну и ту же спин-орбиталь χi, то получим:
Ψ=
1
[χ i (x1 ) ⋅ χ i (x 2 )-χi (x 2 ) ⋅ χi (x1 )] = 0 .
2
(60)
Значит, волновая функция в виде (59) удовлетворяет принципу исключения Паули. С математической точки
зрения волновая функция есть детерминант из спин-орбиталей χi(х), где i = 1 или 2:
Ψ=
1 χ1 (x1 ) χ 2 (x1 )
2 χ1 (x 2 ) χ 2 (x 2 )
(61)
.
Важным свойством детерминанта, известным из курса высшей математики, является то, что он меняет
знак при перестановке двух любых его столбцов или строк; величина его при этом остается неизменной. Это
как раз эквивалентно перемене мест двух электронов: электрон 1 перемещается с орбитали i на орбиталь j, а
электрон 2 – с орбитали j на орбиталь i. Если же два столбца или две строки детерминанта одинаковы (что эквивалентно занятию двумя электронами одной и той же орбитали), то детерминант равен нулю.
Принцип Паули диктует, что две спин-орбитали с одинаковыми пространственными частями (т.е. с одинаковыми квантовыми числами n, l, m) отличаются только спиновыми компонентами. Пренебрегая малым по
величине спин-орбитальным взаимодействием, каждую спин-орбиталь χi(x) можно представить в виде произведения пространственной орбитали ϕi(r) и спиновой функции η(s):
χi(x) = ϕi(r)η(s).
(62)
Рассмотрим атом Не. Из возможных вариантов распределения электронов по 1s и 2s АО атома Не, представленным на рис. 8, основному состоянию соответствует схема а. Волновая функция основного состояния
записывается через определитель следующим образом:
Ψa =
где функция
α (1) α (2)
1 ϕ1 (r1 )α (1) ϕ1 (r2 )α (2)
1
ϕ1 (r1 )ϕ1 (r2 )
=
β (1) β (2)
2 ϕ1 (r1 ) β (1) ϕ1 (r2 ) β (2)
2
,
(63)
ϕ1 (r1 ) - это 1S АО.
Таким образом, функция
Ψa
имеет симметричную пространственную часть и антисимметричную спи-
новую. Полный спин системы равен нулю, поэтому это состояние принято называть синглетным.
20
Рис. 8. Состояния атома Не: а – основное, б, в, г, д, е – возбуждённые
Заполнение электронами орбиталей любого атома происходит в порядке возрастания их энергий. При
этом обязательно выполняется требование принципа Паули и каждая орбиталь может быть заполнена не более
чем двумя электронами. В последнем случае электроны обязательно имеют противоположные спины и называются спаренными. Состояния и соответствующие им АО с одинаковыми квантовыми числами n и l заполняются
так, чтобы суммарный спин электронов был максимален (правило Хунда); другими словами, сначала заполняется максимальное число АО с разными m.
Для нейтральных атомов в основном состоянии последовательность заполнения АО следующая:
1s < 2s < 2p < 3s < 3p < 4s ~ 3d < 4p < 5s ~ 4d < 5p < 6s ~5d ~ 4f < 6p < 7s …
Это ведет к так называемому aufbau-принципу, который служит обоснованием таблицы Менделеева.
Следует подчеркнуть, что порядок заполнения атомных орбиталей обусловлен не отношением их энергий, а
требованием минимума полной энергии атома, которая отлична от суммы одноэлектронных энергий. Здесь
также играют роль спин-орбитальные и другие эффекты. Именно поэтому, в частности, 4s AO заполняются
электронами в атомах калия и кальция раньше, чем 3d АО.
В итоге электроны с одним и тем же значением главного квантового числа n образуют электронные слои.
Эти слои, в свою очередь, построены из оболочек (nl), заполненных электронами с одинаковым значением орбитального квантового числа l. Слои с n = 1, 2, 3, 4,... обозначают буквами K, L, M, N,...
8.
Детерминант Слейтера
Итак, представление многоэлектронной волновой функции в виде детерминанта обеспечивает ее правильные антисимметричные свойства. Кроме того, электроны неразличимы, и их перестановка не меняет свойства системы. Перестановка электронов для волновой функции в виде детерминанта эквивалентна перестановке
местами столбцов (строк), что лишь меняет знак детерминанта. Поскольку волновая функция в принципе определена с точностью до фазового множителя, перемена знака свойств системы не меняет. Приближенная многоэлектронная волновая функция, построенная из ортонормированных спин-орбиталей отдельных электронов, называется детерминантом Слейтера:
Ψ = (N!)
Множитель
( N!)
−
−
1
2
χ1 (x1 )
χ1 (x 2 )
χ 2 (x1 ) .... χ N (x1 )
χ 2 (x 2 ) .... χ N (x 2 )
....
....
....
....
χ1 (x N ) χ 2 (x N ) .... χ N (x N )
(64)
.
1
2 перед детерминантом обеспечивает условия нормировки, здесь N – число электронов.
21
Детерминант Слейтера является единственной функцией, обеспечивающей антисимметричность вол-
новой функции, записанной через орбитали (орбитальное приближение). Следовательно, он дает только одно
решение соответствующих одноэлектронных уравнений.
Хотя электроны неразличимы, в орбитальном приближении каждый электрон описывается "своей" волновой функцией. Системы, в которых все электроны занимают орбитали попарно, называются системами с закрытыми (замкнутыми) электронными оболочками. Для таких систем детерминант Слейтера состоит из дважды занятых электронами орбиталей, число которых равно половине числа электронов. Системы с нечетным
числом электронов называются системами с открытыми (незамкнутыми) оболочками.
9. Метод Хартри-Фока
Аппроксимация многоэлектронной волновой функции единственным детерминантом Слейтера (64) и использование при ее нахождении приближения самосогласованного поля приводят к методу Хартри-Фока
(ХФ). При этом исходное электронное уравнение Шредингера (21) путем довольно громоздких математических вычислений преобразуется в уравнение, где точный гамильтониан H (34) заменен оператором Фока (фокианом):
Fi = −
1 N 2 N Z
∑ ∇ i − ∑i r
2 i
i
N
+
N
⎡
∑ ∑ ⎢∫ χ
i
j
⎣⎢
j
(x j )
2
⎤
1
1
1
dx j − ∫ χ i* ( x j ) χ j ( x j ) dx j ⎥ .
x ij
x ij
2
⎦⎥
(65)
Здесь и далее мы используем принятую в квантовой химии атомную систему единиц: множитель
1/(4πε 0 ) опускается, m = 1, e = 1, ħ =1. Введение атомных единиц делает формулы менее громоздкими.
Различие между F и H состоит в том, что оператор кулоновского электронного взаимодействия
n
n
i
j
∑∑
⎛ e2 ⎞
⎜
⎟
⎜ 4πε r ⎟
0 ij ⎠
⎝
заменен в (65) оператором в квадратных скобках, описывающим взаимодействие каi, j
ждого электрона со средним полем всех остальных электронов с учетом требований принципа Паули. Из условия минимума энергии
Ψ |F|Ψ
возникает набор независимых уравнений для каждой одноэлектронной
орбитали – уравнений Хартри-Фока:
N ⎡⎛
⎤
⎞
⎞
2 1
1⎛
1
hi (xi )χi (xi ) + ∑⎢⎜ ∫ χ j (xj ) −1 dxj ⎟ χi (xi ) − ⎜ ∫ χ*j (xj ) −1 χi (xj )dxj ⎟ χ j (xi )⎥ = εi χi (xi ) . (66)
⎜
⎟
⎟
2 ⎜⎝
rij
rij
j =1 ⎢⎝
⎥⎦
⎠
⎠
⎣
Энергия электрона, находящегося на орбитали
на
χi
χ i , может быть получена умножением слева выражения (66)
и интегрированием по всему пространству:
N
1
е i = h ii + ∑ (J ij − K ij ) ,
2
j=1
(67)
hii = ∫ χ i* ( xi )h i χ i ( xi ) dxi ,
(68)
Здесь
22
1
J ij = ∫ χ i* (xi ) χ i* (xi ) χ j (x j ) χ j (x j )dxi dx j ,
rij
Kij = ∫ χi* (xi ) χ *j (xi )
Одноэлектронный интеграл
h ii
(69)
1
χi (x j ) χ j (x j )dxi dx j .
rij
описывает энергию электрона на орбитали
ных электронов. Двухэлектронный кулоновский интеграл
J ij
(70)
чi
в поле ядра без осталь-
описывает энергию межэлектронного отталкива-
ния при независимом движении электронов. Двухэлектронный обменный интеграл
энергии взаимодействия электронов с параллельными спинами на орбиталях
K ij
отражает понижение
χi и χ j .
Полная энергия атома с замкнутыми электронными оболочками (по 2 электрона на каждой орбитали) в
методе ХФ вычисляется следующим образом:
N
N
N
E = 2∑ εi − ∑∑ (2J ij − K ij ) .
1=1
i =1 j=1
(71)
Подчеркнем, что оператор Фока (65) сам зависит от полного набора одноэлектронных волновых функций
и его решение ищется самосогласованно. По этой причине метод ХФ иногда отождествляют с методом ССП.
Поскольку, однако, общая стратегия самосогласования проявляется в квантовой химии во многих контекстах,
название "метод Хартри-Фока" более точное.
Наличие обменного члена в операторе Фока эквивалентно учету корреляции в движении электронов с
одинаковыми спинами, связанных с разными орбиталями (обменной корреляции). Кулоновская корреляция,
вызванная взаимным отталкиванием электронов независимо от их спинов, в методе ХФ не учитывается: это является следствием приближения независимых частиц и представляет собой существенный недостаток метода.
Уравнения ХФ могут в принципе быть решены численно любым стандартным методом решения интегрально-дифференциальных уравнений.
10. Ограниченный и неограниченный методы Хартри-Фока
Многоэлектронная волновая функция N-электронной системы должна быть антисимметричной и отвечать определенным проекциям N-электрон-ных орбитального углового и спинового моментов. Она также будет
собственной функцией оператора квадрата полного спина системы S2, если она построена из пространственных
спин-орбиталей, занятых парой электронов с противоположными спинами. Когда волновая функция аппроксимируется единственным детерминантом Слейтера, состоящим из таких спин-орбиталей, метод называется ограниченным (по спину) методом Хартри-Фока (ОХФ). Если же требование быть собственной функцией S2 на
волновую функцию не накладывается, метод называется неограниченным методом ХФ (НХФ), и электроны с
противоположными спинами занимают разные АО (рис. 9). В последнем случае достигается большая гибкость
волновой функции. НХФ, как правило, применяется для систем с открытыми оболочками, обеспечивая для них
более низкое значение энергии.
Энергия
23
Ψ ( Li ) = ϕ (1s)αϕ (1s)βϕ (2s)α
Рис. 9. Волновая функция атома Li в ограниченном и неограниченном методе Хартри-Фока
Таблица 5
Собственные значения квадрата спинового момента электронов
для разных спиновых состояний
Спиновое
состояние
Собственное зна-
Спин
2
Спиновое
чение S
состояние
Спин
Собственное значение S2
Синглет
0
0
Квартет
1.5
3.75
Дуплет
0.5
0.75
Квинтет
2
6.0
Триплет
1
2.0
"N-тет"
N −1
2
s(s+1)
Например, энергия дублетного спинового состояния молекулы CH3, рассчитанная ОХФ равна –39.555
a.е. при
S2 = 0.750, тогда как расчет НХФ дает –39.559 a.е. при S2 = 0.761 (табл. 5). Говорят, что волновая
функция НХФ не соответствует чистому спиновому состоянию. Степень чистоты спинового состояния оценивается по величине квадрата полного спина системы S2, которая должно быть равна s(s+1).
Из-за снятия фундаментального требования быть собственной функцией оператора S2 в больших системах иногда получают НХФ-решения с более низкой энергией и с собственными значениями S2, завышенными
более чем на 10% по сравнению с истинными. Решения, у которых симметрия по спину понижена, называются
"НХФ-нестабильными". Такая нестабильность указывает на непригодность описания волновой функции единственным детерминантом. Заметим, однако, что свойства ОХФ-решения часто точны даже в присутствии большой НХФ-нестабильности.
11. Метод Кона-Шэма
Многоэлектронная волновая функция Ψ({x}) связана с электронной плотностью (ЭП) основного состояния
ρ (r ) соотношением:
ρ(r=r1) = |e|
∫ ∫ ..∫ ∫ ∫
r2 r3
s1 s 2 s3
Ψ*({x})Ψ({x})dr2dr3…ds1ds2ds3…
(72)
24
Интегрирование в (72) проводится по пространственным координатам всех электронов, кроме одного, и по
спинам всех электронов. Поскольку электроны неразличимы, ρ(r) есть средняя по системе плотность электронов.
Электронная энергия системы зависит от ρ(r) следующим образом:
E[ρ ] = ∫ VЯД (r )ρ (r)dV +
1 ρ (r)ρ (r')
drdr ′ + G[ρ ] ,
2 ∫ ∫ r − r'
(73)
где Vяд (r ) – потенциал ядра, а G[ρ] – функционал ЭП, представляющий собой сумму кинетической энергии и
неклассической энергии электрон-электронного взаимодействия Exc[ρ], связанной с обменом и корреляцией
электронов. Теорема Хоэнберга-Кона утверждает, что G[ρ] – одинаковый для всех многоэлектронных систем
(универсальный) функционал электронной плотности, причем, точная ЭП основного состояния обеспечивает
минимум функционала (73). Теория, которая изучает способы расчета электронной структуры молекул и кристаллов, основываясь на минимизации функционала (73), называется теорией функционала плотности (ТФП).
Предположим, что ЭП основного состояния взаимодействующих и невзаимодействующих электронов
одинаковы. Выражение для кинетической энергии электронов (7) одинаково для всех систем. Тогда минимизация функционала (73) относительно одноэлектронных функций χi(r) при дополнительном условии их ортонормировки и постоянства числа электронов в системе дает одноэлектронные уравнения Кона-Шэма:
⎡ 1 2
⎤
ρ (r′)
dr′ + υˆxc (r) ⎥ χi (r) = ε i χi (r) ,
⎢ − ∇ + V(r) + ∫
r − r′
⎣ 2
⎦
где ρ =
∑
χi 2,
V( r )
– электростатический потенциал, создаваемый ядром,
i=1,2,3,…
υˆxc (r ) =
i
(74)
δ E xc [ρ ]
δρ
– обменно-
корреляционный потенциал, Exc[ρ] – обменно-корреляционная энергия, зависящая от электронной плотности.
Система уравнений (74) решается методом самосогласованного поля и заменяет в теории функционала
плотности одноэлектронные уравнения Хартри-Фока. Все приближения связаны со способами аппроксимации
обменно–корреляционного потенциала
υˆxc ( r ) :
поиск таких приближений составляет важный раздел теории
функционала плотности.
Энергии одноэлектронных состояний в теории функционала плотности равны
ε i (ТФП ) =
δ E[ρ ]
,
δ ni
(75)
(ni – электронная заселенность состояния i), тогда как в методе Хартри-Фока они равны разности энергий двух
состояний с числом электронов, отличающимся на единицу:
εi(ХФ) = EХФ, n = 1 – EХФ, ni
i
= 0.
(76)
Таким образом, решения одноэлектронных уравнений Кона-Шэма, вообще говоря, отличны от решений, полученных методом Хартри-Фока. Это иногда дает преимущество теории функционала плотности. Например, разложив энергию (76) в ряд Тейлора и приняв во внимание соотношение (75), можно получить для потенциала ионизации достаточно точное выражение
25
1
Ii ≈ – εi(ni = ).
2
(77)
Можно сказать, что потенциал ионизации может быть вычислен с помощью некоего переходного состояния с наполовину уменьшенной электронной заселенностью высшей занятой электронами орбитали.
12. Электронная структура атомов с точки зрения квантовой химии
Однодетерминантная волновая функция многоэлектронного атома, отвечающая определенным проекциям N-электронного орбитального L и спинового S моментов, не обязательно является собственной функцией
операторов квадрата полного орбитального L2 и спинового S2 моментов атома. Тогда собственные функции
этих операторов являются линейными комбинациями детерминантов Слейтера, отвечающими одним и тем же
N
значениям квантовых чисел полных орбитального
L = ∑ li
S=
и спинового
N
∑ Si
моментов в пределах
i
i
некоторой электронной конфигурации. Под электронной конфигурацией атома понимают определенное распределение электронов по оболочкам:
(nl )1k , (nl )k2 , (nl )3k (nl )kb
1
2
3
....,
b
. Каждая (nl)i оболочка представляет
b
собой набор 2(2li+1) спин-орбиталей, из которых kj спин-орбиталей
(∑ k j = N )
заняты электронами и
j
включены в детерминант Слейтера. Эти kj спин-орбитали можно выбрать для каждой (nl)i
оболочки
k
С 2 (j2li +1) способами, как это следует из правил комбинаторики. Например, для электронной конфигурации
атома углерода (1s)2 (2s)2 (2p)2 можно построить 15 детерминантов, а из них составить 15 линейных комбинаций, соответствующих определенным значениям квантовых чисел L и S.
Таким образом, каждой электронной конфигурации соответствует некоторое число однодетерминантных
функций, определяемое числом незамкнутых оболочек. Совокупность этих функций, характеризующаяся одними и теми же значениями квантовых чисел орбитального L и спинового S моментов, называется термом. Отдельные волновые функции терма отличаются квантовыми числами проекций этих моментов Lz и Sz. Если пренебречь спин-орбитальным взаимодействием, то все волновые функции терма отвечают одному и тому же
(2L+1)(2S+1) кратно вырожденному энергетическому уровню атома. Спин-орбитальное взаимодействие расщепляет этот вырожденный уровень на так называемые уровни тонкой структуры. Энергия терма равна средневзвешенному значению энергий уровней тонкой структуры. Именно минимум этой энергии определяет порядок, в котором электроны заполняют атомные оболочки. Отсюда следует, что понятие электронной конфигурации атома не является строгим физическим понятием, а зависит от приближения, в котором рассматривается
атомная электронная структура.
В то же время, оболочечная структура атомов является объективным фактом, определяющим в конечном счете периодические свойства элементов. Полная ЭП оболочечную структуру явно не показывает и важно
иметь способ надежно выделить электронные оболочки атомов. Рассмотрим этот вопрос более подробно.
Электронная плотность атома является скалярной функцией координаты r. Как известно, имеет ли
скалярная функция максимум или минимум в точке экстремума, определяется знаком ее второй производной
или кривизной в этой точке. Кривизна одномерной функции f(x) отрицательна в точке x, если значение f(x)
больше среднего значения этой функции в соседних точках x+dx и x–dx, и положительна в противополож-
26
ном случае.
ЭП атома имеет максимум в положении ядра и спадает по мере удаления от ядра. Рассмотрим одномерную функцию f(x), моделирующую распределение электронов в атоме вдоль радиального направления, а
также ее первую и вторую производные (рис. 10). Угол наклона f(x) в точке x1–Δx больше, чем в точке x1+Δx,
вторая производная по х в точке x1 отрицательна, а это значит, что кривизна f(x) в точке x1 положительна. Наклон f(x) в точке x2+Δx больше, чем в точке x2–Δx, т.е. кривизна f(x) здесь отрицательна. В области от x1 до x2
функция f(x) имеет точку перегиба, в которой кривизна равна нулю.
Говорят, что f(x) концентрируется в областях, в которых
d 2 f (x)
d 2 f (x)
<0 и разрежается при
>0.
dx 2
dx 2
Понятно, что f(x) максимально концентрируется в точке, где отрицательная кривизна максимальна. Наличие
областей локальной концентрации или разрежения f(x) не связано с наличием максимумов или минимумов самой функции. Например, функция f(x), изображенная на рис. 10, имеет максимум при х=0 и других экстремумов не имеет (т.е. отсутствуют точки, в которых f '(x)=0), а области, где
d 2 f (x)
d 2 f (x)
<0
и
>0 имеют место.
dx 2
dx 2
Все сказанное справедливо и в случае трех измерений, в частности, справедливо для функции ρ(r). В
этом случае необходимо рассматривать вторые производные по трем координатам – лапласиан ЭП:
∇ 2 ρ (r ) =
∂2ρ ∂2ρ ∂2ρ
+
+
∂x 2 ∂y 2 ∂z 2
(78)
В направлении от ядра ЭП экспоненциально уменьшается, и в общем случае кривизна ρ(r) вдоль радиального направления от ядра положительна, как и f(x) в точке x1 на рис. 10. Однако две компоненты кривизны,
перпендикулярные радиальной линии, отрицательны. Таким образом, функция ∇2ρ(r) будет принимать разные
знаки по мере удаления от ядра, явно проявляя электронные оболочки атома.
Для каждой квантовой оболочки наблюдается пара областей (одна отрицательная и одна положительная) с внутренней областью, представляющей собой максимум концентрации электронов. Поскольку ЭП концентрируется при ∇2ρ<0, локальный максимум −∇2ρ будет соответствовать максимуму концентрации ЭП, а локальный минимум −∇2ρ указывает на локальное разрежение ЭП.
Распределение лапласиана ЭП в некоторых атомах приведено на рис. 11. Из рисунка видно, что ∇2ρ −
действительно удобная функция для анализа строения электронных оболочек атомов.
27
Рис. 10. График монотонно убывающей функции f(x) и ее первой и второй производной
28
-∇2ρ
Рис. 11. Лапласиан ЭП атомов водорода (а), неона (б), лития (в), бериллия (г), бора (д), углерода (е), азота (ж) и
кислорода (з). Расчет UHF/6-31G(d))
Важно, что лапласиан ЭП определяет концентрацию электронов на валентной оболочке. Эта квантовая оболочка атома также разделена на внутреннюю область, в которой ∇2ρ<0 и внешнюю, в которой ∇2ρ>0. Часть
оболочки, внутри которой ∇2ρ<0, носит название области концентрации заряда валентной оболочки (КЗВО).
При образовании химической связи валентная оболочка искажается, приводя к появлению максимумов, соответствующих числу и относительному расположению электронных пар, отвечающих предсказаниям модели
Льюиса и модели отталкивания электронных пар валентной оболочки Гиллеспи.
Оболочечная структура атома проявляется в радиальной функции распределения (52) по иному. Максимум функции радиального распределения определяет значение r, в котором с наибольшей вероятностью
можно обнаружить ЭП, усредненную по бесконечно малому объему, лежащему между двумя концентрическими сферами радиуса r и r+dr. Однако фактическое распределение ЭП в трехмерном пространстве максимумом в
этой области не обладает.
13. Квантово-химическая трактовка решений одноэлектронных уравнений
Хартри-фоковские энергии орбиталей εi имеют вполне определенный физический смысл. Если удалить с
орбитали χi один электрон (ионизировать атом), изменение энергии системы можно приближенно записать как
N
⎡
1 ⎞⎤
⎛
ΔE ≈ −ε i = − ⎢ h i + ∑ ⎜ J ij − K ij ⎟ ⎥ ,
2 ⎠⎦
j ⎝
⎣
(79)
считая, что после удаления электрона система не перестроится (останется “замороженной”). Этот результат называется теоремой Купманса. Мы видим, что орбитальные ХФ-энергии дают оценку потенциалов ионизации I
29
− энергий, которые необходимо сообщить системе, чтобы удалить какой-либо из ее электронов. По потенциалу ионизации можно судить о прочности связи электрона данной орбитали с атомным остовом.
Первый потенциал ионизации I1 описывает энергию отрыва электрона с высшей занятой атомной орбитали. При этом предполагается, что и исходный атом, и образовавшийся ион находятся в основных (невозбужденных) состояниях. Потенциалы ионизации I2, I3, I4 и т.д. отвечают дальнейшим последовательным отрывам
электронов от ионов. Для атома с N электронами I1 < I2 < I3 <I4 …< IN. Например, для атома углерода первые
потенциалы ионизации с верхней занятой 2p-АО и глубинной 2s-АО составляют соответственно 11.26 и 20.0
эВ. Зависимость потенциала ионизации от порядкового номера элемента имеет ярко выраженный периодический характер, что иллюстрирует рис. 12.
I1, эВ
I1 = E(N) – E(N–1)
Рис. 12. Зависимость потенциалов ионизации элементов от атомного номера
В пределах периода, как правило, с увеличением атомного номера потенциалы ионизации возрастают.
Исключения связаны с устойчивостью замкнутых оболочек (IО < IN; IMg > IAl) с максимальной мультиплетностью. Для экспериментального определения I применяют фотоэлектронную и рентгеноэлектронную спектроскопию.
Более точно потенциал ионизации следует вычислять с помощью выражения (76). Тогда первый потенциал ионизации, вычисляется как разность ХФ энергий атома и образующегося иона:
I1 = E(N) – E(N–1)
(80)
Сродством к электрону Ax называют энергию, которая высвобождается при присоединении к нейтральному атому одного электрона. Наибольшими значениями Ax обладают атомы галогенов. По аналогии со
сказанным выше,
A1 = E(N+1) – E(N),
(81)
т.е. сродство к электрону можно приближенно охарактеризовать энергией низшей свободной (виртуальной)
АО.
30
Приравнивая нулю дифференциал функционала электронной энергии (73) при постоянном числе
электронов и при неизменном ядерном потенциале
d{E[r] - μN[r(r)]}Vяд = const
= 0,
(82)
(N = ∫ ρ(r) dv), можно определить электронный химический потенциал μ:
⎛ ∂E ⎞
μ=⎜
⎟
⎝ ∂ N ⎠ Vяд = const
(83)
При образовании химической системы из атомов химический потенциал выравнивается, при этом происходит переток электронов к атому с большим значением μ. По смыслу это совпадает с эмпирически введенной
Полингом электроотрицательностью.
На рис. 13 представлена зависимость энергии атома от числа электронов. Как видно из рисунка, μ характеризует наклон кривой E = E(N). Представим производную
∂E
∂N
в конечно-разностном виде (предполагая Е
гладкой функцией N):
μ≈
E ( N + 1) − E ( N − 1)
( N + 1) − ( N − 1)
=
[ E ( N + 1) − E ( N )] + [ E ( N ) − E ( N − 1)]
2
(84)
Отсюда следует
μ≈−
I+A
2
,
(85)
что с точностью до знака совпадает с определением электроотрицательности по Малликену (I1 – первый потенциал ионизации, А – сродство к электрону).
Рис. 13. Зависимость энергии атома от числа электронов
Определенная таким образом величина
χ=
I+A
2
(86)
называется абсолютной электроотрицательностью (Парр и др., 1978).
Пирсон (1983), действуя аналогичным образом, показал, что скорость изменения химического потенциала μ при изменении числа электронов N есть абсолютная химическая жесткость
1 ⎛ ∂μ ⎞
η= ⎜
⎟
2 ⎝ ∂N ⎠V
31
I−A
1⎛ ∂ E ⎞
= ⎜ 2⎟
≈
2
= const 2 ⎝ ∂N ⎠V = const
яд
яд
2
(87)
В рамках метода Хартри-Фока атомная жесткость есть
η=
1
(ЕНСАО – ЕВЗАО),
2
(88)
где ЕНСАО – энергия нижней свободной атомной орбитали; ЕВЗАО – энергия высшей занятой атомной орбитали.
Можно также определить индекс электрофильности (Парр и др., 1999)
μ2 χ2
ω=
=
,
2η 2η
(89)
величина которого может использоваться как индекс склонности атома (молекулы) к атаке электрофила. Другой важный индекс – поляризуемость
α , которая представляет собой меру линейного отклика электронной
плотности в присутствии бесконечно малого электрического поля F и представляет собой вторую производную
от энергии:
⎛ ∂2 E ⎞
⎟ ,
⎝ ∂Fa ∂Fb ⎠
α ab = − ⎜
a, b = x, y, z
(90)
Экспериментально наблюдаемое значение поляризуемости является усредненным по координатам:
1
α = (α xx + α yy + α zz )
3
Величины
(91)
χ , η , ω и α являются характеристиками не только для атомов, ионов, молекул, но и для дру-
гих многоэлектронных многоядерных систем. Они лежат в основе широко применяемых в химии концепций
принципа жестких и мягких кислот и оснований, принципа максимальной жесткости Пирсона, и принципа минимальной поляризуемости (Чаттарадж и Поддар, 1998). Все эти характеристики для атомов проявляют явную
периодичность (рис. 14).
Атомы с закрытыми оболочками или подоболочками имеют большую Наиболее электроотрицательные
элементы обладают максимальной электрофильностью ω, и наоборот. Для легких атомов главных подгрупп отношение χ / η приблизительно постоянно и для них электроотрицательность
χ может рассматриваться как
хорошая мера относительной электрофильности в пределах группы. Отметим, что ωF > ωCl и
αF <
χ F > χ Cl, однако
α Cl.
В общем случае поляризуемость
α
уменьшается по периоду и увеличивается в группе с ростом Z. Для η
и ω наблюдается обратная тенденция. В зависимостях η и ω от Z, структура электронных подоболочек проявляется более отчетливо.
Перечисленные характеристики позволяют проводить оценку технологических свойств материалов. Например, как показывает рис. 15, существует явная корреляция между величиной электроотрицательности элемента по Малликену и температурой его перехода в сверхпроводящее состояние Тсп. Наивысшими температурами Тсп обладают металлы с
χ ≈ 3.9 эВ: Nb, Tc, Pb.
32
Рис. 14. Зависимости индексов жесткости η, электрофильности ω и поляризуемости α от порядкового номера
элемента
33
Установлено, что максимальные температуры перехода в сверхпроводящее состояние Тсп бинарных
сплавов наблюдаются для систем со средней электроотрицательностью
χ ≈ 4 эВ.
TСП, К
Рис. 15. Корреляция между температурой перехода в сверхпроводящее состояние и их электроотрицательностью для 36 элементов
Таким образом, по свойствам элементов можно отобрать вещества, потенциально перспективные для поиска материалов с высокими температурами перехода в сверхпроводящее состояние.
Другой важной характеристикой веществ является энергия когезии - энергия, необходимая для разрушения всех межатомных связей. Для чистых металлов в кристаллической фазе эта величина связана с электронной
плотностью свободного атома ρа следующим приближенным соотношением (Тал, 1996):
5
E когез
2
2 ∞
⎤3
− ⎡
3
= (3π 2 ) 3 V0 3 ⎢ ∫ ρ a (r )d 3 r ⎥ ,
10
⎢⎣R 0
⎥⎦
(92)
1
V0 - объем элементарной ячейки кристалла, R0 = ⎛⎜ 3N ⎞⎟ 3 , N – число электронов в атоме, Ω – атомный объем.
⎝ 4πΩ ⎠
Рассчитанные по формуле (92) величины энергий когезии прекрасно согласуются с экспериментальными данными для широкого круга металлов (рис. 16). Аналогичная зависимость установлена также между электронной плотностью свободного атома ρа и модулями сдвига металлов.
Рис. 16. Энергия когезии чистых металлов в кристаллической фазе: ∗ – эксперимент, ο – расчет
Download