книги из ГПНТБ / Хабердитцл, В. Строение материи и химическая связь
.pdf196 Часть II. Введение в квантовохимические расчеты
линейных комбинаций можно построить молекулярные орбитали.
В заключение этой главы следует указать на серьез ные недоразумения по поводу распределения заряда в мно гоэлектронных атомах, которые часто обнаруживаются
у студентов. Рассмотрим в |
качестве |
примера атом N |
с тремя 2д-электронами (рх, |
р у и p z). |
Не следует путать |
симметрию этих орбиталей с симметрией распределения заряда. Последняя характеризуется суммой г|з^ + +
+которая равна cos2(psin2& -f sin2(psin2& -f- cos2&=
=1. Эта сумма не зависит от &и ср, поэтому можно сде лать вывод о сферически симметричном распределении заряда. Это справедливо по отношению ко всем атомам
ссимметрично заполненными оболочками. Следует отме тить, что в 5/73-гибридизованном атоме углерода (см. разд. 6.3.1, часть I) распределение заряда сферически симметричное, а не тетраэдрическое, как часто предпола гают вследствие смешения понятий симметрии орбиталей
исимметрии распределения заряда.
3.2.Введение атомарных единиц
Для упрощения способов записи квантовохимических соотношений очень удобно использовать атомарные еди ницы. В качестве единицы длины служит уже введенный
выше р а д и ус Б о р а а 0 — Л2 = 0,52917-10 8 см (назы
ваемый в последнее время просто «1 Бор»), В качестве единиц массы пг и заряда е используют
массу и заряд электрона; единицей действия служит ве
личина h = |
= 1,0544-10~27 эрг-с. |
В |
качестве еди- |
|
|
|
|
|
Таблица 5 |
Соотношение между |
единицами |
энергии |
||
|
Атомная энергия |
эВ |
|
ккал/моль |
Атомная энергия |
1 |
27,210 |
|
627,71 |
эВ |
3,6752-10“2 |
1 |
|
23,069 |
ккал/моль |
1,5931-10-® |
4,3348-1О-2 |
1 |
4. Введение в рассмотрение иона Н% методом МО |
197 |
в*
ницы энергии выбрана величина — = 2Eh = 27,210 эВ, “О
т. е. удвоенная энергия ионизации атома водорода; для пересчета единиц энергии полезна табл. 5.
4. ВВЕДЕНИЕ В РАССМОТРЕНИЕ ИОНА Hj МЕТОДОМ МО
Чтобы познакомиться с основами метода МО, мы уже проводили краткое качественное обсуждение уравнения
Шредингера для иона Нг+. Теперь познакомимся с осно вами количественной трактовки* при помощи прибли женного вариационного метода. Преобразуем уравнение
(разд. 6.2.4, часть I)
w + T - ( £ - T + 7 T + l r ) ' i' = o |
(|) |
в форму, обычно применяемую в квантовой химии: |
|
H'F = £.'F . |
(2) |
Специфика уравнения (1) заключается в особой форм оператора Гамильтона Н:
Н: |
h2 |
д + 4 — |
га |
— . |
(3) |
|
2т |
R |
гв |
|
Умножим уравнение (2) слева [в левой части уравнения (2) стоит оператор!] на Т и после интегрирования по всему объему получим следующее выражение для энергии:
J* УНУЩ)
Е |
(4) |
* При этом для нас принципиально важно приближение Бор на — Оппенгеймера: движение ядра и электронов должно быть раз делимым, чтобы каждому данному положению ядра соответствовала определенная энергия электронов. Вследствие относительно боль шой массы ядра ошибка очень невелика (например, по ван Флеку, для Н j она составляет < 0,0075 эВ). Эмпирическое правило для порядков величин энергий переходов дает следующие значения: электронные переходы, 10°—101 эВ; колебательные переходы,
Ю-1 эВ; торсионные переходы, 10-а эВ; вращательные переходы, 10-3 эВ.
198 Часть II. Введение в квантовохимические расчеты
Мы представили Y в виде линейной комбинации дейст вительных атомных орбиталей фА и фв. Чтобы можно было воспользоваться вариационным исчислением, введем в
качестве варьируемых параметров коэффициенты линей
ных комбинаций (вместо |
фА, фв запишем ф, срА, фв) |
|
Ф =сАфА+ св<рв. |
(5) |
|
Запишем условие минимизации Е |
|
|
дЕ |
дЕ __0 |
|
дед |
дсв |
|
Сначала подставим уравнение (5) в уравнение (4); тогда
| (сАФа + свФв) Н (сдфд + свФв) dv
Е =
j* (саФа + свФв)2 dv
При умножении получим
с\ j фдНфдйс + 2сдев j фдНфв* + св | ФвНфв^и
САj 4>\dv + 2сАсв JФаФв^о+ Свj q>2Bdv
Введем следующие сокращения и обозначения:
ГФАНфв^о= НАВ= рАВ—резонансный интеграл;
фАНфАйп = НАА — аА—кулоновский интеграл;
фАфв^о = SAB—интеграл перекрывания.
1 '
Имеем
с\НАА + 2сАсвЯАв + св Явв
|
Е = |
(6) |
|
|
СА Saa + 2cacbSAb -f св Sbb |
Запишем частную |
производную по сА и приравняем ее |
|
к нулю |
|
|
1 |
дЕ _ (сдЯаа + сдсвЯдв) (сАSaa + 2cacbSab + св Sbb) |
|
2 |
dcA |
(c A S Aa + 2cacbS ab + св 5 вв)а |
( ca Я д а + 2 с а св Я а в + c b h b b ) (c a ^ a a + св $ а в )
(caSaa + 2cbcaSab + cB Sbb)1
4. Введение в рассмотрение иона Н J методом МО |
199 |
Умножим на выражение, стоящее в знаменателе в скоб ках, проведем сокращение и приведение подобных членов с учетом уравнения (6) и получим:
са ( ^ а а —ESaa) + св (На в —ESAB) —0. |
(7а) |
При проведении аналогичной процедуры с д—- = 0 по лучаем
са (На в — ESAB) + св (Нвв —ESBв) = 0. |
(76) |
Полученная система линейных однородных уравне ний (вековых уравнений) имеет нетривиальные решения только тогда, когда обращается в нуль определитель, составленный из коэффициентов. Учитывая условие нор мировки (5аа = SBB = 1) и тот факт, что здесь аА — = а в = а, имеем
а —Е |
Р а в — E S a b |
P ab ESab |
а —Е |
Отсюда следует два решения:
Es |
g + Р а в |
и Еas |
_ « — Ра в |
(8) |
1+ 5Ав |
1 - S ab ' |
Коэффициенты сА и гв спределяют подстановкой Es или Eas в уравнения (7). Если, креме того, учесть, что
и фА должны быть нормированы, то для симметричного
решения
1
са — с в —
/ 2 ( 1 + S a b )
и для антисимметричного решения
С к — — С о ------ т -...- |
• |
/ 2 ( 1 - S AB)
Физический смысл полученных таким образом решений мы уже обсуждали на стр. 102.
Расчет необходимых интегралов выходит за рамки данного курса. Мы сообщаем здесь только результат. Для этого запишем в атомарных единицах для h, m и е оператор Гамильтона (3)
Н = — f |
Д - - i - |
R |
(9) |
2 |
r A Гв |
200Часть II. Введение в квантовохимические расчеты
иуравнения для разделенных атомов водорода А и В
^ |
^-А--- |
— ^ I 1sa ) — ^ h I 1sa ) ’ |
|
|
( 10) |
^ --- |
Y А ---- |
I 1% ) = Ен I lsB > . |
При этом | lsA> = | lsB > = фА = фв обозначают нор мированные функции основных ls-состояний атомов во дорода (см. табл. 3); Ен — соответствующие им значения энергии. Введем еще следующие обозначения интегра лов:
J a a = J Фа ( ----- Фа^ = < |
1sa |
_1_ |
lsA >, |
||
гв |
|
||||
*^ab = j Фа ^ |
^~|.Фв^у — ( 1sa |
гв |
lsB ) . |
||
Подставим эти соотношения |
в уравнения |
(8): |
|||
£s — Ан |
|
•^аа + ^Ав |
|
||
|
1+ Sab |
|
|||
|
|
|
|
||
|
|
R |
J A A |
J A B |
|
A a s |
— А н |
1 —Sab |
|
Таким образом, для расчета энергии нам нужны сле дующие три интеграла: JАА, J АВ и SAB. Результаты вы числения интегралов мы позаимствуем из соответствую щих литературных источников. Чтобы обеспечить свя зывание, слагаемое в Es, состоящее из интегралов, долж
но компенсировать член —-, характеризующий куло
новское отталкивание двух ядер. Результаты расчета ин тегралов таковы:
Адв = е~R ^ 1 + ■ft+ —j- j .
JAA=— ^ t t ~ e - 2«(l + /?)],
■^А в=-в-Л(1+/г).
4. Введение в рассмотрение иона Н а+ методом МО |
201 |
На рис. 4 показаны кривые энергии для симметрич ного и антисимметричного состояний, полученные при
помощи |
этих |
интегралов. |
|
|
||
Мы |
получили |
энергию, составляющую 60% экспе |
||||
риментального |
значения |
энергии |
связи |
Е — Ен |
||
(64,1 ккал/моль). |
Улучшение нашего очень простого под |
|||||
хода (увеличение |
числа |
варьируемых |
параметров, вве- |
Рис. 4. Энергия Н2+ как функция межъядерного расстояния.
дение эффективных зарядов ядер) приводит к превосход ному совпадению. Так, Финкельштейну и Горовицу в 1928 г. путем вариации величины Z удалось прийти к энергии связи в 51,9 ккал/моль и даже получить правиль ное равновесное расстояние — 1,06 А. Промежуточные значения эффективного заряда ядер, между 1 и 2, можно пояснить следующим образом. Реализуется определенное промежуточное состояние между случаем с очень боль шим межъядерным расстоянием (при котором эффектив но действует только один заряд, Z « 1) и случаем объеди ненного ядра (Z — 2). На рис. 5 показана зависимость 2Эфф от расстояния между ядрами (в атомарных единицах). Для равновесного расстояния получают значение Z3dxb =
=1,24.
Дальнейшим усовершенствованием (Дикинсон) яв
ляется такой вполне допустимый прием, как использо
202 Часть II. Введение в квантовохимические расчеты
вание вместо сферически симметричных ls-функций ато ма водорода, гораздо более близких к реальности функ ций типа
Фа (В) = Фа (в) + ^хфд <в)
(где ось х расположена вдоль связи А—В, а А, — варьи руемый параметр), которые учитывают, если это необхо-
Рис. 5. Зависимость эффективного заряда ядра от расстояния между ядрами для иона H t в методе МО ICouIson С. A ., Trans. Faraday
Soc., 33, 1479 (1937)].
димо, поляризующее влияние других ядер. И действи тельно, такой путь приводит практически к эксперимен тальным значениям, или к величинам, которые получа ются при точном расчете; трактовка последнего здесь
еще преждевременна. Первый точный расчет иона HJ
4. В в е д е н и е в рассм от рение и о н а / / « м ет одом М О |
203 |
принадлежит Хиллераас (1931); она использовала эллип тические координаты (гА ± гв, азимутальный угол отно сительно R), а также подход ЛКАО (Барроу, 1927).
Из этого не следует, однако, что варьировать можно только линейные параметры, как это делается в методе ЛКАО. Так, например, в 1935 г. Джеймс показал, что очень хороших результатов можно достичь, используя функцию
^ = e - M ' A + ' B) [ l + C i (rA- r B)*]
сдвумя варьируемыми параметрами сх и с2.
4.1.Запрет пересечения
Добавим еще одно замечание общего характера, отно сящееся к энергиям, которые находят при помощи вариа ционного метода ЛКАО в виде функций межъядерных расстояний. Исходя из вековых уравнений (7), обозна
чим Еа и |
Ев энергии, соответствующие функциям |
ф |
А |
|||||
и фв: |
|
|
|
|
|
|
|
|
|
|
|
ФаНФа* |
J ФвНфв йт |
( |
И |
) |
|
|
|
|
J |
ф\4т |
£в |
|||
|
|
|
| Фв* |
|
|
|
||
В предположении, что функции нормированы, из |
||||||||
уравнений |
(7) |
получаем |
|
|
|
|
||
|
|
ск (Еа |
Е) + св (НАв—ESab) = 0, |
|
|
|
||
|
|
са (ЕАв—ESab) -(- св (Ев— £) —0. |
(12) |
|||||
Решая |
каждое уравнение |
относительно |
получаем |
|||||
далее |
следующее |
квадратное уравнение: |
|
|
|
|||
|
|
(Е |
Еа)(Е Ев) |
(Яав ESab)2 — 0. |
(13) |
Рассмотрим его левую часть f(E) и построим график, предполагая, что ЕА < Ев (рис. 6). Независимо от кон кретных значений ЕА, Ев, НАВ и 5АВ приходим к сле дующему общему результату: оба значения Е, получен ные при помощи вариационного метода ЛКАО, т. е. нули функции f(E) лежат ниже, чем Ел, и выше, чем Ев.
204 |
Часть II . В в е д е н и е в к ва н т о во хи м и ч ески е расчеты |
Следствием этого оказывается очень важное для двух атомных молекул явление: если рассмотреть ЕА и Ев
как функции R (последний входит в Н в виде |
то |
получается зависимость, изображенная схематически на рис. 7 [E(R) — пунктирная кривая].
f(E)
Если теперь предположить, что при изменении R при определенном его значении может наступить такая ситуация, когда ЕА = Ев и далее Ев < ЕА, то при этом
кривые Еа и Ев пересекутся, тогда как кривые Е = f(R) по-прежнему не пересекаются (рис. 8). Это составляет сущность важного правила непересечения (запрет пере сечения). (Заметим еще, что это правило справедливо только тогда, когда срА и срв обнаруживают определенные свойства симметрии — это условие станет понятным после изучения теории групп.) Запрет пересечения играет большую роль при анализе диссоциации молекул.
4. В в е д е н и е в рассм от рение и о н а Н д мет одом М О |
205 |
4.2. Максимальное перекрывание
Анализ вековых уравнений (7) позволяет получить еще одну характеристику метода МО Л КАО, имеющую общий характер, к которой мы придем здесь чисто ка чественно. Сравнивая уравнения (9) и (10), можносделать вывод, что молекулярный оператор Гамильтона и атом ные операторы Гамильтона для области пространства обоих отдельных атомов подобны. Это означает, что ЕА и Ев близки энергиям отдельных атомов. Перепишем уравнения (2) в форме
(Еа~ Е ) + ^ - ( Я ав- £ 5 ав)= 0 .
|
|
|
|
|
(14) |
|
(^ ав—ESab) + |
(Ев |
Е) = 0. |
||
Теперь покажем, |
что метод Л КАО имеет смысл толь |
||||
ко тогда, |
когда ЕА и Ев не слишком различаются. Пред |
||||
положим, |
что Еа |
Ев, |
тогда по уравнению (11) и срА < |
||
<С фв. Отсюда следует, |
что |
SAB, |
НАВ и квадратичный |
член в уравнении (13) также становятся малыми. Но тогда в уравнении (13) либо величина Е — ЕА, либо величина Е — Ев должна быть также мала. Что отсюда следует? Рассмотрим первый случай (для второго справедливо аналогичное рассуждение) и подставим в уравнение (13) для Е значение ЁА (за исключением только первого чле на в скобках, иначе он обращается в нуль). После преоб разований получим
исоответственно
Е= Е п
Ша в — E AS Ab )2
Ев — Е а
(15)
( Н а в — E b S a b )2
Е в — Е а
Из подстановки первого уравнения (15) в первое урав нение (14) следует, что
£в_ _ — (ЯАв — EASа в )
с а |
Е в — Е а |