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!

# Finite difference Schrodinger equation

Loading...

Similar Threads - Finite difference Schrodinger | Date |
---|---|

A DFT Code for Finite-Temperature Calculations | May 8, 2017 |

I Solving the 3D Poisson Equation Using Finite Difference/Volume | Nov 12, 2016 |

B Finite difference problem | May 16, 2016 |

Finite Differences in Inhomogeneous Media | Apr 21, 2009 |

Nonlinear parabolic equations: finite difference method | Jan 3, 2009 |

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