Строение молекул и квантовая химия

advertisement
Строение молекул и квантовая химия
Версия от 16.01.2008
Москва 2007
Содержание
Глава 1. Основные положения квантовой механики и решение модельных задач.
§1. Основные положения квантовой механики.
§2. Примеры решения стационарного уравнения Шредингера.
§3. Гармонический осциллятор.
§4. Решение уравнения Шрёдингера для водородоподобного атома.
§5. Спин.
Глава 2. Разделение электронного и ядерного движений
§ 6. Адиабатическое приближение и приближение Борна-Оппенгеймера.
Глава 3. Решение электронной задачи в одноэлектронном приближении.
§ 7. Многоэлектронные системы.
§8. Метод Хартри-Фока.
§9. Приближение МО ЛКАО.
§10. Ион H2+ и молекула H2 в приближении МО ЛКАО.
Глава 4. Использование симметрии в квантовой механике
§11. Теория групп и квантовая механика молекул.
§12. π-электронное приближение. Метод Хюккеля. Принципы построения полуэмпирических
вариантов метода МО ЛКАО.
§13. Заряды на атомах, порядок связи, свободная валентность.
Глава 5. Корреляционные эффекты и методы их учета
§14. Корреляционная энергия. Возбужденные состояния молекул.
§15. Метод конфигурационного взаимодействия.
§16. Метод валентных связей или валентных схем.
§17. Многочастичная теория возмущений. Теория возмущений Мёллера-Плессе.
Глава 6. Решение задачи о движении ядер в молекулах
§18. Вращательные и колебательные состояния двухатомных молекул.
§19. Вращательные состояния многоатомных молекул.
§20. Колебательные состояния многоатомных молекул.
Приложение
П1. Атомная система единиц.
П2. Необходимые сведения из теории групп.
П3. Таблицы характеров некоторых точечных групп.
П4.1. Проекторы на неприводимые представления. Построение симметризованных МО.
П4.2. Прямое произведение представлений группы.
П4.3. Классификация электронных состояний атомов и молекул.
П4.4. Правила отбора для колебательных переходов многоатомных молекул.
П5. Построение спиновых функций.
П6. Математическое дополнение.
Литература.
3
7
12
14
22
25
29
31
36
40
47
48
54
57
59
62
65
69
73
78
83
85
93
99
102
103
106
112
116
121
Г л а в а 1. Основные положения квантовой механики и решение модельных задач.
§ 1. Основные положения квантовой механики.
В квантовой механике постулируется, что любое состояние системы, состоящей из
взаимодействующих частиц, полностью определяется комплексной волновой функцией Ψ(q,t), где q –
совокупность пространственных и спиновых координат каждой из частиц, t – время. То есть все
предсказания, которые могут быть сделаны относительно различных свойств системы в данный момент
времени t, определяются значением волновой функции Ψ в этот момент времени. Так, частица с
волновой функцией Ψ(q,t) может быть найдена в любой точке области изменения переменных с
2
плотностью вероятности, определяемой квадратом модуля Ψ (q, t ) = Ψ * (q, t ) ⋅ Ψ (q, t ) .
На волновую функцию накладываются условия однозначности, конечности, непрерывности и
квадратичной интегрируемости во всем пространстве переменных. Выполнение этих условий означает,
что существует интеграл вида ∫ Ψ * (q, t )Ψ (q, t )dq и, более того, этот интеграл имеет конечную величину.
Поясним сказанное на одномерном случае, то есть будем считать, что волновая функция частицы
зависит лишь от одной переменной x. Если область изменения координаты x ( −∞ < x < ∞ или
A < x < B ) разбить на маленькие отрезки Δxi и считать, что в пределах данного отрезка Ψ(x;t)
является постоянной величиной, то вероятность обнаружить частицу в пределах отрезка Δxi будет
определяться, как Ψ * ( x; t ) Ψ ( x; t ) Δxi . Полная вероятность обнаружить частицу где-либо, должна
соответствовать достоверному событию, т.е. быть равной 1 и будет задаваться суммой:
∞
∑ Ψ ( x ; t ) Δx
i
i =1
i
= 1.
(1.1)
Если мы теперь будем уменьшать отрезки Δxi
( Δxi → 0 ) , то, при наложении указанных выше
условий на волновую функцию, получим:
∞
lim ∑ Ψ (xi , t ) Δxi =
Δxi →0
i =1
2
+∞
+∞
*
∫ Ψ (xi , t ) dx = ∫ Ψ (xi , t )Ψ (xi , t )dx =1 .
−∞
2
(1.2)
−∞
Таким образом, приходим к условию, что волновая функция, описывающая одномерное
движение только одной частицы должна быть нормирована на 1. В общем случае волновая функция,
описывающая N частиц, зависит от координат каждой частицы xj, yj, zj и условие ее нормировки должно
быть записано как:
*
(1.3)
∫∫∫ Ψ ( x1 , y1 , z1 ,..., xN , yN , zN , t ) Ψ ( x1 , y1 , z1 ,..., xN , yN , zN , t ) dτ1...dτ N = 1 ,
где dτ 1 = dx1dy1dz1 , dτ 2 = dx2 dy2 dz2 ... dτ N = dxN dy N dz N .
В принципе, волновая функция может быть определена для любой совокупности частиц,
например, для Вселенной или солнечной системы, однако делать это достаточно бессмысленно, так как
если части большой системы практически не взаимодействуют между собой, то полная волновая
функция может быть представлена как произведение волновых функций, описывающих те ее части, где
частицы сильно взаимодействуют друг с другом. Обычно волновые функции используют для описания
микроскопических систем, состоящих из конечного числа сильно взаимодействующих между собой
атомных ядер и электронов, т.е. для атомов, молекул, атомных или молекулярных кластеров и их
ионов. Далее, учитывая взаимодействие микроскопических фрагментов между собой, переходят к
описанию всей макросистемы (вещества).
Рассмотрение системы с учетом ее квантовых свойств необходимо в тех случаях, когда
линейные размеры системы сопоставимы с длиной волны де Бройля. Действительно, если рассмотреть
частицу с массой в 1 грамм (m=1 г), двигающуюся со скоростью v = 1·10-3 м/c, то длина волны де
Бройля λ = h / mv ≈ 7 ⋅ 10 −23 м , будет во много раз меньше размеров частицы. Это означает, что
волновые свойства макротелами проявляться не будут и для их описания можно использовать
классическую механику.
Уравнение, решением которого является волновая функция системы частиц, не может быть
получено путем логических рассуждений. Как и уравнение Ньютона в классической механике, оно
постулируется, и единственным доказательством его правильности служит сравнение теоретических
предсказаний с экспериментальными результатами.
Основное уравнение квантовой механики или уравнение Шрёдингера было постулировано
Эрвином Шрёдингером (E. Schrödinger) в 1926 году:
3
∂Ψ (q, t )
,
∂t
HΨ (q, t ) = i
(1.4)
где q – координаты частицы, t – время, i – мнимая единица (i2=-1),
=
h
(h –постоянная Планка), H
2π
– оператор Гамильтона или гамильтониан.
В общем случае под оператором подразумевается некоторая операция A, которая позволяет
сопоставить каждой функции f в некотором функциональном пространстве другую вполне
определенную функцию g, т.е. Af=g. Например, можно сказать, что производная волновой функции по
времени ∂Ψ/∂t есть результат действия оператора ∂/∂t на волновую функцию Ψ. Соответственно,
действие операторов V(x), являющихся функциями координаты x, сводится к умножению функции на
V(x). В квантовой механике также полагают, что все рассматриваемые операторы A являются
линейными, т.е. выполняются условия: A(cΨ ) = cAΨ и A(Ψ1 + Ψ2 ) = AΨ1 + AΨ2 , где с – произвольная
постоянная, и их собственные функции образуют полные наборы, а следовательно, для любой
волновой функции существует представление в виде разложения по собственным функциям оператора
A:
∞
Ψ = ∑ cn Ψ n .
(1.5)
n
Каждой измеряемой величине (энергии, импульсу, координате и т.д.) сопоставляется свой
оператор A и значение измеряемой на опыте величины определяется как среднее значение этой
величины < A >, определяемой через интеграл вида:
A = ∫∫∫ Ψ * ( x, y , z , t )A(x, y , z , t )Ψ (x, y, z , t )dxdydz ,
(1.6)
Ω
где интегрирование ведется по всей области изменения координат Ω, x, y, z ∈ Ω , которой
сопоставляется оператор A(x,y,z,t). Например, если мы хотим вычислить среднее значение координаты
частицы <x>, описываемой функцией Ψ(x,y,z;t) по всему объему, то необходимо вычислить интеграл
x =
+∞ +∞ +∞
∫ ∫ ∫ Ψ (x, y, z, t )xΨ (x, y, z, t )dxdydz ,
*
(1.7)
− ∞ − ∞− ∞
Оператор H , входящий в уравнение Шрёдингера (1.4), соответствует полной энергии системы E,
которая складывается из кинетической и потенциальной энергии E = T + V, и, соответственно,
гамильтониан есть сумма операторов кинетической и потенциальной энергии: H = T + V . Эти
операторы строятся по аналогии с классической механикой, где кинетическая энергия частицы
mv 2
p2
=
, где p = mv есть импульс частицы с массой m, двигающейся со скоростью
определяется как
2
2m
v.
Оператор импульса p определяется следующим образом:
p = px i + p y j + pz k ,
где i , j , k - единичные векторы (орты) декартовой системы координат, а px = −i
pz = −i
(
∂
.
∂z
)
p = pi p = −
2
Отсюда
2
следует
выражение
для
оператора
∂
∂
; p y = −i
;
∂x
∂y
квадрата
импульса
⎛ ∂2
∂2
∂2 ⎞
2
⎜ 2 + 2 + 2 ⎟ = − Δ. Если мы имеем N частиц, то оператор кинетической
⎝ ∂x ∂y ∂z ⎠
N
энергии будет иметь вид суммы T = ∑
i =1
определяться законом Кулона V =
pi2
. Для частиц с зарядами qi потенциал взаимодействия будет
2mi
1
N
qi q j
∑ r −r
4πε ε
0
i< j
i
, где суммирование идет однократно по всем парам
j
взаимодействующих частиц, то есть взаимодействие между i-ой и j-ой частицами учитывается лишь
один раз. Под ε0 и ε подразумеваются абсолютная и относительная диэлектрические проницаемости
4
вакуума,
ri − rj =
последняя
далее
будет
опущена,
(x − x ) +( y − y ) +(z − z )
2
i
2
j
i
j
i
2
j
поскольку
для
вакуума
равна
единице;
- расстояние между частицами i и j.
Таким образом, оператор гамильтона H для системы из N заряженных частиц, на которые не
действуют внешние силы, можно записать как:
N
1 N qi q j
p2
.
(1.8)
H =∑ i +
∑
4πε 0ε i < j ri − rj
i =1 2 mi
Например, для молекулы воды (H2O) в оператор Гамильтона будут входить импульсы трех ядер
(два ядра водорода и одно кислорода) и 10 электронов, а также потенциалы кулоновского
взаимодействия всех пар частиц. Если индексом α обозначить номер ядра, соответственно Mα, pα и
qα=eZα – его массу, импульс и заряд, а индекс i использовать для идентификации электронов с зарядом
e и массой me, то гамильтониан этой молекулы запишется в следующем виде:
10
pα2
p2
+∑ i +
α =1 2 M α
i =1 2 me
3
H ( H 2O ) = ∑
⎧ 3 10
⎫
3
10
Zα Z β
Zα
1 ⎪
⎪
+∑
+∑
⎨−∑∑
⎬
4πε 0ε ⎪ α =1 i =1 Rα − ri α , β Rα − Rβ
i , j ri − rj ⎪
i< j
α <β
⎩
⎭
e2
(1.9)
Первую сумму в этом уравнении можно интерпретировать как оператор кинетической энергии ядер Tn ,
вторую – как оператор кинетической энергии электронов Te , а три последние суммы соответствуют
потенциалам взаимодействия электронов и ядер Vne , ядер друг с другом Vnn ,и
электронов друг с
другом Vee : H = Tn + Te + Vne + Vnn + Vee . Такая структура оператора Гамильтона сохраняется для
любой молекулы из N ядер и M электронов.
Из соображений удобства в квантовой механике атомов и молекул вместо стандартной системы
единиц (СИ) используют атомную систему единиц, в которой
• за единицу длины принято среднее расстояние электрона от ядра в атоме водорода
a0=0,52917706·10-10 м = 0,52917706·10-8 см = 0,52917706 Å;
• за единицу массы - масса покоя электрона me = 9,1093897·10-31 кг;
• за единицу заряда - элементарный заряд или абсолютное значение заряда электрона
е=1,60217733·10-19 Кл;
•
за единицу действия - постоянная Планка
=
h
= 1, 05457267i10−34 Дж·с.
2π
Единицей энергии в атомной системе единиц является удвоенная энергия основного состояния
электрона в атоме водорода, называемая Хартри. 1 Хартри (или 1 а. е. энергии) = 4,35974812·10-18 Дж.
Скорость света в атомной системе единиц равна 137,04, а единица времени составляет 2,53556·10-17 с. В
приложении П1 атомная система единиц обсуждена более подробно, и приведена таблица перевода
атомной единицы энергии в другие наиболее употребительные в химии и физике единицы измерения.
В атомной системе единиц уравнение Шрёдингера для атомов и молекул записывается в более
простом виде, поскольку ћ=1 и 4πεε0=1. Например, выписанный ранее оператор Гамильтона для
молекулы воды имеет вид:
H
−
( H 2O ) =
1
2
3
∑
α
=1
Δα
1
−
Mα
2
10
∑
i =1
Δi −
3
10
∑
∑
α
=1 i =1
Zα
R α − ri
+
3
∑
α β
,
α <β
Zα Z β
Rα − R β
+
10
∑
i, j
i< j
1
ri − r j
.
Отметим, что введенное правило соответствия для построения операторов квантовой механики
по аналогии с классической механикой нужно использовать только для декартовых координат, так как
оно не инвариантно по отношению к замене переменных, например, к замене декартовых координат на
полярные. Следует, однако, иметь в виду, что существует много физических объектов, которые не
имеют классических аналогов, и поэтому построение для них соответствующих гамильтонианов часто
является наиболее трудной задачей. Отсутствие четкого рецепта для нахождения гамильтониана
объекта является одним из главных недостатков современной квантовой механики.
Если оператор потенциальной энергии V(q,t) не зависит от времени, т.е. V(q,t)≡V(q) для любого
t, то решение уравнения Шрёдингера удобно искать в виде произведения Ψ (q, t ) = Φ(q )Ξ(t ) . Подставив
это произведение в уравнение (1.4) и разделив его правую и левую части на Φ(q )Ξ(t ) , получим
5
⎤
1 ⎡ p 2 Φ (q )
1 ∂Ξ(t )
−
+ V (q )Φ (q )⎥ = i
.
(1.9)
⎢
(t ) ∂t
Φ (q ) ⎣
2m
Ξ
⎦
Так как левая часть уравнения зависит только от переменной q, под которой мы подразумеваем
координаты частиц, составляющих систему, а правая часть – только от времени, то обе части должны
быть постоянной величиной (константой разделения), которую обозначим как E. Функция Ξ ( t ) легко
находится из уравнения:
i ∂Ξ ( t )
= E.
Ξ ( t ) ∂t
(1.10)
Действительно, представив Ξ ( t ) = eα t и подставив в уравнение, получим:
i α eα t = Eeα t ,
разделив обе части на eαt, находим, что α =
Ξ ( t ) = e −iEt .
E
iE
= − , то есть
i
Для функции Φ(q) получаем уравнение:
HΦ (q )
= E или HΦ (q ) = EΦ (q ) ,
(1.11)
Φ (q )
которое называют стационарным уравнением Шрёдингера. Постоянная разделения E имеет физический
смысл полной энергии системы, которая не меняется с течением времени t (случай консервативной
системы в классической механике).
Таким образом, в случае не зависящего от времени гамильтониана волновая функция Ψ(q,t)
имеет вид:
iE
Ψ (q, t ) = Φ(q )e ,
(1.12)
где энергия E находится из решения стационарного уравнения (1.11).
Заметим, что плотность вероятности
−
t
⎛ iE iE ⎞
⎜ − ⎟t
⎠
Ψ * (q, t ) ⋅ Ψ (q, t ) = Φ * (q )Φ (q )e ⎝
= Φ * (q )Φ (q ) также не зависит от времени.
Уравнение HΦ (q ) = EΦ (q ) представляет задачу на собственные значения оператора H .
Функции Φn(q), которые удовлетворяют этому уравнению и соответствующие им энергии En называют
собственными функциями и собственными значениями оператора H . Собственные функции
подразделяют на два типа. Первый класс функций описывает локализованные состояния квантовой
системы, т.е. частицы, входящие в состав системы удерживаются силами в некоторой конечной
области пространства. В этом случае решения уравнения Шрёдингера существуют только для
отдельных значений энергии E = E1, E2 … En …, и мы имеем дискретный спектр собственных значений
гамильтониана. Второй тип функций описывает делокализованные частицы, которым соответствует
непрерывный спектр собственных значений. Если для простоты ограничиться случаем, когда
гамильтониан обладает только дискретным спектром: HΦ n (q ) = En Φ n (q ) , (n=1,..,∞), то произвольная
функция Ψ может быть записана в виде разложения (1.5) по собственным функциям Φn(q), а квадраты
модулей коэффициентов разложения c n
2
= c n* ⋅ c n будут задавать вероятность обнаружить у частицы
значение энергии En.
В 1925 г. Вольфгангом Паули в квантовую механику был введен так называемый принцип
исключения (запрета), который применительно к системам электронов может быть сформулирован
следующим образом: волновая функция системы частиц с полуцелым спином (в частности,
электронов) должна быть антисимметрична относительно перестановки координат любых двух
частиц:
Ψ (q1 ,..., qi ,..., q j ,..., qn , t ) = − Ψ (q1 ,..., q j ,..., qi ,..., qn , t ) .
(1.13)
Смысл этого постулата будет пояснен ниже, при рассмотрении многоэлектронных систем. В
настоящем разделе мы ограничимся лишь его формулировкой.
6
§2. Примеры решения стационарного уравнения Шредингера.
В атомной системе единиц стационарное уравнение Шрёдингера (1.11) для частицы массы m,
двигающейся в потенциале U(x,y,z) имеет вид:
∂2
∂2 ⎞
1 ⎛ ∂2
−
+
+
⎜
⎟ Ψ ( x, y , z ) + U ( x, y , z ) Ψ ( x, y , z ) =
2m ⎝ ∂x 2 ∂y 2 ∂z 2 ⎠
или
E Ψ ( x, y , z )
∂2
∂2
1 ⎛ ∂2
+
+
⎜
2m ⎝ ∂x 2 ∂y 2 ∂z 2
⎞
⎟ Ψ ( x, y , z ) +
⎠
( E − U ( x, y, z ) ) Ψ ( x, y, z ) = 0.
(2.1)
Рассмотрим примеры решения уравнения (2.1) в случае постоянного потенциала U(x,y,z)=const или, в
частном случае, U(x,y,z)=0.
1. Область изменения переменных x,y,z не ограничена, что соответствует свободному движению
частицы.
Заметим, что волновую функцию Ψ(x,y,z) можно представить в виде произведения
Ψ ( x, y, z ) = Ψ (x ) ⋅ Ψ ( y ) ⋅ Ψ ( z ) , так что для каждой функции имеется свое уравнение Шрёдингера
(переменные x, y, z разделяются).
1 ∂2
−
Ψ ( x ) = Ex Ψ ( x ) ;
2m ∂x 2
1 ∂2
−
Ψ ( y ) = Ey Ψ ( y ) ;
2m ∂y 2
2.2
1 ∂2
−
Ψ ( z ) = Ez Ψ ( z ) .
2m ∂z 2
Таким образом, для решения задачи достаточно найти решение одного из уравнений (2.2),
например:
−
1 ∂2
Ψ ( x ) = EΨ ( x ) .
2m ∂x 2
(2.3)
Уравнение (2.3) является дифференциальным уравнением второго порядка, и для нас интересны
его ограниченные, непрерывные и дифференцируемые решения в интервале ( −∞, +∞ ) . Будем решение
Ψ(x) искать в виде Ψ ( x ) = eiα x . Заметим, что если такое решение существует, то всякое другое
решение, получаемое умножением его на постоянный коэффициент, обладает теми же свойствами,
поэтому мы не будем различать решения, отличающиеся постоянным множителем. Если допустимы
два решения, то всякая их линейная комбинация является также решением уравнения. В частности,
если Ψ(x) есть собственная функция H , то ее действительная и мнимая части также являются
собственными функциями.
Подставляя eiα x в уравнение (2.3), получаем:
1 2 iα x
α e = Eeiα x ,
2m
(2.4)
откуда следует, что α 2 = 2mE или α = ± 2mE , причем E > 0 и
Ψ ( x ) = ei
2 mE x
= eipx ,
(2.5)
2
p
для свободного движения частицы. Полное решение запишется как
2m
ip y
Ψ ( x, y, z ) = Aeipx x e y eipz z = Aeipr
(2.6)
так как E =
− ip x
где A – постоянный множитель. Функция e − i 2 mE x = e x также удовлетворяет уравнению (2.3) и
соответствует движению частицы в другую сторону по оси x. Следовательно, общее решение
уравнения (2.1) должно иметь вид:
Ψ ( x, y, z ) = Ae ipr + Be − ipr ,
(2.7)
7
где A и B – произвольные постоянные множители.
Таким образом, в случае свободного движения частицы ее энергия Е может принимать любые
положительные значения, а волновая функция имеет вид (2.7). Такую функцию обычно называют
плоской волной.
Заметим, что разделение переменных возможно и в более сложных случаях, если оператор H
можно представить в виде суммы операторов, зависящих от переменных q1,…,qN, то есть
N
H = ∑ H (qi ) .
(2.8)
i =1
Тогда волновая функция представляется произведением функций:
N
Ψ (q1 , q 2 ,..., q N ) = ∏ Ψi (qi ) = Ψ1 (q1 ) ⋅ Ψ2 (q 2 ) ⋅ ... ⋅ ΨN (q N ) ,
i =1
где функции Ψi(qi) определяются из уравнений HΨi (qi ) = Ei Ψi (qi ) , а полная энергия системы будет
суммой энергий Ei: E = E1 + E2 + ...Ei + E j + Ek ... + EN =
N
∑E .
i =1
i
2. Теперь рассмотрим случай, когда потенциал U ( x, y, z ) = 0 только для определенных значений
переменных в интервалах 0<x<a; 0<y<b; 0<z<c, а вне этих интервалов принимает бесконечные
значения.
⎧ 0 , x ∈ [0 , a ], y ∈ [0 , b ], z ∈ [0 , c ]
U (x , y , z ) = ⎨
(2.9)
⎩ ∞ , x ∉ [0 , a ], y ∉ [0 , b ], z ∉ [0 , c ]
Как и в предыдущем случае, волновую функцию представим в виде произведения одномерных
функций Ψ (x, y, z ) = Ψ ( x ) ⋅ Ψ ( y ) ⋅ Ψ ( z ) , что позволяет перейти к рассмотрению одномерного случая:
1 ∂2
1 ∂2
(
)
(
(
)
)
(
)
Ψ
x
+
E
−
U
x
Ψ
x
=
Ψ (x ) − E x Ψ (x ) = 0 ;
2m ∂x 2
2m ∂x 2
1 ∂2
Ψ( y ) − E y Ψ( y ) = 0 ;
2m ∂y 2
1 ∂2
Ψ(z ) − E z Ψ (z ) = 0 .
2m ∂z 2
(2.10)
Так как движение частицы здесь ограничено потенциальной функцией прямоугольной формы, то эту
задачу обычно называют задачей о движении частицы в «потенциальном ящике с бесконечно
высокими стенками». Общее решение первого из уравнений (2.10) запишется как
Ψ ( x ) = Aeiα x + Be − iα x , или
Ψ ( x ) = A sin (α x ) + B cos (α x ) ,
(2.11)
так как eiα x = cos α x + i sin α x (формула Эйлера).
На решение Ψ(x) в этом случае необходимо наложить граничные условия Ψ ( 0 ) = Ψ ( a ) = 0 ,
так как вероятность найти частицу на стенках и вне такого ящика должна быть равна нулю.
Выполнение граничных условий ведет к следующему:
Ψ ( 0 ) = A sin (α 0 ) + B cos (α 0 ) = 0 ,
(2.12)
так как cos(0)=1 и sin(0)=0, то коэффициент B должен быть равен нулю (B=0). Подставляя x = a ,
получаем A sin (α a ) = 0 , что выполняется, если α a = nπ , где n = 1, 2,3, 4... Значение n = 0 нам не
подходит, так как в этом случае получаем нулевое решение. Таким образом, α =
уравнения (2.10) будет иметь вид
Ψ ( x ) = A sin
nπ
x.
a
(2.13)
Коэффициент А найдем из условия нормировки волновой функции на единицу:
8
nπ
, и тогда решение
a
a
a
∫ Ψ ( x ) Ψ ( x ) dx = 1 = A ∫ sin
*
2
0
0
a
a
∫ Ψ (x )Ψ (x )dx = A ∫ sin
*
2
0
2
0
2
nπ
a
xdx = A2 ;
a
2
⎛ nπ ⎞
x ⎟dx = 1 .
⎜
⎝ a ⎠
(2.14)
⎛ α ⎞ 1 − cos(α )
. Произведя соответствующие
⎟=
2
⎝2⎠
Для взятия интеграла используем формулу sin 2 ⎜
замены
a
∫ sin
0
2
в
⎛ nπ
⎜
⎝ a
a
⎞
x ⎟dx =
nπ
⎠
nπ
a ⎡
⎛ nπ ⎞ 1
x⎟ −
⎢∫ d⎜
2nπ ⎣ 0 ⎝ a ⎠ 2
интеграле
nπ
∫ sin
0
2 nπ
∫
0
2
⎛ nπ
⎜
⎝ a
⎞ ⎛ nπ
x ⎟d ⎜
⎠ ⎝ a
⎞
x⎟ =
⎠
⎛ 2nπ ⎞ ⎛ 2nπ ⎞ ⎤
x⎟d ⎜
x ⎟⎥ =
cos ⎜
⎝ a ⎠ ⎝ a ⎠⎦
=
(2.14),
a
2nπ
получим:
1
⎤ a
⎡
⎢⎣nπ − 0 − 2 (0 − 0)⎥⎦ = 2 .
вычисленное значение интеграла в ту же формулу (2.14), найдем, что A =
Подставив
2
, и, следовательно,
a
волновая функция с учетом нормирующего множителя имет вид:
Ψ (x ) =
2 ⎛ nπ ⎞
x ⎟ , (n,k,l=1,2,3,…,∞).
sin ⎜
a ⎝ a ⎠
(2.15)
Подставив функцию (2.15) в уравнение (2.10), получим выражение для энергии частицы Еx(n):
α
E x (n ) =
n 2π 2
=
.
2m 2ma 2
2
(2.16)
В результате мы получили квантованные значения энергии частицы в одномерном
потенциальном ящике, что изображено на рисунке 2.1, где вместе с уровнями энергии приведены
графики волновых функций для трех нижних состояний частицы.
E3
E2
E1
a
0
U=0
U= ∞
U= ∞
Рисунок 2.1. Вид волновых функций для трех низших состояний частицы в ящике.
Аналогичным способом найдем решения дифференциальных уравнений для переменных y и z.
В окончательном виде энергия и волновая функция для частицы в трехмерном ящике имеют вид:
Ψ ( x, y , z ) =
π 2 ⎛ n2
k2
⎜⎜ 2 + 2
2m ⎝ a
b
8
⎛ kπ ⎞
⎛ lπ
⎛ nπ ⎞
x ⎟ ⋅ sin ⎜
y ⎟ ⋅ sin ⎜
sin ⎜
abc ⎝ a ⎠
⎝ b ⎠
⎝ c
E (n, k , l ) = E x (n ) + E y (k ) + E z (l ) =
+
l2
c2
⎞
⎟⎟ и
⎠
⎞
z ⎟ , (n,k,l=1,2,3,…,∞). (2.17)
⎠
Отметим, что рассмотренная простая задача применима для моделирования движения
электрона в ряде химических задач. В частности, движение π-электрона в цепи полиена CH2═CH–
CH═…CH=CH2 может быть описано как его одномерное движение в ящике с бесконечно высокими
9
стенками, ширина которого определяется длиной цепи полиена. Если N – число атомов углерода в цепи
полиена, то a = R ( N − 1) , где R – межъядерное расстояние между соседними ядрами углерода.
1
3
2
4
N
Рисунок 2.2.
Если рассматривать все π-электроны в полиене как независимые, то они будут заселять
энергетические уровни с разными n по 2 на каждом уровне. Таким образом, последний заселенный
уровень будет иметь энергию
N 2π 2
= EN .
2
4 ⋅ 2 R ( N − 1)
N /2 = n
2=n
1=n
Рисунок 2.2. Энергии уровней полиена.
2
⎛N ⎞ 2
⎜ + 1⎟ π
2
⎠
, и энергия возбуждения электрона на
Первый свободный уровень будет иметь энергию ⎝
2 R ( N − 1)
этот уровень составит:
2
⎛N ⎞ 2
⎜ + 1⎟ π
⎛ ⎛ N ⎞2 N 2 ⎞
N 2π 2
π2
2
⎠
E=⎝
−
=
⎜ ⎜ + 1⎟ −
⎟=
2 R ( N − 1) 4 ⋅ 2 R ( N − 1) 2 R ( N − 1) ⎜⎝ ⎝ 2
4 ⎟⎠
⎠
=
2
⎛ N2
N
N 2 ⎞ π ( N + 1)
2
1
.
+
+
−
=
⎜
⎟
2 R ( N − 1) ⎝ 4
2
4 ⎠ 2 R ( N − 1)
π2
lim E =
N →∞
π2
2R
.
(2.18)
3. Движение частицы с массой m по кругу радиуса r.
r2=x2+y2
r
m
Рисунок 2.4. Круговое движение частицы.
В этом случае уравнение Шрёдингера запишется следующим образом:
∂2 ⎞
1 ⎛ ∂2
2
2
−
⎜ 2 + 2 ⎟ Ψ ( x, y ) + U ( x + y ) Ψ ( x , y ) = E Ψ ( x , y ) .
2m ⎝ ∂x ∂y ⎠
(2.19)
Если r 2 = const , то потенциал U тоже является постоянной величиной, и уравнение выглядит
следующим образом:
10
−
∂2 ⎞
1 ⎛ ∂2
⎜ 2 + 2 ⎟ Ψ ( x, y ) = ( E − U ) Ψ ( x, y ) .
2m ⎝ ∂x ∂y ⎠
(2.20)
Перепишем это уравнение в полярных координатах: x = r cos ϕ ; y = r sin ϕ .
∂ ∂r ∂ ∂ϕ
=
+
∂y ∂y ∂r ∂y
∂ ∂r ∂ ∂ϕ
=
+
∂x ∂x ∂r ∂x
∂
∂ cos ϕ
= sin ϕ +
∂ϕ
∂r
r
∂
∂ sin ϕ
= cos ϕ −
∂ϕ
∂r
r
∂
;
∂ϕ
∂
.
∂ϕ
(2.21)
Подстановка этих выражений в уравнение с учетом того, что
∂Ψ
= 0 , так как r = const , приводит к
∂r
уравнению:
−
1 ∂2
Ψ (ϕ ) = ( E − U ) Ψ (ϕ ) .
2 μ r 2 ∂ϕ 2
(2.22)
Для того чтобы функция была однозначной и непрерывной, необходимо потребовать ее периодичности
Ψ (ϕ ) = Ψ (ϕ + 2π ) . Если функцию искать в виде Aeimϕ , то
m2
m2
1 2 imϕ
или E =
+U .
m e = ( E − U ) eimϕ ; E − U =
2I
2I
2I
im ϕ + 2π )
= Aeimϕ eim 2π , откуда следует eim 2π = 1 , что выполняется только если
При этом Aeimϕ = Ae (
m = 0, ±1, ±2...
Таким образом, получаем квантование энергетических вращательных уровней, которое следует из
граничных условий. Постоянная A определяется из нормировки функции на единицу:
2π
2π
∫ Ψ ( ϕ ) Ψ (ϕ ) d ϕ = A ∫ e
*
2
0
e dϕ = A2 2π = 1 ,
− imϕ imϕ
(2.23)
0
откуда получаем A2 =
1
1
и A=
.
2π
2π
Окончательно, волновые функции и энергии уровней записываются как
1 imϕ
e
2π
m2
Em =
+U,
2I
Ψm =
m = 0, ±1, ±2...
где I = μ r 2 ( момент инерции ).
m= +_ 2
E2=4/2I
m= + 1
E1=1/2I
_
m=0
E0=0
Рисунок 2.5. Энергии уровней для кругового движения частицы.
Эту модельную задачу применим для приближенной оценки энергии возбуждения молекулы
бензола электромагнитным излучением. Будем считать, что шесть π-электронов в бензоле двигаются
независимо друг от друга по окружности с радиусом r = 1,4 Å ≈ 3a0.
11
r
Рисунок 2.6. Схематическое изображение молекулы бензола.
Момент инерции электрона, I = me r 2 = 9 , так как в атомных единицах масса электрона есть 1 и r = 3 .
Получаем следующие уровни энергии:
m= +_ 2
E2=4/(2 9) а.е.
m= +_ 1
E1=1/(2 9) а.е.
m=0
E0=0
Рисунок 2.7. Энергии уровней.
Первый уровень невырожден и может быть заселен двумя электронами с противоположными спинами.
Следующий уровень вырожден дважды и может быть заселен четырьмя электронами. Энергия
возбуждения
E будет определяться как E2 − E1 =
4 1
3
3
− = а.е. = ⋅ 27, 2 эВ = 4.5эВ , что по
18 18 18
18
порядку величины согласуется с экспериментальным значением 4,7 эВ.
§ 3. Гармонический осциллятор.
В классической механике гармоническим осциллятором называют частицу, которая может
двигаться по некоторой оси и на которую действует возвращающая сила, пропорциональная смещению
частицы от начала координат. Если q есть координата частицы на оси, p = mq - ее импульс, m – масса
1 2
kq - потенциальная функция, то функция Гамильтона записывается как:
2
p2 1 2
E = H = T +V =
+ kq .
(3.1)
2m 2
и
В квантовой механике гармоническим осциллятором называют систему, для которой оператор
p
1
H=
+ kq 2 и соответственно волновые функции определяются из решения уравнения
2m 2
Шрёдингера:
2
−
∂2
1
Ψ ( q ) − kq 2 Ψ ( q ) = E Ψ ( q ) .
2
2m ∂q
2
2
(3.2)
Решение этого уравнения разбирается практически во всех книгах, посвященных квантовой
механике, поэтому здесь мы ограничимся лишь формальным описанием схемы его решения и сводкой
получаемых результатов. Подробный вывод (3.2) можно найти в [1,2]. На первом этапе исходное
уравнение (3.2) приводится к виду
d 2 Ψ (q )
+ (α − β 2 q 2 )Ψ (q ) = 0
(3.3)
2
dq
Путем замены α =
2mE
2
и β=
mk
. Это уравнение можно упростить и далее, введя новую
переменную ζ = q β :
12
d 2 Ψ (ζ ) ⎛ α
+ ⎜⎜ − ζ
dζ 2
⎝β
⎞
⎟⎟Ψ (ζ ) = 0 .
⎠
Произведя подстановку Ψ (ζ ) = u (ζ ) ⋅ exp(− ζ
дифференциальному уравнению вида:
⎞
d 2 u (ζ )
du (ζ ) ⎛ α
− 2ζ
+ ⎜⎜ − 1⎟⎟u (ζ ) = 0 .
2
dζ
dζ
⎝β
⎠
2
Если принять, что
(3.4)
2
2 ) , находим, что функция u(ζ) должна удовлетворять
(3.5)
α
= 2ν + 1 , то (3.5), в свою очередь, преобразуется к уравнению Эрмита:
β
d 2 u (ζ )
du (ζ )
− 2ζ
+ 2νu (ζ ) = 0 ,
(3.6)
2
dζ
dζ
Его частное решение при целых ν=n≥0 называется полиномом Эрмита n-степени:
dn
n
u (ζ ) = H n (ζ ) = (− 1) exp(ζ 2 ) n (exp(− ζ 2 )) . Наиболее сложные моменты в процедуре решения
dζ
заключаются в обосновании выбора подстановки для преобразования уравнения (3.4) к (3.5) и
доказательстве существования при целых положительных значениях ν решений (3.4),
удовлетворяющих условиям, накладываемым на волновые функции (см. §2). После этого из выражения
α
= 2n + 1 можно вывести формулу для энергии уровней гармонического осциллятора и найти вид их
β
волновых функций с учетом нормировки.
Окончательно, получим:
1.Собственные значения Ev , определяющие энергетические уровни осциллятора, принимают значения
k⎛
1⎞
1⎞
k
⎛
и n=0,1,2,… (3.7)
⎜ n + ⎟ = ω ⎜ n + ⎟ , где ω =
m⎝
2⎠
2⎠
m
⎝
2.Волновые функции гармонического осциллятора являются произведениями полиномов Эрмита и
гауссовых функций:
En =
1
1
2
ζ
β
⎛
−
⎛
⎞ 2
− q
β ⎞⎟ 2
a
2
2
⎜
(
)
(
)
H
e
ζ
,
Ψ
q
=
⋅
H
q
β
⋅
e
. (3.8)
Ψv ( x) = ⎜
(
)
⎟
n
n
v
1
⎜ 2 n n! π ⎟
v
2
⎝ 2 v !π ⎠
⎝
⎠
Для справки, первые пять полиномов Эрмита имеют вид H 0 (ζ ) = 1 , H 1 (ζ ) = 2ζ , H 3 (ζ ) = 8ζ 3 − 12ζ ,
H 4 (ζ ) = 16ζ 4 − 48ζ 2 + 12 .
2
Рис. 3.1. Волновые функции и уровни энергии гармонического осциллятора.
Из рисунка 3.1 видно, что низший, или как часто говорят, основной уровень энергии при n=0 имеет
ω
. Занимающая этот уровень частица по большей части находится вблизи точки равновесия,
2
так как максимум плотности вероятности ρ (q ) = Ψ0* (q )Ψ (q ) достигается при q=0. Следующий полином
энергию
13
(n=1) пропорционален q, поэтому его произведение на гауссову функцию представляет собой волну с
узлом при q=0. Следующая функция (n=1), как и функция основного состояния имеет максимум при
q=0, но также и еще два в областях с более высоким потенциалом. При дальнейшем возбуждении
гармонического осциллятора главные максимумы волновой функции смещаются к точкам поворота
классического гармонического осциллятора, что в какой-то мере соответствует классическому
результату, так как в точках поворота скорость частицы минимальна, а вероятность ее нахождения
максимальна. Таким образом, при больших возбуждениях среднее распределение плотности
вероятности приближается к классическому, тогда как для n=0 наличие максимума ρ(q) при q=0 не
согласуется с классическим результатом.
§ 4. Решение уравнения Шрёдингера для водородоподобного атома.
Под водородоподобным атомом принято понимать простейшую атомную систему, состоящую
из ядра с зарядом +ZNe и массой M и электрона с зарядом –e, массой me. Если использовать атомную
систему единиц, то в декартовых координатах уравнение Шрёдингера запишется следующим образом:
z
me, -e
(x,y,z)
r
M, e
(X,Y,Z)
R
y
x
⎡ 1
ZN ⎤
∂2
∂2
∂2
1
(
)
(
)
,
Δ
=
+
+
,
Ψ
X
,
Y
,
Z
,
x
,
y
,
z
=
E
Ψ
X
,
Y
,
Z
,
x
,
y
,
z
−
Δ
−
Δ
−
⎥
⎢
N
N
e
∂X 2 ∂Y 2 ∂Z 2
2m e
R Ne ⎦
⎣ 2M
Δe =
∂2
∂2
∂2
+ 2 + 2 ,
2
∂x
∂y
∂z
R Ne = ( X − x ) + (Y − y ) + (Z − z ) ,
(4.1)
где X, Y, Z и x, y, z – координаты ядра и электрона, соответственно. Как мы видим, волновая функция
зависит от шести переменных, которые определяют положение ядра и электрона относительно
неподвижной (лабораторной) системы координат. Разумно полагать, что нас интересует относительное
движение электрона и ядра, а не движение атома водорода как целого. По этой причине перейдем к
новым переменным, в которых можно выделить движение атома как целого и относительное движение
электрона и ядра. За «новые» переменные возьмем координаты центра масс атома (XCM , YCM , ZCM):
MX + me x
MY + me y
MZ + me z
X CM =
; YCM =
; Z CM =
;
M + me
M + me
M + me
2
2
MR + me r
M + me
и относительные координаты:
2
или RCM =
x ' = X − x;
y ' = Y − y; z ' = Z − z
или r ' = R − r
(4.2)
(4.3)
Используя общую формулу для производной при замене переменных, запишем
∂X CM ∂
∂
∂x / ∂
=
+
;
∂X
∂X ∂X CM ∂X ∂x /
∂ ∂X CM ∂
∂x / ∂
=
+
.
(4.4)
∂x ∂X CM
∂x ∂x /
∂x
Аналогичные выражения будут иметь место для производных по производным от Y, y и Z, z.
Продифференцировав новые переменные по исходным, получим
14
∂X CM
∂x /
M
=
;
=1;
∂X
∂X
M + me
∂X CM
me
∂x /
=
;
= −1 .
(4.5)
∂x
M + me ∂x
Подставляя эти выражения в предыдущие равенства, получим запись производных по исходным
переменным через производные по новым переменным:
∂
M
∂
∂
=
+ / ;
∂X M + m e ∂X CM ∂x
me
∂
∂
∂
=
− / .
(4.6)
∂x M + me ∂X CM ∂x
Для вторых производных по исходным переменным получаем следующие выражения:
∂2
∂ ∂ ⎛ M
∂
∂ ⎞⎛ M
∂
∂ ⎞
=
= ⎜⎜
+ / ⎟⎟⎜⎜
+ / ⎟⎟ =
2
∂X
∂X ∂X ⎝ M + me ∂X CM ∂x ⎠⎝ M + me ∂X CM ∂x ⎠
M2
M ⎛ ∂
∂2
∂
∂
∂ ⎞
∂2
⎜⎜ /
⎟
=
+
+
+ /
2
2
/2
M + m e ⎝ ∂x ∂X CM ∂x ∂X CM ⎟⎠
(M + me ) ∂X CM ∂x
и
m e2
∂2
∂2
∂2
=
+
−
2
∂x 2 (M + m e )2 ∂X CM
∂x / 2
−
me
M + me
⎛ ∂
∂
∂
∂
⎜⎜ /
+ /
⎝ ∂x ∂X CM ∂x ∂X CM
⎞
⎟⎟ .
⎠
(4.7)
Аналогичные выражения получим и для частных вторых производных
∂2
∂2
∂2 ∂2
,
и
,
.
∂Y 2 ∂y 2
∂Z 2 ∂z 2
Подставляя выражения для «старых» производных через «новые» в оператор кинетической энергии
1
1
T =−
ΔN −
Δ e , найдем его новый вид:
2M
2m e
T =−
⎛ ∂2
1
∂2
∂2
⎜
+
+
2
2
2
2(M + m e ) ⎜⎝ ∂X CM
∂YCM
∂Z CM
−
⎞
⎟−
⎟
⎠
∂2
∂2
1⎛ 1
1 ⎞⎛ ∂ 2
⎜⎜
⎟⎟⎜⎜ / 2 + / 2 + / 2
+
∂x
∂x
2 ⎝ M me ⎠⎝ ∂x
⎞
⎟⎟ .
⎠
(4.8)
Гамильтониан H в «новых» координатах будет иметь следующий вид:
⎛ ∂2
1
∂2
∂2 ⎞
⎟−
⎜
H =−
+
+
2
2
2
⎟
2(M + me ) ⎜⎝ ∂X CM
∂YCM
∂Z CM
⎠
−
∂2
1 ⎛ ∂2
∂2
⎜⎜ / 2 + / 2 + / 2
2 μ ⎝ ∂x
∂x
∂x
⎞
⎟⎟ −
⎠
ZN
x
/2
+ y/2 + z/2
=
1
1 /
Δ CM −
Δ + U (x / , y / , z / ) = H (RCM ) + H (r / ) .
2(M + m e )
2μ
(4.9)
M ⋅ me
Величину μ =
из выражения (4.9) обычно называют приведенной массой.
M + me
После выполненной замены переменных оператор Гамильтона разделяется на два оператора,
=−
( )
зависящих от переменных центра масс (RCM ) , и относительных переменных r ' , и волновая функция
может быть представлена в этих переменных в виде произведения
15
Ψ (RCM , r / ) = Φ (RCM ) ⋅ Ψ (r / ) .
(4.10)
Подставляя волновую функцию в виде (4.10) в уравнение Шрёдингера и разделив его на Ψ (RCM , r / ) ,
получим:
H (R CM )Φ (R CM ) H (r / )Ψ (r / )
+
− E = 0 (4.11)
Φ (R CM )
Ψ (r / )
или
H (R CM )Φ (R CM )
H (r / )Ψ (r / )
−
+E=
=ε ,
(4.12)
Φ (R CM )
Ψ (r / )
где ε – постоянная величина, так как левая и правая части уравнения равны при произвольных
значениях переменных RCM и r / .
Таким образом, мы приходим к двум уравнениям:
H (R CM )Φ (R CM ) = (E − ε )Φ (R CM ) и
H
(r )Ψ (r ) =
/
/
εΨ
(r )
/
(4.13)
Первое уравнение соответствует свободному движению частицы с массой
как H (RCM ) = −
полученным
( M + m) ,
так
1
Δ CM . Подобная задача уже рассматривалась ранее, и, воспользовавшись
2(M + me )
ранее решением, можно показать, что волновая функция имеет вид
Φ (RCM ) = e
, где E − ε = E / ≥ 0 .
Второе уравнение (4.13) соответствует движению частицы с приведенной массой μ в
центральном поле:
⎡ 1 ⎛ ∂2
⎤
ZN
∂2
∂2 ⎞
⎜⎜ / 2 + / 2 + / 2 ⎟⎟ −
⎢−
⎥ Ψ (x / , y / , z / ) =
/2
/2
/2
μ
2
∂
∂
∂
x
x
x
⎢⎣
x + y + z ⎥⎦
⎝
⎠
= εΨ (x / , y / , z / ) .
(4.14)
Так как масса протона примерно в 1840 раз превосходит массу электрона me, то для не слишком точных
расчетов обычно полагают, что μ = me = 1 , и уравнение Шрёдингера можно представить в компактной
форме:
⎡ 1
Z⎤
⎢− Δ − ⎥ Ψ (r ) = εΨ (r ) ,
r ⎦⎥
⎣⎢ 2
± i RCM
2 ( M + me )( E − ε )
где r - радиус-вектор электрона в системе координат, центрированной на ядре и Z=ZN -заряд ядра.
Решение этого уравнения удобно проводить в сферической системе координат, в которой новые
переменные r, θ и ϕ связаны с декартовыми переменными соотношениями:
x = rSinθCosϕ ,
0≤ r ≤∞;
y = rSinθSinϕ ,
0 ≤θ ≤π ;
z = rCosθ ,
0 ≤ ϕ ≤ 2π .
В сферических координатах оператор Лапласа (лапласиан) будет иметь вид:
Δ=
1 ∂ ⎛ 2 ∂ ⎞ 1
⎜r
⎟ + Λ,
r 2 ∂r ⎝ ∂r ⎠ r 2
16
1 ∂ ⎛
∂ ⎞
1
∂2
где Λ =
.
⎜ sin θ
⎟+
sin θ ∂θ ⎝
∂θ ⎠ sin 2 θ ∂ϕ 2
Для получения этого выражения используют следующую систему уравнений:
⎧ ∂ ∂x ∂ ∂y ∂ ∂z ∂
⎪ ∂r = ∂r ∂x + ∂r ∂y + ∂r ∂z
⎪
⎪ ∂
∂x ∂ ∂y ∂ ∂z ∂
=
+
+
⎨
θ
θ
x
θ
y
∂
∂
∂
∂
∂
∂θ ∂z
⎪
⎪ ∂
∂x ∂ ∂y ∂ ∂z ∂
=
+
+
⎪
⎩ ∂ϕ ∂ϕ ∂x ∂ϕ ∂y ∂ϕ ∂z
(4.15)
Решая эту систему уравнений, найдем выражения для производных по x, y, z через производные по r, θ,
φ. Далее, подставляя их в лапласиан, придем к выше выписанному выражению.
Таким образом, уравнение Шрёдингера для одноэлектронного атома с зарядом ядра Z в
сферических координатах примет следующий вид:
⎡ 1 ∂ ⎛ 2 ∂ ⎞ Λ
Z⎤
⎢− 2r 2 ∂r ⎜ r ∂r ⎟ − 2r 2 − r ⎥ Ψ (r , θ , ϕ ) = εΨ (r , θ , ϕ ) . (4.16)
⎝
⎠
⎣
⎦
Для разделения переменных запишем волновую функцию как произведение радиальной функции на
угловую функцию.
Ψ ( r , θ , ϕ ) = R ( r ) ϒ (θ , ϕ ) .
Подставляя волновую функцию в виде произведения в уравнение (4.16), приходим к двум уравнениям:
Λϒ (θ , ϕ ) = λϒ (θ , ϕ ) ;
(4.17)
⎡ 1 ∂2 1 ∂
λ
Z⎤
(4.18)
⎢− 2 ∂r 2 − r ∂r − 2r 2 − r ⎥ R(r ) = εR(r ) .
⎣
⎦
Заметим, что уравнение (4.17) не включает потенциал взаимодействия электрона с ядром, т.е. угловая
часть ϒ (θ , ϕ ) не зависит от заряда ядра одноэлектронного атома и записывается одинаково для всех
одноэлектронных атомов или ионов.
2
2
Оператор Λ = − M , где M - оператор квадрата момента импульса (количества движения).
Под моментом импульса в механике понимают векторное произведение радиус-вектора r ,
проведенного от некоторой избранной точки к частице, на импульс p :
M = ⎡⎣ r p ⎤⎦ .
(4.19)
Особое значение эта величина имеет для случая движения частицы в поле центральных сил, так как она
является интегралом движения. Операторы проекций момента импульса на оси координат в
соответствии с определениями векторного произведения и оператора импульса есть:
⎛ ∂
∂ ⎞
M x = px y − p y z = i ⎜ z − y ⎟ ;
∂z ⎠
⎝ ∂y
∂ ⎞
⎛ ∂
M y = px z − pz x = i ⎜ x − z ⎟ ;
∂x ⎠
⎝ ∂z
(4.20)
⎛ ∂
∂ ⎞
M y = p y x − px y = i ⎜ y − x ⎟ .
∂y ⎠
⎝ ∂x
Оператор квадрата момента импульса записывается как:
2
M = M x2 + M y2 + M z2 =
2
2
2
⎧⎪⎛ ∂
∂ ⎞ ⎛ ∂
∂ ⎞ ⎛ ∂
∂ ⎞ ⎫⎪
− ⎨⎜ z − y ⎟ + ⎜ x − z ⎟ + ⎜ y − x ⎟ ⎬ .
∂z ⎠ ⎝ ∂z
∂x ⎠ ⎝ ∂x
∂y ⎠ ⎪⎭
⎪⎩⎝ ∂y
При переходе от декартовой системы координат к сферической получаем
17
⎛ 1 ∂ ⎛
∂ ⎞
∂2 ⎞
1
+
M 2 = −Λ = − ⎜
sin
θ
⎟.
⎜
⎟
∂θ ⎠ sin 2 θ ∂ϕ 2 ⎠
⎝ sin θ ∂θ ⎝
(4.21)
2
Решение уравнения на собственные значения для оператора M хорошо известны. Это так
называемые шаровые (или сферические) функции ϒ lm (θ , ϕ ) , которые определяются двумя целыми
числами l и m. При этом l может принимать любые целые значения от 0 до ∞, а m = 0, ±1, ±2... ± l .
2
Собственные значения оператора M принимают значения l(l+1), то есть:
2
M ϒlm (θ ,ϕ ) = l ( l + 1) ϒlm (θ ,ϕ ) ,
где ϒ lm (θ , ϕ ) =
m
Знаком
(4.22)
( l − m )!( 2l + 1) P cosθ e
(
)
( l + m )!4π
m
imϕ
l
.
обозначено абсолютное значение числа m. Функции
Pl m ( cos θ )
называются
присоединенными полиномами Лежандра и определяются как:
Pl
m
( cosθ ) = (1 − cos2 θ )
m/2
dm
Pl ( cosθ ) ,
(d cos θ ) m
(4.23)
l
1
dl
cos 2 θ − 1) .
где Pl ( cos θ ) = l
l (
2 l ! d ( cos θ )
Выражения для шаровых функций с l=0,1,2 имеют следующий вид:
ϒ00 =
1
4π
;
ϒ10 =
3
3
cosθ; ϒ1±1 = ∓
sinθ e±iϕ ;
4π
8π
ϒ20 =
5
(3cos2 θ −1) ; ϒ2±1 = ∓ 815π sinθ cosϕeimϕ ; ϒ2±2 = 3215π sin2 θei2ϕ .
16π
(4.24)
Вернемся к решению уравнения (4.18) для радиальной части волновой функции R(r). С учетом
условия (4.22), это уравнение перепишем в форме:
⎡ 1 ∂ 2 1 ∂ l (l + 1) Z ⎤
−
−
− ⎥ R(r ) = εR (r ) ,
(4.25)
⎢−
2
r ∂r
r⎦
2r 2
⎣ 2 ∂r
где l принимает значения l = 0,1,2,…,∞. При решении этого уравнения необходимо отдельно
рассмотреть два случая:
1. ε>0, что соответствует непрерывному спектру энергий;
2. ε<0, что соответствует дискретному спектру.
.
k
Z2
x и ε = 2 . Тогда уравнение (4.25)
В первом случае сделаем следующие подстановки r =
2Z
2k
принимает вид
⎡ ∂ 2 2 ∂ ⎛ 1 k l (l + 1) ⎞⎤
⎟⎥ R ( x ) = 0 .
⎢ 2 + x ∂x + ⎜ 4 + x −
x 2 ⎠⎦
⎝
⎣ ∂x
Для больших значений x это уравнение будет иметь асимптотический вид
⎛
⎝
Решениями его будут функции R ( x ) = c ⋅ exp⎜ ± i
∂ 2 R(x ) 1
+ R(x ) = 0 .
∂x 2
4
x⎞
⎟.
2⎠
Таким образом, для положительных значений ε существуют решения уравнения Шредингера,
соответствующие
непрерывному
ряду
собственных
значений,
отвечающих
ионизации
водородоподобного атома.
18
Для второго случая сделаем следующие подстановки r =
Z2
n
x и ε = − 2 . Тогда уравнение
2Z
2n
(4.25) принимает вид
⎡ ∂ 2 2 ∂ ⎛ 1 n l (l + 1) ⎞⎤
⎟⎥ R ( x ) = 0 .
⎢ 2 + x ∂x + ⎜ − 4 + x −
2
x
∂
x
⎝
⎠⎦
⎣
⎛ x⎞
⎟ , причем функция u(x) должна
⎝ 2⎠
Будем искать решение этого уравнения в виде: R ( x ) = u ( x )x l exp⎜ −
удовлетворять дифференциальному уравнению:
d 2u
du
+ (2l + 2 − x ) + (n − l − 1)u = 0 .
2
dx
dx
Если положить β = 2l + 1 и α = n + l (n ≥ l + 1) , то это уравнение совпадает с уравнением для
x
присоединенных полиномов Лагерра, поэтому радиальные волновые функции должны иметь вид
⎛ x⎞
R ( x ) = cx l exp⎜ − ⎟ L2nl++l1 ( x ) , где c-нормирующий множитель, а n и l – целые положительные числа.
⎝ 2⎠
Таким образом, собственные значения энергии водородоподобного атома (εn) будут зависеть
только от значения введенного нами числа n, называемого главным квантовым числом, и будут
определяться формулой:
1 Z2
εn = −
, где n=1,2,3,…,∞,
(4.26)
2 n2
Z
а собственные функции R(r ) ≡ Rnl (r ) - от n и l. Если ввести новую переменную σ = r (для атома
n
r
водорода σ = ), то выражения для функций Rnl ( r ) при n = 1, 2,3 ; l = 0,1, 2 будут следующими:
n
3
R1s = R10 = Z 2 2e −σ ;
1 ⎛ 1 ⎞ −σ 2
⎜1 − σ ⎟ e ;
2 2⎝ 2 ⎠
3
1 ⎛ 2
2 2 ⎞ −σ 3
R3 s = R30 = Z 2
⎜1 − σ + σ ⎟ e ;
27 ⎠
3 3⎝ 3
3
1
−σ
R2 p = R21 = Z 2
σe 2;
2 6
3
8
⎛ 1 ⎞ −σ
R3 p = R31 = Z 2
σ ⎜1 − σ ⎟ e 3 ;
27 6 ⎝ 6 ⎠
3
4
−σ
R3d = R32 = Z 2
σ 2e 3 .
81 30
R2 s = R20 = Z
3
2
(4.27)
Плотность вероятности нахождения электрона в точке с координатами (r, θ, ϕ) определяется
как:
ρ (r , θ , ϕ ) = Ψ * (r , θ , ϕ ) ⋅ Ψ (r , θ , ϕ ) = Rnl* (r )Υlm* (θ , ϕ )Rnl (r )Υlm (θ , ϕ )
и вероятность найти электрон в объеме dV = dxdydz = r 2 sin θ drdθ dϕ будет иметь вид:
dw = ρ ( r ,θ , ϕ ) dV .
(4.28)
Вероятность нахождения электрона в конечном объеме будет интегралом от вероятности dw по
заданному элементу объема:
wv =
r2 θ 2 ϕ2
∫ ρ ( r ,θ ,ϕ ) dV = ∫ θ∫ ϕ∫ ρ ( r ,θ ,ϕ ) r
V
r1
1
2
sin θ drdθ dϕ . (4.29)
1
Определим вероятность w(r) обнаружить электрон в шаровом слое толщины dr на расстоянии r от
ядра.
19
r
dr
Для этого проинтегрируем dw по переменным θ и ϕ .
π 2π
w (r ) = ∫
∫ Ψ ( r, θ , ϕ ) Ψ ( r, θ , ϕ ) =
*
0 0
π 2π
= Rnl* ( r ) Rnl ( r ) r 2 dr ∫
∫ ϒ (θ , ϕ ) ϒ (θ , ϕ ) sin θ dθ dϕ = ⎡⎣ R ( r )⎤⎦
*
lm
lm
nl
2
r 2 dr.
0 0
(4.30)
Таким образом, получаем функцию вероятности нахождения электрона на расстоянии r от ядра,
которую принято называть радиальной функцией распределения. На рисунках 4.1 и 4.2 приведены
графики радиальных функций Rnl ( r ) и радиальных функций распределения электронного заряда для
n = 1, 2 ; l = 0,1 .
Рис.4.1 Графики радиальных функций R(r) атома водорода
Рис.4.2. Графики радиальных функций распределения электронной плотности в атоме водорода.
В соответствии с выражением (4.26) для водородоподобного атома энергия уровня зависит
только от значения главного квантового числа, так, что все функции с разными значениями l и m для
заданного значения n относятся к одному и тому же энергетическому уровню, т.е. этот уровень
20
вырожден n раз по l, так как l меняется от 0 до n-1 и 2m+1 раз для каждого значения m ( m = −l ÷ +l ) .
Наличие
вырождения
по
m
imϕ
ϒlm (θ , ϕ ) = N ⋅ Pl
( cos θ ) e
l
к
и m позволяет перейти от комплексных функций
действительным, поскольку любая линейная комбинация
вырожденных волновых функций также является собственной функцией гамильтониана.
Рассмотрим случай l=1, где имеются три функции:
3
Υ10 =
cos θ ;
4π
3
3
sin θ ⋅ e iϕ ; Υ1−1 =
sin θ ⋅ e −iϕ ,
(4.31)
4π
4π
из которых одна вещественная и две комплексные. Если вместо функций ϒ11 и ϒ1−1 перейти к
Υ11 =
1
1
( ϒ11 + ϒ1−1 ) и
( ϒ11 − ϒ1−1 ) , записывая экспоненциальные части через
2
i 2
синусы и косинусы переменной φ с помощью формулы Эйлера: eimϕ = cos ( mϕ ) + i sin ( mϕ ) , то новые
линейным комбинациям
функции, которые будем обозначать как Px , Py и Pz запишутся следующим образом:
Px =
1
(Υ11 + Υ1−1 ) = 6 sin θ ⋅ cos ϕ ;
2
2 2π
Py =
1
(Υ11 − Υ1−1 ) = 6 sin θ ⋅ sin ϕ ;
2 2π
i 2
6
cos θ .
(4.32)
2 2π
Соответственно, с учетом радиальной части волновые функции для l=1 и главного квантового числа n
запишутся как:
Pz = Υ10 =
Rn ( r )
R (r )
⋅ x;
( r sin θ cos ϕ ) =
r
r
R (r )
R (r )
Ψ y = Rn ( r ) sin θ cos ϕ = n
⋅ y;
( r sin θ sin ϕ ) =
r
r
R (r )
R (r )
Ψ z = Rn ( r ) cos θ = n
⋅ z.
( r cos θ ) =
r
r
Ψ x = Rn ( r ) sin θ cos ϕ =
(4.33)
Таким образом, три волновые функции вырожденного уровня l=1 являются произведениями полинома
R (r )
на декартовы координаты x, y, z. Если нарисовать графики зависимости волновых функций Ψ nx ,
r
Ψ ny и Ψ nz от угла θ , то они будут выглядеть следующим образом: орбиталь Ψ nz имеет вид объемной
косинусоиды, вытянутой вдоль оси Z, и знак функции определяется знаком cos θ . При этом Ψ nz
обращается в нуль при θ = 90 , т.е. плоскость xy является ее узловой плоскостью. Сечение Ψ nz в
плоскости zx (или zy) иметь вид, изображенный на рисунке 4.3а. Аналогичные диаграммы для Ψ nx и
Ψ ny представлены на рисунках 4.3b и 4.3с.
z
y
Ψnz
y
Ψny
Ψnx
+
+
x
-
+
x
x
-
-
a
c
b
Рис.4.3 Сечения граничных поверхностей функций Ψ 2 px , y , z .
21
Графическое представление волновой функции позволяет наглядно представить вероятность
найти электрон в некоторой точке пространства, так как плотность вероятности определяется
квадратом модуля функции - Ψ *Ψ . Граничные поверхности волновых функций, т.е. поверхности,
соединяющие точки одинаковых значений функций Ψ 2 px , y , z в плоскостях xy, xz и yz рисуют по
аналогии с изображением высоты местности на географических картах. Например, контуры функции
Ψ 2 pz будут иметь следующий вид в плоскости xz:
0.2
0.2
0.2
0.4
0.4
0.2
Рис.4.3. Граничные поверхности атомных орбиталей типов s, p,d и f.
§5. Спин.
Предположение о наличии у электрона собственного механического момента (спина) впервые
было введено в квантовую механику Гаудсмитом и Уленбеком в 1925 г. для объяснения
экспериментальных результатов.
22
Прямым указанием на существование спина является опыт Штерна и Герлаха (1922 г.) по
расщеплению пучка атомов серебра в магнитном поле, показавший, что атомы серебра имеют
ненулевой магнитный момент в s-состоянии. Компоненты магнитного момента атома в направлении
внешнего магнитного поля имеют дискретные значения, что можно объяснить, если предположить, что
электрон обладает собственным механическим моментом (спином), который проявляет себя так, как
если бы электрон был вращающимся твердым телом. Величина возможной проекции этого
дополнительного момента равна
2
(в атомных единицах 1 ). Соответствующий этому магнитный
2
момент, который определяет энергию взаимодействия электрона с магнитным полем равен μ =
e
.
2mc
При этом отношение магнитного момента к механическому вдвое больше для спина электрона, чем для
орбитального движения.
Гипотеза спина электрона также позволяет объяснить мультиплетную структуру щелочных и
щелочноземельных атомов, эффект Зеемана и т.д.
В квантовой механике вводится линейный самосопряженный векторный оператор спина S с
компонентами S x , S y , S z , для которых выполняются те же коммутационные соотношения, что и для
компонент орбитального (углового) момента Lx , Ly , Lz :
S x S y − S y S x = iS z ;
S y S z − S z S y = iS x ;
(5.1)
S z S x − S x S z = iS y .
По аналогии с угловым моментом, спиновый момент электрона должен характеризоваться
двумя величинами: абсолютным значением вектора спина s и проекцией этого вектора на
какое-либо выбранное направление. Для определения абсолютного значения вводится оператор
S 2 = S ⋅ S = S x2 + S y2 + S z2 ,
(5.2)
а в качестве проекции выбирают Sz.
Обычно в нерелятивистской квантовой механике вводят две спиновые функции α и β,
собственные для оператора Sz с собственными значениями +
2
и −
2
(+
1
1
и −
в атомных
2
2
единицах). Так что
1
S zα (σ ) = α (σ ) ;
2
1
S z β (σ ) = − β ( σ ) .
2
(5.3)
При этом, если ввести формально интегрирование по спиновой переменной, которая в случае
электрона принимает два значения, то получим:
α α = β β = 1,
α β = β α =0.
(5.4)
Вместо операторов Sx и Sy можно ввести два других оператора:
S + = ( S x + iS y ) , S − = ( S x − iS y ) .
(5.5)
Эти операторы называются операторами повышения и понижения, так как они, действуя на функцию с
собственным значением Sz для оператора S z дают собственную функцию S z с собственным значением
на единицу больше или меньше, соответственно.
S+ β = α
и S + α = 0;
S− β = 0 и S− α = β .
(5.6)
23
Действительно, пусть оператор
(S
x
+ iS y ) действует на спиновую функцию β :
(S
x
+ iS y ) β .
Подействуем на эту функцию оператором S z и учтем коммутационные соотношения.
S z S x = iS y + S x S z , S z S y = −iS x + S y S z ;
S z ( S x + iS y ) β = ( S z S x + iS z S y ) β = ( iS y + S x S z + S x + iS y S z ) β =
=
{( S
x
}
+ iS y ) + ( S x + iS y ) S z β = ( S x + iS y ) (1 + S z ) β =
Также можно показать, что
1
⎛ 1⎞
= ( S x + iS y ) ⎜ 1 − ⎟ β = ( S x + iS y ) β .
2
⎝ 2⎠
(S
Т.е.
(S
x
x
+ iS y ) β = α ,
− iS y ) α = β
и
(S
(S
x
x
+ iS y ) α = 0.
− iS y ) β = 0.
В общем случае, если функция является собственной для оператора S z с собственным значением m, то
оператор
(S
x
+ iS y ) переводит ее в собственную функцию S z с собственным значением m+1,
соответственно оператор
(S
x
− iS y ) , действуя на эту функцию, переводит ее в волновую функцию,
собственную для S z с собственным значением m-1. По этой причине эти операторы называют
операторами повышения и понижения, соответственно.
Оператор S 2 можно представить в более удобной форме, без S x2 и S y2 , используя операторы
повышения и понижения - S + и S − . Для этого рассмотрим произведение операторов S+S-:
(S
x
+ iS y )( S x − iS y ) = S x2 + S y2 + iS y S x − iS x S y =
S x2 + S y2 + i ( −iS z ) = S x2 + S y2 + S z .
С другой стороны, S 2 = S x2 + S y2 + S z2 . Подставляя S x2 + S y2 = S + S − − S z , получаем:
S 2 = S + S − + S z2 − S z
или S 2 = S − S + + S z2 + S z . (5.7)
Для многоэлектронных систем вводится соответствующий набор операторов, отвечающих
суммарному спину системы. Если у нас имеется система из N частиц, то операторы будут иметь вид:
S (1, 2,..., N ) = S (1) + S ( 2 ) + ... + S ( N ) и
S z (1...N ) = S z (1) + S z ( 2 ) + ... + S z ( N ) ;
S x (1...N ) = S x (1) + S x ( 2 ) + ... + S x ( N ) ;
S y (1...N ) = S y (1) + S y ( 2 ) + ... + S y ( N ) .
S + (1, 2,..., N ) = S + (1) + S + ( 2 ) + ... + S + ( N )
(5.8)
В общем случае можно показать, что операторы
гамильтонианом
H,
т.е.
выполняется
соотношение
S2 = S ⋅S
и
Sz
[A, H ] = AH − HA = 0 ,
коммутируют с
где
A = Sz , S 2 .
Следовательно, волновая функция, собственная для H , должна быть также собственной функцией
этих операторов.
Таким образом, наряду с тремя пространственными переменными электрона (x,y,z), появляется
еще одна переменная σ, определяющая спин электрона и принимающая дискретные значения +
1
− .
2
24
1
и
2
Download