A Convergence of lattice Ising model

Leonardo Machado
Messages
56
Reaction score
2
Hello everyone.
I'm working on a program to solve 2D Ising model of magnetic materials, using a system with 10x10 spins for simplicity at a temperature of 1E-8 K. I'm using this parameters to get a faster result of m=1 and guarantee it is correct. but...

For now i already pass 300 Monte Carlo's moves and it is in the same position of when it was in 50... is that normal ? I've watch a vídeo on youtube where they took like a 10 million moves to compute a 200x200 system. I'm a bit scared with the time it will cost me.

Thanks.
 
Physics news on Phys.org
Could I check what your initial state was? If you started with an arbitrary random state and immediately put it at such a low temperature, it is highly likely for the system to become stuck in a local equilibrium. Most of the Monte Carlo moves will be rejected because the extremely low temperature essentially prevents moves that increase the energy of the system from occurring.

Also, it is also worth noting that such Monte Carlo simulations are often measured in terms of lattice sweeps (total number of moves / number of lattice sites) - the idea being that in each update one should allow for every lattice spin to have a chance to flip on average.
 
Fightfish said:
Could I check what your initial state was? If you started with an arbitrary random state and immediately put it at such a low temperature, it is highly likely for the system to become stuck in a local equilibrium. Most of the Monte Carlo moves will be rejected because the extremely low temperature essentially prevents moves that increase the energy of the system from occurring.

Also, it is also worth noting that such Monte Carlo simulations are often measured in terms of lattice sweeps (total number of moves / number of lattice sites) - the idea being that in each update one should allow for every lattice spin to have a chance to flip on average.

yes, I've made it totally random and in it first move it was already at such temperature.. So, how should i proceed ? Low the temperature gradually ?
 
I usually prefer the alternative - to heat the system up instead. So instead of starting from a random state, I start with the zero temperature limit of perfectly aligned spins. This avoids the problem of getting trapped in local minima.
 
Fightfish said:
I usually prefer the alternative - to heat the system up instead. So instead of starting from a random state, I start with the zero temperature limit of perfectly aligned spins. This avoids the problem of getting trapped in local minima.

Not sure if I've got it. So, you start your system at a temperature of approximate 0 ( not exactly zero because of the exponential of delta E / T ), then you heat the system until what temperature ? It cannot pass 2.6..something Kelvin because i want to observe the consequence of phase transition.
 
Not sure what it is that you want to achieve but here is a set of pictures of representative states of a Ising system at different temperatures that I've done in the past:
download.png


You can either start with the uniformly aligned state each time or use the final output system state of the lower temperature as the starting state for the next (higher temperature) simulation.

However, it must be noted that these are just representative states of the system. In order to "observe" the phase transition, one has to measure other quantities such as the average magnetization and the susceptibility, which are averaged over the Monte Carlo runs.
 
From the BCS theory of superconductivity is well known that the superfluid density smoothly decreases with increasing temperature. Annihilated superfluid carriers become normal and lose their momenta on lattice atoms. So if we induce a persistent supercurrent in a ring below Tc and after that slowly increase the temperature, we must observe a decrease in the actual supercurrent, because the density of electron pairs and total supercurrent momentum decrease. However, this supercurrent...
Hi. I have got question as in title. How can idea of instantaneous dipole moment for atoms like, for example hydrogen be consistent with idea of orbitals? At my level of knowledge London dispersion forces are derived taking into account Bohr model of atom. But we know today that this model is not correct. If it would be correct I understand that at each time electron is at some point at radius at some angle and there is dipole moment at this time from nucleus to electron at orbit. But how...
Back
Top