计算荧光，优化S1态，结果不收敛

liguangwu

计算一种离子化合物时，优化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

请教大神们，如何解决这个问题。

zjxitcc

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

注意激发态势能面并非总有当前分子的极小点，有可能优化过程中 断键了，说明这个分子光激发后就是会解离，你没法优化出 不解离的激发态局域极小点。另外，如果真的断键了，那就应该用UDFT和U-TDDFT了。

liguangwu

zjxitcc 发表于 2022-12-9 13:28
解决办法都在《解决SCF不收敛问题的方法》http://sobereva.com/61。在基态下解决SCF不收敛问题，然后读取基 ...

这个分子本身就是断键的离子型化合物，不加限域性条件优化时，正负离子会最终成键，加了限域性条件才能优化出一个S0结构，这种情况下优化S1应该如何操作

liguangwu

zjxitcc 发表于 2022-12-9 13:28
解决办法都在《解决SCF不收敛问题的方法》http://sobereva.com/61。在基态下解决SCF不收敛问题，然后读取基 ...

我看了优化的过程结构，并没有产生额外的断键

冰释之川

这年头还有不加色散校正的啊 小心被审稿人臭骂一顿

liguangwu

冰释之川 发表于 2022-12-9 14:03
这年头还有不加色散校正的啊 小心被审稿人臭骂一顿

什么是色散校正？我是个新手

冰释之川

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。

zjxitcc

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

该咋算咋算。遇到SCF不收敛就解决之。

liguangwu

zjxitcc 发表于 2022-12-9 13:28
解决办法都在《解决SCF不收敛问题的方法》http://sobereva.com/61。在基态下解决SCF不收敛问题，然后读取基 ...

“在基态下解决SCF不收敛问题，然后读取基态轨道进行激发态计算。” 具体如何操作，能详细点吗？

zjxitcc

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. [空行]

复制代码