2.2.3. Общая формулировка метода мо

Описание молекулы водорода в рамках метода МО является относительно простой задачей, вследствие малого числа ядер и высокой симметрии объекта. Для более сложных молекул задача существенно осложняется необходимостью вычисления коэффициентов МО вариационным методом. Процедура выглядит следующим образом:

1) выражают каждую МО в виде ЛКАО (в общем виде);

2) строят глобальную волновую функцию в виде определителя Слэтера из МСО;

3) выражают полную энергию молекулы в виде функции от коэффициентов МО: Е=Е(С_ij);

4) дифференцируют эту функцию по всем коэффициентам и приравнивают производные к нулю: Е/С_ij= 0.

Получаемая система уравнений (содержащая nnуравнений), называетсяуравнениями Хартри-Фока-Рутана(ХФР). Она состоит изnодинаковых экземпляров системы, содержащейnуравнений следующего вида:

где F_— матричные элементыоператора Фока, характеризующие энергии взаимодействия пар атомов с номерамии,S_ — интегралы перекрывания для базисных АО с номерамии, а— энергия МО с коэффициентами {С_С_. .С_n }.

Метод решения ХФР-уравнений сводится к следующему. Всякая система однородных линейных уравнений имеет решение только в том случае, если ее определитель равен нулю. Приравнивание определителя системы ХФР к нулю дает "характеристическое уравнение" n-й степени относительно, имеющееnкорней:₁,₂, … ,_n. Подставляя поочередно эти корни в систему, находятnрешений, каждое из которых представляет собой набор коэффициентов для одной из МО: (С₁,С₂, …С_n)₁, (С₁,С₂, …С_n)₂и т.д.

Вычислительные проблемы метода МО связаны с тем, что значения интегралов типа FиSзаранее неизвестны (они зависят от коэффициентов МО, которые являются решением уравнений ХФР). Поэтому приходится прибегать к трудоемкой процедуре самосогласования, аналогичной той, которая используется в методе Хартри-Фока в теории многоэлектронных атомов. По окончании процедуры (после достижения заданной точности) получают самосогласованное решение в виде набора МО и их энергий:

Такой подход к решению системы ХФР носит называние неэмпирического метода (или "метода ab initio") и отличается тем, что в нем все интегралыFиSвычисляются в ходе итерационной процедуры. Известны и альтернативные варианты, которые называютсяполуэмпирическими, так как в них часть интеграловFиS(или даже все из них) находятся из эмпирических данных, например, спектральных или калориметрических. Это позволяет значительно снизить число итераций и ускорить процедуру решения. Полуэмпирических методов существует чрезвычайно много, причем все они дают несколько различные результаты. Самый простой из них (метод Хюккеля) будет подробно рассмотрен ниже.

Следует отметить некоторые важные особенности результатов, получаемых методом МО. Оператор Гамильтона в методе МО строят в виде:

H =  (_i) +  (U_i) +  (U_ij) +  (U_) =  (h_i) +  (U_ij)

где _i= (–²/2m)²_i— оператор кинетической энергииi-го электрона,

U_i= –Z_e²/r_i— оператор потенциальной энергииi-го электрона в кулоновском поле-го ядра,

U_ij=e²/r_ij— операторы межэлектронного отталкивания,

h_i =_i+(U_i) — т.н. "одноэлектронные гамильтонианы".

Поэтому полная энергия молекулы в методе МО задается в виде суммы орбитальных энергий с поправками на межэлектронное взаимодействие:

Е=_i^* + (J_ijK_ij)

где J_ij—кулоновскиеинтегралы, аK_ij—обменныеинтегралы (их физический смысл тот же, что и в теории МЭА).

В большинстве случаев наилучших результатов удается достичь, если учитывать конфигурационное взаимодействие (метод МО-КВ), когда используются многодетерминантные волновые функции, получаемые при "смешивании" нескольких электронных конфигураций. Следует иметь в виду, что КВ — это не особый вид физического взаимодействия, а лишь способ введения поправок в волновую функцию молекулы и ее энергию.