exam03: ROHF 法による三重項メチレン（カルベン）のグラジエント

Restricted open-shell Hartree-Fock (ROHF) 法による三重項メチレン（カルベン）の基底状態の計算です。

! The wavefunction is a pure triplet state (<S**2> = 2), and so has a higher energy than the second example. ! This job tests the ROHF wavefunction and gradient code. It also tests the direct SCF procedure. $CONTRL SCFTYP=ROHF MULT=3 RUNTYP=GRADIENT COORD=CART $END $SYSTEM TIMLIM=1 $END $SCF DIRSCF=.TRUE. $END $BASIS GBASIS=STO NGAUSS=2 $END $GUESS GUESS=HUCKEL $END $DATA Methylene...3-B-1 state...ROHF/STO-2G Cnv 2 Hydrogen 1.0 0.82884 0.7079 0.0 Carbon 6.0 Hydrogen 1.0 -0.82884 0.7079 0.0 $END

$CONTRL group

$CONTRL SCFTYP=ROHF MULT=3 RUNTYP=GRADIENT COORD=CART $END

計算方法はRestricted open-shell Hartree-Fock（制限開殻）法で、三重項状態（MULT=3）のエネルギーとグラジエントを計算する（RUNTYP=GRADIENT）

構造データの入力は、すべての原子の直交座標を記述する（COORD=CART）

$SYSTEM group

$SYSTEM TIMLIM=1 $END

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

$SCF group

$SCF DIRSCF=.TRUE. $END

Direct SCF法を使用する

$BASIS group

$BASIS GBASIS=STO NGAUSS=2 $END

基底関数としてSTO-2Gを用いる

$GUESS group

$GUESS GUESS=HUCKEL $END

最初の分子軌道の係数を求めるのに拡張ヒュッケル法を用いる

$DATA group

$DATA
Methylene...3-B-1 state...UHF/STO-2G     ; タイトル
Cnv 2      ; 対称群   回転軸の次数　この場合はC2vに相当する
     ; C1以外の場合は、通常空行を入れる
Hydrogen 1.0 0.82884 0.7079 0.0    ; 原子の名前（10文字以下）　原子番号　x y z
Carbon 6.0    ;　原子の名前　原子番号　x y z は0
Hydrogen 1.0 -0.82884 0.7079 0.0    ; 原子の名前　原子番号　x y z
$END

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

補足

実際に計算を行う際は、STO-2Gではなくもっと大きな基底関数を用いるべきである。