2010/07/15

[GAMESS] benzyl cation

ベンジルカチオンのπ軌道。
DFT, B3LYP/6-31G(d)で構造最適化&計算。
カチオンの構造最適化なので、入力ファイルの中身に

$CONTRL ICHARG=1 $END

が入っています。

ここから入力ファイルの中身-----
! DFT, B3LYP/6-31G(d), optimization
$CONTRL DFTTYP=B3LYP SCFTYP=RHF RUNTYP=OPTIMIZE $END
$CONTRL ICHARG=1 $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 -5.2029383210 1.8129269588 -0.0002477590
C 6.0 -4.4049708599 2.9581461642 -0.0035676844
C 6.0 -3.0086345385 2.8500015771 -0.0175311352
C 6.0 -4.6121347346 0.5503118155 -0.0002421480
C 6.0 -3.2235334267 0.4327450198 0.0010405267
C 6.0 -2.4245292607 1.5778037838 -0.0019997265
C 6.0 -2.1531967986 4.0816451522 0.0141818071
H 1.0 -4.8794496390 3.9367082607 0.0040439345
H 1.0 -6.2854702705 1.9058876720 0.0053649472
H 1.0 -5.2335421396 -0.3407920442 0.0032410794
H 1.0 -2.7617465824 -0.5507996887 0.0074553494
H 1.0 -1.3423347032 1.4700722064 0.0066902708
H 1.0 -2.3299507141 5.0983747701 0.0076291540
H 1.0 -1.1462179877 3.8588805268 0.0515127073
$END
ここまで-----

最適化構造。全ての原子が同一平面上にあります。
↓HOMO' (-12.51eV)。カチオンの場合、ここまで電子が配置されます。

100715_1.png

↓HOMO (-8.76 eV)。

100715_2.png

↓LUMO (-5.31 eV)。

100715_3.png

モデリング: Avogadro
http://avogadro.openmolecules.net/wiki/Main_Page
軌道計算: GAMESS-US
http://www.msg.chem.iastate.edu/gamess/
軌道表示: MacMolPlt
http://www.scl.ameslab.gov/MacMolPlt/

0 件のコメント:

コメントを投稿