计算化学公社

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

[GROMACS] 木质素与重金属离子作用体系的力场选择问题

[复制链接 Copy URL]

8

帖子

0

威望

433

eV
积分
441

Level 3 能力者

本帖最后由 oorangew 于 2025-7-15 10:49 编辑

各位老师好,我想做的是木质素分子和重金属离子(如铜离子、铅离子)的相互作用分析,后续还需要在基础的木质素分子上进行改性,如图所示,如加入氨基等。查阅了近段时间以来有关于木质素模拟的文章,一般用的都是Josh Vermaas的木质素CHARMM参数:DOI: 10.1039/c8gc03209b。问了作者,说不建议将该参数用于进行后续的改性木质素,因为原参数没有相关基团的参数。
木质素的结构单体与键连接方式比较复杂。所以想问一下各位老师,建议选择什么力场并如何生成后续GROMACS使用的拓扑文件呢?
木质素分子结构示例:(示例分子有近400个原子)

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

2#
发表于 Post on 2025-5-20 06:12:02 | 只看该作者 Only view this author
用sobtop(http://sobereva.com/soft/Sobtop)就可以,基于GAFF或GAFF2力场
北京科音自然科学研究中心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

8

帖子

0

威望

433

eV
积分
441

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-5-20 09:29:25 | 只看该作者 Only view this author
sobereva 发表于 2025-5-20 06:12
用sobtop(http://sobereva.com/soft/Sobtop)就可以,基于GAFF或GAFF2力场

感谢sob老师的回复。这样的话审稿人会不会comment为什么用这个力场呢?

887

帖子

3

威望

1670

eV
积分
2617

Level 5 (御坂)

傻傻的木瓜

4#
发表于 Post on 2025-5-20 09:55:26 | 只看该作者 Only view this author
sobereva 发表于 2025-5-20 06:12
用sobtop(http://sobereva.com/soft/Sobtop)就可以,基于GAFF或GAFF2力场

社长,从结构上看“木质素分子和重金属离子”的相互作用应该涉及水的氧和酚的氧发生配体交换的反应过程吧,这种配位与解离的成键断键过程能用经典分子力场的MD模拟来描述么?
√546=23.36664289109

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

5#
发表于 Post on 2025-5-21 04:05:51 | 只看该作者 Only view this author
oorangew 发表于 2025-5-20 09:29
感谢sob老师的回复。这样的话审稿人会不会comment为什么用这个力场呢?

这是极为主流、稳健的用于有机体系的力场,没什么可comment之处
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

6#
发表于 Post on 2025-5-21 04:07:09 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-5-20 09:55
社长,从结构上看“木质素分子和重金属离子”的相互作用应该涉及水的氧和酚的氧发生配体交换的反应过程吧 ...

除非是反应力场,否则不能
北京科音自然科学研究中心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

8

帖子

0

威望

433

eV
积分
441

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2025-5-22 11:25:59 | 只看该作者 Only view this author
sobereva 发表于 2025-5-21 04:05
这是极为主流、稳健的用于有机体系的力场,没什么可comment之处

好的,谢谢老师!我还有个疑问,木质素的这种结构(图示分子结构有400个原子)用sobtop产生top文件的话可以指认MMFF94力场电荷吗,还是建议计算RESP电荷呢?老师另外想让我做的分子有1000+原子,木质素的结构单元和键连接方式的复杂原因所以整个聚合物的重复的程度较小,所以如果像sobtop主页中提及的套用模型体系来计算各个单元的RESP电荷貌似不太理想。

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

8#
发表于 Post on 2025-5-22 15:02:14 | 只看该作者 Only view this author
oorangew 发表于 2025-5-22 11:25
好的,谢谢老师!我还有个疑问,木质素的这种结构(图示分子结构有400个原子)用sobtop产生top文件的话可 ...

MMFF94可以勉强凑合。也可以尝试用https://molmodel.com/hg/python_DENR_charges/算电荷(我并未实际用过,只是审JCIM稿子时注意到,有问题需咨询作者),比MMFF94的静电势重现性好一些,而且也不依赖于量化计算,但质量明显不及RESP

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
oorangew + 5 非常感谢sob大大!

查看全部评分 View all ratings

北京科音自然科学研究中心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

8

帖子

0

威望

433

eV
积分
441

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2025-5-22 21:59:46 | 只看该作者 Only view this author
sobereva 发表于 2025-5-22 15:02
MMFF94可以勉强凑合。也可以尝试用https://molmodel.com/hg/python_DENR_charges/算电荷(我并未实际用过 ...

好的,感谢sob老师的耐心解答,我先试试用MMFF94

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 04:12 , Processed in 0.176534 second(s), 24 queries , Gzip On.

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