计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First Principle) › AIMD模拟后借助gmx hbond计算水溶剂分子之间氢键数目问 ...

AIMD模拟后借助gmx hbond计算水溶剂分子之间氢键数目问题

查看数: 636 | 评论数: 7 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

uyhvtv

发布时间: 2023-3-21 19:37

正文摘要:

想请问大家一个关于aimd结果用gromacs计算氢键的问题，本人先将O-H做成group后，并将其作为供体和受体，用gmx_mpi hbond计算氢键数目，出现了donor为0的情况，导致没有结果输出，是哪里出现了问题？

回复 Reply

uyhvtv 发表于 Post on 2023-4-11 20:17:37

lyj714 发表于 2023-3-25 00:12
就是top中必须存在[ bonds] 这个字段，并且指定所有的水的氧原子和对应H的成键才行，因为只有成键的氧原 ...

感谢您的回复，我在运行H2O的AIMD过程中，H-O的键长和键角可能是改变的，那[ bonds] 是如何才能添加呢，谢谢

uyhvtv 发表于 Post on 2023-4-11 20:16:18

感谢您的回复，我在运行H2O的AIMD过程中，H-O的键长和键角可能是改变的，那[ bonds] 是如何才能添加呢，谢谢

lyj714 发表于 Post on 2023-3-25 00:12:24

uyhvtv 发表于 2023-3-24 21:49
再次感谢您的回复我查看文件top文献，H和O也是有的，上述是没有截图到，这里关于“成键（O-H之间的)也要 ...

就是top中必须存在[ bonds] 这个字段，并且指定所有的水的氧原子和对应H的成键才行，因为只有成键的氧原子，才能被gmx判断是供体氧，否则孤立的氧只能被认为是受体

lyj714 发表于 Post on 2023-3-22 16:19:59

本帖最后由 lyj714 于 2023-3-22 16:22 编辑

从你发的一系列截图不难看出拓扑中的原子名称都是Ru，这肯定是不行的。gmx hbond判供体和受体只根据拓扑中原子名称的第一个字符和成键关系，而非轨迹gro中的。虽然你用其他工具生成了伪top，但是你必须确保名称正确，并且成键（O-H之间的)也要正确。另外，你也可以用vmd进行氢键分析，可以直接通过gro轨迹中的原子名称

sobereva 发表于 Post on 2023-3-21 23:04:17

gro文件里原子名必须合适，能让GROMACS认得出元素是氧和氢

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register