Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Расчет электронных характеристик молекул полуэмпирическим методом Хюккеля

.pdf
Скачиваний:
135
Добавлен:
08.01.2014
Размер:
879.57 Кб
Скачать

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

) =

 

 

 

 

 

 

 

E1E4

 

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

Соседние файлы в предмете Химия