【作者】 姚丽；

【导师】 宋鹤山； 韩克利；

【作者基本信息】 大连理工大学 ， 理论物理， 2004， 博士

【摘要】 分子反应动力学是在分子水平上研究化学反应微观动态和机理的一门学科。它不仅是宏观化学运动学的基础，同时也为解释基元化学反应中的基本现象提供了丰富的知识。长期以来人们致力于化学反应的精确量子动力学研究，并发展了许多研究方法。目前含时波包量子动力学方法已经成为人们研究三原子和四原子体系强有力的计算工具。采用的量子力学方法能够完全从量子力学基本原理出发，真实而完全的反映原子核在势能面上的运动。相对非含时方法，含时波包方法有许多优点，除了数值计算上效率高之外，它还为动力学提供了物理意义明确而直观的图像，具有经典的直观性又不乏量子力学的准确性，因而，受到人们的喜爱。含时波包方法已经用于研究了许多三原子和四原子反应体系了，但是对反应的含时波包精确计算至今仍是一项非常复杂的工作。 在气相化学反应动力学发展过程中，H+HCl及其同位素衍生物反应起着非常重要的作用；在大气化学和燃烧化学中，F+CH₄及其同位素衍生物的反应是氢提取和高放热反应的一个典型反应原型。所以本文对它们进行了研究。 本文采用含时量子波包方法（理论）对三原子反应H+HCl，H+DCl，D+HCl，D+DCl进行了精确的全维计算，在由Bian和Wemer提出的BW2势能面上研究了它们的动力学行为。当为线性过渡态时，该势能面上提取反应通道的经典势垒为0.184eV，对于交换反应通道也存在一个线性过渡态计算的势垒为0.776eV。计算了这四个反应在给定初始态下的反应几率，积分反应截面及初始态为激发态对反应过程及结果的影响等。在分析结果的时候我们将零点能和隧道效应同时考虑在内。再利用J-shifting方法计算了反应的速率常数。 为了实现量子动力学含时波包方法对多原子反应体系的研究，本文采用张增辉教授提出的一个半刚性振转子（SVRT）的模型，对F+CH₄反应及其同位素衍生物反应进行了含时波包动力学的研究。其势能面用了矫正后的儿势能面已由Corchado，J．C．报道。该势能面的反应势垒0.066eV。实验和理论研究表明在反应中C-H键不影响反应碰撞，它只是一个旁观者。在本文中C-H键作为常数处理使得计算大大简化了。我们计算得到了不同初始振转态的总反应几率，积分反应截面和速率常数。从得到的数值结果可以看出，总反应几率随着碰撞能有一个很尖锐的峰，这一般来讲是与动力学振荡有关系。 并将所有反应的结果和已有的理论以及实验结果进行了比较，发现它们符合得很好，我们又详细的分析了一些细微的差别及这些反应的立体动力学行为。更多还原

【Abstract】 Molecule reaction dynamics is a science of studying microcopic feature and mechanism of chemical reaction in molecular and atomic level. It can provide basic knowledge for interpreting general phenomena in the chemical reactions, and become the base of the macroscopic reaction dynamics.Quantum approach could describe the motion of atomic nuclear on the potential energy surface vividly and completely based on the first principle. Much effort has been devoted to the development of accurate quantum dynamics methods to perform reaction scattering calculation for elementary chemical reactions. At present, time-dependent (TD) quantum wave packet approach has emerged as a powerful computational tool for studying quantum reaction dynamics of triatomic and tetraatomic systems. The main attraction of solving the time-dependent, as opposed to the time-independent, Schrodinger equation stems from its favorable computational scaling with the number of basis functions. This makes the TD wavepacket approach an attractive choice for studying reaction dynamics of large molecular systems. In addition, the TD wavepacket approach is conceptually simple and provides a classical-like interpretation of the numerical results through time propagation of the wavepacket.The gas-phase reaction of H+HC1 and the corresponding to isotopic reactions system have played a major role in the development of chemical kinetics and to the environment in atmospheric chemistry. The reaction of fluorine atom with methane F+CH4HF+CH3 reaction and the corresponding isotopic reactions have been played an important role in atmospheric and combustion chemistry which deserves us to have more study. The reaction systems represent an especially important prototype for highly exothermic hydrogen abstraction reactions of F with hydrocarbons.A theoretical study of the dynamics is given in the thesis of the H+HC1, H+DC1, D+HC1, D+DC1 systems on the potential energy surface (PES) pubilshed by Bian and Werner. In the present work, the four reaction probabilities have been calculated by employing TDWP method for several values of the total angular momentum quantum number J> 0. Reaction probabilities are calculated from various initial rotational states of the reagent. The effects of the zero-point energy and the tunneling effect of the reactions are considered for depicting the behavior of the reactions. Those have then been used to estimate cross sections for exchanged and abstraction channels. The abstraction reaction for the H+HC1 reaction has a collinear transition state, and for the H+HC1 exchange reaction, which also has a collinear transition state. The thermal rate constants of HC1 are calculated by employing the uniform J -shifting method, and a comparison with experimental measurement is provided.In this work, the reaction of fluorine with methane and the isotopic reactions system are calculated through the time-dependent quantum scattering, which has been carried out using the SVRT model, based on the analytic the modified Jl (MJ1) PES which reported by Corchado, J. C. etc. The barrier height about 1.8Kcal/mol, which among the experimental estimate of 0.96-2.39 kcal/mol. The reagent molecule CH4 that is consisted of an H atom and a CH3 fragment which is treated as a semirigid vibrating rotor. The fragment CH3 is fixed in the geometry in the basic SVRT model. For example, since during the reactant to thetransition state the C-H bond in the CH3 group essentially remains constant, which can be treated as a spectator bond, therefore we fixed it at its transition state value-1.090A.In the current study, reaction probability, cross section, and rate constant are calculated for the title reaction from the ground state of the reagent on the MJl PES. Numerical calculation shows that the oscillatory structures in the energy dependence of the calculated reaction probability. Those structures are generally associated with broad dynamical resonances.The calculated results can compared with the results of previous calculations in the references and reaction dynam更多还原