24小时热门版块排行榜    

查看: 1500  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

astringent

铜虫 (著名写手)

[求助] gmx不认小分子PDB残基怎么办??

初学gmx,就遇到问题了,(1)我有一个小分子的结构,用gaussian优化的,然后我就按照帖子上说的用CHEMOFFICE转化成PDB,然后想用pdb2gmx命令得到top文件,但是我使用pdb2gmx时,提示.RTP文件中没有这个小分子的残基(我用opls力场),我该怎么办,是不是需要我把这个小分子的残基手动加到rtp文件中就可以?除了修改.RTP文件外,我还需要该其他的文件吗?麻烦大家指点,如果我自己加残基的话,都需要该那些文件,大致流程是什么?细节我可以看手册。因为初学,我怕我自己写的东西出错,影响整个.RTP文件不能用啊?请问各位,有没有其他方法?
(2)另外我在看gmx自带的ethonal.itp文件时发现,所用的电荷并不是从力场中语原子类型直接对应的电荷?我想问,itp文件中的电荷可以自己随便改?
(3)如果自己在写top文件时,比如写一个angle 2 1 3 和3  1  2是不是对计算结果无影响?另外,我发现我在写top是很容易漏掉某些角度和二面角啊,大家有没有这种经历?有没有什么好方法,一次性将角度,二面角都写全?
问题有点多,请大家指点,不甚感激
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

astringent

铜虫 (著名写手)

引用回帖:
3楼: Originally posted by jawang at 2011-11-26 09:19:23:
你在gmx中用oplsaa力场的话,建议你使用软件acpype
http://code.google.com/p/acpype/
它可以直接产生可以使用的itp文件
对于电荷的话,你可以使用量化方法算出的电荷如resp取代itp中得电荷

我还一个问题请教,我看opls立场中,有电荷分组(cgnr),请问这个分组有什么依据吗?我的小分子是一个弱键化合物(A....B),那么如果我自己定义电荷分组的话,是将AB两部分作为一个整体,使整体电荷为零还是把它定义成两部分,然后A,B部分分别为零呢?我如何判断我定义的CGNR是否合理?谢谢
4楼2011-11-28 04:59:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

yanghua0378

金虫 (小有名气)

★ ★ ★
御剑江湖(金币+2): 谢谢 2011-11-26 19:43:03
jiaoyixiong(金币+1, 模拟EPI+1): 2011-11-28 20:57:13
jiaoyixiong(金币+10): 多谢指教 2011-11-28 20:58:49
cenwanglai(金币-10, 模拟EPI-1): 不能重复奖励。注意EPI授予标准。一般回帖给金币奖励就好。回答的很好的才给予EPI~ 2011-11-28 23:35:56
好像是力场的问题,要自己写才行
2楼2011-11-25 16:07:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jawang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
御剑江湖(金币+4): 谢谢 2011-11-26 19:43:11
astringent(金币+10): 2011-11-28 04:51:21
astringent(金币+5): 2011-11-29 19:51:35
你在gmx中用oplsaa力场的话,建议你使用软件acpype
http://code.google.com/p/acpype/
它可以直接产生可以使用的itp文件
对于电荷的话,你可以使用量化方法算出的电荷如resp取代itp中得电荷
3楼2011-11-26 09:19:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jawang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
ghcacj(金币+3): 谢谢 2011-11-28 12:01:22
astringent(金币+5): 2011-11-29 19:50:59
引用回帖:
4楼: Originally posted by astringent at 2011-11-28 04:59:27:
我还一个问题请教,我看opls立场中,有电荷分组(cgnr),请问这个分组有什么依据吗?我的小分子是一个弱键化合物(A....B),那么如果我自己定义电荷分组的话,是将AB两部分作为一个整体,使整体电荷为零还是把 ...

我觉得将AB作为一个整体比较合理,个人意见,因为这里会有电荷转移的问题。你可以找找相关文献对照一下。
一般是同一种类型的原子才会归为一样的电荷,比如ch3上得三个氢,你具体看自己的体系。
5楼2011-11-28 08:38:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 324分 +7 llllkkkhh 2026-03-18 7/350 2026-03-18 17:18 by 216的故事
[考研] 311求调剂 +3 冬十三 2026-03-18 3/150 2026-03-18 17:15 by 216的故事
[考研] 085600材料与化工 +5 安全上岸! 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 281求调剂(0805) +5 烟汐忆海 2026-03-16 13/650 2026-03-18 14:30 by stone_128
[考研] 085601材料工程专硕求调剂 +6 慕寒mio 2026-03-16 6/300 2026-03-18 14:26 by 007_lilei
[考研] 0703化学调剂 ,六级已过,有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
信息提示
请填处理意见