计算化学公社

标题: ORCA启动TRAH时报错 [打印本页]

作者
Author: pizaipizai@5989 时间: 2021-11-2 13:23
标题: ORCA启动TRAH时报错
我读了下手册，感觉trah是orca用来加速收敛的技术，不过程序准备使用trah时，发生了一个错误，如下。我想求助这个错误是我的orca安装出了问题吗？orca是从帖子各个版本的逆战鲸下载的，501版本linux。
TRAH GRID
---------

General Integration Accuracy    IntAcc    ... 3.467
Radial Grid Type                RadialGrid  ... OptM3 with GC (2021)
Angular Grid (max. ang.)       AngularGrid ... 1 (Lebedev-50)
Angular grid pruning method    GridPruning ... 4 (adaptive)
Weight generation scheme       WeightScheme... Becke
Basis function cutoff          BFCut    ... 1.0000e-11
Integration weight cutoff       WCut       ... 1.0000e-14
Angular grids for H and He will be reduced by one unit
Partially contracted basis set             ... off
Rotationally invariant grid construction    ... off

Total number of grid points                ... 21074
Total number of batches                   ...    338
Average number of points per batch          ...    62
Average number of grid points per atom    ...    878

--------------------------------------------------------------------------------------------
Iter.       energy          ||Error||_2       Shift    TRadius  Mac/Mic       Rej.
--------------------------------------------------------------------------------------------
[file orca_tools/Tool-Numint/qcfxc.cpp, line 14986, Process 41]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

[file orca_tools/Tool-Numint/qcfxc.cpp, line 14986, Process 3]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

[file orca_tools/Tool-Numint/qcfxc.cpp, line 14986, Process 19]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

[file orca_tools/Tool-Numint/qcfxc.cpp, line 14986, Process 13]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

[file orca_tools/Tool-Numint/qcfxc.cpp, line 14986, Process 9]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

[file orca_tools/Tool-Numint/qcfxc.cpp, line 14986, Process 53]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

[file orca_tools/Tool-Numint/qcfxc.cpp, line 14986, Process 49]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

[file orca_tools/Tool-Numint/qcfxc.cpp, line 14986, Process 1]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

[file orca_tools/Tool-Numint/qcfxc.cpp, line 14986, Process 5]: Error: Invalid or unknown value for Exchange in DFT XC-Kernel

作者
Author: zjxitcc 时间: 2021-11-2 14:26
贴输入文件关键词

作者
Author: wzkchem5 时间: 2021-11-2 15:48
参见https://orcaforum.kofo.mpg.de/viewtopic.php?f=11&t=7703

作者
Author: pizaipizai@5989 时间: 2021-11-2 16:58

zjxitcc 发表于 2021-11-2 14:26
贴输入文件关键词

1> ! PWPB95 D3 def2-QZVPP def2/J def2-QZVPP/C RIJCOSX tightSCF noautostart miniprint nopop
2> %maxcore 4000
3> %pal nprocs 64 end
4> * xyz 0 1

作者
Author: pizaipizai@5989 时间: 2021-11-2 17:06

wzkchem5 发表于 2021-11-2 15:48
参见https://orcaforum.kofo.mpg.de/viewtopic.php?f=11&t=7703

感谢！
不过我的qq，校园还有特意新注册的yeah邮箱，人家orca注册都觉得这些邮箱无效。也许需要梯子？

作者
Author: pizaipizai@5989 时间: 2021-11-2 17:10
我手册里发现可以把trah关了，我先关了试试，有结果再贴一波

作者
Author: zjxitcc 时间: 2021-11-2 17:19

pizaipizai@5989 发表于 2021-11-2 17:10
我手册里发现可以把trah关了，我先关了试试，有结果再贴一波

跟TRAH没关系，答案就在ORCA官方论坛那个帖子里

作者
Author: wzkchem5 时间: 2021-11-2 18:26

pizaipizai@5989 发表于 2021-11-2 10:06
感谢！
不过我的qq，校园还有特意新注册的yeah邮箱，人家orca注册都觉得这些邮箱无效。也许需要梯子？

填验证码需要梯子，其他应该不需要。
我贴一下那个帖子的内容吧：
Since TRAH is a second order method, it computes the electronic Hessian and that requires the Kernel (second derivative of the XC functional)
However, this specific kernel for that functional is not natively implemented. You should try using LibXC.
In case that also fails, there is the
! NoTRAH
keyword.

作者
Author: pizaipizai@5989 时间: 2021-11-2 21:36

zjxitcc 发表于 2021-11-2 17:19
跟TRAH没关系，答案就在ORCA官方论坛那个帖子里

感谢！

作者
Author: pizaipizai@5989 时间: 2021-11-2 21:37

wzkchem5 发表于 2021-11-2 18:26
填验证码需要梯子，其他应该不需要。
我贴一下那个帖子的内容吧：
Since TRAH is a second order metho ...

感谢！

作者
Author: pizaipizai@5989 时间: 2021-11-2 21:42
我这个问题好像不算完全解决，但是曲线救国了一下。方法是，把trah关了，然后scf第一遍不收敛，没关系，然后再restart一下scf收敛了。打算先这么凑合一下，细节以后有机会再研究

作者
Author: ABQTrap 时间: 2022-5-26 14:19
本帖最后由 ABQTrap 于 2022-5-26 14:22 编辑

您好，我用ORCA5.0.1计算有机化合物（有一个溴原子，一个磷原子，不含过渡金属）的过渡态，泛函用的wB97X-D3，先用较小基组ma-def2-SVP优化，很顺利就找到过渡态了，收敛没有问题。然后把ma-def2-SVP下的GBW和hess文件作为初猜和频率文件，打算在大基组ma-def2-TZVPP下优化，但是遇到了收敛的问题。用默认的方法时输出文件会显示“WARNING: the maximum gradient error descreased on average only by a factor 1.1 during the last 20 iterations *** Initiating the TRAH-SCF procedure **”。然后启动TRAH，这个计算非常慢，虽然最后收敛了。
我之前计算的其他有机化合物（不含磷原子和过渡金属，但都有溴原子）均不存在这个问题，积分格点是defgrid3，输入文件见图片，请问是什么原因呢？谢谢！

作者
Author: wzkchem5 时间: 2022-5-26 15:22

ABQTrap 发表于 2022-5-26 07:19
您好，我用ORCA5.0.1计算有机化合物（有一个溴原子，一个磷原子，不含过渡金属）的过渡态，泛函用的wB97X-D ...

这个计算碰巧收敛比较困难，所以收敛所需时间长而已。如果不用TRAH，大概率收敛时间更长甚至不收敛，只不过因为这个例子既开了TRAH又算得慢，容易让人误以为开TRAH是计算慢的原因，但其实两者不是因果关系，相反两者都是”这个体系SCF收敛本身就比较困难“的结果。

作者
Author: ABQTrap 时间: 2022-5-26 17:24

wzkchem5 发表于 2022-5-26 15:22
这个计算碰巧收敛比较困难，所以收敛所需时间长而已。如果不用TRAH，大概率收敛时间更长甚至不收敛，只不 ...

好的，谢谢wzkchem5 老师！

作者
Author: 七尺贱 时间: 2022-6-5 21:29
这个一个弱相互作用能的单点计算（wB97M-V/QZVPP），我想问一下出现这种情况还有必要算下去吗？

作者
Author: wzkchem5 时间: 2022-6-5 22:26

七尺贱发表于 2022-6-5 14:29
这个一个弱相互作用能的单点计算（wB97M-V/QZVPP），我想问一下出现这种情况还有必要算下去吗？

应该是有基组线性相关问题。把这个计算停掉，在%scf里设置sthresh 1e-7（如果还不管用就1e-6），然后重试

作者
Author: 七尺贱 时间: 2022-6-5 22:37

wzkchem5 发表于 2022-6-5 22:26
应该是有基组线性相关问题。把这个计算停掉，在%scf里设置sthresh 1e-7（如果还不管用就1e-6），然后重试

谢谢您，我试试

作者
Author: 七尺贱 时间: 2022-6-9 11:54
本帖最后由七尺贱于 2022-6-9 11:56 编辑

wzkchem5 发表于 2022-6-5 22:26
应该是有基组线性相关问题。把这个计算停掉，在%scf里设置sthresh 1e-7（如果还不管用就1e-6），然后重试

老师您好，我已经按照您的方法改成1e-6了，但还是出现这个问题，我看着手册说线性差值小于1.125就会有这个问题，我这个1.1好像也差得不是太多，应该可以继续算下去把
D:\desktop\图片1.png

作者
Author: wzkchem5 时间: 2022-6-9 15:49

七尺贱发表于 2022-6-9 04:54
老师您好，我已经按照您的方法改成1e-6了，但还是出现这个问题，我看着手册说线性差值小于1.125就会有这个 ...

你截图不是报错，注意TRAH启动是完全正常的，只有TRAH不收敛才需要干预。建议去了解一下TRAH是做什么的，这样就知道为什么是这样了。
此外，1.125不是什么“线性差值”，根本没有“线性差值”这个概念，而且这个1.1也不是为了和1.125比的。我不知道你是从哪里看来的1.125，不要看到一个1.125、又因为1.1和1.125接近，就以为应该比较这两者。两者根本没关系。

作者
Author: 七尺贱 时间: 2022-6-9 16:46

wzkchem5 发表于 2022-6-9 15:49
你截图不是报错，注意TRAH启动是完全正常的，只有TRAH不收敛才需要干预。建议去了解一下TRAH是做什么的， ...

对不起老师，是我表述不清楚，TRAH是ORCA让SCF收敛的一种方法。这个是我手册9.7.7里面找到的话，可能我没理解。In cases for which the conventional SCF procedures (DIIS/KDIIS/SOSCF) struggle, we invoke TRAH-SCF automatically (AutoTRAH). For this purpose, we perform a linear interpolation of the norm of the electronic gradient on the log10 scale after a minimum number of SCF iterations AutoTRAHIter (default 20). The number of iterations for interpolations is controlled by AutoTRAHNInter (default 10). If the slope s of the interpolated gradient norm 10-s becomes smaller than AutoTRAHTol (default 1.125), conventional SCF is shut down and TRAH-SCF start from the current set of MOs. Those parameters were optimized for a benchmark set with the purpose to minimize the calculation times even for SCF calculations that are hard to converge.

作者
Author: wzkchem5 时间: 2022-6-9 17:50

七尺贱发表于 2022-6-9 09:46
对不起老师，是我表述不清楚，TRAH是ORCA让SCF收敛的一种方法。这个是我手册9.7.7里面找到的话，可能我没 ...

这里说的是线性外推，和线性相关虽然都有线性两个字，但是半毛钱关系都没有。
建议去看看线性代数课本，了解什么是线性相关，以及它为什么和线性外推没关系。

作者
Author: 七尺贱 时间: 2022-6-9 20:36

wzkchem5 发表于 2022-6-9 17:50
这里说的是线性外推，和线性相关虽然都有线性两个字，但是半毛钱关系都没有。
建议去看看线性代数课本， ...

谢谢老师

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