ORCA算单点能SCF不收敛

Oliviaw · 发表于 Post on 2024-2-21 09:11:26

我的输入行是! 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?)

请问该怎么重新设置？

Oliviaw · 发表于 Post on 2024-2-21 09:12:46

计算的分子是148个原子。我用def2-tzvp+gcp跑过能收敛，但是用def2-QZVP就不收敛。

Strange · 发表于 Post on 2024-2-21 09:44:04

你读一下收敛过的波函数当初猜呗

Oliviaw · 发表于 Post on 2024-2-21 13:46:51

Strange 发表于 2024-2-21 09:44
你读一下收敛过的波函数当初猜呗

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

lucky1999 · 发表于 Post on 2024-2-21 14:02:45

Oliviaw 发表于 2024-2-21 13:46
请问读取gbw文件的关键词是什么呀？

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

zjxitcc · 发表于 Post on 2024-2-21 15:23:49

Oliviaw 发表于 2024-2-21 13:46
请问读取gbw文件的关键词是什么呀？

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

sobereva · 发表于 Post on 2024-2-21 15:36:11

Oliviaw 发表于 2024-2-21 09:12
计算的分子是148个原子。我用def2-tzvp+gcp跑过能收敛，但是用def2-QZVP就不收敛。

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

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

也可以尝试defgrid3，用KDIIS等其它做法，本论坛里有很多讨论

Oliviaw · 发表于 Post on 2024-2-21 17:34:41

sobereva 发表于 2024-2-21 15:36
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员 ...

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

sobereva · 发表于 Post on 2024-2-21 18:38:04

Oliviaw 发表于 2024-2-21 17:34
请问orca3.0.1的grid=5和grid=6，与orca5.0.4的defgrid3是一个档次吗？

没有严格的可比性
只能说grid6和defgrid3的格点质量都已经非常高了

Oliviaw · 发表于 Post on 2024-2-21 19:31:57

sobereva 发表于 2024-2-21 18:38
没有严格的可比性
只能说grid6和defgrid3的格点质量都已经非常高了

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

Oliviaw · 发表于 Post on 2024-2-21 19:57:45

zjxitcc 发表于 2024-2-21 15:23
不需要添加任何关键词，只要inp和gbw文件同名，程序自动会读入轨道和基组，投影到当前目标基组上。换句话 ...

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

wzkchem5 · 发表于 Post on 2024-2-21 20:23:59

Oliviaw 发表于 2024-2-21 12:31
一个很基础的问题，积分格点的设置的精度高是解决了什么问题，为什么要设置defgrid3？

格点精度高，数值积分的精度就高，数值误差就小。有的时候SCF不收敛是因为数值误差太大，计算出来的能量和真实能量之间有误差，虽然误差的幅度小，但是符号随机，导致波函数一直小幅度震荡。这时增大格点就是最有效的方法。
参见https://bdf-manual.readthedocs.i ... tml#scfconvproblems

sobereva · 发表于 Post on 2024-2-21 20:34:32

Oliviaw 发表于 2024-2-21 19:31
一个很基础的问题，积分格点的设置的精度高是解决了什么问题，为什么要设置defgrid3？

积分格点质量不够高可能导致KS矩阵构造精度差，对于某些情况可能导致SCF难收敛。例如之前我用wB97M-V结合很大基组计算18碳环的弱相互作用复合物就遇到这种情况，默认的defgrid2 SCF震荡，改成defgrid3收敛了。

Oliviaw · 发表于 Post on 2024-2-21 21:17:35

wzkchem5 发表于 2024-2-21 20:23
格点精度高，数值积分的精度就高，数值误差就小。有的时候SCF不收敛是因为数值误差太大，计算出来的能量 ...

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

Oliviaw · 发表于 Post on 2024-2-21 21:25:53

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呢？

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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