计算化学公社

标题: 求助：DFT-D3参数设置不识别 [打印本页]

作者
Author: DoubeeTwT 时间: 2020-7-12 22:00
标题: 求助：DFT-D3参数设置不识别
看了sob老师的博文想要在m062x上加上色散校正
没有m062xD3
em=GD3显示没有S6 S8 SR6 参数
查找高斯手册找到对应的数值S6=1.0000 S8=0.0000 SR6=1.6190
使用IOp输入参数发现一个问题IOp(3/175=0)识别不了，那我如何输入S8=0这个参数呢？求教

作者
Author: hebrewsnabla 时间: 2020-7-12 22:04
如果你是G09D或以后的版本，不可能没有M062X的D3参数。如果是之前的版本，当然不识别，因为根本不支持D3。

作者
Author: ldatea 时间: 2020-7-12 22:10
本帖最后由 ldatea 于 2020-7-12 22:20 编辑

sob老师在《DFT-D色散校正的使用》http://sobereva.com/210一文中已经讲的很清楚了，G09 D.01版以前不支持DFT-D3，想要用只能通过外部程序，比如Grimme的DFT-D3程序，目前最新版本3.2.0。sob老师在文中给了此程序的网址，以及已经编译好的3.0.2版本和最基础的使用方法（Linux和Windows都有），程序能直接给出D3的能量和梯度，但没法直接给出Hessian。
M06-2X的零阻尼形式（M06-2X只能用零阻尼）在G09 D.01中已经有内置参数了，请检查G09版本。我感觉你看的是g16的手册，g16才能通过IOP修改D3的参数。老版本的用户手册论坛上就有。

作者
Author: DoubeeTwT 时间: 2020-7-12 22:21

hebrewsnabla 发表于 2020-7-12 22:04
如果你是G09D或以后的版本，不可能没有M062X的D3参数。如果是之前的版本，当然不识别，因为根本不支持D3。

输入：
#p opt=(calcfc,maxcyc=60) freq oniom(m062x/6-31g(d):amber) em=GD3

报错：
Raffenetti 1 integral format.
Two-electron integral symmetry is turned on.
R6DS8: Unable to choose the S8 parameter, IExCor= 0 IXCFnc= 0 ScaHFX= 1.000000 IDFTD=3
Error termination via Lnk1e in /lustre/home/acct-clsst/clsst/Software/g16.A.03/g16/l301.exe at Sun Jul 12 22:20:22 2020.

作者
Author: DoubeeTwT 时间: 2020-7-12 22:22

ldatea 发表于 2020-7-12 22:10
sob老师在《DFT-D色散校正的使用》http://sobereva.com/210一文中已经讲的很清楚了，G09 D.01版以前不支持D ...

我用的是16

作者
Author: ldatea 时间: 2020-7-12 22:41
本帖最后由 ldatea 于 2020-7-12 22:42 编辑

DoubeeTwT 发表于 2020-7-12 22:21
输入：
#p opt=(calcfc,maxcyc=60) freq oniom(m062x/6-31g(d):amber) em=GD3

我试了一下，不用oniom是可以的。我猜这可能跟你用了oniom有关。我用的是C.01版的

作者
Author: DoubeeTwT 时间: 2020-7-12 22:55

ldatea 发表于 2020-7-12 22:41
我试了一下，不用oniom是可以的。我猜这可能跟你用了oniom有关。我用的是C.01版的

原来是这样
现在的问题是我通过IOp确实可以给定参数，但是S8=0我直接设置3/175=0不能识别，但是给3/175=1就可以正常运行，此时S8=0.000001不是严格意义上的0

作者
Author: ldatea 时间: 2020-7-12 23:03
本帖最后由 ldatea 于 2020-7-12 23:11 编辑

DoubeeTwT 发表于 2020-7-12 22:55
原来是这样
现在的问题是我通过IOp确实可以给定参数，但是S8=0我直接设置3/175=0不能识别，但是给3/175= ...

IOp(3/175)
S8 scale factor in Grimme’s D2/D3/D3BJ dispersion.
0       Default (see subroutine R6DS8).
-1       Set S8 to 0.
NNNNNNNN          A value of NNNNNNNN/1000000.
这是论坛上IOP的离线版本，你设0等于什么都没干，要设成0就应该iop(3/175=-1)
http://bbs.keinsci.com/thread-17394-1-1.html这个帖子正好就是2楼发的。oniom我没用过，不排除你用oniom有误的情况。

作者
Author: indec 时间: 2020-7-13 06:50
在 oniom 里使用 D3 校正需要按如下写

#p opt=(calcfc,maxcyc=60) freq oniom(m062x/6-31g(d) em=gd3:amber)

我在 A.03 里如此测试是能算下去的，em=gd3 写在外面，就会报一样的错

作者
Author: DoubeeTwT 时间: 2020-7-13 09:59

indec 发表于 2020-7-13 06:50
在 oniom 里使用 D3 校正需要按如下写

我在 A.03 里如此测试是能算下去的，em=gd3 写在外面，就会报一样 ...

有道理哦~可以写里面这样就只对高层使用，学会了学会了~

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