| 查看: 3086 | 回复: 5 | |||
liangjun44至尊木虫 (著名写手)
|
[求助]
【楼主分享了解决方案】：gromacs中坐标文件与top文件不匹配 已有1人参与
|
|
大家好，刚学gromacs，在找某个教程上的例子练手，遇到个问题想请教一下大家。 构建了一个含有蛋白质、水和sds的系统，然后想要grompp，产生tpr然后在接着往下做时，遇到这样的错误。 Fatal error: number of coordinates in coordinate file (1AKI_solv.gro, 33543)does not match topology (topol.top, 33023)，（差值是520） 我检查过几个文件： 含结晶水的蛋白质gro显示有1555个原子（1321+78*3个结晶水），但topol文件以及蛋白质的itp文件显示只有1321个（top文件里面提示Protein_chain_A 1（应该就是1321这个）和SOL 78（这个SOL是结晶水吧）。 我在加了溶剂之后，下面就多显示了SOL 10376，就是共10454个水分子。加溶剂之前我有插入 -nmol 20个sds分子，sds的itp是prodrg生成的，并在top文件中include进去且手工修改成sds 20 1AKI_solv.gro的原子数33543是对的，但我比较在意的是sds的itp。sds的gro文件里面有43个原子，但当用prodrg生成itp后，这个itp就显示只有17（43-17）*20刚好是520~ 我不知道问题是否出在这里，而且上面蛋白质的itp文件里面都并不包含结晶水（不过水自己有自己的itp被include进top吧？）不知道应该怎样解决这个问题，请大家帮帮忙~ 下面是几个文件的链接 http://pan.baidu.com/share/link?shareid=195959&uk=939815927 解决方案见第4楼。 [ Last edited by 月只蓝 on 2014-4-19 at 14:50 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 MDs-Gromacs |
» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- gromacs主页的蛋白质与配体作用
已经有5人回复
- gromacs中的top文件
已经有16人回复
- gromacs，用pdb文件生成top文件时遇到的问题
已经有4人回复
- GROMACS 小分子top文件
已经有9人回复
- gromacs做全原子模拟时候，突然到MD这一步出错
已经有6人回复
- 【转帖】用Amber+Gaussian做小分子力场
已经有16人回复
水儿CoCo
木虫 (小有名气)
- 应助: 12 (小学生)
- 金币: 2597.7
- 红花: 7
- 帖子: 120
- 在线: 96.8小时
- 虫号: 2718708
- 注册: 2013-10-12
- 性别: MM
- 专业: 药物化学
2楼2014-04-16 08:59:56
3楼2014-04-16 16:52:05
liangjun44
至尊木虫 (著名写手)
- 应助: 5 (幼儿园)
- 金币: 12270.9
- 红花: 1
- 帖子: 1147
- 在线: 58.6小时
- 虫号: 1699935
- 注册: 2012-03-18
- 性别: GG
- 专业: 药剂学
★ ★ ★ ★ ★ ★
ben_ladeng: 金币+1, 谢谢分享经验 2014-04-19 07:27:31
月只蓝: 金币+5, 感谢分享研究成果！ 2014-04-19 14:48:05
ben_ladeng: 金币+1, 谢谢分享经验 2014-04-19 07:27:31
月只蓝: 金币+5, 感谢分享研究成果！ 2014-04-19 14:48:05
|
sds的gro文件里面有43个原子，但当用prodrg生成itp后，这个itp就显示只有17（43-17）*20刚好是520 问题就出在这里，prodrg生成的itp文件只有重原子，H都去掉，我是用了这个itp文件进行操作；但插入sds的gro文件时使用的是全原子的结构文件（包含氢原子），因此出现结构文件和top文件不匹配。我应该在插入sds时选择只有重原子的gro文件，或者自己编写sds的全原子itp文件。 |
4楼2014-04-18 23:58:55
5楼2015-01-12 22:42:10
6楼2018-04-19 10:55:09