6.23(-11) & 7.10(-11)\
\
25 & 26 $^\dag$ & 1.14(-10) & 1.50(-10) & 1.83(-10) & 2.30(-10)\
25 & 27 $^\dag$ & 1.13(-10) & 1.05(-10) & 1.18(-10) & 1.31(-10)\
25 & 28 $^\dag$ & 8.55(-11) & 8.38(-11) & 8.34(-11) & 7.76(-11)\
25 & 29 $^\dag$ & 7.67(-11) & 7.45(-11) & 8.31(-11) & 7.02(-11)\
\
35 & 36 $^\dag$ & 1.16(-10) & 1.34(-10) & 1.73(-10) & 2.32(-10)\
35 & 37 $^\dag$ & 9.63(-11) & 1.12(-10) & 1.21(-10) & 1.11(-10)\
35 & 38 $^\dag$ & 8.33(-11) & 9.20(-11) & 8.56(-11) & 9.19(-11)\
35 & 39 $^\dag$ & 8.77(-11) & 8.51(-11) & 7.65(-11) & 7.51(-11)\
\

For H$_2$, all deexcitation rates $k_{J_1J_1'}(T)$, $J_1\neq J_1'\leq
50$, are fitted with the following formula [@wernli06]: $$\label{eq:fit}
    \log_{10}\left(k_{J_1J_1'}(T)\right)=\sum_{n=0}^{
    4}a^{(n)}_{J_1J_1'} x^n$$ where $x=T^{-1/6}$. As some transitions have zero probability within the QCT approach, the above formula was employed when rates were bigger than 10$^{-12}$ cm$^3$s$^{-1}$ for at least one grid temperature. For these rates, null grid values were replaced by a very small value, namely 10$^{-14}$ cm$^3$s$^{-1}$, to avoid fitting irregularities. All rates not fulfilling this condition are set to zero. Note that below 20 K, QCT rates for low-probability transitions may show a non physical behaviour. All $a^{(n)}_{J_1J_1'}$ 