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
Click For 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?
 

Thread 'For those who ask: "What programming language should I learn?"'

Learn If you want to write code for Python Machine learning, AI Statistics/data analysis Scientific research Web application servers Some microcontrollers JavaScript/Node JS/TypeScript Web sites Web application servers C# Games (Unity) Consumer applications (Windows) Business applications C++ Games (Unreal Engine) Operating systems, device drivers Microcontrollers/embedded systems Consumer applications (Linux) Some more tips: Do not learn C++ (or any other dialect of C) as a...
View full post »

Similar threads

Help with Monte Carlo Wang-Landau JDoS
  • · Replies 1 ·
Replies
1
Views
2K
Problems with Landau Wang algorithm in fortran
  • · Replies 8 ·
Replies
8
Views
2K
A Monte Carlo Wang-Landau_antiferromagnetic perovskite
  • · Replies 3 ·
Replies
3
Views
3K
Phonon density of states and density of states of free electrons
  • · Replies 3 ·
Replies
3
Views
2K
Why is the DOS from Wang-Landau's algoritm non-flat?
  • · Replies 2 ·
Replies
2
Views
2K
Solving the Density of States: Understanding dn/dE
  • · Replies 7 ·
Replies
7
Views
2K
I Density of states in the ideal gas
  • · Replies 2 ·
Replies
2
Views
2K
Derivation of the density of states?
  • · Replies 6 ·
Replies
6
Views
2K
What is the correct formula for the density of states in 2D for normal atoms?
  • · Replies 3 ·
Replies
3
Views
2K
A Interpreting the energy density of QFT vacuum states
  • · Replies 3 ·
Replies
3
Views
1K