计算化学公社

 找回密码 Forget password
 注册 Register
Views: 6333|回复 Reply: 0

[Gaussian/gview] Gaussian有个坑爹的选项UnConAO

[复制链接 Copy URL]

597

帖子

20

威望

4151

eV
积分
5148

Level 6 (一方通行)

发表于 Post on 2019-7-31 20:40:02 | 显示全部楼层 Show all |阅读模式 Reading model
本帖最后由 beefly 于 2019-8-1 02:02 编辑

像Gaussian、Gamess、Orca这类程序处理general contraction类型的基组是非常低效的。比如用HF/ANO-RCC-VDZ算At2的单点能,Molpro、Molcas几分钟就算完了,而前面的几个程序要算几个小时(甚至可能超过一天),甚至还不如用非收缩形式快。Turbomole、Orca开发的def-、def2-、sarc、x2c-等类型基组其实就是一种折衷的办法,仅对1s、2p、3d等最低芯层比较重要的几个基函数(一般小于10)做了收缩,其余的大部分基函数都是非收缩的,这样不仅效率比较快(虽然不如Pople型基组),而且半芯层、价层的精度比较高。

Gaussian使用非收缩基组的一种办法是把基组手工编辑为非收缩形式以后,通过gen关键词读入。另一个办法是设置Int(UnConAO),原始的收缩基组既可以来自内部基组库,也可以通过gen读入。但是UnConAO存在一些问题。


在Gaussian09比较早的几个版本里,UnConAO仅能用于Pople型基组。general contraction基组有大量重复的高斯指数,UnConAO并没有删除它们,导致基组有严重的线性依赖。


在Gaussian09晚期的版本和Gaussian16里,UnConAO仅能保证general contraction基组非收缩后的形式正确,但基函数在程序内部的存储并未优化,所以仍然非常慢。还是前面用HF/ANO-RCC-VDZ算At2单点的例子,用Gaussian16.b,单节点8核。作为极端的测试,这里做了二分量相对论DKH4(ghf/gen int=dkhso),这种情况下程序自动关闭对称性。
  • Gen读入收缩基组:将近1小时。
  • Gen读入收缩基组+UnConAO:算完单电子积分之后没反应,等了半小时后强行终止。
  • Gen读入非收缩基组:3分钟。

评分 Rate

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

查看全部评分 View all ratings

本版积分规则 Credits rule

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

GMT+8, 2023-2-2 22:01 , Processed in 0.179475 second(s), 22 queries .

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