计算化学公社

标题: 关于高斯中LC-BLYP一个参数的设定问题 [打印本页]

作者
Author: ter20 时间: 2015-5-14 16:20
标题: 关于高斯中LC-BLYP一个参数的设定问题
文章审稿意见中，审稿人提到了用LC-BLYP泛函来测试我的体系，我查了他提供的文献，都提到了一个参数"range separating parameter μ"需要设为0.33 bohr-1 (这个单位也很诡异，作者的几篇文献不统一，有的直接不写单位，不知道这里是不是作者写错了，难道不该是无量纲的？)，具体如下图所示：

我查了高斯手册和IOP手册，都没有找到这个参数的设定，求Sob大神指教！不胜感激！

作者
Author: 卡开发发 时间: 2015-5-14 16:41
本帖最后由卡开发发于 2015-5-14 16:58 编辑

IOp(3/107)
Omega for short/long range Hartree-Fock exchange.
0 Standard HF exchange
MMMMMNNNNN Short range HF exchange with NNNNN/10000 and long range exchange with
MMMMM/10000.

IOp(3/108)
Omega for short/long range DFT exchange.
0 Standard DFT exchange or default from functional.
MMMMMNNNNN Short range DFT exchange with NNNNN /10000 and long range DFT exchange
with MMMMM/10000.

IOp(3/109)
Omega for short/long range DFT correlation
0 Standard DFT correlation or default from functional.
MMMMMNNNNN Short range DFT correlation with NNNNN /10000 and long range DFT
correlation with MMMMM/10000.

当中的omega就是μ，估计应该是IOp(3/107)=0330000000和IOp(3/108)=0000003300 IOp(3/109)=0000003300吧。

作者
Author: ter20 时间: 2015-5-14 17:11
我也查到了感觉楼上说的这三个参数，但是感觉不太像，难道是要三个IOp选项中的所有long/short都设成0.3吗？
另外原公式中是μ，这里看起来是ω，不是太懂，求解

作者
Author: 卡开发发 时间: 2015-5-14 17:16

ter20 发表于 2015-5-14 17:11
我也查到了感觉楼上说的这三个参数，但是感觉不太像，难道是要三个IOp选项中的所有long/short都设成0.3吗？ ...

就是这个参数，如J. Chem. P hy s. 136, 134704 (2012)中就是ω表示。但数值上是不是这样设定我不完全确定。

作者
Author: sobereva 时间: 2015-5-14 17:50
ω和μ都有人用，无所谓，是一回事
单位是bohr^-1，否则erf里不是无量纲的。
将3/107和3/108都设为MMMMM00000，就相当于用了ω参数为MMMMM/10000的范围分离泛函。对GGA加LC-时ω默认是0.47/Bohr，等价于0470000000。注意3/107和3/108在IOp手册的字面上的解释有一定误导性，不要管。
3/109不用设。

作者
Author: ter20 时间: 2015-5-14 21:03

sobereva 发表于 2015-5-14 17:50
ω和μ都有人用，无所谓，是一回事
单位是bohr^-1，否则erf里不是无量纲的。
将3/107和3/108都设为MMMMM0 ...

多谢Sob和卡开发发，确认一下，如果我要把μ设为0.3，是否就是用如下关键字：
LC-BLYP IOP(3/107)=0330000000 IOP(3/108)=0330000000

作者
Author: sobereva 时间: 2015-5-14 23:55
应为0300000000

作者
Author: ter20 时间: 2015-5-15 09:37

sobereva 发表于 2015-5-14 23:55
应为0300000000

囧，笔误，真是太感谢Sob了！大赞！

作者
Author: 小范范1989 时间: 2015-9-18 11:49

sobereva 发表于 2015-5-14 23:55
应为0300000000

sob老师好，卡开发发老师好@卡开发发，我也遇到了和楼主一样的问题，而且我也看过老师您写的那个IOP一览的帖子，我想咨询一下几个问题：
1：我也看到了和楼主一样的文献（JPCC2015.119.13535-13544），公式都一样，里面用到了LC-BLYP泛函，默认的u=0.33bohr-1。我想咨询一下老师，这个普通的泛函前面加上LC-是什么意思？
2：同时，我看到另外一篇文献（JCTC2015.11.3851-3858），里面采用的是LC-wPBE，w默认是0.4，如果我想咨询一下，我想调剂w的数值，是不是和楼主这个帖子中u的调剂方式一样IOP（3/107） IOP（3/108）。
3：通过这个平台的帖子，收获太多了。thanks

作者
Author: sobereva 时间: 2015-9-18 13:30

小范范1989 发表于 2015-9-18 11:49
sob老师好，卡开发发老师好@卡开发发，我也遇到了和楼主一样的问题，而且我也看过老师您写的那个IOP一览 ...

1 LC-是对GGA泛函加上长程校正，远程交换部分用100%的HF成份。wB97X实际上也是这样。w参数用来调节近程和远程怎么进行分离。你可以看LC的原文JCP, 115, 3540。

2 应该是。你可以再看看LC-wPBE的原文确认一下是不是指的LZ图中的那个参数。

作者
Author: 小范范1989 时间: 2015-11-28 10:57
本帖最后由小范范1989 于 2015-11-28 11:30 编辑

ter20 发表于 2015-5-14 21:03
多谢Sob和卡开发发，确认一下，如果我要把μ设为0.3，是否就是用如下关键字：
LC-BLYP IOP(3/107)=0330 ...

楼主好，@sobereva sob老师好。按照楼主的方法，我这样计算，可行，但是我想换位LC-pw91，如下的错误 ----------------------------------------------------------
#p lc-pw91/6-31g(d) iop(3/107=0203600000,3/108=0203600000)
----------------------------------------------------------
QPErr --- An ambiguous keyword was detected.
#p lc-pw91/6-31g(d) iop(3/107=0203600000
'
Last state= "DFT"

2：同时我想咨询一下，这样加LC的泛函，除了LC-pbe，LC-wpbe，LC-BLYP，是不是还有LC-pw91，LC-tpss？还有吗？我想自己用另外一个泛函计算，不知道可行不可行，谢谢指点。

作者
Author: sobereva 时间: 2015-11-28 17:05

小范范1989 发表于 2015-11-28 10:57
楼主好，@sobereva sob老师好。按照楼主的方法，我这样计算，可行，但是我想换位LC-pw91，如下的错误 --- ...

PW91得写PW91PW91

作者
Author: 小范范1989 时间: 2015-11-28 17:20

sobereva 发表于 2015-11-28 17:05
PW91得写PW91PW91

偶，确实测试了一下是这样，谢谢老师。

作者
Author: 初初 时间: 2019-5-22 19:03
各位老师好！我在回复审稿意见时也遇到了同样的问题，除了上面提到的μ=0.33设定的问题，审稿人给我的参考文献同时提到了α和β参数的设定，并测试了两种情况α=0.0，β=1.0和α=0.2，β=0.8，并提出后面的比例结果更好，想问下各位老师α和β的参数改如何设定？不胜感激！J. Chem. Theory Comput. 2015, 11, 2199− 2209

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