计算化学公社

标题: 请问在计算激发态时，对应的波长出现负值，这是正常现象吗？ [打印本页]

作者
Author: wangxiaofeng 时间: 2022-11-21 11:53
标题: 请问在计算激发态时，对应的波长出现负值，这是正常现象吗？
输入文件所用关键词和输出文件中出现负值的截图如下

作者
Author: zjxitcc 时间: 2022-11-21 13:05
本帖最后由 zjxitcc 于 2022-11-21 13:15 编辑

这不是正常现象，这表示当前基态UDFT波函数不稳定。scf=qc大多数时候不是明智的选择（尤其是不知道qc为何意、这么写有什么作用时），即使用了scf=qc，也一定要检测波函数稳定性。合理的计算步骤：
第一步，先做基态UDFT单点计算，注意加上关键词stable=opt nosymm检验波函数稳定性、确保波函数稳定；
第二步，读取上一步的轨道，做TDDFT单点计算/结构优化。
（如果需要的话）第三步，读取优化好的激发态极小点结构、基态UDFT轨道，算吸收光谱，nstates可以设置较大。

当前关键词含义较为混乱，难以知晓是想优化第一激发态结构、还是想算吸收光谱（亦或是其他目的）。优化第一激发态的几何结构压根不需要nstates=100，这会浪费大量电费。默认的nstates=3一般就够，（如果你懂检验够不够的话）即使不够，写nstates=5更保险，与100个态计算量有天壤之别。

作者
Author: sobereva 时间: 2022-11-21 14:01
别轻易用SCF=QC，此文明确强调了
解决SCF不收敛问题的方法
http://sobereva.com/61

不知道你算什么体系，如果就5个原子，算100个态毫无意义。而且较高的态基本都是里德堡激发，靠这种泛函和不带弥散的基组根本没谱。而且这么点体系应当用EOM-CCSD等更好级别来算。

作者
Author: wangxiaofeng 时间: 2022-11-21 21:39

zjxitcc 发表于 2022-11-21 13:05
这不是正常现象，这表示当前基态UDFT波函数不稳定。scf=qc大多数时候不是明智的选择（尤其是不知道qc为何意 ...

您好！我是想优化一个第一激发单重态，除了去掉scf=qc,改一下nstates的个数，请问其他关键词还有什么问题吗？

作者
Author: wzkchem5 时间: 2022-11-21 22:08

wangxiaofeng 发表于 2022-11-21 14:39
您好！我是想优化一个第一激发单重态，除了去掉scf=qc,改一下nstates的个数，请问其他关键词还有什么问题 ...

stable=opt是最重要的。相比之下去掉scf=qc和改nstates反倒次要一点

作者
Author: wangxiaofeng 时间: 2022-11-22 10:43

wzkchem5 发表于 2022-11-21 22:08
stable=opt是最重要的。相比之下去掉scf=qc和改nstates反倒次要一点

好的，谢谢您

作者
Author: wangxiaofeng 时间: 2022-11-22 11:11

sobereva 发表于 2022-11-21 14:01
别轻易用SCF=QC，此文明确强调了
解决SCF不收敛问题的方法
http://sobereva.com/61

通过百度网盘分享的文件：jitai-ji…
链接:https://pan.baidu.com/s/1rOH-YJrCRaTYj-5C84913w
提取码:hh2u
复制这段内容打开「百度网盘APP 即可获取」

不是五个原子体系，只是截了一部分，完整的输出文件已上传，请老师解答。

作者
Author: zjxitcc 时间: 2022-11-22 11:18

wangxiaofeng 发表于 2022-11-22 11:11
通过百度网盘分享的文件：jitai-ji…
链接:https://pan.baidu.com/s/1rOH-YJrCRaTYj-5C84913w
提取码:hh ...

2楼，3楼，5楼已经有很多解答了。应该根据上面的解答自己修改后，若有问题才再次提问

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