优化激发态时，总也优化不出来，优化过程中总有化学键的断裂和氢键的参与

xuejunfei1987

优化激发态时，总也优化不出来，优化过程中总有化学键的断裂和氢键的参与，希望sob老师和大家能给我提提意见。下图是基态的结构，
   
刚开始用的输入是下面的
%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.



后来修改了输入为#P td=(nstates=5,root=3) opt(calcfc,maxstep=5,notrust,maxcyc=500) uM062X/6-31+g(d) scf=(Cycles=500) int=fine，自旋多重度还是2，这次收敛了，结构是这样的，这样的结构应该不是正确的吧

我想问优化激发态的时候会有化学键的断裂生成吗，那多重度该怎么定义呢，怎样才能把这个态优化出来，这个态用nto算了是下面的转移态

这样的激发态能优化出来吗，谢谢，问题比较杂，我用的Gaussian 16

sobereva

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

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

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

xuejunfei1987

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

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

sobereva

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

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

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

xuejunfei1987

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

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

sobereva

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

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

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

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