（一）引言：

? 在使用Gaussian進(jìn)行計(jì)算時(shí)，我們根據(jù)所要研究體系的狀態(tài)設(shè)置溫度和壓力等條件，教程二十三介紹了如何設(shè)置體系的溫度和壓力。但在某些體系中，我們想要考察該體系在不同溫度和壓力下的熱力學(xué)量，應(yīng)該怎么辦呢？我們并不需要重新設(shè)定溫度和壓力重新進(jìn)行計(jì)算（重新進(jìn)行頻率計(jì)算需要花費(fèi)大量的計(jì)算時(shí)間）。本文以教程二十三中的甲烷計(jì)算為例，將介紹如何使用freqchk命令產(chǎn)生不同溫度和壓力下的熱力學(xué)數(shù)值。

（二）計(jì)算過(guò)程：

1.計(jì)算設(shè)置：
█ 軟件版本Gaussian 16, GaussView 6。使用方法 # B3LYP/6-31G(d) opt freq
█ 先使用GaussView 6構(gòu)建甲烷模型，然后按照上述方法設(shè)置，保存輸入文件名為methane.gjf，chk文件名為methane.chk，此處省略計(jì)算過(guò)程。

2.Windows (Linux)下使用freqchk命令產(chǎn)生不同溫度和壓力下的熱力學(xué)數(shù)值
█ 首先打開(kāi)Gaussian 16的界面，然后點(diǎn)擊界面上的Utilities功能，選中freqchk命令，然后到chk文件保存目錄，打開(kāi)methane.chk文件。得到如下界面：

█ 注意，Linux 系統(tǒng)下只要在命令行模式下輸入freqchk methane.chk 即可（需要在目標(biāo)chk文件所在目錄下進(jìn)行該命令，如果不在，則輸入絕對(duì)路徑效果一樣）。以下內(nèi)容Windows和Linux都一致。

█ 本文為了能夠更清晰地展示，此處把上圖中的信息以文本形式展現(xiàn)出來(lái)：

█ 設(shè)置完畢參數(shù)后，就會(huì)輸出該溫度和壓力下的熱力學(xué)校正值，如下圖：

█ 值得注意的是，此處并未給出該體系最終的熱力學(xué)數(shù)值。若要得到體系最終的熱力學(xué)數(shù)值，需要使用0K 時(shí)體系的電子能量（即教程二十三中提到的Electronic Energy或EE值）再加上這個(gè)校正值。比如用EE + Thermal Enthalpy Correction 就是該體系在500K，20.0 atm下的焓值；EE + Thermal Free Energy Correction 就是該體系在500K，20.0 atm下的吉布斯自由能值。

█ 最后，還有一個(gè)更為簡(jiǎn)便的命令組合，可以實(shí)現(xiàn)上述過(guò)程。并且把結(jié)果輸入到correct.txt文件中，如下：

freqchk  methane.chk   n   500   20.0    1    yn>  correct.txt