Critical properties of the anisotropic Ising model with competing interactions

The critical properties of the anisotropic Ising model with competing interactions have been investigated by Monte Carlo methods. The region of localization of the Lifshitz point on the phase diagram has been computed. Relations of the finite-size scaling theory are used to calculate the critical exponents of the heat capacity, susceptibility, and magnetization at various values of the competing interaction parameter J ₁. A crossover to a critical behavior characteristic of a multicritical point with increasing parameter J ₁ is shown to be present in the system.     

Exactly soluble planar Ising models on compressible lattices

Two planar Ising models on compressible lattices are considered. The elastic forces act in the horizontal direction only and between nearest-neighbors, but are otherwise arbitrary. The nearest-neighbor exchange interaction is taken as constant for two spins with the same column index and depending on separation for spins on the same row. In the first model (A) the transition remains continuous, and Fisher's theory of renormalized exponents applies; in the second model (B) the additional constraint that spins on a column move as a unit changes the transition to first order.     

Layered Ising models on triangular and honeycomb lattices

We study inhomogeneous Ising models on triangular and honeycomb lattices. The nearest neighbour couplings can have arbitrary strength and sign such that the coupling distribution is translationally invariant in the direction of one lattice axis, i.e. the models have a layered structure. By using a transfer matrix method we derive closed form expressions for the partition functions and free energies. The critical temperatures are calculated. Phase transitions at a finite critical temperature are universally of Ising type. Models with no phase transition may show different behaviour atT=0, which is explicitly shown for fully frustrated models on square, triangular and honeycomb lattices. Finally, generalizations to layered Ising models on more general lattices are discussed.Work performed within the research program of the Sonderforschungsbereich 125 Aachen-Jülich-Köln     

Critical temperatures of Ising models

A general and simple method for calculating critical points of two-dimensional Ising models is presented. As an example we derive Utiyama's critical relation for the generalized square lattice.     

The crossover behaviour of the susceptibility of the Ising model on anisotropic lattices

The behaviour of the reduced susceptibility χ′ of a d-dimensional Ising system in the limit t_d = tanhβJ_d → 0 is investigated. The critical exponent γ_n of (∂ⁿχ′/∂tⁿ_d)_td = 0 is, for n = 2,3, found to satisfy γ_n ? (n + 1)γ₀, and hence, by inequalities found by Liu and Stanley, γ_n = (n + 1) γ₀. For d = 2, the latter identify is valid for all n.     

Ising models on the 2 × 2 × ∞ lattices

Exact analytic solutions are presented for two 2 × 2 × ∞ Ising étagères. The first model has a simple cubic lattice with fully anisotropic interactions. The second model consists of two different types of linear chains and includes noncrossing diagonal bonds on the side faces of the 2 × 2 × ∞ parallelepiped. In both cases, the solutions are expressed through square radicals and obtained by using the obvious symmetry of the Hamiltonians, Z ₂ × C _2v , and the hidden algebraic λλ symmetry of the transfer matrix secular equations. The solution found for the second model is used to analyze the behavior of specific heat in a frustrated many-chain system. The text was submitted by author in English.     

Dimers and the Critical Ising Model on lattices of genus >1

We study the partition function of both Close-Packed Dimers and the Critical Ising Model on a square lattice embedded on a genus two surface. Using numerical and analytical methods we show that the determinants of the Kasteleyn adjacency matrices have a dependence on the boundary conditions that, for large lattice size, can be expressed in terms of genus two theta functions. The period matrix characterizing the continuum limit of the lattice is computed using a discrete holomorphic structure. These results relate in a direct way the lattice combinatorics with conformal field theory, providing new insight to the lattice regularization of conformal field theories on higher genus Riemann surfaces.     

On the critical temperature of non-periodic Ising models on hexagonal lattices

The critical temperature of layered Ising models on triangular and honeycomb lattices are calculated in simple, explicit form for arbitrary distribution of the couplings.     

On the critical temperature of non-periodic Ising models on hexagonal lattices

The critical temperature of layered Ising models on triangular and honeycomb lattices are calculated in simple, explicit form for arbitrary distribution of the couplings.     

Critical Binder cumulant of two-dimensional Ising models

The fourth-order cumulant of the magnetization, the Binder cumulant, is determined at the phase transition of Ising models on square and triangular lattices, using Monte Carlo techniques. Its value at criticality depends sensitively on boundary conditions, details of the clusters used in calculating the cumulant, and symmetry of the interactions or, here, lattice structure. Possibilities to identify generic critical cumulants are discussed.     

Analytical properties of the anisotropic cubic Ising model

We combine an exact functional relation, the inversion relation, with conventional high-temperature expansions to explore the analytic properties of the anisotropic Ising model on both the square and simple cubic lattice. In particular, we investigate the nature of the singularities that occur in partially resummed expansions of the partition function and of the susceptibility.     

Comparison of diluted antiferromagnetic Ising models on frustrated lattices in a magnetic field

We study diluted antiferromagnetic Ising models on triangular and kagome lattices in a magnetic field, using the replica-exchange Monte Carlo method. We observe seven and five plateaus in the magnetization curve of the diluted antiferromagnetic Ising model on the triangular and kagome lattices, respectively, when a magnetic field is applied. These observations contrast with the two plateaus observed in the pure model. The origin of multiple plateaus is investigated by considering the spin configuration of triangles in the diluted models. We compare these results with those of a diluted antiferromagnetic Ising model on the three-dimensional pyrochlore lattice in a magnetic field pointing in the [111] direction, sometimes referred to as the “kagome-ice” problem. We discuss the similarity and dissimilarity of the magnetization curves of the “kagome-ice” state and the two-dimensional kagome lattice.     

Some properties of random Ising models

We consider an Ising model with random magnetic fieldh _i and random nearest-neighbor couplingsJ _ij . The random variablesh _i andJ _ij are independent and identically distributed with a nice enough distribution, e.g., Gaussian. We will prove that (i) at high temperature, infinite volume correlation functions are independent on the boundary conditions and decay exponentially fast with probability 1 and (ii) for any temperature with sufficiently strong magnetic field the correlation functions are again independent on the boundary conditions and decay exponentially fast with probability 1. We also prove that the averaged magnetization of the ground state configuration of the one-dimensional Ising model with random magnetic field is zero, no matter how small is the variance of theh _i .     

Ising model on mixed two-dimensional lattices

We inform results on physical and topological magnitudes related to the ground level of Ising model on mixed two-dimensional lattices of coordination numbers 4 (Kagomé lattices) and 5 (five-point star lattices). We consider little clusters of size N, where N represents the total number of spins, subject to periodic boundary conditions. On these systems we randomly distribute ±J nearest-neighbor interactions (+J: antiferromagnetic, −J: ferromagnetic (F)). Concentration x of F interactions is varied in the interval (0,1). Two different methods are used to obtain results reported here. First, a numerical method related to multiple replicas. Second, an analytical method based on probabilistic analysis of flat and curved plaquettes. Both methods are complementary to each other. Initially, this study is restricted to calculate frustration of plaquettes and bonds, energy and bond order parameter at T=0. The results of magnitudes informed here are compared with the similar ones obtained for honeycomb, square and triangular lattices.     

Cyclic period-3 window in antiferromagnetic potts and Ising models on recursive lattices

The magnetic properties of the antiferromagnetic Potts model with two-site interaction and the antiferromagnetic Ising model with three-site interaction on recursive lattices have been studied. A cyclic period-3 window has been revealed by the recurrence relation method in the antiferromagnetic Q-state Potts model on the Bethe lattice (at Q < 2) and in the antiferromagnetic Ising model with three-site interaction on the Husimi cactus. The Lyapunov exponents have been calculated, modulated phases and a chaotic regime in the cyclic period-3 window have been found for one-dimensional rational mappings determined the properties of these systems.     

The stochastic models method applied to the critical behavior of Ising lattices

The stochastic models (SM) computer simulation method for treating manybody systems in thermodynamic equilibrium is investigated. The SM method, unlike the commonly used Metropolis Monte Carlo method, is not of a relaxation type. Thus an equilibrium configuration is constructed at once by adding particles to an initiallyempty volume with the help of a model stochastic process. The probability of the equilibrium configurations is known and this permits one to estimate the entropy directly. In the present work we greatly improve the accuracy of the SM method for the two and three-dimensional Ising lattices and extend its scope to calculate fluctuations, and hence specific heat and magnetic susceptibility, in addition to average thermodynamic quantities like energy, entropy, and magnetization. The method is found to be advantageous near the critical temperature. Of special interest are the results at the critical temperature itself, where the Metropolis method seems to be impractical. At this temperature, the average thermodynamic quantities agree well with theoretical values, for both the two and three-dimensional lattices. For the two-dimensional lattice the specific heat exhibits the expected logarithmic dependence on lattice size; the dependence of the susceptibility on lattice size is also satisfactory, leading to a ratio of critical exponents/=1.85 ±0.08. For the three-dimensional lattice the dependence of the specific heat, long-range order, and susceptibility on lattice size leads to similarly satisfactory exponents:=0.12 ±0.03,=0.30 ±0.03, and=1.32 ±0.05 (assuming =2/3).