using a suitable Hamiltonian operator Ĥ corresponding to the system. The solution is a complex -valued wavefunction ψ (r, t), a function of the 3D position vector r of the particle at time t, describing the behavior of the system.