How to normalize the density of states in JDoS Wang-Landau?

  • Thread starter Thread starter UFSJ
  • Start date Start date
  • Tags Tags
    Density of states
AI Thread Summary
The discussion centers on the normalization of the Joint Density of States (JDoS) in the Wang-Landau algorithm, particularly in the context of a bi-dimensional Ising model. The user seeks guidance on updating the logarithm of the density of states, ln[Ω(E,m)], noting that for ground states, the density must equal Q (with Q being 2 for the Ising model). They mention that for ground states, there is only one state for maximum and minimum magnetization, which could complicate the normalization process. A key point raised is that normalization may not be necessary for updating the density of states, as Monte Carlo steps primarily consider relative densities, allowing for non-normalized generation of the density of states during the algorithm's execution. Normalization can be performed after the algorithm concludes. The user also seeks clarification on what JDoS entails.
UFSJ
Messages
13
Reaction score
2
Hi guys.

I want some help understanding how I can make the normalization of the JDoS density of states (Ω[E,m]) in the Wang-Landau algorithm. When I am working with DoS (Ω[E]) I use the knowledge that the value of the density of states in the ground states must be
equal to Q (Q = 2 for the Ising model), that is, I make the update ln[Ω(E)] = ln[Ω(E)] - ln[Ω(Eground state )] + ln[2]. However, I don't know how to update the ln[Ω(E,m)] in the JDoS algorithm. I am working with a bi-dimensional space of energy (E) and magnetization (m) in an Ising spin-lattice.
 
Technology news on Phys.org
I assume that you talk about an Ising like model. If you sample in E and m, then the number of states with E_groundstate is one for m=+N and m=-N (where N is the number of spins total and m assumed to be in units of spins pointing up). The number of states with E_groundstate and any other magnetization is zero, which may actually become a problem for you if you do not handle this situation.

I think normalization is irrelevant for updating your estimate for the density of states (DOS). Your Monte-Carlo steps (assuming you do MC) only take the relative DOS into account, and the normalization drops out of the calculation. So you can just generate your DOS estimation non-normalized and normalize after the algorithm is terminated.

What is JDoS?
 
Dear Peeps I have posted a few questions about programing on this sectio of the PF forum. I want to ask you veterans how you folks learn program in assembly and about computer architecture for the x86 family. In addition to finish learning C, I am also reading the book From bits to Gates to C and Beyond. In the book, it uses the mini LC3 assembly language. I also have books on assembly programming and computer architecture. The few famous ones i have are Computer Organization and...
I have a quick questions. I am going through a book on C programming on my own. Afterwards, I plan to go through something call data structures and algorithms on my own also in C. I also need to learn C++, Matlab and for personal interest Haskell. For the two topic of data structures and algorithms, I understand there are standard ones across all programming languages. After learning it through C, what would be the biggest issue when trying to implement the same data...
Back
Top