计算化学公社

标题: Gaussian某些过渡金属体系读取chk文件重算时，没有在一步之内收敛的问题 [打印本页]

作者
Author: yueyue1765 时间: 2023-11-16 15:08
标题: Gaussian某些过渡金属体系读取chk文件重算时，没有在一步之内收敛的问题
各位老师好！
我在用HF/def2-tzvp计算过渡金属配合物体系，先计算单点，并保存了scf收敛的波函数的chk文件。用guess=read geom=check关键词读取波函数进行后续计算时，有一小部分体系并没有在1步内达到收敛，有些甚至到了maxcycle都没有收敛。

附件给出了两次计算的输入卡和log文件，我注意到续算的log文件中，读取初猜波函数得到的第一步能量和初算得到的单点能是相同的，但由于MaxDP=1.56D-5，没达到收敛限，所以进入了下一步迭代，共22步才收敛。虽然这个对计算结果基本没有影响，但想知道这里不收敛的原因是什么呢？

作者
Author: zjxitcc 时间: 2023-11-16 17:01
尝试添加scf(NoIncFock,NoVarAcc)

作者
Author: yueyue1765 时间: 2023-11-19 01:07

zjxitcc 发表于 2023-11-16 17:01
尝试添加scf(NoIncFock,NoVarAcc)

谢谢老师的回复！我的体系并不是没收敛，而是没有在一步之内收敛，我有点好奇原因。刚刚试了添加这两个关键词，仍然是22步收敛。
之前看sobereva老师的博文，这两个关键词应该是用来解决scf不收敛的问题的，好像和我的问题不太一样（？）

作者
Author: hebrewsnabla 时间: 2023-11-19 01:28
本帖最后由 hebrewsnabla 于 2023-11-19 01:32 编辑

2L建议这两个关键词的原因是，VarAcc可以在计算前期降低精度来加快计算，但降低精度就会导致原来收敛的波函数略微偏离收敛标准。NoVarAcc可以关闭之。NoIncFock与之类似但是不太一样，我觉得并不一定有用。

既然排除了这两个原因，还有一种可能的原因是，原本的计算是默认采用了DIIS的，而DIIS会导致最后一步的波函数不只取决于上一步的结果，而是前面若干步的结果。而重新读取波函数续算时，没有前面若干步的信息，就会导致第一步（在有些程序里是第二步）的MaxDP略微偏离0。从这个角度来说如果程序采用的收敛标准是轨道梯度的话就不会有这个问题。

这个问题在做结构优化的时候比较常见，SCF我不确定是否也会有。
如果确实是由这种原因导致的，似乎没有什么解决方法。

作者
Author: zjxitcc 时间: 2023-11-19 12:28
本帖最后由 zjxitcc 于 2023-11-19 12:31 编辑

yueyue1765 发表于 2023-11-19 01:07
谢谢老师的回复！我的体系并不是没收敛，而是没有在一步之内收敛，我有点好奇原因。刚刚试了添加这两个关 ...

它们能帮助解决SCF不收敛问题，但并不局限于此唯一用途，你需要知其然and知其所以然。还有几种测试手段没试：
（1）两个任务都加nosymm
（2）两个任务都加nosymm int=nobasistransform
（3）第二个任务加chkbasis nosymm int=nobasistransform，此时不要写“/def2TZVP”。
然后来说一下各自需要几圈SCF。如果都需要22圈，那原因很可能是4L所说。

作者
Author: yueyue1765 时间: 2023-11-20 18:17

zjxitcc 发表于 2023-11-19 12:28
它们能帮助解决SCF不收敛问题，但并不局限于此唯一用途，你需要知其然and知其所以然。还有几种测试手段没 ...

我按照5L的方法进行了尝试，无论加什么关键词，初算都是29步收敛，读chk文件续算都是22步收敛。打算暂时不管这个问题了，因为对结果影响确实不大。
非常感谢两位老师的回复！@hebrewsnabla @zjxitcc

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