计算化学公社

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

[VMD] 求助:VMD怎么进行多个分子多个frame显示

[复制链接 Copy URL]

7

帖子

0

威望

53

eV
积分
60

Level 2 能力者

请教一下各位老师,我想用VMD显示一个物理过程,这个过程有几十步(几十个frame),其中每一步包括几千个分子(C1,C2,C3组成一个分子结构,A单独组成一个分子结构),并且存在上一步与下一步分子数不相同的情况,请问这个用xyz文件能实现吗?或者有其他更好的实现方法吗?
我目前是用一个xmol文件包括每一步的所有原子,用TK console导入所有xmol文件,一开始mol new,然后在分子数不变的情况下用mol addfile的方式添加多个frame,分子数变换后就mol new,但是这样出来的效果不好,需要手动调整很麻烦,而且无法识别出每个分子结构,无法成键。
希望各位老师能够指点一下!谢谢!

202208202145113391..png (102.4 KB, 下载次数 Times of downloads: 17)

tcl文件

tcl文件

202208202143481222..png (61.31 KB, 下载次数 Times of downloads: 27)

xmol文件

xmol文件

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

8#
发表于 Post on 2022-9-5 22:23:19 | 只看该作者 Only view this author
sobereva 发表于 2022-9-5 20:59
对其中一个体系,把视角调合适了,用ppt里的sv命令保存,然后脚本里每次都调用ppt里的lv命令使视角和之前 ...

我还专门把录音的这几页听了几遍,但还是没有完全重复出来讲义上的效果,请帮忙看一下我的操作在哪里出了问题

https://www.bilibili.com/video/BV1ND4y1q7rw/
又菜又爱玩

6万

帖子

99

威望

5万

eV
积分
120064

管理员

公社社长

7#
发表于 Post on 2022-9-5 20:59:56 | 只看该作者 Only view this author
牧生 发表于 2022-9-5 17:18
非常感谢,在依葫芦画瓢的过程中,整了个半成品,现在可以顺利同时显示氢键和线条了,也有盒子信息了,可 ...

对其中一个体系,把视角调合适了,用ppt里的sv命令保存,然后脚本里每次都调用ppt里的lv命令使视角和之前设的一致




北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

6#
发表于 Post on 2022-9-5 17:18:06 | 只看该作者 Only view this author
本帖最后由 牧生 于 2022-9-5 19:25 编辑
sobereva 发表于 2022-9-4 23:15
用xyz格式也可以传递给Multiwfn盒子信息,需要自己补充,见
使用Multiwfn非常便利地创建CP2K程序的输入 ...

非常感谢,在依葫芦画瓢的过程中,整了个半成品,现在可以顺利同时显示氢键和线条了,也有盒子信息了,可以批量截多个pdb的图。。
但仍有一个问题,整了好久都没弄好。source 脚本得出的图形较小,看不清,如下所示


希望得到的图形是以下的大小,请问改在脚本里加入什么样的语句啊。








脚本如下(本外行写的,一定有很多不合理的地方),我只知道应该设置set current_view,但不会弄

#vmd

display projection Orthographic
color Display Background black
display depthcue off
menu main  on
set allfiles [glob *.pdb]
foreach file $allfiles {
mol new $file

# turn on lights 0 and 1
light 0 on
light 1 on
light 2 off
light 3 off
lappend viewplist [molinfo top]
puts "Rendering $file..."
display nearclip set 0.1
mol delrep 0 top
mol representation hbonds 3.5 35 1
mol addrep top
mol representation lines
mol addrep top
pbc box
display rendermode GLSL
lappend viewplist [molinfo top]
render TachyonInternal $file.bmp
mol delete top
}



又菜又爱玩

6万

帖子

99

威望

5万

eV
积分
120064

管理员

公社社长

5#
发表于 Post on 2022-9-4 23:15:55 | 只看该作者 Only view this author
牧生 发表于 2022-9-4 11:33
请教一下,近日正好遇到了这样一个问题,需要将多个原子总数量不同的pdb文件,做成一个轨迹文件,以便于 ...

用xyz格式也可以传递给Multiwfn盒子信息,需要自己补充,见
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587http://bbs.keinsci.com/thread-21668-1-1.html

每帧作为单独的pdb文件,写VMD脚本,自动载入,渲染成图片,然后再处理下一个。参考

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

4#
发表于 Post on 2022-9-4 11:33:40 | 只看该作者 Only view this author
本帖最后由 牧生 于 2022-9-4 13:35 编辑
sobereva 发表于 2022-8-22 09:32
VMD命令行窗口输入诸如topo readvarxyz test.xyz利用Topotools提供的功能载入,原子数变化的情况Topotools ...

请教一下,近日正好遇到了这样一个问题,需要将多个原子总数量不同的pdb文件,做成一个轨迹文件,以便于用VMD的movie maker转为视频。

比如这里,做水结冰,中间是冰,两边是液态水,我想看液态水分子变少的过程,取了其中10帧液态水为例,其中一帧为
我的方法是,使用脚本调用Multiwfn将这一批pdb文件批量转为xyz文件(有一个问题:Multiwfn不识别四点水的虚原子,MW都变成了Mg,且xyz格式,没有盒子尺寸),再用cat命令合并这一批xyz文件。
然后命令行打开vmd,输入topo readvarxyz Hyd.xyz,确实是10帧,但得到的图像全是点


需要请教解决的问题:
(1)如何才能实现10帧的总原子数量不同的pdb文件,转为一个轨迹?
(2)如果我的上述方法行不通,有没有什么方法,可以将这多个pdb文件按照一定的设置显示方法(如同时显示lines和hbond),批量转成图片,然后再ffmpeg合并为视频?
期待看到的效果是,所显示的分子数量越来越少。


我所用的这10个pdb在附件中

Liquid_frame.rar (217.58 KB, 下载次数 Times of downloads: 6)

又菜又爱玩

7

帖子

0

威望

53

eV
积分
60

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2022-8-22 12:50:43 | 只看该作者 Only view this author
sobereva 发表于 2022-8-22 09:32
VMD命令行窗口输入诸如topo readvarxyz test.xyz利用Topotools提供的功能载入,原子数变化的情况Topotools ...

感谢老师指导,我去试试

6万

帖子

99

威望

5万

eV
积分
120064

管理员

公社社长

2#
发表于 Post on 2022-8-22 09:32:34 | 只看该作者 Only view this author
VMD命令行窗口输入诸如topo readvarxyz test.xyz利用Topotools提供的功能载入,原子数变化的情况Topotools会自动用虚粒子占位
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

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

GMT+8, 2025-8-12 23:08 , Processed in 0.158450 second(s), 24 queries , Gzip On.

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