schroedinger_gaussian
Graham Dennis
Evolve a superposition of the groundstate and first excited state.
t
V1
groundstate_gaussian_break.h5
phi1
excitedstate_gaussian_break.h5
phi
groundstate_gaussian_break.h5
psi
excitedstate groundstate
psi = 1.0/sqrt(2.0) * (phi1 + phi);
500 500
T
T = -0.5*i*ky*ky;
wavefunction
potential
dpsi_dt = T[psi] + V1*psi;
dens
wavefunction
dens = mod2(psi);
T V
wavefunction potential
L
L = i*ky;
// Calculate the potential energy at each point V = (-V1.Im()) * mod2(psi); // Calculate the kinetic energy at each point T = 0.5*mod2(L[psi]);