应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 【求助】VMD画晶胞
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 1324  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

wangyujia

木虫 (正式写手)

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

我是把.xyz文件用VMD打开，这种文件只有坐标，没有晶胞矢量。我想加上晶胞矢量，并多重复几个晶胞来看晶体结构，不知道用VMD可否实现？谢谢！
回复此楼

» 猜你喜欢
1楼 2009-10-16 17:52:17
已阅   回复此楼   关注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
想重复几个自己选
赞一下(11人) 回复此楼
下苦功，三个字，一个叫下，一个叫苦，一个叫功，一定要振作精神，下苦功。
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
赞一下(12人) 回复此楼
下苦功，三个字，一个叫下，一个叫苦，一个叫功，一定要振作精神，下苦功。
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都设置成这种晶胞？
赞一下(10人) 回复此楼
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有用
赞一下(10人) 回复此楼
下苦功，三个字，一个叫下，一个叫苦，一个叫功，一定要振作精神，下苦功。
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
即可
赞一下(11人) 回复此楼
下苦功，三个字，一个叫下，一个叫苦，一个叫功，一定要振作精神，下苦功。
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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 298分 070300求调剂 +13 zwen03 2026-04-02 13/650 2026-04-04 23:13 by 无际的草原
[考研] 化学357分，考研调剂 +8 .Starry. 2026-04-04 9/450 2026-04-04 22:52 by 理想在实践之中
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-04-04 4/200 2026-04-04 22:12 by 啵啵啵0119
[考研] 302求调剂一志愿华中师范大学 +8 小江小江江江 2026-04-02 8/400 2026-04-04 19:50 by 蓝云思雨
[考研] 342求调剂 +3 Liang7111 2026-04-04 5/250 2026-04-04 19:47 by dongzh2009
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 283分求调剂 +5 试试看呗 2026-04-04 5/250 2026-04-04 17:06 by dongzh2009
[考研] 电子信息调剂交叉学科有推荐吗 +6 jhtfeybgj 2026-04-01 8/400 2026-04-04 07:52 by 1753564080
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
[考研] 303求调剂 +9 DLkz1314. 2026-03-30 9/450 2026-04-03 18:34 by ls刘帅
[考研] 266分，求材料相关专业调剂 +13 哇呼哼呼哼 2026-03-30 15/750 2026-04-03 15:24 by arrow8852
[考研] 313求调剂 +3 ～微微凉～ 2026-04-03 3/150 2026-04-03 11:25 by 啵啵啵0119
[考研] 071000生物学调剂 +8 知昭蔓 2026-04-02 8/400 2026-04-03 10:36 by macy2011
[考研] 材料340分调剂 +7 夏夜晚风_long 2026-04-02 9/450 2026-04-02 21:20 by dongzh2009
[考研] 085600，321分求调剂 +13 大馋小子 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[考研] 262求调剂 +9 励志一定发文章 2026-03-31 10/500 2026-04-01 12:22 by sunshine0013
[基金申请] 面上5B能上会吗？ +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 抱歉 +3 田洪有 2026-03-30 3/150 2026-03-30 19:11 by 迷糊CCPs
[考研] 求调剂 +10 家佳佳佳佳佳 2026-03-29 10/500 2026-03-30 18:34 by 544594351
信息提示
关闭
请填处理意见
关闭