orca SCF迷一样的感受

413 · 发表于 Post on 2022-3-10 00:35:27

各位大家好我最近优化了一个结构，跑了一天，然后拿优化好的结构重新做个单点算点性质
奇怪的地方是未优化的结构SCF很快就收敛了
orca自己优化好的结构，SCF竟然报错
我也是也将之前的gbw文件拷贝成新文件名的gbw，但是仍未解决问题
请问大家有这种体验吗？怎么解决这种问题吗？或者我哪里做的不对才会出现这种情况？

wzkchem5 · 发表于 Post on 2023-10-9 16:08:15

mfdsrax2 发表于 2023-10-9 07:56
CCSD(T)的解析梯度有希望吗？

可能可以来得及做完，现在不好说。

mfdsrax2 · 发表于 Post on 2023-10-9 14:56:12

wzkchem5 发表于 2023-10-7 16:07
有很多，我暂时能想到的：sf-X2C；CCSD解析梯度；解析拉曼强度；构象搜索；分子对接（docking）；内转换 ...

CCSD(T)的解析梯度有希望吗？

wzkchem5 · 发表于 Post on 2023-10-7 16:07:47

413 发表于 2023-10-7 02:08
想问下6.0功能更新有啥big news吗？大概啥时候会正式对外？

有很多，我暂时能想到的：sf-X2C；CCSD解析梯度；解析拉曼强度；构象搜索；分子对接（docking）；内转换速率；新的半经验方法（NOTCH）；等等。
估计明年初发布

413 · 发表于 Post on 2023-10-7 09:08:51

本帖最后由 413 于 2023-10-7 09:31 编辑

wzkchem5 发表于 2023-10-6 22:56
DKH2我之前遇到过对于一些比较重的元素（比如Bi）有variational instability，会严重低估Bi-C键长（至少0 ...

想问下6.0功能更新有啥big news吗？大概啥时候会正式对外？

wzkchem5 · 发表于 Post on 2023-10-6 22:56:34

zjxitcc 发表于 2023-10-6 14:00
你的困惑是对的。但现在是2023年，考虑标量相对论应该使用DKH2或sfX2C哈密顿，而不是用完全过时的ZORA（ ...

DKH2我之前遇到过对于一些比较重的元素（比如Bi）有variational instability，会严重低估Bi-C键长（至少0.5埃以上）。所以我现在做这种计算，如果在orca里做的话，一般都是ZORA算opt+freq，sf-X2C算单点（我们内部orca版本已经实现了sf-X2C的单点和解析梯度，但暂时没有解析Hessian），必要情况下单点计算还可以微扰考虑SOC。如果用BDF做就是统一用sf-X2C了，同样的，单点如果有必要的话加SOC。

zjxitcc · 发表于 Post on 2023-10-6 21:00:32

413 发表于 2023-10-6 18:03
我也试过mokit（现在基本是首选项）
但有个小困惑，就是高斯算不了zora

你的困惑是对的。但现在是2023年，考虑标量相对论应该使用DKH2或sfX2C哈密顿，而不是用完全过时的ZORA（不是ORCA的问题，是ZORA的问题）http://bbs.keinsci.com/thread-550-1-1.html

413 · 发表于 Post on 2023-10-6 18:03:47

本帖最后由 413 于 2023-10-6 18:06 编辑

zjxitcc 发表于 2023-10-6 17:23
存在正相关关系。若体系存在基函数线性相关，SCF收敛会变难，即使正确去掉了线性相关，收敛也比常规体系 ...

我也试过mokit（现在基本是首选项）
但有个小困惑，就是高斯算不了zora

413 · 发表于 Post on 2023-10-6 17:48:27

zjxitcc 发表于 2023-10-6 17:21
你有没检验收敛后的波函数稳定性

波函数是稳定的。

   *** CONVERGENCE OF RESIDUAL NORM REACHED ***

-----------------------------
SCF STABILITY ANALYSIS RESULT
-----------------------------

RHF/RKS->UHF/UKS - triplet - external

Root Eigenvalue (au)
   0    0.046522
   1    0.047015
   2    0.047309

Stability Analysis indicates a stable HF/KS wave function.

zjxitcc · 发表于 Post on 2023-10-6 17:23:14

413 发表于 2023-10-6 17:12
请教下，是另外一个zora的小分子吸附在过渡金属的体系，"opt pbe ZORA/RI ZORA-def2-TZVP SARC/J KDIIS N ...

存在正相关关系。若体系存在基函数线性相关，SCF收敛会变难，即使正确去掉了线性相关，收敛也比常规体系要慢。
另外，建议使用fch2mkl小程序传轨道给ORCA，该小程序在传轨道时就考虑了基函数线性相关（如果存在的话），在inp文件里写入了相应关键词，而非等用户在ORCA里摸爬滚打一阵子之后检查才发现。

413 · 发表于 Post on 2023-10-6 17:22:32

zjxitcc 发表于 2023-10-6 17:21
你有没检验收敛后的波函数稳定性

我做下

zjxitcc · 发表于 Post on 2023-10-6 17:21:13

413 发表于 2022-3-14 20:37
我也不太明白。
我试过loose, normal，能量变化走的没样
tight, verytight能分别走到10e-4左右和10e-6 ...

你有没检验收敛后的波函数稳定性

413 · 发表于 Post on 2023-10-6 17:12:25

wzkchem5 发表于 2022-3-14 15:55
试试把thresh、tcut设严一些，此外检查一下有没有基组线性相关

请教下，是另外一个zora的小分子吸附在过渡金属的体系，"opt pbe ZORA/RI ZORA-def2-TZVP SARC/J KDIIS NOSOSCF defgrid3 NOFINALGRIDX notrah"
我检查了“Smallest eigenvalue”，确实如手册那样存在线性相关问题。于是调整了thresh的阈值，现在还在算，还不知道能不能收敛。但我检查了另一个体系，同样存在线性问题，最后也是收敛了。
我想问下线性相关问题与收敛之间存在必然联系吗？

413 · 发表于 Post on 2022-3-14 20:37:19

本帖最后由 413 于 2022-3-16 07:33 编辑

hebrewsnabla 发表于 2022-3-14 18:39
为什么要用verytight呢，这不是增加收敛难度么

我也不太明白。
我试过loose, normal，能量变化走的没样
tight, verytight能分别走到10e-4左右和10e-6左右

坐等大神解释下orca里scf收敛标准是否（和高斯里scf收敛标准一样？）只影响threshold，不影响其他的

orca吊打高斯好几圈
2022.03.16更新体验
Trah超调

hebrewsnabla · 发表于 Post on 2022-3-14 18:39:47

413 发表于 2022-3-14 09:36
请教下，这是不是快接近收敛了？还能加个啥命令呢？
angs PRINTBASIS PRINTMOS pbe ZORA/RI ZORA-def2-T ...

为什么要用verytight呢，这不是增加收敛难度么

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] orca SCF迷一样的感受

浏览过的版块