高精度 W1U 方法的生成焓计算
天泽2018/08/22化学 IP:美国

介绍:

先感谢一下 dracula1429 在计算上的帮助。W1 类理论是计算化学里,气态生成焓计算精度最高的方法之一。W1U 的气态生成焓和实验结果相比,平均误是惊人的 0.3 kcal ,而却完全不带入经验法。

W1 使用个方法使用 coupled cluster (CC) 数值技术。CC 采用基本的Hartree-Fock 分子轨道方法,并使用指数聚类算子构造多电子波函数以解释电子相关性。该计算类似于CBS-X,基函数无限制外推。 

W1U 的高精度也是由于使用了非常大的基组:高达cc-pVQZ + 2d1g和  cc-pV5Z + 2d1f 以及CCSD和CCSD(T)水平,所以如果计算材料超过2个重原子(碳,氧,氮,等等在W1U 里都算是重原子)的材料建议使用超级电脑或者高核心数量工作站。以下计算是使用Gaussian 09在中等配置私人电脑 ( 8 线程Xeon E3-1505M v5 带有 32GB RAM )得出的结果。


计算结果:

以下是7个小型材料的生成焓计算。其中,HN3,拥有3个重原子,和 C2H6 ,2个重原子,暂时是个人最大的计算材料,C2H6 耗时 4小时25 分钟。 1H-四唑跑了8个小时候出现了 diskspace 错误, 肯定是超越我电脑的极限了。

所有气态生成焓是在 NIST 数据库里找出。其中,HN3 的 70.1 气态生成焓不是实验数据,是和俄罗斯人使用超级计算机使用同样 W1U 方法得到的结果。2个结果,拥有 0.18 kcal 的误差,算比较吻合。

俄罗斯人的数据在这里:XXXXXXXXXXXXXXXXXXXX/doi/suppl/10.1021/jp404484q/suppl_file/jp404484q_si_001.pdf

 W1 理论:XXXXXXXXXXXXXXXXXXXXXXXXX/doi/abs/10.1063/1.479454

W1U.png


总结

平均误差稍高于推广的 0.3 kcal 是因为 H-N=O 的气态生成焓和实验结果拥有 1.9kcal 的错误,个人认为这个可能是实验误差。

相对而言,该方法精度特别高。N2H4 的误差只有 0.37 kcal 的误差,在 B3LYP 6-31g(d) 的计算结果里,N2H4 的误差接近 10kcal,而是用同样方法带入到 1999 rice et al 方法后,误差更高。N2H4 在 B3LYP 6-311g+(2d,p) 精度的计算的结果和实验只有 2kcal 差距。

计算时间: W1> G2, CBS-APNO > G4 > G3, G2MP2 > CBS-QB3 > CBS-4M

精度:W1 > CBS-APNO > CBS-QB3, G4 > G3 > G2, G2MP2, CBS-4M

[修改于 6年4个月前 - 2018/08/22 12:30:13]

来自:数理化 / 化学
2
 
已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也

想参与大家的讨论?现在就 登录 或者 注册

所属专业
上级专业
同级专业
天泽
进士 学者 机友 笔友
文章
9
回复
22
学术分
1
2018/07/30注册,3年0个月前活动
暂无简介
主体类型:个人
所属领域:无
认证方式:身份证号
IP归属地:美国
文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
支持的图片格式:jpg, jpeg, png
插入公式
评论控制
加载中...
文号:{{pid}}
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}