应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » gromacs自己添加residue遇到的问题
北京石油化工学院2026年研究生招生接收调剂公告
11/1返回列表
查看: 2962  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 模拟EPI ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jiaoyixiong

荣誉版主 (职业作家)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
chaizhm: 金币+5, 谢谢~ 2013-01-20 10:03:31
zh1987hs: 金币+15, 模拟EPI+1, jiao版真是好娃啊~哈哈 2013-01-20 10:38:01
你还是应该先非常仔细的研究一下手册，真的！

你的意思就是定义了一个新的 分子但是gromacs里面没有这个residue 吧：

下面我以charmm力场下添加 三氧化二铝的模型为例子，说明如何添加：

1、使用MS建立满足自己体系大小的三氧化二铝模型的PDB文件
2、导入到VMD中建立PSF文件
3、通过三氧化二铝的psf和pdb文件，转化成gromacs的 itp 文件，保存成alo_pore_charged.itp，并且放到charmmm27.ff文件夹内（部分内容如下）：

[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 信息根据psf文件中转化，键长，键能自己在这里设定。


4、由于我这里的使用的原子名字，不是charmm力场的标准原子名，因此需要把这些添加到 charmm力场的 atomtypes.atp （charmm27.ff目录内）下面，
; 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文件导入到gromacs中，生成。top文件就可以往下面继续了：
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 ]
中的顺序，一定要和坐标文件（pdb 或者 gro）文件中的顺序一致

ok思路就是这样，具体的多看手册，“atom N not found ”，我想主要就是第四步出现的问题吧

[ Last edited by jiaoyixiong on 2013-1-20 at 09:27 ]
赞一下(1人) 回复此楼
2楼2013-01-20 09:25:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 ilxmc 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +17 超级伊昂大王 2026-03-24 18/900 2026-03-30 21:22 by dick_runner
[考研] 297 地理学070500 复试求调剂 +3 小圆圈圈ooo 2026-03-30 3/150 2026-03-30 21:05 by 余震yz
[考研] 一志愿郑大材料工程290求调剂 +11 Youth_ 2026-03-30 11/550 2026-03-30 21:01 by peike
[考研] 281求调剂 +5 亚克西good 2026-03-26 7/350 2026-03-30 20:42 by dophin1985
[考研] 279求调剂 +12 j的立方 2026-03-29 12/600 2026-03-30 20:30 by dick_runner
[考研] 一志愿郑州大学，080500学硕，总分317分求调剂 +10 举个栗子oi 2026-03-24 11/550 2026-03-30 19:31 by michael2011
[考研] 309求调剂 +15 谁不是少年 2026-03-29 15/750 2026-03-30 19:28 by 2026材料调剂
[考研] 322求调剂 +10 宋明欣 2026-03-27 10/500 2026-03-30 18:47 by 544594351
[考研] 材料专硕调剂 +11 椰椰。 2026-03-29 11/550 2026-03-30 16:21 by wangjy2002
[考研] 环境科学与工程334分求调剂 +6 王一一依依 2026-03-30 8/400 2026-03-30 11:52 by yjolah
[考研] 085404求调剂，总分309，本科经历较为丰富 +6 来财aa 2026-03-25 6/300 2026-03-30 09:48 by 青海小西牛
[考研] 290求调剂 +3 dfffsar 2026-03-29 3/150 2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 343求调剂085601 +3 要努力学习x 2026-03-29 3/150 2026-03-29 18:35 by wxiongid
[考研] 299求调剂 +10 15188958825 2026-03-25 10/500 2026-03-29 17:51 by 王亮_大连医科大
[考研] 一志愿北京理工大学本科211材料工程294求调剂 +8 mikasa的围巾 2026-03-28 8/400 2026-03-29 12:48 by 无际的草原
[考研] 071000生物学求调剂，初试成绩343 +7 小小甜面团 2026-03-25 7/350 2026-03-28 20:25 by 唐沐儿
[考研] 求调剂 +6 芦lty 2026-03-25 7/350 2026-03-28 13:13 by 唐沐儿
[考研] 085701环境工程，267求调剂 +16 minht 2026-03-26 16/800 2026-03-28 12:16 by zllcz
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +8 曼111 2026-03-24 9/450 2026-03-28 07:58 by YYYYX1234
[考研] 一志愿上海理工能源动力（085800）310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
信息提示
关闭
请填处理意见
关闭