计算化学是当代化学领域的一门重要学科。1998年，诺贝尔化学奖授予沃尔特·库恩和约翰·波普尔，以表彰他们在计算化学领域做出的重大贡献。2013年马丁·卡普拉斯、迈克尔·莱维特及亚利耶·瓦谢尔因他们在开发多尺度复杂化学系统模型方面所做的贡献而获得诺贝尔化学奖，标志着计算化学的重要作用愈加为人们所重视。计算化学主要基于量子化学和统计力学等理论来研究化学问题，这两部分内容同时也是物理化学基础课教学中的重点和难点之一。计算化学的研究结果不仅能在分子水平上对实验化学的结果作出合理的解释，还能对分子的化学性质和化学反应的活性、选择性以及产物分布进行预测。随着计算化学的发展、计算软件的普及和计算机硬件水平的提高，许多学生对计算化学产生了浓厚的兴趣。针对这一现象，国内化学教育工作者进行了积极的应对。很多高等院校纷纷开设了计算化学课程^[1-3]，一些从事物理化学教学的教师也有意识地在物理化学基础课中引入计算化学内容，并设计了“甲醛分子的结构和性质的计算化学研究”、“H3反应势能面的构建”等计算化学实验^{[4, 5]}。笔者在物理化学教学实践中也进行了尝试，将计算化学的原理、方法和应用实例引入物理化学教学内容中，并在此基础上设计了一个面向高年级本科生的计算化学实验--${{\text{H}}_{\text{2}}}{\text{(g) + }}\frac{1}{2}{{\text{O}}_{\text{2}}}{\text{(g) }} \to {{\text{H}}_{\text{2}}}{\text{O}}\left( {\text{g}} \right)$反应热力学函数的高水平量化计算研究。该实验具有以下3个特点：(1)使用最新版本的Gaussian 09软件进行高水平量化计算，使学生直接接触到计算化学的前沿领域；(2)通过“构建初始构型→编制输入文件→优化几何构型→计算振动频率→计算热力学函数变化值”的链条式一体化设计思路，使学生通过对一个简单反应的计算了解和掌握Gaussian软件应用的多个重要方面，具有较高的实用性；(3)选取的H₂、O₂和H₂O具有单键、双键两种共价键形式，且分子结构简单、计算量较小、耗时较短，能够保证在规定课时内完成实验内容。

(1) 初步掌握Gaussian 09化学计算软件及配套软件GaussView和UltraEdit的基本使用方法。

(2) 学习使用密度泛函(DFT)方法和热力学组合方法进行构型优化、频率计算和能量计算的基本过程。

Gaussian软件是目前使用极为广泛的一款量子化学计算软件。Gaussian软件最初由诺贝尔奖获得者波普尔等人于1970年开发，后经计算化学家不断补充和完善，该程序历经诸多版本，如Gaussian 94、Gaussian 98、Gaussian 03，目前最新的版本是Gaussian 09。Gaussian一般有两种版本，一种应用于windows操作系统，一种应用于Linux操作系统。一套Gaussian软件一般包括用于计算的Gaussian和用于图形界面处理的GaussView两个组件。Gaussian软件包含了从头算方法、密度泛函方法、半经验方法、分子力学方法等多种计算方法，具有操作简单、功能强大的特点。使用Gaussian 09进行计算的基本过程如下：

(1) 分子构型的构建。通过GaussView或Chem3D软件构建化合物的初始构型。

(2) 计算方法的选取。根据现有的计算条件、模型的大小以及所要解决的问题，选择可行的计算方法，采用UltraEdit程序编制输入文件。

(3) 量化计算。启动Gaussian 09，使用“File”下拉菜单“Open”命令将编制好的输入文件代入Gaussian 09进行构型优化、振动频率计算等，获得分子的平衡几何构型、稳定能态(全部为振动实频而无虚频)、化学反应过渡态(有且只有一个虚频)。

(4) 计算结果的分析和整理。对计算结果进行分析和整理，一般包括构型描述、能量分析、轨道组成分析、电荷和成键分析等，提取有用的信息。

Gaussian 09软件，GaussView 5.0软件，UltraEdit软件，CYLview v1.0软件，联想台式电脑(i5-4460 4G 500G)。

本实验采用Gaussian 09程序中的两种方法：B3LYP/6-31G (d, p)方法和G4方法分别计算${{\text{H}}_{\text{2}}}{\text{(g) + }}\frac{1}{2}{{\text{O}}_{\text{2}}}{\text{(g) }} \to {{\text{H}}_{\text{2}}}{\text{O}}\left( {\text{g}} \right)$反应的热力学函数变化值。

B3LYP是采用广义近似梯度的杂化密度泛函之一。6-31G (d, p)是基组，表示内层轨道用6个高斯函数(GTO)拟合1个Slater函数(STO)，然后用这个STO拟合1个原子轨道；价层轨道则分为内、外轨，内轨用3个GTO拟合1个STO，外轨用1个GTO拟合1个STO，然后用这2个STO拟合1个价层原子轨道。该实验需要用B3LYP/6-31G (d, p)方法对所计算的分子进行构型优化和频率分析，因此需要在计算执行路径行使用关键词opt和freq，其中opt即是对相应的分子进行构型优化，freq即是对相应的分子进行频率分析。

用B3LYP/6-31G (d, p)方法对H2O分子进行计算时的输入文件(通常以gjf为后缀名)和说明如表1所示。

同理，用B3LYP/6-31G (d, p)方法对H2和O2分子进行计算时的输入文件如表2所示。

G4方法是热力学组合方法之一，该类方法是专门用来高精度计算热力学量的方法。该类方法是一套过程的组合：包括中低级别优化构型、振动分析以获得热力学校正量、在高级别下做单点计算并进行适当经验校正或外推，整套过程自动完成，因此热力学组合方法不用指定基组和关键词。

用G4方法对H₂O分子进行计算时的输入文件及说明如表3所示。

同理，用G4方法对H₂、O₂分子进行计算时的输入文件如表4所示。

计算结束后用GaussView或UltraEdit软件打开*.out输出文件查看全部的计算结果。

(1) 在CYLview v1.0中打开*.out格式的输出文件，可以得到B3LYP/6-31G (d, p)方法优化出的H₂、O₂和H₂O分子的三维结构(图1)。

(2) 在UltraEdit软件中打开*.out格式的输出文件，选择热力学函数的理论预测值(表5)。

根据表5中的输出结果，可以计算${{\text{H}}_{\text{2}}}{\text{(g) + }}\frac{1}{2}{{\text{O}}_{\text{2}}}{\text{(g) }} \to {{\text{H}}_{\text{2}}}{\text{O}}\left( {\text{g}} \right)$反应在298.15 K、101325 Pa下的焓变为：

$\begin{array}{l}{\Delta _{\rm{r}}}{H_{\rm{m}}} = \sum {{v_{\rm{B}}}{H_{\rm{m}}}({\rm{B}})} \\ = {H_{\rm{m}}}({{\rm{H}}_{\rm{2}}}{\rm{O,g}}) - {H_{\rm{m}}}({{\rm{H}}_{{\rm{2,}}}}{\rm{g}}) - \frac{1}{2}{H_{\rm{m}}}({{\rm{O}}_{\rm{2}}}{\rm{,g}})\\{\rm{ = - 200573}}{\rm{.99079 - ( - 3058}}{\rm{.86766) - }}\frac{1}{2} \times \left( { - 394646.65810} \right)\\{\rm{ = - 191}}{\rm{.794(kJ}} \cdot {\rm{mo}}{{\rm{l}}^{{\rm{ - 1}}}}{\rm{)}}\end{array}$

类似地，可以计算出该反应的其他热力学函数的变化值分别如下：

${\Delta _{\text{r}}}{U_{\text{m}}}{\text{ =-190}}{\text{.555kJ}} \cdot {\text{mo}}{{\text{l}}^{{\text{-1}}}}{\text{, }}{\Delta _{\text{r}}}{G_{\text{m}}}{\text{ =-178}}{\text{.671kJ}} \cdot {\text{mo}}{{\text{l}}^{{\text{ - 1}}}}{\text{, }}{\Delta _{\text{r}}}{S_{\text{m}}}{\text{ = - 43}}{\text{.976J}} \cdot {\text{mo}}{{\text{l}}^{{\text{ - 1}}}} \cdot {{\text{K}}^{{\text{ - 1}}}}$

在UltraEdit软件中打开*.out格式的输出文件，选择热力学函数的理论预测值列于表6。

根据表6中的输出结果，可以计算${{\text{H}}_{\text{2}}}{\text{(g) + }}\frac{1}{2}{{\text{O}}_{\text{2}}}{\text{(g) }} \to {{\text{H}}_{\text{2}}}{\text{O}}\left( {\text{g}} \right)$反应在298.15 K、101325 Pa下的焓变为：

$\begin{array}{l}{\Delta _{\rm{r}}}{H_{\rm{m}}} = \sum {{v_{\rm{B}}}{H_{\rm{m}}}({\rm{B}})} \\ = {H_{\rm{m}}}({{\rm{H}}_{\rm{2}}}{\rm{O,g}}) - {H_{\rm{m}}}({{\rm{H}}_{{\rm{2,}}}}{\rm{g}}) - \frac{1}{2}{H_{\rm{m}}}({{\rm{O}}_{\rm{2}}}{\rm{,g}})\\ = - 200571.03973 - \left( { - 3057.95923} \right) - \frac{1}{2} \times \left( { - 394547.92880} \right)\\ = - 239.116{\rm{(kJ}} \cdot {\rm{mo}}{{\rm{l}}^{{\rm{ - 1}}}}{\rm{)}}\end{array}$

类似地，可以计算出该反应的其他热力学函数值分别如下：

${\Delta _{\text{r}}}{U_{\text{m}}}{\text{ =-}}237.877{\text{kJ}} \cdot {\text{mo}}{{\text{l}}^{{\text{-1}}}}{\text{, }}{\Delta _{\text{r}}}{G_{\text{m}}}{\text{ =-}}225.983{\text{kJ}} \cdot {\text{mo}}{{\text{l}}^{{\text{ - 1}}}}$

文献给出的量热实验测定结果^[6]为：298.15 K、101325 Pa下反应${{\text{H}}_{\text{2}}}{\text{(g) + }}\frac{1}{2}{{\text{O}}_{\text{2}}}{\text{(g) }} \to {{\text{H}}_{\text{2}}}{\text{O}}\left( {\text{g}} \right)$的△_rH_m=-241.818 kJ·mol^-1。由此可见，G4方法计算出的△_rm^H结果(-239.116 kJ·mol^-1)与实验测定值更接近。

需要说明的是，本实验的计算对象为真空的中分子，因此高水平计算的结果也有可能与实验测定值之间存在微小的差异。另一方面，量子化学计算的精度取决于所用的方法和基组，不同方法和基组得到的计算结果在某些情况下差异会很大，应结合实际体系和所具备的实验条件选择适合的方法和基组。

(1) 初始构型的构建以及对计算方法和基组的选择很大程度上决定了量化计算的准确性和耗时。那么构建合理的初始构型方法有哪些？如何选择适合的计算方法和基组？

(2) 试列举除热力学函数外，Gaussian程序计算还能给出哪些物理化学性质的信息？

本实验是面向大三下学期或大四上学期已学习过物理化学课程中量子化学和统计力学相关知识的学生开设的探索性实验。实验用时5课时，具体为：简要介绍量子化学计算的相关知识，1课时；介绍Gaussian 09、GaussView 5.0、UltraEdit软件的使用方法，如何编制输入文件以及如何查看输出文件，2课时；上机计算并完成数据的计算与整理，2课时。课后学生独立完成实验报告。

笔者在多年的物理化学教学中，一直坚持理论教学和实践教学相结合的教学方法。针对目前物理化学教学中学生普遍存在对量子化学和统计力学部分内容“复杂、难学、实用性低”的看法，笔者进行了在物理化学教学中引入计算化学内容的尝试，如有意识地介绍计算化学与化学其他学科相结合的最新研究进展，并将实验室的计算资源向学生开放，让学生通过计算化学实验了解如何运用物理化学基础知识解决科学问题。实践证明，开设“${{\text{H}}_{\text{2}}}{\text{(g) + }}\frac{1}{2}{{\text{O}}_{\text{2}}}{\text{(g) }} \to {{\text{H}}_{\text{2}}}{\text{O}}\left( {\text{g}} \right)$反应热力学函数的高水平量化计算研究”等计算化学实验，在物理化学教学内容中引入计算化学的原理、方法和应用实例，不但加深了学生对物理化学知识的理解，而且有利于锻炼学生将物理化学基本知识运用到具体科学研究中的能力，提高学生的整体素质。很多学生在进行过计算化学实验后，对物理化学产生了浓厚的兴趣，并在后来的毕业论文(设计)和研究生阶段选择了物理化学和计算化学方向。同时一些选择有机化学、无机化学、材料化学方向的学生也基于计算化学实验的经历，在毕业论文(设计)和研究生阶段有意识地将计算化学方法应用于自身的科研工作中，取得了很好的效果。