关于PW91及PBE的选择问题

sobereva

sigma-jlu 发表于 2023-10-11 19:10
多谢解答。
在用cp2k计算固体过程中，今天碰到一个pbe设置的问题，比如某泛函设置模块：

X跟C的数值没有关联
SCALE_X是交换泛函计算的交换能要乘的系数。诸如你就是纯粹用PBE交换泛函算交换能，那么SCALE_X就是1.0
而如果你用50% PBE交换泛函结合50% HF交换项，自然SCALE_X就是0.5
SCALE_C类似。如果你还要结合MP2或RPA相关能，那么SCALE_C就不是1.0（即只用PBE相关泛函计算相关能）

sigma-jlu

sobereva 发表于 2023-10-11 23:17
X跟C的数值没有关联
SCALE_X是交换泛函计算的交换能要乘的系数。诸如你就是纯粹用PBE交换泛函算交换能， ...

多谢卢老师解答，不知道有没有哪里讲到scale_x与scale_c的系数怎么构成完整的GGA泛函的？我也同步在翻几本参考书。

sigma-jlu

sobereva 发表于 2023-10-11 23:17
X跟C的数值没有关联
SCALE_X是交换泛函计算的交换能要乘的系数。诸如你就是纯粹用PBE交换泛函算交换能， ...

SCALE_X与SCALE_C参数没有关联，那他们的选择与确定，有什么指导原则吗？


96年pbe的原文Phys. Rev. Lett. 77, 3865，给了一对κ=0.967、μ=0.235的经验参数建议，不过该文在PRL被Zhang-Wang comment指出κ=1.245更好，无论是哪个经验参数更好，如果这两个参数跟SCALE_X与SCALE_C之间有某种对应关系，是不是意味着，对于特定的系统与研究的性质，SCALE_X与SCALE_C本身应有一个最合适的值，即肯定不能瞎设，那么在尚未扫描参数进行耗机时人力之前，有没有一些已知的常识需要注意的？

“since the whole Exc is usually developed and parameterized, there is some error compensation between both partners and combination of Ex and Ec functionals issuing from different sources could be dangerous.” J. Chem. Phys. 110, 6158 (1999)

卡开发发

sigma-jlu 发表于 2023-10-11 19:10
多谢解答。
在用cp2k计算固体过程中，今天碰到一个pbe设置的问题，比如某泛函设置模块：

如果是纯PBE的话其实都按1就行，甚至可以不写这些。SCALE_X和SCALE_C其实就是PBE的交换和关联部分的比例，一般情况杂化泛函的情况才需要修改这个比例。例如PBE0的情况是0.75*PBEx+0.25EXX+1.0PBEc，那么PBE下面的SCALE_X就按照0.75，HF的部分的FRACTION则按0.25设置。这里的例子PBE_HOLE_T_C_LR针对的是PBE0_TC_LRC方法的参数。

sobereva

sigma-jlu 发表于 2023-10-12 08:27
SCALE_X与SCALE_C参数没有关联，那他们的选择与确定，有什么指导原则吗？

你用别人现成的泛函的时候，泛函原文怎么定义的就怎么写就完了
你提到的诸如κ这样的纯泛函部分的内在的一些细节参数不是通过这里体现的，是纯泛函部分的代码里体现的（改源代码层面）。你把问题想复杂了。