exam08: RHF + MP2 法によるH₂Oのエネルギーとグラジエント

! This job generates RHF orbitals which should be saved for use with EXAM9. This run, together with EXAM9, ! shows a much more typical MCSCF calculation, which should always be started with some sort of SCF MOs. ! This job also tests the 2nd order Moller-Plesset code, exercizing the parallel code for this, even if p=1. $CONTRL SCFTYP=RHF MPLEVL=2 RUNTYP=GRADIENT $END $SYSTEM TIMLIM=1 MEMDDI=1 PARALL=.TRUE. $END $BASIS GBASIS=N21 NGAUSS=3 $END $GUESS GUESS=HUCKEL $END $DATA Water...RHF/3-21G...exp.geom...R(OH)=0.95781,A(HOH)=104.4776 Cnv 2 OXYGEN 8.0 HYDROGEN 1.0 0.0 0.7572157 0.5865358 $END

$CONTRL group

$CONTRL SCFTYP=RHF MPLEVL=2 RUNTYP=GRADIENT $END

計算方法は非制限Hartree-Fock法（開殻）で、2次のMoller-Plesset摂動計算（MPLEVL=2）により、エネルギーとグラジエントを計算する（RUNTYP=GRADIENT）

COORDが指定されていないためCOORD=UNIQUEとみなされる。構造データの入力は、対称性から一意な原子の直交座標のみを記述する

$SYSTEM group

$SYSTEM TIMLIM=1 MEMDDI=1 PARALL=.TRUE. $END

DDI用のメモリとして1 MWordを指定（MEMDDI=1）し、1コアしかないシステムでも並列MP2計算プログラムの使用を指定する（PARALL=.TRUE.）
計算の時間制限は1分 (TIMLIM=1)　指定しない場合は5,25,600分

$BASIS group

$BASIS GBASIS=N21 NGAUSS=3 $END

基底関数として 3-21Gを用いる

$GUESS group

$GUESS GUESS=HUCKEL $END

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

$DATA group

$DATA
Water...RHF/3-21G...exp.geom...R(OH)=0.95781,A(HOH)=104.4776     ; タイトル
Cnv 2      ; 対称群   回転軸の次数　この場合はC2vに相当する
     ; C1以外の場合は、通常空行を入れる
OXYGEN 8.0    ;　原子の名前（10文字以下）　原子番号　x y z は0
HYDROGEN 1.0 0.0 0.7572157 0.5865358    ; 原子の名前　原子番号　x y z
$END

3番目の原子の位置は対称性から決定される。

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

補足

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