计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 使用tcl脚本计算氢键数结果有误

使用tcl脚本计算氢键数结果有误

查看数: 641 | 评论数: 4 | 收藏 Add to favorites 0

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

帖子模式

zhangfaxue

发布时间: 2024-12-7 13:58

正文摘要:

脚本如下，想要计算两个组分氢键数目随着时间的变化，之前测试了表面活性剂和水分子之间氢键数目发现计算出来结果全是3，为了验证脚本是否存在问题就修改了两个组分都为水分子，结果计算还是氢键数目全部是3，请教各 ...

回复 Reply

zhangfaxue 发表于 Post on 2024-12-8 09:49:31

fhh2626 发表于 2024-12-7 20:39
VMD的extension里面有统计氢键的插件

好的好的谢谢

zhangfaxue 发表于 Post on 2024-12-8 09:49:16

Uus/pMeC6H4-/キ发表于 2024-12-7 14:17
measure hbonds返回的是由给体原子index、受体原子index、氢原子index共3个元素数量相同的子列表构成的列表 ...

谢谢回复！问题已解决。回复是最基本的，很抱歉上一条我疏忽忘记了没有反馈。非常感谢您的帮助。

fhh2626 发表于 Post on 2024-12-7 20:39:51

VMD的extension里面有统计氢键的插件

Uus/pMeC6H4-/キ 发表于 Post on 2024-12-7 14:17:10

measure hbonds返回的是由给体原子index、受体原子index、氢原子index共3个元素数量相同的子列表构成的列表，而llength直接返回的是列表中元素数量而不关心元素是不是还有列表，所以直接把measure hbonds返回结果存到变量$temp里再用llength $temp可不就取成3了。统计氢键数量应该用其中一个子列表来数，比如用第一个子列表的命令写作[llength [lindex $temp 0]]。

有个小小的请求，如果像上帖这样提问得到回答，可否提供个反馈呢？回个帖或者评个分都行，无论回答是很有帮助还是完全不对，提问者的回应可以给回答者一个确认，也可以帮助有相同问题的后来者解惑……

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