Квантовая химия — ее прошлое и настоящее. Развитие электронных представлений о природе химической связи - [41]

Погрешности, вносимые в одноэлектронную матрицу плотности молекулы циклопропана локализацией МО, характеризуются величинами:

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

Анализ заселенности атомных орбиталей

[23]

В квантовой химии функции Ψ, описывающие состояния многоэлектронных систем, обычно выражаются через АО или подобные им более элементарные функции, центрированные на различных ядрах молекулы. Такой выбор одноэлектронного базиса {φ} основан на предположении, что атом, представленный определенным набором орбиталей в таком базисе, сохраняет в молекуле свою индивидуальность. В связи с этим встает задача построения формализма, позволяющего характеризовать состояния химически связанных атомов и установить корреляции между этими состояниями и состояниями невзаимодействующих атомов.

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

атомов. Однако анализ заселенностей осложняется перекрыванием атомных орбиталей в молекуле.

Если последние ортогональны, т. е. не перекрываются, то их заселенности определяются однозначно. При этом заселенность па имеет смысл вероятности нахождения электрона в состоянии, заданном атомной орбиталью φ_а, и может быть выражена как математическое ожидание:

(4.58)

Если атомные орбитали неортогональны, то положение осложняется. Понятие заселенности отдельной АО становится неоднозначным и распадается на дополняющие друг друга понятия полной, неподеленной и аддитивной заселенности и заселенности перекрывания, которые связаны с различными способами ортонормировки атомно-орбитального базиса {φ}.

Полные заселенности (n⁺_а) орбиталей φ_а неортогонального базиса {φ} определяются аналогично заселенностям ортогонального базиса:

(4.59)

Предполагая, что оператор электронной плотности ρ представлен в базисе {φ} матрицей

, определение (4.59) можно переписать также в виде

(4.60)

Детальное исследование заселенностей n⁺_a было проведено Дэвидсоном [37] и Роби [74], которые показали, в частности, что

где n₁ — наибольшая из естественных заселенностей. Это неравенство, как и аналогичные неравенства для определяемых ниже заселенностей n^-_a и n⁰_а, следует из условия антисимметрии многоэлектронной функции Ψ(х₁,..., х_N) относительно перестановок электронных координат.

Неподеленную заселенность (n^-_а) орбитали φ_а можно определить как заселенность ее компоненты, которая ортогональна ко всем прочим орбиталям φ:

(4.62)

где

(4.63)

^aS в формуле (4.63) обозначает матрицу перекрывания, полученную из полной матрицы S вычеркиванием интегралов перекрывания S_ab, включающих рассматриваемую орбиталь φ_а.Такая ортогонализация (аналогичная ортогонализации по методу Шмидта) исключает из полной заселенности n⁺_а ту ее часть, которая принадлежит не только φ_а, но и остальным орбиталям неортогонального базиса (рис. 23).

Рис. 23. Геометрическая иллюстрация к определению неподеленной электронной заселенности

Учитывая отмеченное Галлупом и Норбеком [40] равенство

(4.64)

выражение(4.62) можно привести к чрезвычайно простому виду

(4.65)

В частном случае одноэлектронной системы, состояние которой описывается орбиталью

(4.66)

диагональные элементы матрицы плотности равны

(4.67)

(4.68)

Эта формула, то чиее ее правая часть, приводилась в работе [40], но лишь в качестве промежуточного результата. Окончательное выражение для заселенностей (по Галлупу и Норбеку) получалось путем нормирования n^-_а на единицу:

(4.69)

Обобщение формулы (4.69) на многоэлектронные системы, очевидно, должно осуществляться заменой |С_а|² на Р_аа:

(4.70)

Однако такой подход к проблеме является ошибочным. Расчеты свидетельствуют, в частности, о чрезмерно больших значениях n(GN) для АО внутренних оболочек и неподеленных электронных пар. Например, в молекуле LiH:

Заселенность перекрывания орбитали φ_а с остальными орбиталями неортогонального базиса φ определяется как разность между полной и неподеленной заселенностями:

(4.71)

Заселенность перекрывания

представляет ту долю полной электронной заселенности, которая принадлежит одновременно к рассматриваемой и всем прочим базисным АО. Нетрудно убедиться в том, что величина

равна нулю, если АО φ_а не перекрывается ни с одной из орбиталей базиса φ, т. е. если

(4.72)

для всех b≠a.

Аддитивная заселенность.

Сумма засел енностей n⁺_a или n^-_а по всем базисным орбиталям совпадает с числом электронов (N) в рассматриваемой системе только в том случае, если эти орбитали ортогональны. Иными словами, заселенности орбиталей неортогонального базиса неаддитивны.