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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【楼主分享了解决方案】：gromacs中坐标文件与top文件不匹配
61/1返回列表
查看: 3249  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 303求调剂 +7 SereinQ 2026-04-10 8/400 2026-04-10 21:08 by gong120082
[考研] 材料工程302分求调剂 +18 zyx上岸！ 2026-04-04 18/900 2026-04-10 10:07 by 314126402
[考研] 293求调剂 +5 勇远库爱314 2026-04-08 5/250 2026-04-10 08:46 by vgtyfty
[考研] 297求调剂 +27 GENJIOW 2026-04-07 30/1500 2026-04-09 23:20 by wolf97
[考研] 289求调剂 +4 L1ttleTiger 2026-04-04 4/200 2026-04-09 19:53 by xiayizhi
[考研] 一志愿武理车辆 281 求调剂 +5 上岸研究生. 2026-04-07 5/250 2026-04-09 15:56 by only周
[考研] 368化学求调剂 +13 wwwwabcde 2026-04-07 14/700 2026-04-09 14:47 by heaven_jay
[考研] 化学308分求调剂 +21 你好明天你好 2026-04-07 23/1150 2026-04-08 22:32 by 凯凯要变帅
[考研] 266调剂 +8 daya sun 2026-04-07 9/450 2026-04-08 20:27 by yutian743
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸！ 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 生物学308分求调剂（一志愿华东师大） +8 相信必会光芒万� 2026-04-05 10/500 2026-04-05 12:19 by Hdyxbekcb
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 求生物学专业调剂-332分 +5 云朵遛弯指南 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
信息提示
关闭
请填处理意见
关闭