# Numerical diagonalization of Hamiltonian

Framework

I am trying to diagonalize the Bogoliubov-de Gennes Hamiltonian. The problem is that the Hamiltonian contains a Laplacian. This could be solved by using a discretized Laplacian.

How I tried it

If I use a simple Hamiltonian (just the Laplacian) and use eig(l1dd(...,...)) with l1dd the discretized Laplacian (2 and -1's) as in this example from John Burkardt as FSU, then I get:

What I expect

I expect a linear or quadratic dispersion ($$\omega \propto k \text{ or } k^2$$), but I get neither of the two.

What could be the problem? It seems that eig does a good job, so probably l1dd is the problem?

• We really can't have the first idea what could be the problem. We don't know your code, and the only graph you show doesn't have axes labels and only shows some curve that looks like a sine. I'm afraid you have to provide more information, along with what you have already tried to figure out what is wrong. Feb 22 at 20:52
• I want to know the eigenvalues of the laplacian operator using discretized methods in MATLAB
– MrQ
Feb 22 at 21:30
• @WolfgangBangerth: I have given the code and the graph represents the eigenvalues
– MrQ
Feb 22 at 21:33
• @MrQ It looks like you are the OP and may have accidentally made a second account. If you want to merge these, contact the CMs at this link. You can include your whole Matlab input that you used to produce that plot in your question so that other users can try/edit it to see if they can find/fix the issue. Feb 22 at 22:33
• @MrQ It's worth getting in the habit of labeling all axes. "The graph represents the eigenvalues" is not a useful statement. It could be that each point of the curve represents the real and imaginary parts of eigenvalues in the complex plane. Or the x-axis is the index and the y-axis the absolute value of complex eigenvalues. Or so many more possibilities. Feb 23 at 5:12