计算化学公社

标题: 求助ORCA进行SCF收敛时能量变化很小 [打印本页]

作者
Author: Teresa 时间: 2025-3-20 21:10
标题: 求助ORCA进行SCF收敛时能量变化很小
各位老师，我在进行单点能计算时，SCF的收敛似乎震荡了，这是我的输入文件
! B1LYP D3 def2-TZVPP def2/J RIJCOSX tightSCF noautostart miniprint%maxcore    2000
%pal nprocs 8 end
* xyz 0 6
Ti    0.52110186 0.11889480 0.09456427
Ti    1.59841912 2.28375665 -0.83047139
Ti -0.84051084 2.31446306 0.51074888
Ti -0.59732847 1.34517173 -1.96581057
Ti -2.09809666 0.03480885 -0.01813699
Ti    1.75772708 -0.01907120 -2.19228759
Ti    1.46939454 1.66327692 1.90078577
Ti    3.12436797 0.08982412 0.07275700
Ti    1.57853288 -1.13341671 2.09721151
Ti    1.74522964 -2.12763111 -0.38498059
Ti -0.61704707 -2.09997290 0.98133588
Ti -0.85323750 0.20180149 2.31819566
Ti -0.52312971 -1.47025422 -1.76208672
Na    1.80806481 -3.01869348 -3.55336724
*
输出在进行TRAH时，我发现能量的变化已经非常小了，如下所示
5    -11204.727429931605    3.593258e-05                         (NR MAcro)
6    -11204.727429941937    3.372900e-05                         (NR MAcro)
7    -11204.727429943960    3.532584e-05                         (NR MAcro)
8    -11204.727429948274    2.921376e-05                         (NR MAcro)
9    -11204.727429950606    3.191334e-05                         (NR MAcro)
10    -11204.727429954753    2.856538e-05                         (NR MAcro)
11    -11204.727429957413    3.276574e-05                         (NR MAcro)
12    -11204.727429962280    3.124032e-05                         (NR MAcro)

是否可以取这个能量作为最后计算的能量呢?

作者
Author: sobereva 时间: 2025-3-21 07:26
用B1LYP很奇怪，没什么理由用这个

trah经常耽误事，用notrah，并且搜索论坛找过往关于ORCA SCF不收敛的讨论。

作者
Author: Teresa 时间: 2025-3-21 08:17

sobereva 发表于 2025-3-21 07:26
用B1LYP很奇怪，没什么理由用这个

trah经常耽误事，用notrah，并且搜索论坛找过往关于ORCA SCF不收敛的 ...

好的谢谢sob老师，用B1LYP是因为这个项目需要做一个benchmark，需要计算一下。老师再请教一下用notrah会影响最后的结果吗，因为我之前在计算GGA泛函时没有加

作者
Author: sobereva 时间: 2025-3-21 08:29

Teresa 发表于 2025-3-21 08:17
好的谢谢sob老师，用B1LYP是因为这个项目需要做一个benchmark，需要计算一下。老师再请教一下用notrah会 ...

不会（除非恰巧收敛到了不同波函数）

作者
Author: wzkchem5 时间: 2025-3-21 08:30

Teresa 发表于 2025-3-21 08:17
好的谢谢sob老师，用B1LYP是因为这个项目需要做一个benchmark，需要计算一下。老师再请教一下用notrah会 ...

不会，因为默认情况下程序仅当有必要时才打开TRAH，也就是有可能打开有可能不打开；NoTRAH意味着永远不打开。既然有可能打开有可能不打开的结果，不会自己和自己不可比，那么永远不打开的结果必然和有可能打开有可能不打开可比

作者
Author: Teresa 时间: 2025-3-21 09:32

sobereva 发表于 2025-3-21 08:29
不会（除非恰巧收敛到了不同波函数）

好的谢谢老师

作者
Author: Teresa 时间: 2025-3-21 09:33

wzkchem5 发表于 2025-3-21 08:30
不会，因为默认情况下程序仅当有必要时才打开TRAH，也就是有可能打开有可能不打开；NoTRAH意味着永远不打 ...

好的明白啦谢谢老师

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