search

quantum-visualizations / qmsolve

⚛️ A module for solving and visualizing the Schrödinger equation.
BSD 3-Clause "New" or "Revised" License
931 stars 119 forks source link

In 1D is time double faster than p_x0 in initial gaussians? #20

Closed jesusmontera closed 2 years ago

jesusmontera commented 2 years ago

Hello, when running the sample 1d potential barrier with 0 pot value, if I give 50 AU v0 speed(same that p_x0 because mass is 1), it takes 0.55 AU to walk 50 AU, when it have to take 1 AU. Can this be fixed replacing?: np.exp(p_x0*particle.x*1j) with np.exp(p_x0*particle.x*0.5j) Take note that I've modified your code for plotting AU instead of femtoseconds: time_ax.set_text(u"t = {} au".format("%.3f" % ((animation_data['t'] * TAUFMS)))) and the same for space Amstrong plot

jesusmontera commented 2 years ago

Oh, my mistake in x range plot , 1j is ok, so time have a little near 10% error, sorry