chengmeiling 发表于 2025-6-15 15:53
老师好，我在按照上面步骤做完之后波长还是负值，想问一下之后该怎么解决？

参考态波函数还是不稳定
可能stable=opt确认得到稳定波函数后，TDDFT的时候没读其波函数。没细节没法说

sobereva 发表于 2023-1-7 08:14
(2.1)里对0 3做stable=opt没意义
而且体系本身基态是单重态闭壳层的，(3)里没必要写U，否则计算白白增加很 ...

老师好，我在按照上面步骤做完之后波长还是负值，想问一下之后该怎么解决？

乐平 发表于 2023-1-7 13:45
谢谢 Sob 老师的回复

您的意思是，在 (2.1)步的时候，用 (1.2)步里 T1 的优化后的结构，用 stable=opt ...

对

sobereva 发表于 2023-1-7 08:14
(2.1)里对0 3做stable=opt没意义
而且体系本身基态是单重态闭壳层的，(3)里没必要写U，否则计算白白增加很 ...

谢谢 Sob 老师的回复

您的意思是，在 (2.1)步的时候，用 (1.2)步里 T1 的优化后的结构，用 stable=opt 结合 0 1 计算单点能得到稳定的基态波函数。

然后，还是用  (1.2)步里 T1 的优化后的结构，同时读取刚刚计算得到的稳定波函数，用 0 1 结合 TD=triplet 计算三重态激发态

我这样理解对吗？

(2.1)里对0 3做stable=opt没意义
而且体系本身基态是单重态闭壳层的，(3)里没必要写U，否则计算白白增加很多耗时

应当做的是，对0 1做stable=opt，得到稳定的基态波函数，然后0 1结合TD=triplet算三重态激发态的时候读取刚才的chk

读取某一步的结构，并不意味着要读取同样这一步的波函数。

我看你是想计算T1的结构，那么你的(1.2)已经得到T1的结构和能量了，没必要再用TDDFT。

假如非要用TDDFT再算一遍，应该读取(1.2)的结构，并在此结构上获得单重态的稳定波函数，这个波函数和(1.1)(2.1)都不一样。

当然，(1.2)里面的泛函是可以换的，如要换泛函做TDDFT，则需用同样的泛函做(1.2)，并用同样泛函获取单重态的稳定波函数。

hebrewsnabla 发表于 2023-1-6 15:47
你的TDDFT是基于单重态波函数来做的，应该读取单重态的稳定波函数。

谢谢回复

为了不产生歧义，我在原贴上标记了每一步的编号。
您的意思是，我在 步骤(3) 应该读取 步骤 (1.1) 的 S0 几何优化后的 .chk 文件计算 tddft ，  而不是读取 步骤 (2.1) 的 .chk 文件？

你的TDDFT是基于单重态波函数来做的，应该读取单重态的稳定波函数。