19. Теорема Купмэнса.

Энергия орбитали соответствует потенциалу ионизации при удалении электрона с данной орбитали.

КУПМАНСА ТЕОРЕМА, позволяет рассчитывать потенциал ионизации атома или молекулы в рамках метода молекулярных орбиталей. Состояние электрона в молекуле описывается уравнением Шрёдингера, которое, определяет волновую функцию (орбиталь) и энергию электрона (орбитальную энергию). Согласно Кумпанса теореме, потенциал ионизации равен по абс. величине орбитальной энергии электрона, удаленного из молекулы; система орбиталей в целом при этом приближенно сохраняется. Разность энергий молекулы и иона, полученного удалением одного электрона при сохранении положений ядер, наз. первым вертикальным потенциалом ионизации. Образованию иона в разных электронно-возбужденных состояниях молекулы отвечают разные первые потенциалы ионизации; самый малый из них обычно наз. просто потенциалом ионизации. Кумпанса теорема позволяет также оценивать сродство к электрону по орбитальной энергии низшей незанятой (виртуальной) орбитали. Кумпанса теорема не учитывает согласованности движения электронов в молекуле (электронной корреляции) и перестройки орбиталей иона по сравнению с орбиталями молекулы. Учет этих эффектов возможен методами возмущений теории, однако обычно такие эффекты в значит. степени взаимно компенсируют друг друга. Это позволяет с успехом применять Кумпанса теорему для интерпретации рентгене- и фотоэлектронных спектров молекул, а также при изучении Ожеспектров.

Выбор базисного набора,уравнения Рутаана – Холла

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

1) функции должны вести себя так, чтобы соответствовать физической модели, чтобы при увеличении базисного набора итерационную процедуру самосогласования можно было производить относительно быстро; т.е. функции должны стремиться к нулю при увеличении расстояния между ядрами и электронами до бесконечности;

2) выбранные функции должны позволять легко вычислять все необходимые интегралы.

Первый критерий предполагает предпочтительное использование экспоненциальных функций, локализованных на ядрах, например, функций водородоподобного атома или функций Слэйтера. Но, к сожалению, расчёт интегралов этих функций достаточно трудоёмок. Гауссовы функции достаточно просто позволяют вычислять интегралы, но они достаточно плохо описывают электронную структуру. Однако простота вычисления интегралов в данном случае имеет большее значение.

Гауссовы функции отличаются от функций Слэйтера в своей радиальной части (ξ и α – экспоненциальные множители):

– радиальная часть функции Слэйтера

– радиальная часть функции Гаусса .

Каждую молекулярную орбиталь можно представить в виде разложения по базисным орбиталям. Если в качестве базисных функций взять атомные орбитали (M функций в (3.3.1)), то молекулярные орбитали будут являться линейными комбинациями атомных орбиталей (МО = ЛКАО):

Уравнения Хартри – Фока при этом приобретают следующий вид:

Умножение слева на базисную функцию и интегрирование даёт уравнения Рутаана – Холла (3.3.3) (для систем с закрытой оболочкой). Фактически это уравнения Фока в базисе атомных орбиталей:

FC = SCε , (3.3.3)

где

Матрица S содержит элементы перекрывания между атомными орбиталями, матрица F – элементы матрицы Фока. Каждый матричный элемент F_ij (3.3.4) содержит две чисти оператора Фока: 1) интегралы, соответствующие одноэлектронному оператору; 2) сумму по занятым МО произведений коэффициентов на двухэлектронные интегралы, соответствующие оператору электрон-электронного отталкивания. Последние можно представить как произведение матрицы плотности (3.3.5) и матрицы двухэлектронных интегралов:

где – элементы матрицы плотности. (3.3.5)

Выражение (3.3.6) – более компактная форма матричного элемента матрицы Фока:

или:

F = h + G ^. P . (3.3.7)

Общая энергия в терминах интегралов базисных функций принимает следующий вид (3.3.8):

или (3.3.9), (3.3.10):

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

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

Для построения матрицы Фока необходимо вычислить ряд интегралов:

1) интегралы от двух функций и одноэлектронного оператора (дляМ базисных функций потребуется вычислить таких интегралов) – эти одноэлектронные интегралы называютóстовными, и они описывают взаимодействие электрона со всем набором ядер;

2) интегралы от четырёх функций и двухэлектронного оператора ; их количество пропорционально, они называютсядвухэлектронными. Для двухэлектронных интегралов 4 базисные функции могут быть локализованы на 1, 2, 3 или 4 различных атомных центрах. Соответствующие интегралы называются одно-, двух-, трех- и четырёхцентровыми. Подобные интегралы гораздо проще вычисляются при использовании функций Гаусса.

Таким образом, формально сложность нахождения решений уравнений Рутаана – Холла пропорциональна четвёртой степени размера базиса ().

С увеличением числа базисных функций точность вычисления МО улучшается. В предельном случае, при использовании бесконечного базиса, получается результат, идентичный численному решению уравнений Хартри – Фока. Этот результат называется пределом теории Хартри – Фока. Но этот предел не является точным решением уравнения Шрёдингера. Хартри – фоковский предел соответствует лишь наилучшей однодетерминантной волновой функции.