计算化学公社

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

[GROMACS] packmol填充分子,三明治界面模型,gromacs跑动力学崩溃

[复制链接 Copy URL]

8

帖子

0

威望

83

eV
积分
91

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
对表面活性剂和水溶液的三明治模型进行分子动力学模拟
用Packmol 进行装填分子,z方向一共三层,使表面活性剂分子头基朝水
上中两层或者中下两层合在一起都能正常跑动力学
合在一起就崩溃
调整了水和表活的个数,还有表活和水层的距离,都没有效果
请问怎么解决,谢谢!

这是gromacs报错信息
There are inconsistent shifts over periodic boundaries in a molecule type consisting of 3 atoms.      
The longest distance involved in such interactions is 5.411 nm which is above half the box length.
Either you have excessively large distances between atoms in bonded interactions or your system is exploding.


There is no domain decomposition for 36 ranks that is compatible with the given box and a minimum cell size of 15.9594 nm
Change the number of ranks or mdrun option -rdd or -dds
Look in the log file for details on the domain decomposition


这是packmol 的输入文件
structure bs-12.pdb
  number 1
  atoms 38
  inside box 0. 0. 32. 67. 67. 35
end atoms
  atoms 48
  inside box 0. 0. 17. 67. 67. 18
end atoms
end structure

structure water.pdb
  number 7000
  inside box 0. 0. 35. 69. 69. 105
end structure

structure Na.pdb
  number 120
  inside box 0. 0. 35. 69. 69. 105
end structure

structure Cl.pdb
  number 120
  inside box 0. 0. 35. 69. 69. 105
end structure

structure bs-12.pdb
  number 1
  atoms 38
  inside box 0. 0. 105. 67. 67. 108
end atoms
  atoms 48
  inside box 0. 0. 122. 67. 67. 123
end atoms
end structure

1560

帖子

0

威望

4997

eV
积分
6557

Level 6 (一方通行)

3#
发表于 Post on 2024-12-24 13:34:51 | 只看该作者 Only view this author
大大大概率是某个分子的拓扑参数不合理
又菜又爱玩

4

帖子

0

威望

163

eV
积分
167

Level 3 能力者

2#
发表于 Post on 2024-12-24 11:21:53 | 只看该作者 Only view this author
请问解决了吗?也遇见类似问题了

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

GMT+8, 2026-2-22 05:53 , Processed in 0.156584 second(s), 21 queries , Gzip On.

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