应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 请问应该如何进行VMD配位数着色，
51/1返回列表
查看: 3084  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

楼1970-01-01 08:00:00
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

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这类轨迹文件直接计算
赞一下(2人) 回复此楼
3楼2015-11-24 06:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

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人) 回复此楼
2楼2015-11-23 19:57:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 纯娱乐，不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭