计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 怎么从gibbs.xpm文件中提取自由能数据？

怎么从gibbs.xpm文件中提取自由能数据？

查看数: 4167 | 评论数: 4 | 收藏 Add to favorites 2

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

帖子模式

bingzan

发布时间: 2022-8-17 08:52

正文摘要:

各位老师好，我需要绘制自由能形貌图（横轴是RMSD，纵轴为Rg），不是主成分分析法。算出rmsd和rg后，我用gmx sham输出了gibbs.xpm。我需要提取数据自行绘制。所以我需要把该xpm文件 ...

回复 Reply

starparyer 发表于 Post on 2025-6-21 13:58:05

bingzan 发表于 2022-8-17 17:00
成功了！我发现，如果使用gmx sham -tsham 310 -nlevels 100 -f output.xvg -ls gibbs.xpm -g gibbs.log -l ...

非常感谢楼主提示，找到问题了，链接那个有 -tsham 310 -nlevels 100 -f 特殊设置，所以不能直接copy，可以修改一下那个xpm_to_txt代码，修改之后可以针对双字符xpm进行读取，最后得到结果，感激~

bingzan 发表于 Post on 2022-8-17 17:00:23

成功了！我发现，如果使用gmx sham -tsham 310 -nlevels 100 -f output.xvg -ls gibbs.xpm -g gibbs.log -lsh enthalpy.xpm -lss entropy.xpm 的话，输出的gibbs.xpm是双字符表示颜色的。但是如果按照链接中文献，只用gmx sham -f gsham-file.xvg -ls free-energy-landscape.xpm的话，输出的Gibbs.xpm是单字符表示颜色的，就能用xpm2txt代码了。

bingzan 发表于 Post on 2022-8-17 16:27:11

本帖最后由 bingzan 于 2022-8-17 16:30 编辑

谢谢楼上的解答，但是我不会python，目前刚开始学，不知道怎么改。请问还有别的方式可以提取数据吗？或许gmx energy可以吗？

lyj714 发表于 Post on 2022-8-17 15:28:48

只能自己修改一下脚本，此脚本没有考虑这种双字符表示颜色的xpm文件

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