Abinitio вычисления для нафталина и азулена




1. Оптимизация с помощью MOPAC структуры нафталина и азулена

Нафталин оптимизируется сразу в плоскую структуру без дополнительных параметров.


Азулен же становится плоским только при силовом поле UFF.


2. Оптимизация геометрии средствами GAMESS

Берем файлы nap.out и azu.out
С помощью babel переформатируем координаты в gamin формат и добавим нужные строки: nap_opt.inp и azu_opt.inp
После оптимизации получили nap_opt.log и azu_opt.log

3. Расчёт энергии методом Хартри-Фока и с помощью теории функционала плотности

Составляем файлы:
Хартри-Фок: nap_opt_hf.inp и azu_opt_hf.inp
Теория функционала плотности: nap_opt_dfp.inp и azu_opt_dfp.inp
Рассчеты энергии:
Хартри-Фок: nap_opt_hf.log и azu_opt_hf.log
Теория функционала плотности: nap_opt_dfp.log и azu_opt_dfp.log

4. Таблица энергий

Вещество

E Naptalene

E Azulene

dE, Hartree

dE, kCal/mol

Хартри-Фок

-383.3549061523

-383.2824690219

0.0724371304

45.4549512601736

DFT

-385.6401306541

-385.5857491456

0.0543815085

34.1248860173265



Известно, что энергия изомеризации нафталина в азулен составляет 35.3±2.2 kCal/mol. Видно, что метод теории функционала плотности даёт результат близкий к экспериментальным, выходит, он лучше.


© Alisa Garaeva