计算化学公社

 找回密码 Forget password
 注册 Register

orca SCF迷一样的感受

查看数: 15105 | 评论数: 48 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
413
发布时间: 2022-3-10 00:35

正文摘要:

各位大家好我最近优化了一个结构,跑了一天,然后拿优化好的结构重新做个单点算点性质 奇怪的地方是未优化的结构SCF很快就收敛了 orca自己优化好的结构,SCF竟然报错 我也是也将之前的gbw文件拷贝成新文件名的gbw ...

回复 Reply

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呢,这不是增加收敛难度么

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-17 09:23 , Processed in 0.166998 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list