Приближение МО ЛКАО. Уравнения Рутана.

Уравнения Хартри-Фока для молекул можно, в принципе, решить численно, получив МО в виде таблиц. Извлечение, однако, химической информации из МО, представленных таким образом - нелегкая задача. Существует простое и мощное приближение, позволяющее существенно упростить как решение уравнений ХФ, так и интерпретацию результатов. Суть его в следующем.

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

Коэффициенты разложения с_im определяют “вес” каждой АО в МО; ясно, что одна и та же АО входит в разные МО с разными весами, т.е. эти коэффициенты различны для каждой МО.

Представление (2.12) очень удобно для применения вариационного метода. Применяя его с учетом ортонормированности МО, из условия минимума энергии получают уравнения ХФ в виде:

Элементы матрицы Фока в приближении МО ЛКАО имеют вид:

где

а S_ij –матрица интегралов перекрывания между АО c_m и c_n , а Е_m - одно из решений секулярного уравнения

Уравнения (2.13) называются уравнениями Рутана.

Уравнения Рутана можно переписать в матричном виде:

F C = S C E ,(2.19)

который путем унитарных преобразований F⁺=S^-1/2FS^-1/2 и C⁺=S^-1/2C сводится к стандарнной задаче на собственные значения.

F⁺C⁺= EC⁺. (2.20)

Введем теперь матрицу P c элементами

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

Энергия молекулы с закрытыми оболочками в терминах введенных обозначений записывается следующим образом:

Решение уравнений Рутана осуществляется таким же итерационным методом ССП, как и в случае атома. Блок-схема вычислительного процесса представлена на рис. 2.1.

Двухэлектронные кулоновский и обменный интегралы по АО (или по базисным функциям - см. ниже) являются причиной большинства практических проблем в случае, когда расчет проводится из первых принципов. Число их без учета симметрии равно » N⁴/8, где N - число АО, т.е. очень велико (для сравнения: число одноэлектронных интегралов hm n равно » N²/2 – см. табл. 2.2). Таким образом, вычисление этих интегралов и их запоминание в памяти компьютера представляет собой отдельную проблему, особенно, если учесть, что все четыре орбитали c_m ,c_n , c_l и c_s в общем случае центрированы на разных ядрах. Отсюда ясно, сколь важен для эффективного вычислительного процесса выбор аналитического вида функций, по которым вычисляются интегралы. В так называемых прямых ССП методах (доступных при наличии суперкомпьютеров), эти интегралы не запоминаются, а вычисляются в каждой итерации. При этом малыми по величине интегралами часто пренебрегают, что увеличивает риск ошибки, ибо число их велико, а знаки одинаковы (положительны).

Таблица 2.2 Общее число одноэлектронных и двухэлектронных интегралов в зависимости от базиса

Волновая функция и энергии, получаемые с помощью метода Рутана
(и с помощью метода Хартри-Фока вообще), инвариантны относительно ортогонального преобразования занятых электронами спин-орбиталей. Это означает, в частности, что если вместо АО для построения МО (2.12) будут использованы их линейные комбинации, полученные с помощью ортогональных преобразований, то одновременно и согласовано изменятся как все одно- и двухэлектронные интегралы, так и матрица P. Волновая функция (2.10) и энергия системы (2.22) при этом останутся прежними. Это очень важное свойство метода: используя его всегда можно облегчить концептуальную химическую трактовку результатов, например, перейти к локализованным МО, описывающим электроны связи и неподеленные электронные пары, гибридизованным АО, и т.д.

Ортогональным преобразованием называется преобразование, которое не меняет длин векторов и углов между ними и переводит один ортонормированный базис в другой ортонормированный. Матрица ортогонального преобразования Т обладает свойствами: ТТ' = 1, det T =± 1.

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

Приведем несколько примеров исследования свойств молекул неэмпирическим методом Хартри-Фока.

Таблица 2.3. Длины связи (А) и углы (град.) для циклофосфорамида

Таблица 2.4. Энергии (ккал/моль)вращательные барьеров

После символа // базис, с которым была оптимизирована геометрия.

Таблица 2.5. Энергии (ккал/моль) реакций изомеризации