计算化学公社

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

[Sobtop] sobtop生成普通分子.top文件时出错

[复制链接 Copy URL]

21

帖子

0

威望

328

eV
积分
349

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 shiyunfan 于 2022-6-28 12:15 编辑

分子的.mol2格式如下
1047.mol2 (7.94 KB, 下载次数 Times of downloads: 19)
2246.mol2 (6.55 KB, 下载次数 Times of downloads: 20)
直接在windows环境下运行sobtop到指认原子类型那一步就会闪退
用Cmder运行sobtop时,到指认原子类型那步就会提示没有像(         1)这样的原子类型,当我把.mol2文件中原子下面空格或者是数字和原子之间空格删掉的时候,重新载入.mol2文件就会直接闪退

我把mol2文件中的原子(仅原子不带数字)从上到下输入到了.txt文件中,在生成top文件之前指认了.txt文件中的原子类型(即sobtop中的-2选项),虽然sobtop屏幕上给出的原子类型不对,但接着往下算貌似给出了正确的.top文件和.itp文件
1047.itp (55.21 KB, 下载次数 Times of downloads: 2)
1047.top (334 Bytes, 下载次数 Times of downloads: 0)

1560

帖子

0

威望

5003

eV
积分
6563

Level 6 (一方通行)

2#
发表于 Post on 2022-6-28 15:01:27 | 只看该作者 Only view this author
本帖最后由 牧生 于 2022-6-28 15:07 编辑

1047是对了的,顺利得到itp
2246问题出在分子结构不对,需要使用gv打开时自动补氢,再存为mol2格式,再用sobtop就不会错了
但是电荷还要自己算

又菜又爱玩

21

帖子

0

威望

328

eV
积分
349

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-6-28 15:49:33 | 只看该作者 Only view this author
牧生 发表于 2022-6-28 15:01
1047是对了的,顺利得到itp
2246问题出在分子结构不对,需要使用gv打开时自动补氢,再存为mol2格式,再用s ...

您好,我都是先用protoss补完氢原子然后计算完了RESP2电荷再转.gro .top .itp文件的
现在我这边的问题好像是sobtop它不能识别我.mol2文件下的原子进而没法指认它的原子类型
还有就是您方便让我看看您那边计算完了的1047的itp文件嘛

1560

帖子

0

威望

5003

eV
积分
6563

Level 6 (一方通行)

4#
发表于 Post on 2022-6-28 18:07:01 | 只看该作者 Only view this author
1047.itp (55.21 KB, 下载次数 Times of downloads: 3)

2246.mol2 (3.59 KB, 下载次数 Times of downloads: 8)
又菜又爱玩

21

帖子

0

威望

328

eV
积分
349

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2022-6-28 19:15:13 | 只看该作者 Only view this author

我把原始sdf文件在pymol中保存成了mol2格式,原子前面的空格消失了,现在分子都能正常生成top和itp文件了,我对比了一下咱俩生成的1047.itp,是完全一样的,至于那个2246,我查看了原始结构,还是用protoss补完氢后的mol文件与原始结构一致,而且现在也能正常生成top和itp文件了,所以我想应该是原子前面空格的问题,但仍然非常感谢您

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

6#
发表于 Post on 2022-6-28 19:48:25 | 只看该作者 Only view this author
你的2246.mol2文件格式不规矩
拿其它程序诸如GaussView载入后重新产生mol2就完了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

61

帖子

0

威望

812

eV
积分
873

Level 4 (黑子)

7#
发表于 Post on 2022-6-29 10:01:09 | 只看该作者 Only view this author
sobereva 发表于 2022-6-28 19:48
你的2246.mol2文件格式不规矩
拿其它程序诸如GaussView载入后重新产生mol2就完了

sob老师,我之前发现单个分子结构通过VMD生成的mol2文件载入sobtop就会闪退,通过GW倒是没问题。

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

8#
发表于 Post on 2022-6-30 00:44:12 | 只看该作者 Only view this author
五月雨 发表于 2022-6-29 10:01
sob老师,我之前发现单个分子结构通过VMD生成的mol2文件载入sobtop就会闪退,通过GW倒是没问题。

没有叫GW的程序
VMD生成的mol2文件不太规矩
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

87

帖子

0

威望

1415

eV
积分
1502

Level 5 (御坂)

9#
发表于 Post on 2022-9-14 10:15:24 | 只看该作者 Only view this author
今天发现如果mol2文件的路径太长的话,sobtop生成top的时候也可能会闪退,这各原因可以筛查一下。

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 13:13 , Processed in 0.223974 second(s), 23 queries , Gzip On.

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