计算化学公社

标题: 请教大纯有机体系优化建议（ORCA的SCF不收敛） [打印本页]

作者
Author: ABetaCarw 时间: 2020-8-19 19:33
标题: 请教大纯有机体系优化建议（ORCA的SCF不收敛）
本帖最后由 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

然后发现每次都报错

************************************************************
* Program running with 36 parallel MPI-processes *
* working on a common directory *
************************************************************
------------------------------------------------------------------------------
ORCA SCF GRADIENT CALCULATION
------------------------------------------------------------------------------
Gradient of the Kohn-Sham DFT energy:
Kohn-Sham wavefunction type ... RKS
Hartree-Fock exchange scaling ... 0.000
Number of operators ... 1
Number of atoms ... 453
Basis set dimensions ... 4464
Integral neglect threshold ... 2.5e-11
Integral primitive cutoff ... 2.5e-12
Nuclear repulsion gradient ... done
One Electron Gradient ... done
Pre-screening matrix ... done
RI-J gradient ...
ORCA finished by error termination in SCF gradient
Calling Command: mpirun -np 36 /home/XXX/ORCA/orca421/orca_scfgrad_mpi comp2.scfgrad.inp comp2
[file orca_tools/qcmsg.cpp, line 458]:
.... aborting the run

复制代码

请熟悉ORCA的前辈　指导一下问题出在哪里。

作者
Author: 喵星大佬 时间: 2020-8-19 19:39
求梯度内存不够吧。。。。可以降低并行核数，再就是orca的几何优化其实很蠢的，我算了一个很小的配合物，优化完了振动分析3个虚频。。。。

作者
Author: sobereva 时间: 2020-8-19 20:14
这和SCF不收敛无关，纯粹是受力计算部分的问题

比起BLYP-D3-gCP，我明显更推荐B97-3c，精度更好，耗时差不多
B97-3c也这样的话，加大内存，如果物理内存不够，减少核数再试

作者
Author: ABetaCarw 时间: 2020-8-19 23:59
高斯的报错信息我靠本论坛和谷歌能够99%解决。ORCA的报错信息我是真的无助啊。以前我以为是并行出错了，这次我判定SCF报错的，我以为我找对了……感谢社长和喵星大佬的支持

作者
Author: snljty 时间: 2020-8-20 00:22

ABetaCarw 发表于 2020-8-19 23:59
高斯的报错信息我靠本论坛和谷歌能够99%解决。ORCA的报错信息我是真的无助啊。以前我以为是并行出错了，这 ...

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

作者
Author: winnerwill 时间: 2020-8-20 01:04
体系这么大！%maxcore不给够肯定不行，估计要放大5～10倍

作者
Author: qinzhong605 时间: 2020-8-20 01:31
核数不超过10，内存单核16-20 GB试试，我也刚遇到过类似问题。就是这么解决的

作者
Author: ABetaCarw 时间: 2020-8-20 08:51
本帖最后由 ABetaCarw 于 2020-8-20 09:51 编辑

感谢楼上各位，问题确实去你们所说物理内存不够。按照七楼建议跑了一夜没有报错。
果然计算的问题还是要向靠谱的人求助。我在researchgate上面找了两个讨论该报错的帖子全是在讨论改收敛精度改SCF方法，太坑了。

作者
Author: wzkchem5 时间: 2020-8-20 09:50

ABetaCarw 发表于 2020-8-19 23:59
高斯的报错信息我靠本论坛和谷歌能够99%解决。ORCA的报错信息我是真的无助啊。以前我以为是并行出错了，这 ...

可以在orca官方论坛上发帖求助，这是最靠谱的

作者
Author: shalene 时间: 2020-8-20 10:33
ORCA出现的问题经常都是内存不够导致的，gauss的经验不能照搬

作者
Author: sobereva 时间: 2020-8-20 23:11

ABetaCarw 发表于 2020-8-20 08:51
感谢楼上各位，问题确实去你们所说物理内存不够。按照七楼建议跑了一夜没有报错。
果然计算的问题还是 ...

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

另外，可以在ORCA论坛上用搜索功能搜。ORCA论坛内容用google搜不到

欢迎光临计算化学公社 (http://bbs.keinsci.com/)

Powered by Discuz! X3.3