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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » gmx不认小分子PDB残基怎么办？？
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1528  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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是很容易漏掉某些角度和二面角啊，大家有没有这种经历？有没有什么好方法，一次性将角度，二面角都写全？
问题有点多，请大家指点，不甚感激
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-11-25 05:33:21
已阅   回复此楼   关注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
好像是力场的问题，要自己写才行
赞一下(4人) 回复此楼
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中得电荷
赞一下(1人) 回复此楼
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上得三个氢，你具体看自己的体系。
赞一下(1人) 回复此楼
5楼2011-11-28 08:38:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂，一志愿 南京航空航天大学 ，080500材料科学与工程学硕，总分289分 +9 @taotao 2026-03-29 9/450 2026-03-30 22:29 by 我是小康
[考研] 085602化工求调剂（331分） +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 303求调剂 +7 DLkz1314. 2026-03-30 7/350 2026-03-30 21:07 by peike
[考研] 291求调剂 +8 Y-cap 2026-03-29 11/550 2026-03-30 21:03 by Y-cap
[考研] 抱歉 +3 田洪有 2026-03-30 3/150 2026-03-30 19:11 by 迷糊CCPs
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 317分 一志愿南理工材料工程 本科湖工大 求调剂 +12 芋泥小铃铛 2026-03-28 12/600 2026-03-30 17:06 by wangjy2002
[基金申请] 面上5B能上会吗？ +7 redcom 2026-03-29 7/350 2026-03-30 12:44 by 爱搞研究的小孩
[考研] 283求调剂（080500） +14 A child 2026-03-27 14/700 2026-03-30 12:06 by 探123
[考研] 318求调剂 +7 陈晨79 2026-03-30 7/350 2026-03-30 10:49 by 探123
[考研] 总分293求调剂 +8 加一一九 2026-03-25 11/550 2026-03-29 19:53 by 无际的草原
[考研] 343求调剂 +6 爱羁绊 2026-03-29 6/300 2026-03-29 12:00 by 无际的草原
[考研] 一志愿太原理工安全工程300分，求调剂 +5 0857求调剂. 2026-03-24 6/300 2026-03-28 22:04 by zhq0425
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 085602 307分 求调剂 +7 不知道叫什么！ 2026-03-26 7/350 2026-03-28 09:57 by 神马都不懂
[考研] 085601 材料工程 313分 求调剂 +5 Ong3 2026-03-27 5/250 2026-03-27 12:24 by goldfish51
[考研] 一志愿郑大085600，310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[考研] 一志愿吉大071010，316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 085602化学工程求调剂。 +4 平乐乐乐 2026-03-26 4/200 2026-03-26 17:57 by fmesaito
[考研] 一志愿哈工大，085400，320，求调剂 +4 gdlf9999 2026-03-24 4/200 2026-03-25 23:01 by boxking200
信息提示
关闭
请填处理意见
关闭