| 查看: 2071 | 回复: 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 ] |
回复此楼» 猜你喜欢
313求调剂
已经有4人回复
-
材料调剂
已经有13人回复
-
材料专硕283求调剂
已经有13人回复
-
295求调剂
已经有7人回复
-
266分，一志愿电气工程，本科材料，求材料专业调剂
已经有7人回复
-
307分材料专业求调剂
已经有5人回复
-
一志愿北交大材料工程总分358求调剂
已经有5人回复
-
求调剂
已经有4人回复
-
材料与化工371求调剂
已经有11人回复
-
材料求调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19419.2
- 散金: 10
- 红花: 4
- 帖子: 3560
- 在线: 1035.7小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
4楼2010-01-24 22:17:00
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
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
