计算化学公社

标题: TDDFT优化激发能变成了负值 [打印本页]

作者
Author: shalene 时间: 2017-6-30 14:36
标题: TDDFT优化激发能变成了负值
g09，tddft优化第一三重态，从第三步几何开始，激发能变成了负值。还有继续优化下去的意义么？

谢谢

作者
Author: 我本是个娃娃 时间: 2017-6-30 18:19
请问你前两步做啥了？

作者
Author: 冰释之川 时间: 2017-6-30 19:57
(, 下载次数 Times of downloads: 44)

作者
Author: sobereva 时间: 2017-7-1 01:00
你优化T1直接用0 3在UDFT下优化就完了，别用TDDFT，这里提了
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314

为负说明当前参考态不是真正基态波函数，很可能S0和T1发生了交叉，在第三步的结构开始T1成了真正的基态。

作者
Author: shalene 时间: 2017-7-2 13:51

sobereva 发表于 2017-7-1 01:00
你优化T1直接用0 3在UDFT下优化就完了，别用TDDFT，这里提了
Gaussian中用TDDFT计算激发态和吸收、荧光、 ...

需要和T2比较

作者
Author: sobereva 时间: 2017-7-2 22:40

shalene 发表于 2017-7-2 13:51
需要和T2比较

T2再说T2的事，先把T1弄完了再说

作者
Author: 藍天許 时间: 2019-3-7 19:10
本帖最后由藍天許于 2019-3-7 19:15 编辑

請問如果優化T1，激發能為負值；優化T2時出現l9999.exe 的錯誤；讀進T1結構優化T3，激發能為正，這樣也代表當前的參考態不是真正的S0嗎?

---------------
#P opt td=(triplet, root=3) wb97xd/6-31G(d) int=acc2e=11 scf=(novaracc,noincfock)
---------------

作者
Author: sobereva 时间: 2019-3-8 06:01

藍天許发表于 2019-3-7 19:10
請問如果優化T1，激發能為負值；優化T2時出現l9999.exe 的錯誤；讀進T1結構優化T3，激發能為正，這樣也代表 ...

不同结构下的情况不一样，可能有的结构下默认设置就能收敛到稳定的波函数，但有的结构下可能不能。
只能说在你优化T3最后的结构下，你的参考态波函数已经是这个结构下的真实的S0了。

作者
Author: 藍天許 时间: 2019-3-9 14:41

sobereva 发表于 2019-3-8 06:01
不同结构下的情况不一样，可能有的结构下默认设置就能收敛到稳定的波函数，但有的结构下可能不能。
只能 ...

好的謝謝您

作者
Author: Jameslzy 时间: 2019-6-16 17:01

sobereva 发表于 2017-7-1 01:00
你优化T1直接用0 3在UDFT下优化就完了，别用TDDFT，这里提了
Gaussian中用TDDFT计算激发态和吸收、荧光、 ...

Sob老师好，我最近在进行激发态计算的时候也遇到了激发能变为负值的问题。想比较的是T1和S1的能量差，所以直接用了50-50的选项。关键词如下：
#p B3LYP/TZVP opt td(50-50,nstates=5) em=gd3bj
结果计算过程中，T1的激发能就变成了负的。
想问您一下，比较T1和S1的能量值我这样写关键词合理吗，是否是由于我写了opt，所以对三重态结构进行优化了，导致的三重态能量低于基态呢？

是不是先优化一个基态结构，然后用优化好的基态结构，在计算激发态能量时把opt去掉更合理一点呢？

或者是先优化一个基态结构，用stable=opt找到稳定的波函数，再把关键词改为：#p B3LYP/TZVP opt td(50-50,nstates=5) em=gd3bj guess=read 这样进行操作可以吗？

还有想问您的就是这个作业最后其实并没能正常结束，报了一个没有遇到过的l914错误：
(, 下载次数 Times of downloads: 34)
想问您一下，这个错误是什么原因导致的，怎么解决比较好？

第一次遇到比较S1和T1能量的问题，不是很会操作，辛苦老师解答一下了，谢谢Sob老师

最后我把我之前的输入文件附在后面，或许您能够用的到。
(, 下载次数 Times of downloads: 17)

作者
Author: sobereva 时间: 2019-6-17 01:22

Jameslzy 发表于 2019-6-16 17:01
Sob老师好，我最近在进行激发态计算的时候也遇到了激发能变为负值的问题。想比较的是T1和S1的能量差，所 ...

优化激发态的时候甭用50-50，否则你都不知道被优化的那个态到底是S1还是T1

作者
Author: Jameslzy 时间: 2019-6-17 09:21

sobereva 发表于 2019-6-17 01:22
优化激发态的时候甭用50-50，否则你都不知道被优化的那个态到底是S1还是T1

谢谢老师回复。那么我先优化基态结构，然后用优化好的基态结构，在计算激发态能量时把opt去掉，这时用50-50比较T1和S1的能量，这样操作可以吗？还有就是想问下老师图里的l914错误是什么原因产生的，怎么解决比较好呢？谢谢老师

作者
Author: sobereva 时间: 2019-6-17 21:56

Jameslzy 发表于 2019-6-17 09:21
谢谢老师回复。那么我先优化基态结构，然后用优化好的基态结构，在计算激发态能量时把opt去掉，这时用50- ...

如果你想比基态结构下S1和T1能量，这样可以

没完整的输入输出文件没法说

作者
Author: Jameslzy 时间: 2019-6-17 22:52

sobereva 发表于 2019-6-17 21:56
如果你想比基态结构下S1和T1能量，这样可以

没完整的输入输出文件没法说

谢谢老师回复。如果想比较S1和T1最稳定结构下S1和T1能量就要分别用singlet和triplet关键词以及opt关键词分别交两个作业来比较电子能了是吧？之前说报错的输入和输出文件我放在下面附件里了，麻烦老师帮忙看一下，谢谢老师

(, 下载次数 Times of downloads: 16)
(, 下载次数 Times of downloads: 11)

作者
Author: sobereva 时间: 2019-6-18 02:51

Jameslzy 发表于 2019-6-17 22:52
谢谢老师回复。如果想比较S1和T1最稳定结构下S1和T1能量就要分别用singlet和triplet关键词以及opt关键词 ...

是
对于T1结构尽量用UKS而非TDDFT方式优化

把50-50去了再说，还有问题把nstates增加再试

作者
Author: Jameslzy 时间: 2019-6-18 10:10

sobereva 发表于 2019-6-18 02:51
是
对于T1结构尽量用UKS而非TDDFT方式优化

好的，谢谢老师回复，我这就按照您说的操作尝试一下

作者
Author: liyuanfu7 时间: 2019-10-25 09:41

sobereva 发表于 2017-7-1 01:00
你优化T1直接用0 3在UDFT下优化就完了，别用TDDFT，这里提了
Gaussian中用TDDFT计算激发态和吸收、荧光、 ...

老师您好，我也是遇到了在TDDFT下优化T1，其能量值变为了负值，然后我按照您说的，优化T1直接用了0 3 在UDFT下进行了，输出的结果最后确实收敛到了能量为正的极小点，但是它的S**2并不严格等于2，在此情况下的T1的能量是否可用呢？输入文件及结果文件见附件

作者
Author: sobereva 时间: 2019-10-26 04:24

liyuanfu7 发表于 2019-10-25 09:41
老师您好，我也是遇到了在TDDFT下优化T1，其能量值变为了负值，然后我按照您说的，优化T1直接用了0 3 在U ...

你这个计算完全没有意义
你都没给MM区域原子设置原子电荷，根本没法体现对高层区域电子结构的影响，该怎么考虑仔细看此文
帮助设定Gaussian输入文件中优化flag和MM电荷的小工具
http://sobereva.com/332（http://bbs.keinsci.com/thread-3464-1-1.html）

Gaussian根本就不支持杂化泛函的密度拟合，写/auto毫无意义

不用管<S**2>

作者
Author: liyuanfu7 时间: 2019-10-26 10:11

sobereva 发表于 2019-10-26 04:24
你这个计算完全没有意义
你都没给MM区域原子设置原子电荷，根本没法体现对高层区域电子结构的影响，该怎么 ...

抱歉老师，输入文件传错了，请问我这样计算，得到的T态的能量可以用吗？
还有，请问如何理解这个体系有自旋污染呢

作者
Author: sobereva 时间: 2019-10-27 01:10

liyuanfu7 发表于 2019-10-26 10:11
抱歉老师，输入文件传错了，请问我这样计算，得到的T态的能量可以用吗？
还有，请问如何理解这个体系有 ...

我不知道你要具体要算什么。如果体系基态是单重态的，当前计算毫无意义

DFT计算不用管自旋污染。相关的讨论在公社论坛首页google框搜索，我之前讨论过很多次

作者
Author: liyuanfu7 时间: 2019-10-28 09:27

sobereva 发表于 2019-10-27 01:10
我不知道你要具体要算什么。如果体系基态是单重态的，当前计算毫无意义

DFT计算不用管自旋污染。相关 ...

您好，老师，我是想通过计算来比较S1和T1之间的能量关系。是的，体系基态为单重态，用0 3在UDFT下计算三重态不可以吗？为什么计算没有意义呢？
或者，我用在基态优化好的结构下，直接用50-50关键词直接输出的单三重态的能量这样应该也可以吧

作者
Author: sobereva 时间: 2019-10-29 07:49

liyuanfu7 发表于 2019-10-28 09:27
您好，老师，我是想通过计算来比较S1和T1之间的能量关系。是的，体系基态为单重态，用0 3在UDFT下计算三 ...

基态是单重态，而你写0 3且TD opt，你想优化什么？不管你想算什么，这种组合都明显毫无意义，原理上就是错的
用UDFT干嘛还写TD？

仔细看此文中的三重态计算
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）

作者
Author: liyuanfu7 时间: 2019-10-29 17:02

sobereva 发表于 2019-10-29 07:49
基态是单重态，而你写0 3且TD opt，你想优化什么？不管你想算什么，这种组合都明显毫无意义，原理上就是 ...

嗯嗯。是我记混了，谢谢老师

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