Skip to content

๐Ÿซ Hamiltonian eigenstates

JavaScriptThree.js

This demo visualizes the stationary eigenstates of a 2D wave function:

iโ„โˆ‚ฯˆ(x,y,t)โˆ‚t=[โˆ’โ„22m(โˆ‚2โˆ‚x2+โˆ‚2โˆ‚y2)+V(x,y)]ฯˆ(x,y,t)i\hbar \frac{\partial \psi (x,y,t)}{\partial t}=\left[ -\frac{\hbar ^{2}}{2m}\left(\frac{\partial ^{2}}{\partial x^{2}}+\frac{\partial ^{2}}{\partial y^{2}}\right)+V(x,y) \right]\psi (x,y,t)

Where:

  • ฯˆ(x,y,t)\psi(x, y, t) is the wave function of the system.
  • โˆ‚2/โˆ‚x2+โˆ‚2/โˆ‚y2=โˆ‡2\partial^2/\partial x^2 + \partial^2/\partial y^2 = \nabla^2 is the 2D Laplacian operator.
  • โ„\hbar is the reduced Planckโ€™s constant, set to 11 in the code.
  • mm is the mass of the particle, set to 11 in the code.
  • V(x,y)V(x, y) is the potential energy function, e.g. an isotropic harmonic oscillator potential:
V(x,y)=12k(x2+y2)V(x,y)=\dfrac{1}{2}k(x^2+y^2)
Source

Select different eigenstates to explore how their amplitude and phase form distinct wave-like patterns, including the characteristic lobes and nodes of quantum wavefunctions.

For the isotropic harmonic oscillator, the energies are given by k(nx+ny+1)\sqrt{k}(n_x +n_y+1):

EEnxn_xnyn_y
0.223600
0.447210
0.447201
0.670802
0.670811
0.670820
0.894403
0.894412
0.894421
0.894430
1.118004
1.118013
1.118022
1.118031
1.118040