File:Qm timedomain debroglie wavelength.gif

From Wikimedia Commons, the free media repository
Jump to navigation Jump to search

Qm_timedomain_debroglie_wavelength.gif (360 × 257 pixels, file size: 8.36 MB, MIME type: image/gif, looped, 1,400 frames)

Note: Due to technical limitations, thumbnails of high resolution GIF images such as this one will not be animated. The limit on Wikimedia Commons is width × height × number of frames ≤ 100 million.

Captions

Captions

Gaussian wavepacket accelerating

Summary

[edit]
Description
English: Gaussian wavepacket starting at rest and uniformly accelerating. A particle at rest has no wavelength associated with it, but if its velocity (momentum) starts increasing, its de Broglie wavelength will become shorter and shorter.
Date
Source https://twitter.com/j_bertolotti/status/1108320623120601093
Author Jacopo Bertolotti
Permission
(Reusing this file)
https://twitter.com/j_bertolotti/status/1030470604418428929

Mathematica 11.0 code

[edit]
L = 100; (*System size*)
dx = 0.1; (*Space step*)
dt = 0.01; (*time step*)
k0 = 0; (*Momentum of the incident packet*)
V1 = Table[30 - (30 x)/L, {x, 0, L, dx}]; (*Potential*)
dim = Dimensions[V1][[1]]
(*For a undergraduate level explanation of why this works, see physik.uni-graz.at/~pep/Theses/BachelorThesis_Wachter_2017.pdf *)

\[Alpha] = (I dt)/(2 dx^2);
\[Xi]1 = Table[1 + (I dt)/2 (2/dx^2 + V1[[j]]), {j, 1, dim}];
\[Gamma]1 = Table[1 - (I dt)/2 (2/dx^2 + V1[[j]]), {j, 1, dim}];
U1 = SparseArray[{Band[{1, 1}] -> \[Xi]1, Band[{2, 1}] -> -\[Alpha], Band[{1, 2}] -> -\[Alpha]}, {dim, dim}];
U2 = SparseArray[{Band[{1, 1}] -> \[Gamma]1, Band[{2, 1}] -> \[Alpha], Band[{1, 2}] -> \[Alpha]}, {dim, dim}];
\[Psi]0 = Table[E^(-((x - L/10)^2/(2 2^2))) E^(I k0 x), {x, 0, L, dx}]; (*Initial condition*)
nsteps = 1400; (*Number of time steps*)
U = Inverse[U1].U2; (*Evolution operator*)
evo1 = NestList[Dot[U, #] &, \[Psi]0, nsteps]; (*Apply the evolution operator repeatedly to obtain the soltion at each time step*)
p1 = Table[
    ListPlot[{Re[evo1[[j]] ], Abs[evo1[[j]]]}, Joined -> True, PlotRange -> {-1.1, 1.1}, Axes -> False, PlotStyle -> {Directive[Red, Thick], Directive[Black, Thick]}, PlotLegends -> Placed[{"Re[\[Psi]]", "|\[Psi]|"}, Top] ]
   , {j, 1, nsteps, 1}];
ListAnimate[p1]

Licensing

[edit]
I, the copyright holder of this work, hereby publish it under the following license:
Creative Commons CC-Zero This file is made available under the Creative Commons CC0 1.0 Universal Public Domain Dedication.
The person who associated a work with this deed has dedicated the work to the public domain by waiving all of their rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law. You can copy, modify, distribute and perform the work, even for commercial purposes, all without asking permission.

This file, which was originally posted to https://twitter.com/j_bertolotti/status/1108320623120601093, was reviewed on 21 March 2019 by reviewer Ronhjones, who confirmed that it was available there under the stated license on that date.

File history

Click on a date/time to view the file as it appeared at that time.

Date/TimeThumbnailDimensionsUserComment
current13:49, 21 March 2019Thumbnail for version as of 13:49, 21 March 2019360 × 257 (8.36 MB)Berto (talk | contribs)User created page with UploadWizard

There are no pages that use this file.

Metadata