Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лекции1.pdf
Скачиваний:
210
Добавлен:
30.05.2015
Размер:
7.13 Mб
Скачать

ЛЕКЦИЯ 13

МОДЕЛИ МОЛЕКУЛЯРНЫХ СИСТЕМ

Исследование многочисленных литературных источников показывает, что познание химических особенностей молекул неразрывно связано с познанием их природы. Вопрос о причинах вступления в связь одинаковых, или различных, атомов так, или иначе вставал при анализе результатов химических реакций, но сложность этого вопроса породило много мнений, степень доказанности которых почти всегда конкурировала с простотой и кажущейся близости к истине. Когда эксперимент может быть объяснен с помощью двух различных моделей, появление новой (третьей) модели почти неизбежно. Квантовомеханическая модель атома сразу привлекла к себе внимание ввиду успехов объяснения свойств периодической таблицы Менделеева, и было ясно, что свойства молекул подчиняются тем же законам, что и свойства атомов.

Самой очевидной моделью, появившейся среди химиков еще до квантовой механики, была модель валентных связей (ВС), основанная на примитивном понимании притяжения разноименных, и отталкивания одноименных (по знаку заряда) заряженных частиц. (Такое понимание распространено даже сейчас в смысле взаимодействия отрицательно и положительно заряженных атомов, тогда как одного взгляда на периодическую таблицу достаточно, чтобы убедиться в том, что все атомы электронейтральны). Видимо, более общим названием целой серии методов построения моделей молекул является модель конфигурационного взаимодействия (КВ). Метод молекулярных орбиталей (МО) стал некоторым перекрестком развития этой модели, а также наращиванием собственных особенностей, присущих только квантовохимическому подходу. То, с чем никак не мог справиться метод валентных связей – объяснением двухатомных молекул (из двух одинаковых атомов) логично решил метод молекулярных орбиталей. Однако, этот метод внес и целый ряд новых идей, быстро превратившихся в новые модели. В этом ряду моделей обращает на себя внимание разнообразие типов связей, приведшее к классификации их характеристик – ионные и ковалентные. Это привело к аналогичным результатам в классификации типов молекул, в основу которых были положены типы связей, и из моделей таких молекул родились новые типы связей, которые было трудно отнести к ионным, или ковалентным. Здесь мы ведем речь о «координационных» связях; о молекулах, чья чисто молекулярная природа до сих пор под сомнением – клатратах, ассоциатах, коллоидных частицах, и т.д.; и наконец, о связях частиц в растворах (которые часто

125

характеризуют «просто растворами»). Вряд ли можно критиковать всю эту гамму подходов, применяя ядовитые характеристики и определения. Однако, лишь знания помогают более близко к истине («более правильно») оценить собственную точку зрения.

МОЛЕКУЛА H2 C ПОЗИЦИЙ КВАНТОВОЙ МЕХАНИКИ.

Квантовая химия берет начало с работы немецких ученых В.Гайтлера и Ф.Лондона, опубликованной в 1927 году под названием Взаимодействие нейтральных атомов и гомеополярная связь, которая начинается словами: Взаимодействие между нейтральными атомами до сих пор представляло большие трудности для теоретического рассмотрения. Развитие квантовой механики дало для разработки этой проблемы совершенно новую точку зрения: прежде всего в новой моделираспределение заряда полностью отлично от модели Бора, что уже влечет за собой совершенно новое соотношение сил (Kraftenspiel - буквальный перевод «игру сил») между нейтральными атомами. Этой работой начинается цикл работ под общим названием «метод молекулярных орбиталей как линейной комбинации атомных орбиталей (ЛКАО – МО)».

Не останавливаясь на расчетной части этой статьи, запишем допущения, упрощающие задачу, решенную в терминах квантовой механики тех времен.

1.Электрон не может быть точно локализован на орбитали ( согласно принципу неопределенности), и можно рассчитать только вероятность того, что электрон будет находится в данном элементе объема около ядра.

2.Электрону на орбитали может быть приписана определенная энергия.

3.Орбиталь могут занимать только два электрона.

4.Силы, обуславливающие связи пока принимаем только электрическими.

Рассмотрение начнем с молекулярного иона водорода H 2+ , в

котором реализуется простейшая из всех видов химической связи. Итак, полная энергия представляется суммой потенциальной и кинетической энергии: E = P + K . Если система под-

чиняется

 

законам квантовой механики, то

можно написать, что

ˆ

,

что отражает волновое движение,

записанное через га-

HΨ = EΨ

мильтониан. Существует квадрат волновой функции

(Ψ2 ( x, y, z)dxdydz) ,

которая пропорциональна вероятности нахождения электрона в элементе объема dv = dxdydz декартовых координат.

126

Далее, если

+∞+∞+∞Ψ2dxdydz = 1

(или Ψ2dτ = 1),

то волновая

 

−∞−∞−∞

 

 

функция нормирована. С учетом

комплексных

переменных

ΨΨ dτ = 1.

Каждому состоянию атома водорода 1s, 2s, 2p и т.д. можно приписать соответствующую волновую функцию, и с ее помощью рассчитать плотность вероятности и энергию электрона.

Ясно, что для наших расчетов гамильтониан не зависит от времени.

Свойства H 2+ могут быть рассчитаны в методе молекулярных орбиталей, если представить молекулу в виде двух ядер, охваченных одной молекулярной орбиталью с одним электроном, представленной Ψ - функцией молекулы:

 

ˆ

2

если Ψмол

нормирован.

 

HΨмол = EΨмол,

Ψмолdτ = 1

 

Сами по себе эти выражения не громоздки, трудности возни-

кают при выборе форм гамильтониана и Ψмол

и при использовании

ˆ

иΨ для расчета энергии E.

 

 

H

метода

ЛКАО предполагает,

 

Применение к

расчету H2+

чтоΨмол может быть приближенно представлена в виде линейной

комбинации атомных орбиталей, имеющих свои волновые функции χn . Итак,

Ψмол c1χ1 + c2χ2 .

Можно предположить, что для H2+ коэффициенты c1 и c2 равны, но не будут равными для несимметричной молекулы, скажем LiH . В дальнейшем мы выясним, что число получаемых таким образом молекулярных орбиталей равно в методе ЛКАО числу используемых атомных орбиталей.

Будем считать c1 и c2 параметрами, значения которых мы хотим определить, символ Ψ будем использовать для обозначения волновой функции молекулы, χ1 и χ2 - для обозначения соответствующих атомных орбиталей. E определим в терминах и c1 и c2 и энергий атомных орбиталей. Домик на обозначениях гамильтониана будем упускать, не забывая однако , что это оператор Гамильтона

ˆ

 

 

 

 

H со всеми своими свойствами.

 

 

 

Итак, HΨ = EΨ ,

умножив обе части равенства на Ψ получим

выражение

 

 

 

 

ΨHΨ = EΨ2 .

 

 

 

Интегрирование

по

всему

пространству

да-

ет: ΨHΨdτ = EΨ2dτ ,

или

 

 

 

127

Ψ2dτ

E = ΨHΨdτ .

В последнем уравнении выражение для E получается в такой форме, что проблема координат сильно упрощается. Подставляя Ψ = c1χ1 + c2χ2 , имеем

E = (c1χ1 + c2χ2 )H(c1χ1 + c2χ2 )dτ = (c1χ1 + c2 χ2 )2 dτ

(c1χH1c1χ1 + c1χ1Hc2χ2 + c2χ2 Hc1χ1 + c2 χ2 Hc2χ2 )dτ . (c12χ12 + 2c1c2 χ1χ2 + c22 χ22 )dτ

Можно показать, что для тех значений E , которые имеют физический смысл,

χ1Hχ2dτ = χ2 Hχ1 dτ .

Теперь можно сделать подстановки: H11 = χ1 Hχ1dτ ,

H12 = H21 = χ1Hχ2dτ = χ2 Hχ1dτ , H22 = χ2 Hχ2dτ ,

 

 

S11 = χ12dτ ,

 

S12 = χ1χ2dτ , S22 = χ22dτ .

Следовательно:

 

 

 

 

 

 

 

E =

c2 H

11

+ 2c c

H

12

+ c2 H

22

.

1

 

1 2

 

 

 

2

 

c2 S

 

+ 2c c

S

 

+ S

22

 

 

 

 

 

 

 

1

11

1

2

 

12

 

 

 

Нас интересует минимальное значение энергии. Для этого можно последнее выражение подвергнуть частному дифференциированию по с (отдельно по c1 и c2 и приравнять полученные частные дифференциалы к нулю.

dE

=

(c2 S

11

+ 2c c

S

12

+ c2 S

22

)(2c H

11

+ 2с

H

12

)

 

1

1

2

 

 

2

 

 

1

2

 

 

dc

 

 

(c2 S

 

+ 2c c

S

 

+ c2 S

22

)2

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

11

 

1

2

 

12

2

 

 

 

 

 

(c2 H

11

+ 2c c

 

H

12

+ c2 H

22

)(2c S

 

+ 2с

S

 

)

= 0 .

 

 

 

 

 

 

 

 

 

1

1

2

 

 

 

2

 

 

 

1

11

2

12

 

 

 

 

 

 

 

 

 

 

 

 

(c2 S + 2c c

S + c2 S

22

)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

11

 

 

1

2

12

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2c H

 

 

+ 2c H

 

) =

(c2 H

11

+ 2c c

H

12

+ c2 H

22

)(2c S

 

+ 2с S

)

.

 

 

 

 

 

1

 

1

2

 

 

 

2

 

1 11

2 12

 

 

 

 

 

 

 

 

 

(c2 S

 

+ 2c c S

 

 

 

 

 

 

1

11

 

 

1

 

12

 

 

 

 

 

 

 

+ c2 S

22

)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

11

 

 

1

2

12

2

 

 

 

 

Следовательно: (c1H11 + c1H12 ) = E(c1S11 + c2 S12 )

или c1(H11 ES11 ) + c2 (H12 ES12 ) = 0 .

Таким же образом при dE = 0 dc2

получаем: c1(H12 ES12 ) + c2 (H22 ES22 ) = 0 .

Разрешенные значения E для системы совместных линейных однородных уравнений соответствуют корням векового детерминанта

128

H11-ES11 H12-ES12

= 0

H12-ES12 H22-ES22

Найдя значение E мы можем определить отношения c1 и c2 из системы совместных линейных уравнений. Окончательное значение коэффициентов c1 и c2 должны быть найдены из условия нормировки.

В общем случае, когда Ψ = с1Ψ1 + с2Ψ2 + ... + сnΨn и вековой детерминант принимает вид:

H11

S11 E

H12

S12 E

...

H1n S1n E

 

H12

S12 E

H22

S22 E

...

H2n S2n E

= 0

H13 S13 E

H23 S23 E

...

H3n S3n E

 

...

 

...

...

 

H1n S1n E

H2n S2n E

...

Hnn Snn E

 

 

 

 

 

 

 

 

Такие детерминанты симметричны относительно диагонали (эрмитовы) и имеют n действительных корней.

Для решения детерминантов необходимо теперь вычислить значения величин Hij и Sij

129