- •Лекция 5. Начала компьютерной химии. Математические основы теории мо лкао (Простейшая схема одноэлектронной теории).
- •3. Приближение лкао мо.
- •4. Базис и его свойства.
- •5.1. Эффективный гамидьтониан.
- •6. Процедура расчёта мо лкао
- •6.1. Одноэлектронный принцип (одноэлектронное приближение).
- •6.2. Принцип минимума энергии.
- •6.3. Принцип Паули.
- •6.4. Правило Хунда.
6. Процедура расчёта мо лкао
Молекула состоит из атомов,
которые перечислим, нумеруя в любом порядке, но далее всюду соблюдая принятую нумерацию:
Количество АО у каждого свободного атома не ограничено, и, конечно, наилучшим базисом был бы такой, в котором содержалось максимально возможное количество функций, однако не существует возможности использовать базис бесконечной длины в практических расчётах. Поэтому во всех, даже наиболее совершенных, расчётах используется ограниченное количество базисных функций АО. При выборе базиса используются различные способы ограничения его длины. В более совершенных вариантах теоретических расчётов МО ЛКАО базисные наборы достаточно велики. В более простых приближениях применяют более ограниченные базисные наборы АО.
Наиболее простым и потому распространённым является валентное приближение. Оно состоит в том, для расчёта электронного строения молекулы используют базис, включающий от каждого атома только АО внешнего оптического слоя. Главное квантовое число всех этих АО одно и то же, но при необходимости базис можно расширить за счёт включения и АО близлежащих подуровней с соседними значениями главного квантового числа. Рассчитываемая часть спектра МО заселяется оптическими электронами, находящимися на внешних валентных АО атомов.
У каждого атома в нейтральном зарядовом состоянии число электронов равно его порядковому номеру в Системе Менделеева. В начале расчёта атомы нумеруются. Все внутренние АО, уровни которых лежат ниже оптических АО, считаются атомными, т.е. их вид считается таким же, как и в свободном атоме, и это недалеко от истины. Набор АО последовательно нумеруется.
Согласно математической схеме квантовой механики по Дираку, всякая волновая функция представляется двумя неразрывно связанными образами-«компонентами»: бра- икет-векторами.... Это в чём-то напоминает ситуацию с комплексными числами, каждое из которых само по себе ничего не представляет, обретая смысл лишь в паре со своим комплексно сопряжённым двойником. Похожим образом обстоит дело и сбра-кет- векторами Дирака. Каждый из них неотделим от своего двойника, и лишь вместе они служат основой очень элегантного и компактного построения квантовой механики. Это выглядит следующим образом. Каждой базисной АО отвечают два вектора: бра-вектор и кет-вектор АО, а именно:
«Бра-кет» - символы это не что иное, как «ско- бки»- символы (от англ. brackets-скобки).
Каждый из этих векторов МО представляется в виде разложения в ряд и при этом:
- бра-векторы МО построены только из бра-векторов АО.
- кет-векторы МО построены только из кет-векторов АО.
В наших расчётах приходится часто иметь дело с массивами однотипных объектов. Условимся опускать многократно повторяющиеся символы. Так, производя вычисления, нет смысла повторять всё время один и тот повторяющийся символ в последовательности базисных АО:. С таким же успехом можно ограничиться просто перечислением индексов:1,2,3,...p,...q,..., в нужном контексте идентифицируя с ними сами АО. В итоге получается предельно компактная запись:
АО и их обозначения в виде индексов |
Бракет-символы АО |
Пределы изменения индексов АО |
МО |
Бра-кет- символы МО |
Разложения МО ЛКАО |
; ;
Точно так же, как введены индексные обозначения для базисных АО, ничто
не мешает при необходимости ввести индексные обозначения и для МО...
Свойства МО и составляющих их БАО следующие:-
Введём понятие молекулярного эффективного одноэлектронного гамильтониана.
Молекулярные орбитали являются его собственными функциями.
3 )или эквивалентно, т.е.
Здесь представлено скалярное произведение двух векторов-массивов вида
Математическая справка:
1) Вектором называется упорядоченный массив : .
Вектор-массив, у которого все элементы равны нулю, называется нулевым, например:
2) Из двух массивов равной размерности можно образовать сумму парных произведений элементов с одина-
ковыми индексами. Эта сумма называетсяскалярным произведением массивов*)
Пример:
*) В общем случае читатель может встретить в литературе и другие виды скалярных произведений,
и их структура зависит от природы математических объектов, для которых они введены.
3) Если скалярное произведение равно нулю, то образующие его векторы называютсяортогональными.
Равенство может быть удовлетворено в 2-х случаях:
а) Оба вектора - сомножителя могут оказаться равны нулю: (;)=(0;0)
Для решения систем уравнений этот тривиальный случай не представляет интереса.
б) Один из векторов равен нулю: =0 или =0.
Для решения систем уравнений важен именно этот случай.
4) Таким образом уравнение3)видадолжно быть удовлетворено при условии:
.
Случай А)тривиальный и в качестве решения системы уравнений исключается, поскольку нулевыми быть собственные векторы МО не могут, и потому остаётся
Случай Б)как единственный вариант решения системы уравнений. Это означает, что
hq=0 q(1,2,...n), т.е..
В этой формуле каждому значению индекса qотвечает полный набор возможных значений индекса p, и все
они представлены в одной строке - сумме. Оба индекса изменяются в одних и тех же пределах. Индекс qпробегает все возможные значения от1 до n, и получаем систему изn однотипных линейных уравнений с нулевой правой частью. Согласно теореме Крамера, её нетривиальное решение возможно, если вековой детерминант равен нулю. Раскрытый детерминант представляет из себя алгебраическое уравнение, степень которого равна длине выбранного базиса. Это даёт:
Это алгебраическое уравнение называется вековым уравнениемрешаемой задачи, а в виде ещё нераскрытого детерминанта -вековым детерминантом.
Решение векового уравнения приводит к nкорням, и каждый из них это уровень энергии МО в исследуемой части энергетического спектра МО. Расширяя базис, можно и увеличить число уровней, и отчасти расширить границы исследуемой части спектра энергетических уровней МО. В результате расчёта получаем характеристики, которые можно представить в виде следующей сводки.
Сводка результатов расчёта МО ЛКАО
5)Первая строка полученной сводки результатов это полученные энергетические уровни МО, перечисленные, как правило, в порядке увеличения уровня энергии:. Индексi это номер уровня молекулярной орбитали (МО).
Базисные АО нумеруют в определённой выбранной последовательности, например,
. Коэффициенты при этих АО в составе каждой МО будут перечислены в той же последовательности, образуя вектор . Это один из вполне однозначных способов изображения МО.
Такой вектор называется собственным вектором МО. Иначе собственный вектор называетсяатомно-орбитальным составом МО.
Понятие собственных векторов имеет смысл только в том случае, если определены все базисные единицы разложения, из которых эти векторы состоят.
Располагая собственные векторы в определённом порядке, получаем квадратную матрицу с элементами {cip}, называемуюматрицей составов МО. Если первый индексi это номер уровня МО, то второй индексp или q это уже номер базисной АО. В нашей сводке собственные векторы это столбцы таблицы составов.
Полученные уровни энергии {E 1,E 2,E 3,...E n} в порядке их возрастания заполняют электронами, следуя правилам построения электронной конфигурации.
Они следующие: