Low temperature HD + $ortho$-/$para$-H$_2$ inelastic scattering of astrophysical interest

State-selected total cross sections and thermal rate coefficients are computed for the HD + $ortho$-/$para$-H$_2$ rotational energy transfer collision at low temperatures: 2 K $\lesssim$ T $\lesssim$ 300 K. A modified H$_2$-H$_2$ potential energy surface (PES) devised by Hinde is used for this pure quantum-mechanical dynamical computation. A comparison of the new results for the HD + $ortho$-/$para$-H$_2$ scattering problem and previous calculations computed with the use of other older PESs is presented and discussed.