计算化学公社

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

[Sobtop] 关于Sobtop中bonded_param.dat空行引起bonds和angle指认失败的疑似bug

[复制链接 Copy URL]

3

帖子

0

威望

53

eV
积分
56

Level 2 能力者

本帖最后由 bubaixiji 于 2026-7-12 19:09 编辑

在使用sobtop(包括Version 2026.1.16和更早的dev5的win和linux版本)生成.top文件时,发现即使选择了 3 Prebuilt bonded parameters if possible, missing ones are supplemented by mSeminario method,生成的.itp文件bonds和angle字段中的参数基本全为mSeminario method,而没有prebuild的部分;而选择 4 Same as option 3, but missing ones are arbitrarily guessed,生成的.itp文件bonds和angle均为arbitrarily guessed。

经过排查,该问题是本人利用sobtop生成石英孔隙力场文件时,向bonded_param.dat添加了新的键参数和角参数,而为了区分自己添加的参数和自带的参数,本人空了一行(LJ_param.dat中有空行是不会出现问题的,所以本人想当然以为bonded_param.dat也能空行),这也许导致sobtop无法识别空行后所有的键类型和角的类型,进而只能选择全部基于mSeminario method计算或者 arbitrarily guessed。删去空行后上述问题得到解决。

bonded_param.dat (260.96 KB, 下载次数 Times of downloads: 0) LJ_param.dat (19.4 KB, 下载次数 Times of downloads: 0) 102.mol2 (26.47 KB, 下载次数 Times of downloads: 0)
以上是产生问题的bonded_param.dat文件,以及对应的LJ_param.dat文件,和所使用的.mol2文件,以供大佬们复现,文件主要是添加了INTERFACE力场适于Amber力场的石英参数。

也希望新手在使用sobtop时引以为戒。



本版积分规则 Credits rule

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

GMT+8, 2026-7-13 08:58 , Processed in 0.218445 second(s), 24 queries , Gzip On.

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