计算化学公社

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

[VMD] vmd中怎么显示模拟盒子的边界

[复制链接 Copy URL]

18

帖子

0

威望

306

eV
积分
324

Level 3 能力者

请问各位老师,我把MS中得到的xtd文件转化为xyz格式后导入到VMD中,出现了下面图中显示的情况。想请教各位老师
1、怎么在VMD中显示周期性盒子边界
2、VMD中出现粒子跑散的现象,怎么只显示一个周期内的粒子,也就是说怎么得到和MS中显示的最终构象相似的图像

GHXY%W]7E`XIFP~0LK9P%@1.png (50.12 KB, 下载次数 Times of downloads: 78)

导入VMD后显示的初始构象

导入VMD后显示的初始构象

0(}U(J[UCKCLTM8Z}QOWLFD.png (74.96 KB, 下载次数 Times of downloads: 100)

MS中的初始构象

MS中的初始构象

~U7O1K3NXXKBS}5V9}8SQCE.png (96.13 KB, 下载次数 Times of downloads: 78)

MS中的最终构象

MS中的最终构象

MSQBPS$$EM350(}3[UZWBRB.png (45.07 KB, 下载次数 Times of downloads: 75)

导入VMD后显示的最终构象

导入VMD后显示的最终构象

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

2#
发表于 Post on 2016-6-17 22:29:33 | 只看该作者 Only view this author
没办法,xyz格式不包含盒子信息,VMD里显示不了边界
除非通过特殊手段,写perl脚本,把MS里每一帧盒子信息导出来再利用脚本读入VMD,脚本得你自己写。

通过修改MS的xtd->xyz转化脚本,应该有办法把考虑了PBC的轨迹导出来。可以看看help琢磨琢磨
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

18

帖子

0

威望

306

eV
积分
324

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2016-6-17 23:05:47 | 只看该作者 Only view this author
sobereva 发表于 2016-6-17 22:29
没办法,xyz格式不包含盒子信息,VMD里显示不了边界
除非通过特殊手段,写perl脚本,把MS里每一帧盒子信息 ...

sob老师,我用pbc set  {a b c alpha beta gamma} -all
pbc box -onpbc wrap -all
添加了盒子并将分子wrap到一个盒子中,但出现了一个新问题,周期性边界上的水分子上下连了起来,想请教下该怎么解决

7O[N)501_W@_W8~V6WCG)4I.png (56.12 KB, 下载次数 Times of downloads: 96)

7O[N)501_W@_W8~V6WCG)4I.png

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

4#
发表于 Post on 2016-6-17 23:08:39 | 只看该作者 Only view this author
qw7706885er 发表于 2016-6-17 23:05
sob老师,我用pbc set  {a b c alpha beta gamma} -all
pbc box - ...

drawing method用dynamic bonds
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

18

帖子

0

威望

306

eV
积分
324

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2016-6-18 10:39:51 | 只看该作者 Only view this author
sobereva 发表于 2016-6-17 23:08
drawing method用dynamic bonds

老师,我换成dynamicbonds后,发现很多水分子直接连在一起了,调节distance cutoff得到的结果也不太理想

I1TUQCDAMUTXMYJFQ$1HWYL.png (184.26 KB, 下载次数 Times of downloads: 79)

I1TUQCDAMUTXMYJFQ$1HWYL.png

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

6#
发表于 Post on 2016-6-18 10:51:51 | 只看该作者 Only view this author
qw7706885er 发表于 2016-6-18 10:39
老师,我换成dynamicbonds后,发现很多水分子直接连在一起了,调节distance cutoff得到的结果也不太理想
...

distance cutoff调节不管用那就没其它办法了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

18

帖子

0

威望

306

eV
积分
324

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2016-6-18 11:12:03 | 只看该作者 Only view this author
sobereva 发表于 2016-6-18 10:51
distance cutoff调节不管用那就没其它办法了

老师 我想问下有没有办法分开设置甲烷分子中C-H键和水分子中O-H键的distance cutoff

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

8#
发表于 Post on 2016-6-18 12:40:12 | 只看该作者 Only view this author
qw7706885er 发表于 2016-6-18 11:12
老师 我想问下有没有办法分开设置甲烷分子中C-H键和水分子中O-H键的distance cutoff

点数字键1,然后点不同的分子,如果命令行窗口里显示的residue是不同的,而且结构文件里甲烷和是分子的分开记录的,比如水分子是1-1000,甲烷是1001-1500分子,那么representation你建立两个,一个是residue 1 to 1000,一个是residue 1001 to 1500。分别设distance cutoff
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

18

帖子

0

威望

306

eV
积分
324

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2016-6-18 16:01:43 | 只看该作者 Only view this author
sobereva 发表于 2016-6-18 12:40
点数字键1,然后点不同的分子,如果命令行窗口里显示的residue是不同的,而且结构文件里甲烷和是分子的分 ...

恩恩 知道了 谢谢老师

41

帖子

0

威望

964

eV
积分
1005

Level 4 (黑子)

10#
发表于 Post on 2016-6-21 09:00:24 | 只看该作者 Only view this author
1.如果ms中运行的是nvt
在vmd中pbc wrap -center unitcell -compound resid -all
然后pbc set {a b c}
再pbc box
2. 如果ms中运行的是npt
从ms中把输出盒子尺寸信息修改成xst格式
在vmd中pbc readxst file.xst -all
pbc wrap ...
pbc box

以上两种情况我都测试过,没有问题,关键是pbc wrap 中的-compound选项可以保证在wrap
过程中resid不发生分离

评分 Rate

参与人数
Participants 2
eV +6 收起 理由
Reason
qw7706885er + 4 谢谢
sobereva + 2

查看全部评分 View all ratings

41

帖子

0

威望

964

eV
积分
1005

Level 4 (黑子)

11#
发表于 Post on 2016-6-21 09:03:19 | 只看该作者 Only view this author
顺序写反了
先pbc set 再pbc wrap
手机回复比较麻烦,见谅

18

帖子

0

威望

306

eV
积分
324

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2016-7-1 18:28:51 | 只看该作者 Only view this author
tobeant 发表于 2016-6-21 09:03
顺序写反了
先pbc set 再pbc wrap
手机回复比较麻烦,见谅

谢谢老师

18

帖子

0

威望

306

eV
积分
324

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2016-7-1 18:38:53 | 只看该作者 Only view this author
tobeant 发表于 2016-6-21 09:03
顺序写反了
先pbc set 再pbc wrap
手机回复比较麻烦,见谅

请问老师,在MS里将输出盒子尺寸信息修改成xst格式应该怎么改,刚学不太懂这些

41

帖子

0

威望

964

eV
积分
1005

Level 4 (黑子)

14#
发表于 Post on 2016-7-4 13:59:20 | 只看该作者 Only view this author
qw7706885er 发表于 2016-7-1 18:38
请问老师,在MS里将输出盒子尺寸信息修改成xst格式应该怎么改,刚学不太懂这些

xst就是记录了轨迹中每一帧的盒子信息 格式非常简单 从ms里面复制出来简单修改即可 xst具体格式信息可以看NAMD软件的手册

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

3

帖子

0

威望

79

eV
积分
82

Level 2 能力者

15#
发表于 Post on 2022-3-23 16:28:31 | 只看该作者 Only view this author
怎么修改盒子的框线呀?

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 17:50 , Processed in 0.446694 second(s), 31 queries , Gzip On.

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