We have calculated accurate quantum reactive and elastic cross sections for the prototypical barrierless reaction D + + H 2 (v = 0, j = 0) using a modified hyperspherical scattering method. The considered kinetic energy ranges from the ultracold to the Langevin regimes. A reaction rate coefficient practically constant in no less than eight orders of magnitude is obtained. The availability of accurate results for this system allows one to test the quantum theory by Jachymski et al. [K. Jachymski, M. Krych, P. S. Julienne, and Z. Idziaszek, Phys. Rev. Lett. 110, 213202 (2013)] in a nonuniversal case. The short-range reaction probability is rationalized using statistical model assumptions and related to a statistical factor. This provides a means to estimate one of the parameters that characterizes ultracold processes from first principles. The increasing availability of cold and ultracold samples of atoms and molecules has sprung great interest in chemical reactions at very low temperatures [1][2][3][4][5]. Although new experimental approaches [5] appear highly promising, advances in the field are hampered by technical problems in producing most molecules at low temperatures and high enough densities. In contrast to neutral species, ions can be easily trapped and cooled. The technology of Coulomb crystals in radiofrequency ion traps [6] and the possibility of combining them with traps for neutrals or with slow molecular beams [7,8] promise great progress in the analysis of ion-neutral reactions in the near future.Theoretical simulations employing standard ab initio approaches are not feasible for most of the systems thus far considered. For heavy systems (more convenient experimentally) there are no potential energy surfaces (PESs) accurate enough to describe processes near thresholds. Additionally, most of exact dynamical treatments face insurmountable problems in such regimes. However, in contrast to short-range (SR) chemical interactions, those occurring at long range (LR) can be more easily calculated. Moreover, theoretical approaches based only on the knowledge of the LR part of the PES have been able to describe recent experimental findings nearly quantitatively [1,9]. Indeed, processes at very low collision energies favor LR interactions, leading to the idea of universality in extreme cases [10]: the result of the collision depends exclusively on the LR behavior and not on the details of the PES. In this regard, recently proposed LR parametrization procedures [9,[11][12][13] Fitting experimental data, these models are able to predict nonmeasured values providing some insight into the underlying interactions. In particular, the approach by Jachymski et al., based on multichannel quantum-defect theory (MQDT), provides analytical expressions which can be easily compared with experimental data [9,14,15]. The model has been recently applied to a variety of systems [9,16,17]. In particular, for the Penning ionization of Ar by He( 3 S) [5], the rate coefficients have been fitted in a wide range of collision ener...