计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: 冰释之川
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] ORCA中SCF迭代出现迷之收敛的案例

[复制链接 Copy URL]

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

16#
发表于 Post on 2023-10-17 16:51:28 | 只看该作者 Only view this author
413 发表于 2023-10-17 09:36
手册给了sthresh的区间是1e-5到1e-8。
请问如果我把sthresh设置到1e-3合理吗?因为我试了-4、-9,也有问 ...

其实1e-5已经有些太大了,再大的话结果肯定是不能用的
sthresh只能解决轨道基的线性相关问题,如果调sthresh没用,考虑可能是辅助基线性相关问题,需要用pivoted Cholesky decomposition (PCD)解决,手册里有讲
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

419

帖子

0

威望

1379

eV
积分
1798

Level 5 (御坂)

17#
发表于 Post on 2023-10-17 17:01:14 | 只看该作者 Only view this author
wzkchem5 发表于 2023-10-17 16:51
其实1e-5已经有些太大了,再大的话结果肯定是不能用的
sthresh只能解决轨道基的线性相关问题,如果调sth ...

手册看到了这个命令,想请教下,这个的default是关闭的是吗?
因为手册看到了一句话“This can be requested for any of the basis sets using either the overlap or the Coulomb metric.” 也就是需要指定才打开这个选项?

手册给了示例,如下。只需要完整的copy到输入文件就可以是吗?需要修改啥吗?
%basis
PCDTrimBas Overlap # Trim the orbital basis in the overlap metric
PCDTrimAuxJ Coulomb # Trim the AuxJ basis in the Coulomb metric
PCDTrimAuxJK Coulomb # Trim the AuxJK basis in the Coulomb metric
PCDTrimAuxC Coulomb # Trim the AuxC basis in the Coulomb metric
PCDThresh -1 # Threshold for the PCD: chosen automatically if <0
end

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

18#
发表于 Post on 2023-10-17 17:46:07 | 只看该作者 Only view this author
413 发表于 2023-10-17 10:01
手册看到了这个命令,想请教下,这个的default是关闭的是吗?
因为手册看到了一句话“This can be reque ...

对,默认关闭
全部拷进输入文件是可以的,但为防止误导其他看这个输入文件的人,稍好的做法是用到哪个辅助基组就拷哪一行。例如如果你做的是DFT计算且没有用RIJK,那AuxJK和AuxC就可以不拷,拷了不算错但是容易误导别人以为这两行是必需的
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

419

帖子

0

威望

1379

eV
积分
1798

Level 5 (御坂)

19#
发表于 Post on 2023-10-17 17:49:54 | 只看该作者 Only view this author
wzkchem5 发表于 2023-10-17 17:46
对,默认关闭
全部拷进输入文件是可以的,但为防止误导其他看这个输入文件的人,稍好的做法是用到哪个辅 ...

好的,谢谢指导(手动敬礼)。我试下,看下结构变化。

419

帖子

0

威望

1379

eV
积分
1798

Level 5 (御坂)

20#
发表于 Post on 2023-10-18 00:37:50 | 只看该作者 Only view this author
wzkchem5 发表于 2023-10-17 16:51
其实1e-5已经有些太大了,再大的话结果肯定是不能用的
sthresh只能解决轨道基的线性相关问题,如果调sth ...

您好,我加了PCD的命令,好像还是会出现结构收缩的现象。
如果把sthresh设置的超小,比如-12或-14会有好的影响吗?

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

21#
发表于 Post on 2023-10-18 15:37:27 | 只看该作者 Only view this author
413 发表于 2023-10-17 17:37
您好,我加了PCD的命令,好像还是会出现结构收缩的现象。
如果把sthresh设置的超小,比如-12或-14会有好 ...

sthresh设得太小会导致本来存在的基组相关问题程序没有发现,所以反倒会加重你的问题。sthresh设得太大会导致结果不准,但不会导致SCF能量严重不合理、结构大幅收缩等严重问题。
要不你发一下你的输入输出文件吧
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 00:29 , Processed in 0.164351 second(s), 21 queries , Gzip On.

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