| 查看: 2627 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】向bay__gulf求助！！ 已有4人参与
|
|||
|
目前在做碳纳米管中水的相变，请教一下怎样写.tcl脚本来计算radial density of oxygen 和 orientation distributions of water molecules within CNTs? 虽然这个问题对你可能比较简单，但对我来说却很棘手，望赐教！ 同时恭喜bay__gulf发JPCA！ 也欢迎其他高手指教！ |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
27博士申请
已经有3人回复
-
厉娜
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有12人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有15人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】高校图书馆采购集团和Elsevier公司闹崩了，以后CPL等文献怎么看呢？
已经有21人回复
- 【求助】有没有用DS软件做细胞膜的模拟的？
已经有7人回复
- 【求助】周期边界条件设定
已经有10人回复
- 【求助】有关服务器租用或者借用的事情
已经有15人回复
- 【求助】泊松过程的误差分析
已经有5人回复
- 【求助】请问量子化学的基础是“物理化学”么？
已经有18人回复
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★
jianchaoyv(金币+3): 2010-06-26 08:49:39
lei0736(金币+3):谢谢 2010-06-26 09:17:42
jianchaoyv(金币+3): 2010-06-26 08:49:39
lei0736(金币+3):谢谢 2010-06-26 09:17:42
|
请问vecadd、vecscale、vecsub、veclength这些命令在哪里可以学习一下？根本不知道它们的用途？请给一个关于这类命令的全面的学习资料！！！十分感谢！ ==== vmd ug 的第10章讲的是Vectors and Matrices 这个ug 仅仅是vmd 特有的和常用的tcl 命令, 更深入了解需要Tcl 的资料 lindex $vecd 2 #是O原子的位置矢量 === vecd 是偶极矢量, lindex $vecd 2 是偶极的z 分量 veclength $vecd 是偶极的长度 很荣幸自己的工作被别人分析和引用 |
7楼2010-06-25 16:34:06
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jianchaoyv(金币+10):谢谢指点 2010-06-22 18:29:11
ghcacj(金币+10):谢谢 2010-06-22 19:55:31
jianchaoyv(金币+10):谢谢指点 2010-06-22 18:29:11
ghcacj(金币+10):谢谢 2010-06-22 19:55:31
|
radial density, 也就是径向分布, 很常用, vmd 也提供了一个工具可以直接计算 就是measure gofr, 参考 http://www.ks.uiuc.edu/Research/vmd/current/ug/node133.html, 以及一个gui 界面, 参考 http://www.ks.uiuc.edu/Research/vmd/plugins/gofrgui/ 至于取向分布, 用户较少只能自己写了, vmd 的强大之处不在于内置了多少工具, 更不是做多漂亮的图形 而是提供了一个原子(集团)为数据结构的编程环境, 可以按照自己的需要实现任意的功能 下面是当时用的代码, 现在看来很幼稚了但拿它入门还是不错的 # count the dipolar-z angle and Legendre 2nd order polynomial in a pore, # make sure the psf/dcd be the top molecular set outf [open theta.dat w] set numFrame [molinfo top get numframes] set lgd 0 for { set fr 1 } { $fr < $numFrame } {incr fr } { set wt [atomselect top "water and name OH2 and z < 7.6 and z > -7.6 " frame $fr] set n [$wt num] set theta 0 foreach resid [$wt get resid] { set wto [atomselect top "resid $resid and name OH2" frame $fr] set wth1 [atomselect top "resid $resid and name H1" frame $fr] set wth2 [atomselect top "resid $resid and name H2" frame $fr] set coordo [measure center $wto] set coordh1 [measure center $wth1] set coordh2 [measure center $wth2] set coordh [vecadd $coordh1 $coordh2] set coordh [vecscale 0.5 $coordh] set vecd [vecsub $coordh $coordo] set thetai [expr acos ([expr [lindex $vecd 2] / [veclength $vecd]])] set thetai [expr $thetai * 180.0 / 3.1416] set theta [expr $theta + $thetai] $wto delete $wth1 delete $wth2 delete } set theta [expr $theta/$n ] set costheta [expr cos ([expr $theta*3.1416/180.0])] set lgd [expr $lgd + [expr $costheta * $costheta]] puts $fr puts $outf [format "%d %f" $fr $theta] $wt delete } set lgd [expr $lgd / $numFrame * 1.5 - 0.5] puts "3*<(cosθ)^2> - 1 = $lgd" close $outf |
2楼2010-06-22 16:44:02
3楼2010-06-22 19:41:03
|
请问上面提供的脚本中： set coordh [vecadd $coordh1 $coordh2] set coordh [vecscale 0.5 $coordh] set vecd [vecsub $coordh $coordo] set thetai [expr acos ([expr [lindex $vecd 2] / [veclength $vecd]])] 还有最后set lgd [expr $lgd / $numFrame * 1.5 - 0.5] 这几句看不懂，能不能解释一下？谢谢！ [ Last edited by jianchaoyv on 2010-6-23 at 09:56 ] |
4楼2010-06-23 09:50:22
