计算化学公社

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

[Gaussian/gview] Gaussian使用范围分离泛函时,关键词LC-wPBE和LC-PBEPBE的区别是什么?

[复制链接 Copy URL]

7

帖子

0

威望

255

eV
积分
262

Level 3 能力者

各位老师,大家好:
最近我在学习范围分离泛函相关的内容,并进行了一些尝试,碰到一些问题不太理解:
在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,这些表示的含义我还不理解,希望老师们能解惑。

5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

2#
发表于 Post on 2023-7-5 04:45:15 | 只看该作者 Only view this author
wPBE是指的HSE03交换部分不计短程HF交换项的那部分(为了短程能引入HF,在PBE基础上对泛函做了相应的修改,所以前面有w来注明)。在g09里,wPBEh特指HSE03的交换部分。在09年wPBE被做了些修改,Gaussian里结合LC后叫做LC-wHPBE。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

7

帖子

0

威望

255

eV
积分
262

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-7-5 17:19:40 | 只看该作者 Only view this author
sobereva 发表于 2023-7-5 04:45
wPBE是指的HSE03交换部分不计短程HF交换项的那部分(为了短程能引入HF,在PBE基础上对泛函做了相应的修改, ...

谢谢sob老师的回复。自己之前的理解有一些错误(潜意识地把wPBE的w与ω参数联系在一起了)。

89

帖子

0

威望

1711

eV
积分
1800

Level 5 (御坂)

4#
发表于 Post on 2024-5-2 15:58:06 | 只看该作者 Only view this author
本帖最后由 风起~ 于 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
复制代码

1万

帖子

0

威望

7396

eV
积分
18149

Level 6 (一方通行)

5#
发表于 Post on 2024-5-2 16:45:13 | 只看该作者 Only view this author
风起~ 发表于 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里怎么定义。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

89

帖子

0

威望

1711

eV
积分
1800

Level 5 (御坂)

6#
发表于 Post on 2024-5-2 20:05:50 | 只看该作者 Only view this author
本帖最后由 风起~ 于 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
复制代码

1万

帖子

0

威望

7396

eV
积分
18149

Level 6 (一方通行)

7#
发表于 Post on 2024-5-2 20:15:41 | 只看该作者 Only view this author
风起~ 发表于 2024-5-2 13:05
老师,关键是前面的PBE部分该怎么写呀,其他部分我对照写的,应该没问题
麻烦您了 谢谢~

什么叫“前面的PBE部分”?
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

89

帖子

0

威望

1711

eV
积分
1800

Level 5 (御坂)

8#
发表于 Post on 2024-5-2 20:52:48 | 只看该作者 Only view this author
本帖最后由 风起~ 于 2024-5-3 08:15 编辑
wzkchem5 发表于 2024-5-2 20:15
什么叫“前面的PBE部分”?

就是最开始的泛函部分 !后面的部分

89

帖子

0

威望

1711

eV
积分
1800

Level 5 (御坂)

9#
发表于 Post on 2024-5-3 08:31:09 | 只看该作者 Only view this author
wzkchem5 发表于 2024-5-2 20:15
什么叫“前面的PBE部分”?

可以直接写为LC-PBE么

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 21:39 , Processed in 0.217435 second(s), 20 queries , Gzip On.

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