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

Numerical solution of SE for cylindrical well

  1. Sep 20, 2015 #1
    I am simulating electrons inside a cylindrical well like the one shown on the first figure.
    My current work has been on solving the Schrodinger equation numerically for the above potential and then finding corrections to the solution such that it is consistent with Poissons equation.
    To do so I need to apply the inverse Laplacian to the electron density as shown in figure 2, which will give me the correction to the potential profile (at least in the first iteration). The problem is however, that applying the inverse Laplacian to the electron density shown gives me a correction as shown on figure 3. As seen it makes sense that this is peaked around the origin, since the electron density is highest at this point. What however doesn't make sense, is that it curves upwards in one direction and doesn't go to zero in a spherically symmetric manner. Everything in my problem has circular symmetry, so I have no idea why it doesn't go to zero with circular symmetry. Does anyone have an idea why this can be? Is it a numerical problem? I doubt so since my mesh size is not very big.

    Attached Files:

    • 1.png
      File size:
      7.6 KB
    • 2.png
      File size:
      13.1 KB
    • 3.png
      File size:
      16.7 KB
  2. jcsd
  3. Sep 25, 2015 #2
    Thanks for the post! This is an automated courtesy bump. Sorry you aren't generating responses at the moment. Do you have any further information, come to any new conclusions or is it possible to reword the post?
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook