Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I'm having trouble doing it. I know so far that if we have two Hermitian operators $A$ and $B$ that do not commute, and suppose we wish to find the quantum mechanical Hermitian operator for the product $AB$, then


However, if I have to find an operator equivalent for the radial component of momentum, I am puzzled. It does not come out to be simply


where $\vec{r}$ and $\vec{p}$ are the position and the momentum operator, respectively. Where am I wrong in understanding this?

share|cite|improve this question
up vote 6 down vote accepted

You would have to use the fact that the momentum operator in position space is $\vec{p} = -i\hbar\vec{\nabla}$ and use the definition of the gradient operator in spherical coordinates:

$$\vec{\nabla} = \hat{r}\frac{\partial}{\partial r} + \hat{\theta}\frac{1}{r}\frac{\partial}{\partial\theta} + \hat{\phi}\frac{1}{r\sin\theta}\frac{\partial}{\partial\phi}$$

So the radial component of momentum is

$$p_r = -i\hbar\hat{r}\frac{\partial}{\partial r}$$

However: after a bit of investigation prompted by the comments, I found that in practice this is not used very much. It's more useful to have an operator $p_r'$ that satisfies

$$-\frac{\hbar^2}{2m}\nabla^2 R(r) = \frac{p_r'^2}{2m} R(r)$$

This lets you write the radial component of the time-independent Schrödinger equation as

$$\biggl(\frac{p_r'^2}{2m} + V(r)\biggr)R(r) = E R(r)$$

The action of the radial component of the Laplacian in 3D is

$$\nabla^2 R(r) = \frac{1}{r^2}\frac{\partial}{\partial r}\biggl(r^2\frac{\partial R(r)}{\partial r}\biggr)$$

and if you solve for the operator $p'_r$ that satisfies the definition above, you wind up with

$$p'_r = -i\hbar\biggl(\frac{\partial}{\partial r} + \frac{1}{r}\biggr)$$

This is called the "radial momentum operator." Strictly speaking, it is different from the "radial component of the momentum operator," which is, by definition, $p_r$ as I wrote it above, although I wouldn't be surprised to find people mixing up the terminology relatively often.

share|cite|improve this answer
on second thought, your expression is correct for 1 dimensional space. For D dimensional euclidean space the expression is $$p_r = \frac{\hbar}{i}(\frac{\partial}{\partial r}+\frac{D-1}{2r})$$ Your expression is correct in 1D. Please make a trivial edit so I may atleast remove my downvote. – yayu May 1 '11 at 18:04
@yayu: I did a bit of research and made an edit to reflect that. – David Z May 1 '11 at 18:19

I was able to figure it out, so here goes the clarification for the record.

Classically $$p_r=\hat{D_r} = \hat{r}\cdot\hat{p} = \frac{\hbar}{i}\frac{\partial}{\partial r}$$

However $\hat{D_r}$ is not hermitian. Consider the adjoint $$\hat{D_r^+}= \hat{p}\cdot\hat{r} =\frac{\hbar}{i} \left ( \frac{\partial}{\partial r}+\frac{2}{r} \right ) $$

Now we know from linear algebra how to construct a hermitian operator from an operator and its adjoint: $$\hat{p_r} = \frac{\hat{D_r^+}+\hat{D_r}}{2}=\frac{\hbar}{i} \left ( \frac{\partial}{\partial r}+\frac{1}{r} \right ) $$

And btw, for those who followed my initial question, don't do the following calculational mistake that I committed:

$$p_r = \frac{1}{2}\left(\frac{1}{r}(\vec{r}\cdot\vec{p})+(\vec{p}\cdot\vec{r})\frac{1}{r}\right)\neq \frac{1}{2}\frac{1}{r}(\vec{r}\cdot\vec{p}+\vec{p}\cdot\vec{r})$$

share|cite|improve this answer
Dear @yayu: So your initial guess in the question is correct after all! – Qmechanic May 2 '11 at 20:55
@Qmechanic yes, but I'd thought that this was an ad hoc choice and missed the notion of hermiticity which actually guides the construction. – yayu May 7 '11 at 6:02
@Qmechanic Also thanks, the same notation for operator and unit vectors can be regarded as one of the culprits for my initial mistake in calculation. (i.e making the final inequality an equality) – yayu May 7 '11 at 6:05

Even in one dimension the operator $p_r=-i\partial_r$ on the half line $r>0$ has deficiency indices $(0,1)$. There is thus no way to define it it as a self-adjoint operator. In practical terms this abstract mathematical statement means that there is no set of boundary conditions thta we can impose on the wavefunction $\psi(r)$ that lead to a complete set of eigenfunctions for $p_r$. For example, integration by parts to prove formal hermiticity requires that $\psi(0)=0$ but all potential eigenfunctions are of the form $\psi_k(r)=\exp\{ikr\}$ for some real or complex $k$, and no value of $k$ can make $\psi_k(0)$ be zero.

Since one needs eigenfunctions and eignvalues to assign a probability to the outcome of a measurement of $p_r$ this means that $p_r$ is not an ``observable.''

share|cite|improve this answer
You are just explaining why $\mathrm{i}\partial_r$ isn't Hermitian. The question asks how to construct the correct observable corresponding to radial momentum, not why the naive guess doesn't work. – ACuriousMind Nov 5 '15 at 16:02

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.