Help with Monte Carlo Wang-Landau JDoS

  • Thread starter Thread starter UFSJ
  • Start date Start date
  • Tags Tags
    Monte carlo
Click For Summary
SUMMARY

The discussion focuses on the Wang-Landau algorithm applied to a magnetic perovskite with exchange interactions J1 = 1.66 and J2 = -1.16. It highlights that the joint density of states obtained through the Wang-Landau method does not guarantee the identification of all possible microstates due to the flatness test convergence criterion. The example of a 20x20 Ising model illustrates the impracticality of sampling all 2^400 microstates, emphasizing the limitations of the algorithm in capturing the complete state space within a reasonable timeframe.

PREREQUISITES
  • Understanding of the Wang-Landau sampling algorithm
  • Familiarity with Monte Carlo methods
  • Knowledge of magnetic systems and exchange interactions
  • Basic concepts of statistical mechanics
NEXT STEPS
  • Explore advanced techniques for improving convergence in Wang-Landau sampling
  • Investigate alternative algorithms for sampling microstates in complex systems
  • Learn about the implementation of the Ising model in computational physics
  • Study the implications of microstate sampling limitations on physical properties
USEFUL FOR

Researchers and practitioners in computational physics, particularly those working with Monte Carlo simulations, statistical mechanics, and magnetic materials.

UFSJ
Messages
13
Reaction score
2
Hi, guys.

I have tried to write a Wang-Landau JDoS algorithm to describe a magnetic perovskite with exchange interactions J1 = 1.66 and J2 = -1.16. Then, I have a simple question: in the WL algorithm, the obtained joint density of states must have all possible E x M microstates? Since the convergence criterion in WL is just the flatness test after some Monte Carlo steps (e.g., n * 10^6), it is not guaranteed that all microstates will be identified, correct???
 
Last edited by a moderator:
Technology news on Phys.org
It is not guaranteed that you constructed all microstates. So if that's what you meant by "[not] all microstates will be identified", then you are correct.

Simple example: A 20x20 Ising model has 2^400= 2.5*10^120 microstates. Generating 10^10 microstates per second (10 per nanosecond) would mean that you would need about 10^100 years to sample all microstates. A 2D WL run for such a model will probably take a couple of seconds, maybe minutes.
 
Reply
  • Like
  • Informative
Likes Tom.G and pbuk

Thread 'Who is responsible for the software when AI takes over programming?'

I tried a web search "the loss of programming ", and found an article saying that all aspects of writing, developing, and testing software programs will one day all be handled through artificial intelligence. One must wonder then, who is responsible. WHO is responsible for any problems, bugs, deficiencies, or whatever malfunctions which the programs make their users endure? Things may work wrong however the "wrong" happens. AI needs to fix the problems for the users. Any way to...
View full post »

Similar threads

Graduate Monte Carlo Wang-Landau_antiferromagnetic perovskite
  • · Replies 3 ·
Replies
3
Views
3K
Problems with Landau Wang algorithm in fortran
  • · Replies 8 ·
Replies
8
Views
2K
Why is the DOS from Wang-Landau's algoritm non-flat?
  • · Replies 2 ·
Replies
2
Views
2K
Graduate MOND Theory Primer: Stacy McGaugh's Research Program Explored
  • · Replies 11 ·
Replies
11
Views
4K