计算化学公社

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

[综合交流] 不用色散校正泛函计算

[复制链接 Copy URL]

125

帖子

1

威望

1047

eV
积分
1192

Level 4 (黑子)

本帖最后由 fankui1990129 于 2015-6-14 22:26 编辑

(如果嫌字多就看加粗的部分就行了)最近计算一个体系,用的是B3LYP,最近仔细看了sob老师的帖子,发现这个方法一开始就错了。后来打算用B3LYP-D3,但是这个貌似这个只有G09才开始有。我们学校买的不是09,换言之我就算得到想要的结果也不能用。


我要算的是如图所示的吉布斯自由能。我原先搭建的时候是将原料一起搭建,产物一起搭建。因为不能用B3LYP-D3方法,我想能不能将A、B、C和水用B3LYP分别计算吉布斯自由能。因为我的猜想化学实验已经证明了,量化只是为了增加文章的内容和说服力。所以,请问我的方法可以么?


顺便问问sob老师发过介绍色散校正泛函的帖子么?如果有的话请告诉我一声,我自己去找就好了~~ (现在特别期待sob老师写的书




untitled.gif (1.46 KB, 下载次数 Times of downloads: 50)

反应方程

反应方程

6万

帖子

99

威望

5万

eV
积分
120194

管理员

公社社长

2#
发表于 Post on 2015-6-15 01:55:13 | 只看该作者 Only view this author
好多帖子都涉及色散校正
乱谈DFT-D
http://sobereva.com/83
DFT-D色散校正的使用
http://sobereva.com/210
大体系弱相互作用计算的解决之道
http://sobereva.com/214

分子间相互作用计算没必要非得用B3LYP-D3,用D.01之前版本就支持的M06-2X也能达到很好效果。


我并不清楚A和B的复合物具体是什么相互作用维持的,如果是静电为主,那么B3LYP基本没问题。

至于是每个分子单独计算吉布斯自由能,还是按照复合物来计算吉布斯自由能,这看情况而定,哪个都不能说是错的,意义不一样而已。如果是每个分子单独来计算,B3LYP就更没问题了。
北京科音自然科学研究中心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

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

3#
发表于 Post on 2015-6-15 09:58:37 | 只看该作者 Only view this author
DFT-D是个好东东,简单高效。
另外你可以试试gamess-US-64,可以用b3lyp-d。
当然,首选sob老师的建议
十室之邑,必有忠信如丘焉,不如丘之好学也。

3754

帖子

3

威望

1万

eV
积分
19681

Level 6 (一方通行)

围观吃瓜群众

4#
发表于 Post on 2015-6-15 12:08:47 | 只看该作者 Only view this author
我持一些不同的态度(能力有限仅作为一些讨论):

先做些假设,如对于一个确定的具有两个片段的体系,采用的是GD2或GD3这些经验校正方式:
(1)在真空中的相互作用在DFT-GD2或GD3下相互作用可以分为两个部分,一部分通过SCF得到(SCF中包含一部分色散,但可能不合理),另一部分则来自于色散校正;
(2)加上了溶剂化,DFT-D仍然分为SCF和-D校正得到的两部分能量,因为考虑了溶剂化,SCF的部分包含的色散作用会有一定的变化,溶质和溶剂的色散作用原则上被包含在溶剂化模型中,但需要注意的是-D的部分与电子态无关,这个时候-D校正的能量与真空下计算的结果无差别。

这时候问题就来了,因为分子在溶剂中与溶剂发生相互作用原则上极化率会发生变化,简单根据London公式进行估计色散作用也应发生变化,如果采用DFT-D原有的参数可能无法体现这种趋势(从理论上说可能得调整C6这个参数)。可能采用极化率自洽的色散修正方式(如TS)或是能够较好描述弱作用的泛函如M06-2X再结合溶剂化,这样结果更加合理。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

41

帖子

0

威望

339

eV
积分
380

Level 3 能力者

5#
发表于 Post on 2016-11-18 17:23:33 | 只看该作者 Only view this author
sobereva 发表于 2015-6-15 01:55
好多帖子都涉及色散校正
乱谈DFT-D
http://sobereva.com/83

Sob老师,您好。我看到您的帖子《DFT-D色散校正的使用》中建议使用BLYP-D3(BJ)方法计算柔性大体系构象、大体系之间相互作用。  有没有相关文献证明这个方法应用在大体系中的精度仍可接受。谢谢老师!

6万

帖子

99

威望

5万

eV
积分
120194

管理员

公社社长

6#
发表于 Post on 2016-11-18 18:46:52 | 只看该作者 Only view this author
dy19930508 发表于 2016-11-18 17:23
Sob老师,您好。我看到您的帖子《DFT-D色散校正的使用》中建议使用BLYP-D3(BJ)方法计算柔性大体系构象、 ...

D3的可靠性早就被充分检验了,应用例子很多,诸如DOI: 10.1039/c4cc06722c
北京科音自然科学研究中心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, 2025-8-18 03:01 , Processed in 0.201543 second(s), 23 queries , Gzip On.

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