TwoDMSse
Joe Hope
Schroedinger equation for multiple internal states in two spatial dimensions.
t
phi
phi = j*sqrt(2/sqrt(M_PI/2))*exp(-(x*x+y*y)/4)*exp(i*0.1*x);
U
U=exp(-(x*x+y*y)/4);
V
V=3*(j*(1-k)+(1-j)*k);
VPhi
internalInteraction wavefunction
// Calculate the current normalisation of the wave function. VPhi = V*phi(j => k);
20 100
wavefunction
Lx
Lx = -i*kx*kx*0.5;
Ly
Ly = -i*ky*ky*0.5;
dphi_dt = Lx[phi] + Ly[phi] -i*U*VPhi; spatialInteraction coupling
density
wavefunction
density = mod2(phi);
normalisation
wavefunction
normalisation = mod2(phi);