（一）引言：

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

（二）計算過程：

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

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

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

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

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

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

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

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