24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1381  |  回复: 8
当前主题已经存档。

wangyujia

木虫 (正式写手)

[交流] 【求助】VMD画晶胞

我是把.xyz文件用VMD打开,这种文件只有坐标,没有晶胞矢量。我想加上晶胞矢量,并多重复几个晶胞来看晶体结构,不知道用VMD可否实现?谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

childsliu

捐助贵宾 (正式写手)

二当家的

★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+3,VIP+0):谢谢 10-16 20:39
wangyujia(金币+6,VIP+0):版主们看不过去了,给你加钱了,呵呵!(当然,只是在以小虫之心度版主之腹)我最近囊中羞涩,先把这六个金币给你,等我过几天赚了钱,再给你双倍奉上。你对我的帮助实在是太大了。 10-17 11:16
如果你的晶胞是立方体的
在tkconsole中输入
moinfo top set a x
moinfo top set b y
moinfo top set c z
x , y, z 为你晶胞矢量的数值 单位为埃
vmd mian->Graphical Represention->periodic
想重复几个自己选
下苦功,三个字,一个叫下,一个叫苦,一个叫功,一定要振作精神,下苦功。
2楼2009-10-16 18:33:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)

引用回帖:
Originally posted by childsliu at 2009-10-16 18:33:
如果你的晶胞是立方体的
在tkconsole中输入
moinfo top set a x
moinfo top set b y
moinfo top set c z
x , y, z 为你晶胞矢量的数值 单位为埃
vmd mian->Graphical Represention->periodic
想重复 ...

谢谢啊!很不巧,我的是六方啊!怎么办啊?
ps:你写的命令中少了一个l,应该是molinfo。我先敲了一行,发现错了,然后根据你的命令在说明书中找到了类似的命令,试了一下,OK了!!

[ Last edited by wangyujia on 2009-10-16 at 19:45 ]
3楼2009-10-16 19:34:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

childsliu

捐助贵宾 (正式写手)

二当家的

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+2,VIP+0):谢谢 10-16 20:40
wangyujia(金币+12,VIP+0):终于赚了一点钱,关于VMD的这一点你是手把手的把我教会了,谢谢了! 10-19 08:34
加上
molinfo top set alpha a
molinfo top set beta  b
molinfo top set gamma c
a b c 为三个角 单位为度

如果嫌麻烦
直接用pbc set {x y z a b c} 就行了
我用的是VMD1.8.7
下苦功,三个字,一个叫下,一个叫苦,一个叫功,一定要振作精神,下苦功。
4楼2009-10-16 19:56:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)


yjcmwgk(金币+1,VIP+0): 11-3 20:32
引用回帖:
Originally posted by childsliu at 2009-10-16 19:56:
加上
molinfo top set alpha a
molinfo top set beta  b
molinfo top set gamma c
a b c 为三个角 单位为度

如果嫌麻烦
直接用pbc set {x y z a b c} 就行了
我用的是VMD1.8.7

摸了一段时间,我算是弄明白了。a轴一定得平行于VMD中的x轴,其它的就随便了。但是我以往建的晶胞都设置的是b轴平行于y轴啊!结果a轴和x轴差了30度,请问我这种情况能在VMD中画出合理的晶胞吗?
pbc那个命令我用不了,估计是版本的问题,我以后下一个1.8.7看看。
还有一个问题,这种设置只是对VMD中的一个frame有用,我想看晶胞的movie,能否把所有的frame都设置成这种晶胞?
5楼2009-10-16 20:33:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

childsliu

捐助贵宾 (正式写手)

二当家的

★ ★
yjcmwgk(金币+2,VIP+0): 10-16 21:19
利用transvecinv和move命令 把a轴移动到平行x轴
如果所有的frame都在一个molid里
这种设置当然对所有frame有用
下苦功,三个字,一个叫下,一个叫苦,一个叫功,一定要振作精神,下苦功。
6楼2009-10-16 21:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)

引用回帖:
Originally posted by childsliu at 2009-10-16 21:14:
利用transvecinv和move命令 把a轴移动到平行x轴
如果所有的frame都在一个molid里
这种设置当然对所有frame有用

谢谢啊!可是我的问题还是没有完全解决。
1、move的语法是什么啊?我看到说明书中说是move 4*4 matrix,可是把move和相应的matrix敲进去之后它却说invalid command name "move",对于1.8.7也是这样。
2、VMD中的坐标转换矩阵怎么是四维的啊?第四维是什么?
3、对于加晶胞的操作真的只能对一个frame有用,要不你试试我的movie?MnS的初基单胞,a=b=c=3.738,alpha=beta=gamma=60度。有劳了,谢谢!

[ Last edited by wangyujia on 2009-11-3 at 16:06 ]
7楼2009-10-17 15:19:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

childsliu

捐助贵宾 (正式写手)

二当家的

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+1,VIP+0): 10-18 13:23
对于晶胞我不太熟悉 就先假定a轴矢量为{0 0.5 1.0}为例
tkconsole:
CODE:
set all [atomselect top all]
set nf [molinfo top get numframes]
for {set i 0} {$i < $nf} {incr i} {
     $all move [transvecinv {0 0.5 1}]
}

以上将a轴平行于x轴
如果情况不对 可以继续讨论

对所有frames可以周期性显示的问题
建议安装vmd1.8.7
tkconsole:
pbc set {a b c alpha beta gamma} -all
即可
下苦功,三个字,一个叫下,一个叫苦,一个叫功,一定要振作精神,下苦功。
8楼2009-10-17 22:02:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)

引用回帖:
Originally posted by childsliu at 2009-10-17 22:02:
对于晶胞我不太熟悉 就先假定a轴矢量为{0 0.5 1.0}为例
tkconsole:

[code]
set all [atomselect top all]
set nf [molinfo top get numframes]
for {set i 0} {$i < $nf} {incr i} {
     $all move [ ...

这个问题对我来说算是解决了,自己编了个小程序,把movie.xyz文件改了,呵呵!

Hello!还要麻烦你一下啊!呵呵!

以下命令只能对一个frame有用。我也很奇怪,看你写的命令,应该对所有的frame都有用啊!!
set all [atomselect top all]
set nf [molinfo top get numframes]
for {set i 0} {$i < $nf} {incr i} {
     $all move [transvecinv {0 0.5 1} ] }

还有一个问题:VMD的默认显示是用line型,这可不可以改啊!每次打开VMD都要把graphical representation里的东西改一下。

[ Last edited by wangyujia on 2009-11-5 at 17:36 ]
9楼2009-11-03 16:05:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangyujia 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26年申博自荐-计算机视觉 20+3 光亮我自己 2026-05-01 4/200 2026-05-03 09:01 by lfy8008
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-03 07:22 by jmr1kibhvw
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 qs8kylt1wr 2026-05-01 4/200 2026-05-03 07:19 by jmr1kibhvw
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 d8tbypao1q 2026-05-01 4/200 2026-05-03 05:49 by w5idbsc3mz
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-03 05:07 by jmr1kibhvw
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 4/200 2026-05-03 04:34 by w5idbsc3mz
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 4/200 2026-05-03 04:19 by w5idbsc3mz
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 4/200 2026-05-03 03:49 by w5idbsc3mz
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 4/200 2026-05-03 03:04 by w5idbsc3mz
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k17zww5nna 2026-05-01 5/250 2026-05-03 01:49 by w5idbsc3mz
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 qs8kylt1wr 2026-05-01 4/200 2026-05-03 00:00 by 火炎焱燚炙烤
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:51 by zahslv263y
[博后之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 qs8kylt1wr 2026-05-01 3/150 2026-05-02 15:06 by zahslv263y
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:21 by zahslv263y
[教师之家] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 d8tbypao1q 2026-05-01 3/150 2026-05-02 14:06 by zahslv263y
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k17zww5nna 2026-05-01 5/250 2026-05-02 09:21 by dz7y1p1i98
[考博] 药化及相关博士的申请 5+3 中药学渣 2026-04-29 3/150 2026-05-01 18:14 by babero
[论文投稿] 一篇MDPI论文改变了学习工作和生活 +3 helsinglee 2026-04-29 4/200 2026-05-01 11:45 by bobvan
[考博] 中国地质大学(北京)博士招生补录,数理学院材料科学与工程专业和材料与化工专业 +3 liuxh329 2026-04-29 3/150 2026-04-30 15:42 by 顺利毕业zs
[基金申请] 收到国自然专家邀请后几年才会有本子送过来评 +3 barley5 2026-04-30 3/150 2026-04-30 14:58 by aspect3000
信息提示
请填处理意见