Daily Archives: Wed, 4 January 2017

Free particle in the position basis

References: Shankar, R. (1994), Principles of Quantum Mechanics, Plenum Press. Section 5.1, Exercise 5.1.2.

In quantum mechanics, the free particle has degenerate energy eigenstates for each energy

\displaystyle  E=\frac{p^{2}}{2m} \ \ \ \ \ (1)

where {p} is the momentum. The degeneracy arises because the momentum can be either positive (for a particle moving to the right) or negative (to the left):

\displaystyle  p=\pm\sqrt{2mE} \ \ \ \ \ (2)

Thus the most general energy eigenstate is a linear combination of the two momentum states:

\displaystyle  \left|E\right\rangle =\beta\left|p=\sqrt{2mE}\right\rangle +\gamma\left|p=-\sqrt{2mE}\right\rangle \ \ \ \ \ (3)

This bizarre feature of quantum mechanics means that a particle in such a state could be moving either left or right, and if we make a measurement of the momentum we force the particle into one or other of the two momentum states.

We obtained this solution by working in the momentum basis, but we can also find the solution in the position basis. In that basis, the momentum operator has the form

\displaystyle  P=-i\hbar\frac{d}{dx} \ \ \ \ \ (4)

The matrix elements of this operator in the position basis are

\displaystyle  \left\langle x\left|P\right|x^{\prime}\right\rangle =-i\hbar\delta^{\prime}\left(x-x^{\prime}\right) \ \ \ \ \ (5)

where {\delta^{\prime}\left(x-x^{\prime}\right)} is the derivative of the delta function with respect to the {x}, not the {x^{\prime}}. We can use the properties of this derivative to get a solution in the {X} basis. To be completely formal about it, the derivation of the matrix elements of {P^{2}} in the {X} basis is:

\displaystyle   \left\langle x\left|P^{2}\right|\psi\right\rangle \displaystyle  = \displaystyle  \int\int\left\langle x\left|P\right|x^{\prime}\right\rangle \left\langle x^{\prime}\left|P\right|x^{\prime\prime}\right\rangle \left\langle x^{\prime\prime}\left|\psi\right.\right\rangle dx^{\prime}dx^{\prime\prime}\ \ \ \ \ (6)
\displaystyle  \displaystyle  = \displaystyle  \int\int\left\langle x\left|P\right|x^{\prime}\right\rangle \left(-i\hbar\delta^{\prime}\left(x^{\prime}-x^{\prime\prime}\right)\right)\psi\left(x^{\prime\prime}\right)dx^{\prime}dx^{\prime\prime}\ \ \ \ \ (7)
\displaystyle  \displaystyle  = \displaystyle  -i\hbar\int\left\langle x\left|P\right|x^{\prime}\right\rangle \frac{d\psi\left(x^{\prime}\right)}{dx^{\prime}}dx^{\prime}\ \ \ \ \ (8)
\displaystyle  \displaystyle  = \displaystyle  -i\hbar\int\int\left(-i\hbar\delta^{\prime}\left(x-x^{\prime}\right)\right)\frac{d\psi\left(x^{\prime}\right)}{dx^{\prime}}dx^{\prime}\ \ \ \ \ (9)
\displaystyle  \displaystyle  = \displaystyle  -\hbar^{2}\frac{d^{2}}{dx^{2}}\psi\left(x\right) \ \ \ \ \ (10)

In this basis, the Schrödinger equation is therefore the familiar one:

\displaystyle   \frac{P^{2}}{2m}\left|\psi\right\rangle \displaystyle  = \displaystyle  E\left|\psi\right\rangle \ \ \ \ \ (11)
\displaystyle  \left\langle x\left|\frac{P^{2}}{2m}\right|\psi\right\rangle \displaystyle  = \displaystyle  E\psi\left(x\right)\ \ \ \ \ (12)
\displaystyle  -\frac{\hbar^{2}}{2m}\frac{d^{2}}{dx^{2}}\psi\left(x\right) \displaystyle  = \displaystyle  E\psi\left(x\right)\ \ \ \ \ (13)
\displaystyle  \frac{d^{2}}{dx^{2}}\psi\left(x\right) \displaystyle  = \displaystyle  -\frac{2mE}{\hbar^{2}}\psi\left(x\right) \ \ \ \ \ (14)

This has the general solution

\displaystyle  \psi\left(x\right)=\beta e^{ix\sqrt{2mE}/\hbar}+\gamma e^{-ix\sqrt{2mE}/\hbar} \ \ \ \ \ (15)

[Shankar extracts a factor of {1/\sqrt{2\pi\hbar}} but as he notes, this is arbitrary and can be absorbed into the constants {\beta} and {\gamma} as we’ve done here.]

In this derivation we’ve implicitly assumed that {E>0}, since there is no potential so a free particle can’t really have a negative energy. However, if you follow through the derivation, you’ll see that it works even if {E<0}. In that case, we’d get

\displaystyle  \psi\left(x\right)=\beta e^{-x\sqrt{2m\left|E\right|}/\hbar}+\gamma e^{x\sqrt{2m\left|E\right|}/\hbar} \ \ \ \ \ (16)

That is, the exponents in both terms are now real instead of imaginary. The problem with this is that the first term blows up for {x\rightarrow-\infty} while the second blows up for {x\rightarrow+\infty}. Thus this function is not normalizable, even to a delta function (as was the case when {E>0}), so functions such as these when {E<0} are not in the Hilbert space.