计算化学公社

标题: ORCA6.0中使用wB97M(2)计算单重态开壳层氧气分子单点能出错求助 [打印本页]

作者
Author: MCF 时间: 2024-12-24 20:13
标题: ORCA6.0中使用wB97M(2)计算单重态开壳层氧气分子单点能出错求助
本帖最后由 MCF 于 2024-12-24 20:10 编辑

我用wB97M(2)计算三重态和单重态开壳层氧气分子的单点能的时候发现一个问题，就是算单重态开壳层氧气发现进行第二步wB97M(2)计算的时候ORCA6.0好像忽略了一半的电子数，体系只有8个电子被填充在轨道内，而且能量也差很多，第一步的wB97M-V计算出来是-150多Hartree，但是第二步wB97M(2)算出来是-80多Hartree。三重态计算一切正常，还请各位大佬指点一下是哪里出错了。
下面分别是输入文件和部分计算结果，在算单点能之前两种电子态的氧气分子已经经过B3LYP-D4/def2-TZVP几何优化过了。
单重态开壳层氧气分子单点能输入文件
(, 下载次数 Times of downloads: 19)
单重态wB97M-V计算能量结果
(, 下载次数 Times of downloads: 21)
单重态wB97M(2)计算能量结果
(, 下载次数 Times of downloads: 17)
单重态电子占据轨道情况
(, 下载次数 Times of downloads: 22)
三重态输入文件
(, 下载次数 Times of downloads: 21)
三重态wB97M-V计算能量结果
(, 下载次数 Times of downloads: 19)
三重态wB97M(2)计算能量结果
(, 下载次数 Times of downloads: 16)

作者
Author: zjxitcc 时间: 2024-12-24 20:34
本帖最后由 zjxitcc 于 2024-12-24 20:37 编辑

你自己说“单重态开壳层氧气分子”，但帖子内容似乎没看到你构建了对称破缺初猜，也没看到检验单重态UKS波函数稳定性，也没看到读取对称破缺UKS轨道，如何证明你算的是开壳层单重态？别用双杂化算了半天，发现概念搞错了。
后记：我看到guessmix了，你是想通过它来构建对称破缺初猜？

作者
Author: MCF 时间: 2024-12-24 20:51

zjxitcc 发表于 2024-12-24 20:34
你自己说“单重态开壳层氧气分子”，但帖子内容似乎没看到你构建了对称破缺初猜，也没看到检验单重态UKS波 ...

图中单重态输入文件里面关键词加了guessmix，而且wB97M-V计算S平方算符的期望值接近于1，应该能说明算的是对称性破缺的电子态吧。
(, 下载次数 Times of downloads: 20)

作者
Author: zjxitcc 时间: 2024-12-24 21:03

MCF 发表于 2024-12-24 20:51
图中单重态输入文件里面关键词加了guessmix，而且wB97M-V计算S平方算符的期望值接近于1，应该能说明算的 ...

嗯，这看起来是合理的。建议在双杂化的任务里也加上UKS关键词，看是否还有你帖子里描述的问题。
另外，为何要写MaxIter 1和IgnoreConv 1？

作者
Author: MCF 时间: 2024-12-24 21:14

zjxitcc 发表于 2024-12-24 21:03
嗯，这看起来是合理的。建议在双杂化的任务里也加上UKS关键词，看是否还有你帖子里描述的问题。
另外 ...

MaxIter 1和IgnoreConv 1这俩关键词是手册原文里提到的，应该有特殊作用。
(, 下载次数 Times of downloads: 20)

作者
Author: MCF 时间: 2024-12-24 21:55

zjxitcc 发表于 2024-12-24 21:03
嗯，这看起来是合理的。建议在双杂化的任务里也加上UKS关键词，看是否还有你帖子里描述的问题。
另外 ...

加了uks以后结果正常了，感谢。
(, 下载次数 Times of downloads: 21)

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