基元化学反应态态动力学研究

化学反应动力学是化学领域最基础的学科之一,量子态分辨的基元化学反应动力学在最为基本的原子与分子的层次上对化学反应的机制提供深刻的理解。该领域的科学家们通过精心设计的实验和高精度的理论计算,使得态态反应动力学在过去的半个多世纪中取得了长足的进步,实验和理论的相互结合极大地促进了我们对化学反应本质的认识。本文从实验研究的角度,通过对实验技术的发展和对H₂O光解离、H+H₂、F+H₂、Cl+H₂、OH+H₂、F+CH₄等具体实例的态态动力学研究的简介,概况介绍了过去二十年里态态化学反应动力学研究所取得的进展,希望借此为读者提供对化学反应动力学领域的一个概略认识。     

H+Cl_2→HCl+Cl体系反应几率的量子动力学计算

利用含时波包动力学方法在ＧＨＮＳ势能面上计算了总角动量Ｊ＝０时Ｈ＋Ｃｌ２体系初态确定的累积反应几率,讨论了平动能、Ｃｌ２振动和转动激发对反应几率的影响;并且讨论了Ｈ＋Ｃｌ２体系的波包运动特征．在一定的平动能范围内,Ｃｌ２的振动激发不利于反应的进行;而Ｃｌ２的转动激发特别是第一转动激发能够加快反应．这些结论和反应是早势垒反应的事实是一致的．由于Ｃｌ２是由两个较重的原子组成的,振动和转动态能级密集;同时由于是早势垒反应,能垒较低,波包在运动过程中在两侧的势垒壁之间发生多次反射致使波包发生很大范围内的干涉,需要较多的网格点才能正确地描述体系的波函数．这两个因素使得此体系的计算量较一般三原子体系要大得多．     

三维H+H2(v,j)→H2(v′,j′)+H反应中复合态生成及产物转动态分布的量子散射理论研究

用排列通道线性组合-散射波函数(LCAC-SW,linear combination of arrangement channelsscattering wavefunction)量子反应散射方法计算了H+H2(v,j)→H2(v′,j′)+H三维态-态反应几率,分析了反应体系的复合态生成(或能量共振结构),并由产物的转动态分布解释了能量共振的起源来自于平动态-内态之间的干涉效应.     

小分子光解动力学的理论研究

光解过程是化学中的核心问题之一。量子态分辨的光解动力学可以使人们在原子与分子的层次上深刻理解光解反应的机制。态-态水平的光解动力学在过去四十年中取得了长足的进步,实验和理论的相互结合极大地促进了我们对光解反应本质的认识。本文综述了小分子态-态光解动力学的理论研究进展,总结了H2O和CH3I这两个最具代表性体系的态-态光解动力学研究成果,并提出了该领域未来面临的挑战。     

三维H+H~2(v,j)→H~2(v', j')+H反应中复合态生成及产物转动 态分布的量子 散射理论研究

用排列通道线性组合-散射波函数(LCAC-SW,linearcombinationofarrangementchannels-scatteringwavefunction)量子反应散射方法计算了H+H~2(v,j)→H~2(v',j')+H三维态-态反应几率,分析了反应体系的复合态生成(或能量共振结构),并由产物的转动态分布解释了能量共振的起源来自于平动态-内态之间的干涉效应。     

H+Cl2→ HCl+Cl体系反应几率的量子动力学计算

利用含时波包动力学方法在GHNS势能面上计算了总角动量J=0时H+Cl₂体系初态确定的累积反应几率,讨论了平动能、Cl₂振动和转动激发对反应几率的影响;并且讨论了H＋Cl₂体系的波包运动特征.在一定的平动能范围内,Cl₂的振动激发不利于反应的进行; 而Cl₂的转动激发特别是第一转动激发能够加快反应.这些结论和反应是早势垒反应的事实是一致的.由于Cl₂是由两个较重的原子组成的,振动和转动态能级密集;同时由于是早势垒反应,能垒较低,波包在运动过程中在两侧的势垒壁之间发生多次反射致使波包发生很大范围内的干涉,需要较多的网格点才能正确地描述体系的波函数.这两个因素使得此体系的计算量较一般三原子体系要大得多.     

态-态热力学函数及其应用 I: 双分子交换反应

本文把Levine的态-态热力学函数由振动-振动态推广到振转-振转态, 对A+BC(v, j)→AB(v', j')+C双分子交换反应给出了Gibbs自由能△G^o(v, j→v',j', T)和化学平均常数K(v, j→v', j', T)在谐振子-截锥转子模型下的详细计算公式, 对H+O2(v, j)→HO(v', j')+O作了数值计算, 结果与由实验归纳出的定性规律相符合。利用本文给出的公式不仅可以对化学反应过程描述得更加细致和深刻, 而且可以方便地讨论分子的振动-转动态耦合对化学反应性的影响。     

态-态热力学函数及其应用 I: 双分子交换反应

本文把Levine的态-态热力学函数由振动-振动态推广到振转-振转态, 对A+BC(v, j)→AB(v', j')+C双分子交换反应给出了Gibbs自由能△G^o(v, j→v',j', T)和化学平均常数K(v, j→v', j', T)在谐振子-截锥转子模型下的详细计算公式, 对H+O2(v, j)→HO(v', j')+O作了数值计算, 结果与由实验归纳出的定性规律相符合。利用本文给出的公式不仅可以对化学反应过程描述得更加细致和深刻, 而且可以方便地讨论分子的振动-转动态耦合对化学反应性的影响。     

三维H+H~2(v, j)→H~2(v', j')+H反应中复合态生成及产物转动 态分布的量子 散射理论研究

用排列通道线性组合-散射波函数(LCAC-SW,linearcombinationofarrangementchannels-scatteringwavefunction)量子反应散射方法计算了H+H~2(v,j)→H~2(v',j')+H三维态-态反应几率,分析了反应体系的复合态生成(或能量共振结构),并由产物的转动态分布解释了能量共振的起源来自于平动态-内态之间的干涉效应。     

D_2在Ni(100)表面离解吸附的量子动力学研究

用量子含时波包法研究了Ｄ２在镍表面顶－桥位上离解吸附量子动力学．计算了不同入射动能及初始振转态的离解几率．讨论了分子的同核对称性、转动取向和振动激发对离解几率的影响，并与其他理论计算结果做了比较．     

Quantum wave packet study of the insertion reaction S+H2

S(¹D)+H₂→SH+H reaction at zero total angular momentum has been studied by using a time-dependent quantum real wave packet method. State-to-state and state-to-all reactive scattering probabilities for a broad range of energy are calculated. The probabilities show many sharp peaks that ascribed to reactive scattering resonances. The density plots of the wave function shows that the reaction presents a pure insertion pathway.     

Higher-order split operator schemes for solving the Schrödinger equation in the time-dependent wave packet method: applications to triatomic reactive scattering calculations

The efficiency of the numerical propagators for solving the time-dependent Schr?dinger equation in the wave packet approach to reactive scattering is of vital importance. In this Perspective, we first briefly review the propagators used in quantum reactive scattering calculations and their applications to triatomic reactions. Then we present a detailed comparison of about thirty higher-order split operator propagators for solving the Schr?dinger equation with their applications to the wave packet evolution within a one-dimensional Morse potential, and the total reaction probability calculations for the H + HD, H + NH, H + O(2), and F + HD reactions. These four triatomic reactions have quite different dynamic characteristics and thus provide a comprehensive picture of the relative advantages of these higher-order propagation methods for describing reactive scattering dynamics. Our calculations reveal that the most often used second-order split operator method is typically more efficient for a direct reaction, particularly for those involving flat potential energy surfaces. However, the optimal higher-order split operator methods are more suitable for a reaction with resonances and intermediate complexes or a reaction experiencing potential energy surface with fluctuations of considerable amplitude. Three 4th-order and one 6th-order split operator methods, which are most efficient for solving reactive scattering in various conditions among the tested ones, are recommended for general applications. In addition, a brief discussion on the relative performance between the Chebyshev real wave packet method and the split operator method is given. The results in this Perspective are expected to stimulate more applications of (high-order) split operators to the quantum reactive scattering calculation and other related problems.     

Quantum wave packet study of N(2D) + H2 reactive scattering

N(²D) + H₂ → NH + H reaction at zero total angular momentum is studied by using a time dependent quantum wave packet method. State‐to‐state and state‐to‐all reactive scattering probabilities for a broad range of energy are calculated. The probabilities show many sharp peaks that ascribed to reactive scattering resonances. The probabilities for J > 0 are estimated by using the J‐shifting method. The integral cross sections and thermal rate constants are then calculated. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Moving boundary truncated grid method: Multidimensional quantum dynamics

The moving boundary truncated grid (TG) method is used to study wave packet dynamics of multidimensional quantum systems. As time evolves, appropriate Eulerian grid points required for propagating a wave packet are activated and deactivated with no advance information about the dynamics. This method is applied to the Henon-Heiles potential and wave packet barrier scattering in two, three, and four dimensions. Computational results demonstrate that the TG method not only leads to a great reduction in the number of grid points needed to perform accurate calculations but also is computationally more efficient than the full grid calculations.     

Plane wave packet formulation of atom-plus-diatom quantum reactive scattering

We recently interpreted several reactive scattering experiments using a plane wave packet (PWP) formulation of quantum scattering theory [see, e.g., S. C. Althorpe, F. Fernandez-Alonso, B. D. Bean, J. D. Ayers, A. E. Pomerantz, R. N. Zare, and E. Wrede, Nature (London) 416, 67 (2002)]. This paper presents the first derivation of this formulation for atom-plus-diatom reactive scattering, and explains its relation to conventional time-independent reactive scattering. We generalize recent results for spherical-particle scattering [S. C. Althorpe, Phys. Rev. A 69, 042702 (2004)] to atom-rigid-rotor scattering in the space-fixed frame, atom-rigid-rotor scattering in the body-fixed frame, and finally A+BC rearrangement scattering. The reactive scattering is initiated by a plane wave packet, describing the A+BC reagents in center-of-mass scattering coordinates, and is detected by projecting onto a series of AC+B (or AB+C) plane wave "probe" packets. The plane wave packets are localized at the closest distance from the scattering center at which the interaction potential can be neglected. The time evolution of the initial plane wave packet provides a clear visualization of the scattering into space of the reaction products. The projection onto the probe packets yields the time-independent, state-to-state scattering amplitude, and hence the differential cross section. We explain how best to implement the PWP approach in a numerical computation, and illustrate this with a detailed application to the H+D2 reaction.     

Quantum wave packet calculation of reaction probabilities,cross sections,and rate constants for H+ + LiH → Li + H reaction

The H⁺ + LiH → Li + H reactive scattering has been studied using a quantum real wave packet method. The state‐to‐state and state‐to‐all reaction probabilities for the entitled collision have been calculated at zero total angular momentum. The probabilities for J > 0 are estimated from the J = 0 results by using J‐shifting approximation based on the Capture model. The integral cross sections and thermal rate constants are then calculated. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Analysis of barrier scattering with real and complex quantum trajectories

In this study, an analysis of the one-dimensional Eckart and Gaussian barrier scattering problems is undertaken using approximate quantum trajectories. Individual quantum trajectories are computed using the derivative propagation method (DPM). Both real-valued and complex-valued DPM quantum trajectories are employed. Of interest are the deep tunneling and the higher energy barrier scattering problems in cases in which the scattering barrier is "thick" by comparison to the width of the initial wave packet. For higher energy scattering problems, it is found that real-valued DPM trajectories very accurately reproduce the transmitted probability densities at low orders when compared to large fixed-grid calculations. However, higher orders must be introduced to obtain good probabilities for deep tunneling problems. Complex-valued DPM is found to accurately reproduce transmitted probability densities at low order for both the deep tunneling and the higher energy scattering problems. Of particular note, complex-classical trajectories are found to very nearly give the exact result for the deep barrier tunneling scattering problem, and the complex DPM converges well at high orders for these thick barrier scattering problems. A variety of analyses are performed to elucidate the dynamics of complex-valued DPM trajectories. The complex-extended barrier potentials are examined in detail, including an analysis of the complex force. Of particular interest are initial conditions for complex-valued DPM trajectories known as isochrones. All trajectories launched from an isochrone arrive on the real axis on the transmitted side of the barrier at the same time. The computation and properties of isochrones as well as the behavior of the initial wave packet in the complex plane are also examined.     

Photo-induced desorption of NO from NiO(100): calculation of the four-dimensional potential energy surfaces and systematic wave packet studies

The velocity distributions of the laser-induced desorption of NO molecules from an epitaxially grown film of NiO(100) on Ni(100) have been studied [Mull et al., J. Chem. Phys., 1992, 96, 7108]. A pronounced bimodality of velocity distributions has been found, where the NO molecules desorbing with higher velocities exhibit a coupling to the rotational quantum states J. In this article we present simulations of state resolved velocity distributions on a full ab initio level. As a basis for this quantum mechanical treatment a 4D potential energy surface (PES) was constructed for the electronic ground and a representative excited state, using a NiO5Mg(18+)13 cluster. The PESs of the electronic ground and an excited state were calculated at the CASPT2 and the configuration interaction (CI) level of theory, respectively. Multi-dimensional quantum wave packet simulations on these two surfaces were performed for different sets of degrees of freedom. Our key finding is that at least a 3D wave packet simulation, in which the desorption coordinate Z, polar angle theta and lateral coordinate X are included, is necessary to allow the simulation of experimental velocity distributions. Analysis of the wave packet dynamics demonstrates that essentially the lateral coordinate, which was neglected in previous studies [Klüner et al., Phys. Rev. Lett. 1998, 80, 5208], is responsible for the experimentally observed bimodality. An extensive analysis shows that the bimodality is due to a bifurcation of the wave packet on the excited state PES, where the motion of the molecule parallel to the surface plays a decisive role.