Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Gaussian Wavepacket Tunnel Through double barrier

  1. Oct 3, 2009 #1
    Hello,
    I was trying to design a movie in Matlab. A Gaussian wave packet moves toward double barrier then some wave reflect and some pass out from the barrier . i design the double barrier in Matlab. but dont know how to make movie in Matlab with Gaussian wavepacket. can you help. Do u have the code in Matlab
     
  2. jcsd
  3. Oct 3, 2009 #2

    Pythagorean

    User Avatar
    Gold Member

    There is a programming subforum down below.

    Short answer: the plot command in a for loop. You plot your function (the gaussian wave packet) and then compute the new shape of the function every time through the for loop, and then plot it.

    You can use the pause command (type help pause) to put time between your frame changes.
     
  4. Oct 3, 2009 #3

    Pythagorean

    User Avatar
    Gold Member

    also check out 'help movie'. You can collect each of your plots as a movie frame and build a movie out of it, then once you have the movie array, you can adjust the frame-rate and the number of plays as you want.
     
  5. Oct 3, 2009 #4
    Hej Pythoagorean,
    Thankyou very much for reply. i use the movie frame and it work for some equation. but i was trying move to gaussian wave packet but dont know how implement on matlab mathematically . do u have code for Gaussian wave packet
     
  6. Oct 3, 2009 #5

    Pythagorean

    User Avatar
    Gold Member

    Here's the general analytic structure:

    http://musr.physics.ubc.ca/~jess/hr/skept/GWP/

    Shooting from the hip, I'd use a runge-kutta algorithm to handle the integral, but you may want to see if there's some more experienced replies. I've never coded a gaussian wave packet before.
     
    Last edited by a moderator: Apr 24, 2017
  7. Oct 3, 2009 #6

    f95toli

    User Avatar
    Science Advisor
    Gold Member

    A split-operator FFT metthod is usually more efficient than solving the integral directly for these kinds of problems.
    You should be able to find a description of this method in just about any book on computational physics.
     
  8. Oct 4, 2009 #7
    If I'm understanding these suggestions correctly, then what is being advised is to model the Gaussian wave packet as the superposition of its component frequencies. While this is certainly correct, it's probably overkill for the animation in question. I'm pretty sure you'll get almost the identical picture simply by calculating the reflection/transmission of a monochromatic beam and then just putting a Gaussian envelope over it.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Gaussian Wavepacket Tunnel Through double barrier
  1. Gaussian Wavepackets (Replies: 4)

  2. Gaussian wavepacket (Replies: 1)

  3. Barrier Tunneling (Replies: 4)

Loading...