计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9754|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 求助:DFT-D3参数设置不识别

[复制链接 Copy URL]

138

帖子

0

威望

445

eV
积分
583

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
看了sob老师的博文想要在m062x上加上色散校正
没有m062xD3
em=GD3显示没有S6 S8 SR6 参数
查找高斯手册找到对应的数值S6=1.0000 S8=0.0000 SR6=1.6190
使用IOp输入参数发现一个问题IOp(3/175=0)识别不了,那我如何输入S8=0这个参数呢?求教

上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

138

帖子

0

威望

445

eV
积分
583

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2020-7-13 09:59:01 | 只看该作者 Only view this author
indec 发表于 2020-7-13 06:50
在 oniom 里使用 D3 校正需要按如下写

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

有道理哦~可以写里面这样就只对高层使用,学会了学会了~
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

155

帖子

0

威望

3433

eV
积分
3588

Level 5 (御坂)

9#
发表于 Post on 2020-7-13 06:50:28 | 只看该作者 Only view this author
在 oniom 里使用 D3 校正需要按如下写
#p opt=(calcfc,maxcyc=60) freq oniom(m062x/6-31g(d) em=gd3:amber)

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

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1 正解

查看全部评分 View all ratings

115

帖子

1

威望

1122

eV
积分
1257

Level 4 (黑子)

8#
发表于 Post on 2020-7-12 23:03:40 | 只看该作者 Only view this author
本帖最后由 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有误的情况。

138

帖子

0

威望

445

eV
积分
583

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2020-7-12 22:55:09 | 只看该作者 Only view this author
ldatea 发表于 2020-7-12 22:41
我试了一下,不用oniom是可以的。我猜这可能跟你用了oniom有关。我用的是C.01版的

原来是这样
现在的问题是我通过IOp确实可以给定参数,但是S8=0我直接设置3/175=0不能识别,但是给3/175=1就可以正常运行,此时S8=0.000001不是严格意义上的0
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

115

帖子

1

威望

1122

eV
积分
1257

Level 4 (黑子)

6#
发表于 Post on 2020-7-12 22:41:27 | 只看该作者 Only view this author
本帖最后由 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版的

138

帖子

0

威望

445

eV
积分
583

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2020-7-12 22:22:15 | 只看该作者 Only view this author
ldatea 发表于 2020-7-12 22:10
sob老师在《DFT-D色散校正的使用》http://sobereva.com/210一文中已经讲的很清楚了,G09 D.01版以前不支持D ...

我用的是16
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

138

帖子

0

威望

445

eV
积分
583

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2020-7-12 22:21:42 | 只看该作者 Only view this author
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.
上海交通大学计算化学与分子生物信息学实验室
Shanghai JiaoTong University
Computational Chemistry and Molecular Bioinformatics Laboratory

CCMBI of SJTU
点击访问个人主页

115

帖子

1

威望

1122

eV
积分
1257

Level 4 (黑子)

3#
发表于 Post on 2020-7-12 22:10:21 | 只看该作者 Only view this author
本帖最后由 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的参数。老版本的用户手册论坛上就有。

928

帖子

1

威望

8262

eV
积分
9210

Level 6 (一方通行)

2#
发表于 Post on 2020-7-12 22:04:28 | 只看该作者 Only view this author
如果你是G09D或以后的版本,不可能没有M062X的D3参数。如果是之前的版本,当然不识别,因为根本不支持D3。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-19 15:04 , Processed in 0.467416 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list