| 查看: 2067 | 回复: 13 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 linggang87 的 11 个金币 | |||
linggang87银虫 (小有名气)
|
[交流]
【求助】VASP第一原理分子动力学后的RMSD计算
|
||
|
请教各位一个问题： 用VASP做第一原理分子动力学计算之后，想根据RMSD判断体系是否平衡，但不知如何计算，忘高手不吝指教。 目前本人知道RMSD计算可利用VMD，但可能需要trajectory文件，不知如何得到。另外 VMD手册中提到了 print_rmsd_through_time ，但不知如何运用。忘熟知VMD或有过类似经验的指教。 Thanks in advance！ 参考网址 http://www.vislab.uq.edu.au/users/manuals/vmd_ug/node180.html [ Last edited by linggang87 on 2010-1-24 at 21:51 ] |
回复此楼» 猜你喜欢
333求调剂
已经有9人回复
-
070300化学学硕311分求调剂
已经有3人回复
-
321求调剂
已经有13人回复
-
348求调剂
已经有4人回复
-
0703化学
已经有8人回复
-
材料调剂
已经有11人回复
-
081700学硕，323分，一志愿中国海洋大学求调剂学校
已经有10人回复
-
0703求调剂383分
已经有7人回复
-
0703化学调剂325分
已经有4人回复
-
329求调剂
已经有15人回复
» 本主题相关商家推荐: (我也要在这里推广)
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19419.2
- 散金: 10
- 红花: 4
- 帖子: 3560
- 在线: 1035.7小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
2楼2010-01-24 21:32:17
linggang87
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 48.4
- 散金: 60
- 红花: 1
- 帖子: 198
- 在线: 108.8小时
- 虫号: 635564
- 注册: 2008-10-24
- 专业: 金属材料的制备科学与跨学
3楼2010-01-24 21:48:36
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19419.2
- 散金: 10
- 红花: 4
- 帖子: 3560
- 在线: 1035.7小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
4楼2010-01-24 22:17:00
linggang87
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 48.4
- 散金: 60
- 红花: 1
- 帖子: 198
- 在线: 108.8小时
- 虫号: 635564
- 注册: 2008-10-24
- 专业: 金属材料的制备科学与跨学
|
请高手先来帮我解决这个问题，帖子里面链接中的一个例子的运行： 参照手册中的以下部分： simulation example script Here's a longer script which you might find useful. The problem is to compute the RMSD between each frame of the simulation and the first frame. Usually in a simulation there is no initial global velocity, so the center of mass doesn't move, but because of angular rotations and because of numerical imprecisions that slowly build up, the script aligns the molecule before computing its RMSD. # Prints the RMSD of the protein atoms between each \timestep # and the first \timestep for the given molecule id (default: top) proc print_rmsd_through_time {{mol top}} { # use frame 0 for the reference set reference [atomselect $mol "protein" frame 0] # the frame being compared set compare [atomselect $mol "protein"] set num_steps [molinfo $mol get numframes] for {set frame 0} {$frame < $num_steps} {incr frame} { # get the correct frame $compare frame $frame # compute the transformation set trans_mat [measure fit $compare $reference] # do the alignment $compare move $trans_mat # compute the RMSD set rmsd [measure rmsd $compare $reference] # print the RMSD puts "RMSD of $frame is $rmsd" } } To use this, load a molecule with an animation (for example, $VMDDIR/proteins/alanin.DCD from the VMD distribution). Then run print_rmsd_through_time. Example output is shown here: vmd > print_rmsd_through_time RMSD of 0 is 0.000000 RMSD of 1 is 1.060704 RMSD of 2 is 0.977208RMSD [...] 我首先按照说明 导入$VMDDIR/proteins/alanin. 然后在console 运行了print_rmsd_through_time，但结果显示 measure fit: no atoms selected。 能不能帮我运行一下找到解决方法。不胜感激！ |
5楼2010-01-25 11:09:19
linggang87
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 48.4
- 散金: 60
- 红花: 1
- 帖子: 198
- 在线: 108.8小时
- 虫号: 635564
- 注册: 2008-10-24
- 专业: 金属材料的制备科学与跨学
6楼2010-01-27 14:47:54
freer007
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 1 (幼儿园)
- 金币: 2421.7
- 红花: 1
- 帖子: 676
- 在线: 146小时
- 虫号: 707811
- 注册: 2009-02-24
- 性别: GG
- 专业: 半导体材料
请指点 7楼2010-01-28 13:53:39
linggang87
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 48.4
- 散金: 60
- 红花: 1
- 帖子: 198
- 在线: 108.8小时
- 虫号: 635564
- 注册: 2008-10-24
- 专业: 金属材料的制备科学与跨学
★ ★
wuchenwf(金币+2):xiexie 2010-01-29 09:34
wuchenwf(金币+2):xiexie 2010-01-29 09:34
|
我用的是脚本，感觉挺方便的，发了个资源帖 http://muchong.com/bbs/viewthread.php?tid=1824241&fpage=1 请看一下吧。 |
8楼2010-01-28 16:44:20
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19419.2
- 散金: 10
- 红花: 4
- 帖子: 3560
- 在线: 1035.7小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
9楼2010-01-29 09:36:43
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19419.2
- 散金: 10
- 红花: 4
- 帖子: 3560
- 在线: 1035.7小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
10楼2010-01-29 09:38:24
