Визуализация электронной структуры атома кислорода

Основным состоянием атома кислорода, также как и молекулы, является триплет, т.е. полный электронный спин S=1. Электронная конфигурация атома кислорода:

1s2 ↑↓ 2s2 ↑↓ 2p4 ↑↓

Создаем новое активное окно. В библиотеке элементов выбираем атом кислорода и переносим его в активное окно. На панели меню редактора GaussView выбираем Calculate → Gaussian Calculation Setup. Устанавливаем следующие параметры расчетов: Job Type - energy, Method – DFT B3LYP 6-31(d)G, Spintriplet, Job Title – (например, O_DFT_B3LYP_6-31dG_triplet_energy). Нажимаем Submit, при сохранении файла для удобства лучше указывать то же имя файла, что и во вкладке Title. При запуске расчетов появится окно Gaussian03 (рис.2.2), в окне которого отображается текущее состояние расчетов.

Рис.2.2 Окно программного комплекса Gaussian 03

Информация о распределении электронной плотности атомных и молекулярных орбиталей (АО и МО) содержится в output файле с расширением *.chk. После окончания расчетов открываем файл расширением O_DFT_B3LYP_6-31dG_triplet_energy.chk. На панели команд GaussView выбираем редактор молекулярных орбиталей (табл.1.1). В окне редактора МО отображена диаграмма уровней - распределение электронов по энергетическим уровням молекулярных орбиталей (рис.2.3). Энергия выражена в атомных единицах энергии 1a.u.= 1Хартри=27.2эВ. Alpha и Beta электроны соответствуют противоположным ориентациям электронных спинов.

Рис.2.3 Диаграмма уровней атома кислорода в триплетном состоянии.

Переходим во вкладку Visualize,выделяем мышью интересующие орбитали на диаграмме уровней (выбранные МО подсвечиваются желтым цветом) и нажимаем Update. Визуализированные АО представлены на рис.2.4. Для отображения декартовой системы координат, на панели команд выбираем View → Cartesian Axes.

Одной из характеристик атомов и молекулярных структур является потенциал ионизации- минимальная энергия, необходимая для удаления электрона с внешней оболочки на бесконечность. По теореме Купманса, потенциал ионизации определяется как энергия электронов на высшей занятой орбитали HOMO (High Occupied Molecular Orbital), взятая с противоположным знаком.

(2.12)

Поскольку энергия электрона на 2p↓ орбитали выше, чем на 2p↑, то потенциал ионизации атома кислорода, рассчитанный методом DFT B3LYP 6-31(d)G составляет 0.317 a.u.=8.622 эВ.

Задание 1. Визуализировать электронную плотность атомных орбиталей и соотнести с электронной конфигурацией атома (1s22s12p3) углерода, используя метод DFT B3LYP 6-31(d)G.

№ АО α-орбиталь EАО, a.u. β-орбиталь ЕАО, a.u.
2p↑ -0.369    
2p↑ -0.446    
2p↑ -0.446 2p↓ -0.317
2s↑ -1.006 2s↓ -0.848
  1s↑ -19.291 1s↓ -19.241

Рис. 2.4 Атомные орбитали атома кислорода. Масштаб сохранен