计算化学公社

标题: 关于激发态计算时density关键词使用求助 [打印本页]

作者
Author: zhb 时间: 2021-9-10 13:02
标题: 关于激发态计算时density关键词使用求助
老师您好，我做计算时遇到一个问题，我们知道根据博文《Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法》计算吸收光谱，并且博文里讲述如何产生激发态波函数算电荷，可我按照文章做发现加关键词density后跃迁能和最大吸收峰值都变了，差距很大，那么吸收光谱输入文件用，加关键词density的还是不加density的？（不加关键词density的吸收峰和跃迁能更接近实验值，实验值：跃迁能3.39ev吸收峰367nm）谢谢老师

作者
Author: zjxitcc 时间: 2021-9-10 14:18
两个文件里SCF Done那里的能量是否一样？

作者
Author: zhb 时间: 2021-9-10 14:52

zjxitcc 发表于 2021-9-10 14:18
两个文件里SCF Done那里的能量是否一样？

老师，SCF Done的能量是一样的。是什么原因呢？可以用哪个输入文件。谢谢老师

作者
Author: paramecium86 时间: 2021-9-10 16:07
感觉是不是SCRF的问题。用了density 就让计算从非平衡溶剂变成平衡溶剂了。

作者
Author: zhb 时间: 2021-9-10 16:22

paramecium86 发表于 2021-9-10 16:07
感觉是不是SCRF的问题。用了density 就让计算从非平衡溶剂变成平衡溶剂了。

谢谢，那我计算吸收光谱的话应该用density文件做呢?还是不含关键词的做呢？峰值和能量也对不上了，小白有些懵，请指教

作者
Author: sobereva 时间: 2021-9-11 03:07
一般就两个可能原因
(1)用density时TDDFT的Davidson迭代收敛限会更严，因此会影响激发态计算结果。但通常影响不会大到导致光谱图都能看出显著差异
(2)density影响了溶剂模型的效果，包括对跃迁的响应和溶剂慢部分的考虑
不带density的时候在TD里写上conver=6看看结果，这时候Davidson迭代收敛限将和用density的情况一样。如果这和之前带density的时候相差很大但和不用density时相差很小，就不要用带density时候的数据绘制光谱。

作者
Author: zhb 时间: 2021-9-11 10:33

sobereva 发表于 2021-9-11 03:07
一般就两个可能原因
(1)用density时TDDFT的Davidson迭代收敛限会更严，因此会影响激发态计算结果。但通常 ...

老师，如果带density关键词和不带density的差距很大，我可以用不带density的文件做光谱，然后用带density做激发态ADCH电荷，这样可以吗？因为根据您的博文产生激发态波函数.wfn文件我只能用.fchk文件通过Multiwfn转换产生，我们是送往节点计算的，不能直接产.wfn。如果不可以的话，有什么办法可以用一个文件既能做吸收光谱，又能算激发态ADCH电荷。麻烦老师啦

作者
Author: sobereva 时间: 2021-9-12 05:25

zhb 发表于 2021-9-11 10:33
老师，如果带density关键词和不带density的差距很大，我可以用不带density的文件做光谱，然后用带density ...

可以这样

不直接产生wfn文件，就依靠Multiwfn产生记录了激发态自然轨道的fch文件做激发态波函数分析，下文有例子
在Multiwfn中基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例
http://sobereva.com/403（http://bbs.keinsci.com/thread-9001-1-1.html）

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