计算化学公社

 找回密码 Forget password
 注册 Register
Views: 259|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[VMD] 求助VMD显示模型元素问题

[复制链接 Copy URL]

25

帖子

0

威望

93

eV
积分
118

Level 2 能力者

我建模只有CHONS五种元素,用lammps做模拟,用atomsk将.xyz文件转换成.pdb文件后,VMD却自动识别了一些LI,BE,MG,F,AL元素,若按element显示,其中占比最大的棕色是AL,BE,MG,SC,MN,B,CO,CR这些元素,而不是显示的C元素,我想只显示CHONS五种颜色进行模型的元素显示,请问该怎么做呀

881

帖子

3

威望

1664

eV
积分
2605

Level 5 (御坂)

傻傻的木瓜

2#
发表于 Post on 2025-7-8 12:40:49 | 只看该作者 Only view this author
这不是VMD的问题,根本在于当前操作绕了很多弯路。

1. 最标准的xyz文件/extxyz文件应该用元素名而非元素序号来标记原子,而extxyz的注释行也不是所有软件都能利用,比如VMD不会读注释行的盒子信息但multiwfn可以(http://sobereva.com/477);
2. lammps的dump指令的文档有"You can use the dump_modify element option to change the output from using the (numerical) atom type to an element name (or some other label). This option will help many visualization programs to guess bonds and colors. "这样一句提醒,也就是跑MD途中就要用dump_modify自行把元素改成用符号标记的格式;
3. 实际上VMD直接就能打开(1) lammps导出的.lammpstrj格式的文件(已包含盒子信息),(2) 标准的用元素名来标记原子的.xyz文件(不包含盒子信息),(3) 不标准的拿元素序号来标记原子的xyz文件(也不包含盒子信息);
4. 就算真的只能在已经导出的xyz文件上修改,不清楚atomsk怎么把元素序号映射到元素名上的话,用它还不如用linux的awk等直接做文本操作修改。
√546=23.36664289109

25

帖子

0

威望

93

eV
积分
118

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-7-8 17:12:11 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-7-8 12:40
这不是VMD的问题,根本在于当前操作绕了很多弯路。

1. 最标准的xyz文件/extxyz文件应该用元素名而非元素 ...

感谢回答,.xyz文件已按照您说的加了[size=12.573px]dump_modify[size=12.573px]命令变成标准格式,再用VMD直接打开xyz文件我看好像显示就正确了,且VMD中的color→element也识别这[size=12.573px]CHONS[size=12.573px]五种颜色了,请问这样是对的吧?

881

帖子

3

威望

1664

eV
积分
2605

Level 5 (御坂)

傻傻的木瓜

4#
发表于 Post on 2025-7-8 17:37:04 | 只看该作者 Only view this author
shiyuyang 发表于 2025-7-8 17:12
感谢回答,.xyz文件已按照您说的加了dump_modify命令变成标准格式,再用VMD直接打开xyz文件我看好像显示 ...

不关心盒子信息的话就这样也行,如果要体现三维周期性的话,得自行从别处读取当前帧盒子的三个矢量的长度a, b, c及夹角alpha, beta, gamma,然后在VMD输入指令pbc set {a b c alpha beta gamma}定义。

其实像下面这样直接让lammps输出.lammpstrj文件也可以给VMD读取并识别元素信息(把{element_list}换成按序号排列的元素符号),也不需要像.xyz文件那样另外寻找、补充盒子信息了。
  1. dump             dump1 all custom 100 traj.lammpstrj id element xu yu zu
  2. dump_modify      dump1 sort id element {element_list}
复制代码
√546=23.36664289109

25

帖子

0

威望

93

eV
积分
118

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2025-7-8 18:07:20 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-7-8 17:37
不关心盒子信息的话就这样也行,如果要体现三维周期性的话,得自行从别处读取当前帧盒子的三个矢量的长度 ...

非常感谢,这个格式信息多

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-12 20:27 , Processed in 0.207450 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list