Mean field approximation and entropy

  • Thread starter CAF123
  • Start date
  • #1
CAF123
Gold Member
2,914
88

Homework Statement


Consider a D dimensional Ising model with N sites, defined by the Hamiltonian $$\mathcal H = -J \sum_{\langle i j \rangle} \sigma_i \sigma_j - h \sum_i \sigma_i$$ where the sum extends over nearest neighbours and each spin variable ##\sigma_i = \pm 1##. For a given spin configuration we denote the number of spins up and spins down by ##N_+## and ##N_-## respectively. The magnetisation is defined by ##m=(N_+ - N_-)/N##

a) Using Stirling's approximation, show that the entropy of the system can be written as $$S/N = \log 2 - \frac{1}{2} (1+m) \log(1+m) - \frac{1}{2} (1-m) \log(1-m)$$ and that in the mean field approximation $$f(T,m) = -\frac{1}{2} Jzm^2 + \frac{1}{2} T((1+m) \log(1+m) - (1-m) \log(1-m)) - T\log 2$$

Homework Equations


$$S = \ln \Omega,$$ in units of ##k_B=1##

The Attempt at a Solution


There are N sites and on each site, the spin can be up or down. So isn't the total number of spin configurations (or number of microstates) just ##2^N##? This would give ##S = N \log 2##, the first term in that expansion but I don't see how the other terms come about.

Thanks!
 

Answers and Replies

  • #2
Charles Link
Homework Helper
Insights Author
Gold Member
4,657
2,007
This doesn't seem like an introductory homework, but in any case I got the first part to work by using ## \Omega=N!/((N_+)!(N_-)!) ## . The number of states is apparently given by the number of possible combinations that can create a given ## m ##. It's an exercise in algebra, but their answer for S is correct.
 
  • #3
CAF123
Gold Member
2,914
88
Hi Charles Link,
This doesn't seem like an introductory homework, but in any case I got the first part to work by using ## \Omega=N!/((N_+)!(N_-)!) ## . The number of states is apparently given by the number of possible combinations that can create a given ## m ##. It's an exercise in algebra, but their answer for S is correct.
Ah ok. I tried this expression for ##\Omega## and I did not seem to get the correct result - in particular where would the factor of log 2 come from?

Thanks!
 
  • #4
Charles Link
Homework Helper
Insights Author
Gold Member
4,657
2,007
The factor of ## \ln 2 ## (it's actually ## N \ln 2 ## ) is what is left over after a lot of terms cancel. If I remember correctly ## N_+=N(1+m)/2 ## and ## N_-=N(1-m)/2 ## There are several places where you use ## N=N_+ + N_- ## and when you take ## \ln(\Omega) ## , Stirling's formula is used 3 times. (twice in the denominator). If you are careful with the algebra, I think you will get the same result. (I think the ## N \ln(2) ## comes from the 1/2's in the denominators of ## N_+ ## and ## N_- ##, but I'd need to look again at my calculations to be sure.)
 
Last edited:
  • Like
Likes CAF123
  • #5
CAF123
Gold Member
2,914
88
I see. thanks! The only thing I don't understand is why would the number of states be counted in this way? If I just count the number of spin configurations then the total number of states would be ##2^N## if I can distinguish between the spins. Then this would yield an entropy with only the first term present. Any ideas?
 
  • Like
Likes Charles Link
  • #6
Charles Link
Homework Helper
Insights Author
Gold Member
4,657
2,007
I had Statistical and Thermal physics (we used the book by F. Reif), quite a number of years ago. If the spins are in a lattice, I think your way of counting them might be appropriate, but I am certainly no expert. I picked the method of the combinations simply because it then resulted in the book's answer.
 
  • #7
CAF123
Gold Member
2,914
88
I had Statistical and Thermal physics (we used the book by F. Reif), quite a number of years ago. If the spins are in a lattice, I think your way of counting them might be appropriate, but I am certainly no expert. I picked the method of the combinations simply because it then resulted in the book's answer.
I see, thanks anyway! ok maybe we can await an answer from someone else, @TSny maybe? I thought my spins were on a lattice no? I imagined the Ising model in D dim to be like a hypercubic with spins placed on the nodes of the cubes.
 
  • #8
TSny
Homework Helper
Gold Member
12,751
3,105
The spins are on a lattice.

2N represents the total number of possible micro states. You want the entropy corresponding to a particular value of magnetization m. Only a subset of the 2N possible micro states have the magnetization m. Charles gave you an expression for the number of micro states Ω in the subset. If you use his expression you should get the desired expression for S/N corresponding to a magnetization m.
 
  • Like
Likes CAF123 and Charles Link

Related Threads on Mean field approximation and entropy

Replies
1
Views
571
Replies
3
Views
444
Replies
3
Views
13K
Replies
1
Views
604
Replies
4
Views
3K
Replies
1
Views
483
  • Last Post
Replies
12
Views
1K
  • Last Post
Replies
5
Views
5K
  • Last Post
Replies
6
Views
5K
  • Last Post
Replies
2
Views
1K
Top