Interesting Examples of Computational Physics
  • Welcome to the blog of Computational Physics at the IFM!
  • Dynamical Systems
    • The simple pendulum
    • The spherical pendulum
  • The wave equation
    • The simple discretization
    • Boundary conditions
    • Interfaces between mediums
    • Knots on ropes
    • 2D wave equation
  • Time dependent Schrödinger equation
    • The Crank-Nicolson method
    • Free wave packet
    • Quantum Tunneling
    • Harmonic Oscillator
    • Forced Harmonic Oscillator
    • 2D Harmonic Oscillator
    • Wave-particle duality
Powered by GitBook
On this page

Time dependent Schrödinger equation

Previous2D wave equationNextThe Crank-Nicolson method

Last updated 9 months ago

In this section, illustrative examples are shown of how the wave function evolves under certain scenarios. We present animations of the following scenarios:

  • The

  • The , and harmonic oscillators

  • Examples of the

We start by presenting the equation of movement that rules non-relativistic quantum systems, which is Schrödinger's time dependent equation:

iℏ∂Ψ(x,t)∂t=−ℏ22m∇x2Ψ(x,t)+U(x,t)Ψ(x,t){\rm i} \hbar\frac{\partial \Psi(x,t)}{\partial t} = -\frac{\hbar^2}{2m}\nabla_{x}^2 \Psi(x,t) + U(x,t)\Psi(x,t)iℏ∂t∂Ψ(x,t)​=−2mℏ2​∇x2​Ψ(x,t)+U(x,t)Ψ(x,t)

which, unlike the wave function, describes the evolutiono of complex instead of a real function. However, one computational problem is that the scale at which quantum phenomena are appreciable is really small. For this reason, a change of variables is suggested. One change of dimensions that I came up with is the following:

ξ=xa0,τ=E0ℏt,E0=ℏ2ma02\xi = \frac{x}{a_0},\tau = \frac{E_0}{\hbar}t,E_0 = \frac{\hbar^2}{ma_0^2}ξ=a0​x​,τ=ℏE0​​t,E0​=ma02​ℏ2​

where a_0 is a scale of length appropieate for the problem to be treated. For instance, if we are working with atomic scales, it could be one angstrom or Bohr's radius. These new variables are dimensionless and will be more apropiate for numerical calculations. Along with the axis escalation, the potential and wave functions must also change:

V(ξ,τ)=U(x,t)E0,Ψ(ξ,τ)=a0nψ(x,t),V(\xi,\tau) = \frac{U(x,t)}{E_0},\Psi(\xi,\tau) = \sqrt{a^n_0}\psi(x,t),V(ξ,τ)=E0​U(x,t)​,Ψ(ξ,τ)=a0n​​ψ(x,t),

were n is the dimension in which the wave function is being solved. This is to ensure a dimensionless potential and a normalized wave function. Finally, whenever we want to calculate the momentum of the particle, we can use the dimensionless operator:

pε^=a0ℏp^=−i∇ε\hat{p_{\varepsilon}} = \frac{a_0}{\hbar}\hat{p} = -i\nabla_{\varepsilon}pε​^​=ℏa0​​p^​=−i∇ε​

The resulting dimensionless equation is the following:

i∂Ψ(ξ,τ)∂τ=−12∇ξ2Ψ(ξ,τ)+V(ξ,τ)Ψ(ξ,τ){\rm i} \frac{\partial \Psi(\xi,\tau)}{\partial \tau} = -\frac{1}{2}\nabla_{\xi}^2 \Psi(\xi,\tau) + V(\xi,\tau)\Psi(\xi,\tau)i∂τ∂Ψ(ξ,τ)​=−21​∇ξ2​Ψ(ξ,τ)+V(ξ,τ)Ψ(ξ,τ)

However, to save some writing and have a more comfortable notation, we are going to simply use x and t again, remembering that that we are using code units:

i∂ψ(x,t)∂t=−12∇2ψ(x,t)+V(x,t)ψ(x,t){\rm i}\frac{\partial \psi(x,t)}{\partial t} = -\frac{1}{2}\nabla^2 \psi(x,t) + V(x,t)\psi(x,t)i∂t∂ψ(x,t)​=−21​∇2ψ(x,t)+V(x,t)ψ(x,t)

In the following chapters, we are going to solve this equation with different potentials and initial conditions.

free wave packet
Quantum tunneling
simple
forced
2D
wave-particle duality