 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2940 | 回復(fù): 6 | ||||
| 本帖產(chǎn)生 1 個 模擬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文件時出現(xiàn)這樣的錯誤，我想請問我少了那步呢，不知道還要在什么地方添加?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 |
» 猜你喜歡
求調(diào)劑院校信息
已經(jīng)有4人回復(fù)
-
085600材料與化工306
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有4人回復(fù)
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有8人回復(fù)
-
資源與環(huán)境 調(diào)劑申請(333分)
已經(jīng)有5人回復(fù)
-
280求調(diào)劑
已經(jīng)有12人回復(fù)
-
269專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料學(xué)碩301分求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- gromacs，用pdb文件生成top文件時遇到的問題
已經(jīng)有4人回復(fù)
- topol文件不識別小分子
已經(jīng)有8人回復(fù)
- 【專家會診】bay 討論namd,gmx處理水溶液,蛋白等MD 模擬;討論水及水通道
已經(jīng)有50人回復(fù)
3樓2013-01-20 16:52:54
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
chaizhm: 金幣+5, 謝謝~ 2013-01-20 10:03:31
zh1987hs: 金幣+15, 模擬EPI+1, jiao版真是好娃啊~哈哈 2013-01-20 10:38:01
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
chaizhm: 金幣+5, 謝謝~ 2013-01-20 10:03:31
zh1987hs: 金幣+15, 模擬EPI+1, jiao版真是好娃啊~哈哈 2013-01-20 10:38:01
|
你還是應(yīng)該先非常仔細(xì)的研究一下手冊，真的！ 你的意思就是定義了一個新的 分子但是gromacs里面沒有這個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
4樓2013-01-20 16:58:17
youngfi
木蟲 (小有名氣)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 2419.9
- 散金: 15
- 紅花: 1
- 帖子: 144
- 在線: 100.9小時
- 蟲號: 1850824
- 注冊: 2012-06-07
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
6樓2013-06-19 22:55:41
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +6 | Tyoumou 2026-03-18 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +11 | 妮妮ninicgb 2026-03-15 | 15/750 |
|
|
|
[考研] 326求調(diào)劑 +4 | mlpqaz03 2026-03-15 | 4/200 |
|
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|