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ではなくもっと大きな基底関数を用いるべきである。


Return

2012 Copyright by Hiroshi Kihara