24小时热门版块排行榜    

查看: 2954  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 306求调剂 +6 chuanzhu川烛 2026-03-18 6/300 2026-03-22 19:39 by 公瑾逍遥
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 384求调剂 +3 子系博 2026-03-22 4/200 2026-03-22 11:04 by 搏击518
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666!
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
信息提示
请填处理意见