计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1648|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] [已解决]求助,gromacs进行grompp计算,出现"Atomtype hc not found"报错

[复制链接 Copy URL]

43

帖子

0

威望

551

eV
积分
594

Level 4 (黑子)

本帖最后由 Tanghaoru 于 2024-1-29 16:53 编辑

各位老师好!初学gromacs使用grompp,将文件整合为tpr文件。gromacs报错MIBK.itp文件中的第十六行“Atomtype hc not found”。请问这种情况该如何解决,具体如何操作可以解决。十分感谢!MIBK.itp文件是使用sob老师的Sobtop程序并参考sob老师的方法“http://sobereva.com/635”得到的。

Snipaste_2024-01-27_16-08-32.png (7.67 KB, 下载次数 Times of downloads: 6)

gromacs的报错

gromacs的报错

MIBK.itp

11.63 KB, 下载次数 Times of downloads: 9

ECMIBK.top

235 Bytes, 下载次数 Times of downloads: 12

ECMIBK.gro

341.91 KB, 下载次数 Times of downloads: 4

EC.itp

6.48 KB, 下载次数 Times of downloads: 8

一个做催化实验的并努力向计算转行的在读水货

19

帖子

0

威望

133

eV
积分
152

Level 3 能力者

9#
发表于 Post on 2024-5-14 22:18:31 | 只看该作者 Only view this author
jwmera 发表于 2024-4-22 16:52
我也按照这个方法做的,但是还是报错,Atomtype 1 not found,和你用的一样的力场,不知道咋回事

你的atom type中没有这个原子,自己添加进去参数,或者把这个原子改名字,如果是普通氢的话

1

帖子

0

威望

101

eV
积分
102

Level 2 能力者

8#
发表于 Post on 2024-4-22 16:52:58 | 只看该作者 Only view this author
Tanghaoru 发表于 2024-1-27 22:22
好的谢谢!已排查出问题,将EC.itp和MIBK.itp中的“[ atomtypes ].......”剪切至top文件中,再把重复定 ...

我也按照这个方法做的,但是还是报错,Atomtype 1 not found,和你用的一样的力场,不知道咋回事

43

帖子

0

威望

551

eV
积分
594

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2024-1-28 11:47:28 | 只看该作者 Only view this author
牧生 发表于 2024-1-28 06:51
opls和gaff不兼容的。你应该统一用gaff+amber力场。

好的十分感谢提醒!
一个做催化实验的并努力向计算转行的在读水货

1560

帖子

0

威望

5003

eV
积分
6563

Level 6 (一方通行)

6#
发表于 Post on 2024-1-28 06:51:00 | 只看该作者 Only view this author
opls和gaff不兼容的。你应该统一用gaff+amber力场。
又菜又爱玩

250

帖子

0

威望

1464

eV
积分
1714

Level 5 (御坂)

5#
发表于 Post on 2024-1-27 22:53:03 | 只看该作者 Only view this author
Tanghaoru 发表于 2024-1-27 22:22
好的谢谢!已排查出问题,将EC.itp和MIBK.itp中的“[ atomtypes ].......”剪切至top文件中,再把重复定 ...

不客气

43

帖子

0

威望

551

eV
积分
594

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2024-1-27 22:22:30 | 只看该作者 Only view this author
UTM 发表于 2024-1-27 22:01
拓扑文件里没有“hc”这种氢原子的定义,建议检查一下itp文件。

好的谢谢!已排查出问题,将EC.itp和MIBK.itp中的“[ atomtypes ].......”剪切至top文件中,再把重复定义的原子删去,就没报错了
一个做催化实验的并努力向计算转行的在读水货

250

帖子

0

威望

1464

eV
积分
1714

Level 5 (御坂)

3#
发表于 Post on 2024-1-27 22:01:56 | 只看该作者 Only view this author
拓扑文件里没有“hc”这种氢原子的定义,建议检查一下itp文件。

43

帖子

0

威望

551

eV
积分
594

Level 4 (黑子)

2#
 楼主 Author| 发表于 Post on 2024-1-27 16:26:09 | 只看该作者 Only view this author
此外,给盒子的力场为oplsaa力场
一个做催化实验的并努力向计算转行的在读水货

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-24 11:57 , Processed in 0.196995 second(s), 30 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list