Quantum tunneling isotope exchange reaction H$^2$ + D$^-$ $\rightarrow$ HD + H$^-$

The tunneling reaction H$_2$ + D$^-$ $\rightarrow$ HD + H$^-$ was studied in a recent experimental work at low temperatures (10, 19, and 23~K) by Endres {\it et al.}[1]. An upper limit of the rate coefficient was found to be about 10$^{-18}$ cm$^3$/s. In the present study, reaction probabilities are determined using the ABC program developed by Skouteris {\it et al.}[2]. The probabilities for ortho-H$_2$ and para-H$_2$ in their ground rovibrational states are obtained numerically at collision energies above 50~meV with the total angular momentum $J$ = 0 - 15 and extrapolated below 50~meV using a WKB approach. Thermally averaged rate coefficients for ortho- and para-H$_2$ are obtained; the largest one, for ortho-H$_2$ is about $3.1\times10^{-20}$ cm$^3$/s, which agrees with the experimental results. [1] Endres {\it et al.} PRA 95, 022706 (2017) [2] Skouteris {\it et al.} Comput. Phys. Commun. 133, 128 (2000)