计算化学公社

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

[GROMACS] gmx trjconv命令是否可以使用-select语句?

[复制链接 Copy URL]

262

帖子

0

威望

635

eV
积分
897

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
最近在学习伞形采样,需要根据反应坐标从SMD模拟的xtc文件提取一系列结构进行后续的WHAM模拟和计算,因为距离并不是随时间均匀变化的,所以写了一个python小脚本按0.1 nm的窗口来提取,其中使用了如下命令
  1. gmx trjconv -s pull.tpr -f pull.xtc -o conf{i}.gro -dump {sampling_time:.1f} -quiet -select "group 0"
复制代码
但是发现trjconv好像不能用select语句,这里有什么其他办法吗?


849

帖子

4

威望

1920

eV
积分
2849

Level 5 (御坂)

A Student

2#
发表于 Post on 2025-10-10 08:26:09 | 只看该作者 Only view this author
本帖最后由 student0618 于 2025-10-10 09:27 编辑

The following script use e.g. gmx distance to get distance of each frame then extract https://github.com/Lemkul-Lab/gm ... lin/find_windows.py

In short, use echo 0 | gmx trjconv [options]
敬仰一针见血的指责,厌倦别有用心的赞美。

262

帖子

0

威望

635

eV
积分
897

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-10-10 15:13:38 | 只看该作者 Only view this author
student0618 发表于 2025-10-10 08:26
The following script use e.g. gmx distance to get distance of each frame then extract https://github ...

Thanks for the suggestion. I'll give this method a try.

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 03:41 , Processed in 0.182756 second(s), 20 queries , Gzip On.

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