# 1.5: Electron in a Box

- Page ID
- 1647

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)## 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(x) = 0\) for \(x\) between 0, \(L\) and \(V(x) = \infty\) otherwise—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 = 0\), so the wave equation is just

\[ i\hbar \dfrac{\partial \psi(x,t)}{\partial t} =-\dfrac{\hbar^2}{2m} \dfrac{\partial^2 \psi (x,t)}{\partial x^2}. \tag{1.5.1} \]

A possible plane wave solution is

\[ \psi (x,t) =Ae^{\dfrac{i}{\hbar}(px-Et)}. \tag{1.5.2} \]

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^{\dfrac{i}{\hbar}(px-Et)}\), is *not* a solution to the Schrödinger equation! If we blindly put it into the equation we get

\[E=\dfrac{-p^2}{2m}\]

However, this is an unphysical result.

However, a wave function proportional to \(e^{-\dfrac{i}{\hbar}(px-Et)}\) gives \(E=p^2/2m\), 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^{\dfrac{i}{\hbar}(px-Et)}\) and \(e^{-\dfrac{i}{\hbar}(px-Et)}\) respectively.

**Note that these two solutions have the same time dependence ** \(e^{-\dfrac{iEt}{\hbar}}\).

To decide on the appropriate solution for our problem of an electron in a box, of course we have to bring in the walls—what they mean is that \(\psi=0\) for \(x < 0\) and for \(x > L\) because remember \(|\psi|^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—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\sin \dfrac{px}{\hbar}e^{-\dfrac{iEt}{\hbar}} \tag{1.5.3} \].

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:

\[ \dfrac{pL}{\hbar}=\pi, \, 2\pi, \, 3\pi... \tag{1.5.4} \]

Thus the allowed values of \(p\) are \(hn/2L\), where \(n = 1, 2, 3…\), and from \(E=p^2/2m\) the allowed energy levels of the particle are:

\[ E=\dfrac{p^2}{2m}=\dfrac{1}{2m}\left(\dfrac{h}{2L}\right)^2 \, ,\; \dfrac{4}{2m}\left(\dfrac{h}{2L}\right)^2 \, ,\; \dfrac{9}{2m}\left(\dfrac{h}{2L}\right)^2 \, ,... \tag{1.5.5} \]

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 \(|\psi|^2dx\) or

\[ |A|^2\sin^2 \dfrac{px}{\hbar}dx \tag{1.5.6} \]

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

\[ |\psi|^2dx\]

or

\[ \int\limits_{x=0}^{x=L} |A|^2\sin^2 \dfrac{px}{\hbar}dx=1, \: so\: \dfrac{1}{2}L|A|^2=1 \tag{1.5.7} \]

Hence

\[ \psi(x,t)=\sqrt{\dfrac{2}{L}}\sin \dfrac{px}{\hbar}e^{-\dfrac{iEt}{\hbar}} \tag{1.5.8} \]

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}}\sin \frac{px}{\hbar}e^{-\frac{iEt}{\hbar}} \tag{1.5.9} \]

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

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

## States with Moving Probability Distributions

Recall that the Schrödinger 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: \[ \psi(x,t)=\sqrt{\frac{1}{L}} \left( \sin\frac{\pi x}{L}e^{-\frac{i\pi ht}{4mL^2}}+\sin\frac{2\pi x}{L}e^{-\frac{i\pi ht}{mL^2}} \right) \tag{1.5.10} \]

\[ note: \; h, \; not \; \hbar. \]

(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, \(|\psi(x,t)|^2\) varies in time, so the particle must be moving around—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( \sin\frac{\pi x}{L}+\sin\frac{2\pi x}{L} \right) \tag{1.5.11} \]

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=4mL^2/h\), so \(e^{-\frac{i\pi ht}{4mL^2}}=-1\). At this time, \[ \psi(x,2L^2/h)=\sqrt{\frac{1}{L}} \left( -\sin\frac{\pi x}{L}+\sin\frac{2\pi x}{L} \right) \tag{1.5.11} \]

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 \(|\psi(x,t)|^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^{-\frac{iEt}{\hbar}} \tag{1.5.12} \]

and putting this wave function into the Schrödinger equation we find \[ -\frac{\hbar^2}{2m} \frac{d^2 \psi(x)}{dx^2} +V(x)\psi(x)=E\psi(x) \tag{1.5.13} \]

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\)* 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—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.

## Contributor

- Michael Fowler (Beams Professor, Department of Physics, University of Virginia)