利用高斯进行激发态优化

yinjie

利用高斯进行S1态优化时，收敛标准中前两项都为yes，后两项关于位移的一直不收敛，并且最初始结构的能量竟然是最低的，关于优化时能量的变化放在了图片附件里，请问大家，这样的情况还有必要算下去嘛，要怎样做才会收敛呢？
输入文件如下：
%nprocshared=8
%mem=8GB
%chk=s1.chk
# opt td=nstates=6 cam-b3lyp/6-31+g(d,p) scrf=(solvent=acetonitrile) scf(novaracc,noincfock)int=acc2e=12

s1 opt

0 1
C                  4.68530900    0.48300900   -1.24204000
C                  3.29065900    0.27075000   -1.23293600
C                  2.73721900   -1.04635400   -1.44292300
C                  4.91175000   -1.83072200   -1.69532100
N                  3.57969300   -2.05566500   -1.66712500
N                  5.49807900   -0.58469900   -1.48901100
N                  4.96541000    1.74938300   -0.93619800
N                  3.73859600    2.26998100   -0.75118500
N                  2.71601700    1.45062500   -0.91478600
.....

.....

我本是个娃娃

第一，如果没有特殊情况，一定不要用弥散函数，乱用弥散函数有害无益
第二，可以尝试先不在溶剂情况下优化，气相优化收敛后，在添加溶剂环境优化

liyuanhe211

GaussView 5 对 TD 的 Result-Optimization 的 Total Energy 读的是SCF能量，不是root激发态能量。所以不能看它的曲线，应该自己看输出文件、作图、探究优化趋势。

ggdh

liyuanhe说的对，另外加个opt(maxstep=2)接着算 很有肯能会收敛。

yinjie

我本是个娃娃 发表于 2017-6-21 17:40
第一，如果没有特殊情况，一定不要用弥散函数，乱用弥散函数有害无益
第二，可以尝试先不在溶剂情况下优化 ...

你好，谢谢回复~
因为前面的计算中需要弥散，在这里为了保持一致性，所以加了弥散；对于第二点会进行尝试一下~~多谢指点

yinjie

liyuanhe211 发表于 2017-6-21 19:45
GaussView 5 对 TD 的 Result-Optimization 的 Total Energy 读的是SCF能量，不是root激发态能量。所以不能 ...

你好，谢谢回复~
想问一下，是要看输出文件中的 Total Energy, E(TD-HF/TD-DFT) =  -2952.93659353 这一部分的能量，之后根据这个能量作图，看优化趋势嘛，不知道这么理解对不对，希望帮忙指点一下，谢谢了~~

yinjie

ggdh 发表于 2017-6-21 20:17
liyuanhe说的对，另外加个opt(maxstep=2)接着算 很有肯能会收敛。

你好~
请问是要提取出现在计算的构型，缩小步长，继续进行优化么~~

ggdh

yinjie 发表于 2017-6-22 08:27
你好~
请问是要提取出现在计算的构型，缩小步长，继续进行优化么~~

是的，前提是要用g09的最新版（D或者E)，不然可能出现无法follow某一个态而导致的振荡

yinjie

ggdh 发表于 2017-6-22 09:11
是的，前提是要用g09的最新版（D或者E)，不然可能出现无法follow某一个态而导致的振荡

好的，谢谢指点~~我会试一下~~

yinjie

ggdh 发表于 2017-6-22 09:11
是的，前提是要用g09的最新版（D或者E)，不然可能出现无法follow某一个态而导致的振荡

好的，谢谢指点~~我会试一下的~~

yinjie

ggdh 发表于 2017-6-21 20:17
liyuanhe说的对，另外加个opt(maxstep=2)接着算 很有肯能会收敛。

加入opt(maxstep=2)后优化后，问题得到了解决，谢谢帮忙~~