计算化学公社

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

[xtb] xTB调用ORCA进行几何优化时使用多核心会报错(似乎已解决)

[复制链接 Copy URL]

66

帖子

1

威望

1635

eV
积分
1721

Level 5 (御坂)

本帖最后由 Voidmio 于 2025-12-8 00:30 编辑

今使用xTB调用orca进行联合几何优化,使用命令为:/.../xtb coord --opt tight --orca -I opt.inp |tee opt-part2.out
其中 opt.inp 内容为:

opt.inp (181 Bytes, 下载次数 Times of downloads: 0)

inp 内容为:

inp (251 Bytes, 下载次数 Times of downloads: 0)

可见我在 inp 中规定了 nprocs 1, 这时上述命令是可以正确结束的(输出内容见 opt-part2.out1.txt),可见 xTB 对 ORCA 每轮的单点结果进行了校验,并最终收敛:
opt-part2.out1.txt (429.42 KB, 下载次数 Times of downloads: 0)

但如果我把这个 nprocs 1 改为  nprocs 32(不是因为太多,2 核心也是相同的报错),就会在第一次 ORCA 任务结束后的校验时,xTB报错(输出内容见 opt-part2.out32.txt):
opt-part2.out32.txt (47.67 KB, 下载次数 Times of downloads: 0)
(这是 nprocs 32 的输出末尾)

可见这里第一轮 ORCA 是正确结束的,其不存在 omp 的调度问题(多核心单独运行 ORCA 也可以正确结束);反倒是似乎无关的 xTB 报错了。并且这里nprocs 32第一轮ORCA的输出和nprocs 1应该是完全相同的:
(这是 nprocs 1 的第一次 ORCA 结束,可见 ORCA TERMINATED NORMALLY 后就回到 xTB 了)

环境里写有
export OMP_NUM_THREADS=64
export MKL_NUM_THREADS=64
export OMP_STACKSIZE=4G

均没有改动。单独运行 ORCA 和 xTB 应该都没问题,ORCA、OpenMPI、xTB 的版本分别是5.0.1,4.1.1,6.7.1(前两者也试过6.1.0、4.1.6,没区别)

求教各位老师,这是怎么回事呢





536

帖子

1

威望

5727

eV
积分
6283

Level 6 (一方通行)

2#
发表于 Post on 2025-12-7 21:28:18 | 只看该作者 Only view this author
直接在orca 6.1中用! opt freq xtb2就行了。ORCA 6.1中已经内置了xtb。

66

帖子

1

威望

1635

eV
积分
1721

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2025-12-7 21:56:08 | 只看该作者 Only view this author
niobium 发表于 2025-12-7 21:28
直接在orca 6.1中用! opt freq xtb2就行了。ORCA 6.1中已经内置了xtb。

不中啊,我这个是grimme小玩具中的一环,后续还有程序要读取这次xtb的输出内容,非常头大

66

帖子

1

威望

1635

eV
积分
1721

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2025-12-8 00:32:23 | 只看该作者 Only view this author
现在看来,应该是ompi使用btl时候会莫名其妙warning,其实没报错但给xtb干懵了。
不是我喜欢的字节传输层,直接禁用:export OMPI_MCA_btl=^openib
这样就可以run了,但是不知道有没有什么坏处

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 03:17 , Processed in 0.172475 second(s), 23 queries , Gzip On.

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