| 11 | 1/1 | 返回列表 |
| 查看: 3140 | 回復(fù): 10 | |||
[交流]
【求助】gromacs出現(xiàn)了一個問題
|
|||
|
做GROMACS的時候出現(xiàn)了一個問題: sys.top是這樣的: ; ; File 'sys.top' was generated ; By user: zhu (503) ; On host: apollo2 ; At date: Wed Mar 12 15:48:05 2008 ; ; This is your topology file ; drg ; ; Include forcefield parameters #include "ffG53a6.itp" ; Include surfacetant topology #include "drg.itp" ; Include water topology #include "spce.itp" #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcx fcy fcz 1 1 1000 1000 1000 #endif ; Include generic topology for ions #include "ions.itp" [ system ] ; Name DRG [ molecules ] ; Compound #mols DRG 27 SOL 6490 而建好了盒子,并且插入了水分子的 GRO文件時這樣的: http://pic.www.gaoyang168.com/file.php?id=164 但是當(dāng)grompp -f minim.mdp -c outt.gro -p sys.top -o input.tpr 的時候,出現(xiàn)了這樣的錯誤: number of coordinates in coordinate file (outt.gro, 21279) does not match topology (sys.top, 20226) |
» 搶金幣啦!回帖就可以得到:
+1/288
+2/92
+1/55
+1/40
+1/35
+1/35
+1/31
+1/28
+2/18
+2/14
+1/12
+1/10
+1/8
+1/7
+1/6
+1/5
+1/5
+1/5
+1/2
+1/1
|
itp文件: ; ; ; This file was generated by PRODRG version AA081006.0504 ; PRODRG written/copyrighted by Daan van Aalten ; and Alexander Schuettelkopf ; ; Questions/comments to dava@davapc1.bioch.dundee.ac.uk ; ; When using this software in a publication, cite: ; A. W. Schuettelkopf and D. M. F. van Aalten (2004). ; PRODRG - a tool for high-throughput crystallography ; of protein-ligand complexes. ; Acta Crystallogr. D60, 1355--1363. ; ; [ moleculetype ] ; Name nrexcl DRG 3 [ atoms ] ; nr type resnr resid atom cgnr charge mass 1 CH3 1 DRG CAA 1 0.000 15.0350 2 CH2 1 DRG CAH 2 -0.008 14.0270 3 CH2 1 DRG CAJ 2 -0.007 14.0270 4 CH2 1 DRG CAN 2 -0.008 14.0270 5 CH1 1 DRG CAY 2 0.032 13.0190 6 CH2 1 DRG CAL 2 -0.009 14.0270 7 CH3 1 DRG CAC 3 0.019 15.0350 8 CH2 1 DRG CAP 3 0.142 14.0270 9 OA 1 DRG OAT 3 -0.111 15.9994 10 C 1 DRG CAW 3 0.227 12.0110 11 O 1 DRG OAF 3 -0.387 15.9994 12 CH2 1 DRG CAR 3 0.037 14.0270 13 CH1 1 DRG CBA 3 0.073 13.0190 14 SDMSO 1 DRG SAV 4 0.600 32.0600 15 OM 1 DRG OAE 4 -0.562 15.9994 16 H 1 DRG HBD 4 -0.056 1.0080 17 CH2 1 DRG CAS 4 0.018 14.0270 18 C 1 DRG CAX 5 0.227 12.0110 19 O 1 DRG OAG 5 -0.387 15.9994 20 OA 1 DRG OAU 5 -0.111 15.9994 21 CH2 1 DRG CAQ 5 0.142 14.0270 22 CH1 1 DRG CAZ 5 0.073 13.0190 23 CH2 1 DRG CAM 5 0.037 14.0270 24 CH3 1 DRG CAD 5 0.019 15.0350 25 CH2 1 DRG CAO 6 0.006 14.0270 26 CH2 1 DRG CAK 6 0.006 14.0270 27 CH2 1 DRG CAI 6 0.006 14.0270 28 CH3 1 DRG CAB 6 -0.018 15.0350 [ bonds ] ; ai aj fu c0, c1, ... 2 1 2 0.153 7150000.0 0.153 7150000.0 ; CAH CAA 2 3 2 0.153 7150000.0 0.153 7150000.0 ; CAH CAJ 3 4 2 0.153 7150000.0 0.153 7150000.0 ; CAJ CAN 5 4 2 0.153 7150000.0 0.153 7150000.0 ; CAY CAN 5 6 2 0.153 7150000.0 0.153 7150000.0 ; CAY CAL 5 8 2 0.153 7150000.0 0.153 7150000.0 ; CAY CAP 6 7 2 0.153 7150000.0 0.153 7150000.0 ; CAL CAC 8 9 2 0.143 8180000.0 0.143 8180000.0 ; CAP OAT 10 9 2 0.136 10200000.0 0.136 10200000.0 ; CAW OAT 10 11 2 0.123 16600000.0 0.123 16600000.0 ; CAW OAF 10 12 2 0.153 7150000.0 0.153 7150000.0 ; CAW CAR 13 12 2 0.153 7150000.0 0.153 7150000.0 ; CBA CAR 13 14 2 0.195 4950000.0 0.195 4950000.0 ; CBA SAV 13 17 2 0.153 7150000.0 0.153 7150000.0 ; CBA CAS 14 15 2 0.153 8040000.0 0.153 8040000.0 ; SAV OAE 14 16 2 0.133 8870000.0 0.133 8870000.0 ; SAV HBD 18 17 2 0.153 7150000.0 0.153 7150000.0 ; CAX CAS 18 19 2 0.123 16600000.0 0.123 16600000.0 ; CAX OAG 18 20 2 0.136 10200000.0 0.136 10200000.0 ; CAX OAU 21 20 2 0.143 8180000.0 0.143 8180000.0 ; CAQ OAU 22 21 2 0.153 7150000.0 0.153 7150000.0 ; CAZ CAQ 22 23 2 0.153 7150000.0 0.153 7150000.0 ; CAZ CAM 22 25 2 0.153 7150000.0 0.153 7150000.0 ; CAZ CAO 23 24 2 0.153 7150000.0 0.153 7150000.0 ; CAM CAD 25 26 2 0.153 7150000.0 0.153 7150000.0 ; CAO CAK 26 27 2 0.153 7150000.0 0.153 7150000.0 ; CAK CAI 27 28 2 0.153 7150000.0 0.153 7150000.0 ; CAI CAB [ pairs ] ; ai aj fu c0, c1, ... 1 4 1 ; CAA CAN 2 5 1 ; CAH CAY 3 6 1 ; CAJ CAL 3 8 1 ; CAJ CAP 4 7 1 ; CAN CAC 4 9 1 ; CAN OAT 5 10 1 ; CAY CAW 6 9 1 ; CAL OAT 7 8 1 ; CAC CAP 8 11 1 ; CAP OAF 8 12 1 ; CAP CAR 9 13 1 ; OAT CBA 10 14 1 ; CAW SAV 10 17 1 ; CAW CAS 11 13 1 ; OAF CBA 12 15 1 ; CAR OAE 12 16 1 ; CAR HBD 12 18 1 ; CAR CAX 13 19 1 ; CBA OAG 13 20 1 ; CBA OAU 14 18 1 ; SAV CAX 15 17 1 ; OAE CAS 16 17 1 ; HBD CAS 17 21 1 ; CAS CAQ 18 22 1 ; CAX CAZ 19 21 1 ; OAG CAQ 20 23 1 ; OAU CAM 20 25 1 ; OAU CAO 21 24 1 ; CAQ CAD 21 26 1 ; CAQ CAK 22 27 1 ; CAZ CAI 23 26 1 ; CAM CAK 24 25 1 ; CAD CAO 25 28 1 ; CAO CAB [ angles ] ; ai aj ak fu c0, c1, ... 1 2 3 2 109.5 520.0 109.5 520.0 ; CAA CAH CAJ 2 3 4 2 109.5 520.0 109.5 520.0 ; CAH CAJ CAN 3 4 5 2 109.5 520.0 109.5 520.0 ; CAJ CAN CAY 4 5 6 2 109.5 520.0 109.5 520.0 ; CAN CAY CAL 4 5 8 2 109.5 520.0 109.5 520.0 ; CAN CAY CAP 6 5 8 2 109.5 520.0 109.5 520.0 ; CAL CAY CAP 5 6 7 2 109.5 520.0 109.5 520.0 ; CAY CAL CAC 5 8 9 2 109.5 520.0 109.5 520.0 ; CAY CAP OAT 8 9 10 2 109.5 380.0 109.5 380.0 ; CAP OAT CAW 9 10 11 2 124.0 730.0 124.0 730.0 ; OAT CAW OAF 9 10 12 2 115.0 610.0 115.0 610.0 ; OAT CAW CAR 11 10 12 2 121.0 685.0 121.0 685.0 ; OAF CAW CAR 10 12 13 2 109.5 520.0 109.5 520.0 ; CAW CAR CBA 12 13 14 2 109.5 520.0 109.5 520.0 ; CAR CBA SAV 12 13 17 2 109.5 520.0 109.5 520.0 ; CAR CBA CAS 14 13 17 2 109.5 520.0 109.5 520.0 ; SAV CBA CAS 13 14 15 2 109.5 518.0 109.5 518.0 ; CBA SAV OAE 13 14 16 2 109.5 447.3 109.5 447.3 ; CBA SAV HBD 15 14 16 2 109.5 447.3 109.5 447.3 ; OAE SAV HBD 13 17 18 2 109.5 520.0 109.5 520.0 ; CBA CAS CAX 17 18 19 2 121.0 685.0 121.0 685.0 ; CAS CAX OAG 17 18 20 2 115.0 610.0 115.0 610.0 ; CAS CAX OAU 19 18 20 2 124.0 730.0 124.0 730.0 ; OAG CAX OAU 18 20 21 2 109.5 380.0 109.5 380.0 ; CAX OAU CAQ 20 21 22 2 109.5 520.0 109.5 520.0 ; OAU CAQ CAZ 21 22 23 2 109.5 520.0 109.5 520.0 ; CAQ CAZ CAM 21 22 25 2 109.5 520.0 109.5 520.0 ; CAQ CAZ CAO 23 22 25 2 109.5 520.0 109.5 520.0 ; CAM CAZ CAO 22 23 24 2 109.5 520.0 109.5 520.0 ; CAZ CAM CAD 22 25 26 2 109.5 520.0 109.5 520.0 ; CAZ CAO CAK 25 26 27 2 109.5 520.0 109.5 520.0 ; CAO CAK CAI 26 27 28 2 109.5 520.0 109.5 520.0 ; CAK CAI CAB [ dihedrals ] ; ai aj ak al fu c0, c1, m, ... 5 4 8 6 2 35.3 334.8 35.3 334.8 ; imp CAY CAN CAP CAL 10 9 11 12 2 0.0 167.4 0.0 167.4 ; imp CAW OAT OAF CAR 13 12 17 14 2 35.3 334.8 35.3 334.8 ; imp CBA CAR CAS SAV 14 13 16 15 2 35.3 334.8 35.3 334.8 ; imp SAV CBA HBD OAE 18 17 19 20 2 0.0 167.4 0.0 167.4 ; imp CAX CAS OAG OAU 22 21 23 25 2 35.3 334.8 35.3 334.8 ; imp CAZ CAQ CAM CAO 4 3 2 1 1 0.0 5.9 3 0.0 5.9 3 ; dih CAN CAJ CAH CAA 5 4 3 2 1 0.0 5.9 3 0.0 5.9 3 ; dih CAY CAN CAJ CAH 8 5 4 3 1 0.0 5.9 3 0.0 5.9 3 ; dih CAP CAY CAN CAJ 7 6 5 4 1 0.0 5.9 3 0.0 5.9 3 ; dih CAC CAL CAY CAN 9 8 5 4 1 0.0 5.9 3 0.0 5.9 3 ; dih OAT CAP CAY CAN 5 8 9 10 1 0.0 1.3 3 0.0 1.3 3 ; dih CAY CAP OAT CAW 12 10 9 8 1 0.0 3.8 3 0.0 3.8 3 ; dih CAR CAW OAT CAP 13 12 10 9 1 0.0 1.0 6 0.0 1.0 6 ; dih CBA CAR CAW OAT 17 13 12 10 1 0.0 5.9 3 0.0 5.9 3 ; dih CAS CBA CAR CAW 12 13 14 16 1 0.0 2.9 3 0.0 2.9 3 ; dih CAR CBA SAV HBD 18 17 13 12 1 0.0 5.9 3 0.0 5.9 3 ; dih CAX CAS CBA CAR 13 17 18 20 1 0.0 1.0 6 0.0 1.0 6 ; dih CBA CAS CAX OAU 17 18 20 21 1 0.0 3.8 3 0.0 3.8 3 ; dih CAS CAX OAU CAQ 22 21 20 18 1 0.0 1.3 3 0.0 1.3 3 ; dih CAZ CAQ OAU CAX 25 22 21 20 1 0.0 5.9 3 0.0 5.9 3 ; dih CAO CAZ CAQ OAU 24 23 22 21 1 0.0 5.9 3 0.0 5.9 3 ; dih CAD CAM CAZ CAQ 26 25 22 21 1 0.0 5.9 3 0.0 5.9 3 ; dih CAK CAO CAZ CAQ 27 26 25 22 1 0.0 5.9 3 0.0 5.9 3 ; dih CAI CAK CAO CAZ 28 27 26 25 1 0.0 5.9 3 0.0 5.9 3 ; dih CAB CAI CAK CAO |
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
|
這種問題在GMX中是經(jīng)常遇到的,需要進(jìn)行簡單的手動修改,樓主請仔細(xì)檢查以下問題: 1.你在模擬的過程中是否GENION過?在include下有沒有包含加入離子的mol數(shù)? 2.如果加入離子,仔細(xì)觀察你的top文件,是不是不包含離子的坐標(biāo)?打開gro文件,觀察最下方添加的離子坐標(biāo),copy到top里,替換掉相應(yīng)數(shù)目的水,比如加了2個離子,則替換掉兩個水,即6個原子 3.加入的底物摩爾數(shù)是否正確? 思考以上三個問題,如果還未解決,我們在交流 |

| 11 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|