| 2 | 1/1 | 返回列表 |
| 查看: 512 | 回復(fù): 1 | |||
[交流]
【求助】計(jì)算結(jié)構(gòu)大小 已有1人參與
|
| 我想問一下怎么用gromacs分析聚集結(jié)構(gòu)的大小,也就是如果我想用gromacs算一下我的表面活性劑聚集后每一簇所含的單體數(shù)目,這樣的話怎么辦呢 |
捐助貴賓 (正式寫手)
|
解決這個(gè)問題分兩步: 1. 首先使用top2psf.pl將gromacs的 .top file(或者) .tpr file轉(zhuǎn)換成.psf文件 詳細(xì)見http://www.gromacs.org/Developer_Zone/Programming_Guide/VMD 2. 然后將gromacs的.trr File(或者.trj File)和轉(zhuǎn)換成的.psf文件輸入到VMD中 編寫tcl文件分析,舉個(gè)例子: 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]] #在這里根據(jù)表面活性劑分子之間的距離,生成clusters } puts $outfile "$clusters" #set cluster_num {} foreach a $clusters { #在這里統(tǒng)計(jì)每個(gè)cluster中分子的數(shù)量, } puts $outfile "$cluster_num" } close $outfile mol delete top [ Last edited by itismineok on 2010-7-29 at 09:17 ] |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | Mqqqqqq 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +4 | 西南交通專材355 2026-03-19 | 4/200 |
|
|
[考研] 317求調(diào)劑 +4 | 申子申申 2026-03-19 | 8/400 |
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|