计算化学公社

标题: 激发态T1计算结果为负值 [打印本页]

作者
Author: charles_young 时间: 2023-12-18 18:26
标题: 激发态T1计算结果为负值
最近在计算分子激发态的时候，T1态一直是负值，看了sob老师的帖子，重新检查了自旋多重度并且stable=opt之后，仍然没有解决，我的计算步骤如下：
1、#P B3LYP/6-311G(d,p) opt freq
2、#P b3lyp/6-311g(d,p) stable=opt guess=read geom=check
3、#P b3lyp/6-311g(d,p) opt freq guess=read geom=check
4、#P b3lyp/6-311g(d,p) stable guess=read geom=check
5、#P cam-b3lyp/6-311+g(d,p) td(nstate=10) guess=read geom=check
6、#P cam-b3lyp/6-311+g(d,p) td(triplet,nstate=10) guess=read geom=check
程序都是正常结束，就在第6步T1态计算的时候，结果一直为负
(, 下载次数 Times of downloads: 9)

作者
Author: sobereva 时间: 2023-12-18 20:59
确保第4步提示了波函数稳定。然后第5步使用和第4步完全相同的级别计算，并读第4步的chk当初猜
既然后来改成了cam-b3lyp/6-311+g(d,p)，自然得保证此级别的基态波函数是稳定的，然而你又没在此级别测试
就算stable任务显示波函数不稳定，程序也不会报错，显然得自己看是否提示稳定

作者
Author: charles_young 时间: 2023-12-18 22:28

sobereva 发表于 2023-12-18 20:59
确保第4步提示了波函数稳定。然后第5步使用和第4步完全相同的级别计算，并读第4步的chk当初猜
既然后来改 ...

现在计算结果的第4步的波函数显示是稳定的，The wavefunction is stable under the perturbations considered.
所以我现在是用cam-b3lyp/6-311+g(d,p)重新算一下stable看是否稳定；
或者用之前的b3lyp/6-311g(d,p)算激发态吗？
谢谢sob老师，我先试一下。

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