- #1
Aidan1
- 3
- 0
1. The problem:
Building a (probably very simple) computational model for Ising Spins - particles on a lattice with spin up and spin down, only nearest neighbour interactions. I can't for the life of me figure out the renormalisation constant for the probability of a given particle flipping it's spin.
We are given that the probability of a given particle flipping its spin goes as [tex] P \propto \exp{( -\Delta E/T )} [/tex] for [itex] \Delta E > 0 [/itex] where [itex] \Delta E [/itex] is the change in energy if the particle flips its spin and T is the temperature of the system.
In my code I have set [itex] P = k\exp{( -\Delta E/T )} [/itex] & set k = 1 for now while I check other aspects of the code run okay. I can't figure out an expression for k though - I am aware something's going to need to add up to 1 - tried saying [tex] k\exp{(-\Delta E/T )} + P( \Delta E < 0) + ( 1 - k\exp{(-\Delta E/T )} ) * P(\Delta E > 0) = 1 [/tex] (ie probability of particle flipping spin + probability of not flipping spin = 1 - if [itex] \Delta E < 0 [/itex] the probability of flipping is 1) & just got trivial solutions.
Thanks for the help
Building a (probably very simple) computational model for Ising Spins - particles on a lattice with spin up and spin down, only nearest neighbour interactions. I can't for the life of me figure out the renormalisation constant for the probability of a given particle flipping it's spin.
Homework Equations
We are given that the probability of a given particle flipping its spin goes as [tex] P \propto \exp{( -\Delta E/T )} [/tex] for [itex] \Delta E > 0 [/itex] where [itex] \Delta E [/itex] is the change in energy if the particle flips its spin and T is the temperature of the system.
The Attempt at a Solution
In my code I have set [itex] P = k\exp{( -\Delta E/T )} [/itex] & set k = 1 for now while I check other aspects of the code run okay. I can't figure out an expression for k though - I am aware something's going to need to add up to 1 - tried saying [tex] k\exp{(-\Delta E/T )} + P( \Delta E < 0) + ( 1 - k\exp{(-\Delta E/T )} ) * P(\Delta E > 0) = 1 [/tex] (ie probability of particle flipping spin + probability of not flipping spin = 1 - if [itex] \Delta E < 0 [/itex] the probability of flipping is 1) & just got trivial solutions.
Thanks for the help
Last edited: