anderkve / FYS3150

https://anderkve.github.io/FYS3150
26 stars 14 forks source link

Task 9 plots #101

Closed Norahelgeland closed 1 year ago

Norahelgeland commented 1 year ago

Hello,

We are a bit confused regarding our plots from task 9.

For one slit we get the following plot:

image

For two slits we get the following plot:

image

We expected to observe three peaks for one slit and more peaks for two slits, and we are wondering why it is not symmetric. Is it just because the distance from the slit is too great? Also our heat map for two slits looks like it has three peaks, but this does not show up in the plot for x = 0.8.

Here is our heat map plot for t = 0.002:

image

anderkve commented 1 year ago

Hi @Norahelgeland!

Indeed, these results don't look correct. There may be multiple sources of error, but the most striking problem is how jagged the lines are, corresponding to the strange block-like pattern in the bottom heat map.

My guess would be that there is something wrong with your A and B matrices, or with the code that translates from indices (i,j) to index k. So that is where I'd start looking. A common mistake is to have a code for constructing A and B that works for the (M-2)=3 example matrices in the project description, but doesn't work for general values of M-2. If you suspect this may be the problem, it's worth working/thinking through what these matrices would look like for (M-2)=2 and (M-2)=4, to see the correct pattern.

To test the code I would focus first on the case with zero potential, i.e. just a freely propagating wave packet, and make sure that this looks sensible before you start adding slit potentials.

Norahelgeland commented 1 year ago

Thank you! It turned out that there was a mistake in the sub matrices of our A and be matrices. Now the plots look more correct :)