计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: hdhxx123
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] ORCA 5 ,Yes

[复制链接 Copy URL]

1236

帖子

1

威望

3495

eV
积分
4751

Level 6 (一方通行)

76#
发表于 Post on 2021-7-6 16:36:33 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-6 00:08
ORCA 5发布会部分PPT和视频已开放下载:https://orcaforum.kofo.mpg.de/app.php/dlext/?cat=15
注:如下载 ...

我发现了我的电脑上装了orca5以后高斯跑scf迭代轮数很多的时候莫名其妙的会死掉,注释掉以后发现就好了,怀疑是openmpi的问题
不知道有没有别人遇到过这些问题

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

77#
发表于 Post on 2021-7-6 16:41:34 | 只看该作者 Only view this author
biogon 发表于 2021-7-6 09:36
我发现了我的电脑上装了orca5以后高斯跑scf迭代轮数很多的时候莫名其妙的会死掉,注释掉以后发现就好了, ...

你是说把什么注释掉?以及这个问题可以重复吗?
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

1236

帖子

1

威望

3495

eV
积分
4751

Level 6 (一方通行)

78#
发表于 Post on 2021-7-7 20:16:52 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-6 16:41
你是说把什么注释掉?以及这个问题可以重复吗?

把openmpi环境变量注释掉就没问题了,不注释掉就反复出现

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

79#
发表于 Post on 2021-7-7 20:54:02 | 只看该作者 Only view this author
biogon 发表于 2021-7-7 13:16
把openmpi环境变量注释掉就没问题了,不注释掉就反复出现

那你把设置openmpi环境变量的语句写到调用orca的脚本里就行了,不写到bashrc里,这样就不会影响到高斯
好像orca手册也是建议这样做的
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5万

帖子

99

威望

5万

eV
积分
112349

管理员

公社社长

80#
发表于 Post on 2021-7-7 21:17:48 | 只看该作者 Only view this author
ene 发表于 2021-7-4 15:15
再说一个有关orca_md的问题,不知道为什么运行步数超过1000,任务就自动停掉了……任务停下时会报这样一个 ...

我也发现有bug,每次跑到1203步时候必死,试了几种体系、各种设置,现象都一样。已经发到官方论坛
https://orcaforum.kofo.mpg.de/viewtopic.php?f=11&p=32944#p32944

评分 Rate

参与人数
Participants 2
eV +8 收起 理由
Reason
captain + 5 赞!
ene + 3 好!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

41

帖子

0

威望

3874

eV
积分
3915

Level 5 (御坂)

81#
发表于 Post on 2021-7-8 09:16:59 | 只看该作者 Only view this author
sobereva 发表于 2021-7-7 21:17
我也发现有bug,每次跑到1203步时候必死,试了几种体系、各种设置,现象都一样。已经发到官方论坛
https ...

我也是这种情况,都是1137步就死
select: Bad file descriptor
[file orca_tools/qcmsg.cpp, line 458]:
  .... aborting the run

1236

帖子

1

威望

3495

eV
积分
4751

Level 6 (一方通行)

82#
发表于 Post on 2021-7-8 11:20:33 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-7 20:54
那你把设置openmpi环境变量的语句写到调用orca的脚本里就行了,不写到bashrc里,这样就不会影响到高斯
...

那只能先这么用了,openmpi4.0.5的时候都没有遇到这种问题

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

83#
发表于 Post on 2021-7-8 15:30:06 | 只看该作者 Only view this author
asdf 发表于 2021-7-8 02:16
我也是这种情况,都是1137步就死
select: Bad file descriptor
[file orca_tools/qcmsg.cpp, line 458] ...

这个问题ORCA 5.0.1会修复。本质上是orca_md有一个文件每步MD都会打开,但是忘了关
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

433

帖子

0

威望

2195

eV
积分
2628

Level 5 (御坂)

娃娃儿鱼

84#
 楼主 Author| 发表于 Post on 2021-7-8 17:07:41 | 只看该作者 Only view this author
本帖最后由 hdhxx123 于 2021-7-8 17:32 编辑

TRAH应该是存在bug还是啥的,会收敛着就卡在一个能量上,而显示的能量变化又是在一直降低的,但实际是在震荡。DIIS的veryslowconv特别迅速地就收敛了,没有明显的震荡。另外需要注意的是,TRAH与slowconv及veryslowconv并不会一起用,如果要一直用slowconv或veryslowconv的DIIS时记得用noTRAH把TRAH关掉

temp.zip

10.63 KB, 下载次数 Times of downloads: 8

Gd1_Gau.gbw

749.74 KB, 下载次数 Times of downloads: 0

真·探

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

85#
发表于 Post on 2021-7-8 20:53:25 | 只看该作者 Only view this author
hdhxx123 发表于 2021-7-8 10:07
TRAH应该是存在bug还是啥的,会收敛着就卡在一个能量上,而显示的能量变化又是在一直降低的,但实际是在震 ...

注意一个问题:你把thresh和tcut设得比默认值还松(使用verytightscf时, thresh默认为1e-12,tcut默认为1e-14)。。。这么松的thresh和tcut应该是不足以收敛到verytightscf的标准的。我去掉thresh和tcut的设定,并打开TRAH后,SCF 89步收敛,收敛能量-764.036712959977。你另一个帖子算Tb也有这个问题,本想把积分阈值设严,结果比默认还松,当然不收敛。。。
TRAH能量降低又跳回原值不是bug,是因为TRAH有macroiteration和microiteration,每步macroiteration的最后一步microiteration的能量可以比下一步macroiteration的第一步microiteration低。参见TRAH原始文献https://aip.scitation.org/doi/pdf/10.1063/5.0040798
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

433

帖子

0

威望

2195

eV
积分
2628

Level 5 (御坂)

娃娃儿鱼

86#
 楼主 Author| 发表于 Post on 2021-7-8 21:01:03 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-8 20:53
注意一个问题:你把thresh和tcut设得比默认值还松(使用verytightscf时, thresh默认为1e-12,tcut默认为 ...

了解了,谢谢老师
真·探

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

87#
发表于 Post on 2021-7-8 21:07:03 | 只看该作者 Only view this author
hdhxx123 发表于 2021-7-8 10:07
TRAH应该是存在bug还是啥的,会收敛着就卡在一个能量上,而显示的能量变化又是在一直降低的,但实际是在震 ...

另外顺便一提,TRAH和slowconv/veryslowconv不能结合使用是feature而不是bug。一般只有在能量振荡的时候才会用到slowconv/veryslowconv,而TRAH是自带linear search的,如果能量升高,会自动退回到上一步波函数并缩小步长,也就是该slow的时候才slow,发现上一步迈大了就撤销上一步然后重新走一小步,不需要slow的时候就放开步子走。所以,发现TRAH和slowconv不会一起用就关掉TRAH是没有任何道理的。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

433

帖子

0

威望

2195

eV
积分
2628

Level 5 (御坂)

娃娃儿鱼

88#
 楼主 Author| 发表于 Post on 2021-7-8 21:08:41 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-8 21:07
另外顺便一提,TRAH和slowconv/veryslowconv不能结合使用是feature而不是bug。一般只有在能量振荡的时候 ...

但是我尝试了一下发现veryslowconv的DIIS速度比TRAH的速度收敛快的多。。在波函数质量有保证的前提下,如果DIIS能收敛的话可以选择关掉TRAH
真·探

1万

帖子

0

威望

7390

eV
积分
18143

Level 6 (一方通行)

89#
发表于 Post on 2021-7-8 22:07:45 | 只看该作者 Only view this author
hdhxx123 发表于 2021-7-8 14:08
但是我尝试了一下发现veryslowconv的DIIS速度比TRAH的速度收敛快的多。。在波函数质量有保证的前提下,如 ...

那就是程序自动开TRAH的阈值对于这个体系太松了
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

433

帖子

0

威望

2195

eV
积分
2628

Level 5 (御坂)

娃娃儿鱼

90#
 楼主 Author| 发表于 Post on 2021-7-8 22:25:20 | 只看该作者 Only view this author
wzkchem5 发表于 2021-7-8 22:07
那就是程序自动开TRAH的阈值对于这个体系太松了

老师,想请问一下Tcut和Thresh在orca4.2.1中的设置与scf级别的设置是怎么样的呢?是orca5中有改变吗?我在orca4.2.1手册的9.6节SCF convergence中没有找到具体的定义
真·探

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 05:28 , Processed in 0.216618 second(s), 24 queries , Gzip On.

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