计算化学公社

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

[GROMACS] 求助,使用python脚本进行配体str文件转prm文件时报编码错误

[复制链接 Copy URL]

2

帖子

0

威望

11

eV
积分
13

Level 1 能力者

最近在学习使用gromacs进行分子动力学模拟,在使用python脚本cgenff_charmm2gmx_py3_nx2.py将配体str文件转成prm文件时一直报图中的编码错误,我已经按照教程将脚本文件里f = open(str_filename, 'r')改为f = open(str_filename, 'r',encoding='utf-8')了,但是还是显示这个报错。
请教各位大佬,还有什么方法可以解决这个问题吗?谢谢

1734234008339.jpg (51.38 KB, 下载次数 Times of downloads: 0)

1734234008339.jpg

cgenff_charmm2gmx_py3_nx2.py

36.97 KB, 下载次数 Times of downloads: 0

1240

帖子

6

威望

2553

eV
积分
3913

Level 5 (御坂)

傻傻的木瓜

2#
发表于 Post on 2024-12-15 11:50:28 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2024-12-16 21:36 编辑

是哪个教程说的用GROMACS做分子动力学模拟需要“将配体str文件转成prm文件”?看社长这条评论

编辑:这帖没有达到预期目的,楼主没回答所看教程的名称而是又在上面链接的楼层中问了一遍。
√546=23.36664289109

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

3#
发表于 Post on 2024-12-16 02:50:05 | 只看该作者 Only view this author
GROMACS标准的使用过程中既不牵扯到叫str的东西也不牵扯到叫prm的东西,你去搜GROMACS官方手册也搜不到这种文件。如果某些教程作者自己发明了这些偏门的东西令你搞不明白,且你若非要看这种教程,应当去发邮件问教程作者,其他GROMACS用户没法回答你。极其完整、详细、系统讲授分子动力学模拟知识和GROMACS使用的“北京科音分子动力学与GROMACS培训班”(http://www.keinsci.com/KGMX)里也从未牵扯到这种脚本和文件,而且根本就不需要非得用。sobtop(http://sobereva.com/soft/Sobtop)产生配体的拓扑文件是最佳选择,又普适又方便又极为流行,本论坛里有大量讨论,在本论坛提问时也有人能回答。

我强烈不建议初学者follow某些教程用的一些偏门、绕弯路或是过时的做法去算常规问题(而且很多教程本身也说得不明不白),明明目前有好得多、主流得多的方法。对各种常见问题怎么最好、最顺利地模拟在上述培训里都有充分讲授。
北京科音自然科学研究中心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

1240

帖子

6

威望

2553

eV
积分
3913

Level 5 (御坂)

傻傻的木瓜

4#
发表于 Post on 2024-12-16 21:28:17 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2024-12-16 21:51 编辑
sobereva 发表于 2024-12-16 02:50
GROMACS标准的使用过程中既不牵扯到叫str的东西也不牵扯到叫prm的东西,你去搜GROMACS官方手册也搜不到这种 ...

社长,鉴于楼主这样试图用CGenFF却缺失文件的问题多发,我提议在博文http://sobereva.com/266里讲清楚什么文件才是生成拓扑绝对必要的
√546=23.36664289109

885

帖子

4

威望

2047

eV
积分
3012

Level 5 (御坂)

A Student

5#
发表于 Post on 2024-12-16 22:08:09 | 只看该作者 Only view this author
本帖最后由 student0618 于 2024-12-16 22:12 编辑

真要用charmm+cgenff 的话新版本cgenff网站可一键convert to gromacs,下载unzip 后要用的是xxx_gmx.top xxx_gmx.pdb 以及charmm36.ff文件夹内的xxx_ffbonded.itp,修改一下即可使用。

更简单就是用那个从cgenff下载的包里面版本老一点的chamm.ff 处理蛋白,再改一下xxx_gmx.top include 到体系的top吧。

不过这些都不如用Amber+GAFF (gmx 配合sobtop) 省事。我若非要特别原因用charmm都情愿直接用Amber/GAFF的。
敬仰一针见血的指责,厌倦别有用心的赞美。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2024-12-17 02:31:20 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2024-12-16 21:28
社长,鉴于楼主这样试图用CGenFF却缺失文件的问题多发,我提议在博文http://sobereva.com/266里讲清楚什 ...

没特殊情况一律甭用CGenFF就完了,一律用sobtop产生GAFF/GAFF2拓扑文件
至于什么时候也不该用sobtop产生拓扑文件,前两天我给sobtop主页上加了一个FAQ:http://sobereva.com/soft/Sobtop/#FAQ20
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 02:53 , Processed in 0.250474 second(s), 28 queries , Gzip On.

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