计算化学公社

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

[综合交流] sobEDA计算过程中formchk和unfchk工具进行chk和fch之间的转换问题

[复制链接 Copy URL]

27

帖子

0

威望

235

eV
积分
262

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 ZHUOHAO 于 2025-6-22 11:11 编辑

各位老师,请教一下sobEDA计算过程中formchk和unfchk工具进行chk和fch之间的转换不成功。体系包含了444个原子分成了4个片段去计算。
第一次计算chk和fch转换不成功,错误提示为
Generating fch file of promolecular state via Multiwfn (promol.fch)
Running: unfchk promol.fch promol.chk
Error: segmentation violation
   rax 0000000000000000, rbx 00000000008ceac8, rcx ffffffffffffffff
   rdx 0000000000000f1d, rsp 00007ffe4104d7c8, rbp 00007ffe4104d820
   rsi 000000000000000b, rdi 0000000000000f1d, r8  00002ae44ff28e00
   r9  00002ae44ff28e00, r10 00007ffe4104d220, r11 0000000000000202
   r12 00007ffe4104daa0, r13 00000000008ced21, r14 0000000000000003
   r15 00007ffe4104efe8
  /lib64/libpthread.so.0(+0xf5d0) [0x2ae45010a5d0]
  /lib64/libc.so.6(kill+0x7) [0x2ae45064f4d7]
  unfchk() [0x436020]
  unfchk() [0x41cfd3]
  unfchk() [0x42bc2f]
  unfchk() [0x42d496]
  unfchk() [0x404752]
  unfchk() [0x403839]
  unfchk() [0x40379d]
  /lib64/libc.so.6(__libc_start_main+0xf5) [0x2ae45063b3d5]
  unfchk(sched_setaffinity+0xb9) [0x4036a9]
./sobEDA.sh: 行 159:  3869 已放弃               unfchk promol.fch promol.chk > /dev/null
Generating Gaussian input file of promolecular state via Multiwfn (promol.gjf)
Running: g16 < promol.gjf &> promol.out
Error encountered, please check corresponding output file! Now script exits


根据sob老师提供的解决方案:在用户主目录下的.bashrc里加上比如export GAUSS_MEMDEF=20GB然后重新进入终端。
修改后第二次计算仍然出现计算chk和fch转换不成功,错误提示为:
Generating fch file of promolecular state via Multiwfn (promol.fch)
Running: unfchk promol.fch promol.chk
Error: segmentation violation
   rax 0000000000000000, rbx 00000000008ceac8, rcx ffffffffffffffff
   rdx 000000000002ebd7, rsp 00007ffe71e31478, rbp 00007ffe71e314d0
   rsi 000000000000000b, rdi 000000000002ebd7, r8  00002b2962848e00
   r9  00002b2962848e00, r10 00007ffe71e30ee0, r11 0000000000000206
   r12 00007ffe71e31750, r13 00000000008ced21, r14 0000000000000003
   r15 00007ffe71e32c98
  /lib64/libpthread.so.0(+0xf5d0) [0x2b2962a2a5d0]
  /lib64/libc.so.6(kill+0x7) [0x2b2962f6f4d7]
  unfchk() [0x436020]
  unfchk() [0x41cfd3]
  unfchk() [0x42bc2f]
  unfchk() [0x42d496]
  unfchk() [0x404752]
  unfchk() [0x403839]
  unfchk() [0x40379d]
  /lib64/libc.so.6(__libc_start_main+0xf5) [0x2b2962f5b3d5]
  unfchk(sched_setaffinity+0xb9) [0x4036a9]
./sobEDA.sh: 行 159: 191447 已放弃               unfchk promol.fch promol.chk > /dev/null
Generating Gaussian input file of promolecular state via Multiwfn (promol.gjf)
Running: g16 < promol.gjf &> promol.out
Error encountered, please check corresponding output file! Now script exits

输出的promol.fch文件有轨道信息。看其他帖子sob老师回复可以手动用unfchk,得到chk文件是继续放到当前计算目录下然后重新计算还是其他什么操作

还有个办法是现在的情况是需要继续增大内存重新进入终端吗?如果是继续增大内存,增大为多少合适呢?或者是 还有其他解决办法?




6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

2#
发表于 Post on 2025-6-23 05:17:39 | 只看该作者 Only view this author
先按照下文的FAQ部分排查,避免存在其它可能的问题
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685http://bbs.keinsci.com/thread-39446-1-1.html

有多大空余物理内存GAUSS_MEMDEF就可以尝试设多大

尝试debug而自己unfchk的话,也放到当前目录
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

27

帖子

0

威望

235

eV
积分
262

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-6-23 08:54:20 | 只看该作者 Only view this author
sobereva 发表于 2025-6-23 05:17
先按照下文的FAQ部分排查,避免存在其它可能的问题
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适 ...

好的,谢谢sob老师回复

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 04:45 , Processed in 0.180854 second(s), 20 queries , Gzip On.

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