计算化学公社

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

[ORCA] ORCA启动TRAH时报错

[复制链接 Copy URL]

29

帖子

0

威望

125

eV
积分
154

Level 3 能力者

我读了下手册,感觉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


3806

帖子

4

威望

7996

eV
积分
11882

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2021-11-2 14:26:57 | 只看该作者 Only view this author
贴输入文件关键词
自动做多参考态计算的程序MOKIT

1万

帖子

0

威望

7396

eV
积分
18151

Level 6 (一方通行)

3#
发表于 Post on 2021-11-2 15:48:22 | 只看该作者 Only view this author
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

29

帖子

0

威望

125

eV
积分
154

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2021-11-2 16:58:10 | 只看该作者 Only view this author
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

29

帖子

0

威望

125

eV
积分
154

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-11-2 17:06:35 | 只看该作者 Only view this author
wzkchem5 发表于 2021-11-2 15:48
参见https://orcaforum.kofo.mpg.de/viewtopic.php?f=11&t=7703

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

29

帖子

0

威望

125

eV
积分
154

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2021-11-2 17:10:03 | 只看该作者 Only view this author
我手册里发现可以把trah关了,我先关了试试,有结果再贴一波

3806

帖子

4

威望

7996

eV
积分
11882

Level 6 (一方通行)

MOKIT开发者

7#
发表于 Post on 2021-11-2 17:19:05 | 只看该作者 Only view this author
pizaipizai@5989 发表于 2021-11-2 17:10
我手册里发现可以把trah关了,我先关了试试,有结果再贴一波

跟TRAH没关系,答案就在ORCA官方论坛那个帖子里
自动做多参考态计算的程序MOKIT

1万

帖子

0

威望

7396

eV
积分
18151

Level 6 (一方通行)

8#
发表于 Post on 2021-11-2 18:26:29 | 只看该作者 Only view this author
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.
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

29

帖子

0

威望

125

eV
积分
154

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2021-11-2 21:36:45 | 只看该作者 Only view this author
zjxitcc 发表于 2021-11-2 17:19
跟TRAH没关系,答案就在ORCA官方论坛那个帖子里

感谢!

29

帖子

0

威望

125

eV
积分
154

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2021-11-2 21:37:55 | 只看该作者 Only view this author
wzkchem5 发表于 2021-11-2 18:26
填验证码需要梯子,其他应该不需要。
我贴一下那个帖子的内容吧:
Since TRAH is a second order metho ...

感谢!

29

帖子

0

威望

125

eV
积分
154

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2021-11-2 21:42:25 | 只看该作者 Only view this author
我这个问题好像不算完全解决,但是曲线救国了一下。方法是,把trah关了,然后scf第一遍不收敛,没关系,然后再restart一下scf收敛了。打算先这么凑合一下,细节以后有机会再研究

75

帖子

0

威望

1721

eV
积分
1796

Level 5 (御坂)

12#
发表于 Post on 2022-5-26 14:19:38 | 只看该作者 Only view this author
本帖最后由 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,输入文件见图片,请问是什么原因呢?谢谢!

input file.JPG (43.62 KB, 下载次数 Times of downloads: 43)

input file.JPG

TRAH-2.JPG (202.09 KB, 下载次数 Times of downloads: 31)

TRAH-2.JPG

TRAH-1.JPG (255.64 KB, 下载次数 Times of downloads: 29)

TRAH-1.JPG

1万

帖子

0

威望

7396

eV
积分
18151

Level 6 (一方通行)

13#
发表于 Post on 2022-5-26 15:22:35 | 只看该作者 Only view this author
ABQTrap 发表于 2022-5-26 07:19
您好,我用ORCA5.0.1计算有机化合物(有一个溴原子,一个磷原子,不含过渡金属)的过渡态,泛函用的wB97X-D ...

这个计算碰巧收敛比较困难,所以收敛所需时间长而已。如果不用TRAH,大概率收敛时间更长甚至不收敛,只不过因为这个例子既开了TRAH又算得慢,容易让人误以为开TRAH是计算慢的原因,但其实两者不是因果关系,相反两者都是”这个体系SCF收敛本身就比较困难“的结果。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

75

帖子

0

威望

1721

eV
积分
1796

Level 5 (御坂)

14#
发表于 Post on 2022-5-26 17:24:56 | 只看该作者 Only view this author
wzkchem5 发表于 2022-5-26 15:22
这个计算碰巧收敛比较困难,所以收敛所需时间长而已。如果不用TRAH,大概率收敛时间更长甚至不收敛,只不 ...

好的,谢谢wzkchem5 老师!

312

帖子

2

威望

2302

eV
积分
2654

Level 5 (御坂)

15#
发表于 Post on 2022-6-5 21:29:55 | 只看该作者 Only view this author
这个一个弱相互作用能的单点计算(wB97M-V/QZVPP),我想问一下出现这种情况还有必要算下去吗?

202206052128388662..png (348.62 KB, 下载次数 Times of downloads: 18)

202206052128388662..png

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 01:05 , Processed in 0.596590 second(s), 23 queries , Gzip On.

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