The infinite square well (particle in a box)

Required math: calculus

Required physics: time-independent Schrödinger equation

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

Shankar, R. (1994), Principles of Quantum Mechanics, Plenum Press. Section 5.2, Exercise 5.2.5.

To get the feel of how to solve the time-independent Schrödinger equation in one dimension, the most commonly used example is that of the infinite square well, sometimes known as the ‘particle in a box’ problem. First, recall the Schrödinger equation itself:

\displaystyle   -\frac{\hbar^{2}}{2m}\frac{\partial^{2}\Psi}{\partial x^{2}}+V(x)\Psi \displaystyle  = \displaystyle  i\hbar\frac{\partial\Psi}{\partial t} \ \ \ \ \ (1)

Remember that the ‘time-independent’ bit refers to the potential function {V} which is taken to be a function of position only; the wave function itself, which is the solution of the equation, will in general be time-dependent.

The infinite square well is defined by a potential function as follows:

\displaystyle  V(x)=\begin{cases} 0 & 0<x<a\\ \infty & \mathrm{otherwise} \end{cases} \ \ \ \ \ (2)

An area with an infinite potential means simply that the particle is not allowed to exist there. In classical physics, you can think of an infinite potential as an infinitely high wall, which no matter how much kinetic energy a particle has, it can never leap over. Although examples from classical physics frequently break down when applied to quantum mechanics, in this case, the comparison is still valid: an infinitely high potential barrier is an absolute barrier to a particle in both cases.

We saw in our study of the time independent Schrödinger equation that separation of variables reduces the problem to solving the spatial part of the equation, which is

\displaystyle   -\frac{\hbar^{2}}{2m}\frac{d^{2}\psi}{dx^{2}}+V(x)\psi \displaystyle  = \displaystyle  E\psi \ \ \ \ \ (3)

where the constant {E} represents the possible energies that the system can have. It is important to note that both {E} and {V(x)} are unknown before we solve the equation. In classical physics, we would be allowed to specify {E} since it is just the kinetic energy that the particle has inside the well. Classically, {E} can be any positive quantity, and the particle would just bounce around inside the well without ever changing its speed (assuming the walls were perfectly elastic and there was no friction). In quantum physics, as we will see, {E} can have only certain discrete values, and these values arise in the course of solving the equation.

In an infinite square well, the infinite value that the potential has outside the well means that there is zero chance that the particle can ever be found in that region. Since the probability density for finding the particle at a given location is {|\Psi|^{2}}, this condition can be represented in the mathematics by requiring {\psi(x)=0} if {x<0} or {x>a}. This condition is forced from the Schrödinger equation since if {V(x)=\infty}, any non-zero value for {\psi(x)} would result in an infinite term in the equation. However, it is certainly not rigorous mathematics, since multiplying infinity by zero can be done properly only by using a limiting procedure, which we haven’t done here. A proper treatment of the infinite square well is as a limiting case of the finite square well, where {V(x)} can have a large but finite value outside the well. However, the mathematics for solving the finite square well is considerably more complicated and tends to obscure the physics. Readers who are worried, however, can be reassured that the energy levels in the finite square well do become those in the infinite square well when the proper limit is taken.

Inside the well, {V(x)=0} so the Schrödinger equation becomes

\displaystyle   -\frac{\hbar^{2}}{2m}\frac{d^{2}\psi}{dx^{2}} \displaystyle  = \displaystyle  E\psi\ \ \ \ \ (4)
\displaystyle  \frac{d^{2}\psi}{dx^{2}} \displaystyle  = \displaystyle  -k^{2}\psi\ \ \ \ \ (5)
\displaystyle  \mathrm{with\;}k \displaystyle  \equiv \displaystyle  \frac{\sqrt{2mE}}{\hbar} \ \ \ \ \ (6)

This differential equation has the general solution

\displaystyle   \psi(x) \displaystyle  = \displaystyle  A\sin(kx)+B\cos(kx) \ \ \ \ \ (7)

for unspecified (yet) constants {A} and {B}. If you don’t believe this, just substitute the solution back into the equation.

How can we determine {A} and {B}? To do this, we need to appeal to Born’s conditions on the wave function. Born’s first condition is clearly satisfied here: {\psi} is single-valued. The second condition of {\psi} being square integrable we’ll leave for a minute. The third condition is that {\psi} must be continuous. We have argued above that {\psi=0} outside the well, so in particular, this means that at the boundaries {x=0} and {x=a} we must have {\psi=0}. If we impose that condition on our general solution above, we get:

\displaystyle   \psi(0) \displaystyle  = \displaystyle  0\ \ \ \ \ (8)
\displaystyle  A\sin(0)+B\cos(0) \displaystyle  = \displaystyle  0\ \ \ \ \ (9)
\displaystyle  B \displaystyle  = \displaystyle  0\ \ \ \ \ (10)
\displaystyle  \psi(a) \displaystyle  = \displaystyle  0\ \ \ \ \ (11)
\displaystyle  A\sin(ka) \displaystyle  = \displaystyle  0\ \ \ \ \ (12)
\displaystyle  ka \displaystyle  = \displaystyle  n\pi\ \ \ \ \ (13)
\displaystyle  \frac{\sqrt{2mE}}{\hbar} \displaystyle  = \displaystyle  \frac{n\pi}{a}\ \ \ \ \ (14)
\displaystyle  E \displaystyle  = \displaystyle  \frac{n^{2}\pi^{2}\hbar^{2}}{2ma^{2}} \ \ \ \ \ (15)

where {n} is an integer. The useful values of {n} are just the positive integers. To see this, note that if {n=0}, then {\psi(x)=0} everywhere which besides being very boring, is also no good as a probability density since its square modulus cannot integrate to 1. Negative integers don’t really give new solutions, since {\sin(-x)=-\sin x,} so the negative sign can be absorbed into the (still undetermined) constant {A}. Also, the energies depend only on the square of {n} so the sign of {n} doesn’t matter physically.

We can now return to the square-integrable condition and use it to determine {A}. Remember that the integral is over all space in which the particle can be found, so in this case we are interested in {0\le x\le a}.

\displaystyle   \int_{0}^{a}|\psi|^{2}dx \displaystyle  = \displaystyle  1\ \ \ \ \ (16)
\displaystyle  A^{2}\int_{0}^{a}\sin^{2}kx \displaystyle  = \displaystyle  1\ \ \ \ \ (17)
\displaystyle  \displaystyle  = \displaystyle  A^{2}\frac{1}{2}\int_{0}^{a}\left[1-\cos\left(2kx\right)\right]dx\ \ \ \ \ (18)
\displaystyle  \displaystyle  = \displaystyle  \frac{A^{2}}{2}\left[x-\frac{1}{2k}\sin\left(2kx\right)\right]_{0}^{a}\ \ \ \ \ (19)
\displaystyle  \displaystyle  = \displaystyle  A^{2}\frac{a}{2}\ \ \ \ \ (20)
\displaystyle  A \displaystyle  = \displaystyle  \pm\sqrt{\frac{2}{a}} \ \ \ \ \ (21)

where we used 13 in 19 to eliminate the sine term.

Since it is only the square modulus of the wave function that has physical significance, we can ignore the negative root and take the final form of the wave function as

\displaystyle   \psi(x) \displaystyle  = \displaystyle  \sqrt{\frac{2}{a}}\sin\frac{n\pi x}{a} \ \ \ \ \ (22)

Notice what has happened here. Applying the first boundary condition at {x=0} allowed us to eliminate {B}. But the other boundary condition at {x=a} ended up giving us a condition on {E} rather than {A}. (Well, ok, we could have used the second boundary condition to set {A=0} but then we would have {\psi(x)=0} everywhere again.) Not only that, but the energy levels are discrete; thus the infinite square well is the first case in which the Schrödinger equation has actually predicted quantization in a system.

So the boundary conditions on the differential equation have put restrictions on the allowable energies. The acceptable solutions for {\psi} are determined by the condition {ka=n\pi} and so the various {\psi} functions are just lobes of the sine function. The lowest energy, called the ground state, occurs when {n=1} and {\psi} is half a sine wave, consisting of the bit between {x=0} and {x=\pi}. The next state at {n=2} corresponds to a single complete cycle of the sine wave; {n=3} contains 1.5 cycles and so on.

Eagle-eyed readers will have noticed that in all the excitement over discovering quantization, we have neglected to look at Born’s fourth condition: that of continuous first derivatives. Clearly this condition is violated, since the derivative of {\psi} outside the well is 0 (since {\psi=0} outside the well), but inside the derivative is {kA\cos(kx),}which is {kA} at {x=0} and {\pm kA} at {x=a} (the sign depends on whether {n} is odd or even). Neither of these derivatives is zero.

The reason is, of course, because of the infinite potential function which is not physically realistic. In fact, what happens in the (real-world) finite square well is that the wave function inside the potential barrier (that is, just off the ends of the well) is not zero, but a decaying exponential which tends to zero the further into the barrier you go. In that case, it is possible to make both the wave function and its first derivative continuous at both ends of the well (and it is precisely that condition which makes the mathematics so much more complicated in the finite square well).

In fact, this effect happens in any potential where the energy of the particle is less than that of a (finite) potential barrier: the particle’s wave function extends into the barrier region. So does that mean that the particle has a probability of appearing inside a barrier? Technically yes, but in practice it usually doesn’t do the particle much good, since the probability of being outside the barrier is usually a lot greater. However, there is one case where this barrier penetration effect does occur, and that is if the barrier is thin enough for the wave function to have a significant magnitude on the other side of the barrier. That is, if we have a particle in a finite well, but the wall of the well is fairly thin and there is another well (or just open space) on the other side of the barrier, then the wave function starts off with a respectable magnitude inside the main well, extends into the barrier (but gets attentuated exponentially in doing so), but before the attenuation gets so severe that the wave function becomes very small, it bursts through to the other side of the barrier. That means that, yes, there is a definite probability that the particle can spontaneously appear outside the well without having to jump over the barrier. In effect, it tunnels through the barrier and escapes. The effect, not surprisingly, is known as quantum tunneling and is one of the main causes of some forms of radioactive decay. But that’s a topic for another post.

57 thoughts on “The infinite square well (particle in a box)

  1. Pingback: The time-independent Schrödinger equation – general solutions « Physics tutorials

  2. Pingback: The wave function – making energy measurements « Physics tutorials

  3. Pingback: The free particle « Physics tutorials

  4. Pingback: Harmonic oscillator – summary « Physics tutorials

  5. Pingback: Energy states: bound and scattering states « Physics tutorials

  6. Pingback: Delta-function well – bound state « Physics tutorials

  7. Pingback: Finite square well – bound states & even wave functions « Physics tutorials

  8. Pingback: Wave function: Born’s conditions « Physics tutorials

  9. Pingback: Index – Physics – Quantum mechanics « Physics tutorials

  10. Pingback: The Schrödinger equation « Physics tutorials

  11. Pingback: Infinite square well – minimum energy « Physics tutorials

  12. Pingback: Infinite square well – uncertainty principle « Physics tutorials

  13. Pingback: Infinite square well – combination of two lowest states « Physics tutorials

  14. Pingback: Infinite square well – phase difference « Physics tutorials

  15. Pingback: Infinite square well – triangular initial state « Physics tutorials

  16. Pingback: Infinite square well – particle in left half « Physics tutorials

  17. Pingback: Infinite square well – average energy « Physics tutorials

  18. Pingback: Complex exponentials and trig functions « Physics tutorials

  19. Pingback: Infinite square well – centered coordinates « Physics tutorials

  20. Pingback: Infinite square well – cubic sine initial state « Physics tutorials

  21. Pingback: Infinite square well – change in well size « Physics tutorials

  22. Pingback: Quantum revival time « Physics tutorials

  23. Pingback: Hybrid infinite-finite square well « Physics tutorials

  24. Pingback: Infinite square well with delta function barrier « Physics tutorials

  25. Pingback: Infinite square well: momentum « Physics tutorials

  26. Pingback: Infinite square well: momentum space wave functions « Physics tutorials

  27. Pingback: Infinite square well in three dimensions « Physics tutorials

  28. Pingback: Infinite spherical well – numerical solutions « Physics tutorials

  29. Pingback: Infinite square well: 2 particle systems | Physics tutorials

  30. Pingback: Exchange force: infinite square well | Physics tutorials

  31. Pingback: Electron gas: a crude model of a solid | Physics tutorials

  32. Pingback: Statistical mechanics in quantum theory: counting states | Physics tutorials

  33. Pingback: First order non-degenerate perturbation theory | Physics tutorials

  34. Pingback: Virial theorem | Physics tutorials

  35. Pingback: Quantum dots | Physics pages

  36. Pingback: WKB approximation | Physics pages

  37. Pingback: Time-dependent perturbation uniform in space | Physics pages

  38. Pingback: Time-dependent perturbation of the infinite square well | Physics pages

  39. Pingback: The adiabatic approximation in quantum mechanics | Physics pages

  40. Pingback: Infinite square well with variable delta function barrier: ground state energy | Physics pages

  41. Pingback: Second quantizing the tight-binding hamiltonian | Physics pages

  42. Pingback: The free particle | Physics pages

  43. Joseph McKeown

    Can somebody clarify how to get from step (18) to step (19) and then step (20) please? I don’t understand how the sin^2 has been integrated to get the expression in (19) and I also don’t understand the apparent simplification from (19) to (20)…

    1. gwrowe Post author

      A very good question; I don’t know what I was thinking when I wrote that. I’ve modified the derivation to (hopefully) make it clearer.

  44. Joseph McKeown

    Okay I now understand your workings, however…
    If one substitutes a=1, x=0.5, it appears to me that the probability of finding the particle becomes 2… which is clearly wrong. Do you have anything on this website involving individual examples of finding probabilities in this context?

    1. gwrowe Post author

      Remember that {\left|\psi\right|^{2}} is a probability density, not a probability. That is, it has meaning only when integrated over a range of {x} values, so you can’t say that {\left|\psi\left(x\right)\right|^{2}} is the probability of finding the particle at a precise location {x}. All you can say is that the probability of finding the particle in the interval {\left[0.5,0.5+dx\right]} is {\left|\psi\left(0.5\right)\right|^{2}dx}. Since {dx} is an infinitesimal, this quantity is definitely less than 1. The condition {\int_{0}^{a}\left|\psi\right|^{2}dx=1} says that the particle has to be found somewhere in the interval, but it doesn’t say where.

  45. hitesh

    two situations were given:-
    if (1) bottom of the well is not at zero potential, i.e., V(x) = V_0 for 0<x<L
    (2) for symmetric well, i.e., V(x)=0 for -L/2 < x < L/2, and V(x) = infinity for x L/2.

    Can you comment on this

  46. Pingback: Multiplicity of a 2-dim ideal gas | Physics pages

  47. Pingback: Continuity of the wave function – Born’s conditions revisited | Physics pages

  48. Pingback: Vector spaces: definitions and examples | Physics pages

  49. Pingback: Vibrating string – normal mode analysis | Physics pages

  50. Pingback: Infinite square well – expanding well | Physics pages

  51. Pingback: Infinite square well – force to decrease well width | Physics pages

  52. Pingback: Fermions and bosons in the infinite square well | Physics pages

  53. D. Krause

    Since infinite potential wells are not realistic, can we say that we *never* can confine a particle in a well? Alternatively, can we say that we cannot completely isolate a particle ever? What about Dehmelt experiments?
    Another question: how to describe this situation in QFT? Are we confining the whole infinite particle-field in the well? In this case, how to explain other particles of the same kind outside the well, say in Saturn?
    Thanks for answering.


Leave a Reply

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