计算化学公社

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

[GROMACS] 进行pull处理时提示distance与box size不合理报错

[复制链接 Copy URL]

97

帖子

0

威望

634

eV
积分
731

Level 4 (黑子)

进行拉伸过程中,提示如下错误:

然后尝试将盒子z方向扩大到20nm,在进行pull,错误一样:

查看前面的warning发现:

但是group1 2 相隔距离可视化所示并没有那么远(盒子尺寸5 5 20):

mdp设置如下:

请问:
1、pull时计算两个分子距离是怎么计算的,提示的6.525nm怎么获得的?质心距离吗?(我的大分子大概最长边2.8nm)
2、出现上述错误可能是哪些原因?该怎么处理?

29

帖子

0

威望

854

eV
积分
883

Level 4 (黑子)

2#
发表于 Post on 2025-1-18 08:48:10 | 只看该作者 Only view this author
本帖最后由 FrancisCho 于 2025-1-18 08:59 编辑

1、6.525nm应该是您当前模拟状态下两个组的质心的距离(我翻看了一下我自己模拟中的nohup.out中初始distance是0.06nm,和用gmx distance给出的值很接近)
2、一般GMX中SMD的拖动长度不应超过该方向上盒子尺寸的一半,如要超过,则应使用pull_coord1_geometry    = direction-periodic、pull_coord1_vec         = 0.0 0.0 -1.0以及pull-pbc-ref-prev-step-com        = yes,在这种情况下不能进行NPT,但不会报错(我做US就是这样设置的)

Snipaste_2025-01-18_08-46-27.png (167.16 KB, 下载次数 Times of downloads: 1)

GMX官网对于pull_coord1_geometry = direction-periodic描述

GMX官网对于pull_coord1_geometry    = direction-periodic描述

Snipaste_2025-01-18_08-58-30.png (94.04 KB, 下载次数 Times of downloads: 1)

gmx distance

gmx distance

Snipaste_2025-01-18_08-57-04.png (111.57 KB, 下载次数 Times of downloads: 0)

pull-code

pull-code

97

帖子

0

威望

634

eV
积分
731

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-1-20 18:30:23 | 只看该作者 Only view this author
FrancisCho 发表于 2025-1-18 08:48
1、6.525nm应该是您当前模拟状态下两个组的质心的距离(我翻看了一下我自己模拟中的nohup.out中初始distanc ...

谢谢指导,确实是质心的距离,我再检查一下

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 07:22 , Processed in 0.148373 second(s), 23 queries , Gzip On.

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