I Tunneling with Gaussian Wave Packet

Tertius
Messages
57
Reaction score
10
TL;DR Summary
I have a simulation built of an initial wave packet approaching a barrier and tunneling.
The goal is to have accurate 1D numerical results for tunneling probabilities through an arbitrary barrier without relying on analytic approximations such as WKB. If there is a more ideal approach to this, I am happy to change tactics. Time independent, for example, but I am not sure how to set up the boundaries at the edges so it doesn't become a bound system.

Most of the resources i've found detail way to solve for bound systems (oscillators, potential wells, etc), but I haven't found one that produces tunneling probabilities from a numerical method.

My initial simulation is a Gaussian wavepacket approaching a barrier, but I am finding that the numerical results are of course dependent on the initial location of the wavepacket. This is expected because the time evolution spreads out the wavepacket as it approaches the barrier.

I am attaching a snapshot of the simulation (both real and imaginary parts of ##\psi## are shown). I am computing the tunneling probability as the integrated probability after the barrier divided by the integrated probability before the barrier (computed at each time step, and taking the maximum).

I am concerned this approach is dependent on initial position. Is there a better general approach?

previously consulted resources:
https://arxiv.org/html/2403.13857v1#S3 https://www.reed.edu/physics/faculty/wheeler/documents/Quantum Mechanics/Miscellaneous Essays/Gaussian Wavepackets.pdf

Screenshot 2024-06-12 at 07.36.41.png
 
Physics news on Phys.org
Tertius said:
My initial simulation is a Gaussian wavepacket approaching a barrier, but I am finding that the numerical results are of course dependent on the initial location of the wavepacket. This is expected because the time evolution spreads out the wavepacket as it approaches the barrier.
You should try a higher central momentum for the wave packet then. You want to be in a regime where the spreading of the wave packet is negligible (compared to ##d\braket{x}/dt##) over the entire simulation.

Tertius said:
I am attaching a snapshot of the simulation (both real and imaginary parts of ##\psi## are shown). I am computing the tunneling probability as the integrated probability after the barrier divided by the integrated probability before the barrier (computed at each time step, and taking the maximum).
In the conditions I mentioned above, you will get a splitting of the wave packet into separated reflected and transmitted parts, from which you can calculate the transmission probability, which will be independent of time once the interaction with the barrier is over (assuming that the wave packet never reaches the end of the grid). The maximum might not correspond to the long-time (converged) result.

Also, for square barriers you can get exact transmission probabilities for a plane wave (see, e.g., Sakurai's Modern Quantum Mechanics). The numerical result can of course differ as you are using a wave packet, but it is possible to make the wave packet narrow enough in momentum such that the two results are practically indistinguishable.
 
That makes sense. And then, from what i've found in literature, it is common to have the spreading parameter ##\sigma = \frac{h}{2 \lambda}##. That should, in combination with the momentum being high enough, give physically real and consistent results, I would hope.
 
Not an expert in QM. AFAIK, Schrödinger's equation is quite different from the classical wave equation. The former is an equation for the dynamics of the state of a (quantum?) system, the latter is an equation for the dynamics of a (classical) degree of freedom. As a matter of fact, Schrödinger's equation is first order in time derivatives, while the classical wave equation is second order. But, AFAIK, Schrödinger's equation is a wave equation; only its interpretation makes it non-classical...
I am not sure if this falls under classical physics or quantum physics or somewhere else (so feel free to put it in the right section), but is there any micro state of the universe one can think of which if evolved under the current laws of nature, inevitably results in outcomes such as a table levitating? That example is just a random one I decided to choose but I'm really asking about any event that would seem like a "miracle" to the ordinary person (i.e. any event that doesn't seem to...
Back
Top