计算化学公社

标题: 求助Q-Chem 4如何修改长程校正泛函的omega值 [打印本页]

作者
Author: ChunLinX 时间: 2020-5-7 23:12
标题: 求助Q-Chem 4如何修改长程校正泛函的omega值
本帖最后由 ChunLinX 于 2020-5-11 23:59 编辑

各位老师好，我想在Q-Chem 4中使用如ωB97XD这类的长程校正泛函，但对于长程校正泛函一般需要进行ω值tuning。
我查看了Q-Chem 4的说明书，里面有一个关键词为Omega，可以用来修改ω值，但这个关键词好像只能适用于LRC-wPBE
和BNL。我试过在使用ωB97XD时加入Omega命令，但是输出文件的ω并没有变化。我看了之前有在Q-Chem中使用

ωB97XD并修改ω值的文章。

请问应该用什么关键词来在Q-Chem中修改像ωB97XD这类长程校正泛函的ω值呢？谢谢老师！

作者
Author: myzpku 时间: 2020-5-11 10:33
wB97X-D这一泛函由于在training的时候已经fix omega = 0.2 a0^-1, 并不是非常适合用来做omega tuning。用LRC-wPBE或者LC-BLYP会比较好。

如果一定要改变wB97X-D的omega值，可以用如下方法:
在$rem section设置以下四个关键字
exchange = gen
lrc_dft = true (range-separated funcitonal)
omega = xxx (1000*omega值, 设成200即为标准的wB97X-D）
dft_d = empirical_chg (对应的dispersion correction)

然后再加上一个自定义xc functional的section
$xc_functional
X HF 0.222036
X wB97X-D 1.0
C wB97X-D 1.0
$end

在此基础上改变omega值即可实现tuning。另外我发现如果用最新的Q-Chem 5.2，在method = wB97X-D时直接改omega值也是有效的。

作者
Author: ChunLinX 时间: 2020-5-12 00:17

myzpku 发表于 2020-5-11 10:33
wB97X-D这一泛函由于在training的时候已经fix omega = 0.2 a0^-1, 并不是非常适合用来做omega tuning。用LR ...

谢谢老师！但是我试了一下，输出文件报错了“Invalid functional created by AddX” ，这是不是说明要把wB97X-D完全拆开来写入$xc_functional部分呢？但感觉这样又太麻烦，而且容易出错。

作者
Author: myzpku 时间: 2020-5-12 04:13
我研究了一下，的确用qchem 4似乎做不了wB97X-D的omega tuning。我昨天说的在当前的q-chem 5.2下是可以正常运行的（见附件） (, 下载次数 Times of downloads: 16)

(, 下载次数 Times of downloads: 18)

作者
Author: ChunLinX 时间: 2020-5-12 09:44

myzpku 发表于 2020-5-12 04:13
我研究了一下，的确用qchem 4似乎做不了wB97X-D的omega tuning。我昨天说的在当前的q-chem 5.2下是可以正常 ...

应该就是版本的问题，我就用其它泛函吧。谢谢老师的回答！

作者
Author: GF1995 时间: 2021-9-6 22:59

myzpku 发表于 2020-5-11 10:33
wB97X-D这一泛函由于在training的时候已经fix omega = 0.2 a0^-1, 并不是非常适合用来做omega tuning。用LR ...

你好，wB97XD泛函不应该是0.22短程HF交换+0.78短程wB97XD交换+1.0长程HF交换+1.0wB97XD相关势吗？为什么在X wB97X-D 后面写1.0呢？

作者
Author: myzpku 时间: 2021-9-8 07:12

GF1995 发表于 2021-9-6 22:59
你好，wB97XD泛函不应该是0.22短程HF交换+0.78短程wB97XD交换+1.0长程HF交换+1.0wB97XD相关势吗？为什么 ...

wB97X-D短程交换项的系数(0.78)应该是被吸收到泛函本身了，关于这一点可以看作者在此之前wB97X的文章中的Eq. 28: https://aip.scitation.org/doi/full/10.1063/1.2834918

作者
Author: GF1995 时间: 2021-9-8 11:46

myzpku 发表于 2021-9-8 07:12
wB97X-D短程交换项的系数(0.78)应该是被吸收到泛函本身了，关于这一点可以看作者在此之前wB97X的文章中的 ...

好的，谢谢您，我再研究一下

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