Suppose I am given some 1D Hamiltonian:(adsbygoogle = window.adsbygoogle || []).push({});

H = ħ^{2}/2m d^{2}/dx^{2}+ V(x) (1)

Which I want to solve on the interval [0,L]. I think most of you are familiar with the standard approach of discretizing the interval [0,L] in N pieces and using the finite difference formulas for V and the second derivative in (1), which can then be formulated as a matrix equation, which may be diagonalized for the eigenvector and eigenvalues.

Now for all this to work one has to assume that the wave-function goes to zero outside the interval [0,L], which follows if one makes the effort of writing up the finite difference expressions. My question is: Is there a way to enforce another boundary condition with this method? I am solving a problem, where it would be beneficial to enforce the wave function to take a non-zero value on the boundaries of the interval [0,L]. Is this possible with the standard finite difference method or should I look at more advanced methods?

**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!

# I Finite difference Hamiltonian

Have something to add?

Draft saved
Draft deleted

Loading...

Similar Threads - Finite difference Hamiltonian | Date |
---|---|

I Finite universe, quantized linear momentum, and the HUP... | Sep 7, 2017 |

I Finite square well | Mar 19, 2017 |

I Finite difference method for Schrödinger equation | Oct 23, 2016 |

Some questions about the Finite Potential Well | Aug 24, 2016 |

A Approximating a QF with finite-dimensional Hilbert space | Jul 14, 2016 |

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