2010/07/09

[GAMESS] 1,3-butadiene

1,3-ブタジエンのπ軌道。
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 -5.1239939214 1.2238560269 0.1059725242
C 6.0 -3.9908907316 1.8938823772 -0.1321032023
C 6.0 -2.6840225261 1.3140294664 0.0535620469
C 6.0 -1.5502889350 1.9832495223 -0.1848541474
H 1.0 -2.6276783633 0.2853193784 0.4036097812
H 1.0 -1.5473427010 3.0102376992 -0.5349133184
H 1.0 -0.5871068582 1.5062066961 -0.0310976074
H 1.0 -6.0874912222 1.7003021720 -0.0475774738
H 1.0 -5.1263768737 0.1966558159 0.4554607173
H 1.0 -4.0474723492 2.9227138823 -0.4816780598
$END
ここまで-----

↓最も準位の低いπ軌道 (-8.60 eV)。

100709_1.png

↓HOMO (-6.15 eV)。

100709_2.png

↓LUMO (-0.52 eV)。

100709_3.png

↓最も準位の高いπ軌道 (+2.42 eV)。

100709_4.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/

2010/07/08

[GAMESS] ethylene

エチレンのπ軌道。
DFT, B3LYP/6-31G(d)で構造最適化&計算。一瞬で終了。

ここから入力ファイルの中身-----
! DFT, B3LYP/6-31+G(d)
$CONTRL DFTTYP=B3LYP SCFTYP=RHF RUNTYP=ENERGY MAXIT=100 $END
$SYSTEM TIMLIM=600000 MEMORY=80000000 $END
$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 DIFFSP=.TRUE. $END
$GUESS GUESS=HUCKEL $END
$SCF DIRSCF =.TRUE. FDIFF =.FALSE. $END

$DATA

C1
C 6.0 -5.7096904870 3.9866513087 0.0000000000
C 6.0 -4.5373174627 3.3464410169 0.0000000000
H 1.0 -3.6001887723 3.8946627780 -0.0000000000
H 1.0 -4.4913366989 2.2621372591 0.0000000000
H 1.0 -5.7556435437 5.0710436024 0.0000000000
H 1.0 -6.6467596626 3.4385008169 -0.0000000000
$END
ここまで-----

↓HOMO (-7.46 eV)。

100708_1.png

↓LUMO (-0.14 eV)。LUMOと言えども準位が正の値ではないのです。

100708_2.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/

2010/07/07

[GAMESS] benzene

ベンゼンのπ軌道。
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 -2.1948388776 2.4300766609 1.9586799938
C 6.0 -3.5627489207 2.5576768934 1.7173704112
C 6.0 -4.1730092146 1.7793859170 0.7340566694
C 6.0 -1.4373995572 1.5241684676 1.2166003291
C 6.0 -2.0477117313 0.7457226721 0.2331100462
C 6.0 -3.4154113550 0.8733398883 -0.0081240342
H 1.0 -5.2387893116 1.8787699275 0.5460102749
H 1.0 -4.1529758331 3.2636585766 2.2956933732
H 1.0 -0.3715676048 1.4246847191 1.4045563757
H 1.0 -1.7193639208 3.0366049843 2.7249551908
H 1.0 -3.8910176198 0.2668446897 -0.7744005102
H 1.0 -1.4574078021 0.0398130751 -0.3451162192
$END
ここまで-----

↓最も準位の低いπ軌道 (-9.69 eV)。黄色と緑色の境目が1つしかないのに注目。

100707_1.png

↓HOMO (-6.61 eV)。2つあります。

100707_2.png

100707_3.png

↓LUMO (+0.19 eV)。これも2つ。

100707_4.png

100707_5.png

↓最も準位の高いπ軌道 (+4.54 eV)。黄色と緑色の境目がたくさん。

100706_6.png

↓お気に入り (+2.56 eV)。

100707_7.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/

2010/07/06

[GAMESS] trimethylsilyldiazomethane

ジアゾメタンの代替として有機合成に使われるトリメチルシリルジアゾメタン。

↓モデリングして分子力場計算で構造を最適化すると…ジアゾメタンと同じく、曲がってました。

100706_0.png

得られた座標データに計算パラメータを加えて、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 -8.8566933285 1.5002915810 -0.0509405067
C 6.0 -7.0363636916 1.9081656074 -0.1162174194
C 6.0 -9.8472844728 3.0841892299 0.0838504429
C 6.0 -9.1991354679 0.4179617066 1.4384486976
C 6.0 -9.3408792803 0.5898979621 -1.6137373542
N 7.0 -6.7259076815 2.6701355531 1.0481225144
H 1.0 -6.7926242282 2.4969697649 -1.0049368740
H 1.0 -6.4253042278 1.0002200468 -0.1245381427
H 1.0 -9.6293599827 3.6021095454 1.0225274977
H 1.0 -10.9203442337 2.8721164454 0.0544418481
H 1.0 -9.6103791731 3.7606947647 -0.7425696000
H 1.0 -8.9522469484 0.9420102949 2.3665742329
H 1.0 -8.6036809143 -0.4989887235 1.3955403038
H 1.0 -10.2556329613 0.1367158888 1.4770449761
H 1.0 -9.1346052910 1.2009508699 -2.4975321076
H 1.0 -10.4082040768 0.3494320808 -1.6069717817
H 1.0 -8.7809034367 -0.3453256006 -1.7070955760
N 7.0 -6.4678666386 3.2546503813 1.9288267484
$END
ここまで-----

↓HOMO (-5.69 eV)。炭素−窒素−窒素のところがまっすぐになりました。

100706_1.png

↓LUMO (-0.87 eV)。HOMO, LUMO共にジアゾメタン(5.90 eV/1.20 eV)よりも高く、HOMO-LUMO準位差(4.82 eV/4.70 eV)は大きくなってます。

100706_2.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/

2010/07/05

[GAMESS] diazomethane

メチル化剤等として有機合成に使われる(ただし毒性あり)ジアゾメタン。

↓モデリングして分子力場計算で構造を最適化すると…折れ曲がってます(たぶん間違った構造)。

100705_0.png

得られた座標データに計算パラメータを加えて、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.5062850794 3.5210793468 0.0408353170
N 7.0 -2.5781507499 2.8632400298 0.5634773115
N 7.0 -1.6346723394 3.3535608528 0.8902702351
H 1.0 -3.4317177932 4.6092667170 -0.1184894116
H 1.0 -4.4219991822 3.0015812836 -0.2686474484
$END
ここまで-----

↓HOMO (-5.90 eV)。やはり直線状の構造でした。

100705_1.png

↓LUMO (-1.20 eV)。

100705_2.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/

2010/07/04

浜焼市場

しらすかきあげ丼

2010/07/03

[GAMESS] MePC

1-ホスファプロピンの構造異性体のCH3PC。アセトニトリルとメチルイソシアニドの関係と同じです。構造最適化はいつも通り、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 -6.3128585037 1.6437576812 0.0000000000
P 15.0 -4.7408818568 2.3213369961 0.0000000000
H 1.0 -7.0615670543 2.4081760916 0.0000000000
H 1.0 -6.4376658323 1.0098654853 0.8529383417
H 1.0 -6.4269680629 1.0641990942 -0.8921831946
C 6.0 -2.9193852014 3.0491078035 0.0000000000
$END
ここまで-----

↓HOMO (-6.20 eV)。

100702_1.png

↓LUMO (-2.10 eV)。1-ホスファプロピンよのものより、HOMOは高く、LUMOは低くなっています。

100702_2.png

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