计算化学公社

标题: 高斯中可调控的LC-PBE的α和β参数该如何设定，在输入文件中该怎么写? [打印本页]

作者
Author: wenxiao 时间: 2020-3-6 00:07
标题: 高斯中可调控的LC-PBE的α和β参数该如何设定，在输入文件中该怎么写?
各位老师好，请教一个问题，高斯中可调控的LC-PBE的α和β参数该如何设定，在输入文件中该怎么写，不胜感激！

作者
Author: xp47 时间: 2020-3-6 14:02
本帖最后由 xp47 于 2020-3-6 14:03 编辑

具体用法可以在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)

作者
Author: wenxiao 时间: 2020-3-6 20:22

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)？谢谢！

作者
Author: xp47 时间: 2020-3-7 16:33

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

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

作者
Author: wenxiao 时间: 2020-3-8 22:10

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

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

作者
Author: xp47 时间: 2020-3-9 10:02
本帖最后由 xp47 于 2020-3-9 10:27 编辑

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函数下的参数手动改掉然后编译下.

作者
Author: wenxiao 时间: 2020-3-9 13:55

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

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

作者
Author: xp47 时间: 2020-3-9 14:05

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

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

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