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