近年来，随着量子化学理论的不断完善和计算机技术与计算化学软件的发展，计算化学已逐渐与实验化学、理论化学一起并称化学科学的三个重要支柱，在材料、催化、生物和医药等领域都已获得广泛应用。许多化学教育工作者都希望借助计算化学方法来解决化学实验教学中遇到的各种问题，因此使学生学习并掌握一定计算化学实验技能变得愈发必要与迫切。目前，国内很多高校纷纷开设了计算化学实验课程^[1–3]，设计了一系列面向高年级本科生的计算化学实验，这些教学实践不仅有助于学生了解和掌握计算化学的基本原理、方法和应用，还有助于学生加深对一些复杂抽象的化学理论知识的理解。

计算化学实验是在计算机上开展的以分子模拟为工具，通过“理论、计算、实验”协同作用来解决各种科学问题的实验。计算化学实验具有教学成本低和实验安全性高的优点，可以帮助学生直观、深入地学习与掌握微观抽象的物理化学基本概念。“电动势的测定”实验为南开大学化学国家级实验教学示范中心面向化学、材料及环境等学科专业本科生开设的一个经典的物理化学实验，采用补偿法(对消法)原理^[4]线下开展实验。借助教育信息化的机遇，通过在Cu-Zn原电池电动势的测定物理化学实验中应用信息化技术手段，以计算机作为实验工具，学生通过操作计算机上的Gaussian 16 ^[5]软件，结合物理化学基本理论自主完成实验。本实验可实现信息化技术与实验课程教学高度融合，突破时间与空间的限制，有效帮助学生加深对热力学循环、吉布斯自由能、电极电势、能斯特方程和电池电动势等基本概念的理解，对学生后续的课程学习及日后从事科研工作都有所帮助，为高端技术人才培养提供基础支撑。

(1) 掌握Gaussian 16软件的基本使用方法，能用密度泛函方法进行理论计算，用基于密度的溶剂模型(Solvation Model Based on Density，SMD)对溶剂化自由能进行计算。

(2) 设计合理的热力学循环来求解热力学量。通过模拟得到的溶剂化自由能，结合升华能、电离能实验值，计算出Cu-Zn原电池的标准电动势。

(3) 加深对标准摩尔吉布斯自由能、电极电势、电池标准电动势等概念的理解。

Cu电极和Zn电极都具有较高的交换电流密度值，由它们构成的原电池可逆性比较好。按文献^[6]数据，Cu/Cu²⁺和Zn/Zn²⁺的交换电流密度均为0.2 A·m⁻²，两个反应相对较大的交换电流密度，意味着它们的可逆性都比较好，在有相对小的电流通过电极时，电极电位基本保持不变，因而Cu-Zn原电池的电动势比较稳定。本文中Cu-Zn原电池由热力学平衡态的电极构成，可视为可逆电池，其电动势的计算遵循能斯特方程^[7]：

(1) $E=E^{\Theta}-\frac{R T}{z F} \ln \left(\prod\limits_{\mathrm{B}} a_{\mathrm{B}}^{v_{\mathrm{B}}}\right)$

在所有组分活度都为1 mol·dm⁻³的标准状态下，E = E^Ө。电池的标准电动势为正、负极的标准电极电势之差，可利用电池反应的标准摩尔吉布斯自由能变与电池标准电动势之间的关系^[8]进行计算，其计算式为：

(2)$ \begin{aligned}E^{\Theta} & =E_{\mathrm{m}, \mathrm{c}, \mathrm{r}}^{\Theta}-E_{\mathrm{m}, \mathrm{a}, \mathrm{r}}^{\Theta}=\left(E_{\mathrm{m}, \mathrm{c}, \mathrm{abs}}^{\Theta}-E_{\mathrm{SHE}, \mathrm{abs}}^{\Theta}\right)-\left(E_{\mathrm{m}, \mathrm{a}, \mathrm{abs}}^{\Theta}-E_{\mathrm{SHE}, \mathrm{abs}}^{\Theta}\right)=E_{\mathrm{m}, \mathrm{c}, \mathrm{abs}}^{\Theta}-E_{\mathrm{m}, \mathrm{a}, \mathrm{abs}}^{\Theta} \\& =\left[\left(-\frac{\Delta_{\mathrm{r}} G_{\mathrm{m}, \mathrm{c}}}{z F}\right)-\left(-\frac{\Delta_{\mathrm{r}} G_{\mathrm{m}, \mathrm{a}}^{\Theta}}{z F}\right)\right]=-\frac{\Delta_{\mathrm{r}} G_{\mathrm{m}, \mathrm{c}}^{\Theta}-\Delta_{\mathrm{r}} G_{\mathrm{m}, \mathrm{a}}^{\Theta}}{z F}=-\frac{\Delta_{\mathrm{r}} G_{\mathrm{m}}^{\Theta}}{z F}\end{aligned} $

式(2)中z为1 mol电极反应中得失电子的数目，F为法拉第常数，$ {E}_{\text{m, c, r}}^{\text{Ө}} $和$ {E}_{\text{m, a, r}}^{\text{Ө}} $分别表示正、负极相对于标准氢电极SHE的标准电极电势，$ {E}_{\text{m, c, abs}}^{\text{Ө}} $和$ {E}_{\text{m, a, abs}}^{\text{Ө}} $分别表示正、负极的标准绝对电极电势。

本实验的温度规定为25 ℃，$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c}}^{\text{Ө}} $和$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, a}}^{\text{Ө}} $分别为正、负极电极反应的标准摩尔吉布斯自由能变。

(I) $\text { 正极: } \mathrm{Cu}^{2+}(\mathrm{aq})+2 \mathrm{e}^{-} \rightarrow \mathrm{Cu}(\mathrm{s})$

(II) $\text { 负极: } \mathrm{Zn}(\mathrm{s})-2 \mathrm{e}^{-} \rightarrow \mathrm{Zn}^{2+}(\mathrm{aq})$

(III) $\text { 电池总反应: } \mathrm{Cu}^{2+}(\mathrm{aq})+\mathrm{Zn}(\mathrm{s}) \rightarrow \mathrm{Cu}(\mathrm{s})+\mathrm{Zn}^{2+}(\mathrm{aq})$

由式(2)可知，欲求Cu-Zn原电池的标准电动势$ {{E}}^{\text{Ө}} $，只需求得电极反应(Ⅰ)和(Ⅱ)的$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c}}^{\text{Ө}} $和$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, a}}^{\text{Ө}} $。以$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c}}^{\text{Ө}} $的计算为例($ {{-}{\Delta }}_{\text{r}}{{G}}_{\text{m, a}}^{\text{Ө}} $的计算方法同理)，设计了一个热力学循环，如图 1所示。

在图 1中Cu电极一侧，过程(1)表示从Cu(g)转化为Cu(s)的过程。Cu(s)的升华吉布斯自由能$ {{\Delta }}_{\text{sub}}{G}\left(\text{Cu}\right) $可通过Δ_subG = $ {{\Delta }_{\mathrm{f}}{H}}_{\text{gas}}^{\text{Ө}} $– T($ {{S}}_{\text{g}\text{, }\text{1bar}}^{\text{Ө}} $− $ {{S}}_{\text{s}\text{, }\text{1bar}}^{\text{Ө}} $)计算得到，其中$ {{\Delta }_{\mathrm{f}}{H}}_{\text{gas}}^{\text{Ө}} $即为Cu(g)的标准生成焓，$ {{S}}_{\text{g}\text{, }\text{1bar}}^{\text{Ө}} $与$ {{S}}_{\text{s}\text{, }\text{1bar}}^{\text{Ө}} $^[9]分别为Cu(g)和Cu(s)的标准熵。则过程(1)的吉布斯自由能变为$ {{-}{\Delta }}_{\text{sub}}{G} $。过程(2)为Cu²⁺(g)转化为Cu(g)的过程。此过程的吉布斯自由能变值即为Cu(g)的电离能相反值$ {-}{{\Delta }}_{\text{ion}}{G}\left(\text{Cu}\right) $^[10]。过程(3)为Cu²⁺(aq)转变为Cu²⁺(g)的过程。此过程的吉布斯自由能变为Cu²⁺(g)的溶剂化自由能$ {\Delta }{{G}}^{\text{*}} $(Cu²⁺)的相反值$ {-}{\Delta }{{G}}^{\text{*}}\left({\text{Cu}}^{\text{2+}}\right) $^[11]。

可得：

(3)$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c}}^{\text{Ө}}\text{=}{-}\text{(}{{\Delta }}_{\text{sub}}{G}\left(\text{Cu}\right)\text{+}{{\Delta }}_{\text{ion}}{G}\left(\text{Cu}\right){+ \Delta }{{G}}^{\text{*}}\text{(}{\text{Cu}}^{2+}\text{))} $

结合式(2)可知，通过$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c}}^{\text{Ө}} $来求解电池的标准电动势与将正、负极的标准电极电势相减来计算电池标准电动势的方法是等价的。

根据图 1中的热力学循环，将实验得到的热力学数据Δ_subG、Δ_ionG和$ {\Delta }{{G}}_{\text{exp}}^{\text{*}} $按照式(3)进行计算，得到Cu、Zn电极发生还原反应时的标准摩尔吉布斯自由能变$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c, exp}}^{\text{Ө}} $和$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, a, exp}}^{\text{Ө}} $ (表 1)。

利用表 1中的$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c, exp}}^{\text{Ө}} $与$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, a, exp}}^{\text{Ө}} $，根据式(2)可求出依据热力学循环计算的Cu-Zn原电池的标准电动势为：

${E}_{\text{exp}}^{\text{Ө}}\text{=}\text{}{-}\frac{{{\Delta }}_{\text{r}}{{G}}_{\text{m, exp}}^{\text{Ө}}}{zF}\text{}\text{=}\text{}{-}\frac{{{\Delta }}_{\text{r}}{{G}}_{\text{m, c, exp}}^{\text{Ө}}{-}\text{}{{\Delta }}_{\text{r}}{{G}}_{\text{m, a, exp}}^{\text{Ө}}}{zF}= 1.0958\text{ V }$

查文献^[12]中的标准电极电势表得到Cu、Zn电极在25 ℃时标准电极电势分别为$ {E}_{{\text{Cu}}^{\text{2+}}\text{/Cu}}^{\text{Ө}}\text{}\text{=}\text{}\text{0.34 V} $和$ {E}_{{\text{Zn}}^{\text{2+}}\text{/Zn}}^{\text{Ө}}\text{=}\text{}{-}\text{0.7626 V} $，故查表得到的Cu-Zn原电池标准电动势为：

$ {E}_{\text{ref}}^{\text{Ө}}\text{}\text{=}\text{}{E}_{{\text{Cu}}^{\text{2+}}\text{/Cu}}^{\text{Ө}}{-}\text{}{E}_{{\text{Zn}}^{\text{2+}}\text{/Zn}}^{\text{Ө}}\text{}\text{=}\text{}\text{1.1026 V} $

通过上述计算，可知依据热力学循环计算的Cu-Zn原电池标准电动势的相对误差为：

$ {RE}\text{}\text{=}\text{}\frac{{E}_{\text{exp}\text{}}^{\text{Ө}}{-}{\text{}\text{E}}_{\text{ref}}^{\text{Ө}}}{{E}_{\text{ref}}^{\text{Ө}}}\text{}\text{×}\text{}\text{100%}\text{}\text{=}\text{}{-}\text{0.62%} $

对比文献^[9–11]中实验测得的热力学数据实验值(表 1)计算得到的Cu-Zn原电池的标准电动势值(1.0958 V)与文献^[12]中查到的Cu-Zn原电池的标准电动势的值可知，相对误差仅为−0.62%，表明通过热力学循环求解Cu-Zn原电池标准电动势的方法具有正确性。

Gaussian软件是目前使用较广泛的一款量子化学计算软件，其最新版本是Gaussian16。Gaussian软件通常包括用于计算的Gaussian和用于图形界面处理的GaussView两个组件。其中Gaussian软件涵盖了从头算方法、密度泛函方法和半经验方法等多种计算方法^[13]。密度泛函理论基于Hohenberg-Kohn定理，将3N维多电子波函数问题简化为3维电子密度问题，用电子密度分布函数求解薛定谔方程，从而得到系统总能量及相关性质。Gaussian中，可通过对优化好的结构进行频率分析，进而计算分子的内能、焓、吉布斯自由能等热力学量。Gaussian中吉布斯自由能的计算过程为G = E_ele + G_corr，其中E_ele为单点能，即分子的电子能量，包括电子的动能、电子与原子核间的势能、电子之间的相互作用能以及原子核之间的相互作用能；G_corr则是温度与压力的影响所带来的吉布斯自由能修正项。由于实验上许多构成电池的分子与离子的溶剂化自由能较难直接测得，故本实验采用Gaussian计算方法得到离子的溶剂化自由能，进而计算出Cu-Zn原电池的标准电动势。

实验中Cu²⁺的溶剂化自由能$ {\Delta }{{G}}_{\text{cal}}^{\text{*}} $(Cu²⁺)是通过Gaussian计算的Cu²⁺溶液相的吉布斯自由能(G_aq)与气相吉布斯自由能(G_g)相减得到。因为气相和溶液相的G_corr通常相差不大，故不需要做溶剂模型下的频率分析来获得$ {\Delta }{{G}}_{\text{cal}}^{\text{*}} $(Cu²⁺)，而只做单点能计算，用溶液相单点能$ {E}_{\text{ele(aq)}} $与气相单点能$ {E}_{\text{ele(g)}} $相减即得到$ {\Delta }{{G}}_{\text{cal}}^{\text{*}} $(Cu²⁺)。本实验中，使用Gaussian16计算Cu-Zn原电池标准电动势的过程如图 2所示。

Gaussian 16软件，GaussView 5.0软件，DELL台式电脑(处理器：i5-6500；核心数：4；运行内存：32GB)。

构建初始结构，打开GaussView，点击Cu原子，保存为输入文件(*.gjf)，采用M06-2X/6-31G(d)密度泛函方法，对Cu²⁺进行单点能的计算。M06-2X属于Truhlar课题组开发的明尼苏达泛函系列，是采用广义近似梯度的杂化密度泛函之一^[14]。6-31G(d)是基组的名称，表示内层轨道用6个Gauss型函数(GTO)拟合1个Slater型函数(STO)，这1个STO拟合1个原子轨道，价层轨道分为内(I)、外(O)轨道，内层轨道用3个GTO拟合1个STO(I)，外层轨道用1个GTO拟合一个STO(O)，然后1个STO(I)和STO(O)拟合一个价层原子轨道；在6-31G基组的基础上给重原子添加了d轨道极化函数。

编辑输入文件，Gaussian中默认的计算项目即单点能计算，因此计算单点能时关键词sp可不写。Gaussian默认状态为气相的计算，不需额外增加关键词，溶液相的计算则需要添加溶剂模型。溶剂模型包括隐式溶剂模型和显式溶剂模型，前者不考虑溶剂分子的实际结构，只把溶剂视为具有一定介电常数的连续介质，后者则是计算时在体系中直接增加溶剂分子，考虑溶质与溶剂分子间的相互作用以及溶剂分子的不同排布方式，这会导致计算量大幅增加。本实验采用的SMD溶剂模型^[15]，是由Truhlar及其合作者发展的一种参数化的基于SCRF的隐式溶剂化模型，该模型已经用2821个溶剂化数据的训练集进行了参数化，包括112个水性离子溶剂化自由能。Gaussian中使用scrf关键词来调用溶剂模型，因此在计算执行命令行关键词中需添加“scrf=SMD”，因SMD模型默认的溶剂为水，故“solvent=water”可不写。在电荷与多重度行，Cu²⁺的电荷为+2，有一个单电子，因此电荷和多重度设置为“2 2”。Cu²⁺的坐标直接设定为“0.0 0.0 0.0”，气相单点能计算输入文件如表 2所示。溶液相单点能计算的输入文件只需将表 2中执行命令行关键词中加上“scrf=SMD”。

输入文件编辑完成后，点击GaussView中“Calculation”下拉菜单“Gauss Calculation Setup”命令进入设置Gaussian计算选项，点击“Submit”选项，提交计算任务，也可通过在Gaussian软件中打开已建立的输入文件，进行计算任务提交。计算完成后会得到单点能计算的输出文件(*.log，若用Gaussian软件直接提交计算，则输出文件格式为*.out)。用文本编辑器打开输出文件，在“HF=”处读取Cu²⁺和Zn²⁺气相和溶液相的单点能数据。

由4.1小节中提取的Cu²⁺和Zn²⁺气相和溶液相单点能数值，计算得到它们的溶剂化自由能$ {\Delta }{{G}}_{\text{cal}}^{\text{*}} $。因为SMD溶剂化自由能是指从298.15 K的理想气体到298.15 K理想溶液的自由能变化，并且在气相和溶液相中使用相同的浓度(1 mol·L⁻¹)。但本实验所选取的状态是标准状态，为气相$ {{p}}^{\text{Ө}} $下溶质到1 mol·L⁻¹溶液相状态过程的自由能的变化，因此需要加上修正项，修正项的计算为：

$ {{\Delta }}_{\text{aq}}{{G}}_{\text{corr}\text{}}\text{=}\text{}{RT}\text{ln}\frac{{RT}}{{{p}}^{\text{Ө}}}\text{}\text{=}\text{}\text{1.89 kcal}\cdot{\text{mol}}^{{-}\text{1}}\text{}\text{=}\text{}\text{0.0820 eV} $

所以，

$ {\Delta }{{G}}_{\text{cal}}^{\text{*}}\text{}\text{=}\text{}{{G}}_{\text{aq}}{-}\text{}{{G}}_{\text{g}}\text{}\text{=}\text{}{E}_{\text{ele}\left(\text{aq}\right)}{-}\text{}{E}_{\text{ele}\left(\text{g}\right)}\text{}\text{+}\text{}{{\Delta }}_{\text{aq}}{{G}}_{\text{corr}} $

将Gaussian计算得到的溶剂化自由能与查到的原子升华自由能^[9]和电离自由能变^[10]按式(3)进行计算，可以得到Cu电极发生还原反应时的标准吉布斯自由能$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c}}^{\text{Ө}} $。同理，Zn负极侧$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, a}}^{\text{Ө}} $的计算过程如同$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c}}^{\text{Ө}} $，进而计算得到Cu-Zn原电池反应的标准摩尔吉布斯自由能变化$ {{\Delta }}_{\text{r}}{{G}}_{\text{m}}^{\text{Ө}} $，再由式(2)可计算出Cu-Zn原电池的标准电动势。

由于Cu²⁺与Zn²⁺所带电荷量相同，离子半径相近，实验测得二者的溶剂化自由能相近，这与本实验的Gaussian计算结果($ {\Delta }{{G}}_{\text{cal}}^{\text{*}} $，见表 3)相一致。两种离子溶剂化自由能的Gaussian计算值$ {\Delta }{{G}}_{\text{cal}}^{\text{*}} $与实验测得的溶剂化自由能值($ {\Delta }{{G}}_{\text{exp}}^{\text{*}} $，见表 1)相比有所偏差，可能是由于Gaussian计算中选用的SMD溶剂模型，忽略了离子的真实溶剂化结构造成的。

根据表 3中Gaussian计算的溶剂化自由能，可计算得到Cu-Zn原电池电极反应的$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, c}}^{\text{Ө}} $和$ {{\Delta }}_{\text{r}}{{G}}_{\text{m, a}}^{\text{Ө}} $值，进而得到电池总反应的$ {{\Delta }}_{\text{r}}{{G}}_{\text{m}}^{\text{Ө}} $，通过式(2)可得到其标准电动势为：

$ {E}_{\text{cal}}^{\text{Ө}}\text{}\text{=}\text{}{-}\frac{{{\Delta }}_{\text{r}}{{G}}_{\text{m}}^{\text{Ө}}}{zF}\text{}\text{=}\text{}\text{1.4447 V} $

$ {E}_{\text{cal}}^{\text{Ө}} $与文献^[11]得到的$ {E}_{\text{ref}}^{\text{Ө}} $的相对误差为：

$ {RE}\text{}\text{=}\text{}\frac{{E}_{\text{cal}}^{\text{Ө}}{-}\text{}{E}_{\text{ref}}^{\text{Ө}}}{{E}_{\text{ref}}^{\text{Ө}}}\text{}\text{×}\text{}\text{100%}\text{}\text{=}\text{}\frac{\text{1.4447}\text{}{-}\text{}\text{1.1026}}{\text{1.1026}}\text{}\text{×}\text{}\text{100%}\text{}\text{=}\text{}\text{31.03%} $

可见$ {E}_{\text{cal}}^{\text{Ө}} $与$ {E}_{\text{ref}}^{\text{Ө}} $之间的相对误差较大，一方面是因为在构成Cu-Zn原电池的Cu、Zn电极的固/液界面处会分别形成带相反电荷的双电层。在Cu电极一侧会形成电极表面带正电、电极附近溶液相带负电的双电层，在Zn电极一侧则会形成电极表面带负电、其附近溶液相带正电的双电层。因而电化学反应的本质与固/液界面的场强密切相关，且固/液界面的场强会对体系的过电位产生影响，而本文的计算中没有考虑过电位，这会造成一些计算误差。另一方面，实际电极过程，比如电极反应(Ⅰ)，包括了Cu²⁺去溶剂化过程和跨界面的电子(2e)转移过程，在图 1的热力学循环中，这些过程的自由能都按照零来考虑的。在实际测量中，Cu、Zn电极会有不同程度的极化，产生的过电位对实验结果有所影响。而固/液界面的场强、电子逸出功以及金属离子由液相进入固相时的去溶剂化效应都会影响体系的过电位，因此会造成一些计算误差。此外，本实验中采用SMD隐式溶剂模型计算离子的溶剂化能，忽略了离子的溶剂化结构，因而也会造成一定的计算误差。Cu²⁺和Zn²⁺的电荷量较大且较为集中，Cu²⁺的基态还存在单电子，易与溶剂水分子之间产生强静电相互作用。因此，条件允许的情况下可选用考虑了离子的溶剂化结构的显式溶剂模型来计算离子的溶剂化自由能，但此模型虽然可提高计算的准确性，却必须要优化大量的溶剂分子的结构，致使计算量会大幅增加。考虑到计算成本，本实验选用了一种简便高效且计算量较小的隐式溶剂模型计算了离子的溶剂化自由能，此方法还可用于计算其他离子的溶剂化自由能，结合热力学循环，可求得不同电池的标准电动势。

根据文献^[16]，考虑到Cu²⁺添加一些显式的溶剂水分子，再对Cu²⁺的水合团簇使用SMD模型计算团簇的溶剂化自由能，利用热力学循环便能求出Cu²⁺的溶剂化自由能。此方法称为团簇连续化溶剂模型，是一种结合了显隐式溶剂模型的混合溶剂模型，考虑了溶质分子真实的溶剂化结构，使溶剂化自由能的计算结果更精确。因此，在使用量化软件进行计算时，不仅需要正确的理论方法，还需要选用合适的模型。

1) 根据计算结果得知本实验的误差主要来自于离子溶剂化自由能的计算过程，请思考原因。尝试更换计算方法，看能否解决这个问题。

2) 请阅读文献^[16]，尝试用团簇连续化溶剂模型计算Cu²⁺和Zn²⁺溶剂化自由能，并与实验值进行比较。

本实验以计算机作为实验工具，发展了经典的“电动势的测定”这一物理化学实验。在实验中，学生可通过操作Gaussian 16软件线上开展Cu-Zn原电池标准电动势的测定实验，可与线下实际操作互补结合。学生通过自主设计热力学循环，计算出原电池反应的标准摩尔吉布斯自由能变，进而得到Cu-Zn原电池的标准电动势。此实验可实现现代化信息化技术与实验教学有效融合，帮助学生加深对标准摩尔吉布斯自由能、电极电势以及电池电动势等基本概念的理解，引导学生进行发散思维，感悟计算方法中严密的逻辑关系，有利于提高学生分析解决问题能力、科学研究能力与动手实践能力。