PHZ3400-11 Problem Set 3

Problem 1 - Critical Behavior of a Ferromagnet

a) Show that the mean-field equation of state for a ferromagnetic Ising model is given by (use units where $k_{B}=1$; $\beta =1/T$):

$m=\tanh \left(\beta Jzm+\beta h\right)$.

b) Consider $h=0$, and solve this equation numerically (using Mathematica, MAPLE, or FORTRAN), and plot $m(T)$. At $T<T_{c}=Jz$, you will find three solutions. Pick the one with positive magnetization and plot it as a function of temperature. Show the printout of your computer plot as a part of the homework.

c) By expanding the RHS in powers of $m$, to cubic order, solve analytically for the magnetization close to $T_{c}$. Plot the expression you obtain on the same graph as the full numerical solution (dashed lines). The two curves should agree only close to $T_{c}$

d) Now choose a small but finite external magnetic field $h>0$. Solve the above equation numerically, again picking the positive solution. You should find that transition is "rounded" and the magnetization remains finite at any temperature.

e) Now chose $T=T_{c}=Jz$, and expand the above equation to linear order in $h$, and to cubic order in $m$. Solving for $m$ as a function of $h$, show that at the critical point

$m\propto h^{1/\delta }$, where the critical exponent $\delta =3$.

Problem 2 Gap to low temperature excitations

At $T=0$, our equation of state in zero external field gives fully saturated magnetization $m=1$ (consider the positive solution). At very small but finite temperature, there is a very small temperature correction of the form

$\delta m(T)=1-m(T)\propto \exp\{-E_{g}/T\}$. Determine the energy of the gap $E_{g}$to low-temperature excitations, in terms of the interaction $J$, and the coordination number $z$. How does it compare to the energy to perform a single spin-flip from the ground state (all spins aligned)?

Problem # - Hysteresis

Consider the same equation os state for $T=0.5Jz<T_{c}$. Solve numerically the equation os state starting at large field (say $h=3$, and plot the solution as a function of $h$. You should find that the magnetization, which is close to $m=1$ for large field, decreases as the field is rediced, but remains finite at $h=0+$. Continue plotting it by further reducing the field, until you reach the (negative) coercive field $h_{-}$, where the positive solution can no longer be found. Repeat the calculation starting at very large negative fields and finding the negative magnetization solution, which can be found all the way to the (positive) coercive field $h_{+}$. In the field range between the coercive fields, there are two solutions and we find the hysteresis loop, which is typically seen at any first-order transition.