An example of a two-dimensional crystal in which each lattice point interacts with six neighbors is the so-called \triangular Ising lattice. In this work the transfer matrix method is used to nd the thermodynamic properties of the triangular Ising model with a limited number of rows, in the absence or presence of a magnetic eld. In the absence of a magnetic eld the results have been used to nd an analytical expression for the partition function of a model with any number of rows, from which the thermodynamic properties of the triangular planar lattice model are obtained as a polynomial in terms of x = e 2j tanh 2j=2 cosh2j, where J = jkT is the interaction energy of the two nearest neighbor spins. The degree of the polynomial depends on the number of rows in the model. The results are compared with the exact values. The calculated free energy is found to be very accurate, but the calculated internal energy and heat capacity show some small deviations from the exact values. In the presence of a magnetic eld the solution to the matrices is too complicated for a general analytical expression to be proposed for a model with any number of rows. However, our exact solution for the model with n = 4, 5, and 7 reveals an important point: that such results are independent of n when the eld-spin interaction energy is almost equal to or larger than that of the spin-spin interaction.
An example of a two-dimensional crystal in which each lattice point interacts with six neighbors is the so-called \triangular Ising lattice. In this work the transfer matrix method is used to nd the thermodynamic properties of the triangular Ising model with a limited number of rows, in the absence or presence of a magnetic eld. In the absence of a magnetic eld the results have been used to nd an analytical expression for the partition function of a model with any number of rows, from which the thermodynamic properties of the triangular planar lattice model are obtained as a polynomial in terms of x = e 2j tanh 2j=2 cosh2j, where J = jkT is the interaction energy of the two nearest neighbor spins. The degree of the polynomial depends on the number of rows in the model. The results are compared with the exact values. The calculated free energy is found to be very accurate, but the calculated internal energy and heat capacity show some small deviations from the exact values. In the presence of a magnetic eld the solution to the matrices is too complicated for a general analytical expression to be proposed for a model with any number of rows. However, our exact solution for the model with n = 4, 5, and 7 reveals an important point: that such results are independent of n when the eld-spin interaction energy is almost equal to or larger than that of the spin-spin interaction.
We have used the two-dimensional Ising model with a limited number of rows, but with the coordination number of four for each site, to set up the transfer matrix for the model. From the solution of such a matrix, the exact thermodynamic properties have been obtained for the model with a definite number of rows, n . We have solved the matrix for n ≤ 7 and n ≤ 10 in the presence and absence of a magnetic field, respectively. On the basis of such solutions, we have proposed an analytical expression for the partition function of the model with any number of rows in the absence of a magnetic field. The proposed expression becomes more accurate when n is larger, in such a way that it becomes very accurate for n ≥ 8 and is exact for n → ∞. Our results show that the singularity of the specific heat occurs only for the model with infinite number of rows. In the presence of a magnetic field, the solution to the matrixes is too complicated to propose a general analytical expression for the partition function of the model with any number of rows. However, our exact solution for the model with n = 4, 5, 6, and 7 reveals an important point that such results are independent of n when the field−spin interaction energy is almost equal to or larger than that of the spin−spin interaction.
For generalized 2D Ising model in an external magnetic field with the interaction of nearest neighbors, next nearest neighbors, all kinds of triple interactions and the quadruple interaction the formulas for finding free energy per lattice site in the thermodynamic limit were derived on a certain set of exact disordered solutions depending on seven parameters. Lattice models are considered with boundary conditions with a shift (similar to helical ones), and with cyclic closure of the set of all points in natural ordering. The elementary transfer matrix with nonnegative matrix elements are constructed. On the set of disorder solutions the largest eigenvalue of the transfer matrix is constant for every size of considering planar lattice, and, in particular, in the thermodynamic limit. Free energy per lattice site in the thermodynamic limit is expressed through the natural logarithm of the largest eigenvalue of transfer matrix. This largest eigenvalue can be found for a special form of eigenvector with positive components. The numerical example show the existence of nontrivial solutions of the resulting systems of equations. The system of equations and the value of free energy in the thermodynamic limit will remain the same for 2D generalized Ising models with Hamiltonians, in which the values of two (out of four) neighboring maximal spins in the natural ordering are replaced by the values of the spins at any other two lattice points adjacent in the natural ordering, this significantly expands the set of models having disordered exact solutions. The high degree of symmetry and inductive construction of the components of the eigenvectors, which disappear when going beyond the framework of the obtained set of exact solutions, is an occasion to search for phase transitions in the vicinity of this set of disordered solutions.
The partition function and the spontaneous magnetization of a two-dimensional (2D) Ising model with the next nearest interactions at zero magnetic field are derived analytically. At first, the transfer matrices are analyzed in three representations, i.e., Clifford algebraic representation, transfer tensor representation, and schematic representation, to inspect nontrivial topological structures in this system. The system is equivalent to a triangular Ising model plus an interaction along the z axis, so that the approaches developed for the 3D Ising models are modified to be appropriable for analytically solving the solution of the 2D Ising model with the next nearest interactions. The comparison with the exact solutions of other Ising lattices reveals that either the increase of the number of interactions in a unit cell or the presence/increase of topological contributions enhances the critical point of the Ising lattices. The results obtained in this work are helpful for understanding the physical properties of the 2D magnetic materials.
To solve the controversy, regarding the existence of an analytic solution to the 1-D Ising model with nearest-neighbor (NN) and next-nearest-neighbor (NNN) interactions in the presence of a magnetic field, we apply the transfer matrix method to solve the 1-D Ising model in the presence of a magnetic field, taking both NN and NNN interactions into account. We show that it is possible to write a transfer matrix only if the number of sites is even. Even in such a case, it is impossible to diagonalize the transfer matrix in an analytic form. Therefore, we employ a numerical method to obtain the eigenvalues of the transfer matrix. Moreover, the heat capacity, magnetization, and magnetic susceptibility versus temperature for different values of the competition factor (the ratio of NNN to NN interactions) are shown.
The Ising model at magnetic field H=0 is one of the most important exactly solved models in statistical mechanics, solved on the square lattice in 1944 by Onsager, and later by others on the triangular lattice; its magnetic susceptibility at H=0 continues to be an unsolved aspect of the model, however. The susceptibility can either be viewed as a sum over all correlation functions of the integrable model at H=0, or else as a second derivative of the free energy of the non-integrable Ising model in a field. Therefore, its analytic properties, though derived from the non-integrable model, can be studied through series expansion of the integrable correlation functions. We begin this process by analyzing the first four terms in the form factor expansion of the diagonal correlation functions, and after summing over the diagonal form factor expansion, the first four terms of the diagonal susceptibility expansion. We have been able to reduce the form factor and susceptibility expansion terms, given as multi-dimensional integrals, to closed-form functions in all cases. Under the limit of H going to infinity with the interaction energy E going to minus infinity, the isotropic Ising model becomes a hard particle lattice gas model, with exclusion of nearest neighbor lattice sites. In this limit, the triangular lattice Ising model becomes the exactly solved hard hexagon model, while the square lattice becomes the non-integrable hard squares model. We study in detail these models for finite lattice sizes in order to understand the differences between integrable and non-integrable models. We consider the partition functions zeros and their density for different boundary conditions, and find notable differences in the density which is attributed to an extra factorization in the transfer matrices of hard hexagons which is absent in hard squares. We also study the special point at fugacity z=-1 of hard squares where all eigenvalues of the transfer matrix are equimodular and where the grand partition function's value depends on boundary conditions. | 300 pages