计算化学公社

 找回密码 Forget password
 注册 Register
Views: 4372|回复 Reply: 37
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] ORCA算单点能SCF不收敛

[复制链接 Copy URL]

210

帖子

0

威望

564

eV
积分
774

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
我的输入行是! PW6B95 D3 def2-QZVP def2/J def2-QZVPP/C RIJCOSX tightSCF noautostart miniprint nopop NoTRAH
跑了125圈,2天时间,报错如下:
     ERROR
     This wavefunction IS NOT CONVERGED!
     And it is 'too far' from meeting the convergence thresholds.
     Aborting ...
     Please restart calculation (with larger maxiter?)


请问该怎么重新设置?

210

帖子

0

威望

564

eV
积分
774

Level 4 (黑子)

2#
 楼主 Author| 发表于 Post on 2024-2-21 09:12:46 | 只看该作者 Only view this author
计算的分子是148个原子。我用def2-tzvp+gcp跑过能收敛,但是用def2-QZVP就不收敛。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
Keys + 1 同样用TZVP收敛而QZVP不收敛

查看全部评分 View all ratings

71

帖子

0

威望

1187

eV
积分
1258

Level 4 (黑子)

3#
发表于 Post on 2024-2-21 09:44:04 | 只看该作者 Only view this author
你读一下收敛过的波函数当初猜呗

210

帖子

0

威望

564

eV
积分
774

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2024-2-21 13:46:51 | 只看该作者 Only view this author
Strange 发表于 2024-2-21 09:44
你读一下收敛过的波函数当初猜呗

请问读取gbw文件的关键词是什么呀?

74

帖子

0

威望

1857

eV
积分
1931

Level 5 (御坂)

5#
发表于 Post on 2024-2-21 14:02:45 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-21 13:46
请问读取gbw文件的关键词是什么呀?

http://bbs.keinsci.com/thread-31423-1-1.html

4289

帖子

4

威望

9541

eV
积分
13910

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2024-2-21 15:23:49 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-21 13:46
请问读取gbw文件的关键词是什么呀?

不需要添加任何关键词,只要inp和gbw文件同名,程序自动会读入轨道和基组,投影到当前目标基组上。换句话说,你只要复制一份gbw文件就行。当然,如果你非要绕远路,添加一些关键词,也是可以实现的。
自动做多参考态计算的程序MOKIT

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

7#
发表于 Post on 2024-2-21 15:36:11 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-21 09:12
计算的分子是148个原子。我用def2-tzvp+gcp跑过能收敛,但是用def2-QZVP就不收敛。

有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员必读贴里明确说了。

写def2-QZVPP/C多余,又不是双杂化泛函

也可以尝试defgrid3,用KDIIS等其它做法,本论坛里有很多讨论
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

210

帖子

0

威望

564

eV
积分
774

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2024-2-21 17:34:41 | 只看该作者 Only view this author
sobereva 发表于 2024-2-21 15:36
有别人回复之前若需要对帖子进行修改、补充,应直接编辑原帖,不要通过回帖进行补充,这点在置顶的新社员 ...

请问orca3.0.1的grid=5和grid=6,与orca5.0.4的defgrid3是一个档次吗?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

9#
发表于 Post on 2024-2-21 18:38:04 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-21 17:34
请问orca3.0.1的grid=5和grid=6,与orca5.0.4的defgrid3是一个档次吗?

没有严格的可比性
只能说grid6和defgrid3的格点质量都已经非常高了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

210

帖子

0

威望

564

eV
积分
774

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2024-2-21 19:31:57 | 只看该作者 Only view this author
sobereva 发表于 2024-2-21 18:38
没有严格的可比性
只能说grid6和defgrid3的格点质量都已经非常高了

一个很基础的问题,积分格点的设置的精度高是解决了什么问题,为什么要设置defgrid3?

210

帖子

0

威望

564

eV
积分
774

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2024-2-21 19:57:45 | 只看该作者 Only view this author
zjxitcc 发表于 2024-2-21 15:23
不需要添加任何关键词,只要inp和gbw文件同名,程序自动会读入轨道和基组,投影到当前目标基组上。换句话 ...

有简单的方法肯定直接用,我复制了gbw,改了名字和inp一样。感觉读取收敛的波函数是最靠谱的方法。

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

12#
发表于 Post on 2024-2-21 20:23:59 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-21 12:31
一个很基础的问题,积分格点的设置的精度高是解决了什么问题,为什么要设置defgrid3?

格点精度高,数值积分的精度就高,数值误差就小。有的时候SCF不收敛是因为数值误差太大,计算出来的能量和真实能量之间有误差,虽然误差的幅度小,但是符号随机,导致波函数一直小幅度震荡。这时增大格点就是最有效的方法。
参见https://bdf-manual.readthedocs.i ... tml#scfconvproblems
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

13#
发表于 Post on 2024-2-21 20:34:32 | 只看该作者 Only view this author
Oliviaw 发表于 2024-2-21 19:31
一个很基础的问题,积分格点的设置的精度高是解决了什么问题,为什么要设置defgrid3?

积分格点质量不够高可能导致KS矩阵构造精度差,对于某些情况可能导致SCF难收敛。例如之前我用wB97M-V结合很大基组计算18碳环的弱相互作用复合物就遇到这种情况,默认的defgrid2 SCF震荡,改成defgrid3收敛了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

210

帖子

0

威望

564

eV
积分
774

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2024-2-21 21:17:35 | 只看该作者 Only view this author
wzkchem5 发表于 2024-2-21 20:23
格点精度高,数值积分的精度就高,数值误差就小。有的时候SCF不收敛是因为数值误差太大,计算出来的能量 ...

defgrid3是只在能力计算的时候设置,还是说优化和频率也要统一用defgrid3呢?
如果是已经优化收敛的结构,要补充加defgrid3,我只需读取优化收敛的结构的gwb,然后inp里重新设置defgrid3和opt就可以吗?
还是说我就用优化收敛的结构的坐标,然后重新加defgrid3,opt呢?

210

帖子

0

威望

564

eV
积分
774

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2024-2-21 21:25:53 | 只看该作者 Only view this author
sobereva 发表于 2024-2-21 20:34
积分格点质量不够高可能导致KS矩阵构造精度差,对于某些情况可能导致SCF难收敛。例如之前我用wB97M-V结合 ...

不知道18碳环的文章会不会探究lewis acidity/basicity和charge transfer?

我的体系还要讲不同金属配位对弱相互作用的影响,相关的文章都会提到酸性和电荷转移,但是没见到可以计算酸性的强度的文章。charge transfer有办法可以计算出一个电荷变化的数值,也不知道charge transfer 会不会影响interaction?

那个能量分解的方法有个静电相互作用的项,跟charge transfer有没有关系啊?

您的文章有没有酸性的讨论和charge transfer呢?

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 23:19 , Processed in 0.287804 second(s), 21 queries , Gzip On.

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