Project 4: Difference between revisions

From phys660
Jump to navigationJump to search
No edit summary
No edit summary
Line 4: Line 4:


<math> E = - J \sum_{\langle i,j \rangle} S_i S_j - H \sum_i S_i </math>
<math> E = - J \sum_{\langle i,j \rangle} S_i S_j - H \sum_i S_i </math>
where <math> \langle i,j \rangle </math> means that the sum runs over all pairs of nearest neighbor spins (attached to the lattice sites), and <math> J </math> is the exchange coupling. Positive  <math> J >0 </math> interaction between the spins leads to ferromagnetism, will negative <math> J >0 </math>  gives rise to antiferromagnetism in the Ising model. Choose units so that <math> |J|=1 </math> and the Boltzmann constant is <math> k_B=1 </math>. The free parameters of the model are then just the temperature <math> T </math> and the applied field <math> H </math>.
Parth I

Revision as of 13:13, 2 April 2012

First and second order phase transition in the Ising model of ferromagnetism

Write a program for the 2D Ising model on a 64×64 square lattice with periodic boundary conditions which is placed in an applied external magnetic field H. The energy of a particular state of the system is given by:

E=Ji,jSiSjHiSi

where i,j means that the sum runs over all pairs of nearest neighbor spins (attached to the lattice sites), and J is the exchange coupling. Positive J>0 interaction between the spins leads to ferromagnetism, will negative J>0 gives rise to antiferromagnetism in the Ising model. Choose units so that |J|=1 and the Boltzmann constant is kB=1. The free parameters of the model are then just the temperature T and the applied field H.

Parth I