计算化学公社

 找回密码 Forget password
 注册 Register
Views: 13352|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 求问:如何使用gmx自带xpm2ps工具来作图分析xpm文件内容

[复制链接 Copy URL]

129

帖子

0

威望

955

eV
积分
1084

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
1. 使用命令gmx hbond -s md.tpr -f md.xtc -hbm得到的hbmap.xpm文件,文件内容如下(加粗字体):
/* XPM */
/* This file can be converted to EPS by the GROMACS program xpm2ps */
/* title:   "Hydrogen Bond Existence Map" */
/* legend:  "Hydrogen Bonds" */
/* x-label: "Time (ps)" */
/* y-label: "Hydrogen Bond Index" */
/* type:    "Discrete" */
static char *gromacs_xpm[] = {
"21 13   2 1",
"   c #FFFFFF " /* "None" */,
"o  c #FF0000 " /* "Present" */,
/* x-axis:  290000 290500 291000 291500 292000 292500 293000 293500 294000 294500 295000 295500 296000 296500 297000 297500 298000 298500 299000 299500 300000 */
/* y-axis:  0 1 2 3 4 5 6 7 8 9 10 11 12 */
"oo o        ooooo oo ",
"  o  o o            o",
"                 o   ",
"   o                 ",
"     o   oo oo       ",
"                  o  ",
"  o o                ",
"o o    o             ",
"            o    o   ",
"                 o   ",
"           o o       ",
"          o          ",
"oooo o     ooo       "


这文件里的圈圈应该是表示氢键,我想问下为什么会有13列的圆圈而且每一列上圆圈的位置也不一样?这代表的是什么意思。

2. 命令gmx mdmat 来寻找蛋白质内两两氨基酸距离最小值,生成的dm.xpm矩阵文件手册里说是可通过工具xpm2ps可以转换成Postscript picture进而可以作图来查看,但是有以下出错消息,想知道为什么没法用xpm2ps来分析?


Command line:
  gmx xpm2ps -f dm.xpm

There is 1 matrix in dm.xpm
Matrix 0 is 348 x 348

-------------------------------------------------------
Program:     gmx xpm2ps, version 2019

Standard library logic error (bug):
(exception type: St11logic_error)
basic_string::_S_construct null not valid

5

帖子

0

威望

75

eV
积分
80

Level 2 能力者

2#
发表于 Post on 2021-9-30 16:38:54 | 只看该作者 Only view this author
请问第二个问题怎么解决呀?我出现了一样的问题

129

帖子

0

威望

955

eV
积分
1084

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-2-17 08:13:43 | 只看该作者 Only view this author
程序正常安装使用一般不会有问题

5

帖子

0

威望

75

eV
积分
80

Level 2 能力者

4#
发表于 Post on 2022-3-28 20:35:32 | 只看该作者 Only view this author
why99 发表于 2021-9-30 16:38
请问第二个问题怎么解决呀?我出现了一样的问题

使用18版可以解决第二个问题

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 01:20 , Processed in 0.149833 second(s), 20 queries , Gzip On.

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