the small-$\rho$ asymptotics is governed by the large-$k$ behavior of the integrand, which gives a scaling as $\propto \rho^{(2-1.9)}\log\rho$, and numerics indicates that this is indeed the case: