| 查看: 2070 | 回复: 6 | ||
[求助]
dump 输出原子名称
|
|
各位好： lammps定义原子时用数字1,2,3，导致在vmd中显示的name也为1,2,3，请问怎么才能使lammps结果在vmd显示的name为真实名字名称N,H,O?  8660473.png |
回复此楼» 猜你喜欢
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有7人回复
-
求调剂
已经有20人回复
-
一志愿郑州大学085600求调剂
已经有20人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问能用DFT方法计算一个体系中原子由于所处位置不同而势能不同，能输出原子势能吗
已经有4人回复
- 从高斯激发态优化输出文件中提取原子坐标
已经有9人回复
- 我用Gaussview打开输出文件，为甚么只显示原子，键都不显示，怎么样就显示键了
已经有6人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【求助】关于LAMMPS 运行的问题
已经有6人回复
- 【转帖】极化连续介质模型（PCM）
已经有12人回复
2楼2012-10-17 19:04:50
luxujie3
金虫 (著名写手)
- 应助: 54 (初中生)
- 金币: 4964.3
- 散金: 27
- 红花: 4
- 帖子: 1216
- 在线: 266.4小时
- 虫号: 756922
- 注册: 2009-04-24
- 性别: GG
- 专业: 高分子合成化学
3楼2012-10-18 16:37:51
blueybz
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 3023.5
- 散金: 359
- 红花: 1
- 帖子: 469
- 在线: 247.7小时
- 虫号: 497342
- 注册: 2008-02-02
- 性别: GG
- 专业: 理论和计算化学
4楼2012-10-25 02:42:58
【答案】应助回帖
★ ★ ★ ★ ★
chaizhm: 金币+5, 谢谢~ 2012-10-25 09:23:05
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
【答案】应助回帖
★ ★ ★
zh1987hs: 金币+3, 谢谢 2012-10-25 19:09:03
zh1987hs: 金币+3, 谢谢 2012-10-25 19:09:03
|
昨天正好研究了下这个。你可以看下 dump_modify 。可以手动修改的。dump custom 里面可以选择输出元素 element 。 不过不修改的话都是默认的 C .(如dump 1 all custom 1 MCM.mol id x y z element )。而 dump_modify 可以修改。如有3种原子 Si O H， 可以写 dump_modify 1 element Si O H . 不过这个会不会写到自动生成的输出格式我就没试过了。另外补充一点， 自定义输出的格式定义，实际除了 %d %g 也可以用 %s 就是用来处理个，如对应上面那个dump 可以写 dump_modify 1 format "%5d %9.5g %9.5g %10s"。 如果对文件格式比较熟的话，应该可以用这个方法实现。 |
6楼2012-10-25 19:03:16
7楼2013-04-14 12:43:16
