关于TDDFT优化三重态激发态T1/T2，激发能为负值如何消除

qingxin

zjxitcc 发表于 2022-8-28 18:10
否。我说的循环操作是指目标电子态的合理计算->出现负激发能->stable=opt算基态单点->目标电子态的结构优 ...

谢谢老师，这个操作看懂了，已经在算，只是不太明白多次stable=opt算基态单点，怎么会不同

qingxin

hebrewsnabla 发表于 2022-8-28 03:09
td单点就是只写td不写opt

老师，只算td基态单点的时候，得到波函数，opt时读取该波函数，得到的结构应该既不是基态，也不是T1结构。
按2L说的再次回到td基态单点，核验波函数，两次td之间结构一直是基态才对，这样理解对吗？

hebrewsnabla

qingxin 发表于 2022-8-28 21:22
老师，只算td基态单点的时候，得到波函数，opt时读取该波函数，得到的结构应该既不是基态，也不是T1结构 ...

看不懂你在说什么

qingxin

hebrewsnabla 发表于 2022-8-28 22:58
看不懂你在说什么

老师，可能我没表达清楚，就是问，您说结构变了，指的是什么结构，和我们再次td单点能有什么关系呢？

LittlePupil

我想你可能需要搞清楚计算过程中“基态”和“激发态”的概念。

很明显，结构优化过程中结构是在（至少程序是在尽力让结构）不断接近激发态的稳定构型的，这个时候做的波函数稳定性计算只是对于特定的结构算这个结构的基态稳定波函数，而实际上这个结构只是优化过程中的一个中间结构，既不是基态的稳定构型，也不是激发态的稳定构型，或者说，既不是常规实验上的基态构型，也不是常规实验上的激发态构型。

18L说的“结构变了”应该是指你两次做基态波函数稳定性计算的结构不一样，因为你在完成一次基态波函数稳定性计算后是要读取波函数后做激发态结构优化的。

qingxin

LittlePupil 发表于 2022-8-28 23:58
我想你可能需要搞清楚计算过程中“基态”和“激发态”的概念。

很明显，结构优化过程中结构是在（至少程 ...

谢谢，

是否可以理解这样：每次要换不同的结构（即前一次优化后的结构）来优化激发态，同时再次读取核验的基态波函数（这个不停核验，结果会变化吗？）

hebrewsnabla

14L说的循环操作的适用条件是：操作正确，但是td opt之后就是有负的激发能；

但是你的情况是，td opt的第一步（等同于td 单点）的激发能就是负的，说明操作不正确，还没有到考虑循环操作的时候。

所以我让你先把td 单点搞对再说（即没有负的激发能）。


“换不同的结构（即前一次优化后的结构）来优化激发态，同时再次读取核验的基态波函数（这个不停核验，结果会变化吗？）” -- 会。

qingxin

hebrewsnabla 发表于 2022-8-29 13:48
14L说的循环操作的适用条件是：操作正确，但是td opt之后就是有负的激发能；

但是你的情况是，td opt的 ...

老师辛苦了，您看我这样td 单点能对吗？

%chk=SQ1-THF-02.chk
%mem=12000MB
#p td=(triplets, nstates=6) ub3lyp/6-31g(d,p) scrf=(solvent=thf) stable=opt g09defaults

hebrewsnabla

qingxin 发表于 2022-8-29 17:41
老师辛苦了，您看我这样td 单点能对吗？

%chk=SQ1-THF-02.chk

我不记得高斯支不支持一个任务里面做stable=opt和td了。如果不支持你要分两步做。

LittlePupil

老师辛苦了，您看我这样td 单点能对吗？

http://bbs.keinsci.com/thread-26563-1-1.html

建议看看这篇帖子里z老师的回复，希望对你有帮助。

LittlePupil

qingxin 发表于 2022-8-29 12:22
谢谢，

是否可以理解这样：每次要换不同的结构（即前一次优化后的结构）来优化激发态，同时再次读取核 ...

“每次要换不同的结构（即前一次优化后的结构）来优化激发态”是，因为原理上正确操作的话结构应该越来越接近激发态稳定结构。

“同时再次读取核验的基态波函数”也对，但是需要读取的是本轮优化起点，也就是前次结构优化过程中激发能出现负值的那个结构（称为结构A）的稳定波函数，也就是说在此之前你需要先对结构A进行波函数稳定性计算。

“不停核验，结果会变化吗”，如果你是指“核验”结果，那确实可能不会变——理想情况下应该都是提示波函数不稳定后优化到稳定波函数，但我们真正关心的是结构的稳定波函数，这个当然会变——理论上，每次“结构A”的结构都不会和前次相同，对应的稳定波函数也会不同。

qingxin

LittlePupil 发表于 2022-8-29 19:35
“每次要换不同的结构（即前一次优化后的结构）来优化激发态”是，因为原理上正确操作的话结构应该越来越 ...

谢谢，那个贴也拜读了

我正在计算，到时再来回复结果

15689900879@163

zjxitcc 发表于 2022-8-27 19:39
读取UDFT波函数，加stable=opt进行单点计算、检验波函数稳定性、获得稳定波函数。然后读取稳定波函数，继续 ...

老师您好，在进行UDFT stable=opt 计算时，我在自旋多重度的定义上有疑问。
1. 对于S1态结构，出现负的激发能，自旋多重度设为1吗？
2. 对于T1态结构，出现负的激发能，自旋多重度设为3吗？

另外，我在查阅相关帖子时，发现有的在UDFT下进行stable=opt计算，有的在RDFT下进行stable=opt计算，使用U/RDFT有何要求吗？
在RDFT stable=opt计算时，S1态和T1态的自旋多重度又分别设置多少合适呢？

wzkchem5

15689900879@163 发表于 2024-3-27 10:05
老师您好，在进行UDFT stable=opt 计算时，我在自旋多重度的定义上有疑问。
1. 对于S1态结构，出现负的 ...

如果你做RDFT的stable=opt计算，且不存在任何instability，那么结果和不加stable=opt一样，这种情况下你应该已经知道该怎么做了；
如果你做RDFT的stable=opt计算，且只存在singlet instability（这种情况比较罕见），那么仍然是RDFT计算，自旋多重度写法都和不加stable=opt的RDFT计算一致；
如果你做RDFT的stable=opt计算，且存在triplet instability，那么程序会自动切换到UDFT，所以一切适用于UDFT的结论在这里都自动适用。
自旋多重度的写法参见https://mattermodeling.stackexch ... -triplets/9357#9357

15689900879@163

wzkchem5 发表于 2024-3-27 17:43
如果你做RDFT的stable=opt计算，且不存在任何instability，那么结果和不加stable=opt一样，这种情况下你 ...

谢谢老师的回复！您提供的链接内容我看了，不确定理解的对不对，麻烦您帮我看看。
1. 对于基态是S0，且为single-referential，使用RDFT TD(triplet) opt 0 1计算T1态，此时寻找T1的稳定基态波函数，RDFT stable=opt 0 1；
2. 对于基态是S0，且为multireferential，使用UDFT opt 0 3计算T1态，此时寻找T1的稳定基态波函数，UDFT stable=opt 0 3。
上述1和2得到的稳定基态波函数，均使用RDFT TD(triplet) opt guess=read 0 1 继续优化。

您提到，当进行RDFT的stable=opt计算，且存在triplet instability，那么程序会自动切换到UDFT，我就不确定1对不对了。