| 查看: 1961 | 回复: 6 | |||
[求助]
dump 输出原子名称
|
|
各位好： lammps定义原子时用数字1,2,3，导致在vmd中显示的name也为1,2,3，请问怎么才能使lammps结果在vmd显示的name为真实名字名称N,H,O?  8660473.png |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有16人回复
-
求个博导看看
已经有19人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问能用DFT方法计算一个体系中原子由于所处位置不同而势能不同，能输出原子势能吗
已经有4人回复
- 从高斯激发态优化输出文件中提取原子坐标
已经有9人回复
- 我用Gaussview打开输出文件，为甚么只显示原子，键都不显示，怎么样就显示键了
已经有6人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【求助】关于LAMMPS 运行的问题
已经有6人回复
- 【转帖】极化连续介质模型（PCM）
已经有12人回复
2楼2012-10-17 19:04:50
luxujie3
金虫 (著名写手)
- 应助: 54 (初中生)
- 金币: 4957.7
- 散金: 27
- 红花: 4
- 帖子: 1214
- 在线: 266.2小时
- 虫号: 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