Metropolis sampling of quantum hall wavefunction

I thought I’d whip up a fast little script to show the radial density of the fluid from a laughlin style wavefunction.

Doing it algebraically and exactly I’ve run into some hiccups that may or may no iron out.

Using monte carlo sampling is easy enough.

Using the metropolis algorithm i suggest a new configuration of positions of the electrons, then accept or reject depending on whether a random acceptance ratio is greater that the ratio of the probabilities of the configurations (probabilities being given by the square of the wavefunction).

Then I make a histogram.

There are a ton of inefficiencies in this implementation, some that would be very very easy to fix. For example, I could put the binning in the loop. That’s save me that huge data structure and let the thing run for hours. But whatever. Works good enough for now.

laughlindensity

Not totally sure everything is right, but at least we are getting a filling of 1/3 for the laughlin state. So that’s something.

I believe that peak is an actual effect. I think I’ve seen it in plots before.

Leave a Reply

Your email address will not be published. Required fields are marked *