发展新的波包传播、坐标变换方法研究态态分子反应动力学

State-to-state quantum reaction dynamics study provides the most detailed observable of reaction dynamics experiments. But the theoretical calculations are still a great challenge of amount of the computation.We will develop two methods to raise the efficiencies of the wave-packet propagation and transformation between the reactant and product coordinates, respectively. A new split-operator method is introduced in wave-packet propagation and decreases communication among the computational nodes. Therefore, it solves the bottleneck of the parallelism. A new method separates the coordinate transformation into two parts, and turns the five-loop calculation to two four-loop calculations. It decreases the amount of computation and needs no extra memory. A state-to-state dynamics code will developed with the new methods and parallelly implemented on GPUs/CPUs. The methods will be tested and compared with the previous works in the amount of computation and applicability in different systems. We will discuss the feasibility to extend the new methods to tetratomic systems. At last the code will be applied on several prototype reactions to obtain accurately quantitative results, such as the state-to-state reaction probabilities, integral and differential cross section, to study the mechanism of reaction.

态态量子动力学研究给出了最精细的反应动力学实验可观测量。但是，理论计算的计算量仍是一个很大的挑战。我们将发展两种新的方法分别提高态态动力学研究中波包传播效率和反应物产物坐标系之间的变换效率。新的波包传播方法将引入一个分裂算符，显著减少计算节点间的通信，解决并行瓶颈。新的坐标变换方法将变换分成两部分，把一个五重循环变成两个四重循环，减少变换所需的计算量的同时不消耗系统更多的内存。利用新方法，编写一整套态态动力学程序，并实现GPUs/CPUs并行化。与前人的方法进行测试比较，比较各种方案之间的计算量区别以及应用于各种反应体系时不同的适用性。讨论新方法推广到四原子体系的可行性。最后应用程序计算几个典型体系的态态反应几率、积分截面、微分截面等精确定量的动力学结果，研究反应机理。

态态量子动力学计算给出了反应几率、微分散射截面等信息，为深入了解小分子化学反应过程提供了途径。我们发展了新的波包传播方法和坐标变换方法，并编写了一整套态态动力学程序实现态态量子动力学研究。新的波包传播方法将引入一个分裂算符，显著减少计算节点间的通信，解决并行瓶颈。新的坐标变换方法将变换分成两部分，把一个五重循环变成两个四重循环，减少变换所需的计算量的同时不消耗系统更多的内存。程序实现GPUs/CPUs并行化。测试表明在两个K20m GPU上的平均运算速度是CPU串行的运算速度的83.5倍，最大的加速比是122。在此基础上，我们应用程序实现了Ne + D2+，H(2S) + NH (X3Σ-)，S + HD 和 Li + H2+等体系的研究，得到了态态反应几率、积分截面、微分截面等精确定量的动力学结果，研究了反应机理。

内容获取失败，请点击重试