I am simulating electrons inside a cylindrical well like the one shown on the first figure.(adsbygoogle = window.adsbygoogle || []).push({});

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.

**Physics Forums - The Fusion of Science and Community**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Numerical solution of SE for cylindrical well

Loading...

Similar Threads - Numerical solution cylindrical | Date |
---|---|

I Numerical solution of Schrödinger equation | Dec 9, 2016 |

Numerical solution of Schrodinger equation | Sep 5, 2015 |

Schrodinger equation numerical solution | Aug 23, 2015 |

Numerically solutions with periodic boundary conditions | Sep 22, 2014 |

Numerical solution of one dimensional Schrodinger equation | Apr 18, 2014 |

**Physics Forums - The Fusion of Science and Community**