# previous index next

# 5. Electron in a Box

*Michael Fowler, University of Virginia *

### Plane Wave Solutions

The best way to gain understanding of Schrödinger’s equation is to solve it for various potentials. The simplest is a one-dimensional “particle in a box” problem. The appropriate potential is $V\left(x\right)=0$ for $x$ between $0,L$ and $V\left(x\right)=\infty $ otherwise$\u2014$that is to say, there are infinitely high walls at $x=0$ and $x=L,$ and the particle is trapped between them. This turns out to be quite a good approximation for electrons in a long molecule, and the three-dimensional version is a reasonable picture for electrons in metals.

Between $x=0$ and $x=L$ we have $V\left(x\right)=0,$ so the wave equation is just

$i\hslash \frac{\partial \psi (x,t)}{\partial \text{}\text{\hspace{0.05em}}\text{}\text{}t}=-\frac{{\hslash}^{2}}{2m}\frac{\partial {}^{2}\psi (x,t)}{\partial \text{\hspace{0.05em}}{x}^{2}}$.

A possible plane wave solution is

$\psi (x,t)=A{e}^{\frac{i}{\hslash}(px-Et)}$.

On inserting this into the zero-potential Schrödinger equation above we find $E={p}^{2}/2m,$ as we expect.

It is very important to notice
that the complex conjugate, proportional to ${e}^{-\text{\hspace{0.17em}}\frac{i}{\hslash}(px-Et)}$, is *not* a solution to the
Schrödinger equation! If we blindly put
it into the equation we get

$E=-{p}^{2}/2m,$

an unphysical result.

However, a wave function
proportional to ${e}^{-\text{\hspace{0.17em}}\frac{i}{\hslash}(px+Et)}$ gives *E* = *p*^{2}/2*m*,
so this plane wave *is* a solution to
the equation.

Therefore, the two allowed plane-wave solutions to the zero-potential Schrödinger equation are proportional to ${e}^{\frac{i}{\hslash}(px-Et)}$ and ${e}^{-\text{\hspace{0.17em}}\frac{i}{\hslash}(px+Et)}$ respectively.

**Note that these two solutions have the same time
dependence **${e}^{-\text{\hspace{0.17em}}\frac{iEt}{\hslash}}$.

To decide on the appropriate
solution for our problem of an electron in a box, of course we have to bring in
the walls$\u2014$what they mean is that $\psi \left(x,t\right)=0$ for $x<0$ and for $x>L$ because remember ${\left|\psi \left(x,t\right)\right|}^{2}$ tells us the probability of finding the
particle anywhere, and, since it’s in the box, it’s trapped *between* the
walls, so there’s zero probability of finding it outside.

The condition $\psi =0$ at $x=0$ and $x=L$ reminds us of the vibrating string with two fixed ends$\u2014$the solution of the string wave equation is standing waves of sine form. In fact, taking the difference of the two permitted plane-wave forms above gives a solution of this type:

$\psi (x,t)=A\mathrm{sin}\frac{px}{\hslash}{e}^{-\text{}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\frac{iEt}{\hslash}}$.

This wave function satisfies the Schrödinger equation between the walls, it vanishes at the $x=0$ wall, it will also vanish at $x=L$ provided that the momentum variable satisfies:

$\frac{pL}{\hslash}=\pi ,2\pi ,3\pi ,\dots $

Thus the allowed values of $p$ are $hn/2L,$ where $n=1,2,3\dots $, and from $E={p}^{2}/2m$ the allowed energy levels of the particle are:

$E=\frac{{p}^{2}}{2m}=\frac{1}{2m}{\left(\frac{h}{2L}\right)}^{2},\text{}\frac{4}{2m}{\left(\frac{h}{2L}\right)}^{2},\text{}\frac{9}{2m}{\left(\frac{h}{2L}\right)}^{2},\mathrm{...}$

Note that these energy levels become more and more widely spaced out at high energies, in contrast to the hydrogen atom potential. (As we shall see, the harmonic oscillator potential gives equally spaced energy levels, so by studying how the spacing of energy levels varies with energy, we can learn something about the shape of the potential.)

What about the overall multiplicative constant $A$ in the wave function? This can be real or complex. To find its value, note that at a fixed time, say $t=0,$ the probability of the electron being between $x$ and $x+dx$ is ${\left|\psi \left(x,t\right)\right|}^{2}dx$ or

${\left|A\right|}^{2}{\mathrm{sin}}^{2}\frac{px}{\hslash}dx.$

The *total *probability of
the particle being *somewhere* between $0,L$ must be unity:

$\underset{x=0}{\overset{x=L}{\int}}|A{|}^{2}{\mathrm{sin}}^{2}\frac{px}{\hslash}}dx=1,\text{so}\frac{1}{2}L|A{|}^{2}=1.$

Hence

$\psi (x,t)=\sqrt{\frac{2}{L}}\mathrm{sin}\frac{px}{\hslash}{e}^{-\text{\hspace{0.17em}}\frac{iEt}{\hslash}}\text{\hspace{0.17em}}.$

When $A$ is fixed in this way, by demanding that the
total probability of finding the particle somewhere be unity, it is called the *normalization
constant*.

### Stationary States

Notice that at a later time the probability distribution for the wave function

$\psi (x,t)=\sqrt{\frac{2}{L}}\mathrm{sin}\frac{px}{\hslash}{e}^{-\text{\hspace{0.17em}}\frac{iEt}{\hslash}}$

is *the same*, because *time
only appears as a phase factor* in this time-dependent function, and so does
not affect ${\left|\psi \right|}^{2}.$

*A state with a time-independent probability distribution is called
a* ** stationary state**.

### States with Moving Probability Distributions

Recall that the Schrodinger equation is a *linear* equation, and *the sum of any two solutions is also a
solution* to the equation. That means
that we can add two solutions having *different*
energies, and still have a legal wave function.
We shall establish that in this case, *the probability distribution varies in time*.

The simplest way to see how this must be is to look at an example. Let’s add the ground state to the first excited state, and normalize the sum:

$\begin{array}{l}\psi (x,t)=\sqrt{\frac{1}{L}}\left(\mathrm{sin}\frac{\pi x}{L}{e}^{-\text{\hspace{0.17em}}\frac{i\pi ht}{4m{L}^{2}}}+\mathrm{sin}\frac{2\pi x}{L}{e}^{-\text{\hspace{0.17em}}\frac{i\pi ht}{m{L}^{2}}}\right)\\ \text{note:}h\text{,not}\hslash .\end{array}$.

(You can check the normalization
constant at $t=0.$ ) For
general $x,$ the two terms in the bracket rotate in the
complex plane at different rates, so their sum has a time-varying magnitude.
That is to say, ${\left|\psi \left(x,t\right)\right|}^{2}$ varies in time, so the particle must be moving
around$\u2014$this is *not *a stationary
state.

*Exercise*: To see this, note that at $t=0$ the wave function is:

$\psi (x,0)=\sqrt{\frac{1}{L}}\left(\mathrm{sin}\frac{\pi x}{L}+\mathrm{sin}\frac{2\pi x}{L}\right)$

and sketch this function: the particle is more likely to be found in the left-hand half of the box.

Now, suppose the time is $t=4m{L}^{2}/h,$ so ${e}^{-\text{\hspace{0.17em}}\frac{i\pi ht}{4m{L}^{2}}}=-1.$ At this time,

$\psi (x,2{L}^{2}/h)=\sqrt{\frac{1}{L}}\left(-\mathrm{sin}\frac{\pi x}{L}+\mathrm{sin}\frac{2\pi x}{L}\right)$

and it’s easy to see that the
particle is more likely to be found in the *right*-hand
half.

That is to say, this wave
function, a linear sum of wave functions corresponding to different energies,
has a probability distribution that sloshes back and forth in the box: and, any
attempt to describe a classical-type particle motion, bouncing back and forth, *necessarily* involves adding quantum wave
functions of different energies. Note
that the frequency of the sloshing motion depends on the *difference* of the two energies: how constructively the two
components interfere depends on the difference of the phases in the energies at
the time. A single energy wave function
always has a static probability distribution.

Of course, the *total*
probability of finding the particle *somewhere* in the box remains unity:
the normalization constant is time-independent.

### The Time-Independent Schrödinger Equation: Eigenstates and Eigenvalues

The only way to prevent ${\left|\psi \left(x,t\right)\right|}^{2}$ varying in time is to have all its parts
changing phase in time at the same rate. *This means they all correspond to
the same energy*. If we restrict our considerations to such *stationary
states*, the wave function can be factorized

$\psi (x,t)=\psi (x){e}^{-\text{\hspace{0.17em}}\frac{iEt}{\hslash}}$

and putting this wave function into the Schrödinger equation we find

$-\frac{{\hslash}^{2}}{2m}\frac{d{\text{}}^{2}\psi (x)}{d\text{}{x}^{2}}+V(x)\psi (x)=E\psi (x).$

This is the *time-independent
Schrödinger equation*, and its solutions are the spatial wave functions for
stationary states, states of definite energy. These are often called*
eigenstates* of the equation.

The values of energy
corresponding to these eigenstates are called the *eigenvalues*.

### An Important Point: What, Exactly, Happens at the Wall?

Consider again the wavefunction for the lowest energy state
of a particle confined between walls at $x=0$ and $x=L.$ The reader should sketch the wavefunction from
some point to the left of $x=0$ over to the right of $x=L.$ To the left of $x=0,$ the wavefunction is exactly zero, then at $x=0$ it takes off to the right (inside the box) as
a sine curve. In other words, at the origin the slope of the wavefunction $\psi $ is zero to the left, nonzero to the right. *There
is a discontinuity in the slope at the origin: this means the second derivative
of **$\psi \left(x\right)$** is *infinite* at the origin. *On
examining the time-independent Schrödinger equation above, we see the equation
can *only* be satisfied at the origin because the potential becomes
infinite there$\u2014$the wall
is an infinite potential. (And, in fact, since $\psi $ becomes zero on approaching the origin from
inside the box, the limit must be treated carefully.)

It now becomes obvious that if the box does not have infinite walls, but merely high ones, $\psi $ describing a confined particle cannot suddenly go to zero at the walls: the second derivative must remain finite. For non-infinite walls, $\psi $ and its derivative must be continuous on entering the wall. This has the important physical consequence that $\psi $ will be nonzero at least for some distance into the wall, even if classically the confined particle does not have enough energy to “climb the wall”. (Which it doesn’t, if it’s confined.) Thus, in quantum mechanics, there is a non-vanishing probability of finding the particle in a region which is “classically forbidden” in the sense that it doesn’t have enough energy to get there.