“第25届北京科音初级量子化学培训班” 将于6月11至14日于北京举办,是初学者真正从头一次性完整、系统性学习量子化学计算,从而能很快上手量子化学研究的重要机会,比起自己摸索、鼓捣能少走无数弯路。报名正在进行中,请点击此链接查看,欢迎参加!

计算化学公社

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

[GROMACS] 研究小分子跨膜,Gromacs加墙体(WALL)模拟报错求助

[复制链接 Copy URL]

7

帖子

0

威望

95

eV
积分
102

Level 2 能力者

我研究的内容是小分子跨膜,所以我想利用pbc=xy来设置在Z方向上不周期防止小分子通过边界跨越膜,通过附件中nvt.mdp的设置,在Z方向上添加了碳原子墙(C来源于小分子的atomtype)实现添加墙体,但是不管本地运行还是服务器运行都没能跑通,本地运行的log文件与服务器运行的log文件也在附件内,Gromacs所用版本为2020.6,本地运行显示如图,nvt模拟确实运行了一段时间,跑了一定的步数,但是突然就停止了,没有gromacs运行结束的结束语也没有报错,log里也没有相关信息,只有服务器日志有报错显示,但是不知道如何解决,请各位大佬帮忙看看,是不是我哪里配置错了,top、mdp和log文件已上传,但是gro的构型文件太大不够传

c95ab7e0a09f2eae6a55d1f82815f54.png (36.14 KB, 下载次数 Times of downloads: 12)

本地运行结束图

本地运行结束图

4d7511feac30ca12300640e2dd44303.png (23.48 KB, 下载次数 Times of downloads: 6)

服务器运行结束图

服务器运行结束图

nvt.mdp

2.34 KB, 下载次数 Times of downloads: 4

配置文件

nvt.log

27.47 KB, 下载次数 Times of downloads: 2

本地运行日志

topol.top

772 Bytes, 下载次数 Times of downloads: 2

top文件

nvt-服务器.log

22.7 KB, 下载次数 Times of downloads: 1

服务器运行日志

382

帖子

0

威望

950

eV
积分
1332

Level 4 (黑子)

2#
发表于 Post on 3 day ago | 只看该作者 Only view this author
去掉墙,看能不能正常跑

6万

帖子

99

威望

6万

eV
积分
127221

管理员

公社社长

3#
发表于 Post on 前天 01:29 | 只看该作者 Only view this author
如置顶的新社员必读贴和论坛首页的公告栏所示,上传超过500KB的gro等文本型文件前必须先压缩再上传,以节约论坛空间、节约下载耗时。别不压缩就试图直接上传,论坛规则故意让这种行为无法成功。更何况还有叫网盘的东西

按照http://sobereva.com/soft/Sobtop#FAQ8的逻辑排查。显然遇到问题应当先简化体系,诸如想到先去掉wall尝试是否能跑成功判断和wall是否有关。以及增加轨迹输出频率后,检查崩溃后轨迹的最后一帧的结构看是否能看出什么
北京科音自然科学研究中心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

7

帖子

0

威望

95

eV
积分
102

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 前天 08:19 | 只看该作者 Only view this author
KazusaT 发表于 2026-5-29 22:45
去掉墙,看能不能正常跑

去掉WALL,设置PBC=xyz可以跑通,我确认问题和wall有关

1618

帖子

0

威望

5246

eV
积分
6864

Level 6 (一方通行)

5#
发表于 Post on 前天 08:58 | 只看该作者 Only view this author
Rororozoroni 发表于 2026-5-30 08:19
去掉WALL,设置PBC=xyz可以跑通,我确认问题和wall有关

你把Z轴再延长一点,让wall距离盒子内物质远一些。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
SharkYYX2025 + 2 正解

查看全部评分 View all ratings

又菜又爱玩

7

帖子

0

威望

95

eV
积分
102

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 前天 13:14 | 只看该作者 Only view this author
牧生 发表于 2026-5-30 08:58
你把Z轴再延长一点,让wall距离盒子内物质远一些。

好的,我尝试一下

175

帖子

0

威望

981

eV
积分
1156

Level 4 (黑子)

开心小鲨鱼 ORCA的饲养员/夜宵

7#
发表于 Post on 前天 13:34 | 只看该作者 Only view this author
楼上说得对,很可能是墙重叠了实际原子。
我也遇到过这种问题。我是用XY周期性,Z方向依次是:Z=0 墙,晶格模型,吸附的小分子,真空区10埃,Z=max 墙。晶格和最下面的墙需要隔开1埃,packmol建模时候坐标就要设定留点距离(晶格不能太接近盒子底部)。因为墙是mdp top后加的,gro里没有实体原子坐标,很容易晶格建模时候忘记留点空隙。第二,真空区也要足量。

7

帖子

0

威望

95

eV
积分
102

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 5 hour ago | 只看该作者 Only view this author
牧生 发表于 2026-5-30 08:58
你把Z轴再延长一点,让wall距离盒子内物质远一些。

我延长之后模拟的最后一帧是这样的,不太对劲,水分子跑出盒子了,然后水分子还进入到了膜内部,正常模拟是不可能进入到内部的,相比正常不加墙的模拟就是这两处区别

733a35df2e5fe126b7a620f28b4d32c.png (411.99 KB, 下载次数 Times of downloads: 0)

NVT.trr的最后一帧

NVT.trr的最后一帧

本版积分规则 Credits rule

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

GMT+8, 2026-6-1 15:57 , Processed in 0.491974 second(s), 24 queries , Gzip On.

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