CC BY
6
Метод улучшения сходимости одноцентрового разложения молекулярной
орбитали
А. С. Каспржицкий, В.А. Явна Ростовский государственный университет путей сообщения
Аннотация: Исследование взаимодействия электромагнитного излучения с молекулярными системами дает большую часть информации об их структуре и свойствах. Интерпретация экспериментальных данных напрямую определяется знанием структуры энергетических уровней и ее изменением при переходе данных систем в возбужденное состояние. Ключевой задачей методов расчета молекулярных орбиталей возбужденных состояний является точное описание возникающих вакансий молекулярного остова, приводящих к радиальной релаксации электронной плотности. В работе предложен метод точного описания электронной плотности молекулярной системы вблизи ее ядер, обеспечивающий решения задачи медленной сходимости молекулярных орбиталей в одноцентровом представлении. Реализация вычислительной процедуры рассмотрена на примере двухатомной молекулы фтороводорода. Проведена оценка энергетических характеристик основного и ионизированных состояний молекулы.
Ключевые слова: одноцентровый метод, молекулярная орбиталь, возбужденное состояние, ионизированное состояние, глубокая оболочка.
Теоретическое исследование процессов поглощения, фотоионизации, упругого и неупругого рассеяния рентгеновского фотона молекулярной системой является сложной задачей, требующей точного описания возбужденных состояний с возникающими вакансиями в глубоких и субвалентных электронных оболочках[1-4]. Существующие методы расчета волновых функций можно классифицировать по форме их представления как многоцентровые (МЦ) [5] и одноцентровые (ОЦ) [6].
В МЦ представлении молекулярные орбитали (МО) в разных областях пространства описываются относительно центров, связанных с ядрами молекулярной системы. В таком подходе основное ограничение связанно с выбором базисного набора волновых функций и влияет на точность и объём вычислений [7,8]. При проведении расчетов набольшее распространение получили построения базисных функций в виде суперпозиции слейтеровских [9] или гауссовых [10] орбиталей. Наиболее известная реализация МЦ подхода является метод линейной комбинации атомных орбиталей (МО ЛКАО) [11, 12] уравнения, которого были впервые записаны Рутаном [5]. Задача минимизации функционала энергии, в данном случае, сводится к нахождению коэффициентов разложения МО по базисным атомным орбиталям. Достаточно полный обзор работ по использованию МО ЛКАО, представленный в [13], показывает, что основной недостаток метода - принципиальные трудности в получении возбужденных
состояний молекулярных систем, обусловленных проблемой выбора оптимального базисного набора.
Иной подход заключается в распространении методов расчета волновых функций атомных электронов на молекулярные системы, получивших свое развитие в ОЦ методе расчета МО [6, 14,15]. ОЦ подход основан на представлении МО в виде разложения по базису функций, центрированных относительно одной точки пространства. Применение ОЦ метода имеет ряд преимуществ, особенно при изучении процессов электрон-фотонного взаимодействия в многоэлектронных системах [15-17]. С вычислительной точки зрения, простота использование ОЦ по сравнению с МЦ обусловлена развитой теорией и методами атомных расчетов. Переход к молекулярным системам определяется введением суммирования по орбитальным моментам, поскольку формально одноэлектронная атомная волновая функция есть составляющая ОЦ разложения МО. Данный подход позволяет избежать проблем, связанных с избыточностью базиса, или по крайней мере контролировать их в отличие от МЦ подхода. Кроме того, отсутствие многоцентровых интегралов позволяет существенно упростить процедуру проведения расчетов.
Исторически первый способ реализации ОЦ метода заключался в представлении парциальных гармоник ОЦ разложения в виде линейной комбинации орбиталей слейтеровского типа и непосредственной минимизации выражений для полной энергии [6,11]. В ином подходе, подробно описанном в работах [16,18], предложен другой вариант ОЦ метода, который сводится к решению секулярного уравнения. При этом, в качестве базисных функций используются как функции дискретного, так и функции непрерывного спектра, которые являются решением уравнений Хартри-Фока. Данная реализация ОЦ подхода несвободна от проблем подбора оптимальных базисных наборов, а также от медленной сходимости рядов ОЦ разложения МО и сильного влияния на энергию состояния парциальных гармоник ОЦ разложения МО с большими значениями орбитального момента при низком вкладе этих гармоник в норму функции. Последние два недостатка также можно отнести к ОЦ методу в целом.
Применение вариационной принципа для получения волновых функций электронов в случае атома приводит к широко известной системе интегро-дифференциальных уравнений Хартри-Фока. В молекулярных системах ОЦ представление позволяет получить систему интегро-дифференциальных уравнений, связывающих разные парциальные гармоники МО. Такой подход является наиболее точным благодаря
ликвидации проблемы подбора базисных наборов для описания МО. В работах [17,19] данный метод был реализован при расчете волновой функции возбужденного фотоэлектрона дискретного и непрерывного спектров на основе системы связанных интегро-дифференциальных уравнений. Вместе с тем, существенное влияние на точность расчищаемых значений оказывает медленная сходимость МО в ОЦ представлении. Согласно оценки, выполненной в [19], для описания одноэлектронных энергий МО молекул с лигандами второго периода таблицы Менделеева с точностью ~3% в ОЦ разложении МО необходим учет парциальных гармоник с орбитальными моментами I >50. Рост числа слагаемых ОЦ разложения вызывает увеличение системы дифференциальных уравнений, обусловливающее расширение требуемых вычислительных мощностей с незначительным вкладом, получаемых высших гармоник разложения в наблюдаемые величины, но важным для энергетических характеристик систем. В работах [17,20,21] проблема медленной сходимости ОЦ разложения МО решена «эффективным» учетом влияния парциальных гармоник с большими значениями орбитального момента, путем увеличения коэффициента участия последней из учитываемых парциальной гармоники в ограниченном ОЦ разложении при расчете МО фотоэлектрона. Вместе с тем, такой «действенный» прием при расчете отдельных спектроскопических величин, включая сечение фотопоглощения, не позволяет правильно описать межэлектронное кулоновское взаимодействие.
В работе предложен метод точного описания электронной плотности молекулярной системы вблизи ее ядер. Численная процедура основана на решении системы связанных интегро-дифференциальных уравнений для всех МО, входящих в электронную конфигурацию состояния.
Интегро-дифференциальные уравнения, лежащие в основе итерационной схемы получения радиальных частей МО, являются обобщением уравнений Хартри-Фока впервые записанных для атома [22]. Уравнения получены на основе минимизации функционала энергии в приближении Борна-Оппенгеймера [23] для волновой функции Ы-электронной молекулярной системы ОТ'^ййЦ^И, представляемой в виде антисимметризованного произведения одноэлектронных МО фпу. После процедуры варьирования мы получаем набор систем связанных
интегро-дифференциальных уравнений Хартри-Фока, каждая из которых соответствует определенной МО пу электронной конфигурации молекулы. При этом, сама одноэлектронная МО в нерелятивистском ОЦ
представлении имеют вид [6]:
(1)
где (г) — радиальная и У^ (А, у) —угловая части парциальной гармоники ; —симметрии МО с фиксированным (для линейной молекулы) значением Ц — проекции момента импульса на ось квантования 02 (ось молекулы), п — главное квантовое число, у — симметрия МО и г^^ - сферические координаты.
С учетом (1) система связанных интегро-дифференциальных уравнений для МО пу имеет вид:
(2)
где (г) —локальный потенциал, описывающий взаимодействие
электронов МО пу с ядрами молекулярной системы и прямое кулоновское взаимодействие с электронами других МО, ^ (г) —нелокальный потенциал, описывающий обменное кулоновское взаимодействие электронов МО пу с электронами других МО.
Для точного описания электронной плотности в близи ядер молекулярной системы и решения задачи медленной сходимости функционального ряда (1) МО пу при решении системы связанных дифференциальных уравнений (2), представлена в виде:
= А,-... (3)
^^ = (г) Р? (т) + 2ЯГСг),
СЙ-2
где А— нормировочный коэффициент, Л.™^ — коэффициент участия волновой функции ф^^(г) — с главным квантовым числом п и симметрией у1 - ого лиганда молекулярной системы, имеющая вид:
= (4)
Функции (г) представляют собой разложенные в ряд по
сферическим гармоникам волновых функции электронных оболочек атомов, входящих в молекулярную систему в качестве лигандов (атомов не находящихся в начале координат).
В данном случае, парциальные гармоники МО в (3) со значениями орбитального момента получены решением системы уравнений (2).
Члены разложения (3) с орбитальными моментами 1> получены на основе линейной комбинации (г) с коэффициентами участия Л'Д". Такое
представление МО обусловлено тем, что форма парциальных сферические гармоники с большими значениями орбитального момента с точностью до коэффициента определяется потенциалом лиганда молекулярной системы.
Выбор значения 10 в (3) осуществляется в результате численного эксперимента и зависит от характера молекулярной системы и количества парциальных гармоник, участвующих в образовании химической связи соединения.
Для определения коэффициентов участия в (3) осуществляется
сшивка парциальных гармоник, используя равенство:
= (5)
То есть, радиальная часть 10 —ой парциальной гармоники пу —молекулярной орбитали представляется в виде суперпозиции радиальных парциальных гармоник функций лигандов с весовыми коэффициентами .
:
Выражение (5) позволяют составить систему алгебраических
ayL fta
уравнений для определения коэффициентов участия Афункций лигандов
с симметрией у в виде:
Е^'ГК'^^МСМ' (6)
Определенная, таким образом,МО нормируются на единицу.
В качестве объекта для проведения апробации численной процедуры выбрана молекула фтороводорода (ИБ). При расчетах выбрана система отсчета, связанная с атомом Б. Для систематизации одноэлектронных состояний использованы квантовые числа, характерные для молекул точечной группы симметрии Расчеты МО молекулы НБ выполнены
дляравновесного расстояния между ядрами атомов водорода и фторайр = 1.733 а.и. [24].Сравнительные расчеты энергетических
характеристик (одноэлектронные энергии МО и полная энергия) основного состояния 1 (7г гЗсрг1?ЕМолекулы НБ, представленные в Таблице №1.
Таблица № 1
Абсолютные значения
одноэлектронныеэнергииМОиполнаяэнергиямолекулыНР (в атомных единицах)
Тип 1g 2g 3g 1n Полная энергия
a 26.2923 1.5990 0.7674 0.6464 100.0577
b 26.2905 1.5995 0.7656 0.6492 100.0594
c 26.2915 1.5995 0.7656 0.6492 100.0616
В пункте (а) Таблицы №1 представлен расчет, выполненный МЦ методом МО ЛКАО по программе General Atomic and Molecular Electronic Structure
System [12]. Результаты расчетов, представленные в пункте (b) Таблицы №1, получены в данной работе с использованием МО разложенных в ряд (1) и
ограниченных t <1 ft, показывают достаточную точность, которая при оценке
энергетических характеристик не уступает МЦ методам расчета. Использование высших гармоник волновых функций лигандов для точного описания МО в области ядер молекул дает дополнительный вклад в полную энергию и закономерно меняет энергетический уровень 1о-состояния молекулы HF. Наблюдаемые в пункте (с) Таблицы №1 изменения энергетических величин обусловлены, прежде всего, вкладом 1s-состоянияатома H, являющимся лигандом в данном случае. Построение МО
определяется выбором Iq, включающим в структуру МО функции лигандов.
Существенное влияние высших гармоник ОЦ разложения на энергетические характеристики молекулярной системы наблюдается при росте заряда лигандаб?£.. Для оценки данного влияния проведены расчета полной энергии электронной конфигурации З^'З^'Зйр' It;^ зависимости от заряда который вместо случая ZL= 1 менялся в диапазоне SL =2+5. Результаты расчета полных энергий представлены в Таблице № 2.
Таблица № 2
Полнаяэнергияэлектронного состояния Icp' 1'ге4, (в атомных
единицах)
Полная энергия (a) Полная энергия(Ь)
2 -101.0902 -101.0951
3 -102.8639 -102.9189
4 -105.1593 -105.4553
5 -107.7773 -108.5403
В столбце (а) Таблица № 2приведен расчет полных энергий на основе МО с ОЦ разложением ограниченным I 4. Расчет, учитывающий функции лиганда с разложением до I = 60, представлен в столбце (Ь) Таблица № 2. Видно, что с ростом заряда ZL высшие гармоники ОЦ разложения МО начинаю существенно влиять на энергетические характеристики молекулярной системы. Результаты расчета МО электронного состояния показывают, что данная итерационная схема решает основную проблему ОЦ методов, а именно, описания поведения МО в области ядер лиганда. При этом видно, что с ростом заряда лиганда существенно увеличивается вклад высших гармоник МО, учитываемых данным ОЦ подходом.
Таким образом, представленная итерационная схема ОЦ расчета МО становится эффективным при исследовании молекул с неводородными лигандами. Результаты исследования показывают, что учет высших гармоник разложения волновых функций электронных оболочек лигандов является необходимым условием при проведении расчетов в ОЦ представлении для получения точных энергетических характеристик молекулярных систем.
Работа выполнена при поддержке ФГБОУ ВО «Ростовский государственный университет путей сообщения» (Договор о предоставлении гранта ФГБОУ ВО «РГУПС» № 30/1 от 27.01.2017 года).
Литература
1. Dutoi A.D., Leone S.R. Simulation of X-ray transient absorption for following vibrations in coherently ionized F2 molecules// Chem. Phys. 2017. V. 482.pp. 249-264.
2. Смирнова Ю.О., Положенцев О.Е., Леонтьева Д.В., Чайников А.П., Сучкова С.А., Гуда А.А., Ломаченко К.А., Смоленцев Н.Ю., Подковырина Ю.С., Солдатов М.А., Кравцова А.Н., Солдатов А.В. Разработка нового
комплексного метода определения параметров 3D наноразмерной атомной и электронной структуры материалов на основе методик XAFS, XRD и Raman
ivdon.ru/magazine/archive/n4p1y2012/1268.
3. Yavna V., Hopersky A., Nadolinsky A., Yavna S. Orientation effects in elastic scattering of polarized X-rays by linear molecules// J. Synchrotron Radiat. 2001.V. 8.pp. 240-242.
4. Кременная М.А., Солдатов М.А., Чайников А.П., Подковырина Ю.С., Бугаев А.Л., Ломаченко К.А., Кравцова А.Н. Рентгеноспектральное исследование и компьютерное моделирование локальной атомной структуры центра связывания иона меди в бета амилоиде //Инженерный вестник Дона, 2013, №2URL: ivdon.ru/magazine/archive/n2y2013/1640.
5. Roothaan C.C.J. Self-Consistent Field Theory for Open Shell of Electronic Systems// Rev. Mod. Phys. 1960.V. 32. pp. 179-185.
6. Bishop D.M. Single-Center Molecular Wave Function. In Advances in Quantum Chemistry. Holland: Elsevier, 1967.408 p.
7. Pople J.A. Molecular orbital methods in organic chemistry // Acc. Chem. Res. 1970.V. 3.pp. 217-223.
8. Wilson S. Electron Correlation in Molecules. USA: Dover Publications, 2007. 304 p.
9. Bunge C.F., Barrientos J.A., Bunge, A.V., Cogordan J.A. Hartree-Fock and Roothaan-Hartree-Fock energies for the ground states of He through Xe // Phys. Rev. A. 1992.V 46.pp. 3691-3696.
10. Huzinaga S., Gaussian-Type Functions for the Polyatomic Systems//J. Chem. Phys. 1965.V. 42. pp. 1293-1302.
11. Roothaan C.C.J. New Development in Molecular Orbital Theory // Rev. Mod. Phys. 1951.V.23.pp. 69-89.
// Инженерный
вестник
Дона, 2012, №4 URL:
12. Schmidt M.W., Baldridge K.K., Boatz J.A., et al. General atomic and molecular electronic structure system// J. Comput. Chem. 1993. V. 14. pp. 1347— 1363.
13. Lowe J.P. Quantum Chemistry. USA: Academic Press, 1978. 599 p.
14. Faisal F.H.M. Electron-Molecular Interactions: I-Single Center Wave-Functions and Potentials// J. Phys. B: At. Mol. Phys. 1970.V. 3.pp. 636-640.
15. Yavna V.A., Nadolinsky A.M., Hopersky A.N. Theoretical study of processes of multiple excitation/ionization in 2s-photoabsorption of the CO molecule // J. Electr. Spectr. Relat. Phen. 1998. V. 94.pp. 49-57.
16. Yavna V.A., Nadolinsky A.M., Demekhin V.F. Theoretical study of inner shell photoabsorption spectra of simple molecules// J. Electr. Spectr. Relat. Phen. 1994. V. 68.pp. 267-275.
17. Demekhin Ph. V., Omel'yanenko D. V., Lagutin B. M.,et al. Investigation of photoionization and photodissociation of an oxygen molecule by the method of coupled differential equations // Opt. Spectrosc. 2007. V. 103.pp. 318-329.
18. Yavna V.A., Nadolinsky A.M., Hopersky A.N. Theoretical study of processes of multiple excitation/ionization in 2s-photoabsorption of the CO molecule// J. Electr. Spectr. Relat. Phen. 1998. V. 94.pp. 49-57.
19. Lagutin B. M., Migal Yu. F. Single-center method of calculation for clusters and molecules other than hydrides//Theor. Exp. Chem. 1987.V. 25. pp.1017
20. Yavna V.A., Nadolinsky A.M., Demekhin V. Ph. Theoretical study of photoabsorption cross sections in the regions of CK- and NK-ionization thresholds of CO and N2 molecules// Opt. Spectrosc. 1990. V. 69.pp. 1278-1284.
21. Yavna V.A., Nadolinsky A.M., Demekhin V. Ph. Theoretical study of XANES 2o-shells CO // Opt. Spectrosc. 1992. V. 73.pp. 1081-1085.
22. Cacelli I. The Hartree-Fock Method. In Reference Module in Chemistry, Molecular Sciences and Chemical Engineering.Holland: Elsevier, 2015.14 p.
23. Sobel' Man I.I. Introduction to the Theory of Atomic Spectra.USA: Pergamon Press, 1972. 626 p.
24. Kane P.P. Inelastic scattering of X-rays and gamma rays by inner shell electrons // Phys. Reports. 1992.V. 218.pp. 67-139.
References
1. Dutoi A.D., Leone S.R. Chem. Phys. 2017. V. 482. pp. 249-264.
2. Smirnova Ju.O., Polozhencev O.E., Leont'eva D.V., Chajnikov A.P., Suchkova S.A., Guda A.A., Lomachenko K.A., Smolencev N.Ju., Podkovyrina Ju.S., Soldatov M.A., Kravcova A.N., Soldatov A.V. Inzhenernyj vestnik Dona (Rus), 2012, №2. URL: ivdon.ru/magazine/ archive/ n4p1y2012/1268.
3. Yavna V, Hopersky A, Nadolinsky A, Yavna S. J. Synchrotron Radiat. 2001. V. 8. pp. 240-242.
4. Kremennaja M.A., Soldatov M.A., Chajnikov A.P., Podkovyrina Ju.S., Bugaev A.L., Lomachenko K.A., Kravcova A.N. Inzhenernyj vestnik Dona (Rus), 2013, №2. URL:ivdon.ru/magazine/ archive/n2y2013/1640.
5. Roothaan C.C.J.Rev. Mod. Phys. 1960. V. 32. pp. 179-185.
6. Bishop D.M. Single-Center Molecular Wave Function. In Advances in Quantum Chemistry. Holland: Elsevier, 1967. 408 p.
7. Pople J.A. Acc. Chem. Res. 1970. V. 3. pp. 217-223.
8. Wilson S. Electron Correlation in Molecules. USA: Dover Publications, 2007. 304 p.
9. Bunge C.F., Barrientos J.A., Bunge, A.V., Cogordan J.A. Phys. Rev. A. 1992. V 46. pp. 3691-3696.
10. Huzinaga S.J. Chem. Phys. 1965. V. 42. pp. 1293-1302.
11. Roothaan C.C.J. Rev. Mod. Phys. 1951. V.23. pp. 69-89.
12. Schmidt M.W., Baldridge K.K., Boatz J.A., et al. J. Comput. Chem. 1993. V. 14. pp. 1347-1363.
13. Lowe J.P. Quantum Chemistry. USA: Academic Press, 1978. 599 p.
14. Faisal F.H.M. J. Phys. B: At. Mol. Phys. 1970. V. 3. pp. 636-640.
15. Yavna V.A., Nadolinsky A.M., Hopersky A.N. J. Electr. Spectr. Relat. Phen. 1998. V. 94. pp. 49-57.
16. Yavna V.A., Nadolinsky A.M., Demekhin V.F. J. Electr. Spectr. Relat. Phen. 1994. V. 68. pp. 267-275.
17. Demekhin Ph. V., Omel'yanenko D. V., Lagutin B. M., et al. Opt. Spectrosc. 2007. V. 103. pp. 318-329.
18. Yavna V.A., Nadolinsky A.M., Hopersky A.N. J. Electr. Spectr. Relat. Phen. 1998. V. 94. pp. 49-57.
19. Lagutin B. M., Migal Yu. F. Theor. Exp. Chem. 1987. V. 25. pp.10-17
20. Yavna V.A., Nadolinsky A.M., Demekhin V. Ph. Opt. Spectrosc. 1990. V. 69. pp. 1278-1284.
21. Yavna V.A., Nadolinsky A.M., Demekhin V. Ph. Opt. Spectrosc. 1992. V. 73. pp. 1081-1085.
22. Cacelli I. The Hartree-Fock Method. In Reference Module in Chemistry, Molecular Sciences and Chemical Engineering. Holland: Elsevier, 2015. 14 p.
23. Sobel' Man I.I. Introduction to the Theory of Atomic Spectra. USA: Pergamon Press, 1972. 626 p.
24. Kane P.P. Phys. Reports. 1992. V. 218. pp. 67-139.
