 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2944 | 回復(fù): 6 | ||||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ，點(diǎn)擊這里進(jìn)行查看 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[交流]
gromacs自己添加residue遇到的問題 已有3人參與
|
||||
|
大家好，最近一直在研究gromacs的residue這塊，我根據(jù)prodrg給出的top file自己在.rtp中添加了新的residue，具體步驟如下，在選擇的force field文件夾中的.rtp中添加residue，在ffbonded.itp里增加新的鍵的類型，最后更新residuetypes.dat，我認(rèn)為添加的內(nèi)容都挺對的，但用pdb2gmx生成top文件時(shí)出現(xiàn)這樣的錯(cuò)誤，我想請問我少了那步呢，不知道還要在什么地方添加?xùn)|西呢？ Fatal error: atom N not found in buiding block 2MOL while combining tdb and rtp For more information and tips for troubleshooting, please check the GROMACS |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 高分子學(xué)習(xí)--qlh |
» 猜你喜歡
323求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有4人回復(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ù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- gromacs，用pdb文件生成top文件時(shí)遇到的問題
已經(jīng)有4人回復(fù)
- topol文件不識別小分子
已經(jīng)有8人回復(fù)
- 【專家會診】bay 討論namd,gmx處理水溶液,蛋白等MD 模擬;討論水及水通道
已經(jīng)有50人回復(fù)
4樓2013-01-20 16:58:17
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
chaizhm: 金幣+5, 謝謝~ 2013-01-20 10:03:31
zh1987hs: 金幣+15, 模擬EPI+1, jiao版真是好娃啊~哈哈 2013-01-20 10:38:01
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
chaizhm: 金幣+5, 謝謝~ 2013-01-20 10:03:31
zh1987hs: 金幣+15, 模擬EPI+1, jiao版真是好娃啊~哈哈 2013-01-20 10:38:01
|
你還是應(yīng)該先非常仔細(xì)的研究一下手冊，真的！ 你的意思就是定義了一個(gè)新的 分子但是gromacs里面沒有這個(gè)residue 吧： 下面我以charmm力場下添加 三氧化二鋁的模型為例子，說明如何添加： 1、使用MS建立滿足自己體系大小的三氧化二鋁模型的PDB文件 2、導(dǎo)入到VMD中建立PSF文件 3、通過三氧化二鋁的psf和pdb文件，轉(zhuǎn)化成gromacs的 itp 文件，保存成alo_pore_charged.itp，并且放到charmmm27.ff文件夾內(nèi)（部分內(nèi)容如下）： [moleculetype] ; molname nrexcl al2o3Nanopore 1 [atoms] ; id type resnr residu atom cgnr charge mass 1 OAL 1 ALO O8 1 -1.05 15.9994 2 OAL 1 ALO O9 2 -1.05 15.9994 3 OAL 1 ALO O10 3 -1.05 15.9994 4 AL 1 ALO AL7 4 1.601296 26.9815 5 OAL 1 ALO O2 5 -1.05 15.9994 6 OAL 1 ALO O13 6 -1.05 15.9994 7 OAL 1 ALO O14 7 -1.05 15.9994 8 OAL 1 ALO O15 8 -1.05 15.9994 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 20062 AL 1 ALO AL4 20062 1.601296 26.9815 20063 AL 1 ALO AL5 20063 1.601296 26.9815 20064 AL 1 ALO AL6 20064 1.601296 26.9815 [bonds] ; i j funct length force.c. [angles] ; i j k funct angle force.c. 這里具體的bond angle 信息根據(jù)psf文件中轉(zhuǎn)化，鍵長，鍵能自己在這里設(shè)定。 4、由于我這里的使用的原子名字，不是charmm力場的標(biāo)準(zhǔn)原子名，因此需要把這些添加到 charmm力場的 atomtypes.atp （charmm27.ff目錄內(nèi)）下面， ; atom types below are under my definition ALS 26.9815 ; the atoms on the AL2O3 pore's surface AL 26.9815 ; AL atoms in the AL2O3 OAL 15.9994 ; O atoms in the AL2O3 5、把建立好的pdb文件導(dǎo)入到gromacs中，生成。top文件就可以往下面繼續(xù)了： genbox -cp sys_box.pdb -cs spc216.gro -o sys_water.gro -p sys.top （.top文件可以手動修改） 最后的sys.top如下： ; Include forcefield parameters #include "charmm27.ff/forcefield.itp" ; Include your topologies #include "charmm27.ff/alo_pore_charged.itp" ; Include chain topologies #include "dsdna45_DNA_chain_A.itp" #include "dsdna45_DNA_chain_B.itp" ; Include water topology #include "charmm27.ff/spc.itp" #ifdef POSRES_WATER ; Position restraint for each water oxygen [ position_restraints ] ; i funct fcx fcy fcz 1 1 1000 1000 1000 #endif ; Include topology for ions #include "charmm27.ff/ions.itp" [ system ] ; Name dsDNA in nanopore in water [ molecules ] ; Compound #mols DNA_chain_A 1 DNA_chain_B 1 al2o3Nanopore 1 SOL 88342 K 1803 CL 1847 注意注意注意： .top 中的 [ molecules ] 中的順序，一定要和坐標(biāo)文件（pdb 或者 gro）文件中的順序一致 ok思路就是這樣，具體的多看手冊，“atom N not found ”，我想主要就是第四步出現(xiàn)的問題吧 [ Last edited by jiaoyixiong on 2013-1-20 at 09:27 ] |
2樓2013-01-20 09:25:24
3樓2013-01-20 16:52:54
youngfi
木蟲 (小有名氣)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 2419.9
- 散金: 15
- 紅花: 1
- 帖子: 144
- 在線: 100.9小時(shí)
- 蟲號: 1850824
- 注冊: 2012-06-07
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
6樓2013-06-19 22:55:41
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +11 | 吃吃吃才有意義 2026-03-19 | 11/550 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|