Chapter 8: Q. 8.19 (page 345)
The critical temperature of iron is . Use this value to make a rough estimate of the dipole-dipole interaction energy , in electron-volts.
Short Answer
The dipole interaction energy of iron =
Chapter 8: Q. 8.19 (page 345)
The critical temperature of iron is . Use this value to make a rough estimate of the dipole-dipole interaction energy , in electron-volts.
The dipole interaction energy of iron =
All the tools & learning materials you need for study success - in one app.
Get started for freeIn this problem you will use the mean field approximation to analyse the behaviour of the Ising model near the critical point.
(a) Prove that, when
(b) Use the result of part (a) to find an expression for the magnetisation of the Ising model, in the mean field approximation, when T is very close to the critical temperature. You should find (not to be confused with 1/kT) is a critical exponent, analogous to the f defined for a fluid in Problem 5.55. Onsager's exact solution shows that in two dimensions, while experiments and more sophisticated approximations show that in three dimensions. The mean field approximation, however, predicts a larger value.
(c) The magnetic susceptibility is defined as . The behaviour of this quantity near the critical point is conventionally written as , where y is another critical exponent. Find the value of in the mean field approximation, and show that it does not depend on whether T is slightly above or slightly below Te. (The exact value of y in two dimensions turns out to be 7/4, while in three dimensions .)
Show that, if you don't make too many approximations, the exponential series in equation includes the three-dot diagram in equation . There will be some leftover terms; show that these vanish in the thermodynamic limit.
Modify the ising program to compute the average energy of the system over all iterations. To do this, first add code to the initialise subroutine compute the initial energy of the lattice; then, whenever a dipole is flipped, change the energy variable by the appropriate amount. When computing the average energy, be sure to average over all iterations, not just those iterations in which a dipole is actually flipped (why?). Run the program for a 5 x 5 lattice for T values from 4 down to l in reasonably small intervals, then plot the average energy as a function of T. Also plot the heat capacity. Use at least 1000 iterations per dipole for each run, preferably more. If your computer is fast enough, repeat for a 10x 10 lattice and for a 20 x 20 lattice. Discuss the results. (Hint: Rather than starting over at each temperature with a random initial state, you can save time by starting with the final state generated at the previous, nearby temperature. For the larger lattices you may wish to save time by considering only a smaller temperature interval, perhaps from 3 down to 1.5.)
Modify the ising program to simulate a one-dimensional Ising model.
(a) For a lattice size of 100, observe the sequence of states generated at various temperatures and discuss the results. According to the exact solution (for an infinite lattice), we expect this system to magnetise only as the temperature goes to zero; is the behaviour of your program consistent with this prediction? How does the typical cluster size depend on temperature?
(b) Modify your program to compute the average energy as in Problem 8.27. Plot the energy and heat capacity vs. temperature and compare to the exact result for an infinite lattice.
(c) Modify your program to compute the magnetisation as in Problem 8.28. Determine the most likely magnetisation for various temperatures and sketch a graph of this quantity. Discuss.
The Ising model can be used to simulate other systems besides ferromagnets; examples include anti ferromagnets, binary alloys, and even fluids. The Ising model of a fluid is called a lattice gas. We imagine that space is divided into a lattice of sites, each of which can be either occupied by a gas molecule or unoccupied. The system has no kinetic energy, and the only potential energy comes from interactions of molecules on adjacent sites. Specifically, there is a contribution of to the energy for each pair of neighbouring sites that are both occupied.
(a) Write down a formula for the grand partition function for this system, as a function of , T, and p.
(b) Rearrange your formula to show that it is identical, up to a multiplicative factor that does not depend on the state of the system, to the ordinary partition function for an Ising ferromagnet in the presence of an external magnetic field B, provided that you make the replacements and . (Note that is the chemical potential of the gas while uB is the magnetic moment of a dipole in the magnet.)
(c) Discuss the implications. Which states of the magnet correspond to low density states of the lattice gas? Which states of the magnet correspond to high-density states in which the gas has condensed into a liquid? What shape does this model predict for the liquid-gas phase boundary in the P-T plane?
What do you think about this solution?
We value your feedback to improve our textbook solutions.