s i s j µ B H Figure 3.12: A possible spin conguration for an Ising model on a square lattice (in two dimensions).
|
|
- Stephanie Webster
- 5 years ago
- Views:
Transcription
1 s i which can assume values s i = ±1. A 1-spin Ising model would have s i = 1, 0, 1, etc. We now restrict ourselves to the spin-1/2 model. Then, if there is also a magnetic field that couples to each spin, and assuming only nearest-neighbour interactions, we have the Hamiltonian H ({s k }) = J nn s i s j µ B H i s i. (3.51) nn stands for nearest neighbours. The number of nearest neighbours is called coordination number, z. In a two-dimensional square lattice z = 4. The coupling constant between spins is J, which here is assumed to be constant (i.e. not to depend on the location of the spins). In turn, we may have J > 0 or < 0. In Fig a possible configuration for a spin-1/2 Ising model on a square lattice is represented. Many variations of this model are possible. For example, the coupling constant J can be site-dependent. The Heisenberg model introduces a sophistication by way of a vector spin vector: H ({s k }) = J ij s i s j µ B H s i. (3.52) s i are freely-rotating vectors in the 3D space. One variation of this is the XY model, which is Heisenberg s model on the xy plane. The Gaussian, spherical, Potts and percolation models are also frequently used. The spin-1/2 Ising model with J > 0 is a toy model for i Figure 3.12: A possible spin conguration for an Ising model on a square lattice (in two dimensions). the ferromagnetic transition. The ferromagnetic transition occurs in ferromagnets, and is associated with spontaneous (i.e. at zero magnetic field) ferromagnetic order: above the critical temperature T c the magnetisation of the material is zero; below the critical temperature a spontaneous magnetisation arises: 0, T < T c, m(t) = 0, T T c. (3.53) 68
2 The transition is of second order at zero field, but it can also be reached at T < T c by varying the magnetic field from H < 0 to H > 0, and in this case it is of first order. Fig shows a schematic phase diagram in the T m and H T plane; note that the phase diagram is similar to that of the liquid-vapour phase transition, with H playing the role of pressure and the magnetisation per spin, m = M/N, m = N s i N = s i, (3.54) playing the role of density. The response function for a magnetic system is the magnetic susceptibility at zero field: ξ = m. (3.55) H H=0 As corresponds to a second-order transition, ξ [a second derivative of the free energy G(H, T)] diverges as T T c + at zero field. First derivatives, for example, the magnetisation M = ( G/ H) T, are finite and continuous at H = 0. Another quantity of interest is the correlation function, the analogue of the total correlation function h(r) for a fluid: Γ ij = Γ i j = s i s j s i s j = s i s j m 2. (3.56) This function remain short-ranged at H = 0 (it is an exponentially-decaying function of the separation i j ) except at T = T c, where it becomes a long-ranged function. It is easy to demonstrate that ξ = β N ( s i s j s i s j ) = β ij j ( s1 s j m 2) (3.57) [again this is analogous to the relation between κ T and the integral of h(r) for a fluid]. Figure 3.13: (Colour on the Web). Schematic phase diagrams in the m T and H T plane of the spin-1/2 Ising model. The shaded area represents the region of twophase coexistence. 69
3 3.4.1 Mean-field approximation for the 2D-Ising model: Weiss theory The Ising model in 2D is a cornerstone in the history of statistical mechanics as it is exactly solvable and contains a non-trivial phase transition (the ferromagnetic transition). This was shown by Onsager. The derivation of the exact solution is a bit cumbersome though. Therefore, we will only mention the exact solution without explicitely obtaining it. Here we resort to the mean-field approximation, which comes in various disguises. One of them is the Weiss or effective-field theory. The Weiss approximation consists of assuming that each spin (for instance the i-th spin) interacts not with the neighbouring spins in each system s conguration, {s j }, but with their averaged values, s j = m (Fig. 3.14). The energy of the i-th spin, ǫ i, is therefore written as ǫ i = J j i s i s j µ B Hs i = (Jzm + µ B H)s i µ B H eff s i. (3.58) The total energy is E = NJzm 2 /2 Nµ B Hm (the factor 1/2 accounts for the double counting of pair interactions). The underlying assumption is that the system is ideal (i.e. spins are non-interacting): every spin, with average magnetisation m, interacts with an external effective eld given by µ B H eff = Jzm + µ B H, equal to the real external field H plus an internal averaged field due to the neighbouring spins. But it is not an ideal approximation completely, since the effective field depends in turn on the average value of the spin, m: it is a self-consistent approximation, which expresses itself by means of the trascendental equation in m that we are about to derive now. Since spins are in principle ideal, the average value of s i can be obtained using its associated reduced Boltzmann factor exp (βǫ i ): m = s i e βǫ i s i =±1 e βǫ i s i =±1 = eβ(jzm+µ BH) e β(jzm+µ BH) e β(jzm+µ BH) + e β(jzm+µ BH) = tanh [β(jzm + µ BH)]. (3.59) This is a trascendental (self-consistent) equation in m [compare with Eqn. (1.83) where spin coupling was neglected]. Note that this equation does not explicitely depend on the dimensionality of the system but through the coordination number z. This is a typical feature of the mean-field approximation Mean-field approximation for the 2D-Ising model: Weiss theory Eqn. (3.59) has to be solved numerically. We consider the case H = 0. A simple means of obtaining the solution is based on inspection of Fig In the figure the straight line f(m) = m and the curve f(m) = tanh(βjzm) are plotted; the solution of Eqn. (3.59) is the point(s) where the two functions cross. It is immediately evident that, depending on the value of the temperature, two types of solution arise: 70
4 Figure 3.14: (Colour on the Web). A fluctuating spin s 0 embedded in a medium of spins with average value of their spin m. If the slope of the function f(m) = tanh (βjzm) at the origin is less than unity, then only one solution, m = 0, exists: The stable phase is the disordered phase (zero magnetisation). If the slope of the function f(m) = tanh (βjzm) at the origin is larger than unity, then two solutions, m = ±m 0 0, exist, corresponding to a degenerate ordered phase. The slope is: ( ) d Jz [1 dm [tanh (βjzm)] m=0 = tanh 2 (βjzm) ] kt = Jz m=0 kt 1. (3.60) The critical temperature T c, separating the two regimes of solutions, is therefore: kt c J = z. (3.61) As the number of neighbours increases, the critical temperature increases; this was to be expected, since more neighbours means stronger interactions and increasing ability of the system to counteract the disordering effect of temperature via fluctuations. For example, in the triangular lattice (2D), we have z = 6, and T c will be larger than that for the square lattice, with z = 4. For H 0 the solution of Eqn. (3.59) is as follows. We can write { [ ( f(m) = tanh [β(jzm + µ B H)H] = tanh βjz m µ BH Jz )]}, (3.62) so that the hyperbolic-tangent function in Fig is displaced to the left if H > 0 or to the right if H < 0, giving one solution with m > 0 and possibly two solutions with m < 0 in the first case and one solution with m < 0 and possibly two solutions with m > 0 in the second. This implies the equation of state H = H(m, T) shown in Fig. 3.16, which can in fact be obtained using Eqn. (3.59) and solving for H: βµ B H = arctanh m βjzm. (3.63) 71
5 Three temperatures are shown. The equation of state for T > T c is similar to that of a paramagnetic material. For T < T c a loop appears; this is similar to the van der Wals loop in the case of fluids. Here again we can interpret that the true phase transition proceeding at H = 0 between two states with magnetisation ±m(t) (T < T c ) is in fact given by the meanfield theory but a Maxwell construction must be applied (equality of free energies G 1 = G 2 in the two phases). In this case, due to the symmetry of the function H = H(m), the construction leads immediately to H = 0 for the location of the phase transition. There is a region of unstable states with H/ m < 0 (which implies a negative magnetic susceptibility, i.e. mechanical instability) and two regions of mechanically stable but thermodynamically metastable states. Figure 3.15: Graphical solution of Eqn. (3.59): points where the straight line f(m) = m (dashed line) and the curve f(m) = tanh (βjzm) (continuous lines) cross for three temperatures: kt/jz = 1/2 (T < T c, with solution at ±m), 1 (T = T c also with solution at m = 0) and 2 (T > T c, also with solution at m = 0). Solutions are indicated with filled circles. average value of their spin m. Based on all these considerations, we can represent the phase diagram in the T m plane, shown in Fig The dashed line is the line of spinodal points (where the magnetic susceptibility becomes infinite). The region inside the coexistence curve is the two-phase region, where the two phases coexist at composition given by the analogue of the lever rule. In this region we can distinguish between the region of spinodal decomposition inside the spinodal line, and the two regions of nucleation, between the coexistence and 72
6 spinodal lines (analogous regions exist for the liquid-vapour transition). These words refer to the different kinetic mechanisms by which the phase transition proceeds, a theme not covered by our equilibrium mean-field theory. In the nucleation region domains of the material with reversed magnetisation appear; they have to overcome an activation barrier but, once overcome, the domains grow without limit until all the material has reversed the magnetisation. In the spinodal-decomposition region there is no activation for the formation of domains and there is a fluctuation in the whole system, not just in local areas, that drives the transition. The prediction of meanfield theory is that there is a Figure 3.16: Equation of state H = H(m, T) for the spin1/2 Ising model, in meanfield approx- imation, for three temperatures: kt/jz = 2.2 (T < T c ), 1 (T = T c ) and 6 (T > T c ). phase transition at any dimensionality. This is wrong: there is no such transition in a 1D ferromagnet. To see this using a simple argument, let us consider the ground state (T = 0) of a row of spins (all pointing in one direction, say up ). The energy of this state is E 0 = NJ. We now consider a simple fluctuation from the ground state where half of the spins have flipped from up to down (thus creating a grain boundary). The energy is now E 1 = E 0 +2J and the excess energy per spin is e = (E 1 E 0 )/N = 2J/N. But for large N this is such a small energy that even for very low T these type of fluctuations, which destroy the ground state, will very likely happen: the ordered state is only possible at T = 0 or, if we want, the phase transition occurs at T = 0. In fact, including the entropy S = k log (N 1) [there are N 1 ways to locate the grain boundary], the associated free energy is G = 2J kt log (N 1), which is negative in the thermodynamic limit for all 73
7 temperatures. However, mean-field theory predicts a phase transition at kt c /J = z = 2! Incidentally, the same type of argument in 2D (i.e. the creation of grain boundaries) leads to a prediction of a finite critical temperature, actually quite close to the exact result. Let us now discuss the spin-1/2 Ising model with antiferromagnetic coupling. At H = 0, Figure 3.17: Phase diagram of the spin-1/2 Ising model, in meanfield approximation, in the temperature-magnetisation per spin plane, T m, at zero magnetic field. Temperature is measured in units of the critical temperature. The critical point is depicted by a circle. The shaded region is the twophase region, where the system is not thermodynamically stable and separates into two equivalent phases with opposite magnetisation. The dashed lines indicate the spinodal line where the magnetic susceptibility is infinity. the completely ordered phase on the square lattice consists of two interpenetrating square lattices, a and b, with opposite magnetisations m a = m b = 1. In the disordered phase m a = m b = 0 (Fig. 3.18). The condition m a + m b = 0 is in fact true for any T because of symmetry. Spins in a given sublattice do not in fact interact (for nearest-neighbour interactions) and are only subject to interactions with spins of the other sublattice. In mean-field approximation the Hamiltonian can be written H = J s i s j = n.n. i s (a) i H (a) + j s (b) j H (b), H (a) = Jzm a, H (b) = Jzm b. (3.64) Then: m a = s (a) i = s (a) i =±1 s (a) i =±1 s (a) i e βh(a) e βh(a) = tanh (βjzm b )tanh (βjzm a ), 74
8 m b = s (b) s (b) i = i =±1 s (b) i =±1 s (b) i e βh(b) e βh(b) = tanh (βjzm a ) tanh (βjzm b ). (3.65) Therefore, in mean-field approximation, each of the sublattice magnetisations, m a and m b, depend on T following the same law as in the ferromagnetic case. In particular, there is a phase transition between the disordered phase and the ordered phase (where the magnetisation of the sublattices are non-zero and with opposite signs) at a critical temperature kt c /J = z. This is ultimately based on the symmetry of H with respect to the two sublattices, a symmetry which is lost when H 0. In this case the transition temperature decreases until it becomes zero beyond some value of the magnetic field (Fig. 3.18). Figure 3.18: (Colour on the Web). On the left, phase diagram of the antiferromagnetic spin-1/2 Ising model, in the H T plane. The dashed line is a line of critical points (second-order phase transitions). On the right, depiction of the ordered and disordered phases, indicating the population of the two sublattices a and b. Let us calculate the critical exponents. The exponent β can be easily calculated from the trascendental equation (3.19). Setting H = 0 and recognising that close to the critical point m 1, we can expand the hyperbolic-tangent function: Then, since kt c = Jz: 1 T c T 1 3 m = tanh (βjzm) = βjzm 1 3 (βjzm)3 + O ( m 5). (3.66) ( ) 3 Tc m 2 m ± ( ) T 3/2 ( ) Tc T 1/2 3 (T c T) 1/2, (3.67) T Tc T so that β = 1/2, the same result obtained with the van der Waals equation for a fluid. The analogue of the compressibility for a fluid is here the magnetic susceptibility at zero field, ξ: ξ = m (T T H c ) γ, T > T c. (3.68) H=0 75
9 Differentiating (3.59) with respect to H: m H = { 1 tanh 2 [β (Jmz + µ B H)] } ( βjz m ) H + µ BH ( = (1 m 2 ) βjz m ) H + µ BH. (3.69) Solving for the derivative, and calculating at H = 0 (which means m = 0 for T > T c ): m = βµ B(1 m 2 ) H H=0 1 βjz(1 m 2 = βµ BT (T T ) c ) 1, (3.70) m=0 T T c which means γ = 1, again as in van der Wals theory. All other critical exponents can be seen to be equal to those derived from van der Waals theory (and in fact from any mean-field theory). The conclusion is that, at the level of mean-field theory, fluids and magnets are in the same universality class. 3.5 Improvement of Weiss theory: Bethe theory In Bethe theory correlations of two spins are taken into account. This is realised by considering fluctuations of one spin and of its z neighbours, i.e. a cluster of z + 1 spins rather than just one-spin clusters, as in Weiss theory. The z neighbouring spins are still coupled to their neighbours through the average magnetisation per spin m. For a 2D lattice in zero external field (H = 0) the energy associated to such a cluster is (see Fig. 3.19): z ǫ cluster = Js 0 s i (z 1)mJ Now the partition function of the cluster is Z cluster = e βj(s 0+(z 1)m)(s 1 +s 2 + +s z) s 0,s 1,,s z=±1 z s i. (3.71) = 2 z cosh z {βj[1 + (z 1)m]} + 2 z cosh z {βj[1 (z 1)m]}. (3.72) The average value of the central spin is s0 = 1 Z cluster [2 z cosh z {βj[1 + (z 1)m]} 2 z cosh z {βj[1 (z 1)m]}], (3.73) whereas the average value of one of its neighbours is s j = 1 Z cluster [ 2 sinh {βj[1 + (z 1)m]}2 z 1 cosh z 1 {βj[1 + (z 1)m]} 2 sinh {βj[1 (z 1)m]}2 z 1 cosh z 1 {βj[1 (z 1)m]} ]. (3.74) Since these two have to be the same for consistency, s 0 = s j, and we arrive at the equation e 2βJ(z 1)m = coshz 1 {βj[1 + (z 1)m]} cosh z 1 {βj[1 (z 1)m]} 76 (3.75)
10 Analysis of this equation can be done graphically, or else as follows: m = 0 is always a solution, so the two functions cross at this point. Since both functions are monotonically increasing but the exponential diverges whereas the ratio of hyperbolic cosines has a horizontal asymptote, they can only cross at an additional point (symmetric with respect to the origin in view of the symmetry m m of the equation) if the slope of the latter is larger or equal to the slope of the former, which is 2βJ(z 1); the condition is then: 2βJ(z 1) 2 tanhβj 2βJ(z 1) β c J = arctanh (z 1) 1. (3.76) For the square lattice (z = 4) this yields β c J = or kt c /J = , to be compared with the exact result kt c /J = and with the mean-field result ktc/j = 4. This is a considerable improvement over Weiss theory. However, it can be shown that the critical exponents derived from Bethe s theory have the classical (mean-field) values: Bethe s theory is a mean-field theory, more sophisticated than Weiss, but still mean field. Figure 3.19: (Colour on the Web). A fluctuating cluster of spins (in red) embedded in a medium of spins with average value of their spins m. 3.6 Mean-field approximation: Bragg-Williams theory The mean-field theory of Weiss can be obtained in a completely different way by writing an effective free energy that is minimised with respect to the magnetisation per spin m. This procedure is based on the result that the thermodynamic free energy is a minimum, at fixed values of the thermodynamic external parameters, with respect to the internal, microscopic variables. Here the internal variable is taken to be m. We construct the free energy, G(H, T; m) = E(H, T; m) TS(H, T; m) as follows. First, we have seen that the energy, in mean-field approximation, is E(H, T; m) = 1 2 NJzm2 Nµ B Hm. (3.77) 77
11 Now we write an entropy. Within mean-field approximation spins are uncorrelated, so that the entropy can be written as S = k log Ω, where Ω is the number of ways to impose a given value of the average magnetisation per spin, m, in a system of independent spins. If N + and N are the number of up and down spins, respectively, and if M = Nm is the magnetisation: Then: M = N + N, N = N + + N. (3.78) N + = 1 2 (N + M) = N 2 (1 + m), N = 1 2 (N M) = N (1 m), (3.79) 2 so that Ω = N! N +!N! = N! [ ] [ ]. (3.80) N N 2 (1 + m)! 2 (1 m)! After short manipulations, the entropy can be written as S Nk = 1 + m ( ) 1 + m log 1 m ( ) 1 m log. (3.81) The free energy is then: G(H, T; m) = 1 2 NJzm2 Nµ B Hm [ ( ) 1 + m 1 + m + NkT log m 2 ( 1 m log 2 )]. (3.82) The equilibrium state (value of m) of the system is obtained from: G = 0. (3.83) m Then: Rearranging the equation: eq. G m = NJzm Nµ BH + NkT m 1 m = e2β(jzm+µ BH) log ( ) 1 + m = 0. (3.84) 1 m m = tanh [β (Jzm + µ B H)], (3.85) which is the same trascendental equation for m as that obtained in the Weiss theory. 3.7 Ising model in 1D (Ising chain) In 1D the Ising model can be solved analytically. Let us deal with the case H = 0 first. The Hamiltonian is H = J N 1 78 s i s i+1, (3.86)
12 Figure 3.20: (Colour on the Web). Left: free boundary conditions (spins s 1 and s N only have one neighbour). Right: periodic boundary conditions (all spins have two neighbours). where we used free boundary conditions (i.e. the s 1 and s N spins only have one neighbour, s 2 and s N 1 respectively). The partition function is: Z N = s 1 =±1 s 2 =±1 s N 1 =±1 s N =±1 βj e N 1 s i s i+1. (3.87) The sum over s N can be done explicitely since s N only appears once in the sum: s N =±1 independent on the value of s N 1. Then we can write: e βjs N 1s N = 2 cosh βj, (3.88) Z N = (2 cosh βj)z N 1 Z N = (2 cosh βj) N 2 Z 2, (3.89) with Z 2 = e βjs 1s 2 = 4 cosh βj, (3.90) s 1 =±1 s 2 =±1 so finally Z N = 2(2 coshβj) N 1, (3.91) and the free energy is G = NkT log Z N = kt [log 2 + (N 1) log (2 coshβj)] = NkT log (2 cosh βj). (3.92) The last equality applies in the thermodynamic limit. In the case H 0 it is more useful to write the Hamiltonian using periodic boundary conditions: N H = J s i s i+1 1 N 2 µ BH (s i + s i+1 ). (3.93) 79
13 The partition function is then: Z N = s 1 =±1 s 2 =±1 = s 1 =±1 s 2 =±1 s N =±1 β e N N e β s N =±1 [ Js i s i ] 2 µ BH(s i + s i+1 ) [ Js i s i µ BH(s i + s i+1 )]. (3.94) At this point the use of transfer matrices facilitates manipulations. We define the matrix: ( ) ( P11 P P = 1, 1 e β(j+µ B H) e βj ), (3.95) P 1,1 P 1, 1 e βj e β(j µ BH) in terms of which the partition function can be written: Z N = { N } P si s i+1 = P s1 s 2 P s2 s 3 P sn s 1 s 1 =±1 s 2 =±1 s N =±1 s 1 =±1 s 2 =±1 s N =±1 = ( ) P N = Tr P N. (3.96) s 1 =±1 Now λ 1, λ 2 are the eigenvalues of P, then λ N 1, λn 2 are the eigenvalues of PN and The eigenvalues λ 1, λ 2 are: The free energy is λ 1,2 = e βj cosh (βµ B H) ± = e βj cosh (βµ B H) ± G = kt log ( λ N 1 + λ N 2 Z N = Tr P N = λ N 1 + λn 2. (3.97) e 2βJ cosh 2 (βµ B H) + 4 cosh βj sinh βj e 2βJ sinh 2 (βµ B H) + e 2βJ. (3.98) ) = kt N log λ 1 + log ( ) N λ2 1 + λ 1 NkT log λ 1 as N (3.99) since λ 1 > λ 2. The free energy is then: [ ] G = NkT log e βj cosh (βµ B H) + e 2βJ sinh 2 (βµ B H) + e 2βJ. (3.100) This expression reduces to Eqn. (3.92) at H = 0. The magnetisation is: m = s i = 1 N G H = sinh (βµ B H). (3.101) sinh 2 (βµ B H) + e 4βJ 80
14 From here we see that, at zero field, there can be no spontaneous magnetisation. The magnetic susceptibility is ξ = m β cosh (βµ B H) = = βe H H=0 (1 + e 4βJ sinh (βµ B H)) sinh 2βJ. (3.102) 2 (βµ B H) + e 4βJ ξ never diverges at finite T, but becomes arbitrarily large as T 0, indicating that we can assume a kind of phase transition to take place at this temperature (even though the ordered state can never be reached!). Next we calculate the correlation function: Γ j = s i s i+j s i s i+j = s i s i+j (3.103) (we used the fact that s i = 0 for all temperatures). First, we consider the more general function s i s i+j, and assume i and i + j to be far from both ends of the chain. Second, we take the ferromagnetic coupling constant of the pair l and l + 1 to be J l, a constant that will be set to J at the end of the calculations. The partition function will be We have s i s i+j l = 1 Z N = 1 Z N s 1 =±1 s 2 =±1 s 1 =±1 s 2 =±1 s N =±1 N 1 Z N = 2 (2 cosh βj l ). (3.104) s N =±1 l=1 β s i s i+j e N 1 J l s i s i+1 β (s i s i+1 )(s i+1 s i+2 ) (s i+2 s i+3 ) (s i+j 1 s i+j ) e The last equality follows from the property s 2 i = 1. Now consider the derivative: N 1 j Z N j β J l s i s i+1 = e J i J i+1 J i+j 1 J i J i+1 J i+j 1 s 1 =±1 s 2 =±1 s N =±1 = β j Therefore: s 1 =±1 s 2 =±1 s i s i+j l = s N =±1 ( 2 ) Z N β j β (s i s i+1 ) (s i+1 s i+2 ) (s i+2 s i+3 ) (s i+j 1 s i+j )e j Z N J i J i+1 J i+j 1 N 1 N 1 J l s i s i+1. J l s i s i+1. (3.105) (3.106) 2 Z N β j {cosh βj 1 cosh βj 2 (β sinh βj i ) (β sinh βj i+1 ) (β sinh βj i+j 1 ) β sinh βj N 1 } = tanhβj i tanhβj i+1 tanhβj i+j 1, (3.107) 81
15 where Eqn. (3.91) has been used in the last step. Now specifying at J l = J: The correlation function is s i s i+j = (tanh βj) j. (3.108) Γ j = s i s i+j s i s i+j = s i s i+j = (tanh βj) j. (3.109) Since tanh βj < 1, the correlation function is a decreasing function of distance j (as it should be). We define the correlation length, χ, as tanh βj = e 1/χ χ = [log (tanh βj)] 1, (3.110) so that the correlation function is an exponential: Γ j = e j/χ. (3.111) χ is a positive quantity, and quantifies the exponential decay of Γ. As T 0 +, we have tanh βj 1, and log (tanhβj) 0, so that χ, as would be expected if we were approaching a critical point. 3.8 The lattice gas and lattice-alloy models Here we consider some lattice models that can be formulated to treat a number of problems different from the ferromagnetic systems but that can be seen to map exactly onto a spin-1/2 Ising model. Binary alloy We begin with the binary alloy. Here we have two species A and B, and each lattice site can be occupied by either species. N A and N B are the (in principle fixed) number of particles of each species (no vacancies are allowed). Clearly mapping onto the Ising model can be made by assigning A species to up spin and B species to down spin (or the other way round). The Hamiltonian would be ǫ AA, s i = +1, s j = +1, H = n.n.ǫ si,s j s i s j, ǫ si,s j = ǫ AB, s i = +1, s j = 1 or s i = 1, s j = +1, (3.112) ǫ BB, s i = 1, s j = 1. In the Ising model we would have ǫ AA = ǫ BB = J and ǫ AB = J, so the binary-alloy model seems a bit more general. To see more clearly the possible configurations, we express the energy of the i-th configuration of the system, E i, as follows: E i = ǫ AA N (i) + ǫ AA ABN (i) + ǫ AB BBN (i), (3.113) BB where N kl are the number of kl pairs in the lattice, and kl are the energies associated with kl pairs. If z is the coordination number, it is easy to see that zn A = 2N (i) AA + N AB, (i) zn B = 2N (i) BB + N(i) AB (3.114) 82
16 by the following procedure: take A sites and draw a line to all its z nearest neighbours; then zn A lines will have been drawn, and there will be double lines between AA pairs and single lines between AB pairs (and no lines between BB pairs). Hence the first line of Eqns. (3.114) follows. Doing the same with B sites we have the second line. This means that only one of the three variables N kl are independent, and we choose N (i) AB as the independent one. Then: where E i = 1 2 ǫ AA ( ) zna N (i) AB + ǫabn (i) AB ǫ BB ( ) znb N (i) AB = z 2 ǫ AAN A + z 2 ǫ BBN B + ǫn (i) AB, (3.115) ǫ ǫ AB ǫ AA 2 ǫ BB 2 (3.116) is the mixing energy. The mixing energy gives the tendency for the system to keep different species apart (ǫ > 0), i.e. to macroscopically segregate into two phases (one rich in A s and the other in B s) or together (ǫ < 0), giving a uniform mixture. Note that the same construction can be made for the Ising model; in that case ǫ AA = ǫ BB = J, ǫ AB = J, and ǫ = 2J, so that the ferromagnetic Ising model gives phase separation into +1-rich and 1- rich domains. In the antiferromagnetic Ising model J < 0 and ǫ < 0 which corresponds to antiferromagnetic ordering; this case corresponds to a lattice of alternating A and B species in the binary-alloy case. Therefore, in the alloy language we may have: ǫ > 0: liquid binary mixtures where two phases, one rich in A and the other rich in B phase-separate below a critical (consolute) temperature T c. An experimental T x phase diagram of a mixture of sodium fluoride and aluminium fluoride is shown in Fig. 3.21(a). We define the composition or mole fraction x as x = N A N A + N B = N A N. (3.117) In the case of the figure, x corresponds to the fraction of the complex AlF 3 with respect to the total number. In the figure various transitions, involving liquid and solid phases, are visible, but here it is the one between two liquids, liq 1 and liq 2, that is of interest to us. ǫ < 0: solid binary mixtures where a lattice with compositional disorder (sites occupied with A or B more or less at random) becomes ordered (with sites of one sublattice mostly occupied by a given species in a regular arrangement) below some temperature; this is the so-called order-disorder transition, an example of which is the β-brass phase of CuZn, the ordered phase here being two interpenetrating simple-cubic lattices of Cu and Zn atoms, which are occupied at random above a given temperature, Fig. 3.21(b). In each of these two cases results from the Ising model can be used, since there is an exact mapping between the two models. In particular, the mean-field approximation of Weiss for a ferromagnet can be used to predict phase behaviour in binary fluid mixtures 83
17 Figure 3.21: (a) Phase diagram, in the T x plane, of the mixture NaF/AlF 3, with x being the mole fraction of the second. (b) Structure of β-brass, CuZn. or binary solid alloys. But let us rewrite the theory in terms of the new variable (number of AB pairs N AB or equivalently mole fraction x). Since particles are taken to be independent in mean-field theory, the number of N AB pairs can be obtained as follows: since zn/2 is the total number of bonds, the probability of having an AB pair in the i-th conguration is N (i) AB/(zN/2). We approximate this as N (i) AB 1 2 zn 2 ( NA N ) ( NB N ) = 2 N AN B N 2, (3.118) i.e. as the product of probabilities of having a site occupied by A and the other by B (in fact we know this is not quite correct because of correlations). Then: E i = z 2 ( ǫ AA N A + ǫ BB N B + 2ǫ ) N N AN B = Nz 2 [ǫ AAx + ǫ BB (1 x) + 2ǫx(1 x)]. (3.119) The entropy can be obtained from Eqn. (3.81) directly, taken into account the equivalence between magnetisation m and mole fraction x: m = N + N N N A N B N = 2N A N N = 2x 1. (3.120) Then, taking ǫ AA = ǫ BB for simplicity (which symmetrises the mixture), the grand potential is Ω NkT zǫ = x log x + (1 x) log (1 x) + kt x(1 x) µ A kt x µ B (1 x), (3.121) kt where an unimportant constant term zǫ AA /2kT has been omitted, and contributions from chemical potentials µ A and µ B (which play the role of the magnetic field in the Ising model), giving contributions µ A N A and µ B N B, added. Also, all energies can be measured in units of kt, so we can write ω = x log x + (1 x) log (1 x) + zǫ x(1 x) µ x, (3.122) 84
18 where ω = Ω NkT, ǫ = ǫ kt, and a constant term µ B /kt has been deleted. µ = µ A µ B, (3.123) kt Let us take ǫ > 0 (liquid binary mixture) and analyse the free energy with respect to x. The entropic term, called mixing entropy, S = Nk [x log x + (1x) log (1x)], (3.124) always favours mixing (disorder), having a maximum at x = 1/2 (complete mixing), and discouraging formation of pure phases (x = 0 or x = 1), for which entropy is minimum. The energetic term, E = NzkTǫ x(1 x), ǫ > 0, (3.125) have minima at x = 0 and x = 1, and a maximum at x = 1/2, favouring formation of pure phases (i.e. segregation into Arich and Brich phases). In Fig we plot ω(x) for µ = 0 (making the mixture fully symmetric with respect to x = 1/2) for various temperatures. For ǫ > ǫ c there is a region of instability; the critical point verifies ω x = 0, 2 ω x 2 = 0, log x 1 x + zǫ µ 2zǫ x = 0, 1 x x 2zǫ = 0, (3.126) with solution x c = 1/2 and ǫ c = 2/z or kt c /ǫ = z/2. For T < T c there are three solutions, one at x = 0 (unstable) and two symmetrically placed with respect to x c, at x = 1/2 ± x 0, corresponding to the two (stable) coexisting phases, one rich in A particles, the other rich in B particles. For ǫ < 0 the x variable is not a good one since we need something that measures the sublattice occupation. Once this is identified, one can follow an equivalent number of steps and treat the case of a solid binary alloy in the same vein. Lattice gas Consider now the lattice gas, a model meant to approximately represent gases and liquids. We introduce occupancy variables c i for each lattice site, so that c i = 0 if the lattice site is empty and c i = 1 if it is occupied. Exclusion effects typical of fluids at high density arise implicitely by not permitting more than one particle to occupy any given site. We consider an attractive interaction between nearest-neighbour particles, and write the Hamiltonian as H = ǫ c i c j µn = ǫ c i c j µ n.n. n.n. i c i. (3.127) To see that this model is isomorphic to the spin-1/2 Ising model, we introduce spin variables s i as s i = 2c i 1. (3.128) 85
19 Figure 3.22: Free energy ω as a function of composition x for reduced inverse temperature zǫ = 2.5, 2.0 and 1.5, which correspond, respectively, to temperatures below, equal to and above the critical temperature T c. The dashed line is the two-phase region of coexistence for the temperature. The Hamiltonian is then: H = ǫ 1 n.n. 2 (s i + 1) 1 2 (s j + 1) µ (s i + 1) 2 i = N ( ) zǫ µ ǫ i s j 4 n.n.s 1 ( ) zǫ µ i s i, (3.129) which corresponds to the Ising model with the identification J ǫ and µ B H (zǫ/2 + µ)/2. We can then say that there occurs a first-order liquid-vapour phase transition at a chemical potential µ = zǫ/2 at temperatures below kt c /ǫ = z/4, which is the critical temperature. 3.9 Landau theory As we have seen, in a second-order phase transition second derivatives of the free energy, such as the specific heat, the isothermal compressibility or the magnetic susceptibility, exhibit divergences at the critical point, while the first derivatives (entropy, magnetisation or density) remain continuous. From this point of view, all second-order phase transitions have something in common. In 1937 Landau extended this idea and formulated a theory which was meant to explain all second-order phase transitions from a common viewpoint. The Landau theory first identifies an order parameter, η, a parameter which distinguishes the high symmetry phase, where η = 0, from the low-symmetry phase, for which η 0. The order parameter may be a scalar, vector or tensor with respect to the symmetry 86
20 operations associated to the problem. For the moment we will assume it to be a scalar. Since the transition is continuous, the order parameter is continuous at the critical point (T = T c ). Then the dependences of the free energy are written as G = G(h, T; η) (h is an ordering field: a magnetic field in the case of spins, a chemical potential or difference in chemical potentials in a fluid, etc.) In the case of a magnetic field, the obvious order parameter is the magnetisation per spin, m. The second assumption of Landau s theory is that the free energy of the system can be expanded in powers of the order parameter. From this expansion, all critical properties of the system may be obtained. In the absence of an ordering field, h = 0, where the transition is to take place (external fields are usually coupled to the first power of the order parameter), Landau argued that only even powers of η can be included in the expansion if a symmetry (invariance of F) under the operation η η is assumed. Then we write: G(T; η) = A(T) + B(T)η 2 + C(T)η 4 + (h = 0). (3.130) This expansion should in principle be valid sufficiently close to the critical point, since in that region the order parameter is either zero (in the high-symmetry phase) or very small. For the sake of convenience, we express the temperature in terms of the scaled temperature Now we assume the validity of the expansions t = T T c. (3.131) T A(t) = A 0 + A 1 t + A 2 t 2 +, B(t) = B 0 + B 1 t + B 2 t 2 +, C(t) = C 0 + C 1 t + C 2 t 2 + (3.132) For the time being, we truncate the expansion (3.130) after the fourth power in η. This implies that the parameter C(t) > 0 for all t, so that the free energy is bounded from below (implying stability of the system). The equilibrium value of the order parameter is obtained by minimising the free energy with respect to η at a fixed value of the reduced temperature t: [ 2B(t)η + 4C(t)η 3 = 0 η = 0, η = ± B(t) ] 1/2. (3.133) 2C(t) Clearly, if B(t) > 0, the only possible solution is η = 0, i.e. the high-symmetry phase (stable for T T c ). In order for the low-symmetry phase to be stable, we must have B(t) < 0 for T < T c. The simplest possibility is to have B 0 = 0, C 0 > 0 so that close to the transition we can write G(T; η) = G 0 (t) + B(t)η 2 + C 0 η 4, B(t) = B 1 t, (3.134) 87
3 Phase transitions. 3.1 Introduction. 3.2 Thermodynamics of phase transitions
3 Phase transitions 3.1 Introduction Phase transitions are ubiquitous in Nature. We are all familiar with the different phases of water (vapour, liquid and ice), and with the change from one to the other.
More informationThe Ising model Summary of L12
The Ising model Summary of L2 Aim: Study connections between macroscopic phenomena and the underlying microscopic world for a ferromagnet. How: Study the simplest possible model of a ferromagnet containing
More informationS i J <ij> h mf = h + Jzm (4) and m, the magnetisation per spin, is just the mean value of any given spin. S i = S k k (5) N.
Statistical Physics Section 10: Mean-Field heory of the Ising Model Unfortunately one cannot solve exactly the Ising model or many other interesting models) on a three dimensional lattice. herefore one
More informationLectures 16: Phase Transitions
Lectures 16: Phase Transitions Continuous Phase transitions Aims: Mean-field theory: Order parameter. Order-disorder transitions. Examples: β-brass (CuZn), Ferromagnetic transition in zero field. Universality.
More informationStatistical Physics. Solutions Sheet 11.
Statistical Physics. Solutions Sheet. Exercise. HS 0 Prof. Manfred Sigrist Condensation and crystallization in the lattice gas model. The lattice gas model is obtained by dividing the volume V into microscopic
More informationSpontaneous Symmetry Breaking
Spontaneous Symmetry Breaking Second order phase transitions are generally associated with spontaneous symmetry breaking associated with an appropriate order parameter. Identifying symmetry of the order
More information5.1 Mean-Field Treatments of the Ising Model
Chapter 5 Mean-Field Theory 5.1 Mean-Field Treatments of the Ising Model In this Section, we are discussing various approaches to obtain a mean-field solution to the Ising model. In fact, several of the
More information3. General properties of phase transitions and the Landau theory
3. General properties of phase transitions and the Landau theory In this Section we review the general properties and the terminology used to characterise phase transitions, which you will have already
More informationPhase transitions and critical phenomena
Phase transitions and critical phenomena Classification of phase transitions. Discontinous (st order) transitions Summary week -5 st derivatives of thermodynamic potentials jump discontinously, e.g. (
More informationPhase Transitions and Critical Behavior:
II Phase Transitions and Critical Behavior: A. Phenomenology (ibid., Chapter 10) B. mean field theory (ibid., Chapter 11) C. Failure of MFT D. Phenomenology Again (ibid., Chapter 12) // Windsor Lectures
More informationJ ij S i S j B i S i (1)
LECTURE 18 The Ising Model (References: Kerson Huang, Statistical Mechanics, Wiley and Sons (1963) and Colin Thompson, Mathematical Statistical Mechanics, Princeton Univ. Press (1972)). One of the simplest
More informationClusters and Percolation
Chapter 6 Clusters and Percolation c 2012 by W. Klein, Harvey Gould, and Jan Tobochnik 5 November 2012 6.1 Introduction In this chapter we continue our investigation of nucleation near the spinodal. We
More informationVI. Series Expansions
VI. Series Expansions VI.A Low-temperature expansions Lattice models can also be studied by series expansions. Such expansions start with certain exactly solvable limits, and typically represent perturbations
More informationPhysics 212: Statistical mechanics II Lecture XI
Physics 212: Statistical mechanics II Lecture XI The main result of the last lecture was a calculation of the averaged magnetization in mean-field theory in Fourier space when the spin at the origin is
More informationPH4211 Statistical Mechanics Brian Cowan
PH4211 Statistical Mechanics Brian Cowan Contents 1 The Methodology of Statistical Mechanics 1.1 Terminology and Methodology 1.1.1 Approaches to the subject 1.1.2 Description of states 1.1.3 Extensivity
More informationWORLD SCIENTIFIC (2014)
WORLD SCIENTIFIC (2014) LIST OF PROBLEMS Chapter 1: Magnetism of Free Electrons and Atoms 1. Orbital and spin moments of an electron: Using the theory of angular momentum, calculate the orbital
More informationPhysics 127b: Statistical Mechanics. Second Order Phase Transitions. The Ising Ferromagnet
Physics 127b: Statistical Mechanics Second Order Phase ransitions he Ising Ferromagnet Consider a simple d-dimensional lattice of N classical spins that can point up or down, s i =±1. We suppose there
More informationStatistical Thermodynamics Solution Exercise 8 HS Solution Exercise 8
Statistical Thermodynamics Solution Exercise 8 HS 05 Solution Exercise 8 Problem : Paramagnetism - Brillouin function a According to the equation for the energy of a magnetic dipole in an external magnetic
More informationStatistical mechanics, the Ising model and critical phenomena Lecture Notes. September 26, 2017
Statistical mechanics, the Ising model and critical phenomena Lecture Notes September 26, 2017 1 Contents 1 Scope of these notes 3 2 Partition function and free energy 4 3 Definition of phase transitions
More informationRenormalization Group for the Two-Dimensional Ising Model
Chapter 8 Renormalization Group for the Two-Dimensional Ising Model The two-dimensional (2D) Ising model is arguably the most important in statistical physics. This special status is due to Lars Onsager
More informationGinzburg-Landau Theory of Phase Transitions
Subedi 1 Alaska Subedi Prof. Siopsis Physics 611 Dec 5, 008 Ginzburg-Landau Theory of Phase Transitions 1 Phase Transitions A phase transition is said to happen when a system changes its phase. The physical
More informationPhase transitions and finite-size scaling
Phase transitions and finite-size scaling Critical slowing down and cluster methods. Theory of phase transitions/ RNG Finite-size scaling Detailed treatment: Lectures on Phase Transitions and the Renormalization
More informationPhase Transitions in Condensed Matter Spontaneous Symmetry Breaking and Universality. Hans-Henning Klauss. Institut für Festkörperphysik TU Dresden
Phase Transitions in Condensed Matter Spontaneous Symmetry Breaking and Universality Hans-Henning Klauss Institut für Festkörperphysik TU Dresden 1 References [1] Stephen Blundell, Magnetism in Condensed
More informationThe Equilibrium State
Materials Science & Metallurgy Part III Course M16 Computation of Phase Diagrams (Revision) H. K. D. H. Bhadeshia The Equilibrium State Equilibrium is a state in which no further change is perceptible,
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
1 Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 2, 24 March 2006 1 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationMetropolis Monte Carlo simulation of the Ising Model
Metropolis Monte Carlo simulation of the Ising Model Krishna Shrinivas (CH10B026) Swaroop Ramaswamy (CH10B068) May 10, 2013 Modelling and Simulation of Particulate Processes (CH5012) Introduction The Ising
More informationPhenomenological Theories of Nucleation
Chapter 1 Phenomenological Theories of Nucleation c 2012 by William Klein, Harvey Gould, and Jan Tobochnik 16 September 2012 1.1 Introduction These chapters discuss the problems of nucleation, spinodal
More informationChapter 4 Phase Transitions. 4.1 Phenomenology Basic ideas. Partition function?!?! Thermodynamic limit Statistical Mechanics 1 Week 4
Chapter 4 Phase Transitions 4.1 Phenomenology 4.1.1 Basic ideas Partition function?!?! Thermodynamic limit 4211 Statistical Mechanics 1 Week 4 4.1.2 Phase diagrams p S S+L S+G L S+G L+G G G T p solid triple
More informationc 2007 by Harvey Gould and Jan Tobochnik 28 May 2007
Chapter 5 Magnetic Systems c 2007 by Harvey Gould and Jan Tobochnik 28 May 2007 We apply the general formalism of statistical mechanics developed in Chapter 4 to the Ising model, a model magnetic system
More informationCritical Behavior I: Phenomenology, Universality & Scaling
Critical Behavior I: Phenomenology, Universality & Scaling H. W. Diehl Fachbereich Physik, Universität Duisburg-Essen, Campus Essen 1 Goals recall basic facts about (static equilibrium) critical behavior
More informationPhysics 127b: Statistical Mechanics. Renormalization Group: 1d Ising Model. Perturbation expansion
Physics 17b: Statistical Mechanics Renormalization Group: 1d Ising Model The ReNormalization Group (RNG) gives an understanding of scaling and universality, and provides various approximation schemes to
More informationPhysics 127b: Statistical Mechanics. Landau Theory of Second Order Phase Transitions. Order Parameter
Physics 127b: Statistical Mechanics Landau Theory of Second Order Phase Transitions Order Parameter Second order phase transitions occur when a new state of reduced symmetry develops continuously from
More information8.334: Statistical Mechanics II Spring 2014 Test 2 Review Problems
8.334: Statistical Mechanics II Spring 014 Test Review Problems The test is closed book, but if you wish you may bring a one-sided sheet of formulas. The intent of this sheet is as a reminder of important
More informationThe 1+1-dimensional Ising model
Chapter 4 The 1+1-dimensional Ising model The 1+1-dimensional Ising model is one of the most important models in statistical mechanics. It is an interacting system, and behaves accordingly. Yet for a variety
More informationNumerical Analysis of 2-D Ising Model. Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011
Numerical Analysis of 2-D Ising Model By Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011 Contents Abstract Acknowledgment Introduction Computational techniques Numerical Analysis
More informationThe (magnetic) Helmholtz free energy has proper variables T and B. In differential form. and the entropy and magnetisation are thus given by
4.5 Landau treatment of phase transitions 4.5.1 Landau free energy In order to develop a general theory of phase transitions it is necessary to extend the concept of the free energy. For definiteness we
More informationCritical Exponents. From P. Chaikin and T Lubensky Principles of Condensed Matter Physics
Critical Exponents From P. Chaikin and T Lubensky Principles of Condensed Matter Physics Notice that convention allows for different exponents on either side of the transition, but often these are found
More informationIntroduction to Phase Transitions in Statistical Physics and Field Theory
Introduction to Phase Transitions in Statistical Physics and Field Theory Motivation Basic Concepts and Facts about Phase Transitions: Phase Transitions in Fluids and Magnets Thermodynamics and Statistical
More informationThe Phase Transition of the 2D-Ising Model
The Phase Transition of the 2D-Ising Model Lilian Witthauer and Manuel Dieterle Summer Term 2007 Contents 1 2D-Ising Model 2 1.1 Calculation of the Physical Quantities............... 2 2 Location of the
More informationPHASE TRANSITIONS IN SOFT MATTER SYSTEMS
OUTLINE: Topic D. PHASE TRANSITIONS IN SOFT MATTER SYSTEMS Definition of a phase Classification of phase transitions Thermodynamics of mixing (gases, polymers, etc.) Mean-field approaches in the spirit
More informationPhase Transition & Approximate Partition Function In Ising Model and Percolation In Two Dimension: Specifically For Square Lattices
IOSR Journal of Applied Physics (IOSR-JAP) ISS: 2278-4861. Volume 2, Issue 3 (ov. - Dec. 2012), PP 31-37 Phase Transition & Approximate Partition Function In Ising Model and Percolation In Two Dimension:
More information4 Critical phenomena. 4.1 Ginzburg Landau theory of fluctuations
4 Critical phenomena 4.1 Ginzburg Landau theory of fluctuations We know that, in the critical region, fluctuations are very important but, as we have realised from the previous discussion, what is missing
More informationChapter 6. Phase transitions. 6.1 Concept of phase
Chapter 6 hase transitions 6.1 Concept of phase hases are states of matter characterized by distinct macroscopic properties. ypical phases we will discuss in this chapter are liquid, solid and gas. Other
More informationPhase transition and spontaneous symmetry breaking
Phys60.nb 111 8 Phase transition and spontaneous symmetry breaking 8.1. Questions: Q1: Symmetry: if a the Hamiltonian of a system has certain symmetry, can the system have a lower symmetry? Q: Analyticity:
More informationS j H o = gµ o H o. j=1
LECTURE 17 Ferromagnetism (Refs.: Sections 10.6-10.7 of Reif; Book by J. S. Smart, Effective Field Theories of Magnetism) Consider a solid consisting of N identical atoms arranged in a regular lattice.
More informationKrammers-Wannier Duality in Lattice Systems
Krammers-Wannier Duality in Lattice Systems Sreekar Voleti 1 1 Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7. (Dated: December 9, 2018) I. INTRODUCTION It was shown by
More informationSuggestions for Further Reading
Contents Preface viii 1 From Microscopic to Macroscopic Behavior 1 1.1 Introduction........................................ 1 1.2 Some Qualitative Observations............................. 2 1.3 Doing
More informationImperfect (non-ideal) Gases
Chapter 2 Imperfect non-ideal) Gases Before getting to a completely quantum treatment of particles, however, let s step back for a moment and consider real gases. Virtually all textbooks on statistical
More informationMonte Carlo study of the Baxter-Wu model
Monte Carlo study of the Baxter-Wu model Nir Schreiber and Dr. Joan Adler Monte Carlo study of the Baxter-Wu model p.1/40 Outline Theory of phase transitions, Monte Carlo simulations and finite size scaling
More informationLinear Theory of Evolution to an Unstable State
Chapter 2 Linear Theory of Evolution to an Unstable State c 2012 by William Klein, Harvey Gould, and Jan Tobochnik 1 October 2012 2.1 Introduction The simple theory of nucleation that we introduced in
More informationSimple lattice-gas model for water
Simple lattice-gas model for water A. Ciach, W. T. Góźdź and A. Perera Motivation H-bonds and packing of molecules Construction of the lattice model Results in mean-field approximation Water molecule and
More informationPhase Transitions in Multicomponent Systems
Physics 127b: Statistical Mechanics Phase Transitions in Multicomponent Systems The ibbs Phase Rule Consider a system with n components (different types of molecules) with r s in equilibrium. The state
More informationChapter 2 Ensemble Theory in Statistical Physics: Free Energy Potential
Chapter Ensemble Theory in Statistical Physics: Free Energy Potential Abstract In this chapter, we discuss the basic formalism of statistical physics Also, we consider in detail the concept of the free
More informationPrinciples of Equilibrium Statistical Mechanics
Debashish Chowdhury, Dietrich Stauffer Principles of Equilibrium Statistical Mechanics WILEY-VCH Weinheim New York Chichester Brisbane Singapore Toronto Table of Contents Part I: THERMOSTATICS 1 1 BASIC
More information(1) Consider the ferromagnetic XY model, with
PHYSICS 10A : STATISTICAL PHYSICS HW ASSIGNMENT #7 (1 Consider the ferromagnetic XY model, with Ĥ = i
More information8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization
8.334: Statistical Mechanics II Problem Set # 4 Due: 4/9/14 Transfer Matrices & Position space renormalization This problem set is partly intended to introduce the transfer matrix method, which is used
More informationUNDERSTANDING BOLTZMANN S ANALYSIS VIA. Contents SOLVABLE MODELS
UNDERSTANDING BOLTZMANN S ANALYSIS VIA Contents SOLVABLE MODELS 1 Kac ring model 2 1.1 Microstates............................ 3 1.2 Macrostates............................ 6 1.3 Boltzmann s entropy.......................
More informationLecture 11: Long-wavelength expansion in the Neel state Energetic terms
Lecture 11: Long-wavelength expansion in the Neel state Energetic terms In the last class we derived the low energy effective Hamiltonian for a Mott insulator. This derivation is an example of the kind
More informationIntroduction to Disordered Systems and Spin Glasses
Introduction to Disordered Systems and Spin Glasses Gautam I. Menon and Purusattam Ray Abstract This chapter reviews some basic material required for the study of disordered systems, establishing some
More informationSlightly off-equilibrium dynamics
Slightly off-equilibrium dynamics Giorgio Parisi Many progresses have recently done in understanding system who are slightly off-equilibrium because their approach to equilibrium is quite slow. In this
More informationThe glass transition as a spin glass problem
The glass transition as a spin glass problem Mike Moore School of Physics and Astronomy, University of Manchester UBC 2007 Co-Authors: Joonhyun Yeo, Konkuk University Marco Tarzia, Saclay Mike Moore (Manchester)
More informationChemistry C : Polymers Section. Dr. Edie Sevick, Research School of Chemistry, ANU. 5.0 Thermodynamics of Polymer Solutions
Chemistry C3102-2006: Polymers Section Dr. Edie Sevick, Research School of Chemistry, AU 5.0 Thermodynamics of Polymer Solutions In this section, we investigate the solubility of polymers in small molecule
More informationMonte Carlo Simulation of the Ising Model. Abstract
Monte Carlo Simulation of the Ising Model Saryu Jindal 1 1 Department of Chemical Engineering and Material Sciences, University of California, Davis, CA 95616 (Dated: June 9, 2007) Abstract This paper
More informationA Monte Carlo Implementation of the Ising Model in Python
A Monte Carlo Implementation of the Ising Model in Python Alexey Khorev alexey.s.khorev@gmail.com 2017.08.29 Contents 1 Theory 1 1.1 Introduction...................................... 1 1.2 Model.........................................
More information1. Scaling. First carry out a scale transformation, where we zoom out such that all distances appear smaller by a factor b > 1: Δx Δx = Δx/b (307)
45 XVIII. RENORMALIZATION GROUP TECHNIQUE The alternative technique to man field has to be significantly different, and rely on a new idea. This idea is the concept of scaling. Ben Widom, Leo Kadanoff,
More informationUniversità degli Studi di Ferrara. Critical Properties of the Potts Glass with many states
Università degli Studi di Ferrara Dottorato di Ricerca in Fisica Ciclo XXIII Coordinatore Prof. F. Frontera Critical Properties of the Potts Glass with many states Settore Scientifico Disciplinare FIS/03
More informationExamination paper for TFY4245 Faststoff-fysikk, videregående kurs
Side 1 av 6 Department of Physics xamination paper for TFY445 Faststoff-fysikk, videregående kurs Academic contact during examination: Ragnvald Mathiesen Phone: 976913 xamination date: 0.06.015 xamination
More informationMagnetism at finite temperature: molecular field, phase transitions
Magnetism at finite temperature: molecular field, phase transitions -The Heisenberg model in molecular field approximation: ferro, antiferromagnetism. Ordering temperature; thermodynamics - Mean field
More informationPhases of matter and phase diagrams
Phases of matter and phase diagrams Transition to Supercritical CO2 Water Ice Vapor Pressure and Boiling Point Liquids boil when the external pressure equals the vapor pressure. Temperature of boiling
More informationScaling Theory. Roger Herrigel Advisor: Helmut Katzgraber
Scaling Theory Roger Herrigel Advisor: Helmut Katzgraber 7.4.2007 Outline The scaling hypothesis Critical exponents The scaling hypothesis Derivation of the scaling relations Heuristic explanation Kadanoff
More informationSpontaneous magnetization of the square 2D Ising lattice with nearest- and weak next-nearest-neighbour interactions
Phase Transitions Vol. 82, No. 2, February 2009, 191 196 Spontaneous magnetization of the square 2D Ising lattice with nearest- and weak next-nearest-neighbour interactions H.J.W. Zandvliet a * and C.
More informationStatistical Mechanics
Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2
More informationIsing model and phase transitions
Chapter 5 Ising model and phase transitions 05 by Alessandro Codello 5. Equilibrium statistical mechanics Aspinsystemisdescribedbyplacingaspinvariableσ i {, } at every site i of a given lattice. A microstate
More information1 Phase Transitions. 1.1 Introductory Phenomenology
1 Phase Transitions 1.1 Introductory Phenomenology Consider a single-component system in thermodynamic equilibrium. Let us describe the system with the set of (independent) variables T, p, and N. The appropriate
More informationNATURAL SCIENCES TRIPOS. Past questions. EXPERIMENTAL AND THEORETICAL PHYSICS Minor Topics. (27 February 2010)
NATURAL SCIENCES TRIPOS Part III Past questions EXPERIMENTAL AND THEORETICAL PHYSICS Minor Topics (27 February 21) 1 In one-dimension, the q-state Potts model is defined by the lattice Hamiltonian βh =
More informationPhysics 115/242 Monte Carlo simulations in Statistical Physics
Physics 115/242 Monte Carlo simulations in Statistical Physics Peter Young (Dated: May 12, 2007) For additional information on the statistical Physics part of this handout, the first two sections, I strongly
More informationFerromagnetism. Iron, nickel, and cobalt are ferromagnetic.
Ferromagnetism Technische Universität Graz Institute of Solid State Physics Ferromagnetism elow a critical temperature (called the Curie temperature) a magnetization spontaneously appears in a ferromagnet
More informationEvaporation/Condensation of Ising Droplets
, Elmar Bittner and Wolfhard Janke Institut für Theoretische Physik, Universität Leipzig, Augustusplatz 10/11, D-04109 Leipzig, Germany E-mail: andreas.nussbaumer@itp.uni-leipzig.de Recently Biskup et
More informationREVIEW: Derivation of the Mean Field Result
Lecture 18: Mean Field and Metropolis Ising Model Solutions 1 REVIEW: Derivation of the Mean Field Result The Critical Temperature Dependence The Mean Field Approximation assumes that there will be an
More information16 Singular perturbations
18.354J Nonlinear Dynamics II: Continuum Systems Lecture 1 6 Spring 2015 16 Singular perturbations The singular perturbation is the bogeyman of applied mathematics. The fundamental problem is to ask: when
More informationV. SCALING. A. Length scale - correlation length. h f 1 L a. (50)
V. SCALING A. Length scale - correlation length This broken symmetry business gives rise to a concept pretty much missed until late into the 20th century: length dependence. How big does the system need
More informationFinite-size analysis via the critical energy -subspace method in the Ising models
Materials Science-Poland, Vol. 23, No. 4, 25 Finite-size analysis via the critical energy -subspace method in the Ising models A. C. MAAKIS *, I. A. HADJIAGAPIOU, S. S. MARTINOS, N. G. FYTAS Faculty of
More informationNotes on Spin Operators and the Heisenberg Model. Physics : Winter, David G. Stroud
Notes on Spin Operators and the Heisenberg Model Physics 880.06: Winter, 003-4 David G. Stroud In these notes I give a brief discussion of spin-1/ operators and their use in the Heisenberg model. 1. Spin
More informationElectromagnetism II. Instructor: Andrei Sirenko Spring 2013 Thursdays 1 pm 4 pm. Spring 2013, NJIT 1
Electromagnetism II Instructor: Andrei Sirenko sirenko@njit.edu Spring 013 Thursdays 1 pm 4 pm Spring 013, NJIT 1 PROBLEMS for CH. 6 http://web.njit.edu/~sirenko/phys433/phys433eandm013.htm Can obtain
More informationMagnets, 1D quantum system, and quantum Phase transitions
134 Phys620.nb 10 Magnets, 1D quantum system, and quantum Phase transitions In 1D, fermions can be mapped into bosons, and vice versa. 10.1. magnetization and frustrated magnets (in any dimensions) Consider
More informationComplex Systems Methods 9. Critical Phenomena: The Renormalization Group
Complex Systems Methods 9. Critical Phenomena: The Renormalization Group Eckehard Olbrich e.olbrich@gmx.de http://personal-homepages.mis.mpg.de/olbrich/complex systems.html Potsdam WS 2007/08 Olbrich (Leipzig)
More informationarxiv:quant-ph/ v2 24 Dec 2003
Quantum Entanglement in Heisenberg Antiferromagnets V. Subrahmanyam Department of Physics, Indian Institute of Technology, Kanpur, India. arxiv:quant-ph/0309004 v2 24 Dec 2003 Entanglement sharing among
More informationEquilibrium, out of equilibrium and consequences
Equilibrium, out of equilibrium and consequences Katarzyna Sznajd-Weron Institute of Physics Wroc law University of Technology, Poland SF-MTPT Katarzyna Sznajd-Weron (WUT) Equilibrium, out of equilibrium
More informationExamination paper for TFY4245 Faststoff-fysikk, videregående kurs
Side 1 av 7 Department of Physics Examination paper for TFY445 Faststoff-fysikk, videregående kurs Academic contact during examination: Ragnvald Mathiesen Phone: 976913 Examination date: 08.06.017 Examination
More informationRecitation: 10 11/06/03
Recitation: 10 11/06/03 Ensembles and Relation to T.D. It is possible to expand both sides of the equation with F = kt lnq Q = e βe i If we expand both sides of this equation, we apparently obtain: i F
More informationwhere we have recognised A(T) to be the free energy of the high-symmetry phase, G 0 (t). The order parameter is then t 1/2, t 0, (3.
where we have recognised A(T) to be the free energy of the high-symmetry phase, G 0 (t). The order parameter is then [ η = ± B ] 1/2 1t t 1/2, t 0, (3.1) 2C 0 giving β = 1/2, as in mean-field theory. In
More informationMetropolis, 2D Ising model
Metropolis, 2D Ising model You can get visual understanding from the java applets available, like: http://physics.ucsc.edu/~peter/ising/ising.html Average value of spin is magnetization. Abs of this as
More informationThe Second Virial Coefficient & van der Waals Equation
V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.v.33). As discussed before, the two-body potential is characterized by
More informationLuigi Paolasini
Luigi Paolasini paolasini@esrf.fr LECTURE 7: Magnetic excitations - Phase transitions and the Landau mean-field theory. - Heisenberg and Ising models. - Magnetic excitations. External parameter, as for
More informationf(t,h) = t 2 g f (h/t ), (3.2)
Chapter 3 The Scaling Hypothesis Previously, we found that singular behaviour in the vicinity of a second order critical point was characterised by a set of critical exponents {α,β,γ,δ, }. These power
More informationQuantum Phase Transitions
Quantum Phase Transitions Subir Sachdev Department of Physics Yale University P.O. Box 208120, New Haven, CT 06520-8120 USA E-mail: subir.sachdev@yale.edu May 19, 2004 To appear in Encyclopedia of Mathematical
More informationPhase Transitions of Random Binary Magnetic Square Lattice Ising Systems
I. Q. Sikakana Department of Physics and Non-Destructive Testing, Vaal University of Technology, Vanderbijlpark, 1900, South Africa e-mail: ike@vut.ac.za Abstract Binary magnetic square lattice Ising system
More informationfiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES
Content-Thermodynamics & Statistical Mechanics 1. Kinetic theory of gases..(1-13) 1.1 Basic assumption of kinetic theory 1.1.1 Pressure exerted by a gas 1.2 Gas Law for Ideal gases: 1.2.1 Boyle s Law 1.2.2
More informationRenormalization Group analysis of 2D Ising model
Renormalization Group analysis of D Ising model Amir Bar January 7, 013 1 Introduction In this tutorial we will see explicitly how RG can be used to probe the phase diagram of d > 1 systems, focusing as
More informationExsolution and Ordering in Solid Solutions
Exsolution and Ordering in Solid Solutions 1.1 Introduction to solid solutions A solid solution is a single phase which exists over a range of chemical compositions. Some materials are able to tolerate
More information