The reactions between ground-state H + 2 (X 2 + g (v = 0, J = 0)) and D 2 forming HD + 2 + H and H 2 D + + D were investigated in the range of collision energies E coll between E coll /k B = 0 and 10 K using a merged-beam approach. The reaction rates measured experimentally are compared to those obtained for the reaction between H + 2 and H 2 forming H + 3 + H under similar experimental conditions. Below 1 K, a clear enhancement of the reaction rate coefficient compared to the Langevin rate measured at higher collision energies was observed in both reaction systems. This enhancement is interpreted as originating from the interaction between the charge of H + 2 and the quadrupole of para D 2 and ortho H 2 molecules in the J = 1 rotational level. The enhancement of the reaction with D 2 was found to be significantly less than that of the reaction with H 2 , reflecting the relative population of the J = 1 rotational level of H 2 (75%) and D 2 (33%) in natural samples at low temperatures. Simulations of the experimental results based on the theoretical predictions of the reaction cross sections by Dashevskaya et al. [J. Chem. Phys. 145, 244315 (2016)] reveal agreement within the experimental uncertainties. The branching ratio η of the reaction involving H + 2 and D 2 and forming H 2 D + and D (η) near E coll = 0 was determined to be 0.341(15). Time-of-flight measurements of the velocity distributions of the reaction products are compatible with an isotropic product emission with an average total kinetic energy of 0.45(5) eV for both channels, representing about 30% of the total energy released by the reaction.