Infinite square well – change in well size

Required math: calculus

Required physics: Schrödinger equation

Reference: Griffiths, David J. (2005), Introduction to Quantum Mechanics, 2nd Edition; Pearson Education – Chapter 2, Post 38.

Suppose a particle is in the ground state in an infinite square well in the interval {x\in[0,a]}. If the right wall suddenly moves to {x=2a}, what effect does this have on the allowable energies?

The ground state for an infinite square well of width {a} is

\displaystyle \psi_{1}=\sqrt{\frac{2}{a}}\sin\left(\frac{\pi x}{a}\right)

The stationary states for a well of width {2a} are

\displaystyle \psi_{n}=\frac{1}{\sqrt{a}}\sin\left(\frac{n\pi x}{2a}\right)

The initial wave function can be expanded as a series using the new states:

\displaystyle \psi_{1}=\sum_{n=1}^{\infty}c_{n}\psi_{n}

where the {c_{n}} are found using the orthogonal property of the stationary states:

\displaystyle c_{n}=\frac{1}{\sqrt{a}}\sqrt{\frac{2}{a}}\int_{0}^{a}\sin\left(\frac{\pi x}{a}\right)\sin\left(\frac{n\pi x}{2a}\right)dx

For {n\ne2}, this integral evaluates to

\displaystyle c_{n}=\frac{4\sqrt{2}}{\pi\left(4-n^{2}\right)}\sin\left(\frac{n\pi}{2}\right)

For {n=2}, we have

\displaystyle c_{2} \displaystyle = \displaystyle \frac{1}{\sqrt{a}}\sqrt{\frac{2}{a}}\int_{0}^{a}\sin^{2}\left(\frac{\pi x}{a}\right)dx
\displaystyle \displaystyle = \displaystyle \frac{1}{\sqrt{2}}

The first few probabilities for stationary states are then

\displaystyle \left|c_{1}\right|^{2} \displaystyle = \displaystyle \frac{32}{9\pi^{2}}
\displaystyle \displaystyle = \displaystyle 0.36
\displaystyle \left|c_{2}\right|^{2} \displaystyle = \displaystyle 0.5
\displaystyle \left|c_{3}\right|^{2} \displaystyle = \displaystyle \frac{32}{25\pi^{2}}
\displaystyle \displaystyle = \displaystyle 0.13

All higher coefficients are smaller, with all even coefficients from 4 upwards being zero.

The most probable energy is therefore

\displaystyle E_{2} \displaystyle = \displaystyle \frac{2^{2}\pi^{2}\hbar^{2}}{2m\left(2a\right)^{2}}
\displaystyle \displaystyle = \displaystyle \frac{\pi^{2}\hbar^{2}}{2ma^{2}}

This happens to be the ground state
of the original well, so, unsurprisingly, the most probable energy for the particle is the one it starts with.

The next most probable energy is

$latex \displaystyle E_{1}=\frac{\pi^{2}\hbar^{2}}{8ma^{2}}

The average energy could be found if we could sum the infinite series

\displaystyle \left\langle E\right\rangle \displaystyle = \displaystyle \frac{1}{2}E_{2}+\sum_{\mathrm{odd}\; n}\frac{32}{\pi^{2}\left(4-n^{2}\right)^{2}}E_{n}
\displaystyle \displaystyle = \displaystyle \frac{\pi^{2}\hbar^{2}}{8ma^{2}}\left[2+\sum_{\mathrm{odd}\; n}\frac{32n^{2}}{\pi^{2}\left(4-n^{2}\right)^{2}}\right]

Clearly this isn’t an easy sum, but we can calculate the expectation value of the Hamiltonian directly as an integral to get the result

\displaystyle \left\langle E\right\rangle \displaystyle = \displaystyle \int_{0}^{a}\psi_{1}(x)\left(-\frac{\hbar^{2}}{2m}\right)\frac{d^{2}}{dx^{2}}\psi_{1}(x)dx
\displaystyle \displaystyle = \displaystyle -\frac{\hbar^{2}}{ma}\int_{0}^{a}\sin\left(\frac{\pi x}{a}\right)\frac{d^{2}}{dx^{2}}\sin\left(\frac{\pi x}{a}\right)dx
\displaystyle \displaystyle = \displaystyle \frac{\hbar^{2}\pi^{2}}{ma^{3}}\int_{0}^{a}\sin^{2}\left(\frac{\pi x}{a}\right)dx
\displaystyle \displaystyle = \displaystyle \frac{\hbar^{2}\pi^{2}}{2ma^{2}}
\displaystyle \displaystyle = \displaystyle E_{2}

If you’re interested, a side effect of this result is that we can evaluate the sum above:

\displaystyle \sum_{\mathrm{odd}\; n}\frac{32n^{2}}{\pi^{2}\left(4-n^{2}\right)^{2}}=2

This series doesn’t converge all that quickly (it goes as {1/n^{2}}). The sum of the first 100 terms gives 1.991893832.

A slight rearrangement gives yet another way of calculating {\pi}:

\displaystyle \sum_{\mathrm{odd}\; n}\frac{16n^{2}}{\left(4-n^{2}\right)^{2}}=\pi^{2}

2 thoughts on “Infinite square well – change in well size

  1. G

    Hi, thanks for your explanation! But I have got a few questions regarding this, what if it is a symmetric well of -a/2<x<a/2 and it is expanded suddenly to -3a/2<x<3a/2 ? how will the probability of its energy state vary? Say, it is initially in the ground state.

    1. growescience

      If all you’re interested in are the energy levels, the location of the well doesn’t matter, so you could analyze your problem with a well that starts out at $latex {x\in\left[0,a\right]}&fg=000000&bg=e5e4e8$ and then expands to $latex {x\in\left[0,3a\right]}&fg=000000&bg=e5e4e8$. You’d then just have to run through the same calculation replacing $latex {2a}&fg=000000&bg=e5e4e8$ by $latex {3a}&fg=000000&bg=e5e4e8$. I don’t think that would give rise to any particular difficulties, but I haven’t checked it.


Leave a Reply

Your email address will not be published. Required fields are marked *

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong>