激发态T1计算结果为负值

charles_young

最近在计算分子激发态的时候，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态计算的时候，结果一直为负

sobereva

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

charles_young

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老师，我先试一下。

az￥

charles_young 发表于 2023-12-18 22:28
现在计算结果的第4步的波函数显示是稳定的，The wavefunction is stable under the perturbations consid ...

您好，我现在也碰到了这个问题，我计算T1态是分为两步，先用 0 3 加 opt freq wb97xd/def2tzvp 优化结构，之后在这个结构上用 td(triplet,nstates=6) wb97xd/def2tzvp 自旋多重度为0 1计算了T1，结果out文件中第一个激发能是负的  我看了您这个帖子之后，我在第一步的结构基础上加了stable=opt看看波函数稳不稳定，结果是稳定的，然后在这个稳定的结构上又加opt freq进行了优化，之后再次进行了第二步的计算  但out文件中第一个激发能还是负值。能请教一下您我这个计算流程有啥问题吗？或者您后面是怎么解决了这个问题？

wal

az￥ 发表于 2025-9-10 22:16
您好，我现在也碰到了这个问题，我计算T1态是分为两步，先用 0 3 加 opt freq wb97xd/def2tzvp 优化结构 ...

说明你的结构在当前级别基态是三重态 无需解决

az￥

wal 发表于 2025-9-10 22:36
说明你的结构在当前级别基态是三重态 无需解决

那如果我要是想做一个关于SOC能级图的话，该怎么去画这个能级图呢？

wal

az￥ 发表于 2025-9-11 09:36
那如果我要是想做一个关于SOC能级图的话，该怎么去画这个能级图呢？

应该是要用DFT/ROCIS来算SOC

az￥

wal 发表于 2025-9-11 11:12
应该是要用DFT/ROCIS来算SOC

我是在ORCA上用TDDFT算SOC，但是如果要是看SOC的话，不是得先有一个能级图吗，这个能级图的话，t2态的能量是比s0的能量还要低的

wal

az￥ 发表于 2025-9-12 16:12
我是在ORCA上用TDDFT算SOC，但是如果要是看SOC的话，不是得先有一个能级图吗，这个能级图的话，t2态的能 ...

怎么了

az￥

wal 发表于 2025-9-12 16:20
怎么了

正常的不是应该t2态的能量比s0高吗

wal

az￥ 发表于 2025-9-15 20:45
正常的不是应该t2态的能量比s0高吗

你体系特殊的话有可能T2低这样