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

计算化学公社

 找回密码
 现在注册!
查看: 1931|回复: 3

[量化理论] 高性价比热力学组合方法G4(MP2)-6X的计算方法

[复制链接]

2万

帖子

25

威望

2万

eV
积分
45687

管理员

公社社长

发表于 2018-5-27 00:04:18 | 显示全部楼层 |阅读模式
高性价比热力学组合方法G4(MP2)-6X的计算方法

文/Sobereva @北京科音
First release: 2018-May-26  Last update: 2019-Jul-30

在2011年,J. Chem. Theory Comput., 7, 112-120中,几个非Gn系列方法的官方人员提出了G4(MP2)-6X热力学组合方法,号称耗时与G4(MP2)相仿佛,而精度很接近昂贵得多的G4。其相对于G4(MP2)最主要变化是把MP2换成了SCS-MP2,把CCSD(T)的CCSD和(T)部分的相关能乘了系数,把优化和振动分析的泛函从B3LYP改为BMK。此方法竟然一直到当下最新的G16 B.01里都仍然没有被加入,好在文章的补充材料里给出了基于Gaussian做此方法计算的Perl脚本,使用很简单,这里简单说一下用法。原作者给的脚本只能输出H(0)和H(T),没法输出指定温度下的自由能和内能,因此我对脚本进行了一些修改使之能够输出。

首先需要编辑模板.gjf文件,内容如下,也可以直接从这里下: G4MP2_6x.gjf (686 Bytes, 下载次数: 24)

评分

参与人数 8eV +50 收起 理由
aqhuangry + 5 赞!
Shine剪水 + 5 谢谢
wangyj + 5 谢谢
zsu007 + 15 赞!
978142355 + 5 好物!
captain + 5 谢谢
一颗赛艇 + 5 GJ!不知什么时候能出个评测
我本是个娃娃 + 5 大博士赛高

查看全部评分

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

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

81

帖子

1

威望

271

eV
积分
372

Level 3 能力者

发表于 2019-7-25 15:35:39 | 显示全部楼层
谢谢sob老师的分享 有一问题请教 我的小分子液相和气相优化结构有较大的差异 想问能否用热力学的方式来算?
我的想法是这样

%chk=mol.chk
# BMK/6-31+G(2df,p) Opt SCRF=SMD

A molecule G4(MP2)-6X calculation

0 1
C                  0.00000000    0.00000000   -0.56221066
H                  0.00000000   -0.92444767   -1.10110537
H                 -0.00000000    0.92444767   -1.10110537
O                  0.00000000    0.00000000    0.69618930

--Link1--
%chk=mol.chk
# Geom=AllCheck Guess=Read BMK/6-31+G(2df,p) Freq SCRF=SMD

--Link1--
%chk=mol.chk
# Geom=AllCheck Guess=Read CCSD(T,FrzG4)/GTBas1

--Link1--
%chk=mol.chk
# Geom=AllCheck Guess=Read MP2(FrzG4)/GTMP2LargeXP

--Link1--
%chk=mol.chk
# Geom=AllCheck Guess=Read HF/GFHFB3

--Link1--
%chk=mol.chk
# Geom=AllCheck Guess=Read HF/GFHFB4

--Link1--
%chk=mol.chk
# Geom=AllCheck Guess=Read M052X/6-31G*

--Link1--
%chk=mol.chk
# Geom=AllCheck Guess=Read M052X/6-31G* SCRF=SMD


最后加上溶解自由能

2万

帖子

25

威望

2万

eV
积分
45687

管理员

公社社长

 楼主| 发表于 2019-7-25 21:18:58 | 显示全部楼层
wuzhiyi 发表于 2019-7-25 15:35
谢谢sob老师的分享 有一问题请教 我的小分子液相和气相优化结构有较大的差异 想问能否用热力学的方式来算? ...

不是“热力学的方式”而是“热力学组合方法的方式”

这么做从实际角度上可行,但原理上不十分严格,毕竟热力学组合方法里的参数是对于气相优化的情况拟合的
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

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

81

帖子

1

威望

271

eV
积分
372

Level 3 能力者

发表于 2019-7-30 16:02:02 | 显示全部楼层
本帖最后由 wuzhiyi 于 2019-7-30 17:05 编辑

https://pubs.acs.org/doi/10.1021/acs.jctc.9b00449
这里好像有人提出了一个将G4(MP2)-6X拓展到第五行以内的所有的主族元素
作者说在有些测试集表现比原本的G4(MP2)-6X好
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2019-9-16 09:03 , Processed in 0.166569 second(s), 27 queries .

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