计算化学公社

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

[GROMACS] 求助GROMACS进行NVT模拟报错,结构会断裂

[复制链接 Copy URL]

2

帖子

0

威望

45

eV
积分
47

Level 2 能力者

本帖最后由 hezhichao 于 2023-9-1 09:10 编辑

各位老师,使用GROMACS做分子动力学模拟时,做好能量最小化后,再进行NVT模拟会报错,结构会发生断裂,如图所示,求助各位老师该如何解决该问题,谢谢

我的体系中,基底是石墨烯,跑的分子结构如下图所示,断裂部分是共价键,模型是用高斯建的,所用力场是amber力场,用sob的那个软件直接生成的。

微信图片_20230831090305.jpg (136.63 KB, 下载次数 Times of downloads: 11)

微信图片_20230831090305.jpg

1693529820979.jpg (69.95 KB, 下载次数 Times of downloads: 10)

1693529820979.jpg

160

帖子

7

威望

786

eV
积分
1086

Level 4 (黑子)

2#
发表于 Post on 2023-8-31 15:10:04 | 只看该作者 Only view this author
又见提问一图流。。。细节呢?

石墨烯(疑似)上面具体是啥分子?
模型怎么建的?
断是怎么个断法?
非共价键断了那不是正常?
如果是共价键断了,那检查力场有没有错?

建议参考此文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://bbs.keinsci.com/forum.php ... 87&fromuid=6031
(出处: 计算化学公社)

160

帖子

7

威望

786

eV
积分
1086

Level 4 (黑子)

3#
发表于 Post on 2023-9-11 11:06:56 | 只看该作者 Only view this author
帖子编辑后大概懂了, 未必是分子断了只是你没显示出来,你用vmd显示出xy两个方向看看就知道了。
vmd控制台里pbc unwrap显示就正常了

2

帖子

0

威望

45

eV
积分
47

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2023-9-11 15:00:36 | 只看该作者 Only view this author
本帖最后由 hezhichao 于 2023-9-12 15:17 编辑
tomwong4253 发表于 2023-9-11 11:06
帖子编辑后大概懂了, 未必是分子断了只是你没显示出来,你用vmd显示出xy两个方向看看就知道了。
vmd控制 ...

感谢您的答复,但是不止是显示分子断裂的问题,只要是会报错,“fatal error”,无法进行继续计算。

160

帖子

7

威望

786

eV
积分
1086

Level 4 (黑子)

5#
发表于 Post on 2023-9-16 13:06:31 | 只看该作者 Only view this author
hezhichao 发表于 2023-9-11 15:00
感谢您的答复,但是不止是显示分子断裂的问题,只要是会报错,“fatal error”,无法进行继续计算。

那就是你虚线部分(盲猜是C=S和N-H之间)的参数没有给好,所以报错了,仔细检查你的itp文件

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 20:34 , Processed in 0.295508 second(s), 30 queries , Gzip On.

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