| 查看: 2010 | 回复: 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 ] |
回复此楼» 猜你喜欢
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
-
多组分精馏求助
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
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