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

1.Подготовка структур

SMILES записи нафталина и азулена:

C1=CC=C2C=CC=C2C=C1 Azulene 
C1CCC2CCCCC2C1 Naphtalene
          

Плоские структуры были получены при использования поля UFF для obgen и оптимизации в Mopac на уровне PM6

Нафталин	Азулен

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

 $CONTRL COORD=CART UNITS=ANGS SCFTYP=RHF RUNTYP=OPTIMIZE $END
 $BASIS  GBASIS=N31 NGAUSS=6  $END
 $SYSTEM MWORDS=2 $END
 $DATA
	  

3.Расчет энергии методами Хартри-Фока и DFT в GAMESS

Для расчета энергии на основе выходных данных предыдущего пункта подготовлены файлы для расчета энергии двумя методами для каждой молекулы.
Для метода Хартри-Фока нафталин, азулен c заголовком:

 $CONTRL COORD=CART UNITS=ANGS SCFTYP=RHF RUNTYP=ENERGY $END
 $BASIS GBASIS=N31 NGAUSS=6 POLAR=POPN31 NDFUNC=1 $END
 $GUESS GUESS=HUCKEL $END
 $SYSTEM MWORDS=2 $END
 $DATA

Для метода DFT нафталин, азулен с заголовком:

 $CONTRL COORD=CART UNITS=ANGS DFTTYP=B3LYP RUNTYP=ENERGY $END
 $BASIS GBASIS=N31 NGAUSS=6 POLAR=POPN31 NDFUNC=1 $END
 $GUESS GUESS=HUCKEL $END
 $SYSTEM MWORDS=2 $END
 $DATA
	  

4.Сравнение методов расчета энергии

Согласно экспериментальным данным, энергия изомеризации нафталина в азулен составляет 35.3±2.2 kCal/mol, что очень хорошо соответствует значению, полученному методу DFT. Метод Хартри-Фока, в данном случае, дает слишком большую ошибку.

5.HOMO/LUMO

Орбиталь / Вещество	HOMO	LUMO
Нафталин
Азулен