第11届北京科音分子动力学与GROMACS培训班将于8月13~16日于北京举办,是系统性学习分子动力学与GROMACS使用的最佳机会,请点击此链接查看培训详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register
Views: 183|回复 Reply: 6

[GROMACS] 求助:为什么gro和top文件原子数对不上

[复制链接 Copy URL]

27

帖子

0

威望

51

eV
积分
78

Level 2 能力者

发表于 Post on 2022-7-4 21:33:39 | 显示全部楼层 Show all |阅读模式 Reading model
求助各位找了半天没发现原子数对不上啊
1111111.png

mix3.gro

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

topol.top

1.86 KB, 下载次数 Times of downloads: 7

20

帖子

2

威望

605

eV
积分
665

Level 4 (黑子)

发表于 Post on 2022-7-4 21:43:54 | 显示全部楼层 Show all
水的力场参数错了,tip4p是四点水模型,而你gro用的是三点水

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
lonemen + 5 赞!

查看全部评分 View all ratings

675

帖子

0

威望

2059

eV
积分
2734

Level 5 (御坂)

发表于 Post on 2022-7-4 21:44:40 | 显示全部楼层 Show all
看起来像是将单胞扩展以后,再加入了水和别的分子。你确认一下各个分子数量,以及gro里面是三点水模型还是四点水模型

27

帖子

0

威望

51

eV
积分
78

Level 2 能力者

 楼主 Author| 发表于 Post on 2022-7-5 11:49:29 | 显示全部楼层 Show all
casea 发表于 2022-7-4 21:43
水的力场参数错了,tip4p是四点水模型,而你gro用的是三点水

谢谢解决了,确实是四点水模型没改

27

帖子

0

威望

51

eV
积分
78

Level 2 能力者

 楼主 Author| 发表于 Post on 2022-7-5 11:51:35 | 显示全部楼层 Show all
本帖最后由 吧唧爱吃糖 于 2022-7-5 11:54 编辑
牧生 发表于 2022-7-4 21:44
看起来像是将单胞扩展以后,再加入了水和别的分子。你确认一下各个分子数量,以及gro里面是三点水模型还是 ...

改了四点水模型,但是grompp的时候提示这个
The optimal PME mesh load for parallel simulations is below 0.5
  and for highly parallel simulations between 0.25 and 0.33,
  for higher performance, increase the cut-off and the PME grid spacing.

我没管他继续mdrun之后提示


Energy minimization has stopped, but the forces have not converged to the
requested precision Fmax < 1 (which may not be possible for your system). It
stopped because the algorithm tried to make a new step whose size was too
small, or there was no change in the energy since last step. Either way, we
regard the minimization as converged to within the available machine
precision, given your starting configuration and EM parameters.
You might need to increase your constraint accuracy, or turn
off constraints altogether (set constraints = none in mdp file)

Steepest Descents converged to machine precision in 3152 steps,
but did not reach the requested Fmax < 1.
Potential Energy  = -8.80181406162373e+04
Maximum force     =  4.90022566027328e+02 on atom 1882
Norm of force     =  1.71190334880236e+01


请问一下师兄这个是怎么回事啊?以下是我的mdp文件和gro文件

min.mdp

532 Bytes, 下载次数 Times of downloads: 4

mix3_new.gro

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

27

帖子

0

威望

51

eV
积分
78

Level 2 能力者

 楼主 Author| 发表于 Post on 2022-7-5 15:41:02 | 显示全部楼层 Show all
忽略上面的内个note修改emtol为500后发现min.gro变成了这样,不知道哪里有问题
11111.png

27

帖子

0

威望

51

eV
积分
78

Level 2 能力者

 楼主 Author| 发表于 Post on 2022-7-5 17:09:02 | 显示全部楼层 Show all
吧唧爱吃糖 发表于 2022-7-5 15:41
忽略上面的内个note修改emtol为500后发现min.gro变成了这样,不知道哪里有问题

解决了,打扰各位了

本版积分规则 Credits rule

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

GMT+8, 2022-8-12 01:04 , Processed in 0.460890 second(s), 28 queries .

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