| 查看: 569 | 回复: 1 | |||
[交流]
【求助】计算结构大小 已有1人参与
|
| 我想问一下怎么用gromacs分析聚集结构的大小，也就是如果我想用gromacs算一下我的表面活性剂聚集后每一簇所含的单体数目，这样的话怎么办呢 |
回复此楼» 猜你喜欢
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
-
面上意见出来了
已经有23人回复
itismineok
捐助贵宾 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 537.8
- 散金: 50
- 帖子: 372
- 在线: 97.5小时
- 虫号: 425200
- 注册: 2007-07-22
- 性别: GG
- 专业: 生物大分子结构与功能
★ ★ ★
雪狼乖乖(金币+3):谢谢 希望能完整化 2010-07-28 21:08:00
雪狼乖乖(金币+3):谢谢 希望能完整化 2010-07-28 21:08:00
|
解决这个问题分两步： 1. 首先使用top2psf.pl将gromacs的 .top file（或者） .tpr file转换成.psf文件 详细见http://www.gromacs.org/Developer_Zone/Programming_Guide/VMD 2. 然后将gromacs的.trr File(或者.trj File)和转换成的.psf文件输入到VMD中 编写tcl文件分析，举个例子： set psf your.psf set dcd your.dcd mol load psf $psf dcd $dcd set outfile [open cluster.dat w]; set nf [molinfo top get numframes] # Select cluster. #set sel [atomselect top "resname polar"] # polar 表面活性剂的残基名 for {set i 1 } {$i < $nf } { incr i } { # Select cluster number. set cluster [atomselect top "resname polar and name only" frame $i] set crd [$cluster get {x y z}] #set clusters {}{} foreach a $crd b $crd { set diffr [veclength [vecsub $a $b]] #在这里根据表面活性剂分子之间的距离，生成clusters } puts $outfile "$clusters" #set cluster_num {} foreach a $clusters { #在这里统计每个cluster中分子的数量, } puts $outfile "$cluster_num" } close $outfile mol delete top [ Last edited by itismineok on 2010-7-29 at 09:17 ] |
2楼2010-07-28 17:07:52
30