  We have calculated quantum reactive and elastic cross-sections for D$^{+}+$
para-H$_2$($v$=0, $j$=0) $\rightarrow$ H$^+$ + HD collisons using the
hyperspherical quantum reactive scattering method [Chem. Phys. Lett. 1990,169,
473]. The H$_{3}^{+}$ system is the prototype of barrierless ion-molecule
reactions, apart from its relevance in astrochemistry. The considered collision
energy ranges from the ultracold regime, where only one partial wave is open,
up to the Langevin regime, where many of them contribute. At very low kinetic
energies, both an accurate description of the long-range (LR) region in the
potential energy surface (PES), and long dynamical propagations, up to
distances of 10$^{5}$ a0, are required. Accordingly, calculations have been
carried out on the PES by Velilla et al. [J. Chem. Phys. 2008, 129,084307]
which accurately reproduces the LR interactions. Hyperspherical methodology was
recently modified in order to allow the accurate inclusion of LR interactions
while minimizing the computational expense. Such implementation is shown
particularly suitable for systems involving ions, where the $R^{-4}$ behavior
largely extends the range of the potential. We find a reaction rate coefficient
which remains almost constant in a kinetic energy range of more than ten orders
of magnitude. In particular, the value reached in the Wigner regime, where only
one partial wave is open, is paradoxically the Langevin classical value (within
a few percent) expected at high energies. Results are discussed in terms of
{\em universality} and related to the recently published quantum defect theory
[Phys. Rev. Lett. 2013, 110, 213202]. Since the system has small exothermicity
and a low number of channels, it is possible to test such model in a case where
the loss probability at short range is appreciably far from unity.
