Gaussian09 CBS-QB3计算输出文件的能量计算公式

清微

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



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

sobereva

CBS-QB3(0K)已经包含了ZPE，E(thermal)也包含ZPE，加起来显然不是T下的内能

0.019855-0.016991-39.069832=-39.066968

清微

这是我分析得到的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：这种组合方法和博文“谈谈能量的基组外推”中，分步计算再外推有什么区别？

GMW

清微 发表于 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是计算的能量吗？

TYCHEN

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

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

TYCHEN

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

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

sobereva

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

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