The "Ising model" is central to many problems and systems in statistical mechanics. It first came to the fore in the study of ferromagnetic systems. It was found that the use of such **simplified models** paved the way for greater understanding via modeling of more complex systems. Let's now look at an elementary, statistical mechanical spin system. A fairly mundane example is a two-dimensional Ising model for ferromagnetic matter. It contains magnetic domains for which the individual spin magnets can be subject to sudden reversals. For a simple example, think of the 2D model of 4 x 4 elementary spin magnets as shown below:

The degree of order, as well as information, is determined from what is called "

**", or the**

*the spin excess***net spin difference**(up minus down or vice versa). The larger this number, the greater the degree of order, and the

**lower the entropy**of the system. Obviously, since 0 denotes an extremely low number, we can deduce large entropy.

Consider the system S(2) in more detail, noting the right side orientations of the elementary spin magnets. Here we get: 14 spin ups - 2 spin downs = 12 spin ups, or in other words the spin excess = 12. This system, S(2), has much higher degree of order (less entropy) than the system S(1). (We should add here that

*higher entropy*- as in S(1) - corresponds to

*the most probable state*, defined by the minimal spin excess of zero.)

Accessing such simple systems allows us to infer fundamental measures applicable to the systems, for example the

**"magnetic moment**" of a state, as well as the "

*degeneracy function*". Consider an N= 2 model system with either 2 ups (two up arrows) or 2 downs.

Then, if m denotes the magnetic permeability we can have:

M = + 2 m or M = -2 m

where the first is the magnetic moment for two spin- up particles, and the second for two spins down. One can also, of course, have the mixed state inclusive of one spin up plus one spin down, then:

M = O m or O

Meantime, the degeneracy function computes the number of states having the same value of m (individual spins) or M such that:

g(N,m) = N!/ (½N + m)! (½N - m)!

_{av}]

where [M

_{av}]

_{s})

__Bulletin of the American Mathematical Society__, Vol. 53, No. 2, p. 187). In this paper the author presents us with the 2D ice crystal Ising model shown below:

The most basic probability for any such system is "thermal equilibrium". Thus, at some temperature T if the system attains thermal equilibrium then the probability of any particular configuration decays exponentially with the energy of C, which is analogous to E =

_{m}B in the ferromagnetic case. The probability of any particular configuration dependent on T is then:

Where k is Boltzmann's constant, 1.38 x 10

^{-23}J

**/**K.

One will also make use of the "partition function":

Z(T) =

_{C}exp (- Energy (C)/ kT)

which as Okounkov notes, really functions as a "normalization factor" given that it "makes the probabilities sum to 1". In this Ising ice crystal model, then, the energy is "the sum of interactions of all adjacent squares." Since the total number of squares is fixed (see stipulation (1)) then the energy must be proportional to the total length of the contours separating white from blue.

To identify the contours is easy. If the energetic reader will run off a copy of the image of the 2D rectangle, then take a black magic marker and trace around each ice crystal region as it appears, he will have generated the contours. The normalization for energy is then (

*op. cit*.):

As in the case of the ferromagnetic system entropy competes with order (energy). In the Ising ice crystal energy is saved via clumping. If we designate an "order parameter" such that

**then in the ice crystal case the larger b the stronger tendency for order. Interestingly, as Okounkov notes there is a critical temperature T**

^{-1}**> 0 above which entropy wins, it is:**

_{c}b = 0.5 ln (

**Ö**2 + 1)

Below T

**and for ice crystal concentrations above a certain threshold a crystal will form as the size of the container goes to infinity.**

_{c}**Suggested Problems:**

1) Quantify the magnetic energy for the system at time t(o) compared to time t(o) + t, if the magnetic energy of one spin magnet can be written:

M = - m

**B**cos Θ

where m is the magnetic moment (-eL/2m, L = 1) and assume Θ = +/- π, and

**B**= 0.1T.

2) Say that S = log (g) determines the entropy for a simple statistical mechanical system, where g denotes the number of accessible states. Then estimate S for the 2D ice crystal model - including any errors that might enter.

## No comments:

Post a Comment