[GAMESS] phosphabenzene

ベンゼンの炭素原子を１つリン原子に置き換えたフォスファベンゼンの構造最適化。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
C 6.0 -3.2319609783 4.2474004795 0.0005148532
C 6.0 -4.3493730452 3.4117748214 -0.0023482711
C 6.0 -4.2746495965 2.0198651773 0.0009464602
C 6.0 -3.0383187494 1.3747397225 0.0074783662
C 6.0 -1.8156417002 2.0446899055 0.0111350406
P 15.0 -1.6348652520 3.7179553475 0.0083025766
H 1.0 -5.1889108059 1.4315590665 -0.0016005813
H 1.0 -5.3355376129 3.8750566004 -0.0074490582
H 1.0 -0.8953229088 1.4661159780 0.0162078390
H 1.0 -3.0261579769 0.2848993782 0.0098644932
H 1.0 -3.3772786511 5.3242458695 -0.0023989107
$END
ここまで-----

…重い原子が入ると計算時間が長くなります。
炭素−リン原子結合が炭素−炭素原子よりも長いので、ベンゼンのような正六角形ではなくなっていますが、平面性は保たれています。



↑最適化構造のHOMO(-6.56 eV)。リン原子のp軌道、でかいです。しかも中心部に何かいます。



↑LUMO(-1.25 eV)。炭素原子のp軌道とリン原子のp軌道の形がかなり異なっているのが解るかと。