Hover over the menu bar to pick a physics animation.

uantum Particle Wave Packet in a Parabolic Potential

In this version I include two features that were not included previously.
1. For a gaussian wave packet envelope, I plot the algebraic expression that agrees with the Time Dependent Schrodinger Equation (TDSE) solution of the width variation as a function of time.
2. The modeling of a three initial envelope shape wave packets in terms of the stationary states of the parabola.
The motivation for including these items was a book "Quantum Mechanics" 3rd Edition, 1998 John Wiley and Sons Inc. by Eugen Merzbacher The animation shows the propagation of a quantum wave packet

`psi_0(x,0)="envelope"(x)exp(ikx)`

in a parabolic potential well. The wave packet that I use is a simple complex sinusoid embedded in various envelopes (see "Starting Wave Packet Envelope Shapes"). Its propagating function changes are most accurately computed by solving the time dependent Schrodinger equation (TDSE)

`iℏ(dPsi)/dt=-ℏ^2/(2m)(del^2Psi)/(delx^2)+V(x)Psi`

using finite element analysis (FEA) but this animation just uses finite difference (FD) methods to get adequate accuracy. The time dependence of a given starting packet can also be obtained by use of a series of the stationary states with their respective time factor

`psi(x,t)=sum_(n=0)^ooc_n*H_n(x)exp(-i(n+1/2)*omega*t) `

where `omega` is defined below, `H_n(x)` are the Hermite polynomials, and `c_n` are the coefficients needed to fit the initial wave packet. For reasonably small initial wave vectors, an upper limit of n=100 or so is sufficient. For this analysis we will let the Planck constant, ℏ, and the mass, m, of the particle both equal 1. FEA clearly shows that the wavepacket's propagation is strictly periodic with a period defined by the coefficient of the quadratic potential `V(x)=(b/2)x^2` so that the period is `2*pi/omega` where `omega=sqrt(b)`. Note that `V(x)` is the potential energy associated with a mass and linear spring system where the spring constant is `b`. The behavior of the wave packet is another manifestation of the Bohr Correspondence Principle Just as a classsic particle moves in a parabolic potential well, the centroid of the magnitude squared of the wave packet moves sinusoidally with function `x_c(t)=x_(Max)*sin(omega*t)` where `x_(Max)=k_0/omega` is the turn around value of x and where `k_0` is the initial wave vector, `k`, of the packet. The Schrodinger equation requires that the integral of the wave packet's magnitude be constant but its amplitude and width can vary and these will be shown by the animation. A very important distinction between what may be expected from classic physics and what happens in quantum physics is the effect of the initial width, `sigma_0`, of the wave packet envelope. To show how the time dependent solutions relate to the stationary state solutions, I have provided a check box "Show Stationary State Plots" which plots the steady state solution, eigenfunction, of the parabolic well for the particle energy, `En=1/2 k_0^2` which is the same as `(n+1/2)ℏomega`. Note that the spatial frequencies of the moving wave packet and the eigenfunction are reasonably equal over the entire range, `|x_(Max)|`. I let `sigma_0=r*x_(Max)` where r is less than 1 and is selectable by a slider. I also provide sliders to select Initial Wave Velocity, Parabolic Potential Coefficient, Time Increment for solving the TDSE, and the total number of grid points (Gridding Range). In cases where the quality of the plots seems less than desirable, the learner might want to experiment with these latter two sliders. After a full period of propagation, I provide an option to plot the position and momentum variances for the present set of parameters. As a check on the solution accuracy, included in these plots is that of the integral of the magnitude of `psi`. The equation that provides the width variation of a gaussian packet is given as

`deltax^2=(cdx_0)^2+(sdp_0)^2`

where `c=cos(omega*t)`, `s=sin(omega*t)`, `dx_0^2` is the position variance of the initial packet and `dp_0^2=1/(bdx0^2)` is the momentum variance of the initial packet. When the "gaussian" starting shape is chosen, this expression is plotted as a black trace along with the iterative or series time dependent solutions. Using a Chrome browser, the reader is welcome to press F12 on a Windows computer in order to examine the Javascript code that I have written.