A simple Python program for numerically solving the Schrödinger equation in one dimension.
This was originally developed as part of my Bachelor's thesis. For the original code used for the thesis, see branch bachelor's-thesis.
When the program is run it will open an animated plot that visualizes the evolvement of the chosen initial wave function under the influence of a chosen potential according to the Schrödinger equation.
The currently available algorithms for solving the Schrödinger equation are the finite difference methods:
To run the program, in the src/ directory simply run simulate.py, e.g.
./simulate.pyThe default setup is an arbitrary harmonic oscillator with a gaussian wave packet used as the wave function.
- More dimensions
- More solving algorithms
- More potentials