24小时热门版块排行榜    

查看: 2053  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yiha

新虫 (初入文坛)

[求助] dump 输出原子名称

各位好:
      lammps定义原子时用数字1,2,3,导致在vmd中显示的name也为1,2,3,请问怎么才能使lammps结果在vmd显示的name为真实名字名称N,H,O?

8660473.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
chaizhm: 金币+5, 谢谢~ 2012-10-25 09:23:05
现成的没有,不过思路简单:
1、把lammps轨迹载入VMD
2、做循环,直到你的轨迹全部循环完毕
3、保存成PDB,或者XYZ 格式的文件,其中就有坐标信息

去看看VMD所使用的几个 tcl 命令
for
set sel [atomselect top "your atoms that will be saved"]
$sel writepdb yourfilename.pdb

或者直接输出到另一个文件中
set outfile [open $yourfilename.dat w]
$sel get {x y z}
puts $outfile $sel

思路如此,去查手册搞懂各个命令,完善此脚本即可。
5楼2012-10-25 07:51:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
zh1987hs: 金币+2, 谢谢 2012-10-18 11:13:22
这个好像是lammps的弊端,因为lammps输入的时候就没有这样的定义。
2楼2012-10-17 19:04:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luxujie3

金虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流, lammps输入都是把名字数字话了,输出肯定没有,不过在把轨迹载入到VMD的时候,可以自己写一个tcl脚本,把数字对应的原子名修改过来 2012-10-18 16:55:51
我使用那么久好像也没有听说过有输出确切原子名称的lammps命令,不过在VMD查看的时候可能会有,你可以多看下VMD的使用方面的介绍或许会得到答案的
3楼2012-10-18 16:37:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

blueybz

木虫 (正式写手)

引用回帖:
2楼: Originally posted by jiaoyixiong at 2012-10-17 19:04:50
这个好像是lammps的弊端,因为lammps输入的时候就没有这样的定义。

版主,请教个问题。如何把dump文件中的坐标分别提取出来呢?比如dump中储存有1000个snapshot,有没有脚本可以把坐标分别提取出来,
4楼2012-10-25 02:42:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 十三加油 2026-03-21 4/200 2026-03-21 16:20 by barlinike
[考研] 286求调剂 +6 Faune 2026-03-21 6/300 2026-03-21 16:09 by ColorlessPI
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666!
[考研] 310求调剂 +3 baibai1314 2026-03-16 3/150 2026-03-21 03:56 by JourneyLucky
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
信息提示
请填处理意见