The correct answer was given by Cosmas Zachos, however, I want to add a remark concerning the assertion that "$\hat p$ converges to $0$ in the classical limit"
One can actually look directly at the limit of the impulsion operator in a more mathematically precise way: the momentum operator is defined by $\hat p = -i\hbar\nabla$. Notice first that it is an unbounded operator on $L^2$, so one should be careful about what "convergence to $0$" means (in operator norm, $\hbar\cdot\infty = ∞$ does not converge to $0$).A way to look at the convergence of unbounded operator is by using weak convergence. An example is to look at the limit of its mean value
$\lim_{\hbar\to 0}\mathrm{Tr}(\hat p\,\rho)$ for any nice density matrix $\rho$.
However, when doing the classical limit of quantum mechanics, one should forget that the size of operators also depend on $\hbar$.
A classical way to perform the limit $\hbar\to 0$ is to introduce the Wigner transform
$$
W_{\rho}(x,p) = \frac{1}{\hbar^d}\int_{\mathbb R^3} e^{-i\,y\cdot p/\hbar}\,\rho(x+y/2,x-y/2)\,\mathrm d y
$$
which is indeed an object that converges to a classical distribution of the phase space. A quick computation then shows that
$$
\mathrm{Tr}(\hat p\,\rho) = \int_{\mathbb R^6} p\, W_{\rho}(x,p)\,\mathrm d x\,\mathrm d p.
$$
In particular, if $W_{\rho}(x,p)$ converges to a classical function of the phase space $f$ in a suitable topology (which is the assumption that indicates that $\rho$ is $\hbar$ dependent), then
$$
\mathrm{Tr}(\hat p\,\rho) \to \int_{\mathbb R^6} p\, f(x,p)\,\mathrm d x\,\mathrm d p
$$
which is not $0$.