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

计算化学公社

 找回密码
 现在注册!
查看: 1115|回复: 2

[量化理论] 给基组以even-tempered方式增加弥散函数的工具adddiffuse

[复制链接]

1万

帖子

25

威望

1万

eV
积分
35687

管理员

公社社长

发表于 2016-9-23 20:41:27 | 显示全部楼层 |阅读模式
给基组以even-tempered方式增加弥散函数的工具adddiffuse

文/Sobereva @北京科音  2016-Sep-23


有的基组本身没弥散函数,或者虽然有弥散函数但只有一层,比如aug-cc-pVTZ,而对于精确计算超极化率、里德堡激发态的目的一层弥散函数不够。

有一种构建基组的方式称为even-tempered,即相邻GTF的指数是成比例的。若当前基组中某个角动量函数中指数最小的是ζ(n),第二小的是ζ(n-1),则这个角动量的弥散函数的指数应当为ζ(n)*[ζ(n)/ζ(n-1)]。

例如aug-cc-pVTZ对碳的三层D极化函数的指数:
1.0970000
0.3180000
0.1000000
指数最小的两个的比例是0.1/0.318=0.314,因此,再增加一层弥散函数(构成d-aug-cc-pVTZ),指数应当为0.1*0.314=0.0314,如果再增加一层(构成t-aug-cc-pVTZ),就是0.1*0.314^2=0.00986。

笔者写的adddiffuse是用来按照上述方式增加弥散函数的。下载地址:http://sobereva.com/soft/adddiffuse.rar,其中含编译好的windows版和源程序。

启动后需要输入含有Gaussian格式的基组定义的文件。可以从EMSL上拷出来放到文本文件里。文件中可以包含多个元素,每个元素末尾都要有****。压缩包里的aug-cc-pVTZ.txt就是个例子。启动程序后,输入这个文件名,然后选择要增加几层弥散函数,以及对什么角动量增加弥散函数(比如只增加S、P、D弥散就输入SPD,如果输入A就对所有角动量都增加弥散)。然后程序就会在当前目录下产生new.txt,这是添加弥散函数后的Gaussian格式的基组定义。

new.txt里面每个元素前头多了个负号,这样做的好处是,用gen使用自定义基组的时候,若把一堆元素的基组定义放在坐标后头,即使其中有的元素没有在当前计算中出现Gaussian也不会报错。adddiffuse不支持SP壳层。如果某个角动量只含一层基函数,则adddiffuse也不会对之增加弥散函数。

或许有人觉得这程序没多大用,因为用计算器点几下就得到弥散函数的指数了。但是,当你同时需要对多个基组、多个元素都增加一层乃至更多层弥散的时候,手工去弄就很痛苦了,还容易出错。

此程序一个重要用处是基于EMSL上的aug-的Dunning相关一致性基组产生d-aug-和t-aug-版本。虽然EMSL上也有d-aug-的定义,但是撰文时只对前两周期有定义,而t-aug-则在EMSL上根本没有。另外,虽然Gaussian里面内置了d-aug-的定义,比如写CCSD/daug-cc-pVTZ,但是其最外层弥散函数的指数和EMSL上的不符。EMSL上的d-aug-的最外层指数是通过上述even-tempered方式构建的,而Gaussian里的不知道怎么来的,大抵没even-tempered方式构建的合理。

评分

参与人数 4eV +20 收起 理由
978142355 + 5 牛!
captain + 5 谢谢
jliang + 5 好物!
我本是个娃娃 + 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!

2363

帖子

9

威望

4027

eV
积分
6570

Level 6 (一方通行)

首席卖萌官

发表于 2016-9-23 22:10:51 | 显示全部楼层
大宝博士知识的渊博,我等实在是无法望其项背
She doesn't love me.
Even so,
my heart has been taken away by her.

39

帖子

0

威望

957

eV
积分
996

Level 4 (黑子)

发表于 2016-9-24 09:29:14 | 显示全部楼层
sob老师真伟大!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2018-11-15 08:57 , Processed in 0.108216 second(s), 26 queries .

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