24小时热门版块排行榜    

查看: 1594  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:24 by ca3DVAx4OB19
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 9/450 2026-08-02 18:23 by ca3DVAx4OB19
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 10/500 2026-08-02 18:18 by ca3DVAx4OB19
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 AS6fX1tr1on4 2026-08-02 12/600 2026-08-02 18:14 by ca3DVAx4OB19
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 13/650 2026-08-02 18:14 by ca3DVAx4OB19
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 18:12 by ca3DVAx4OB19
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 18:11 by ca3DVAx4OB19
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 18:11 by ca3DVAx4OB19
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cBGhfJFuNMca 2026-08-02 20/1000 2026-08-02 18:04 by ca3DVAx4OB19
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 20/1000 2026-08-02 18:04 by ca3DVAx4OB19
[基金申请] 准备明年的基金了 +7 Tide man 2026-07-27 7/350 2026-08-02 16:43 by gltch
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 9/450 2026-08-02 16:32 by jHcrHMhbKCFV
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[基金申请] 面上提前没消息,有中的吗 +9 archvillain 2026-08-02 11/550 2026-08-02 15:44 by alexon
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见