Technically, if we want to prove an operator $\hat{A}$ to be an hermitian,we should prove $ \left< \psi \hat{A}|\psi \right> =\left< \psi |\hat{A}\psi \right>$. It works well in Cartesian coordinate when proving the momentum operator $\hat{p}$.
However when dealing with $\hat{p}_r\equiv-i\hbar\left(\frac{\partial}{\partial r}+\frac{1}{r}\right)$ in Spherical coordinate.
In Spherical coordinate system, we have:
$$ \langle\psi \hat{A} \mid \psi\rangle=\int_{0}^{\infty} \int_{0}^{\pi} \int_{-\frac{\pi}{2}}^{\frac{\pi}{2}}(\hat{A} \psi)^{\dagger} \psi r^{2} \sin \theta \mathrm{d} \varphi \mathrm{d} \theta \mathrm{d} r=4 \pi \int_{0}^{\infty}(\hat{A} \psi)^{\dagger} \psi r^{2} \mathrm{~d} r $$
Imply $\hat{p}_r$ , then we have: $$ \begin{aligned} & 4 \pi \int_{0}^{\infty}\left(-i \hbar\left(\frac{\partial}{\partial r}+\frac{1}{r}\right) \psi\right)^{\dagger} \psi r^{2} \mathrm{~d} r \\ =& 4 \pi i \hbar \int_{0}^{\infty}\left(\frac{1}{r} \frac{\partial}{\partial r}(r \psi)\right)^{\dagger} \psi r^{2} \mathrm{~d} r \\ =& 4 \pi i \hbar \int_{0}^{\infty}\left(\frac{\partial}{\partial r}(r \psi)\right)^{\dagger}(\psi r) \mathrm{d} r \\ =& 4 \pi i h\left(\left.(r \psi)^{\dagger}(\psi r)\right|_{0} ^{\infty}-\int_{0}^{\infty}(r \psi)^{\dagger} \frac{\partial}{\partial r}(\psi r) \mathrm{d} r\right) \end{aligned} $$
When proving $\hat{p}$ to be an hermitian, the first term now vanish because there's no extra $r$ in the equation and the wave function equals 0 when it goes to $\pm \infty$.But now, this step cannot be used here.
So how to prove it?