计算化学公社

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

[GROMACS] 求助gromacs在有机溶剂体系的平衡相模拟计算过程中出现吐核错误

[复制链接 Copy URL]

9

帖子

0

威望

61

eV
积分
70

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师好,我在不同有机溶剂体系中模拟沥青质的聚集行为,但在平衡相模拟过程中断,终端中出现“段错误(吐核)”的字样,具体警告内容如下所示。想请教一下为什么会出现“吐核”错误,该错误出现后应该从哪个文件入手进行修改呢?谢谢。
WARNING: Using the slow plain C kernels. This should not happen during routine usage on supported platforms.starting mdrun 'asp-box2'
500000 steps,   1000.0 ps.
step 2700, will finish Wed Mar 30 12:13:53 2022
Step 2761, time 5.522 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.000011, max 0.000107 (between atoms 2537 and 2573)
bonds that rotated more than 30 degrees:
atom 1 atom 2  angle  previous, current, constraint length
   2189   2237   30.4    0.1090   0.1090      0.1090
   2190   2238   31.0    0.1090   0.1090      0.1090

Step 2769, time 5.538 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.270460, max 37.492767 (between atoms 2190 and 2238)
bonds that rotated more than 30 degrees:
atom 1 atom 2  angle  previous, current, constraint length
   2189   2237   90.0    0.1090   3.6836      0.1090
   2190   2238   90.0    0.1090   4.1957      0.1090
Wrote pdb files with previous and current coordinates

WARNING: Listed nonbonded interaction between particles 2184 and 2238
at distance 4.286 which is larger than the table limit 2.009 nm.

This is likely either a 1,4 interaction, or a listed interaction inside
a smaller molecule you are decoupling during a free energy calculation.
Since interactions at distances beyond the table cannot be computed,
they are skipped until they are inside the table limit again. You will
only see this message once, even if it occurs for several interactions.

IMPORTANT: This should not happen in a stable simulation, so there is probably something wrong with your system. Only change the table-extension distance in the mdp file if you are really sure that is the reason.



Step 2770, time 5.54 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.000658, max 0.119843 (between atoms 2189 and 2237)
bonds that rotated more than 30 degrees:
atom 1 atom 2  angle  previous, current, constraint length
   2189   2237   90.0    3.6836   0.1221      0.1090
   2190   2238   90.0    4.1957   0.1110      0.1090

Step 2771, time 5.542 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.244306, max 33.750931 (between atoms 2190 and 2238)
bonds that rotated more than 30 degrees:
atom 1 atom 2  angle  previous, current, constraint length
   2189   2237   90.0    0.1221   3.3524      0.1090
   2190   2238   90.0    0.1110   3.7879      0.1090
Wrote pdb files with previous and current coordinates

Step 2772, time 5.544 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.000011, max 0.000082 (between atoms 4518 and 4571)
bonds that rotated more than 30 degrees:
atom 1 atom 2  angle  previous, current, constraint length
   2185   2235   62.9    0.1094   0.1090      0.1090
   2188   2236   58.0    0.1094   0.1090      0.1090
   2190   2238   32.7    3.8030   0.1090      0.1090

Step 2773, time 5.546 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.103036, max 18.976585 (between atoms 2190 and 2238)
bonds that rotated more than 30 degrees:
atom 1 atom 2  angle  previous, current, constraint length
   2185   2235   55.2    0.1090   0.1090      0.1090
   2188   2236   49.7    0.1090   0.1090      0.1090
   2189   2237   70.7    0.1090   0.1090      0.1090
   2190   2238   90.0    0.1090   2.1774      0.1090
Wrote pdb files with previous and current coordinates

Step 2774, time 5.548 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 13846.613281, max 2526945.750000 (between atoms 2192 and 2239)
bonds that rotated more than 30 degrees:
atom 1 atom 2  angle  previous, current, constraint length
   2185   2235   90.0    0.1090   0.6697      0.1090
   2189   2237   90.1    0.1090 17540.2715      0.1090
   2190   2238   89.7    2.1774 15725.8203      0.1090
   2192   2239   90.0    0.1080 272910.2500      0.1080
   2208   2244   90.0    0.1080 28832.7930      0.1080
Wrote pdb files with previous and current coordinates
段错误(吐核)


1560

帖子

0

威望

5001

eV
积分
6561

Level 6 (一方通行)

2#
发表于 Post on 2022-3-30 07:46:00 | 只看该作者 Only view this author
极大概率是因为top文件不合理
又菜又爱玩

9

帖子

0

威望

61

eV
积分
70

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2022-3-30 13:23:17 | 只看该作者 Only view this author
牧生 发表于 2022-3-30 07:46
极大概率是因为top文件不合理

你好,我不太理解,想请问一下可能存在的不合理的因素有哪些呢

1560

帖子

0

威望

5001

eV
积分
6561

Level 6 (一方通行)

4#
发表于 Post on 2022-3-30 16:07:24 | 只看该作者 Only view this author
Snown 发表于 2022-3-30 13:23
你好,我不太理解,想请问一下可能存在的不合理的因素有哪些呢

http://bbs.keinsci.com/thread-25599-1-1.html

第四楼
又菜又爱玩

9

帖子

0

威望

61

eV
积分
70

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2022-3-30 18:36:20 | 只看该作者 Only view this author
牧生 发表于 2022-3-30 16:07
http://bbs.keinsci.com/thread-25599-1-1.html

第四楼

多谢指导

6万

帖子

99

威望

6万

eV
积分
125153

管理员

公社社长

6#
发表于 Post on 2022-3-31 15:01:00 | 只看该作者 Only view this author
北京科音自然科学研究中心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, 2026-2-23 22:42 , Processed in 0.152932 second(s), 20 queries , Gzip On.

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