计算化学公社

标题: 如何使用Gaussian计算双自由基电离能 [打印本页]

作者
Author: Trevor77 时间: 2024-1-23 22:42
标题: 如何使用Gaussian计算双自由基电离能
本帖最后由 Trevor77 于 2024-1-23 23:04 编辑

使用B3lyp优化结构，再使用G4计算电荷为0和1时能量，在计算普通分子时很正常，但计算自由基尤其是双自由基时，常常报错或结果偏差很大，怀疑双自由基不适用于这种方法，读了很多帖子没有非常相关的，尤其是在电荷为1时不知该如何设置，请教各位老师正确的算法。计算的分子如图。

优化关键词为：#p opt freq ub3lyp/6-311++g(d,p)
能量关键词为：#p g4
neutron和ion电荷和自旋分别为：0 1 和 1 2（如果报错我会改自旋多重度）

作者
Author: zjxitcc 时间: 2024-1-23 22:46
你先展示一下你“使用B3lyp优化结构”输入文件关键词，让大家看看你结构优化有没做对

作者
Author: Trevor77 时间: 2024-1-23 23:03

zjxitcc 发表于 2024-1-23 22:46
你先展示一下你“使用B3lyp优化结构”输入文件关键词，让大家看看你结构优化有没做对

已补充

作者
Author: sobereva 时间: 2024-1-23 23:19
G4没法拿来算双自由基。算普通自由基没有问题
此体系也没必要用G4，直接用CCSD(T)算电子能量部分就完了
弥散函数加得完全多余，毫无意义

算双自由基也不是光写个U就了事的，仔细看
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

作者
Author: zjxitcc 时间: 2024-1-24 03:26
本帖最后由 zjxitcc 于 2024-1-24 03:28 编辑

Trevor77 发表于 2024-1-23 23:03
已补充

两种任务的关键词都是错的，没有构造对称破缺初猜，花了约2倍机时算出闭壳层单重态结果，只能全重算了。重算的时候注意，在结构优化中，初始结构除了需要构造对称破缺初猜，最好还应检验波函数稳定性。

作者
Author: Trevor77 时间: 2024-1-24 21:56

sobereva 发表于 2024-1-23 23:19
G4没法拿来算双自由基。算普通自由基没有问题
此体系也没必要用G4，直接用CCSD(T)算电子能量部分就完了
...

谢谢您的指导，我再重新计算一下

作者
Author: Trevor77 时间: 2024-1-24 21:56

zjxitcc 发表于 2024-1-24 03:26
两种任务的关键词都是错的，没有构造对称破缺初猜，花了约2倍机时算出闭壳层单重态结果，只能全重算了。 ...

谢谢您的指导，我再重新计算一下

作者
Author: Trevor77 时间: 2024-1-27 17:22

sobereva 发表于 2024-1-23 23:19
G4没法拿来算双自由基。算普通自由基没有问题
此体系也没必要用G4，直接用CCSD(T)算电子能量部分就完了
...

请教一下老师，我用新的算法算出来数值和之前差别不大，是哪里用错了吗
优化
# opt freq def2tzvp um062x guess=mix nosymm
优化波函数
# def2tzvp um062x guess=mix nosymm stable=opt
计算
# ccsd(t) def2tzvp guess=read

作者
Author: sobereva 时间: 2024-1-28 02:46

Trevor77 发表于 2024-1-27 17:22
请教一下老师，我用新的算法算出来数值和之前差别不大，是哪里用错了吗
优化
# opt freq def2tzvp um06 ...

“优化波函数”意义不明
几何优化过程中必须通过查看<S**2>确保得到了对称破缺波函数，并且最好在几何优化后，用guess=read读取优化任务的chk做波函数稳定性测试。只有优化过程的波函数是稳定的，优化的结果才有意义

说清楚“偏差很大”是相对于什么，当前结果是多少、参照数据是多少、参照数据哪里来的，搞清楚和当前计算有无可比性

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