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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 【求助】VMD画晶胞
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 1325  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 +9 A child 2026-04-04 9/450 2026-04-04 23:35 by lqwchd
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 10/500 2026-04-04 23:12 by 无际的草原
[考研] 材料334求调剂 +15 Eecho# 2026-04-03 15/750 2026-04-04 23:05 by 无际的草原
[考研] 070300化学学硕311分求调剂 +5 梁富贵险中求 2026-04-04 5/250 2026-04-04 22:51 by 无际的草原
[考研] 295求调剂 +4 A你好研究生 2026-04-04 5/250 2026-04-04 22:46 by yu221
[考研] 一志愿南航，数一英一学硕317求调剂！！ +3 Acaciad 2026-04-04 3/150 2026-04-04 22:29 by lbsjt
[考研] 材料调剂 +12 一样YWY 2026-04-02 13/650 2026-04-04 20:49 by 蓝云思雨
[考研] 285求调剂 +4 AZMK 2026-04-04 5/250 2026-04-04 16:45 by cql1109
[考研] 本9一志愿2 0854低分专硕286求调剂 +9 芒种111 2026-04-04 9/450 2026-04-04 11:01 by tangruihua
[考研] 一志愿江南大学085501机械工程专硕326分，本科佳木斯大学 +4 顾若浮生 2026-04-03 8/400 2026-04-03 22:31 by zhq0425
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
[考研] 275求调剂 +14 waltzh 2026-04-01 14/700 2026-04-03 11:24 by qqq850810
[考研] 材料340分调剂 +7 夏夜晚风_long 2026-04-02 9/450 2026-04-02 21:20 by dongzh2009
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 085600，321分求调剂 +13 大馋小子 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
[考研] 085410人工智能 初试316分 求调剂 +3 残星拂曙 2026-03-31 3/150 2026-04-01 11:09 by 小熊raider
[考研] 一志愿西电085401数一英一299求调剂 六级521 +4 爱吃大鸭梨 2026-03-31 4/200 2026-03-31 11:51 by 搏击518
[考研] 一志愿华中师范化学332分求调剂 +3 Lyy930824@ 2026-03-29 3/150 2026-03-30 20:15 by DHUSHUAI
[考研] 296求调剂 +10 彼岸t 2026-03-29 10/500 2026-03-30 10:50 by 探123
信息提示
关闭
请填处理意见
关闭