计算化学公社

标题: 计算荧光，优化S1态，结果不收敛 [打印本页]
作者
Author: liguangwu    时间: 2022-12-9 13:22
标题: 计算荧光，优化S1态，结果不收敛
计算一种离子化合物时，优化S0时添加了限域性条件，后面在优化S1时，去掉限域性条件后，计算很难收敛，请教下有没有好的解决办法，

S0时，# opt=modredundant cam-b3lyp/6-31g(d,p) geom=connectivity
S1时，# opt td=(singlets,nstates=3,root=1) cam-b3lyp/6-31g(d,p) geom=connectivity

输出文件显示：
>>>>>>>>>> Convergence criterion notmet.
SCF Done:  E(RCAM-B3LYP) =  -4056.65654202     A.U. after 129 cycles
            NFock=128  Conv=0.14D-04     -V/T= 2.0078
Convergence failure -- runterminated.
Error termination via Lnk1e in/app1/centos6.3/gnu/apps/gaussian/g16a6/g16/l502.exe
请教大神们，如何解决这个问题。
作者
Author: zjxitcc    时间: 2022-12-9 13:28
本帖最后由 zjxitcc 于 2022-12-9 13:31 编辑

解决办法都在《解决SCF不收敛问题的方法》http://sobereva.com/61。在基态下解决SCF不收敛问题，然后读取基态轨道进行激发态计算。

注意激发态势能面并非总有当前分子的极小点，有可能优化过程中 断键了，说明这个分子光激发后就是会解离，你没法优化出 不解离的激发态局域极小点。另外，如果真的断键了，那就应该用UDFT和U-TDDFT了。
作者
Author: liguangwu    时间: 2022-12-9 13:48
zjxitcc 发表于 2022-12-9 13:28
解决办法都在《解决SCF不收敛问题的方法》http://sobereva.com/61。在基态下解决SCF不收敛问题，然后读取基 ...

这个分子本身就是断键的离子型化合物，不加限域性条件优化时，正负离子会最终成键，加了限域性条件才能优化出一个S0结构，这种情况下优化S1应该如何操作
作者
Author: liguangwu    时间: 2022-12-9 13:51
zjxitcc 发表于 2022-12-9 13:28
解决办法都在《解决SCF不收敛问题的方法》http://sobereva.com/61。在基态下解决SCF不收敛问题，然后读取基 ...

我看了优化的过程结构，并没有产生额外的断键
作者
Author: 冰释之川    时间: 2022-12-9 14:03
这年头还有不加色散校正的啊 小心被审稿人臭骂一顿
作者
Author: liguangwu    时间: 2022-12-9 14:05
冰释之川 发表于 2022-12-9 14:03
这年头还有不加色散校正的啊 小心被审稿人臭骂一顿

什么是色散校正？我是个新手
作者
Author: 冰释之川    时间: 2022-12-9 14:14
liguangwu 发表于 2022-12-9 14:05
什么是色散校正？我是个新手

《谈谈“计算时是否需要加DFT-D3色散校正？”》http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）
对于主要只体现化学键作用的体系，加D3是没必要的，虽然加了也没坏处。诸如优化有机小分子、简单配合物（如二茂铁）、简单阴阳离子体系（比如乙酸钠、硫酸钙）等等，加D3基本不会对结果产生任何影响。一方面是D3里面有阻尼函数，阻尼函数会使得校正函数在近程区域衰减为0，从而避免影响原先DFT泛函就已经能描述良好的范畴，另一方面是对于这些较强的相互作用，微小的色散作用的影响本身就是可以忽略的。但如果是有机大分子、配体比较大的配合物（此时配体间容易有弱相互作用）、比较大的阴阳离子对等等情况，弱相互作用的存在不可忽略，当使用对色散作用描述烂的泛函的时候建议加上D3。


作者
Author: zjxitcc    时间: 2022-12-9 14:17
liguangwu 发表于 2022-12-9 13:48
这个分子本身就是断键的离子型化合物，不加限域性条件优化时，正负离子会最终成键，加了限域性条件才能优 ...

该咋算咋算。遇到SCF不收敛就解决之。
作者
Author: liguangwu    时间: 2022-12-9 14:53
zjxitcc 发表于 2022-12-9 13:28
解决办法都在《解决SCF不收敛问题的方法》http://sobereva.com/61。在基态下解决SCF不收敛问题，然后读取基 ...

“在基态下解决SCF不收敛问题，然后读取基态轨道进行激发态计算。” 具体如何操作，能详细点吗？
作者
Author: zjxitcc    时间: 2022-12-9 15:38
liguangwu 发表于 2022-12-9 14:53
“在基态下解决SCF不收敛问题，然后读取基态轨道进行激发态计算。” 具体如何操作，能详细点吗？

一种合理示例（但不仅限于这一种）:
  1. %mem=48GB
  2. %nprocshared=24
  3. %chk=xxx.chk
  4. #p RCAM-B3LYP/6-31G(d,p) scf(maxcycle=128,xqc)
  5. [空行]
  6. 标题行
  7. [空行]
  8. 电荷 自旋多重度
  9. 坐标
  10. [空行]
  11. --Link--
  12. %mem=48GB
  13. %nprocshared=24
  14. %chk=xxx.chk
  15. #p opt freq TD RCAM-B3LYP chkbasis scf(maxcycle=256) guess=read geom=allcheck
  16. [空行]
复制代码





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3