计算化学公社

标题: Gaussian09 CBS-QB3计算输出文件的能量计算公式 [打印本页]

作者
Author: 清微 时间: 2017-12-27 22:57
标题: Gaussian09 CBS-QB3计算输出文件的能量计算公式
本帖最后由清微于 2017-12-27 23:04 编辑

高斯的CBS计算输出文件中最后几项：

(, 下载次数 Times of downloads: 31)

E(Thermal)应该等于Etrans+Erot+Evib，从而CBS-QB3 Energy= CBS-QB3 (0 K) + E(Thermal)，但实际上不等，难道后一项也要乘以校正因子？

作者
Author: sobereva 时间: 2017-12-28 03:26
CBS-QB3(0K)已经包含了ZPE，E(thermal)也包含ZPE，加起来显然不是T下的内能

0.019855-0.016991-39.069832=-39.066968

作者
Author: 清微 时间: 2017-12-28 12:57
这是我分析得到的CBS-QB3各项能量值：
CBS-QB3计算共有5步：I.B3LYP/CBSB7 opt    II.B3LYP/CBSB7 freq    III.CCSD(T)/6-31+G(d')    IV.MP4SDQ/CBSB4    V.MP2/CBSB3 CBSExtrap=(NMin=10,MinPop)
左边是CBS计算的总结部分，右边标注了对应的值出现在第几步
E(SCF)=HF                                        V
DE(MP2)=MP2-HF                               V
DE(CBS)=E2(CBS)-DE(MP2)             V
DE(MP34)=MP4SDQ-MP2                   IV
DE(CCSD)=CCSD(T)-MP4SDQ          III
DE(Int)=CBS-Int-E2(CBS)                   V
DE(Empirical)可能与OIii这一项有关    V

最终，Eelec=E(SCF)+DE(MP2)+DE(CBS)+DE(MP34)+DE(CCSD)+DE(Int)+DE(Empirical)

CBS-QB3 (0 K)=Eelec+E(ZPE)
CBS-QB3 Energy=Eelec+E(Thermal)
CBS-QB3 Enthalpy=CBS-QB3 Energy+RT
CBS-QB3 Free Energy=CBS-QB3 Enthalpy-TS

问题1：这样的总结对么？
问题2：文献中给出DE(Empirical)计算公式（见附图），m，Eh分别是什么含义？后面两项在输出文件哪个部分？
问题3：这种组合方法和博文“谈谈能量的基组外推”中，分步计算再外推有什么区别？

作者
Author: GMW 时间: 2023-4-4 21:54

清微发表于 2017-12-28 12:57
这是我分析得到的CBS-QB3各项能量值：
CBS-QB3计算共有5步：I.B3LYP/CBSB7 opt II.B3LYP/CBSB7 freq ...

你好我在CBS-QB3计算吉布斯自由能的时候遇到一些问题
问题1第一步是CBSB7是怎么选上的啊？Gaussion View5.0.9里没找到这个基组
问题2第III-V中Job Type是计算的能量吗？

作者
Author: TYCHEN 时间: 2023-9-18 21:55

GMW 发表于 2023-4-4 21:54
你好我在CBS-QB3计算吉布斯自由能的时候遇到一些问题
问题1第一步是CBSB7是怎么选上的啊？Gaussion Vie ...

我也想问，CBSB7还有CBS-QB3到底指的是啥啊？具体是啥基组

作者
Author: TYCHEN 时间: 2023-9-18 21:57

GMW 发表于 2023-4-4 21:54
你好我在CBS-QB3计算吉布斯自由能的时候遇到一些问题
问题1第一步是CBSB7是怎么选上的啊？Gaussion Vie ...

http://sobereva.com/172 你看看这个

作者
Author: sobereva 时间: 2023-9-18 23:53

TYCHEN 发表于 2023-9-18 21:55
我也想问，CBSB7还有CBS-QB3到底指的是啥啊？具体是啥基组

这种事自己看原文便知
没有叫Gaussion View，正经用Gaussian的内行人也没人用GaussView产生输入文件，问怎么用这种菜鸟界面产生输入文件毫无意义

欢迎光临计算化学公社 (http://bbs.keinsci.com/)