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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【楼主分享了解决方案】：gromacs中坐标文件与top文件不匹配
61/1返回列表
查看: 3243  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-01-14 21:03:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水儿CoCo

木虫 (小有名气)
ben_ladeng: 见4楼 2014-04-19 07:26:51
你好，我也遇到了这种问题，请问你是怎么解决的呢
赞一下 回复此楼
2楼2014-04-16 08:59:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

870277182

铁虫 (小有名气)

【答案】应助回帖

★ ★ ★
月只蓝: 金币+1, 应助指数+1, 鼓励交流！ 2014-04-16 18:29:13
liangjun44: 金币+2, ★有帮助, 不是这样的，但感谢你的应助 2014-04-18 23:51:48
我猜你是不是在加水，加离子的时候没有-p topol.top 啊
赞一下 回复此楼
3楼2014-04-16 16:52:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangjun44

至尊木虫 (著名写手)
★ ★ ★ ★ ★ ★
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文件。
赞一下(1人) 回复此楼
4楼2014-04-18 23:58:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

背包旅客

新虫 (初入文坛)
内容已删除
回复此楼
5楼2015-01-12 22:42:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Wineze

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by 背包旅客 at 2015-01-12 22:42:10
你好，弱弱地问一句，sds是什么东西啊？我现在在做蛋白质—配体复合物的那个例子，也是用prodrg生成的itp文件，现在也遇到这个问题“Fatal error:number of coordinates in coordinate file (solv.gro, 33049） do ...

您好，请问您最后解决这个问题了吗？应该怎么做呢？打扰了。
赞一下 回复此楼
6楼2018-04-19 10:55:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liangjun44 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408 366分，本科211，一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 0858求调剂 5+5 Gky09300550， 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 一志愿0703化学招61最终排名62化学求调剂 +24 招61排名62 2026-04-07 28/1400 2026-04-10 16:15 by yx54321
[考研] 电子信息279求调剂，有书读就行 +4 wwwooden 2026-04-08 7/350 2026-04-10 11:34 by wwwooden
[考研] 初试分332，一志愿报考西北工业大学， +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[考研] 一志愿郑州大学 22408 305分求调剂 +3 安小满zzz 2026-04-08 3/150 2026-04-09 16:16 by wp06
[考研] 299求调剂 +6 ZVVZ13 2026-04-08 6/300 2026-04-09 15:48 by wp06
[考研] 1U盾记得记得就 +9 sanjin020722 2026-04-08 10/500 2026-04-09 14:11 by 诗与自由
[考研] 349学科化学045106求调剂，化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 材料科学与工程320求调剂，080500 +12 黄瓜味薯片 2026-04-06 12/600 2026-04-08 16:26 by luoyongfeng
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 求调剂 +11 wwwwabcde 2026-04-07 11/550 2026-04-07 23:16 by JourneyLucky
[考研] 319分085702安全工程求调剂 +6 rious 2026-04-05 6/300 2026-04-07 09:42 by jp9609
[考研] 一志愿青科085500，初试295分，公共课213分 +3 遇到的人愿望都� 2026-04-05 3/150 2026-04-05 18:45 by 蓝云思雨
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 358求调剂 +7 秋gk 2026-04-04 7/350 2026-04-05 13:29 by huangmoli
[考研] 调剂 +4 是可乐不是可乐 2026-04-04 4/200 2026-04-04 19:41 by 唐沐儿
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
信息提示
关闭
请填处理意见
关闭