 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 975 | 回復(fù): 2 | ||
zxy964777869銅蟲 (正式寫手)
|
[求助]
editconf 和pdb2gmx坐標(biāo)轉(zhuǎn)換錯誤
|
|
各位MD大俠， 最近準(zhǔn)備在gromacs中建聚合物的top文件，于是采用一個CF3CF2CF2CF3的結(jié)構(gòu)按照在rtp文件中添加殘基的方式建立力場，這一塊倒是成功了，初始pdb坐標(biāo)為： HETATM 1 C1 PSFB 1 -5.110 2.179 -2.411 C HETATM 2 C2 PSFB 1 -3.581 2.024 -2.520 C HETATM 3 F11 PSFB 1 -5.771 0.846 -2.466 F HETATM 4 F12 PSFB 1 -5.595 3.012 -3.547 F HETATM 5 F13 PSFB 1 -5.463 2.856 -1.130 F HETATM 6 F21 PSFB 1 -2.964 3.378 -2.469 F HETATM 7 F22 PSFB 1 -3.253 1.384 -3.825 F HETATM 8 C1 PSFE 2 -3.033 1.154 -1.366 C HETATM 9 F11 PSFE 2 -3.627 -0.210 -1.445 F HETATM 10 F12 PSFE 2 -3.386 1.771 -0.058 F HETATM 11 C2 PSFE 2 -1.498 1.031 -1.450 C HETATM 12 F21 PSFE 2 -1.099 0.398 -2.737 F HETATM 13 F22 PSFE 2 -1.019 0.186 -0.321 F HETATM 14 F23 PSFE 2 -0.871 2.379 -1.344 F END 執(zhí)行 editconf -f psf.pdb -o psf.gro 后得到的gro文件如下： S C A M O R G 14 0PSFB C1 1 -0.511 0.000 0.900 0PSFB C2 2 -0.358 0.100 0.400 0PSFB F11 3 -0.577 0.100 0.600 0PSFB F12 4 -0.559 0.500 0.200 0PSFB F13 5 -0.546 0.300 0.600 0PSFB F21 6 -0.296 0.400 0.800 0PSFB F22 7 -0.325 0.300 0.400 0PSFE C1 8 -0.303 0.300 0.400 0PSFE F11 9 -0.362 0.700 0.000 0PSFE F12 10 -0.338 0.600 0.100 0PSFE C2 11 -0.149 0.800 0.100 0PSFE F21 12 -0.109 0.900 0.800 0PSFE F22 13 -0.101 0.900 0.600 0PSFE F23 14 -0.087 0.100 0.900 0.00000 0.00000 0.00000 ~ 可以看到pdb和gro文件中y和z坐標(biāo)值完全不一致，現(xiàn)在也不知道錯誤出在哪里？ 執(zhí)行pdb2gmx的屏幕輸出為 Reading test3.pdb... WARNING: all CONECT records are ignored Read 14 atoms Analyzing pdb file Splitting chemical chains based on TER records or chain id changing. There are 1 chains and 0 blocks of water and 0 residues with 14 atoms chain #res #atoms 1 ' ' 2 14 WARNING: there were 3 atoms with zero occupancy and 9 atoms with occupancy unequal to one (out of 14 atoms). Check your pdb file. Opening force field file ./oplsaa_nafion.ff/atomtypes.atp Atomtype 1 Reading residue database... (oplsaa_nafion) Opening force field file ./oplsaa_nafion.ff/aminoacids.rtp Residue 64 Sorting it all out... Opening force field file ./oplsaa_nafion.ff/aminoacids.hdb Opening force field file ./oplsaa_nafion.ff/aminoacids.n.tdb Opening force field file ./oplsaa_nafion.ff/aminoacids.c.tdb Back Off! I just backed up topol.top to ./#topol.top.1# Processing chain 1 (14 atoms, 2 residues) Warning: Starting residue PSFB0 in chain not identified as Protein/RNA/DNA. Warning: Starting residue PSFE0 in chain not identified as Protein/RNA/DNA. Problem with chain definition, or missing terminal residues. This chain does not appear to contain a recognized chain molecule. If this is incorrect, you can edit residuetypes.dat to modify the behavior. 8 out of 8 lines of specbond.dat converted successfully Checking for duplicate atoms.... Now there are 2 residues with 14 atoms Making bonds... Warning: Long Bond (1-2 = 0.323042 nm) Warning: Long Bond (1-3 = 0.861571 nm) Warning: Long Bond (1-4 = 0.425705 nm) Warning: Long Bond (1-5 = 0.532362 nm) Warning: Long Bond (5-6 = 0.503829 nm) Warning: Long Bond (8-9 = 0.568754 nm) Warning: Long Bond (8-10 = 0.425705 nm) Warning: Long Bond (8-11 = 0.603089 nm) Warning: Long Bond (11-12 = 0.708237 nm) Warning: Long Bond (11-13 = 0.512156 nm) Warning: Long Bond (11-14 = 1.06482 nm) Number of bonds was 14, now 13 Generating angles, dihedrals and pairs... Before cleaning: 27 pairs Before cleaning: 27 dihedrals Keeping all generated dihedrals Making cmap torsions...There are 27 dihedrals, 0 impropers, 24 angles 27 pairs, 13 bonds and 0 virtual sites Total mass 238.028 a.m.u. Total charge -0.162 e Writing topology Back Off! I just backed up posre.itp to ./#posre.itp.1# Writing coordinate file... --------- PLEASE NOTE ------------ You have successfully generated a topology from: test3.pdb. The Oplsaa_nafion force field is used. --------- ETON ESAELP ------------ gcq#121: "Step On the Brakes" (2 Unlimited) 附件上傳建立聚合物力場的的一個初始步驟，也供初學(xué)者參考 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : gromacs模擬筆記.docx
2013-02-28 15:08:17, 14.79 K
» 猜你喜歡
289材料與化工（085600）B區(qū)求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿武理材料工程348求調(diào)劑
已經(jīng)有6人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑
已經(jīng)有6人回復(fù)
-
323求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有5人回復(fù)
-
317求調(diào)劑
已經(jīng)有16人回復(fù)
-
287求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
308求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- gromacs 運行pdb2gmx中出錯，pdb2gmx: error while loading shared libraries
已經(jīng)有3人回復(fù)
- [新手求教]安裝Gromacs后，查看pdb2gmx命令，顯示無此命令
已經(jīng)有9人回復(fù)
zxy964777869
銅蟲 (正式寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 19.5
- 散金: 2246
- 紅花: 8
- 帖子: 500
- 在線: 900.3小時
- 蟲號: 1355217
- 注冊: 2011-07-25
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
，發(fā)現(xiàn)為文件內(nèi)部坐標(biāo)位置的錯誤，在將pdb文件的第五列和第六列之間減少一個空格后，生成的gro文件就正常了，還是第一次遇到這種無厘頭的錯誤�？磥鞰S生成的pdb文件也應(yīng)該小心 
2樓2013-03-01 13:48:12
zxy964777869
銅蟲 (正式寫手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 19.5
- 散金: 2246
- 紅花: 8
- 帖子: 500
- 在線: 900.3小時
- 蟲號: 1355217
- 注冊: 2011-07-25
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
3樓2013-03-01 13:59:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zxy964777869 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +5 | Zhao anqi 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 17/850 |
|
|
|
[考研] 考研調(diào)劑 +4 | 來好運來來來 2026-03-21 | 4/200 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|