Thomas Phillip S, Somers Mark F, Hoekstra Anne W, Kroes Geert-Jan
Leiden Institute of Chemistry - Theoretical Chemistry, Postbus 9502, Leiden 2300 RA, Netherlands.
Phys Chem Chem Phys. 2012 Jun 28;14(24):8628-43. doi: 10.1039/c2cp40173h. Epub 2012 May 18.
Gas phase and surface reactions involving polyatomic molecules are of central importance to chemical physics, and require accurately fit potential energy surfaces describing the interaction in their systems. Here, we propose a method for generating a High Dimensional Model Representation (HDMR) of a multidimensional potential energy surface (PES) and apply it to reactive molecule-surface scattering problems. In the HDMR treatment, all N degrees of freedom (DOF) of an N-dimensional PES are represented but only n < N are explicitly coupled. The HDMR is obtained from Chebyshev polynomial expansions for each degree of freedom, where expansion coefficients are efficiently evaluated using discrete cosine transform (DCT) algorithms and properties of Chebyshev polynomials. HDMR surfaces constructed for the reactive scattering of H2 from Pt(111) and Cu(111) are used in quantum dynamics simulations; the resultant state-resolved reaction and scattering probabilities are compared to those from simulations using full (6D) PESs and n-mode PESs from previous work. The results are encouraging, and suggest that this method may be applicable to "late barrier" reactive systems for which the previously-used n-mode representation fails.
涉及多原子分子的气相和表面反应对于化学物理至关重要,并且需要精确拟合描述其系统中相互作用的势能面。在此,我们提出一种生成多维势能面(PES)的高维模型表示(HDMR)的方法,并将其应用于反应性分子-表面散射问题。在HDMR处理中,N维PES的所有N个自由度(DOF)都有表示,但只有n < N个自由度被明确耦合。HDMR是通过对每个自由度进行切比雪夫多项式展开得到的,其中展开系数使用离散余弦变换(DCT)算法和切比雪夫多项式的性质进行有效评估。为H2从Pt(111)和Cu(111)的反应性散射构建的HDMR表面用于量子动力学模拟;将得到的态分辨反应和散射概率与之前工作中使用全(6D)PES和n模PES的模拟结果进行比较。结果令人鼓舞,并表明该方法可能适用于先前使用的n模表示失效的“晚势垒”反应系统。