计算化学公社

 找回密码 Forget password
 注册 Register

gmx dipoles 计算氨基酸偶极矩时出错

查看数: 2610 | 评论数: 7 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-9-16 21:47

正文摘要:

    各位老师好,今天利用gmx dipoles 计算某个氨基酸偶极矩时出现了图中所示的错误。好像是顺序对应问题,但是我检查过了index与gro文件,该氨基酸的原子序号与index文件完全对应的。想问下老师还有没有 ...

回复 Reply

新月之痕 发表于 Post on 2023-3-19 09:45:00
uenh1998 发表于 2023-3-16 17:03
在index文件里添加要计算的氨基酸的序号  dipoles计算时候带上相应的index文件

非常感谢您的解答,我就是用的这个方法,结果会报同样的错误。
后来我改用社长的vmd方法计算这些氨基酸的偶极矩了,和gmx的误差在0.5%以内,很可靠
uenh1998 发表于 Post on 2023-3-16 17:03:00
新月之痕 发表于 2023-3-16 12:03
是的,但是我只想要计算蛋白质分子中的部分氨基酸的偶极矩应该怎么做呢?请问你当时是怎么解决的呢?

在index文件里添加要计算的氨基酸的序号  dipoles计算时候带上相应的index文件
新月之痕 发表于 Post on 2023-3-16 12:03:33
uenh1998 发表于 2023-3-16 11:00
应该是用gromacs计算偶极矩,定义的组必须是完整的分子,否则会出错

是的,但是我只想要计算蛋白质分子中的部分氨基酸的偶极矩应该怎么做呢?请问你当时是怎么解决的呢?
uenh1998 发表于 Post on 2023-3-16 11:00:30
新月之痕 发表于 2023-3-16 10:25
你好,我也遇到了同样的问题,请问你是如何解决的呢?

应该是用gromacs计算偶极矩,定义的组必须是完整的分子,否则会出错
新月之痕 发表于 Post on 2023-3-16 10:25:20
uenh1998 发表于 2022-9-17 14:57
感谢老师,已经尝试成功了!

你好,我也遇到了同样的问题,请问你是如何解决的呢?
uenh1998 发表于 Post on 2022-9-17 14:57:08
sobereva 发表于 2022-9-17 00:19
没细节不好说
也可以按下文把原子电荷载入VMD后,用自带的Extensions - Visualization - Dipole Moment Wa ...

感谢老师,已经尝试成功了!
sobereva 发表于 Post on 2022-9-17 00:19:02
没细节不好说
也可以按下文把原子电荷载入VMD后,用自带的Extensions - Visualization - Dipole Moment Watcher计算指定选区的偶极矩,还能可视化
将GROMACS的原子电荷信息读入VMD的方法
http://sobereva.com/365http://bbs.keinsci.com/thread-5417-1-1.html

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
Jus + 5 谢谢

查看全部评分 View all ratings

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

GMT+8, 2026-2-25 09:34 , Processed in 0.172839 second(s), 26 queries , Gzip On.

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