计算化学公社

标题: 优化激发态时，总也优化不出来，优化过程中总有化学键的断裂和氢键的参与 [打印本页]

作者
Author: xuejunfei1987 时间: 2019-10-19 16:39
标题: 优化激发态时，总也优化不出来，优化过程中总有化学键的断裂和氢键的参与
优化激发态时，总也优化不出来，优化过程中总有化学键的断裂和氢键的参与，希望sob老师和大家能给我提提意见。下图是基态的结构，
(, 下载次数 Times of downloads: 17)
刚开始用的输入是下面的
%NPROCSHARED=14
%mem=50GB
%chk=GGT-M06.chk
#P td=(nstates=5,root=3) opt(calcfc,maxcyc=500) uM062X/6-31+g(d) scf=(maxcycles=500)
自旋多重度是设的2，因为初始是一个自由基，但是总是收敛失败，结构是下面这样的：

>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UCAM-B3LYP) =  -307.959188481    A.U. after  501 cycles
         NFock=500  Conv=0.13D-02    -V/T= 2.0084
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7560 S= 0.5030
<L.S>= 0.000000000000E+00
KE= 3.053836119891D+02 PE=-1.221145931873D+03 EE= 3.567684123386D+02
Annihilation of the first spin contaminant:
S**2 before annihilation    0.7560, after    0.7500
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/lzu/g16/l502.exe at Sat Oct 19 06:49:20 2019.

(, 下载次数 Times of downloads: 23)

后来修改了输入为#P td=(nstates=5,root=3) opt(calcfc,maxstep=5,notrust,maxcyc=500) uM062X/6-31+g(d) scf=(Cycles=500) int=fine，自旋多重度还是2，这次收敛了，结构是这样的，这样的结构应该不是正确的吧
(, 下载次数 Times of downloads: 25)
我想问优化激发态的时候会有化学键的断裂生成吗，那多重度该怎么定义呢，怎样才能把这个态优化出来，这个态用nto算了是下面的转移态
(, 下载次数 Times of downloads: 34)
这样的激发态能优化出来吗，谢谢，问题比较杂，我用的Gaussian 16

作者
Author: sobereva 时间: 2019-10-20 00:12
写scf=(Cycles=500)毫无意义，仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

本来扭曲的结构SCF就难收敛，还用对当前体系计算根本没意义的弥散函数，是自讨苦吃。

写int=fine更是雪上加霜，本来明尼苏达系列泛函比起其它泛函就更有可能因为g09默认的fine的积分格点精度不够而难收敛，作为g16用户还不如用默认的ultrafine。

激发态出现成键方式变化没有什么特别的。

作者
Author: xuejunfei1987 时间: 2019-10-20 12:49

sobereva 发表于 2019-10-20 00:12
写scf=(Cycles=500)毫无意义，仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（h ...

谢谢老师，那我该怎么优化这个态啊那个优化出来的结构明显不对啊，频率分析出好几个虚频。
您能给一些具体的建议吗，比如泛函，基组，算法之类的，怎样才能有一个好的效果。激发态成键方式变化，那多重度不用管吗，我现在的进展就卡在了这儿，始终优化不出来这个态，这个态对我很重要。

作者
Author: sobereva 时间: 2019-10-20 19:46

xuejunfei1987 发表于 2019-10-20 12:49
谢谢老师，那我该怎么优化这个态啊那个优化出来的结构明显不对啊，频率分析出好几个虚频。
您能给一些 ...

Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

基组把弥散函数去了就完了
用PBE0

作者
Author: xuejunfei1987 时间: 2019-10-21 15:51

sobereva 发表于 2019-10-20 19:46
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.ke ...

sob老师，我按照您的建议做了一下激发态优化，还是收敛不了，这个转移态能用cam-b3lyp吗，我也试过这个，也收敛不了，是不是存在势能面交叉呢，如果我想做个势能面扫描，应该扫哪个键合适

作者
Author: sobereva 时间: 2019-10-21 22:08

xuejunfei1987 发表于 2019-10-21 15:51
sob老师，我按照您的建议做了一下激发态优化，还是收敛不了，这个转移态能用cam-b3lyp吗，我也试过这个 ...

仔细看此文，结合实际输出进行分析

谈谈势能面交叉对激发态优化的影响
http://sobereva.com/468（http://bbs.keinsci.com/thread-12318-1-1.html）

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