计算化学公社

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

[新手求助] 如何获取大分子体系的非溶剂环境下原子电荷?

[复制链接 Copy URL]

5

帖子

0

威望

73

eV
积分
78

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
希望探究蛋白小分子体系非溶剂环境下原子电荷与溶剂环境下原子电荷有多大区别。
通过sob老师的培训班了解到,AMBER03力场之后蛋白质体系的原子电荷是结合隐式溶剂模型用RESP方法拟合的(作为非极化力场,预先考虑了溶剂环境对原子电荷的影响),不是真空环境下的原子电荷。
对于2000+个原子的体系,使用Gaussian做真空下结构优化计算单点能得到波函数用Multiwfn做RESP电荷拟合似乎计算量较大。
这样的情况下,该如何获取该体系非溶剂环境下的原子电荷呢?

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

6#
发表于 Post on 2025-12-15 23:45:09 | 只看该作者 Only view this author
ZDcccA 发表于 2025-12-15 10:42
谢谢社长耐心解答!想再追问一下,如果用AMOEBA这样的可极化力场在不加水盒子的情况下赋予蛋白质体系原子 ...

看不懂你的意思
可极化力场完全不是用来算原子电荷的
北京科音自然科学研究中心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

5

帖子

0

威望

73

eV
积分
78

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2025-12-15 10:42:31 | 只看该作者 Only view this author

谢谢社长耐心解答!想再追问一下,如果用AMOEBA这样的可极化力场在不加水盒子的情况下赋予蛋白质体系原子电荷是不是也可以起到一样的效果,即得到气相条件下蛋白质原子电荷?

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

4#
发表于 Post on 2025-12-14 21:38:42 | 只看该作者 Only view this author
ZDcccA 发表于 2025-12-14 11:59
谢谢sob老师!是气相蛋白质体系的原子电荷,分子力场的蛋白质原子电荷似乎是考虑了溶剂极化效应的(至少A ...

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

5

帖子

0

威望

73

eV
积分
78

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-12-14 11:59:11 | 只看该作者 Only view this author
本帖最后由 ZDcccA 于 2025-12-14 12:03 编辑
sobereva 发表于 2025-12-13 21:37
说清楚具体是什么体系
如果是蛋白质,其整体的原子电荷是分子力场定义的各个残基的原子电荷拼接得到的,根 ...

谢谢sob老师!是气相蛋白质体系的原子电荷,分子力场的蛋白质原子电荷似乎是考虑了溶剂极化效应的(至少AMBER03及之后是)。如果不是要极化后的蛋白质原子电荷而是要气相条件下的,是需要自己重新算每种氨基酸的气相的RESP电荷吗?谢谢指点!

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

2#
发表于 Post on 2025-12-13 21:37:13 | 只看该作者 Only view this author
说清楚具体是什么体系
如果是蛋白质,其整体的原子电荷是分子力场定义的各个残基的原子电荷拼接得到的,根本不需要自己算。只有存在非标准残基的情况下才需要自己算非标准残基部分的原子电荷,做法下文说了
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441http://bbs.keinsci.com/thread-10880-1-1.html
北京科音自然科学研究中心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

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

GMT+8, 2026-1-23 21:45 , Processed in 0.165131 second(s), 22 queries , Gzip On.

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