## Abstract

Though the idea to use numerical techniques, in order to solve complex three-dimensional problems, has become quite old, computational techniques have gained immense importance in past few decades because of the advent of new generation fast and efficient computers and development of algorithms as parallel computing. Many mathematical problems have no exact solutions. Depending on the complexity of the equations, one needs to use approximate methods. But there are problems, which are beyond our limits, and need support of computers. Ernst Ising published his PhD dissertation in the form of a scientific report in 1925. He used a string of magnetic moments; spin up (+1/2) and spin down (−1/2), and applied periodic boundary conditions to prove that magnetic phase transition does not exist in one dimensions. Lars Onsager, latter, exactly solved the phase transition problem in two dimensions in 1944. It is going to be a century-old problem now. A variety of potential applications of Ising model are possible now a days; classified as Ising universality class models. It has now become possible to solve phase transition problems in complex three-dimensional geometries. Though the area of spinotronics still needs more engagements of computational techniques, its limited use have provided good insights at molecular scale in recent past. This chapter gives a brief introduction to Ising model and its applications, highlighting the developments in the field of magnetism relevant to the area of solid state physics.

### Keywords

- surface-directed phase separation
- wetting-dewetting
- Monte Carlo simulation

## 1. Introduction

Ernst Ising (Figure 1) was born on May 10, 1900, in Loe Koln. He started schooling in 1907 and obtained his diploma at the gymnasium there in the year 1918. After brief military training, he studied mathematics and physics at Gottingen University in the year 1919. After a short gap, he continued his studies and learnt astronomy apart of other subjects. He got focused to theoretical physics at the suggestion of Professor W. Lenz. He started investigating ferromagnetism under supervision of W. Lenz by the end of the year 1922. Ising published short paper in 1925 as a summary of his doctoral thesis [1, 2]. He exactly calculated partition function for one-dimensional lattice system of spins. Ising had first proven that no phase transition to a ferromagnetic ordered state occurs in one dimension at any temperature.

His argument in the favor of his mathematical note was very simple. Suppose, if one of the spins get flipped at a random position because of thermal agitation, there is no force available, which can stop the neighboring spins to flip in the same direction. And this process will go on and on, and completely ordered state will not remain stable at a finite temperature. Thus no phase transition will occur at a finite temperature. Ideally speaking, any ordered state will always remain like a metastable state at finite temperature and nothing more. Molecular motion seizes at absolute zero temperature. So, one may expect that no spin fluctuations may occur at absolute zero temperature. Henceforth, the stable ordered state is a natural outcome at absolute zero temperature. But, it cannot be said to be a critical temperature in true sense. The existence of phase transition at this temperature has no physical meaning, because there is no temperature below it. After going through some approximate calculations, Ising purportedly showed that his model could not exhibit a phase transition in two and three dimensions, either. Latter, his conclusion was proven to be erroneous [1, 2] (Figure 2).

Barry Simon has quoted it very well “This model was suggested to Ising by his thesis advisor, Lenz. Ising solved the one-dimensional model, and on the basis of the fact that the one-dimensional model had no phase transition; he asserted that there was no phase transition in any dimension. As we shall see, this is false. It is ironic that on the basis of an elementary calculation and erroneous conclusion, Ising’s name has become among the most commonly mentioned in the theoretical physics literature. But history has had its revenge. Ising’s name, which is correctly pronounced “E-zing”, is almost universally mispronounced “I-zing”.”

Ising’s paper credited Wilhelm Lenz for his original idea, who had first proposed it in the year 1920. W. Lenz was Ising’s research supervisor. It has been often rendered as Lenz-Ising model in many citations. Lenz suggested that dipolar atoms in crystals are free to rotate in quantized manner. He proposed quantum treatment of dipole orientations, though in its classical version, Ising considered only two spin states, i.e., S = ±½. Ising discussed his results with Professor Lenz and Dr. Wolfgang Pauli, who was teaching at Hamburg at that time. Ising’s work was first cited by famous contemporary scientist Heisenberg. Heisenberg was first one to realize the failure of Lenz-Ising model. In order to explain ferromagnetism, he developed his own theory, using complicated interactions of spins. There are more scientists in the list, whose contribution to Lenz-Ising model or simply say Ising model must be cited here, because of their historical relevance. They have greatly enriched and contributed to this new model. This list includes scientists like Gorskly (1928), R. H. Fowler (1930), Bragg and Williams (1934), R. Peierls (1936), J. G. Krikwood (1938), Hens Bethe (1939), Kramers and Wannier (1941), and Onsager (1942). They further extended Ising model to a new class of problems.

## 2. Application of Ising model

Ising model has been extensively used for solving a variety of problems [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. Some of the problems are discussed, here, with appropriate examples.

### 2.1 Phase separation and wetting/dewetting

Ising model was first exploited for investigating spontaneous magnetization in ferromagnetic film (i.e. magnetization in the absence of external magnetic field). An example case of Ising model using metropolis algorithm is shown in Figure 3. Transition temperature depends on the strength of the inter-spin exchange coupling; the dominating term governs the kinetics, when long-range interactions are introduced in the calculations. Latter, it was used to study phase separation in binary alloys and liquid-gas phase transitions (i.e., condensation of molecule in one region of space of the box). Binary alloys constitute of two different atoms. At temperature T = 0, Zn-Cu alloy; known as brass, gets completely ordered. This state is said to be β-brass. In β-brass state, each Zn atom is surrounded by eight copper atoms, placed at the corners of the unit cell of the body-centered cubic structure and vice versa. The occupation of each site can be represented by:

The interaction energy between A-A, B-B, and A-B type of atoms are represented by ε_{AA}, ε_{BB}, and ε_{AB}, respectively. Phase separation has been studied vastly, using Ising model [4, 5, 6]. A phase is simply a part of a system, separated from the other part by the formation of an interface; that essentially means that two components aggregate and form rich regions of A and B type of molecules with an interface in between them. The evolution of two distinct phases, when an initial random but homogeneous mixture is annealed below a definite temperature, is known as phase separation. Phase separation leads to discontinuity and inhomogeneity in the systems. This happens because the phase-separated regions are energetically more stable. Phase separation has been an old problem and has been extended to study diverse phenomena ranging from magnetic liquid-liquid phase separation to protein-protein phase separation in biological systems. This process has also been studied in the presence of external surfaces having affinity to one type of atom or molecule (Figure 4a). Both theoretical and experimental methods have been exploited and have been found in close agreement. Formation of long ridges and circular drops has been reported numerous occasions using lattice-based Ising model. For example, one may look into John W. Cahn research paper published in *The Journal of Chemical Physics* in the year 1965. The TEM image taken for Vycor, in which one phase had been leached out and the voids were filled with lead (Figure 4b).

### 2.2 Lattice-based liquid-gas model

Yang and Lee first coined the term lattice gas in the year 1952. A lattice should have larger volume (V) than the number of lattice molecules (N), so that some of the nodes or lattice vertices are left empty (i.e., N < V). No lattice vertex can be occupied by more than one particle. The interaction potential between two atoms at lattice sites i and j is given by Eq. (2):

For surface affinity of lower surface to i^{th} liquid molecule, we chose:

The occupation number (n_{i}) of a lattice site i is given by:

One example case is shown in Figure 5. Here, we chose lattice size of 128 × 128 × 48. The fluid-fluid molecule and wall-liquid molecule interactions are defined, respectively, in Eqs. (2) and (3). In canonical ensemble, the three-dimensional lattice is swept one by one; by choosing sites regularly with one of its nearest neighbor (i.e., i = n and i + 1 = n + 1). Change in energy is calculated during exchanging of these two sites; the exchange move is accepted, if Exp [−ΔE/k_{B}T] is found to be greater than or equal to a random number generated between [0, 1]. For all cases of studies here, ε = 1 and J_{0} = 12.0, and only the lower surface is functional, while the upper surface has only hard-sphere interaction with the fluid molecules. Average number density for liquid-like molecules is taken as 0.25 [16].

Figure 6 shows micrograph of self-aligned liquid columns. The system evolves from an initial homogeneous mixture of liquid- and gas-like molecules obtained by annealing the system at high temperature for few thousand MC cycles. Dynamic Monte Carlo simulation has been used with continuous but random trial movements of the molecules. The lattice-based Ising model using Eqs. (2) and (3) is also supposed to give same results, at least qualitatively.

### 2.3 Spin glasses

Crystalline solids possess short- and long-range order along its crystal axes and maintain its periodicity in three dimensions. Liquids possess only short-range order, and its molecules have no long-range correlation. Liquid molecules retain only short-range order. Gases possess neither of the two. These are the three phases, in which any matter may exist. What are the glasses then? Glasses are solids, possessing no long-range order. Molecules may only locally arrange themselves to minimize its free energy. If the molecular arrangement is completely random, then a term “random media” is assigned to that. Glasses are understood as supercooled liquids. If a liquid is frozen abruptly, so that the molecules do not get sufficient time to organize themselves, some local order can be retained inside the frozen liquid. Glasses have one peculiar property. These retain relatively higher entropy even at quite low temperatures. One example is Mn doped in metals as impurity. Mn atoms interact with other Mn so (i.e. impurity atom) via RKKY interaction

Lenz-Ising model did not remain limited to above problems only, but it was extensively used to study liquid mixtures, ternary and quaternary alloys, polymer and their mixtures, random walk problem, and many others. The important aspect of Ising model is that a variety of problems (including some problems mentioned above) can be investigated by the similar kind of modeling and approach all together. It is no longer necessary to develop a different kind of theory for each type of cooperative phenomenon. Despite of all the above, it has been ironical that the inventor of the model, Ernst Ising, gave up the idea on working it, any further presuming that his model has no physical significance. He realized after two decades that he had become famous for his model because of the results obtained by other scientist based on his model, rather by his own work. It has been a queer sensation that the results of Ising model matched with any experimental data or the model was bit artificial. As for as the exponents were concerned, they were of universal nature, and a wide variety of systems have the same Ising exponents. The experimental evidence in favor of it remained a challenge, for many decades. In the year 1974, an alloy was found, which first showed that its magnetic behavior exactly matched with the Onsager result.

## 3. Mathematical formulation in one dimension

Various textbooks are available nowadays, which discuss Ising model and its applications in greater details [19, 20, 21, 22]. Here, brief theory of one-dimensional Ising model is presented. H, Q, and A stands for Hamiltonian, partition function, and free energy of the system, respectively:

Some thermodynamic functions are defined as follows:

### 3.1 Case A: free boundary with zero field

Partition function is given by:

Here, K = β J.

We now define a new variable:

Then we can assign σ to two values, i.e., ±1:

In order to consider contributions from all possible configurations {S_{1}, S_{2}, S_{3}………S_{N}}, we need to provide the set of numbers {σ_{1}, σ_{2}, σ_{3}……..σ_{N-1}}; here each S_{i} can take two values as ±1. Configuration in a lattice description means a particular set of values of all spins; if there are N numbers of vertices, there will be 2^{N} different configurations as a result of permutation and combination of spins. The space, thus formed with these configurations, is called configuration space. Here, summing over σ_{i} will give only half value of Q, henceforth, we can write:

### 3.2 Case B: periodic boundary with zero field

Now, the partition function is given by:

Here, S_{N+1} = S_{1}

Since (S_{i})^{2} = 1, we can write S_{1}S_{N}=S_{1}. S_{2}. S_{2}. S_{3}. S_{3}………... S_{N-1}. S_{N-1}. S_{N}

Here, second part in exponential has been converted into a summation series:

It can be shown that in thermodynamic limit, (i.e. N → ∞), the free energy of the system converge to a finite value. Readers are left with the exercise. So, periodic boundary condition, as shown in Figure 7 (invented by Ising), really helps one to get rid of constructing infinitely large systems. Using appropriate boundary conditions, one may obtain realistic results using large but finite number of spins.

## 4. Critical phenomena

A lot of research work has been dedicated to observe system behavior near critical points [23, 24, 25, 26, 27]. The relevant thermodynamic variables exhibit power-law dependences on the parameter (T − T_{c}) specifying the distance away from the critical point. The critical points are marked by the fact that different physical quantities pertaining to the system pose singularities at the critical point. These singularities are expressed in terms of power laws of (T − T_{c}) characterized by critical exponents. As, for example, magnetization <M > identified as an order parameter in magnetism, shows dependence on critical temperature (T_{c}), with exponent β as follows other exponents are also listed below.

Reduced temperature *t* ≡ (*T* − *T* _{c })/*T* _{c }.

α: specific heat c (*t*) ∼ *t* ^{−α}; B ≡ h = 0.

β: spontaneous magnetization M (t) ∼ (−*t*)^{β}, *T* ≤ *T* _{c }, B ≡ h = 0.

γ: magnetic susceptibility χ = ∂M/∂ℎ, *T* ∼ |*t*|−^{γ}, B ≡ h = 0.

δ: critical Isotherm M (h) ∼ |ℎ|^{1/δ} sgn (ℎ), *t* = 0.

ν: correlation length, ξ ∼ |*t*|−^{ν}, B ≡ h = 0.

η: correlation function G (*r*) ∼ *r* ^{(−d+2−η)}, *t* = 0, B ≡ h = 0.

### 4.1 Scaling hypothesis and renormalization group theory

Kadanoff first suggested that, when a system is near critical temperature, individual spins may be grouped into blocks of spins [23]. It is possible because of the fact that the spin-spin correlation length becomes exceedingly large near T_{c} and details of individual spins no longer remain important. In transformed system, each block plays the role of a single spin. Now, the spin variable associated with a single block is denoted by symbol σ_{i}. σ_{i} can take values ±1. The new system is composed of N′ spins (Figure 8).

Lattice constant:

In order to preserve the spatial density of the degrees of freedom of spins in the system, the spatial distances are rescaled by the factor l.

Now, the partition function can be updated as follows:

This idea was first propounded by Kadanoff, and was later developed by Wilson. This process is also referred as decimation process. A new exchange coupling constant is assigned for interaction between σ_{i.} This new construction of lattice does not alter the free energy of the system, and it remains the same as obtained by the original method. The rescaling process helps to find relations between various exponents. More detailed discussion on this topic can be found in standard textbooks of *Statistical Mechanics* by Patharia, Huang, etc. Since this process involves length transformation or a change of scale, Wilson introduced the concept of renormalization group theory after removing certain deficiencies in Kadanoff’s scaling hypothesis. A greater detail of this is omitted here, because that is beyond the scope of the chapter.

## 5. Physical realization: simulation results based on Ising model

We now discuss some of the simulation results obtained using Ising model. Figure 9 shows spontaneous magnetization for a simple cubic crystal (i.e., scc lattice). As the strength of exchange coupling between spin-up and spin-down (J_{AB}) decreases, the critical temperature lowers down. Lower values of J_{AB} weaken the spin flip-flop mechanism; henceforth the system requires further cooling, so that the spin-spin correlation overcomes the fluctuations. Spontaneous magnetization occurs in the absence of external magnetic field [28]. The confirmation of spontaneous process is further confirmed in Figure 10. Figure 10 is plotted for spin correlation function vs. temperature of the system [28]. The critical temperature is marked by the presence of discontinuity in it. Above critical temperature, the magnetization abruptly falls to zero, which is an indication of paramagnetic state. The critical temperature in ferromagnetic thin film is known as Curie temperature. We observe similar kind of behavior with antiferromagnetic films, though below critical point (also known as Neel temperature), the net average magnetization becomes zero, because opposite spins are energetically favored in this case. The schematic diagram is shown in Figure 11 [28]. Magnetization vs. external magnetic field curves are plotted in Figure 12(a)–(d) for different sets of parameters [28].

Simulation results obtained for a magnetically striped system as schematically shown in Figure 13 are reported in Figures 14**–**17 [29]. One or two alternate rectangular regions are created, using external field. Figure 14 shows the gradual transition at the interface, where a definite value of external field suddenly gets zero. The spin polarizations in two regions show sharp boundary. The magnetized film, in presence of magnetic field, induces the magnetic zones in proximity where its close external field is zero. Micrograph also indicates for spin-spin phase separation. The corresponding average magnetization vs. temperature and spin correlation function vs. temperature are also plotted in Figures 15 and 16, respectively, but these studies are done using Monte Carlo simulation with semi-infinite free boundary conditions. It has been observed that these systems have relatively high critical transition temperatures. Figure 17 shows the magnetization process with two alternate magnetized zones [29].

Low-dimensional magnetic heterostructures play vital role in spinotronics. Ferromagnets can induce magnetic ordering through a 40-nm-thick amorphous paramagnetic layer, when placed in its close proximity. One has to reconcile with long-range magnetic interaction to correctly measure the extent of induced magnetization. Readers may go through the *Nature Communications* article of F. Magnus et al. published in the year 2016 [17]. The magnetic properties of ferromagnetic materials with reduced dimensions get altered; when the thickness of a film is reduced below a critical value, the ferromagnetic to paramagnetic transition disappears [18]. Finite-size effects may also weaken or enhance magnetic interactions at the boundaries, as well as restrict the evolution of spin-spin correlation length. Extension of these ideas to model magnetic heterostructures, comprising of multiple magnetic and/or nonmagnetic layers, gives insight into interfacial phenomena. Many current and emerging technologies are based on this central problem. This may be very useful in understanding and exploring problems as metal-insulator transition, which is at the core of many state-of-the-art technologies. Henceforth, computational techniques, especially Ising model, can now be extended to develop and enrich science, for making new technologies. Though, its use can be said at the nascent stage, but with the advancement in computer hardware and efficient algorithms, it’s applications in areas related to spinotronics appears to be bright.