计算化学公社

标题: 关于Gaussian-DFT自定义泛函公式的参数 [打印本页]

作者
Author: brothers 时间: 2015-10-20 10:45
标题: 关于Gaussian-DFT自定义泛函公式的参数
本帖最后由 brothers 于 2015-10-20 20:05 编辑

研究中要用到自定义泛函，所以参考了gaussian09的帮助文档以及A chemist's guide to DFT 那本书，作为对比，看了B3LYP的定义。
Gaussian09中是如此定义的
(, 下载次数 Times of downloads: 29)
书中是如此定义的
(, 下载次数 Times of downloads: 33)
问题如下：(因为自己不是特别懂，所以将理解的东西尽量写下来，有错误的话请各位不吝赐教)
关于交换部分：
Gaussian09的定义中，P2部分是non-local的exact HF exchange，P1*P4是local的Slater交换。杂化的目的，就是通过引入non-local的exchange来修正local的exchange对hole的描述。那么，要求P2+P1*P4=1.0？假若我调整HF比例，是否要满足这一限制？
书中B88的交换泛函直接用的E_B88，为什么Gaussian却用了个deltaE来表示？deltaE与E等价吗？（看有的书中，E=E_lsda+deltaE_B88）

关于相关部分：
看书里的定义，local和nonlocal的E_C，其系数加和为1.0，可是Gaussian中的定义看起来却没有这种限制（P5=0.81，P6=1.0），这是怎么回事呢？

这些关系实在是把我看糊涂了，local啊，non-local啊，感觉也是模模糊糊。请大神指点一下迷津，同时推荐一些相关资料我再读读看~~万分感激~~

作者
Author: sobereva 时间: 2015-10-21 03:50
local就是指LSDA，这类泛函一般都是用SVWN，即指交换部分是Slater(S)，相关部分是VWN。
non-local就是指LDA以外的那些东西，比如GGA、HF杂化部分。

注意无论对于交换还是相关部分，都是E_nonlocal=ΔE_nonlocal+E_local。因此E_B88=E_VWN+ΔE_B88，E_LYP=E_VWN+ΔE_LYP。

交换部分，P2+P1*P4必须等于1.0，调整的时候必须满足此限制。由于一般都是让P1=1.0，所以只需要P2+P4=1.0就够了。

相关部分，P6必须要求是1.0。你也可以这么写，就容易和6-28式对照了：
P6*E_VWN+P5*ΔE_LYP
=P6*E_VWN+P5*(E_LYP-E_VWN)
=(P6-P5)*E_VWN+P5*E_LYP
可见，P5相当于6-28中的c，是可调的。

作者
Author: brothers 时间: 2015-10-21 08:43

sobereva 发表于 2015-10-21 03:50
local就是指LSDA，这类泛函一般都是用SVWN，即指交换部分是Slater(S)，相关部分是VWN。
non-local就是指LD ...

非常感谢！我会继续研究研究相关知识，不能白费这么精彩的回复

欢迎光临计算化学公社 (http://bbs.keinsci.com/)