 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 3196 | 回復(fù): 5 | |||
liangjun44至尊木蟲 (著名寫手)
|
[求助]
【樓主分享了解決方案】：gromacs中坐標(biāo)文件與top文件不匹配 已有1人參與
|
|
大家好，剛學(xué)gromacs，在找某個(gè)教程上的例子練手，遇到個(gè)問題想請(qǐng)教一下大家。 構(gòu)建了一個(gè)含有蛋白質(zhì)、水和sds的系統(tǒng)，然后想要grompp，產(chǎn)生tpr然后在接著往下做時(shí)，遇到這樣的錯(cuò)誤。 Fatal error: number of coordinates in coordinate file (1AKI_solv.gro, 33543)does not match topology (topol.top, 33023)，（差值是520） 我檢查過幾個(gè)文件： 含結(jié)晶水的蛋白質(zhì)gro顯示有1555個(gè)原子（1321+78*3個(gè)結(jié)晶水），但topol文件以及蛋白質(zhì)的itp文件顯示只有1321個(gè)（top文件里面提示Protein_chain_A 1（應(yīng)該就是1321這個(gè)）和SOL 78（這個(gè)SOL是結(jié)晶水吧）。 我在加了溶劑之后，下面就多顯示了SOL 10376，就是共10454個(gè)水分子。加溶劑之前我有插入 -nmol 20個(gè)sds分子，sds的itp是prodrg生成的，并在top文件中include進(jìn)去且手工修改成sds 20 1AKI_solv.gro的原子數(shù)33543是對(duì)的，但我比較在意的是sds的itp。sds的gro文件里面有43個(gè)原子，但當(dāng)用prodrg生成itp后，這個(gè)itp就顯示只有17（43-17）*20剛好是520~ 我不知道問題是否出在這里，而且上面蛋白質(zhì)的itp文件里面都并不包含結(jié)晶水（不過水自己有自己的itp被include進(jìn)top吧？）不知道應(yīng)該怎樣解決這個(gè)問題，請(qǐng)大家?guī)蛶兔 下面是幾個(gè)文件的鏈接 http://pan.baidu.com/share/link?shareid=195959&uk=939815927 解決方案見第4樓。 [ Last edited by 月只藍(lán) on 2014-4-19 at 14:50 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 MDs-Gromacs |
» 猜你喜歡
一志愿070300浙大化學(xué)358分，求調(diào)劑！
已經(jīng)有4人回復(fù)
-
一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有5人回復(fù)
-
323求調(diào)劑
已經(jīng)有6人回復(fù)
-
352求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有8人回復(fù)
-
277材料科學(xué)與工程080500求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有18人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
280分求調(diào)劑 一志愿085802
已經(jīng)有7人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- gromacs主頁的蛋白質(zhì)與配體作用
已經(jīng)有5人回復(fù)
- gromacs中的top文件
已經(jīng)有16人回復(fù)
- gromacs，用pdb文件生成top文件時(shí)遇到的問題
已經(jīng)有4人回復(fù)
- GROMACS 小分子top文件
已經(jīng)有9人回復(fù)
- gromacs做全原子模擬時(shí)候，突然到MD這一步出錯(cuò)
已經(jīng)有6人回復(fù)
- 【轉(zhuǎn)帖】用Amber+Gaussian做小分子力場(chǎng)
已經(jīng)有16人回復(fù)
水兒CoCo
木蟲 (小有名氣)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 2597.7
- 紅花: 7
- 帖子: 120
- 在線: 96.8小時(shí)
- 蟲號(hào): 2718708
- 注冊(cè): 2013-10-12
- 性別: MM
- 專業(yè): 藥物化學(xué)
2樓2014-04-16 08:59:56
3樓2014-04-16 16:52:05
liangjun44
至尊木蟲 (著名寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 12270.9
- 紅花: 1
- 帖子: 1147
- 在線: 58.6小時(shí)
- 蟲號(hào): 1699935
- 注冊(cè): 2012-03-18
- 性別: GG
- 專業(yè): 藥劑學(xué)
★ ★ ★ ★ ★ ★
ben_ladeng: 金幣+1, 謝謝分享經(jīng)驗(yàn) 2014-04-19 07:27:31
月只藍(lán): 金幣+5, 感謝分享研究成果！ 2014-04-19 14:48:05
ben_ladeng: 金幣+1, 謝謝分享經(jīng)驗(yàn) 2014-04-19 07:27:31
月只藍(lán): 金幣+5, 感謝分享研究成果！ 2014-04-19 14:48:05
|
sds的gro文件里面有43個(gè)原子，但當(dāng)用prodrg生成itp后，這個(gè)itp就顯示只有17（43-17）*20剛好是520 問題就出在這里，prodrg生成的itp文件只有重原子，H都去掉，我是用了這個(gè)itp文件進(jìn)行操作；但插入sds的gro文件時(shí)使用的是全原子的結(jié)構(gòu)文件（包含氫原子），因此出現(xiàn)結(jié)構(gòu)文件和top文件不匹配。我應(yīng)該在插入sds時(shí)選擇只有重原子的gro文件，或者自己編寫sds的全原子itp文件。 |
4樓2014-04-18 23:58:55
5樓2015-01-12 22:42:10
Wineze
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1
- 帖子: 3
- 在線: 1.9小時(shí)
- 蟲號(hào): 8503259
- 注冊(cè): 2018-04-11
- 性別: MM
6樓2018-04-19 10:55:09
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 liangjun44 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 310求調(diào)劑 +4 | baibai1314 2026-03-16 | 4/200 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|