计算化学公社

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

[ORCA] 请教大纯有机体系优化建议(ORCA的SCF不收敛)

[复制链接 Copy URL]

541

帖子

3

威望

6622

eV
积分
7223

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 ABetaCarw 于 2020-8-19 19:36 编辑

我是用ORCA优化一个453原子的纯有机体系,我已经用MOPAC初优化过了结构,然后用以下关键进行ORCA计算
! BLYP D3 GCP(DFT/SVP) def2-SVP def2/J opt
%maxcore  2000
%pal nprocs   36 end
%base "comp2"
* xyz 0 1


然后发现每次都报错




  1.            ************************************************************
  2.            *        Program running with 36 parallel MPI-processes    *
  3.            *              working on a common directory               *
  4.            ************************************************************
  5. ------------------------------------------------------------------------------
  6.                          ORCA SCF GRADIENT CALCULATION
  7. ------------------------------------------------------------------------------

  8. Gradient of the Kohn-Sham DFT energy:
  9. Kohn-Sham wavefunction type      ... RKS
  10. Hartree-Fock exchange scaling    ...    0.000
  11. Number of operators              ...    1
  12. Number of atoms                  ...  453
  13. Basis set dimensions             ... 4464
  14. Integral neglect threshold       ... 2.5e-11
  15. Integral primitive cutoff        ... 2.5e-12

  16. Nuclear repulsion gradient       ... done
  17. One Electron Gradient            ... done
  18. Pre-screening matrix             ... done
  19. RI-J gradient                    ...
  20. ORCA finished by error termination in SCF gradient
  21. Calling Command: mpirun -np 36  /home/XXX/ORCA/orca421/orca_scfgrad_mpi comp2.scfgrad.inp comp2
  22. [file orca_tools/qcmsg.cpp, line 458]:
  23.   .... aborting the run

复制代码

请熟悉ORCA的前辈 指导一下问题出在哪里。
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

2#
发表于 Post on 2020-8-19 19:39:09 | 只看该作者 Only view this author
求梯度内存不够吧。。。。可以降低并行核数,再就是orca的几何优化其实很蠢的,我算了一个很小的配合物,优化完了振动分析3个虚频。。。。

6万

帖子

99

威望

6万

eV
积分
125189

管理员

公社社长

3#
发表于 Post on 2020-8-19 20:14:10 | 只看该作者 Only view this author
这和SCF不收敛无关,纯粹是受力计算部分的问题

比起BLYP-D3-gCP,我明显更推荐B97-3c,精度更好,耗时差不多
B97-3c也这样的话,加大内存,如果物理内存不够,减少核数再试
北京科音自然科学研究中心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

541

帖子

3

威望

6622

eV
积分
7223

Level 6 (一方通行)

4#
 楼主 Author| 发表于 Post on 2020-8-19 23:59:57 | 只看该作者 Only view this author
高斯的报错信息 我靠本论坛和谷歌能够99%解决。ORCA的报错信息我是真的无助啊。以前我以为是并行出错了,这次我判定SCF报错的,我以为我找对了……感谢社长和喵星大佬的支持
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

5#
发表于 Post on 2020-8-20 00:22:56 | 只看该作者 Only view this author
ABetaCarw 发表于 2020-8-19 23:59
高斯的报错信息 我靠本论坛和谷歌能够99%解决。ORCA的报错信息我是真的无助啊。以前我以为是并行出错了,这 ...

不是写着finished by error terminated in SCF gradient呢嘛

82

帖子

4

威望

2540

eV
积分
2702

Level 5 (御坂)

6#
发表于 Post on 2020-8-20 01:04:57 | 只看该作者 Only view this author
体系这么大!%maxcore不给够肯定不行,估计要放大5~10倍

113

帖子

0

威望

3156

eV
积分
3269

Level 5 (御坂)

7#
发表于 Post on 2020-8-20 01:31:20 | 只看该作者 Only view this author
核数不超过10,内存单核16-20 GB试试,我也刚遇到过类似问题。就是这么解决的

541

帖子

3

威望

6622

eV
积分
7223

Level 6 (一方通行)

8#
 楼主 Author| 发表于 Post on 2020-8-20 08:51:59 | 只看该作者 Only view this author
本帖最后由 ABetaCarw 于 2020-8-20 09:51 编辑

感谢楼上各位,问题确实去你们所说 物理内存不够。按照七楼建议 跑了一夜没有报错。
果然 计算的问题还是要向靠谱的人求助。我在researchgate上面找了两个讨论该报错的帖子 全是在讨论 改收敛精度 改SCF方法,太坑了。
恍惚月余,深谙人与人之间的差距。以后还应努力学习,才能与强者比肩。

1万

帖子

0

威望

9892

eV
积分
22140

Level 6 (一方通行)

9#
发表于 Post on 2020-8-20 09:50:30 | 只看该作者 Only view this author
ABetaCarw 发表于 2020-8-19 23:59
高斯的报错信息 我靠本论坛和谷歌能够99%解决。ORCA的报错信息我是真的无助啊。以前我以为是并行出错了,这 ...

可以在orca官方论坛上发帖求助,这是最靠谱的
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
本团队长期招收研究生,有意者可私信联系

227

帖子

0

威望

2800

eV
积分
3027

Level 5 (御坂)

10#
发表于 Post on 2020-8-20 10:33:28 | 只看该作者 Only view this author
ORCA出现的问题经常都是内存不够导致的,gauss的经验不能照搬

6万

帖子

99

威望

6万

eV
积分
125189

管理员

公社社长

11#
发表于 Post on 2020-8-20 23:11:56 | 只看该作者 Only view this author
ABetaCarw 发表于 2020-8-20 08:51
感谢楼上各位,问题确实去你们所说 物理内存不够。按照七楼建议 跑了一夜没有报错。
果然 计算的问题还是 ...

RG上如果是那种引用次数、发文章数目很多的人的回答,大概率靠谱。如果没怎么发过文章、没什么引用的人的回答,不要轻易相信,上面坑爹回复、瞎讨论的特别多,毕竟没人把关。

另外,可以在ORCA论坛上用搜索功能搜。ORCA论坛内容用google搜不到
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 09:44 , Processed in 0.181001 second(s), 20 queries , Gzip On.

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