找回密码
 立即注册
查看: 19|回复: 4

求助:gromacs模拟MXene膜固体材料

[复制链接]

1

主题

2

回帖

17

积分

新手上路

积分
17
发表于 2026-3-26 15:49:31 | 显示全部楼层 |阅读模式
我做的模型大致如下(海水淡化):石墨烯| 海水区 | MXene | 纯水区 | 石墨烯,海水区模拟的是0.5M的海水(2770H2O+23Na+和23Cl-),纯水区添加了717个H2O,MXene材料用的是Ti3C2,在组装盒子把pdb文件转为gro文件之后(第一张图),做完em我用VMD查看了em.gro文件,出现了如2图所示(其中比较大的蓝色的球代表MXene材料里的C,粉色的代表Ti原子,石墨烯里的C原子也是蓝色的),可以看到固体材料原子已经飞散了,还未到后续施加压力步骤盐离子已经提前扩散,求大佬指教!
屏幕截图 2026-03-26 154552.png 屏幕截图 2026-03-26 154359.png

2

主题

56

回帖

155

积分

版主

积分
155
发表于 2026-3-27 08:57:58 | 显示全部楼层
石墨烯应该冻结住

1

主题

2

回帖

17

积分

新手上路

积分
17
 楼主| 发表于 3 天前 | 显示全部楼层
lauthirteen 发表于 2026-3-27 08:57
石墨烯应该冻结住

是通过写索引文件来冻结吗,我试着在em的时候写入了冻结的参数,可是还是会原子飞散,求指教!

2

主题

56

回帖

155

积分

版主

积分
155
发表于 3 天前 | 显示全部楼层
索引文件不是冻结文件,可以在mdp里面加入:
freezegrps = MOL
freezedim = Y Y Y
这种关键字。MOL改为你想冻结的残基名称。如果有多个,可以继续追加,类似:
freezegrps = MOL  UNK
freezedim = Y Y Y Y Y Y

1

主题

2

回帖

17

积分

新手上路

积分
17
 楼主| 发表于 前天 08:25 | 显示全部楼层
lauthirteen 发表于 2026-4-1 16:47
索引文件不是冻结文件,可以在mdp里面加入:
freezegrps = MOL
freezedim = Y Y Y

非常感谢!我这就去试试
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|模拟之家-SimuHome ( 蜀ICP备2023003244号-3 )

GMT+8, 2026-4-4 11:25 , Processed in 0.028408 second(s), 24 queries .

快速回复 返回顶部 返回列表