计算化学公社

标题: 求助：采用Gaussian做TD计算在输出完组态系数后卡住 [打印本页]

作者
Author: Novice 时间: 2019-7-5 10:23
标题: 求助：采用Gaussian做TD计算在输出完组态系数后卡住
本帖最后由 Novice 于 2019-7-5 13:56 编辑

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

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

作者
Author: sobereva 时间: 2019-7-5 12:51
贴具体关键词

作者
Author: Novice 时间: 2019-7-5 13:36

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

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

作者
Author: sobereva 时间: 2019-7-5 17:02
尝试CPHF=grid=fine。还不行换个基组或泛函再试

作者
Author: Novice 时间: 2019-7-5 18:07
本帖最后由 Novice 于 2019-7-5 18:21 编辑

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

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

作者
Author: sobereva 时间: 2019-7-6 17:16

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)