Physlets run in a Javaenabled browser on the latest Windows & Mac operating systems.
If Physlets do not run, click here for help on updating Java and setting Java security.
Section 10.2: The Quantummechanical Infinite Square Well
Please wait for the animation to completely load.
In the infinite square well potential, a particle is confined to a box of length L by two infinitely high potential energy barriers:
V = ∞ x ≤ 0 , V = 0 0 < x < L , V = ∞ x ≥ L .
We begin with the timeindependent Schrödinger equation in one dimension for 0 < x < L:
−(ħ^{2}/2m)(d^{2}/dx^{2}) ψ(x) = E ψ(x) .
The solution to this differential equation is a combination of sines and cosines (or complex exponentials; here, because of the boundary conditions, we choose sines and cosines). The general solution is
ψ(x) = Asin(kx) + Bcos(kx) ,
where k^{2 }= 2mE/ħ^{2}.
We still need to satisfy the boundary conditions, and thereby determine A, B, and k (and therefore E). For the wall on the left we need ψ(x) = 0, for x ≤ 0 and for the wall on the right we need ψ(x) = 0 for x ≥ L, which means that
ψ(0) = 0 and ψ(L) = 0 .
The cosine part of the general solution does not vanish at the x = 0 boundary since, cos(0) ≠ 0, and therefore we must have B = 0. We are left with ψ(x) = Asin(kx) and determining both A and k. This situation is shown in the animation. The boundary condition at x = 0 is already solved and you can vary the energy to see the effect on the energy eigenfunction (this is the shooting method). In the animation, ħ = 2m = 1. Restart. For what values of the energy (and therefore k) is the boundary condition at x = L satisfied?
At the boundary x = L we must have a k such that ψ(L) = Asin(kL) = 0. Satisfying this boundary condition can be accomplished by requiring that kL = nπ, where n = 0, ±1, ±2, ±3,…. We can eliminate the negative values of n as these energy eigenfunctions are different from the positive n values by just an overall phase factor of −1. The n = 0 value must be considered more carefully. If we allow n = 0, this yields k = 0 and amounts to a zerocurvature solution to the timeindependent Schrödinger equation. A zerocurvature solution is of the form Ax + B and cannot satisfy the boundary conditions and be nonzero in the well. Therefore, k = 0 is not a valid possibility^{2 }and we have that:
ψ_{n}(x) = Asin(nπx/L) for 0 < x < L with n = 1, 2, 3,… , (10.2)
and zero otherwise. Since k = (2mE/ħ^{2})^{1/2}, the energy becomes:
E_{n} = n^{2}π^{2}ħ^{2}/2mL^{2}. (10.3)
But what about A? The energy eigenfunction must be normalized to satisfy Born's probabilistic interpretation, so
∫ ψ_{n}^{*}(x)ψ_{n}(x) dx = 1 . [integral from −∞ to +∞] (10.7)
Since ψ_{n}(x) is a real function and most of the spatial integral vanishes because the energy eigenfunction is zero everywhere except between 0 and L, we are simply left with calculating
A^{2 }∫ sin^{2}(nπx/L) dx → A^{2 }∫ sin^{2}(nπx/L) dx = A^{2}L/2 = 1 , [integral from −∞ to +∞] (10.8)
which tells us that A = (2/L)^{1/2}. Therefore the energy eigenfunction
ψ_{n}(x) = (2/L)^{1/2}sin(nπx/L) for 0 < x < L (10.4)
satisfies the normalization condition. We also can consider the integral
∫ ψ_{m}^{*}(x)ψ_{n}(x) dx = Ln cos(nπ)sin(mπ)/((m^{2 }− n^{2})π) + Lmcos(mπ)sin(nπ)/((n^{2 }− m^{2})π) = 0 [integral from 0 to L]
for m ≠ n. Therefore we can represent these two equations together as
∫ψ_{m}^{*}(x)ψ_{n}(x) dx = δ_{mn} [integral from 0 to L]
where δ_{mn }is the Kronecker delta which is defined such that δ_{nn }= 1 and δ_{m ≠ n }= 0. Hence the solutions to the infinite square well are orthogonal and normalized, or orthonormal.

Please wait for the animation to completely load.
In the second animation, the first 10 normalized energy eigenfunctions are shown for a box with L = 1, along with the energy spectrum. In the animation ħ = 2m = 1. You can clickdrag in the energy spectrum on the left to change the energy state. As you do so, the displayed energy turns from green to red.
Note that these energy eigenfunctions, ψ_{n}(x), are only nonzero in the spatial region 0 < x < L and are zero everywhere else.^{3} These energy eigenfunctions are not so simple after all. In fact, these energy eigenfunctions have a kink (a discontinuous first derivative) at x = 0 and x = L. Normally this is not acceptable for a energy eigenfunction, but here the potential energy function for the infinite square well is so badly behaved, these energy eigenfunctions are actually acceptable.
We can now also calculate expectation values of position and momentum. We find that: <x> = L/2 and <p> = 0, as expected. We also find that: <x^{2}> = L^{2} (1/3 − 1/(2n^{2}π^{2})), <p^{2}> = n^{2}π^{2}ħ^{2}/L^{2}, and hence <H> = n^{2}π^{2}ħ^{2}/(2mL^{2}), again as expected.
^{2}Such an argument is also made by M. A. Morrison in Understanding Quantum Physics: A User's Manual. The correct derivation of why the k = 0 case cannot exist in the infinite square well was originally stated by M. Bowen and J. Coster, "Infinite Square Well: A Common Mistake," Am. J. Phys. 49, 8081 (1980) with followup discussion in R. C. Sapp, "Ground State of the Particle in a Box," Am. J. Phys. 50, 11521153 (1982) and L. Yinji and H. Xianhuai, "A Particle Ground State in the Infinite Square Well," Am. J. Phys. 54, 738 (1986).
^{3}These energy eigenfunctions can also be written as ψ_{n}(x) = (2/L)^{1/2} sin(nπx/L) Θ(x) Θ(L − x). This representation uses two Heaviside step functions, Θ(ξ), to explicitly show the region in which the energy eigenfunction is valid. This can be accomplished because the step function, Θ(ξ), is zero for ξ < 0 and is 1 for ξ > 0.
« previous
next »