[GAMESS] silabenzene

ベンゼンの炭素原子がひとつケイ素原子に置き換わったシラベンゼン。
いつものようにDFT, 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
Si 14.0 -7.4577545904 0.4017921552 0.0131771064
C 6.0 -7.4496545260 -1.0694611592 1.1227846618
C 6.0 -6.2857046226 -1.7194162094 1.5486571923
C 6.0 -4.9830704958 -1.3330507560 1.1666995446
C 6.0 -4.6806302733 -0.2559986456 0.3061555292
C 6.0 -5.6457134914 0.5752623141 -0.2736846777
H 1.0 -5.3505513706 1.3788225517 -0.9385513271
H 1.0 -3.6307594385 -0.0672518076 0.0765974584
H 1.0 -4.1473810447 -1.9165389721 1.5590110288
H 1.0 -6.3875722222 -2.5807660387 2.2106734366
H 1.0 -8.4237530558 -1.4231484041 1.4404544231
H 1.0 -8.1421420789 0.0711870719 -1.2602812986
$END
ここまで-----

↓HOMO (-5.63 eV)。ベンゼンのものとも、ホスファベンゼンのものとも、微妙に違ってます。ケイ素原子に結合している水素原子が、ベンゼン環とほぼ同じ平面にあるのが驚き。
ケイ素原子の中心に何かいますね。



↓LUMO (-0.63 eV)。軌道や準位他、他のベンゼン類縁体と比べるとおもしろいかも。