计算化学公社

标题: 利用高斯进行激发态优化 [打印本页]

作者
Author: yinjie 时间: 2017-6-21 17:18
标题: 利用高斯进行激发态优化
利用高斯进行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
.....

.....

作者
Author: 我本是个娃娃 时间: 2017-6-21 17:40
第一，如果没有特殊情况，一定不要用弥散函数，乱用弥散函数有害无益
第二，可以尝试先不在溶剂情况下优化，气相优化收敛后，在添加溶剂环境优化

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

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

作者
Author: yinjie 时间: 2017-6-21 20:45

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

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

作者
Author: yinjie 时间: 2017-6-21 20:50

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

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

作者
Author: yinjie 时间: 2017-6-22 08:27

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

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

作者
Author: ggdh 时间: 2017-6-22 09:11
本帖最后由 ggdh 于 2017-6-22 09:12 编辑

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

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

作者
Author: yinjie 时间: 2017-6-22 09:22

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

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

作者
Author: yinjie 时间: 2017-6-22 09:22

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

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

作者
Author: yinjie 时间: 2017-6-27 15:11

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

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

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