计算化学公社

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

[新手求助] 用oniom计算重组能(考虑外重组能)的情况

[复制链接 Copy URL]

71

帖子

0

威望

379

eV
积分
450

Level 3 能力者

我想用绝热势能面法(四点法)计算重组能,考虑周围分子对其的影响(相当于考虑外重组能),以晶体结构为基础,选取中心分子用QM处理,周围分子用MM处理,计算0电荷下优化态的能量及以上述优化结构为基础给中心分子一个正电荷的单点能,Gaussian路径部分如下:


任务可以正常结束,但计算的HF能量值为0.34左右,是关键字部分设定出错 了吗?而且用四点法的话计算重组能的值为负值,该如何计算呢?谢谢

6万

帖子

99

威望

5万

eV
积分
120201

管理员

公社社长

2#
发表于 Post on 2016-5-5 10:26:55 | 只看该作者 Only view this author
具体用的什么计算公式?如果直接能量相减,内重组能肯定也会掺进去,已经将之扣除了么?
如果你觉得0.34很小而异常,分子力场的能量和从头算给出的不是一个数量级的,因为参考能量不同,分子力场是以所有力场参数都在平衡值为能量0点的。而且还要注意结果取没取错,是取了UFF算的整个体系的还是取了算model体系的。
北京科音自然科学研究中心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

71

帖子

0

威望

379

eV
积分
450

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2016-5-5 15:04:10 | 只看该作者 Only view this author
sobereva 发表于 2016-5-5 10:26
具体用的什么计算公式?如果直接能量相减,内重组能肯定也会掺进去,已经将之扣除了么?
如果你觉得0.34很 ...

公式1是计算重组能的公式,我的理解是如果按照oniom分层,不固定任何分子,计算上述4点按公式1得到的值就是总重组能,如果只计算外重组能,则是在0,1电荷下优化整个体系,在计算单点能时,固定中心QM处理的分子,然后这4点按公式1得到的是外重组能,这样的理解正确吗?尤其是绿色线标记部分,我不太理解是什么意思,求解答!还有一个问题是,计算ONIOM任务后,能量值用那个值更合理呢?





6万

帖子

99

威望

5万

eV
积分
120201

管理员

公社社长

4#
发表于 Post on 2016-6-1 09:37:20 | 只看该作者 Only view this author
lp2014 发表于 2016-5-5 15:04
公式1是计算重组能的公式,我的理解是如果按照oniom分层,不固定任何分子,计算上述4点按公式1得到的值就 ...


假设C是+1态,基本流程
a 在中性下优化整体
b 把中心分子电荷设成1,用a的结构计算单点能
c 把中心分子电荷设成1,用a的结构当初始结构优化整体,优化的时候将中心分子坐标冻住
d 在中性下,用c的结构计算单点能

外重组能=(1/2)*(b-c+d-a)

具体怎么操作看你用什么理论级别。诸如分子力场,中心分子电荷是0还是1,就在相应电荷设定下计算拟合静电势电荷,给中心分子赋上。ONIOM(QM:MM)貌似无法直接实现,因为Gaussian貌似不允许冻住QM所有原子,除非自己改代码,或者问客服有什么奇淫技巧实现。

ONIOM能量取最后extrapolated energy
北京科音自然科学研究中心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

71

帖子

0

威望

379

eV
积分
450

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2016-6-7 11:01:36 | 只看该作者 Only view this author
Gaussian中ONIOM(QM:MM)的确不允许冻住QM所有原子,自己改代码能力有限,问题只能暂时搁置了,还是 谢谢细心解答!

6万

帖子

99

威望

5万

eV
积分
120201

管理员

公社社长

6#
发表于 Post on 2016-6-7 12:31:07 | 只看该作者 Only view this author
lp2014 发表于 2016-6-7 11:01
Gaussian中ONIOM(QM:MM)的确不允许冻住QM所有原子,自己改代码能力有限,问题只能暂时搁置了,还是 谢谢 ...

而且直接用高斯QM:MM也不准确,哪怕是电子嵌入,一方面UFF力场本身不很准,和MM3、MMFF94差距不小,而且MM原子电荷没法与QM区域电子结构达到极化效果的自洽,最好类似这篇文章加入一个迭代过程:J. AM. CHEM. SOC. 2008, 130, 12377–12384
北京科音自然科学研究中心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

204

帖子

0

威望

2717

eV
积分
2921

Level 5 (御坂)

7#
发表于 Post on 2017-6-5 23:35:32 | 只看该作者 Only view this author
sobereva 发表于 2016-6-7 12:31
而且直接用高斯QM:MM也不准确,哪怕是电子嵌入,一方面UFF力场本身不很准,和MM3、MMFF94差距不小,而且M ...

搭车顺问,对于文献中的超晶胞的选择有什么习惯或原则
比如



文献JACS, 2008, 130, 12377中也是这样



根本没有解释问啥这样选择超晶胞,求帮助,谢谢

6万

帖子

99

威望

5万

eV
积分
120201

管理员

公社社长

8#
发表于 Post on 2017-6-6 08:54:49 | 只看该作者 Only view this author
赵云跳槽 发表于 2017-6-5 23:35
搭车顺问,对于文献中的超晶胞的选择有什么习惯或原则
比如


1 不能让体系原子数太多,否则算不动
2 必须足够大,使得边界效应对结果影响小。这和具体问题,体系结构有关。凭经验拿不准的话可以对结果随晶胞平移复制次数进行收敛性测试
北京科音自然科学研究中心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

204

帖子

0

威望

2717

eV
积分
2921

Level 5 (御坂)

9#
发表于 Post on 2017-6-6 09:44:54 | 只看该作者 Only view this author
谢谢sob老师

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 05:09 , Processed in 0.564137 second(s), 23 queries , Gzip On.

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