计算化学公社

标题: Gaussian使用范围分离泛函时，关键词LC-wPBE和LC-PBEPBE的区别是什么？ [打印本页]
作者
Author: symbolhone    时间: 2023-7-4 10:51
标题: Gaussian使用范围分离泛函时，关键词LC-wPBE和LC-PBEPBE的区别是什么？
各位老师，大家好：
最近我在学习范围分离泛函相关的内容，并进行了一些尝试，碰到一些问题不太理解：
在sob老师的这篇文章中：在Gaussian中自定义范围分离泛函的方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)
里面有解读过LC-PBE，我也进行了尝试，输入关键词LC-PBEPBE，计算时泛函的参数与sob老师讲到的一致：
IExCor=11009 DFT=T Ex=LC-PBE+HF Corr=PBE ExCW=0 ScaHFX=  1.000000
ScaDFX=  1.000000  1.000000  1.000000  1.000000 ScalE2=  1.000000  1.000000
IRadAn=      5 IRanWt=     -1 IRanGd=            0 ICorTp=0 IEmpDi=  4
HFx  wShort=  0.000000 wLong=  0.470000 cFull=  0.000000 cShort=  0.000000 cLong=  1.000000
DFx  wShort=  0.000000 wLong=  0.470000 cFull=  0.000000 cShort=  0.000000 cLong=  1.000000

另外，我尝试了LC-wPBE的关键词，参数为：
IExCor=32609 DFT=T Ex+Corr=LC-wPBE ExCW=0 ScaHFX=  1.000000
ScaDFX=  1.000000  0.000000  1.000000  1.000000 ScalE2=  1.000000  1.000000
IRadAn=      5 IRanWt=     -1 IRanGd=            0 ICorTp=0 IEmpDi=  4
HFx  wShort=  0.000000 wLong=  0.400000 cFull=  0.000000 cShort=  0.000000 cLong=  1.000000
DFx  wShort=  0.000000 wLong=  0.400000 cFull=  0.000000 cShort=  0.000000 cLong=  1.000000

其中ScaDFX的第二项为0，即c=0，也就是GGA交换项对LDA交换项的梯度修正部分系数为0，显然二者是不同的。
我尝试的第三个关键词LC-whPBE，参数为：
IExCor=33909 DFT=T Ex+Corr=LC-wHPBE ExCW=0 ScaHFX=  1.000000
ScaDFX=  1.000000  1.000000  1.000000  1.000000 ScalE2=  1.000000  1.000000
IRadAn=      5 IRanWt=     -1 IRanGd=            0 ICorTp=0 IEmpDi=  4
HFx  wShort=  0.000000 wLong=  0.400000 cFull=  0.000000 cShort=  0.000000 cLong=  1.000000
DFx  wShort=  0.000000 wLong=  0.400000 cFull=  0.000000 cShort=  1.000000 cLong=  0.000000

回看PBE泛函，其b=0 c=d=1 e=f=1，结合LC后，先定义为全局100% HF交换+100% DFT交换的状态，之后再定义不同范围如何杂化，则bcdef均变为1，
那么LC-wPBE明显是在此基础上又进行了修改的，因为它的c为0。


因此我有以下问题想请教各位老师：
1、LC-wPBE和LC-PBEPBE二者的区别是什么，将c定义为0是基于什么原因或者他们的应用有什么不同。
2、我阅读一些文献中提到的LRC-wPBEh，LC-whPBE，比如：“The Journal of chemical physics. 2009, 130(5): 54112.”那是否是我尝试的第三个关键词LC-whPBE？
如果不是的话，那么LC-whPBE关键词运行的计算又代表了什么含义呢？
3、对于泛函的参数： IExCor=32609 DFT=T Ex+Corr=LC-wPBE ExCW=0。
其中IExCor、Ex+Corr=LC-wPBE、Ex=LC-PBE+HF、Corr=PBE、ExCW=0，这些表示的含义我还不理解，希望老师们能解惑。
作者
Author: sobereva    时间: 2023-7-5 04:45
wPBE是指的HSE03交换部分不计短程HF交换项的那部分（为了短程能引入HF，在PBE基础上对泛函做了相应的修改，所以前面有w来注明）。在g09里，wPBEh特指HSE03的交换部分。在09年wPBE被做了些修改，Gaussian里结合LC后叫做LC-wHPBE。
作者
Author: symbolhone    时间: 2023-7-5 17:19
sobereva 发表于 2023-7-5 04:45
wPBE是指的HSE03交换部分不计短程HF交换项的那部分（为了短程能引入HF，在PBE基础上对泛函做了相应的修改， ...

谢谢sob老师的回复。自己之前的理解有一些错误（潜意识地把wPBE的w与ω参数联系在一起了）。
作者
Author: 风起~    时间: 2024-5-2 15:58
本帖最后由 风起~ 于 2024-5-2 16:22 编辑

我查到CAM-B3LYP在orca中为如下构造，可是不知道如何自定义LC-whPBE呀，还请老师指教@wzkchem5   
  1. ! RKS CAM-B3LYP DEF2-SVP
  2. # default parameters for CAM-B3LYP:
  3. %method
  4. RangeSepEXX True # must be set
  5. RangeSepMu 0.33 # should not be set to 0.0 or below
  6. RangeSepScal 0.46 # should sum to 1 with ACM-A and ACM-B
  7. ACM 0.19, 0.35, 0.81 # ACM-A, ACM-B, ACM-C(same as B3LYP)
  8. end
复制代码
作者
Author: wzkchem5    时间: 2024-5-2 16:45
风起~ 发表于 2024-5-2 08:58
我查到CAM-B3LYP在orca中为如下构造，可是不知道如何自定义LC-whPBE呀，还请老师指教@wzkchem5

一楼（以及LC-whPBE的原文）已经给了定义需要的参数了，填进去即可。
如果不知道哪个参数该填哪，就去看CAM-B3LYP的原文。你看了CAM-B3LYP的原文，看了CAM-B3LYP在orca里怎么定义，又看了LC-whPBE的原文，必然知道LC-whPBE在orca里怎么定义。
作者
Author: 风起~    时间: 2024-5-2 20:05
本帖最后由 风起~ 于 2024-5-2 20:07 编辑
wzkchem5 发表于 2024-5-2 16:45
一楼（以及LC-whPBE的原文）已经给了定义需要的参数了，填进去即可。
如果不知道哪个参数该填哪，就去看 ...

老师，关键是前面的PBE部分该怎么写呀，其他部分我对照写的，应该没问题
麻烦您了 谢谢~
  1. ! RKS PBE DEF2-SVP
  2. %method
  3. RangeSepEXX True # must be set
  4. RangeSepMu 0.4 # should not be set to 0.0 or below
  5. RangeSepScal 1.0 # should sum to 1 with ACM-A and ACM-B
  6. ACM 0.0, 0.0, 1.0 # ACM-A, ACM-B, ACM-C(same as B3LYP)
  7. end
复制代码


CAM-B3LYP
  1. IExCor=20419 DFT=T Ex+Corr=CAM-B3LYP ExCW=0 ScaHFX=  1.000000
  2. ScaDFX=  1.000000  1.000000  1.000000  0.810000 ScalE2=  1.000000  1.000000
  3. IRadAn=      5 IRanWt=     -1 IRanGd=            0 ICorTp=0 IEmpDi=141
  4. HFx  wShort=  0.000000 wLong=  0.330000 cFull=  0.190000 cShort=  0.000000 cLong=  0.460000
  5. DFx  wShort=  0.000000 wLong=  0.330000 cFull=  0.190000 cShort=  0.000000 cLong=  0.460000
复制代码
LC-wHPBE
  1. IExCor=33909 DFT=T Ex+Corr=LC-wHPBE ExCW=0 ScaHFX=  1.000000
  2. ScaDFX=  1.000000  1.000000  1.000000  1.000000 ScalE2=  1.000000  1.000000
  3. IRadAn=      5 IRanWt=     -1 IRanGd=            0 ICorTp=0 IEmpDi=  4
  4. HFx  wShort=  0.000000 wLong=  0.400000 cFull=  0.000000 cShort=  0.000000 cLong=  1.000000
  5. DFx  wShort=  0.000000 wLong=  0.400000 cFull=  0.000000 cShort=  1.000000 cLong=  0.000000
复制代码
作者
Author: wzkchem5    时间: 2024-5-2 20:15
风起~ 发表于 2024-5-2 13:05
老师，关键是前面的PBE部分该怎么写呀，其他部分我对照写的，应该没问题
麻烦您了 谢谢~

什么叫“前面的PBE部分”？
作者
Author: 风起~    时间: 2024-5-2 20:52
本帖最后由 风起~ 于 2024-5-3 08:15 编辑
wzkchem5 发表于 2024-5-2 20:15
什么叫“前面的PBE部分”？

就是最开始的泛函部分 !后面的部分
作者
Author: 风起~    时间: 2024-5-3 08:31
wzkchem5 发表于 2024-5-2 20:15
什么叫“前面的PBE部分”？

可以直接写为LC-PBE么
作者
Author: lycheeho    时间: 2026-1-14 21:14
sobereva 发表于 2023-7-5 04:45
wPBE是指的HSE03交换部分不计短程HF交换项的那部分（为了短程能引入HF，在PBE基础上对泛函做了相应的修改， ...

你好Sob老师，请问LC-wHPBE泛函可以计算找到分子开壳层单重态基态吗？我用guess=mix计算后，看log文件里面没有S2项信息，且轨道是闭壳层的。谢谢！
作者
Author: sobereva    时间: 2026-1-14 23:59
lycheeho 发表于 2026-1-14 21:14
你好Sob老师，请问LC-wHPBE泛函可以计算找到分子开壳层单重态基态吗？我用guess=mix计算后，看log文件里 ...

怎么算跟用什么泛函没直接关系，做法下文该说的都说了。虽然HF成份高的泛函相对更容易得到对称破缺态
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82
作者
Author: lycheeho    时间: 2026-1-15 10:20
sobereva 发表于 2026-1-14 23:59
怎么算跟用什么泛函没直接关系，做法下文该说的都说了。虽然HF成份高的泛函相对更容易得到对称破缺态
谈 ...

好的，谢谢sob老师耐心回复！



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3