(一)引言：

? 在使用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ì)算過程：

1.計(jì)算設(shè)置：

█ 軟件版本Gaussian 16, GaussView 6。使用方法 # B3LYP/6-31G(d) opt freq

█ 先使用GaussView 6構(gòu)建甲烷模型，然后按照上述方法設(shè)置，保存輸入文件名為methane.gjf，chk文件名為methane.chk，此處省略計(jì)算過程。

2.Windows (Linux)下使用freqchk命令產(chǎn)生不同溫度和壓力下的熱力學(xué)數(shù)值

█ 首先打開Gaussian 16的界面，然后點(diǎn)擊界面上的Utilities功能，選中freqchk命令，然后到chk文件保存目錄，打開methane.chk文件。得到如下界面：

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

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

█ 設(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)上述過程。并且把結(jié)果輸入到correct.txt文件中，如下：

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

總結(jié)

以上是生活随笔為你收集整理的高斯用服务器能算chk文件么,教程（二十四）使用freqchk命令进行热力学校正的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。