[GAMESS] アニソールの構造最適化

アニソールの構造最適化。B3LYP/6-31G(d)。

-----ここから入力ファイルの中身
! DFT, B3LYP/6-31G(d), optimization
$CONTRL DFTTYP=B3LYP SCFTYP=RHF RUNTYP=OPTIMIZE $END
$SYSTEM TIMLIM=600000 MEMORY=80000000 $END
$STATPT OPTTOL=0.0001 NSTEP=200 PROJCT=.FALSE. $END
$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 $END
$GUESS GUESS=HUCKEL $END
$SCF DIRSCF =.TRUE. $END
$DATA

C1
C 6.0 -2.3788041087 3.7700777149 -0.0045667344
C 6.0 -3.6097404063 3.1173033583 -0.0111821668
C 6.0 -3.6553150440 1.7255140675 -0.0045611798
C 6.0 -1.1897169321 3.0334600326 0.0090878101
C 6.0 -1.2297952595 1.6386660305 0.0163387948
C 6.0 -2.4686884916 0.9914506840 0.0089193666
H 1.0 -4.6122211549 1.2099047574 -0.0098104120
H 1.0 -4.5322427059 3.6924496090 -0.0216236257
H 1.0 -0.2523531413 3.5793100841 0.0134149355
H 1.0 -2.3420204316 4.8569149158 -0.0100954429
H 1.0 -2.5073527911 -0.0952663449 0.0140747427
O 8.0 -0.1475765734 0.8015185110 0.0299363804
C 6.0 1.1357179043 1.4126240688 0.0496353749
H 1.0 1.3028009127 2.0097752490 -0.8529194671
H 1.0 1.8853060344 0.6153991373 0.0641982435
H 1.0 1.2729840585 2.0146669949 0.9540052888
$END
ここまで-----

電子供与基がベンゼン環上にあるため、HOMO, LUMO準位がベンゼンのもの(-6.64 eV, +0.16 eV)よりも高く、HOMO-LUMO準位差は小さくなっています。
↓最適化構造のHOMO (-5.77 eV)。酸素−メチル炭素結合がベンゼン環と同じ平面上にあることに注目。酸素原子のp軌道がHOMOに関与しているのもわかります。



↓一方LUMO (+0.19 eV)。こちらには酸素原子のp軌道は関与していません。