exam23: PM3法によるSi3の構造最適化

! semiempirical calculation, using the MOPAC/GAMESS combo
$CONTRL SCFTYP=RHF RUNTYP=OPTIMIZE COORD=ZMT ICHARG=-1 $END
$SYSTEM TIMLIM=1 $END
$BASIS GBASIS=PM3 $END
$DATA
Silyl anion...comparison of semiempirical models
Cnv 3

Si
H 1 rSiH
H 1 rSiH 2 aHSiH
H 1 rSiH 2 aHSiH 3 aHSiH -1

rSiH=1.15
aHSiH=110.0
$END

$CONTRL group

 $CONTRL SCFTYP=RHF RUNTYP=OPTIMIZE COORD=ZMT ICHARG=-1 $END

計算方法は制限Hartree-Fock法（閉殻）で、構造最適化を行う（RUNTYP=OPTIMIZE）

構造データの入力にはGaussianスタイルの内部座標を用いる (COORD=ZMAT)　系の電荷は -1 ( ICHARG=-1)

$SYSTEM group

 $SYSTEM TIMLIM=1 $END

計算の時間制限は1分 (TIMLIM=1)　指定しない場合は5,25,600分

$BASIS group

 $BASIS GBASIS=PM3 $END

PM3ハミルトニアンを用いる

$DATA group

 $DATA
Silyl anion...comparison of semiempirical models      ; タイトル
Cnv 3      ; 対称群   回転軸の次数　この場合はC3v に相当する
     ; C1以外の場合は、通常空行を入れる
Si    ; Zマトリックスで構造を指定する。ケイ素を原点に置く
H 1 rSiH    ; 2番目の水素と1番目の原子の間の結合の長さはrSiH
H 1 rSiH 2 aHSiH   ; 3番目の水素と1番目の原子の間の結合の長さはrSiHで結合角はaHSiH
H 1 rSiH 2 aHSiH 3 aHSiH -1   ; 4番目の水素と1番目の原子の間の結合の長さはrSiHで、結合角 ∠4-1-2および∠4-1-3はaHSiH

rSiH=1.15
aHSiH=110.0
$END

距離の単位は$CONTR groupのUNITSで指定する。指定しない場合はÅ単位となる。角度の単位は度である。

Return

2012 Copyright by Hiroshi Kihara