Расчет электронных характеристик молекул полуэмпирическим методом Хюккеля
.pdf• |
kC-B. = 0.7 |
|
h( B ) = -1 |
|
|
• |
kC-N. = 1 |
|
h( N )= 0.5 |
|
|
•• |
kC-N. = 0.8 |
|
h( N ) = 1.5 |
|
|
•+ |
|
kC-N. = 0.7 |
h( N ) = 2 в NO2 группе |
|
|
_________________ |
|
|
• |
kC-O. = 1.41 |
|
h( O ) = 1 |
||
•• |
kC-O. = 0.8 |
|
h(O ) = 2 |
|
|
________________ |
|
|
h(F) = 3 |
kC-F. = 0.7 |
|
h(cl) = 2 |
kC-Cl. = 0.4 |
|
h(Br) = 1.5 |
kC-Br. = 0.3 |
|
_________________ |
|
|
h(Х) = 2 сверх сопряжение с (СН3 = Х) |
kC-Xr. = 0.7 |
|
h(С) = - 0.5 С- связанный с (СН3 = Х) |
|
1.6.1Расчет молекулы формальдегида
Вмолекуле формальдегида атом кислорода поставляет в π систему 1 электрон (гетероатом 1 рода, связь С=О – двойная). Параметры Стрейтвизера равны:
• |
α2 = α0 + 1β |
|
||
h(O ) = 1, |
(93) |
|||
kC-O. = 1.41 |
β12 = |
|
β |
|
2 |
|
Секулярные уравнения в приближении Хюккеля запишутся следующим образом:
С1х + С2 |
|
= 0 |
|
|||
2 |
|
|||||
С1 |
|
|
+ С 2(х +1) = 0 |
(94) |
||
|
2 |
|||||
Х2 |
+ Х − 2 = 0 |
|
||||
Расставляя в порядке возрастания все Х, получим: |
|
|||||
Х1 |
= − 2, Х2 = 1 |
|
Подставляя полученные значения в секулярные уравнения определим набор коэффициентов для каждого Х. Структура двух найденных М.О. будет выглядеть следующим образом:
Ψ1 = 0,46χ1 + 0,885χ2 |
(95) |
|
Ψ2= 0,885χ1 − 0,46χ2 |
|
|
ε1 |
= α + 2β |
(96) |
ε2 |
= α − β |
|
Электронные плотности на атомах и заряды на каждом атоме равны:
23
P11 |
= 2(0,46)2 = 0,42 |
q1 |
= + 0,58 |
|
P22 |
= 2(0,885)2 = 1,58 |
q2 |
= −0,58 |
(97) |
P21π = 2(0,46)(0,885) = 0,814 |
|
|
|
Согласно этому расчету свободная валентность для любого атома может быть определена:
Fµ = 3 − nµ, где nµ - сумма порядков связей окружающих атом µ.
Еполн= ∑ ni • εi = 2(α +2β) = 2α + 4β |
|
(98) |
||||
i |
|
|
|
|
|
|
Индекс свободной валентности Fµ на каждом атоме составляет: |
||||||
F1 = F2 = |
|
− 0,814 = 0, 916 |
|
|
(99) |
|
3 |
|
|
||||
1.6.1.1 Расчет дипольного момента формальдегида |
|
|||||
µХ = ∑ |
qνХν , |
µУ = ∑ qνУν , |
µZ = |
∑ qνZν |
||
ν |
|
|
|
ν |
|
ν |
(100) |
|
|
µХ2 |
|
µУ2+ |
µZ2 |
µ2 |
= |
+ |
||||
(101) |
|
|
|
|
|
|
ZC = 0, RCO = 1,26 Å
µZ = µ = - 0,58 e • 0,16 = -0,67 e Å • 4,8 = -3,12 D
1.6.2 Расчет системы (NR2)2CO на примере молекулы мочевины
(NH2)2CO
В молекуле (NH)2CO каждый атом азота поставляет в π систему по 2 электрону (гетероатом 2 рода, связь С N – оодинарная). Параметры Стрейтвизера в этом случае равны:
•• |
α2 = |
α0 + 2β |
||
h( N ) = 2 |
||||
(102) |
|
|
|
|
kC-N. = 0.8, β12 (CO) = 1,5, |
β13(CN) |
= 1,5β |
||
Секулярные уравнения в приближении Хюккеля запишутся следующим |
||||
образом: |
|
|
||
С1х + С2 |
|
+ С3 + С4 |
= 0 |
|
2 |
|
24
С1 |
2 |
+ С2(х +2) |
= 0 |
(103) |
|
С1 |
+ 0 |
+ С3(х +2) |
= 0 |
|
|
С1 |
+ 0 |
+ 0 + С3(х +2) = 0 |
|
Применяем элементы симметрии (ось С2 проходящая через связь С = О). Соотношение для коэффициентов запишется:
С3 = С4 (симм. относительно Х) С3 = -С4 (антисимм. относительно Х)
С1 = -С1= 0 |
(104) |
С2 = -С2= 0
С учетом этих условий секулярные уравнения перепишутся в виде:
Симметричный случай: |
|
|
|
|||||
С1х + С2 |
2 |
+ 2С3 |
= 0 |
|
(105) |
|||
С1 |
|
2 |
+ С2(х +2) |
= 0 |
|
|
||
С1 |
|
|
+ 0 + С3(х +2) |
= 0 |
|
|
||
Х2 + 2Х − 4 = 0, |
Х12 = − 1 ± |
|
(106) |
|||||
5 |
||||||||
Х1 |
= -3,23, |
|
|
Х2 = -3,23, |
|
|||
Антисимметричный случай: |
|
|||||||
С3(х +2) |
= 0 |
|
|
(107) |
||||
Х3 |
= -2 |
|
|
|
|
|
Расставляя в порядке возрастания все Х, получим:
Х1 = − 3,23, Х2 = − 2, Х3 = + 1,23 (108)
Энергии М.О. равны:
ε1 = α + 3,23β
ε2 = α + 2β (109)
ε2 = α − 1,23β
Структура двух найденных М.О. будет выглядеть следующим образом:
Ψ1 = 0,57χ1 + 0,875χ2+ 0,46(χ3 + χ4) |
|
|
|
|
Ψ2= + 0,707χ2 − 0,5(χ3 − χ4) |
|
|
(110) |
|
Ψ2= + 0,707(χ3 − χ4) |
|
|
|
|
ε1 |
= |
α |
+ |
2β |
(111)
ε2 = α − β
25
Электронные плотности на атомах и заряды на каждом атоме равны:
q1(C) = 2(0,57)2 = 0,33 • 2 = 0,66 |
|
|
|
|||
q2(O) |
= |
2(0,5)2 |
+ |
2(0,707)2 |
= |
1,5 |
(112)
q3(N) = 2(0,46)2 + 2(0,5)2 + 2(0,707)2 = 0,21 • 2 + 1,5 = 1,92
Таким образом, на атомах азота локализован заряд равный 2 − 1,92 =
0,08 e
Общее количество электронов в π- системе данной молекулы 6, поэтому
∑ qν = 6
ν
1.6.2.1 расчет дипольного момента молекулы мочевины (NH2)2CO
RCO = 1,16 Å, |
RCN = 1,4 Å |
(113) |
Направление оси Z совпадает со связью С=О, тогда координату атома N по оси Z можно вычислить
ZN |
= |
1,4 Å • cos |
60 = 0,7 Å |
(114) |
|
|
|
Проекции дипольного момента µХ = µУ = 0, проекция дипольного |
|||
момента по оси Z отлична от нуля. |
|
||
µZ |
= ∑ |
qνZν = 0,08 e • 2(0,7) Å − 1,16 • 0,5 e Å = − 0,692 • 4,8 |
|
|
ν |
|
|
= − 3,32 D |
|
|
|
µ2 = µХ2 + µУ2+ µZ2 |
|
||
ZC = 0, |
RCO = 1,26 Å |
(115) |
µZ = µ = - 0,58 e • 0,16 = -0,67 e Å • 4,8 = -3,12 D
1.7 Индуктивный π эффект в альтернантных углеводородах, атом-атомная поляризация
Индуктивное влияние гетероатома на π электроны (Iπ индуктивный эффект) вызывает поляризацию заряда между помеченными и не помеченными
атомами АУ. То есть заряд qµ* на отмеченных звездами атомах увеличивается, в то время как заряд на не отмеченных атомах qµ - уменьшается. Эффект
26
можно пояснить вычислив изменение заряда δqµ на каждом атоме воспользовавшись теорией возмущений.
Пусть возмущающее действие выражается в появлении дополнительного
малого слагаемого |
в гамильтониане H`. |
|
Возмущение |
H` возникает вследствие замещения атома углерода на |
|
другой гетероатом. Тогда |
|
|
Н = Н0 + Н` |
(116) |
возмущенную М.О. можно записать как не возмущенную волновую функцию и поправку первого порядка
|
∫ϕ |
H`ϕ |
dV |
|
|
|
Ψ`i =ϕ i + ∑ |
i |
j |
|
ϕ j |
(117) |
|
Ei −E j |
||||||
j |
|
|
Индекс i означает номер занятой М.О., индекс j номер не занятой электронами (свободной) М.О. Каждую из этих М.О. можно представить как линейную комбинацию А.О.
Для того чтобы не потерять перекрестные члены в выражениях, где встречается произведение вида ϕ i • ϕ i обе функции необходимо
представить в виде суммы по разным индексам. |
|
|||||||||||||
ϕ |
i |
= ∑ |
С |
µχµ |
ϕ |
j |
= ∑ |
С |
νχν |
|
|
(118) |
||
|
µ |
i |
|
|
|
|
ν |
j |
|
|
|
|
||
Индексы µ, ν означают номер А.О. Сiµ, Сiνχν - есть коэффициенты на |
||||||||||||||
А.О. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя (118) в (117) получим: |
|
|
|
|
||||||||||
|
` |
|
|
|
∫ |
∑C |
χ |
µ |
H ′∑C |
χ |
dV |
|
||
|
|
|
|
µ iµ |
|
|
ν |
jν ν |
|
|
|
|||
Ψ i =ϕ i + |
∑j |
|
|
|
|
|
|
ϕ j |
(119) |
|||||
|
|
|
Ei −E j |
|
|
Необходимо помнить, что возмущение, описываемое дополнительным слагаемым в гамильтониане H` касаются только атома µ. В методе Хюккеля
каждый атом углерода вносит в π систему только 1 электрон, находящийся на 2Р А.О. не участвующей в образовании гибридной связи. Только по этой причине в методе Хюккеля число атомов равно числу А.О. Изменение
кулоновского интеграла αсвязанного с возмущающим действием H` на атоме µ запишутся:
В связи с этим оператор H` будет действовать только А.О. µ находящуюся на атоме µ. На все другие χν оператор Н` не действует, их можно вынести
27
перед оператором и в числителе (119) знаки суммирования исчезнут. Тогда (119) можно переписать в виде:
Ψ`i =ϕ i + ∑ |
CiµC jµδαµ |
• ∑ Сjνχν |
(120) |
|||||||
|
|
|
||||||||
где |
|
j |
|
Ei −E j |
ν |
|
||||
|
|
|
|
|
|
|
|
|
|
|
∫ |
χµ Н`χµdV = δαµ |
|
|
(121) |
||||||
Перепишем еще раз выражение (120) в более удобном виде: |
|
|||||||||
Ψ`i =ϕ i + ∑ |
∑ |
|
CiµC jµC jνδαµ |
•χν |
(122) |
|||||
|
|
|
||||||||
|
|
ν |
j |
|
|
Ei −E j |
|
|
|
|
Введем обозначение |
|
|
|
|
|
|
|
|||
∑ |
|
CiµC jµC jνδαµ |
= Сiν` |
|
|
(123) |
||||
|
|
|
|
|||||||
j |
|
Ei −E j |
|
|
|
|
|
|
|
|
С учетом (122) и (123) выражение (120) перепишется в виде: |
|
|||||||||
Ψ`i =ϕ i + ∑ Сiν` χν |
|
|
(124) |
|||||||
|
|
ν |
|
|
|
|
|
|
|
|
Согласно Хюккелю электронная плотность на атоме вычисляется следующим образом::
qµ = ∑ ni Сiν2 |
(125) |
i |
|
Напомним, что индекс i пробегает занятые М.О. Для закрытой оболочки
ni =2
Изменение заряда на атоме µ можно записать:
δqν = 2 |
∑ (Сiν + Сiν`) 2 |
− 2 |
∑ |
Сiν2 |
(126) |
|||||
|
|
|
i |
|
|
|
i |
|
|
|
Раскрывая скобки, получим: |
|
|
|
|
|
|||||
δ |
qν = 2 |
∑ |
` |
+ 2 |
∑ |
|
` |
) |
2 |
(127) |
|
i |
2Сiν Сiν |
i |
(Сiν |
|
|||||
|
|
|
|
|
|
|
|
|
28
Вторым слагаемым в (127) пренебрегаем в виду порядка малости и с учетом (123) выражение (127) перепишется в виде:
δ |
|
δα |
µ ∑ |
∑ |
CiµCiνC jµC jν |
|
(128) |
|||
|
|
|
|
|||||||
qν = 4 |
|
|
j |
|
Ei −E j |
|||||
|
|
|
|
i |
|
|
|
|||
Отсюда |
|
|
|
|
|
|
|
|||
δ |
δα |
µ = 4 |
∑ |
∑ |
|
CiµCiνC jµC jν |
(129) |
|||
|
|
|
||||||||
qν / |
|
|
j |
|
Ei −E j |
|||||
|
|
|
|
|
i |
|
|
|||
Величина δqν / |
δαµ носит название атом атомной поляризации и |
обозначается πµ,ν. Таким образом, окончательно получим:
πµ,ν.= |
4 i j |
|
∑ ∑ |
(130)
CiµCiνC jµC jν Ei −E j
Воспользуемся полученным результатом и рассчитаем π1,2 π1,3 и π1,1 для молекулы бутадиена.
1.7.1Расчет атом-атомной поляризации в бутадиене
Вмолекуле бутадиена 4 углеродных центра. Используя метод Хюккеля, мы получим 4 М.О. две из которых заняты и две свободны.(см. выше). Следовательно i = 1,2, а j = 3,4.
|
|
|
|
|
C |
µC |
νC |
µC |
ν |
|
|
|
|
|
C µC |
νC jµC jν |
|
|
||||||
π |
|
|
∑ ∑ |
i |
|
i |
|
j |
j |
|
|
|
π |
|
∑ |
1 1 |
|
|
|
|
||||
µ,ν.= 4 |
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
+ |
|
|||||||
|
i |
j |
|
Ei −E j |
|
1,2. = 4( |
j |
|
|
E1 −E j |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
C |
µC νC jµC jν |
|
|
C µC νC µC ν |
|
|
C µC νC |
µC ν |
|
|||||||||||||
+ ∑j |
|
2 |
2 |
|
|
|
|
) = 4( |
|
1 |
1 |
3 |
3 |
|
+ |
1 |
1 |
4 |
|
4 |
+ |
|||
|
|
E2 −E j |
|
|
|
|
|
E1 −E3 |
|
|
|
E1 −E4 |
|
29
|
C2µC2νC3µC3ν |
+ |
C2µC2νC4µC4ν |
) = 4( |
C11C12C31C32 |
+ |
|||||||
|
|
|
|
|
|
||||||||
|
E2 −E3 |
|
|
E2 −E4 |
|
|
|
E1 |
−E3 |
||||
C11C12C41C42 |
+ |
C21C22C31C32 |
+ |
C21C22C41C42 |
) = |
|
|
||||||
|
|
|
|
|
|||||||||
E1−E4 |
|
E2 −E3 |
E2 |
−E4 |
|
|
= 4[(0,37 • 0,6 • 0,6 • (-0,37)) / 2,24β + (0,37 • (-0,37) • 0,6 • (-0,6)) / 3,24β + + (0,6 • 0,6 • 0,37 • (-0,37)) / 1,24β + (0,6 • 0,37 • 0,37 • (-0,6)) / 2,24β] = -
0,402 / β
π1,3. = 4 • 0,05(-2 / 2,24 + 1 / 3,24 + 1 / 1,24) / β = 4 • 0,05(-0,893 + 0,806 + 0,309) / β = 0,044 / β
π1,4. = 4(0,05 • 2 / 2,24 − 0,13 / 1,24 − 0,019 / 3,24) / β = 4(-0,066) = 0,264 /
β
π1,1. = 4[((0,37)2 • (0,6)2 ) / 2,24 + ((0,37)2 • (0,37)2 ) / 3,24 + + ((0,6)2 • (0,6)2 ) / 1,24 +((0,6)2 • (0,37)2 ) / 2,24] / β = 4(0,045 + 0,006 +
0,105) = 4(0,156) = 0,624 / β
Выпишем полученные результаты.
π1,1. = 0,624 / β π1,2. = -0,402 / β π1,3. = 0,044 / β π1,4. = 0,264 / β
1.7.2 Расчет индуктивного π эффекта в молекуле пропеналя
СН2СНСНО
Замещение одного из концевых атомов углерода в бутадиене кислородом приводит к возмущению системы, которое проявляется как индуктивный эффект в пропенале.
Изображая только углеродные центры, возмущенную молекулу можно изобразить схематически следующим образом.
О*1 •________• С2 ________ С*3 •________• С4
Электронная плотность на каждом атоме с учетом возмущающего действия атома кислорода определяется следующим образом:
qµ = 1+ πµ,х.(αх − αс).
30
Поскольку .(αО − αС) = β. Отсюда .αО = αС + β)
1.8 Практические задания для самостоятельной работы студентов
На основе теоретического материала и приведенных примеров рассчитать методом Хюккеля структуру МО, электронные, спектральные характеристики и описать химические свойства и реакционную способность следующих систем:
-циклопропан, анион (катион);
-метиленциклопропен, анион (катион);
-циклобутан, анион (катион);
-циклопентан, анион (катион);
-циклогексан, анион (катион);
-гексатриен, анион (катион);
-нафталин, анион (катион);
-антрацен, анион (катион);
-бензильный радикал, анион (катион);
-фенол, анион (катион);
-анилин, анион (катион);
-фульвен, анион (катион);
-пиррол, анион (катион);
-пиридин, анион (катион);
31
2 Основная рекомендуемая литература
1.Мелешина А.М. Курс квантовой механики для химиков. – М.: Высшая школа, 1980.-215с.
2.Голованов И. Б., Пискунов А.К., Сергеев Н.М. Элементарное введение в
квантовую биохимию. – М.: Наука, 1969.-235с.
3.Слетер Дж. Электронная структура молекул. – М.: Мир, 1965.-345с.
4.Л. Цюлике. Квантовая химия: Основы и общие методы. - Т.1.- М.: Мир, 1976.-512с.
5.Фларри Р. Квантовая химия. – М.: Мир, 1985.-472с.
6.Минкин В. И., Симкин Б.Я., Миняев Р.М.. Теория строения молекул: Электронные оболочки. – М.: Высшая школа, 1979.-408с.
7.Жидомиров Г. М, Багатурьянц А.А., Абронин Н.Д. Прикладная квантовая химия: Расчеты реакционной способности и механизмов химических реакций. – М.: Химия, 1979.-292с.
8.Флайгер У. Строение и динамика молекул. - Т.1.- М.: Мир, 1982.- 407с.
9.Клопман Г., Иванс Р. Методы пренебрежения дифференциальным перекрыванием. Полуэмпирические методы расчета электронных
структур. - Т.2.- М.: Мир, 1980.-437с.
10.Эйринг Г., Уолтер Дж., Кимбалл Дж. Квантовая химия. - М.: ИЛ, 1948, 597с. .
3 Дополнительная литература
11.Базилевский М.Б. Метод молекулярных орбиталей и реакционная способность органических молекул. – М.: Химия, 1969.-303с.
Список использованных источников
12.Минаев Б.Ф., Минаева В.А., Лыженкова И.И. Введение в квантовую химию. – Черкассы.: ЧХТИ, 1998.-161с.
13.Дьюар М. Теория молекулярных орбиталей в органической химии. - М.:
Мир, 1967.-438с.
14.Хигаси К., Баба Х., Рембаум А. Квантовая органическая химия. – М.:
Мир, 1967.-195с.
15.Стрейтвизер Дж. Теория молекулярных орбиталей для химиков органиков. – М.: Мир, 1965.-384с.
16.Маррел Дж., Кеттл С., Теддер Дж. Химическая связь. – М.: Мир, 1980.- 382с.
17.Заградник Р., Полак Р. Квантовая химия. – М.: Мир, 1979.-504с. 18.П. Эткинс. Кванты: Справочник концепций. – М.: Мир, 1977.-496с.
19.Сигал Дж. Полуэмпирические метода расчета электронной структуры. -
Т.2.- М.: Мир, 1980. 327с.
32