Abinitio


1. Оптимизация нафталена и азулена MOPAC.


Нафталин сразу оптимизировался, а для азулена пришлось применить силовое поле UFF
Наталин

Азулен

2. Оптимизация нафталена и азулена GAMESS.



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

3. Рассчет энергии.



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

4. Рассчет четырех систем.



Вещество E Naptalene E Azulene дельтаE Hartree дельтаE kCal/mol
Хартри-Фок -383.3549061523 -383.2825146536 0,07243713 45,45495126
DFT -385.6401306541 -385.5857339240 0,05439673 34,13443771


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