求助：采用Gaussian做TD计算在输出完组态系数后卡住

Novice

Gaussian 09E进行# td=(50-50,nstates=13) m062x/tzvp maxdisk=100GB iop(9/40=3) density 计算，在输出完各跃迁的组态系数之后
...........
    225 <- 245       -0.00206
SavETr:  write IOETrn=   770 NScale= 10 NData=  16 NLR=1 NState=   26 LETran=     478.
          Differentiating once with respect to electric field.
                with respect to dipole field.
          Electric field/nuclear overlap derivatives assumed to be zero.
          There are     1 degrees of freedom in the 1st order CPHF.  IDoFFX=0 NUNeed=     1.
LinEq1:  Iter=  0 NonCon=     1 RMS=2.38D-04 Max=3.90D-02 NDo=     1
AX will form     1 AO Fock derivatives at one time.
LinEq1:  Iter=  1 NonCon=     1 RMS=2.31D-04 Max=3.24D-02 NDo=     1
...........
LinEq1:  Iter=888 NonCon=     1 RMS=2.47D-10 Max=4.36D-09 NDo=     1

如上，本来应该较快完成的任务，在 LinEq1后 Iter=888次，这是为什么？再继续等下去是否已经没有意义，是不是应改停了任务重新就算了？
------------

补充，一个本该数个小时左右完成的任务，在进行了二十多个小时后报错结束了
LinEq1:  Iter=999 NonCon=     1 RMS=2.24D-10 Max=4.71D-09 NDo=     1
LinEq1:  Iter=*** NonCon=     1 RMS=2.37D-10 Max=4.67D-09 NDo=     1
CPHF failed to converge in LinEq1.
Error termination via Lnk1e in /g16/l1002.exe at Fri Jul  5 11:18:01 2019.
Job cpu time:       9 days  5 hours  5 minutes 14.8 seconds.
Elapsed time:       0 days 18 hours 39 minutes 46.1 seconds.
File lengths (MBytes):  RWF=  25119 Int=      0 D2E=      0 Chk=    713 Scr=      1

个人感觉应该是积分精度不够。欢迎指正。

sobereva

贴具体关键词

Novice

sobereva 发表于 2019-7-5 12:51
贴具体关键词

已补充，并增加了最新的报错描述。

sobereva

尝试CPHF=grid=fine。还不行换个基组或泛函再试

Novice

sobereva 发表于 2019-7-5 17:02
尝试CPHF=grid=fine。还不行换个基组或泛函再试

多谢社长。有个疑问，这一步是不是和加入了density关键词有关呢？此步的作用是什么？

sobereva

Novice 发表于 2019-7-5 18:07
多谢社长。有个疑问，这一步是不是和加入了density关键词有关呢？此步的作用是什么？

这就是density关键词导致的
为了得到激发态弛豫密度，需要迭代求解CPHF方程，恰好当前这一步难收敛