[GAMESS] 久しぶりに出しました

↓はい、私の無知のせいです。

LUCY: DO NOT WORRY, CHARLIE BROWN, WE LEARN MORE
FROM LOSING THAN WE DO FROM WINNING.
CHARLIE: THEN THAT MAKES ME THE SMARTEST PERSON IN
THE WORLD!

ILLEGAL EXTENDED BASIS FUNCTION REQUESTED.

咲きました

他のつぼみも咲きかけで桜色

桜のつぼみが

最近暖かいせいか、つぼみがふくらんできました。

[GAMESS] SCF DID NOT CONVERGE...

【症状】
GAMESSでTDDFT計算していたら

" SCF DID NOT CONVERGE...NO TDDFT=EXCITE CALCULATION"

というエラーが出てABNORMALLYに計算終了。


【解決までの道のり】
検索すると"($SCF) FDIFF =.FALSE. ($END)"を加えるといいかもということが書かれていて、早速追加、再計算。たいていはこれで解決してるのだが、再び同じメッセージが。

とりあえず($SCF) DUMP =.TRUE. ($END)を加えたものを夜中に走らせて様子見。

スクラッチファイルを作らないように"($SCF) DIRSCF =.TRUE. ($END)"とパラメータ指定するとデフォルトでFDIFF=.TRUE.となるようで。これをFALSEに変えてもエラーが出るとは。
解決法はきっとどこかにあると思う。

"SCF IS UNCONVERGED, TOO MANY ITERATIONS"がとても匂う(ITERが30のところで止まっている)。

これがビンゴだとすると、SCF計算の繰り返しの最大値がデフォルトで30(MAXIT=30)なので、30回では少なすぎるから100回に増やせばいいのではと。
で、

($CONTRL) MAXIT=100 ($END)

とすると、そんなパラメータを$CONTROLグループに入れるなと門前払い。
PukiWiki - GAMESS
を見て

($MCSCF) MAXIT=100 ($END)

とするといいかもしれないと思い、トライ。門前払いは回避できたが、SCF計算は30回でストップ。logのパラメータをチェックするとMAXIT=30のまま。


もう一度$CONTROLグループにMAXIT=100を加えると門前払い。あれ？$CONTROLの行の$ENDが正常に読み込まれていない。1行の文字数オーバーかと考え、あえて記述する必要のないUNITS=ANGSを削除。再トライすると門前払いは回避、logにもMAXIT=100となっている記述が。

$CONTRL OPTIONS
---------------
SCFTYP=RHF RUNTYP=ENERGY EXETYP=RUN
MPLEVL= 0 CITYP =NONE CCTYP =NONE VBTYP =NONE
DFTTYP=B3LYP TDDFT =EXCITE
MULT = 1 ICHARG= 0 NZVAR = 0 COORD =UNIQUE
PP =NONE RELWFN=NONE LOCAL =NONE NUMGRD= F
ISPHER= -1 NOSYM = 0 MAXIT = 100 UNITS =ANGS
PLTORB= F MOLPLT= F AIMPAC= F FRIEND=
NPRINT= 7 IREST = 0 GEOM =INPUT
NORMF = 0 NORMP = 0 ITOL = 20 ICUT = 9
INTTYP=BEST GRDTYP=BEST QMTTOL= 1.0E-06

どうやら「1行の文字数オーバー」のようです(笑)。ほんまかいな。

さて、肝心のSCF計算はというと…

DIIS CONVERGED
--------------

31回で収束(ぉぃ)。無事垂直励起エネルギー計算に進みました。

【解決＆計算を通常終了することができたパラメータ】※まだグダグダです。
! TDDFT, B3LYP/6-31+G(d) (by kyoroski)
$CONTRL DFTTYP=B3LYP TDDFT=EXCITE 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
$TDDFT NSTATE=5 $END
$SCF DIRSCF =.TRUE. FDIFF =.FALSE. DAMP=.TRUE. $END

"($SCF) DAMP=.TRUE. ($END)"を除いてもいいかもしれないと思っていますが未確認。

PS
TDDFT計算をする時は、"($BASIS) DIFFSP=.TRUE."としておいた方が無難と有機化学者のための計算化学に記述されています。

ガーリック・ジョー

業務連絡: 明日まではニンニク臭いので私には近づかないように。

[GAMESS] パラメータをブログで表示できるかどうかのテスト

グダグダなパラメータでアレですが。
"＄"が表示できるかどうか。

---

! TDDFT, B3LYP/6-31+G(d)
$CONTRL DFTTYP=B3LYP TDDFT=EXCITE SCFTYP=RHF RUNTYP=ENERGY UNITS=ANGS $END
$SYSTEM TIMLIM=600000 MEMORY=80000000 $END
$BASIS GBASIS=N31 NGAUSS=6 NDFUNC=1 DIFFSP=.TRUE. $END
$GUESS GUESS=HUCKEL $END
$TDDFT NSTATE=5 $END
$SCF DIRSCF =.TRUE. FDIFF =.FALSE. $END

魚河岸まる天

沼津まで。帰宅途中に訃報、一度帰宅して帰省。