计算化学公社

标题: 单重态IRC验证后，三重态需要重新跑IRC吗？ [打印本页]

作者
Author: archer 时间: 2021-3-4 09:49
标题: 单重态IRC验证后，三重态需要重新跑IRC吗？
研究单重和三重态的反应路径，在单重态下找到过渡态并进行IRC验证。三重态的过渡态结构和单重态的结构非常接近，还需要再跑IRC吗？

另外，计算单重态时用unstricted DFT，同时也计算了三重态，还需要计算波函数的稳定性吗？直接根据能量判断基态是否合理。

作者
Author: zjxitcc 时间: 2021-3-4 09:55
毫无疑问一定要检验波函数稳定性，这看能量看不出什么（要说看自旋布局分析还有点可能性）。如果波函数不稳定，要先优化波函数至稳定（stable=opt），然后读取波函数继续做结构优化、找过渡态，因为波函数变了，力也变了。

三重态最好也跑一下IRC，尤其是如果三重态TS与单重态TS能量相差不多的时候，就更要跑了。

作者
Author: archer 时间: 2021-3-4 11:19

zjxitcc 发表于 2021-3-4 09:55
毫无疑问一定要检验波函数稳定性，这看能量看不出什么（要说看自旋布局分析还有点可能性）。如果波函数不稳 ...

谢谢！三重态和单重态的反应路径我都要研究。波函数稳定性，我只需要检验单重态是UHF还是RHF更稳定就可以吗？

作者
Author: zjxitcc 时间: 2021-3-4 11:27

archer 发表于 2021-3-4 11:19
谢谢！三重态和单重态的反应路径我都要研究。波函数稳定性，我只需要检验单重态是UHF还是RHF更稳定就可以 ...

单重态、三重态都要检验波函数稳定性。我不知道你这里的“UHF”是泛指UHF/UDFT 还是真的仅指UHF，你用UDFT算的，就要在UDFT下检验波函数稳定性，不能跑到UHF下去检验（这是错的），计算级别要一致。

作者
Author: archer 时间: 2021-3-4 11:41

zjxitcc 发表于 2021-3-4 11:27
单重态、三重态都要检验波函数稳定性。我不知道你这里的“UHF”是泛指UHF/UDFT 还是真的仅指UHF，你用U ...

好的，谢谢。我用的UDFT。

作者
Author: archer 时间: 2021-3-12 15:38

zjxitcc 发表于 2021-3-4 11:27
单重态、三重态都要检验波函数稳定性。我不知道你这里的“UHF”是泛指UHF/UDFT 还是真的仅指UHF，你用U ...

你好，我用UDFT优化单重态，波函数检查不稳定。优化波函数，复制.chk，geom=check, guess=read，重新优化无法收敛。请问这种情况如何解决，谢谢！

作者
Author: zjxitcc 时间: 2021-3-12 15:45

archer 发表于 2021-3-12 15:38
你好，我用UDFT优化单重态，波函数检查不稳定。优化波函数，复制.chk，geom=check, guess=read，重新优化 ...

建议使用“优化波函数”、“优化结构”、“SCF不收敛”，“结构优化不收敛”等术语，这样区分度明显。在这个帖子里，说“用UDFT优化单重态”不知道你是在优化波函数，还是优化结构。

照我的理解（不知是否有误）：你先用UDFT在单重态下做结构优化，完了发现波函数不稳定，然后优化波函数至稳定，接着读取该波函数继续做结构优化。然后碰到了“不收敛”问题，但是你没讲是“SCF不收敛”，“结构优化不收敛”，需对症下药：
《解决SCF不收敛问题的方法》http://sobereva.com/61
《量子化学计算中帮助几何优化收敛的常用方法》http://sobereva.com/164

作者
Author: archer 时间: 2021-3-12 17:03

zjxitcc 发表于 2021-3-12 15:45
建议使用“优化波函数”、“优化结构”、“SCF不收敛”，“结构优化不收敛”等术语，这样区分度明显。在 ...

谢谢回复，我的描述就是这个意思，我用优化后的波函数作为初始猜测，重新用UDFT优化单重态结构。我有一些波函数不稳定的结构，重新优化基本都会遇到SCF不收敛。

作者
Author: 嘤嘤嘤 时间: 2021-9-11 16:47
请问你三重态下的过渡态怎么找的啊？和单重态的思路一样吗？只不过用三重态去优化结构？

作者
Author: archer 时间: 2021-9-11 20:39

嘤嘤嘤发表于 2021-9-11 16:47
请问你三重态下的过渡态怎么找的啊？和单重态的思路一样吗？只不过用三重态去优化结构？

是的，自旋为3进行优化

作者
Author: 嘤嘤嘤 时间: 2021-9-12 16:47
那要是不收敛怎么办啊？我优化好多都是link9999

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