24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2979  |  回复: 2

lirry1991

铜虫 (初入文坛)

[求助] 请问应该如何进行VMD配位数着色, 已有2人参与

RT
回复此楼

» 收录本帖的淘帖专辑推荐

哎哎哎

» 猜你喜欢

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

luftf

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
参加这个帖子,里面有详细解释:
http://simulation.haotui.com/viewthread.php?tid=20716
========================================
最近论坛里许多朋友都想使用vmd来显示lammps运算中原子的辅助性质,比如将原子根据其应力、能量、温度等辅助性质进行着色,下面的帖子进行了很详细的讨论http://www.mdbbs.org/viewthread. ... ;amp;highlight=user
在这个帖子的基础上我分享一下具体的实现过程。

1>首先要在lammps的input script里设定要输出的辅助性质并将其输出。这里我们以中心对称参数csp为例
.................
compute   csp   all   centro/atom
.................
dump        pe   all   custom   1000   dump.pe   id type x y z c_csp

2>将得到的dump.csp文件使用附件里的csp2pdb.awk工具转换为PDB格式
gawk -f csp2pdb.awk dump.csp > dump.csp.pdb

3>为了使用vmd里的user选项来使原子根据csp值进行着色,需要通过附件里的pdbofactor.tcl工具将PDB格式里的bfactor转换成vmd里的user。操作如下:
打开vmd后,在vmd的命令行里依次敲入以下两行命令即可读入转换后的PDB文件:
source pdbofactor.tcl
pdbofactor dump.csp.pdb
(注意:需要先进入到包含pdbofactor.tcl和dump.csp.pdb的路径)

4>接下来就是使用user来使原子根据其csp值显示颜色了。
Graphical Representations窗口里,draw style里的coloring method选择user,这个时候图形里的原子就会根据csp值来显示不同颜色了。
你也可以在Trajectory里调整Color scale data range来调整color map以达到最好的显示效果。
下图就是我使用vmd显示的Ni3Al纳米压痕过程的一个图片。
2楼2015-11-23 19:57:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbb20121991

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
lirry1991: 金币+5, ★★★很有帮助 2015-12-04 11:16:20
楼上说的这个之外,我之前还找到了另两个处理办法跟大家分享。
http://dxli75.blog.163.com/blog/static/106768289201032004249763/
这个需要把你着色那行参数单独提取出来
http://blog.sciencenet.cn/blog-566436-823784.html
这个对应则可以直接对于lammpstrj这类轨迹文件直接计算
3楼2015-11-24 06:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lirry1991 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856专硕求调剂 希望是a区院校 +24 好好休息好不好 2026-04-09 27/1350 2026-04-13 22:22 by pies112
[考研] 22专硕求调剂 +8 haoyun上岸 2026-04-11 10/500 2026-04-13 22:15 by zhq0425
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 材料299专硕求调剂 +13 +21 2026-04-09 13/650 2026-04-13 14:16 by 张zhihao
[考研] 335求调剂 +18 想上岸呀!! 2026-04-12 19/950 2026-04-13 14:15 by 张zhihao
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +8 慕绝cc 2026-04-09 8/400 2026-04-13 14:08 by 张zhihao
[考研] 考研二轮调剂 +10 故人?? 2026-04-09 10/500 2026-04-13 09:55 by szhize
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 270求调剂 +14 杨乐369 2026-04-11 14/700 2026-04-11 20:16 by 蓝云思雨
[考研] 275求调剂 +9 1624447980 2026-04-08 10/500 2026-04-11 10:20 by Delta2012
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 一志愿京区985,085401,与本科专业一致,电子信息工程, +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 085501机械英二77总分294求调剂,接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 求调剂 +3 猪肉墩粉条cc 2026-04-08 4/200 2026-04-09 10:05 by 猪肉墩粉条cc
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 263分B区求调剂 +6 李nihao 2026-04-08 6/300 2026-04-08 09:38 by 南开小綦
信息提示
请填处理意见