A Python code to simulate Schrodinger’s Time Dependent Wave Equation Simulation for Particle in a Box (1D).
Schrodinger wave equation is a mathematical expression describing the energy and position of the electron in space and time, taking into account the matter wave nature of the electron inside an atom. It is based on three considerations which are listed below:
- Classical plane wave equation
- Broglie’s Hypothesis of matter-wave
- Conservation of Energy
The initial ψ(x, 0) is generated by a complex gaussian wave packet.
The first step is to normalize the wave function given by the equation below,
Now, the energy Eigen-state wave functions are then found. Again, the energy Eigen-state wave function need to be normalized and the following equation is derived,
Discrete quantized wave vectors are required to solve the Time-Dependent equation, ψ(x, t). The energy quantization for the system becomes,
Finally, the last item the time dependent Schrödinger equation depends on is Cn, the expansion coefficient that changes the probability amplitude. The indefinite integral of the complex conjugate of eigenstate wave function and the initial wave function with respect to the spatial dimension x,
Thus the final result that we can obtain by combining the expressions for ψn(x), En and Cn is,
Now, as t increases the wave equation is determined.
Following are the conclusions that we can derive after gaining the simulation and mathematical solution of the wave equation,
- The mathematical result for the solution of the wave equation turned out to be a complex value consisting of a real and an imaginary part.
- We plotted the absolute values, real values and imaginary values on the plot plane and used pause() function available in matplotlib Python Library to obtain a simulation out of continuous plotting of the results.