计算化学公社

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

[Gaussian/gview] 单点任务出错Error: segmentation violation

[复制链接 Copy URL]

62

帖子

0

威望

379

eV
积分
441

Level 3 能力者

本帖最后由 Sunca 于 2025-1-24 08:03 编辑

Gaussian单点任务时报错界面如下:
g16 com.gjf
Error: segmentation violation
   rax 0000000000000000, rbx ffffffffffffffff, rcx 000009bc7644275b
   rdx 0000000000000000, rsp 00007ffe99d94878, rbp 00007ffe99d94df0
   rsi 000000000000000b, rdi 000000000000cd8c, r8  000009bb763ff000
   r9  0000000000000401, r10 000009bc76412f60, r11 0000000000000206
   r12 00007ffe99d94e38, r13 0000000000000000, r14 0000000000000000
   r15 00000000000003e6
  /lib/x86_64-linux-gnu/libc.so.6(+0x42520) [0x9bc76442520]
  /lib/x86_64-linux-gnu/libc.so.6(kill+0xb) [0x9bc7644275b]
  /home/sunca/Gauss16/g16/l302.exe() [0x4755d9]
  /home/sunca/Gauss16/g16/l302.exe() [0x48e5aa]
  /home/sunca/Gauss16/g16/l302.exe() [0x4e50bc]
  /home/sunca/Gauss16/g16/l302.exe() [0x56f4e2]
  /home/sunca/Gauss16/g16/l302.exe() [0x549498]
  /home/sunca/Gauss16/g16/l302.exe() [0x500851]
  /home/sunca/Gauss16/g16/l302.exe() [0x4599c9]
  /home/sunca/Gauss16/g16/l302.exe() [0x454389]
  /home/sunca/Gauss16/g16/l302.exe() [0x450beb]
  /home/sunca/Gauss16/g16/l302.exe() [0x450af4]
  /lib/x86_64-linux-gnu/libc.so.6(+0x29d90) [0x9bc76429d90]
  /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0x80) [0x9bc76429e40]
  /home/sunca/Gauss16/g16/l302.exe() [0x4509e9]
Aborted
另外上传了gjy和log文件,我用的基组文件ma-TZVPP.txt来自http://sobereva.com/509。希望有大佬回复


加了下面两句后跑起来了,暂时没报错。可能是我原子太多了?然后我520个原子大概要跑多久啊?
%mem=32000MB
%nprocshared=8

com.gjf

23.95 KB, 下载次数 Times of downloads: 3

com.log

327.46 KB, 下载次数 Times of downloads: 6

72

帖子

0

威望

637

eV
积分
709

Level 4 (黑子)

2#
发表于 Post on 2025-1-24 09:13:34 | 只看该作者 Only view this author
来教你

先看这个http://bbs.keinsci.com/thread-4829-1-1.html

然后打开你的log文件,Ctrl+F termination ,找到

  1. Out-of-memory error in routine ShPair-CalcS2-2 (IEnd=     373609315 MxCore=     316521421)
  2. Use %mem=615MW to provide the minimum amount of memory required to complete this step.
  3. Error termination via Lnk1e in /home/sunca/Gauss16/g16/l302.exe at Fri Jan 24 07:38:48 2025.
复制代码


65

帖子

0

威望

533

eV
积分
598

Level 4 (黑子)

3#
发表于 Post on 2025-1-24 10:18:01 | 只看该作者 Only view this author
好奇楼主是什么样的牛逼机器,520原子的体系,用TZVPP这么大的基组,估计没认真学过卢老师的贴子,这么大体系为什么要这么大的基组,为什么要加弥散?

13176 basis functions, 19632 primitive gaussians, 14840 cartesian basis functions
  1130 alpha electrons     1130 beta electrons

2260个电子,13176个基函数,Oh my God, 这得用啥服务器可以算出来,8个CPU够吗?

62

帖子

0

威望

379

eV
积分
441

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2025-1-24 12:26:19 | 只看该作者 Only view this author
本帖最后由 Sunca 于 2025-1-24 12:29 编辑
savorw 发表于 2025-1-24 10:18
好奇楼主是什么样的牛逼机器,520原子的体系,用TZVPP这么大的基组,估计没认真学过卢老师的贴子,这么大体 ...

我刚入门,想做单点任务得到fch文件然后算氢键结合能,参考的http://sobereva.com/522例子中用的基组,我不太懂,您推荐我应该用什么基组?

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

5#
发表于 Post on 2025-1-24 14:31:44 | 只看该作者 Only view this author
Sunca 发表于 2025-1-24 12:26
我刚入门,想做单点任务得到fch文件然后算氢键结合能,参考的http://sobereva.com/522例子中用的基组,我 ...

恰当简化体系,只保留要用522博文里的做法算氢键键能附近的原子
北京科音自然科学研究中心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

62

帖子

0

威望

379

eV
积分
441

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-1-24 20:38:49 | 只看该作者 Only view this author
sobereva 发表于 2025-1-24 14:31
恰当简化体系,只保留要用522博文里的做法算氢键键能附近的原子

1.卢老师,我想计算两条长链分子总的氢键结合能,两个分子间大概有10个氢键,我需要计算每个氢键的键能然后加和,这个计算量大吗,我没什么概念

2.简化体系是直接删除要计算氢键处以外的其他原子吗,大概像下图一样,保留多少个原子合适?

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

7#
发表于 Post on 2025-1-25 10:47:13 | 只看该作者 Only view this author
Sunca 发表于 2025-1-24 20:38
1.卢老师,我想计算两条长链分子总的氢键结合能,两个分子间大概有10个氢键,我需要计算每个氢键的键能然 ...

总氢键作用能直接用整体减两个片段能量算就完了,不需要用AIM拓扑分析的方式
北京科音自然科学研究中心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

62

帖子

0

威望

379

eV
积分
441

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2025-1-25 19:08:57 | 只看该作者 Only view this author
本帖最后由 Sunca 于 2025-1-25 22:48 编辑
sobereva 发表于 2025-1-25 10:47
总氢键作用能直接用整体减两个片段能量算就完了,不需要用AIM拓扑分析的方式

谢谢卢老师的回复,我还有下面3个问题想要问下您:
1.为什么整体减两个片段能量得到的是氢键相互作用,那非氢键的静电相互作用和范德华作用呢?因为我的目的就是单独计算出氢键作用
2.我的结构是从我做的分子动力学模拟中截取的一帧,计算单点能然后算氢键作用能的话,应该保持这样的结构不变还是需要再做优化?
3.我一共520个原子,用Gaussian,16个cpu核心算单点能用什么方法和基组算合适,比如B3LYP-D3(BJ)+什么基组好呢

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

9#
发表于 Post on 2025-1-26 13:55:37 | 只看该作者 Only view this author
Sunca 发表于 2025-1-25 19:08
谢谢卢老师的回复,我还有下面3个问题想要问下您:
1.为什么整体减两个片段能量得到的是氢键相互作用,[ ...

1 当前体系中氢键作用能占分子间总相互作用能的主要成分。如果你非要和范德华作用分得那么清楚,那就需要用Multiwfn做AIM拓扑分析。

2 应当优化

3 这种计算条件、这么大体系,用xtb里的GFN2-xTB或MOPAC里的PM6-D3H4X算相互作用能。DFT太费劲
北京科音自然科学研究中心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

62

帖子

0

威望

379

eV
积分
441

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2025-2-7 20:33:49 | 只看该作者 Only view this author
sobereva 发表于 2025-1-26 13:55
1 当前体系中氢键作用能占分子间总相互作用能的主要成分。如果你非要和范德华作用分得那么清楚,那就需要 ...

谢谢卢老师,这么大体系用计算BCP来单独算氢键键能的那个方法用不了是吧。那我这个体系有没有别的目前硬件条件能用的方法计算氢键键能?

6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

11#
发表于 Post on 2025-2-8 18:46:40 | 只看该作者 Only view this author
Sunca 发表于 2025-2-7 20:33
谢谢卢老师,这么大体系用计算BCP来单独算氢键键能的那个方法用不了是吧。那我这个体系有没有别的目前硬 ...

氢键位点及周围一圈原子用和我的JCC 2019文章相同的基组,距离氢键位点远的地方用小基组如3-21G,能算得动
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 18:03 , Processed in 0.197918 second(s), 23 queries , Gzip On.

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