 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (1754)
- >考研 (119)
- >導(dǎo)師招生 (103)
- >蟲(chóng)友互識(shí) (71)
- >基金申請(qǐng) (36)
- >休閑灌水 (23)
- >碩博家園 (17)
- >博后之家 (14)
- >教師之家 (13)
- >微米和納米 (11)
- >論文投稿 (11)
- >考博 (9)
- >找工作 (9)
- >文獻(xiàn)求助 (8)
- >公派出國(guó) (7)
- >經(jīng)濟(jì)學(xué) (6)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 974 | 回復(fù): 2 | ||
zxy964777869銅蟲(chóng) (正式寫(xiě)手)
|
[求助]
editconf 和pdb2gmx坐標(biāo)轉(zhuǎn)換錯(cuò)誤
|
|
各位MD大俠， 最近準(zhǔn)備在gromacs中建聚合物的top文件，于是采用一個(gè)CF3CF2CF2CF3的結(jié)構(gòu)按照在rtp文件中添加殘基的方式建立力場(chǎng)，這一塊倒是成功了，初始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)在也不知道錯(cuò)誤出在哪里？ 執(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) 附件上傳建立聚合物力場(chǎng)的的一個(gè)初始步驟，也供初學(xué)者參考 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : gromacs模擬筆記.docx
2013-02-28 15:08:17, 14.79 K
» 猜你喜歡
一志愿華中科技大學(xué)071000，求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有4人回復(fù)
-
085600材料與化工306
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有4人回復(fù)
-
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有8人回復(fù)
-
資源與環(huán)境 調(diào)劑申請(qǐng)(333分)
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- gromacs 運(yùn)行pdb2gmx中出錯(cuò)，pdb2gmx: error while loading shared libraries
已經(jīng)有3人回復(fù)
- [新手求教]安裝Gromacs后，查看pdb2gmx命令，顯示無(wú)此命令
已經(jīng)有9人回復(fù)
zxy964777869
銅蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 19.5
- 散金: 2246
- 紅花: 8
- 帖子: 500
- 在線: 900.3小時(shí)
- 蟲(chóng)號(hào): 1355217
- 注冊(cè): 2011-07-25
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
，發(fā)現(xiàn)為文件內(nèi)部坐標(biāo)位置的錯(cuò)誤，在將pdb文件的第五列和第六列之間減少一個(gè)空格后，生成的gro文件就正常了，還是第一次遇到這種無(wú)厘頭的錯(cuò)誤�？磥�(lái)MS生成的pdb文件也應(yīng)該小心 
2樓2013-03-01 13:48:12
zxy964777869
銅蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 19.5
- 散金: 2246
- 紅花: 8
- 帖子: 500
- 在線: 900.3小時(shí)
- 蟲(chóng)號(hào): 1355217
- 注冊(cè): 2011-07-25
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
★ ★
jiaoyixiong: 金幣+2, 感謝分享經(jīng)驗(yàn)！以 我自己的經(jīng)驗(yàn)，別的軟件生成的PDB 文件，再用VMD 處理一下，比如修改你所說(shuō)的殘基名等，再倒入到gromacs中，一般就不會(huì)出錯(cuò)了。VMD 生成的PDB 文件是標(biāo)準(zhǔn)的。 2013-03-01 14:46:09
jiaoyixiong: 金幣+2, 感謝分享經(jīng)驗(yàn)！以 我自己的經(jīng)驗(yàn)，別的軟件生成的PDB 文件，再用VMD 處理一下，比如修改你所說(shuō)的殘基名等，再倒入到gromacs中，一般就不會(huì)出錯(cuò)了。VMD 生成的PDB 文件是標(biāo)準(zhǔn)的。 2013-03-01 14:46:09
|
主要錯(cuò)誤時(shí)在對(duì)pdb文件的殘基進(jìn)行修改的時(shí)候?qū)OL換成PSFE等四個(gè)字符的殘基，從而導(dǎo)致了pdb文件的數(shù)據(jù)格式出現(xiàn)混亂，造成了錯(cuò)誤。誤說(shuō)成了MS的原因，在此改正 |
3樓2013-03-01 13:59:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zxy964777869 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +6 | Tyoumou 2026-03-18 | 9/450 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +18 | 大米飯！ 2026-03-15 | 22/1100 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|