计算化学公社

标题: 第一激发态振子强度很小或为负 [打印本页]

作者
Author: yangjunfang 时间: 2021-4-7 10:17
标题: 第一激发态振子强度很小或为负
求助：您好，我想问一下我激发态优化正常结束，但是第一激发态振子强度很小或为负，原因是什么，有什么办法解决吗？

作者
Author: 小范范1989 时间: 2021-4-7 13:27
换高HF的泛函试试

作者
Author: zjxitcc 时间: 2021-4-7 14:09
很小可以。为负是什么样的？-0.0001这种？理论上是不允许为负的

作者
Author: 冰释之川 时间: 2021-4-7 14:17
(, 下载次数 Times of downloads: 49)

作者
Author: niobium 时间: 2021-4-7 15:51

冰释之川发表于 2021-4-7 14:17

楼主说的是振子强度

作者
Author: 喵星大佬 时间: 2021-4-8 00:38

niobium 发表于 2021-4-7 15:51
楼主说的是振子强度

振子强度还能为负的嘛？

作者
Author: niobium 时间: 2021-4-8 09:02

喵星大佬发表于 2021-4-8 00:38
振子强度还能为负的嘛？

我错了。振子强度的符号应该和激发能是一样的。所以4楼的回复是正确的。

作者
Author: yangjunfang 时间: 2021-4-8 10:00
好的，明白了，谢谢各位，我先试试

作者
Author: 冰释之川 时间: 2021-4-8 13:17

niobium 发表于 2021-4-8 09:02
我错了。振子强度的符号应该和激发能是一样的。所以4楼的回复是正确的。

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

作者
Author: niobium 时间: 2021-4-8 15:02

冰释之川发表于 2021-4-8 13:17

谢谢

作者
Author: yangjunfang 时间: 2021-4-10 19:08
按照之前各位所说的先找了真正的基态波函数，然后写guess=read进行激发态计算，发现在激发态优化的过程中能量越来越大，这是什么情况啊

作者
Author: pyyu 时间: 2021-4-13 17:43
正常，激发态能量本来就比基态高啊，只要到后面收敛就行了

作者
Author: yangjunfang 时间: 2021-5-8 15:31
您好，我按照您说的先稳定了波函数，然后再算激发态，也加了关键词guess=read，但是结果还是激发能为负，SCF不收敛，请问这是什么原因啊

作者
Author: wzkchem5 时间: 2021-5-8 16:14

yangjunfang 发表于 2021-5-8 08:31
您好，我按照您说的先稳定了波函数，然后再算激发态，也加了关键词guess=read，但是结果还是激发能为负，SC ...

你确定TDDFT计算和stable计算的理论级别严格一致吗？如果严格一致，应该迭代一步就收敛了，即使考虑到数值误差，3、4步也收敛了，不可能不收敛。

作者
Author: yangjunfang 时间: 2021-5-10 19:16
本帖最后由 yangjunfang 于 2021-5-10 19:18 编辑

wzkchem5 发表于 2021-5-8 16:14
你确定TDDFT计算和stable计算的理论级别严格一致吗？如果严格一致，应该迭代一步就收敛了，即使考虑到数 ...

理论级别是一致的，有可能是我的关键词有问题吗，输入guess=read时，%chk对于chk文件的命名是基态优化之后chk文件的名称吗，还是得重新命名

作者
Author: wzkchem5 时间: 2021-5-10 19:40

yangjunfang 发表于 2021-5-10 12:16
理论级别是一致的，有可能是我的关键词有问题吗，输入guess=read时，%chk对于chk文件的命名是基态优化之 ...

stable计算的分子结构是激发态结构优化收敛后的结构吗？

作者
Author: yangjunfang 时间: 2021-5-10 19:41

wzkchem5 发表于 2021-5-10 19:40
stable计算的分子结构是激发态结构优化收敛后的结构吗？

是基态的

作者
Author: wzkchem5 时间: 2021-5-10 19:45

yangjunfang 发表于 2021-5-10 12:41
是基态的

所以你到底要检验基态结构下的稳定性，还是激发态结构下的稳定性？
如果你的chk文件是激发态结构优化的chk文件的话，里面的波函数是结构优化最后一步的波函数，不是基态结构下的波函数。

作者
Author: yangjunfang 时间: 2021-5-10 19:49
本帖最后由 yangjunfang 于 2021-5-10 20:12 编辑

wzkchem5 发表于 2021-5-10 19:45
所以你到底要检验基态结构下的稳定性，还是激发态结构下的稳定性？
如果你的chk文件是激发态结构优化的c ...

我是激发态优化结束后，激发能为负数，那这个应该是稳定基态波函数还是激发态波函数呢？我一直以为是要稳定基态波函数，然后再在这个基础上读取基态的chk文件来重新优化激发态，不知道这样做对吗？还是应该怎么做？

作者
Author: wzkchem5 时间: 2021-5-10 21:43

yangjunfang 发表于 2021-5-10 12:49
我是激发态优化结束后，激发能为负数，那这个应该是稳定基态波函数还是激发态波函数呢？我一直以为是要稳 ...

是要在激发态结构下找到稳定的基态波函数。
逻辑是这样的：你的目的是计算激发态结构下的激发态，但是只有在同一个结构下的基态波函数稳定的时候，这个结构下的激发态计算结果才有意义。所以要保证激发态结构下的基态波函数是稳定的。基态结构下的基态波函数是否稳定，和激发态结构下的结果是否有意义毫无关系。

作者
Author: yangjunfang 时间: 2021-5-11 10:20
您好，我还想问一下，怎么在激发态结构下找到稳定的基态波函数呢？
是将现在的激发态结构取出来，然后加关键词stable=opt重新优化，再在这个基础上重新进行激发态优化吗

作者
Author: wzkchem5 时间: 2021-5-11 15:00

yangjunfang 发表于 2021-5-11 03:20
您好，我还想问一下，怎么在激发态结构下找到稳定的基态波函数呢？
是将现在的激发态结构取出来，然后加关 ...

对，如果重新优化结构后结构有了比较大的差别，那最好再stable=opt一次

作者
Author: yangjunfang 时间: 2021-5-11 15:22
好的，谢谢您，我试试

作者
Author: yangjunfang 时间: 2021-5-12 09:31
老师您好，我按照您说的现在原来激发态结构的基础上稳定了对应的基态波函数，也检查过结果了，显示基态波函数已经稳定，然后加关键词guess=read读取当前的结构重新进行激发态优化，但是在out文件中前几步的激发能还是正值，可是到后面几步的时候激发能又成为负值了（还没优化完），这是什么原因？

作者
Author: sobereva 时间: 2021-5-18 01:28

yangjunfang 发表于 2021-5-12 09:31
老师您好，我按照您说的现在原来激发态结构的基础上稳定了对应的基态波函数，也检查过结果了，显示基态波函 ...

把步长上限减小再试

或者用后来激发能已经变负的结构重新用stable=opt得到稳定波函数，再继续优化

作者
Author: yangjunfang 时间: 2021-5-19 13:36

sobereva 发表于 2021-5-18 01:28
把步长上限减小再试

或者用后来激发能已经变负的结构重新用stable=opt得到稳定波函数，再继续优化

嗯嗯，好的，我试试

作者
Author: yangjunfang 时间: 2021-5-20 20:51

wzkchem5 发表于 2021-5-10 19:40
stable计算的分子结构是激发态结构优化收敛后的结构吗？

是的

作者
Author: yangjunfang 时间: 2021-5-21 09:07

yangjunfang 发表于 2021-5-20 20:51
是的

就很奇怪一直优化不出来，不知道哪里出问题了

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