Motion of Electron Wave Packet in a 1D Square Potential Well of Infinite Depth

The modes are the eigenvalues and eigenvectors (psi(x))of the second order Schrodinger equation (SE) where h with the bar (ℏ) is the Planck constant, m is the mass and E is the energy of the particular mode whose wave function is the eigenvector. For the square well of infinite depth, the modes are ψ_n(x)=sin(nπx/w) where n is an integer from 1 to infinity and is the mode number, w is the width of the well and x runs between 0 and w. This leads to energies E_n=(nℏπ/w)²/(1*mass) For this program I have chosen to set the Planck constant (divided by 2*pi), ℏ, equal to its actual value, 1.056e-34 Joule seconds and the mass, m, equal to that of an electron 0.911e-30. The animation shows the propagation of a quantum wave packet in a infinitely deep square potential well. The wave packet that I use is a simple complex sinusoid embedded in various envelopes (see "Starting Wave Packet Envelope Shapes"). Here the time dependence of a given starting packet will be obtained by use of a series of the stationary states with their respective time factor ψ(x,t)=sum[cn*ψ_n(x)exp(-i*(E_n/hbar)*t), {n,0,infinity}] where En are the eigenvalues of the potential, ψ_n(x) are the stationary state solutions, and cn are the coefficients needed to fit the initial wave packet. For reasonably small initial wave vectors, an upper limit of n=20 or so is sufficient. I give the learner the choice of 3 wave packet envelope shapes:
1: Gaussian (standard bell curve used in statistics)
2: Triangular
3: Square
The evolution of the Gaussian is the slowest but it eventually gets spread out over the full potential well width. The Triangle, which has a slope discontinuity, remains intact for a classic period or two but it definitely loses its resemblance to its original shape. The Square, which has a value discontinuity loses its identity after only 1 classic period. The evolution of these shapes is in stark contrast to that seen in the case of a parabolic potential.
Since these are all composed of higher mode numbers, one should expect, by the correspondence principle, that their behavior should emulate that of a particle bouncing between two walls but it does not. For the parabolic potential, the motion of the wave packet centroid exactly matches the motion of a classic particle in the same potential and is therefor periodic at the same period. For the square well potential, the motion of the packet centroid starts out approximately matching the particle motion but with a lag. However, since the time dependence of the packet is exp(-i(En/ℏ)t) it turns out that the packet motion is indeed periodic with a period tP=2*π*ℏ/E₁ where E₁=1/2*(π/w)². To show that the packet returns to it original shape and position after this period I have given the option of either Quantum Period or Classic Period. The Quantum Period is much longer than the Classic Period so the learner will have to ' be patient to see that the Quantum Wave Function is exactly periodic. Thus, although the wave packet propagation in a square well is not periodic with the same period as an equivalent classic particle, it is periodic with the period 2*π*ℏ/E₁ where E₁ is the energy of the lowest stationery state. This period, which I call a quantum period, is usually much longer than that of a classic particle. At times distant from the start and end of a quantum period the centroid of the packet stays reasonably well centered in the square well so its location is unknown. The packet spends most of its time on the flat portion of the potential. On this portion of the domain, the Schrodinger equation reads dψ/dt=-i*ℏ/(2*mass)*d2ψ/dx2. With our choice of ψ=exp(-i*k*x), the second derivative is just d2ψ/dx2=-k²ψ. Therefore decreasing the value of ℏ/mass causes the change of ψ (or wave packet speed) to be decreased by the same factor. It turns out that this decreases both the rate of motion of the centroid of the packet as well as the rate of spreading of the width of the packet by the same factor. i.e. v_packet=v_packet1*hbar/mass and w_packet=sqrt(w_packet0^2+i*hbar/mass*time) where v_packet1 is the rate of speed of the packet when hbar/mass=1 and w_packet is the initial width of the packet which becomes a constant width when ℏ goes to zero. Recall that the velocity of a particle packet is related to its wave vector, k, as (ℏ/mass)k. For a square well of width w, the lowest value of k is pi/w. Therefore the speed of all of the states of a square well scales as ℏ/mass. The width of the packet grows as σ=√(σ0^2+i*ℏ/mass*t) where i=√-1. For the lowest mass "regular" particle, the electron, this ratio is about 0.000116 m^2/sec. Reasonable choices of square well width would be that for an electron with one electron volt (eV) of energy. This results in widths that are commensurate with the orbitals of the Hydrogen atom. The average wave vector, k, for the wave packet should result in many waves within the square well width. A choice of average energy of the wave packet should then be a few hundred eV and that is what I have provided here.