| 2 | 1/1 | 返回列表 |
| 查看: 513 | 回復(fù): 1 | |||
[交流]
【求助】計算結(jié)構(gòu)大小 已有1人參與
|
| 我想問一下怎么用gromacs分析聚集結(jié)構(gòu)的大小,也就是如果我想用gromacs算一下我的表面活性劑聚集后每一簇所含的單體數(shù)目,這樣的話怎么辦呢 |
捐助貴賓 (正式寫手)
|
解決這個問題分兩步: 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文件分析,舉個例子: 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)計每個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ā)表 | |
|---|---|---|---|---|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +5 | 77chaselx 2026-03-20 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 生物學(xué)調(diào)劑招人!。 +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|