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](http://lh3.ggpht.com/_rjhP34DYJXE/TD79X3qhPEI/AAAAAAAAARo/OHdxG5yb620/100715_1.png?imgmax=800)
↓HOMO (-8.76 eV)。
![100715_2.png](http://lh5.ggpht.com/_rjhP34DYJXE/TD79vg-coDI/AAAAAAAAARs/qW31oBFoqHc/100715_2.png?imgmax=800)
↓LUMO (-5.31 eV)。
![100715_3.png](http://lh6.ggpht.com/_rjhP34DYJXE/TD7-CfWy33I/AAAAAAAAARw/q2F8hvaVi5A/100715_3.png?imgmax=800)
モデリング: 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 件のコメント:
コメントを投稿