Expectation Value of Momentum in H-Atom

[quantumdude]

Here's a silly question. I'm sure I should know the answer, but alas my most recent QM course was 9 years ago.

I sat down to calculate the expectation value of momentum in the H-atom today, because some kid on another forum wanted to know how fast an electron in an atom is. I was going to calculate <p> and divide it by m for him.

So I wrote down <p>=\int_V\psi_{100}^*(-i\hbar\nabla)\psi_{100} d^3x (\psi_{100} being the ground state H-atom wavefunction). Upon noticing that I've got a solitary i sandwiched between two purely real functions, I didn't even bother to continue. It's going to come out imaginary.

What am I doing wrong?

 

[StatMechGuy]

Not evaluating the integral. Just compute it and see.

Then, calculate <p^2> and take its root mean square and tell the kid that.

 

[Meir Achuz]

The vector momentum =0.

 

[vanesch]

So I wrote down &lt;p&gt;=\int_V\psi_{100}^*(-i\hbar\nabla)\psi_{100} d^3x (\psi_{100} being the ground state H-atom wavefunction). Upon noticing that I've got a solitary i sandwiched between two purely real functions, I didn't even bother to continue. It's going to come out imaginary.

What am I doing wrong?

It should come out 0, and in fact you just proved this:
the p operator being hermitean, its expectation value must be real for any vector. You showed that it must be a purely imaginary number. Hence it can only be 0.

And that's normal too: the expectation of momentum for any stationary bound state must be 0: given that a stationary state will not have the momentum change in time, and given that the fact that it is bound will confine (most of) the wavefunction to a bounded region, the C.O.G. cannot leave that bounded region, which it should, if there was a non-zero expectation value of momentum, after some time.

 

[quantumdude]

The integral is not zero. The ground state wavefunction is:

\psi_{100}(r)=\frac{1}{\sqrt{\pi}}\frac{1}{a_0^{3/2}}\exp(-r/a_0)

If we integrate \psi^*_{100}(-i\hbar\nabla)\psi_{100} over all space we get:

\frac{-i\hbar}{\pi a_0^3}\int_0^{2\pi}\int_{-1}^1\int_0^{\infty}\exp(-r/a_0)\frac{d}{dr}\exp(-r/a_0)r^2drd(\cos(\theta))d\phi

Since the wavefunction doesn't depend on the angular coordinates, this integral is proportional to the following:

i\int_0^{\infty}r^2\exp(-2r/a_0)dr

The integrand is not identically zero, and is nonnegative everywhere. So it's integral cannot vanish. So, it seems this result is purely imaginary and not zero.

There's got to be a mistake somewhere, but I don't see where.

 

[smallphi]

The nabla is a vector operator. It produces something proportional to r_vector. Integrating that over spherically symmetrical space gives zero vector.

 

[nrqed]

The integral is not zero. The ground state wavefunction is:

\psi_{100}(r)=\frac{1}{\sqrt{\pi}}\frac{1}{a_0^{3/2}}\exp(-r/a_0)

If we integrate \psi^*_{100}(-i\hbar\nabla)\psi_{100} over all space we get:

\frac{-i\hbar}{\pi a_0^3}\int_0^{2\pi}\int_{-1}^1\int_0^{\infty}\exp(-r/a_0)\frac{d}{dr}\exp(-r/a_0)r^2drd(\cos(\theta))d\phi

Since the wavefunction doesn't depend on the angular coordinates, this integral is proportional to the following:

i\int_0^{\infty}r^2\exp(-2r/a_0)dr

The integrand is not identically zero, and is nonnegative everywhere. So it's integral cannot vanish. So, it seems this result is purely imaginary and not zero.

There's got to be a mistake somewhere, but I don't see where.

I agree with Smallphi.
You forgot to include \hat{r}. Integrating this over all angles gives zero.

 

[olgranpappy]

One simple way to get an appropriate estimate of the velocity is to divide 13.6eV by the electron mass, multiply by 2 and then take the square-root thus finding that the electron travels at a bit less than one hundredth of the speed of light. This is correct because the expectation value of the kinetic energy is equal to the negative of the expectation value of the Hamiltonian for a Coulomb potential. I.e. -(-13.6eV).

Also, you can see that <\vec p> must be zero for any real wave-function (btw, eigenfunctions of this Hamiltonian can always be *chosen* as real) by just integrating by parts... you end up with the equation <\vec p>=-<\vec p> which implies that <\vec p>=0

 

[Hans de Vries]

To get rid of the " i " you might want to use following momentum
density operator:

{ i\hbar \over 2} \left( \frac{\partial \psi^*}{\partial x_i}\psi\ -\ \psi^*\frac{\partial \psi}{\partial x_i }\right)

Where the psi is the relativistic wave function related to the Schrödinger
wave function as follows:

\psi\ =\ \Psi \ e^{-imc^2t/\hbar}, \qquad \mbox{$\Psi$ is the Schroedinger wave function}

This gives the same result without the imaginary part. It's always real.
This is the Klein-Gordon momentum density operator, as well as the Dirac
momentum density operator, (after the spin contributions are split off
via the Gordon decomposition).Regards, Hans

 

[Hans de Vries]

divide 13.6eV by the electron mass, multiply by 2 and then take the square-root thus finding that the electron travels at a bit less than one hundredth of the speed of light.

A little bit less than one hundredth of the speed of light is 1/137.036 to be
a bit more exact, a value also known as alpha.


Regards, Hans