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

A Solving the Harper equation

  1. Apr 18, 2017 #1
    For a single particle in a 2D square lattice in the presence of an Abelian magnetic field Schroedinger's equation transforms into Harper's equation

    [tex]g(m+1) + g(m-1) = [E - 2 cos(2\pi m \alpha- \nu)]g(m)[/tex]
    [tex]\psi(x,y)=\psi(ma,na) = e^{i\nu n} g(m) \\ \alpha= \frac{e a^2 B}{h c}[/tex]

    I am familiar with a solution that involves matrix multiplication and the condition about the trace of wilson loop.
    I can also plot Hofstatder butterfly by constructing hamiltonian matrix in some gauge (e.g Landau) and diagonalizing it.

    What I don't know is how to solve Harper equation in a different way or how to get butterfly by simply finding points [itex](E, \alpha) [/itex]. Solving this kind of equations is new to me. I don't know what to do if I put some particular values of E and [itex]\alpha[/itex] into equation. Do I need to assume how g(0) looks like? Any help?
  2. jcsd
  3. Apr 23, 2017 #2
    Thanks for the thread! 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? The more details the better.
  4. Jun 6, 2017 #3
    Are you asking how to numerically solve this? Here is a nice paper going over matrix expansion of the Harper equation:
    But to numerically solve you need to write the equation in matrix form and loop over possible parameter values. You should just take the raw Hofstandter Hamiltonian, expand it in matrix form, and diagonalize it numerically (MATLAB and Python are good tools for this).
    Code (Text):
    for k = -pi:pi
           H = %some cell with k plugged in
           E(k,:) = eig(H)
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted