“第25届北京科音初级量子化学培训班” 将于6月11至14日于北京举办,是初学者真正从头一次性完整、系统性学习量子化学计算,从而能很快上手量子化学研究的重要机会,比起自己摸索、鼓捣能少走无数弯路。报名正在进行中,请点击此链接查看,欢迎参加!

计算化学公社

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

[Amber] Amber中CPPTRAJ中氢键分析相关问题求助

[复制链接 Copy URL]

102

帖子

0

威望

632

eV
积分
734

Level 4 (黑子)

蛋白设计爱好者

跳转到指定楼层 Go to specific reply
楼主
https://amberhub.chpc.utah.edu/h ... d-a-small-molecule/

以上是Amber官方氢键分析教程

给出的例子中

parm vph.prmtop
trajin vph.nc 1 last 1
hbond contacts :1-154 avgout avg.dat series uuseries hbond.gnu nointramol
go
lifetime contacts[solutehb] out contacts-lifetime.dat
go


请问各位,命令行中载入轨迹的“1 last 1”是什么意义呢?


310

帖子

0

威望

1525

eV
积分
1835

Level 5 (御坂)

2#
发表于 Post on 2022-4-17 01:13:19 | 只看该作者 Only view this author
本帖最后由 Frozen-Penguin 于 2022-4-17 11:16 编辑

读取范围是轨迹的第一帧到最后一帧,每隔一帧读取一次,上述写法就是读取整个轨迹。
如果写 1 100 10,就会读取轨迹中的10帧,如果什么都不写则是读取完整的轨迹。
类似gromacs中的-b-e-dt的用法,不过amber用的是帧数,gromacs用的是时间。

102

帖子

0

威望

632

eV
积分
734

Level 4 (黑子)

蛋白设计爱好者

3#
 楼主 Author| 发表于 Post on 2022-4-17 12:28:48 | 只看该作者 Only view this author
Frozen-Penguin 发表于 2022-4-17 01:13
读取范围是轨迹的第一帧到最后一帧,每隔一帧读取一次,上述写法就是读取整个轨迹。
如果写 1 100 10,就 ...

感谢,明白了

本版积分规则 Credits rule

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

GMT+8, 2026-6-3 06:02 , Processed in 0.431923 second(s), 25 queries , Gzip On.

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