请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 2651|回复: 12

[综合交流] 给ahlrichs的def2系列基组加弥散的方法

[复制链接]

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

发表于 2016-6-8 10:58:03 | 显示全部楼层 |阅读模式
给ahlrichs的def2系列基组加弥散的方法

文/Sobereva @北京科音   2016-Jun-8


ahlrichs等人2005年搞的def2系列基组(PCCP,7,3297)如今比较流行,构建得不错,比较普适,前六周期完全覆盖(镧系锕系除外),同等计算量下精度秒杀pople基组,比Dunning相关一致性基组明显更适合DFT计算,还有对应的密度拟合基组从而在ORCA/turbomole程序中跑得速度如飞,是DFT计算的首选基组之一。但是开发者并没专门提出def2系列带弥散的版本,而有些问题又必须有弥散函数(讨论见《谈谈弥散函数和“月份”基组》http://sobereva.com/119),经常有人问怎么给def2系列基组加弥散,这里就简单说一下。

1 Dunning相关一致性基组的弥散函数可以直接挪到Alhrichs基组上,例如加到def2-QZVP上可以称为aug-def2-QZVP,例如Phys. Chem. Chem. Phys., 13, 6670-6688 (2011)就这么用。也可以按照月份基组的思路,砍掉aug上不重要的高角动量弥散函数以降低耗时。如果是Gaussian用户不知道具体怎么做的话,仔细阅读和理解此文《详解Gaussian中混合基组、自定义基组和赝势基组的输入》(http://sobereva.com/60)。另外,Frank Jensen搞的aug-pc-n系列基组的弥散函数也可以借到def2上,由于aug-pc-n基组是专给DFT计算优化的,故对于DFT计算来说可能结果比借用aug-cc-pVnZ的弥散效果更好。

2 在J. Chem. Phys., 133, 134105 (2010)中作者提出def2-SVPD、TZVPD、TZVPPD、QZVPD、QZVPPD,分别是对def2-SVP、TZVP、TZVPP、QZVP、QZVPP加上弥散函数,弥散函数的指数是通过优化原子的HF极化率得到(同样适合DFT下极化率计算)。这些基组用于其它需要弥散函数的任务,比如算弱相互作用,预期也能起到不错效果,但未必比本文其它的加弥散方式做法更好。这些带D的def2基组可以在EMSL上获得。

3 Truhlar等人提出一种通用的给出原本不含弥散函数的基组以最低限度的弥散函数的策略,这样的基组以ma-开头,含义是minimal augmentation,详见J. Chem. Theory Comput., 7, 3027、Theor. Chem. Acc., 128, 295。也就是将原先基组中指数最小的s和p的指数除以3作为弥散函数的指数,但不对氢加弥散函数。这种处理后的def2基组可以直接从这里获得Gaussian格式的定义:http://comp.chem.umn.edu/basissets/basis.cgi。如ma-def2-TZVP在里面被简写为ma-TZVP。对于大多数需要弥散函数的问题,实际上像这样只给重原子加一层s和p弥散就已经解决绝大部分问题了,因此如果想给def2加弥散但又不想花费过高代价的话,这种ma-方式加弥散是比较理想的。

4 若当前基组中某个角动量函数中指数最小的是ζ(n),第二小的是ζ(n-1),则这个角动量的弥散函数的指数应当为ζ(n)*[ζ(n)/ζ(n-1)]。这称为Even-tempered方法,见Int. J. Quantum Chem., 113, 21-34 (2013)式3,类似于靠几何级数关系来确定更弥散的函数的指数。这是比较普适的“自助”方式给原本没弥散函数的基组加各角动量弥散的做法。


顺带一提,在Int. J. Quantum Chem., 116, 1084 (2016)中,作者提出了R-ORP基组,是在def2基础上加了一层s和p弥散同时调节了指数,目的是让基组大小在aug-cc-pVDZ的级别但计算(超)极化率性能有明显提升。实际测试表明此基组算极化率效果一般,但算第一超极化率比较碉,根据作者的测试达到>=daug-cc-pVDZ的水准。不过此基组目前只对C、H、O、N、F有定义。

另外,def2系列基组对弥散函数的要求低于pople系列基组。特别是在3-zeta级别,def2加ma后改进远比pople系列基组加+后的改进要小。所以像计算势垒这种加弥散有益但不是必须的情况,如果对精度要求不是很高的话def2不加弥散也无大碍。

评分

参与人数 7eV +34 收起 理由
leebo + 5 好物!
bianwenbo + 4 牛!
bomsaude + 5 谢谢
978142355 + 5 相当的给力
xpyp + 5 赞!
我本是个娃娃 + 5 好物!
北纬18° + 5 好物!

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

2331

帖子

9

威望

3887

eV
积分
6398

Level 6 (一方通行)

首席卖萌官

发表于 2016-6-8 13:01:30 | 显示全部楼层
抢沙发,鼓掌!
为雪国耻身先去,重整河山待后生。

122

帖子

0

威望

2762

eV
积分
2884

Level 5 (御坂)

发表于 2016-6-8 13:03:02 | 显示全部楼层
计算弱相互作用,用了def2-qzvp, 是否还需要加弥散函数呢?

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

 楼主| 发表于 2016-6-8 13:07:59 | 显示全部楼层
xpyp 发表于 2016-6-8 13:03
计算弱相互作用,用了def2-qzvp, 是否还需要加弥散函数呢?


如果是DFT-D3校正的泛函级别来算的话,不用加,而且D3本身拟合参数的时候就是在这个基组下的。如果用高级后HF方法要求贼精确,当金标准,而且是色散主导的作用的话,可以再加个aug-或者考虑BSSE
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

12

帖子

0

威望

103

eV
积分
115

Level 2 能力者

发表于 2016-6-8 18:30:47 | 显示全部楼层
本帖最后由 bianwenbo 于 2016-6-8 18:32 编辑

@sobereva 我感觉加了弥散的基组有这样的排序(都是乞丐级别):aug-cc-pVDZ > ma-def2-SVP ≈ 6-31+G*
至少我在计算单点能的时候使用同样参数,前面Dunning的耗时大概是中间Ahlrichs这个的三倍。

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

 楼主| 发表于 2016-6-8 19:01:40 | 显示全部楼层
bianwenbo 发表于 2016-6-8 18:30
@sobereva 我感觉加了弥散的基组有这样的排序(都是乞丐级别):aug-cc-pVDZ > ma-def2-SVP ≈ 6-31+G*
至 ...


因为aug-cc-pVDZ对d角动量、对氢都有弥散函数,ma就只有sp的弥散,所以耗时会高一截
你可以用月份基组
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

283

帖子

2

威望

2744

eV
积分
3067

Level 5 (御坂)

发表于 2016-6-11 10:31:07 | 显示全部楼层

还是玩点真东西吧

240

帖子

1

威望

1901

eV
积分
2161

科音成员

发表于 2016-8-1 16:38:36 | 显示全部楼层
本帖最后由 helpme 于 2016-8-24 11:25 编辑

感谢sob介绍的方法。

针对使用ma-def2tzvp基组,提供一种具体的操作方式(另外的方式比如使用genecp当然也可以)。

1. 首先得到def2tzvp基组的具体形式(可以在网站下载,或者用高斯的def2tzvp + gfinput来输出)。把s和p的指数系数最小的两项系数除以3,作为弥散函数。(这个也可以到http://comp.chem.umn.edu/basissets/basis.cgi 下载,但是可能有些元素不全)。

2. 编辑一个基组文件 ma-def2tzvp.gbs,放在$g09root/basis下面内容包含所研究的所有原子的弥散函数信息(这个文件可以慢慢积累,以后还可以用)。包含C和O的文件的ma-def2tzvp.gbs文件如下(注意元素符号前面的-号,这样即使体系不含该元素,程序也不会报错退出):
-C 0
S   1 1.00       0.000000000000
      0.3172148001D-01  0.1000000000D+01
P   1 1.00       0.000000000000
      0.3352274557D-01  0.1000000000D+01
****
-O 0
S   1 1.00       0.000000000000
      0.6168178787D-01  0.1000000000D+01
P   1 1.00       0.000000000000
      0.5826140423D-01  0.1000000000D+01
****

3. 使用的时候,大致是这个样子(考虑到弥散函数常常用于弱相互作用的计算,这里把B3LYP-D3(BJ)的形式也写进去了):

# b3lyp/def2tzvp extrabasis empiricaldispersion=gd3 opt freq

test

0 1
C                  0.00000000    0.42348200    0.00000000
O                 -1.12068900   -0.21083400    0.00000000
。。。

@/home/soft/g09/basis/ma-def2tzvp.gbs /N

上面最后一行如果不加 /N ,那么会把ma-def2tzvp.gbs 文件里面的所有内容写到log文件里面,看上去不好看。                                                                                                                    

评分

参与人数 1eV +5 收起 理由
sobereva + 5

查看全部评分

华北电力大学数理学院,理论与计算化学,团簇结构与反应机理。欢迎报考硕士研究生。

186

帖子

0

威望

1927

eV
积分
2113

Level 5 (御坂)

发表于 2017-8-29 11:30:45 | 显示全部楼层
helpme 发表于 2016-8-1 16:38
感谢sob介绍的方法。

针对使用ma-def2tzvp基组,提供一种具体的操作方式(另外的方式比如使用genecp当然 ...

丁老师,学生有个疑问,您给出的ma-TZVP文件内容和网站上下载的不太一样,网站上给出的信息更多一些,写了更多的s和p轨道的定义,甚至给出了d和f轨道的定义,这跟您写的有什么不同,对于计算有什么影响?

ma-tzvp.png
站在宇宙中心呼唤爱

240

帖子

1

威望

1901

eV
积分
2161

科音成员

发表于 2017-9-20 08:24:26 | 显示全部楼层
kevin 发表于 2017-8-29 11:30
丁老师,学生有个疑问,您给出的ma-TZVP文件内容和网站上下载的不太一样,网站上给出的信息更多一些,写 ...

我给的例子里面,关键词加了 extrabasis,表示在原有基组的基础上,额外加一些基函数(用来弥散的),所以就不用再把def2tzvp本身包含的基函数写一遍了。你写的很多很长的部分,实际上就是原来def2tzvp包含的。

回复晚了不好意思。最近只是浏览网站没有登陆,这才看到。
华北电力大学数理学院,理论与计算化学,团簇结构与反应机理。欢迎报考硕士研究生。

186

帖子

0

威望

1927

eV
积分
2113

Level 5 (御坂)

发表于 2017-9-21 17:34:33 | 显示全部楼层
本帖最后由 kevin 于 2017-9-22 09:05 编辑
helpme 发表于 2017-9-20 08:24
我给的例子里面,关键词加了 extrabasis,表示在原有基组的基础上,额外加一些基函数(用来弥散的),所 ...

谢谢丁老师,您的意思是,您的方法是 原有的TZVP + ma-从网站http://comp.chem.umn.edu/basissets/basis.cgi下载下来的是整个的ma-TZVP
如果用网站里的内容编辑成ma-TZVP.gbs文件,输入文件的写法是下面这样么?
# b3lyp  extrabasis empiricaldispersion=gd3 opt freq



test


0 1
坐标


@/home/soft/g09/basis/ma-def2tzvp.gbs /N
站在宇宙中心呼唤爱

240

帖子

1

威望

1901

eV
积分
2161

科音成员

发表于 2017-9-25 15:34:58 | 显示全部楼层
# b3lyp/def2tzvp extrabasis
empiricaldispersion=gd3 opt freq

就是命令行要把 def2tzvp和 extrabasis都写上。前者表示正常版本的def2tzvp,后者表示还要加点东西——要加的东西就是文件末尾那个了。
华北电力大学数理学院,理论与计算化学,团簇结构与反应机理。欢迎报考硕士研究生。

5

帖子

0

威望

144

eV
积分
149

Level 2 能力者

发表于 2018-4-19 11:01:35 | 显示全部楼层
本帖最后由 熊琳 于 2018-4-19 11:02 编辑
helpme 发表于 2017-9-25 15:34
# b3lyp/def2tzvp extrabasis
empiricaldispersion=gd3 opt freq

初学者刚接触Gaussian中的def2tzvp这个基组,不知道针对其的关键词如何使用,之前写的是“def2-tzvp”,出现语法报错的消息。然后刚好看到这个帖子,去掉“-”就OK了。嘿嘿,谢谢老师
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-9-26 13:47 , Processed in 0.206567 second(s), 27 queries .

快速回复 返回顶部 返回列表