Gamallo Pablo, Akpinar Sinan, Defazio Paolo, Petrongolo Carlo
Departament de Ciència de Materials i Química Física & Institut de Química Teòrica i Computacional, Universitat de Barcelona, C/Martí i Franquès 1, 08028 Barcelona, Spain.
Department of Physics, Firat University, 23169 Elazig, Turkey.
Phys Chem Chem Phys. 2017 Feb 8;19(6):4454-4461. doi: 10.1039/c6cp07756k.
Following previous studies on the O(P) + H(XΣ) collisions, we present the nonadiabatic quantum dynamics of the reactions OH(XΠ) + H' → OH'(XΠ) + H, exchange (e), → OH(XΣ) + H'(S), quenching (q), and → OH' (XΣ) + H(S), exchange-quenching (eq). The reactants and products correlate via the ground X[combining tilde]A'' and first excited ÃA' electronic states of OH, which are the degenerate components of linear Π species. Therefore, they are strongly perturbed by nonadiabatic Renner-Teller (RT) effects, opening the (q) and (eq) channels that are closed in the Born-Oppenheimer approximation. Using accurate potential energy surfaces (PESs) and RT matrix elements, initial-state-resolved reaction probabilities, real-time dynamics, cross sections, and rate constants of the product channels are obtained through the time-dependent real wavepacket (WP) method and full coupled-channel calculations. Owing to the nonadiabatic couplings, the WP jumps from the excited ÃA' surface to the X[combining tilde]A'' ground PES, avoiding any barrier, opening the quenching channels, and giving many collision complexes into the deep minima of both PESs, as it is clearly shown by the oscillations of the reaction probabilities and by the time-dependent WP dynamics. All the results show that the nonadiabatic-RT channels (q) and (eq) are highly reactive, much more than the adiabatic one (e), pointing out large RT effects. The reactivity of the quenching channels is similar, accounting for 97% of the overall reactivity. In fact, the maximum values of the (q) and (eq) cross sections σ and σ are equal to 31.6 Å, whereas the maximum σ value equals 1.34 Å, and the maximum values of the rate constants k, k, and k are 2.07 × 10, 2.45 × 10, and 0.23 × 10 cm s. Some calculations show that the centrifugal-sudden and the truncated coupled-channel approximations cannot be employed for the (q) channel. After a sharp increase at the threshold, σ and σ decrease at larger collision energies while σ and the rate constants depend slightly on the collision energy and temperature, respectively. These findings are consistent with the barrierless nature of both PESs and the exoergicity of the quenching products, with the small role played by the centrifugal and RT barriers in the reactant channel, and with the large RT couplings in the OH intermediates. Finally, we contrast the present results with those for the opposite reactions O + H and for the nonadiabatic-induced quenchings NH + H' and OH + H'.
继之前关于O(P) + H(XΣ)碰撞的研究之后,我们展示了反应OH(XΠ) + H' → OH'(XΠ) + H(交换反应 (e))、→ OH(XΣ) + H'(S)(猝灭反应 (q))以及→ OH'(XΣ) + H(S)(交换 - 猝灭反应 (eq))的非绝热量子动力学。反应物和产物通过OH的基态X[波浪线]A''和第一激发态ÃA'电子态相关联,它们是线性Π物种的简并组分。因此,它们受到非绝热伦纳 - 泰勒(RT)效应的强烈扰动,开启了在玻恩 - 奥本海默近似中关闭的(q)和(eq)通道。使用精确的势能面(PESs)和RT矩阵元,通过含时实波包(WP)方法和全耦合通道计算获得了产物通道的初态分辨反应概率、实时动力学、截面和速率常数。由于非绝热耦合,波包从激发的ÃA'表面跃迁到X[波浪线]A''基态势能面,避开任何势垒,开启猝灭通道,并在两个势能面的深阱中产生许多碰撞复合物,反应概率的振荡和含时波包动力学清楚地表明了这一点。所有结果表明,非绝热 - RT通道(q)和(eq)具有高反应活性,远高于绝热通道(e),这表明存在较大的RT效应。猝灭通道的反应活性相似,占总反应活性的97%。事实上,(q)和(eq)截面σ和σ的最大值等于31.6 Å,而最大σ值等于1.34 Å,速率常数k、k和k的最大值分别为2.07×10、2.45×10和0.23×10 cm³/s。一些计算表明,对于(q)通道不能采用离心突然近似和截断耦合通道近似。在阈值处急剧增加后,σ和σ在较大碰撞能量下减小,而σ和速率常数分别对碰撞能量和温度的依赖性较小。这些发现与两个势能面的无势垒性质以及猝灭产物的放能性一致,与反应物通道中离心和RT势垒所起的小作用一致,也与OH中间体中的大RT耦合一致。最后,我们将当前结果与相反反应O + H以及非绝热诱导猝灭NH + H'和OH + H'的结果进行了对比。