HCN与HNC异构化反应活化能的量化计算实验

doi:10.3866/PKU.DXHX201912058

摘要/Abstract

摘要：

结合计算化学教学内容，介绍一个面向高年级本科生的与化学动力学相关的计算化学实验。运用量子化学软件Gaussian 09及其可视化软件GaussView和文本编辑器UltraEdit，通过量子化学方法搜索$ \text{HCN}\rightleftharpoons \text{HNC} $异构化反应的过渡态。计算反应物（HCN）、过渡态（TS）和产物（HNC）的能量，从而求算反应的活化能。旨在通过本实验，使学生初步掌握运用量子化学计算方法求算化学反应活化能的思路和方法。

关键词: 计算化学实验, 量化计算, Gaussian 09软件, 活化能

Abstract:

On the basis of teaching content of computational chemistry course, we introduce a computational chemistry experiment to senior undergraduates. This experiment uses Gaussian 09, GaussView and UltraEdit softwares to search the transition state of $ \text{HCN}\rightleftharpoons \text{HNC} $ isomerization reaction, calculate the energies of reactant (HCN), transition state (TS) as well as the product (HNC), and consequently to obtain the reaction activation energy. Through this experiment, students will learn how to perform the energy calculation with the quantum chemistry calculation methods to obtain activation energy of a reaction.

Key words: Computational chemistry experiment, Quantum chemistry calculation, Gaussian 09 software, Activation energy

许秀芳. HCN与HNC异构化反应活化能的量化计算实验[J]. 大学化学, 2021, 36(2): 1912058.

Xiufang Xu. Quantum Chemistry Calculation Experiment on Activation Energy for Isomerization of HCN and HNC[J]. University Chemistry, 2021, 36(2): 1912058.

链接本文: https://www.dxhx.pku.edu.cn/CN/10.3866/PKU.DXHX201912058

https://www.dxhx.pku.edu.cn/CN/Y2021/V36/I2/1912058

图/表 10

图1

表1

表2

表3

表4

图2

图3

图4

表5

图5

参考文献 16

1	刘国杰; 黑恩成. 大学化学, 2013, 28 (5), 73.
2	罗渝然; 俞书勤; 张祖德; 姚天扬; 高盘良. 大学化学, 2010, 25 (3), 35.
3	邱晓航; 李一峻; 韩杰; 尚贞锋. 基础化学实验, 第2版北京: 科学出版社, 2017, 324
4	甘礼华; 陈龙武; 李天旻. 大学化学, 1996, 11 (3), 38.
5	陈九顺. 大学化学, 1995, 10 (3), 48.
6	王贵昌; 赵学庄. 大学化学, 1997, 12 (1), 24.
7	Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; et al. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford CT, 2013.
8	Gong J. B. ; Ma A. ; Rice S. A. J. Chem. Phys. 2005, 122, 204505. doi: 10.1063/1.1913398
9	Sun Z. P. ; Zheng Y. J. J. Phys. Chem. A 2015, 119, 2982. doi: 10.1021/jp511440w
10	Bowman J. M. ; Gazdy B. J. Phys. Chem. A 1997, 101, 6384.
11	胡建明; 李奕; 李俊; 章永凡; 周立新. 化学学报, 2003, 61 (4), 476. doi: 10.3321/j.issn:0567-7351.2003.04.004
12	Lasorne B. ; Gatti F. ; Baloitcha E. ; Meyer H. -D. ; Desouter-Lecomte M. J. Chem. Phys. 2004, 121, 644. doi: 10.1063/1.1760713
13	Koch D. M. ; Toubin C. ; Xu S. C. ; Peslherbe G. H. ; Hynes J. T. J. Phys. Chem. C 2007, 111, 15026. doi: 10.1021/jp076220h
14	Khalouf-Rivera J. ; Carvajal M. ; Santos L. F. ; Pérez-Bernal F. J. Phys. Chem. A 2019, 123, 9544. doi: 10.1021/acs.jpca.9b07338
15	Fukui K. Acc. Chem. Res. 1981, 14, 363. doi: 10.1021/ar00072a001
16	许秀芳. 大学化学, 2018, 33 (1), 49.

输入文件	各部分的名称	各部分的含义
#p B3LYP/6-311+G(d, p) opt freq	计算执行路径行	指定计算方法、基组、计算项目行首以#开始，段后加空行
HCN	标题行	作业的简要描述，段后加空行
0 1	电荷与多重度	确定电荷与电子自旋状态，段后无空行
C -2.73936174 1.56028366 0.00000000 H -3.80936174 1.56028366 0.00000000 N -1.59276174 1.56028366 0.00000000	分子说明部分 (此处用的是笛卡尔坐标格式，用内坐标也可以)	确定分子中各原子坐标，段后通常加空行

输入文件	各部分的名称	各部分的含义
#p B3LYP/6-311+G(d, p) opt freq	计算执行路径行	指定计算方法、基组、计算项目
HNC	标题行	作业的简要描述
0 1	电荷与多重度	确定电荷与电子自旋状态
C -3.37765962 1.59574466 0.00000000 N -2.23105962 1.59574466 0.00000000 H -1.23105962 1.59574466 0.00000000	分子说明部分(此处用的是笛卡尔坐标格式)	确定分子中各原子坐标

输入文件	各部分的名称	各部分的含义
#p B3LYP/6-311+G(d, p) opt=(ts, calcfc, noeigen) freq	计算执行路径行	指定计算方法、基组、计算项目
TS	标题行	作业的简要描述
0 1	电荷与多重度	确定电荷与电子自旋状态
C 0.00000000 0.11010811 -0.62544218 H 0.00000000 -0.99614153 -0.20995858 N 0.00000000 0.03041465 0.55797751	分子说明部分(此处用的是笛卡尔坐标格式)	确定分子中各原子坐标

输入文件	各部分的名称	各部分的含义
#p B3LYP/6-311+G(d, p) irc=(CalcFC, Recalc=10, maxpoints=200, stepsize=15, LQA)	计算执行路径行	指定计算方法、基组、计算项目
IRC	标题行	作业的简要描述
0 1	电荷与多重度	确定电荷与电子自旋状态
C 0.08091600 0.61884200 0.00000000 H -1.05191100 0.24372000 0.00000000 N 0.08091600 -0.56525400 0.00000000	分子说明部分(此处用的是笛卡尔坐标格式)	确定分子中各原子坐标

Species	E₀/(kJ·mol^-1)	H/(kJ·mol^-1)	G/(kJ·mol^-1)
HCN	-245321.81032	-245312.66570	-245372.64525
TS	-245136.11395	-245126.19744	-245191.24683
HNC	-245263.27479	-245253.40291	-245314.42741