计算化学公社

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

[GROMACS] gromacs中获取原子坐标的命令

[复制链接 Copy URL]

72

帖子

0

威望

297

eV
积分
369

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
各位gromacs大神,
      大家好!有没有谁知道gromacs中如何获取原子的坐标的命令?我在手册上没有找到。
      谢谢大家!

294

帖子

0

威望

3533

eV
积分
3827

Level 5 (御坂)

2#
发表于 Post on 2017-9-26 20:12:49 | 只看该作者 Only view this author
gromacs中如何获取原子的坐标?问题不是很理解,你是需要在轨迹中提取坐标信息吗?
如果是的话:gmx trjconv 可以转换成 .pdb和.gro格式,直接使用文本编辑器打开

72

帖子

0

威望

297

eV
积分
369

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2017-9-27 11:01:51 | 只看该作者 Only view this author
ruanyang 发表于 2017-9-26 20:12
gromacs中如何获取原子的坐标?问题不是很理解,你是需要在轨迹中提取坐标信息吗?
如果是的话:gmx trjcon ...

谢谢您!我开始表达错了,我想问的是,如何得到一个分子质心的坐标。

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

4#
发表于 Post on 2017-9-27 12:00:36 | 只看该作者 Only view this author
用VMD最灵活省事
先用atomselect命令选择指定的区域,假设叫atomselect0,然后运行
measure center atomselect0 weight mass
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

72

帖子

0

威望

297

eV
积分
369

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2017-9-27 17:09:43 | 只看该作者 Only view this author
sobereva 发表于 2017-9-27 12:00
用VMD最灵活省事
先用atomselect命令选择指定的区域,假设叫atomselect0,然后运行
measure center atoms ...

gromacs不行么,因为我这个是gromacs的分析脚本,不想再用VMD了

115

帖子

0

威望

3892

eV
积分
4007

Level 6 (一方通行)

6#
发表于 Post on 2017-9-28 16:28:33 | 只看该作者 Only view this author
lujun 发表于 2017-9-27 17:09
gromacs不行么,因为我这个是gromacs的分析脚本,不想再用VMD了

seq 0 1999 | g_traj -f  AA/trajpbc.xtc -s AA/topol.tpr -oxt cg.gro -n mapping.ndx -com -ng 2000 -b 100000
可以用这种粗粒化映射的方法取出来
seq 0 1999 是取2000个
-com就表示计算的是质心
mapping里面一个质心放一组

评分 Rate

参与人数
Participants 2
eV +6 收起 理由
Reason
lonemen + 5 谢谢
sobereva + 1

查看全部评分 View all ratings

5

帖子

0

威望

138

eV
积分
143

Level 2 能力者

7#
发表于 Post on 2021-9-23 12:21:08 | 只看该作者 Only view this author
先定义好索引,然后 gmx traj -f -s -com -n -ox

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 02:28 , Processed in 0.146497 second(s), 21 queries , Gzip On.

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