高斯中可调控的LC-PBE的α和β参数该如何设定，在输入文件中该怎么写?

wenxiao

各位老师好，请教一个问题，高斯中可调控的LC-PBE的α和β参数该如何设定，在输入文件中该怎么写，不胜感激！

xp47

具体用法可以在IOP手册中寻找：

比如：
C2IOp(76)
C     IOp(76) ... Mixing of HF and DFT.  Negative values correspond
C           0 ... Default: pure HF, DFT or mixed in accord with IOp(76)
C                 implied by IOp(74) or set by IOp(77).
...

C          -8 ... B97-1 coefficients.
C          -7 ... HCTH coefficints.
C          -6 ... B98 coefficients.
C          -5 ... mPW91PW91 coefficients.
C          -4 ... Becke 3 coefficients:  aLSD + (1-a)HF + b(dBx) +
C                 VWN + c(LYP-VWN), with a=0.8 b=0.72 c=0.81
C                 Note that Becke actually used Perdew correlation
C                 rather than LYP.
C          -3 ... Becke "Half and Half" 0.5 HF + 0.5 Xc + Corr
C          -2 ... Coefficients of 0 and 0 (no exchange).
C          -1 ... Coefficients of 0.0 and 1.0 for DFT and HF,
C                 respectively.
C           0 ... Default: pure HF, DFT or mixed in accord with IOp(76)
C  MMMMMNNNNN ... Mixture of MMMMM/10000 DFT exchange and NNNNN/10000
C                 HF exchange.  The DFT exchange factor multiplies any
C                 implied by IOp(74) or set by IOp(77).
C

C2IOp(77)
C     IOp(77) ... Mixing of local and non-local exchange:
C                 by IOp(76)
C2IOp(78)
C     IOp(78) ... Mixing of local and non-local correlation:
C                 by IOp(76)
此外 对于带Omega的杂化泛函可以用以下关键词调控
C     IOp(107) ... Omega for short/long range Hartree-Fock exchange
C     IOp(108) ... Omega for short/long range DFT exchange
C     IOp(109) ... Omega for short/long range DFT correlation
C     IOp(126) ... Omega for short/long range 1/r operator in E(2,AA)
C     IOp(127) ... Omega for short/long range 1/r operator in E(2,AB)

wenxiao

xp47 发表于 2020-3-6 14:02
具体用法可以在IOP手册中寻找：

比如：

感谢回复，如果设置omega=0.146, α=0.1, β=0.9是这样写吗：LC-PBEPBE IOP(3/107=0146000000,3/108=0146000000，3/119=1000001000,3/120=0000009000)？谢谢！

xp47

wenxiao 发表于 2020-3-6 20:22
感谢回复，如果设置omega=0.146, α=0.1, β=0.9是这样写吗：LC-PBEPBE IOP(3/107=0146000000,3/108=0146 ...

应该是可以的（我之前是直接改重新编译 你这么一说才想起来还有IOp

wenxiao

xp47 发表于 2020-3-7 16:33
应该是可以的（我之前是直接改重新编译 你这么一说才想起来还有IOp

我按这样的IOP设置计算的结果和文献（J. Chem. Phys. 147, 044301 (2017)）不符，可能这样设置是不正确的。您说的直接改重新编译具体怎样操作呢？谢谢

xp47

wenxiao 发表于 2020-3-8 22:10
我按这样的IOP设置计算的结果和文献（J. Chem. Phys. 147, 044301 (2017)）不符，可能这样设置是不正确的 ...

我建议还是把5个omega 都设置下，如高斯16c的现有定义omega如下(我觉着你得改下这里)：
        Omega(2,IHFX) = Pt40
        Omega(3,IHFX) = Zero
        Omega(5,IHFX) = One
        Omega(2,IDFX) = Pt40
        Omega(3,IDFX) = Zero                 
        Omega(5,IDFX) = One


（DFX=DFT_x HFX=HF_x）

C     Omega(1,...) : w for short-range;
C     Omega(2,...) : w for long-range;
C     Omega(3,...) : coefficient of full range;
C     Omega(4,...) : coefficient of short range;
C     Omega(5,...) : coefficient of long range;
比如IOp130 131


重新编译的话就是将原来GETDFT函数下的参数手动改掉然后编译下.

wenxiao

xp47 发表于 2020-3-9 10:02
我建议还是把5个omega 都设置下，如高斯16c的现有定义omega如下(我觉着你得改下这里)：
        Omega(2 ...

请问您有具体的例子可以参考下吗？谢谢

xp47

wenxiao 发表于 2020-3-9 13:55
请问您有具体的例子可以参考下吗？谢谢

有的 https://www.researchgate.net/pos ... M-b3lyp_in_Gaussian 参考下CAM的