A. B. Babaev

Critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder

A Monte Carlo method is applied to simulate the static critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder. Numerical results are presented for the spin concentrations of p = 1.0, 0.95, 0.9, 0.8, 0.6 on L × L × L lattices with L = 20–60 under periodic boundary conditions. The critical temperature is determined by the Binder cumulant method. A finite-size scaling technique is used to calculate the static critical exponents α, β, γ, and ν (for specific heat, susceptibility, magnetization, and correlation length, respectively) in the range of p under study. Universality classes of critical behavior are discussed for three-dimensional diluted systems.

Tricritical point of the three-dimensional Potts model (q = 4) with quenched nonmagnetic disorder

The effect of quenched nonmagnetic impurities on the phase transitions in the three-dimensional Potts model with the number of spin states q = 4 for the case of the simple cubic lattice is studied using the Monte Carlo method. The phase transitions in this model are studied for spin density p ranging from 1.0 to 0.70. The position of the tricritical point at the phase diagram is determined.

Phase transition properties of three-dimensional systems described by diluted potts model

Monte Carlo simulations are performed to analyze phase transitions in three-dimensional systems described by the 3-state Potts model with nonmagnetic impurities. Numerical results are presented for systems with spin concentrations p = 1.00, 0.95, 0.90, 0.80, 0.70, and 0.65 on lattices of size L varying between 20 and 44. Binder’s cumulant analysis shows that the introduction of quenched disorder in the form of non-magnetic impurities induces a crossover from first-order to second-order phase transition. The finite-size scaling method is used to calculate the static critical exponents α, γ, β, and ν for specific heat, susceptibility, magnetization, and correlation length, respectively.

Phase transitions in the two-dimensional ferro- and antiferromagnetic potts models on a triangular lattice

The phase transitions in the two-dimensional ferro- and antiferromagnetic Potts models with q = 3 states of spin on a triangular lattice are studied using cluster algorithms and the classical Monte Carlo method. Systems with linear sizes L = 20–120 are considered. The method of fourth-order Binder cumulants and histogram analysis are used to discover that a second-order phase transition occurs in the ferromagnetic Potts model and a first-order phase transition takes place in the antiferromagnetic Potts model. The static critical indices of heat capacity (α), magnetic susceptibility (γ), magnetization (β), and correlation radius index (ν) are calculated for the ferromagnetic Potts model using the finite-size scaling theory.

Phase transitions in a three-dimensional diluted Potts model with 4 spin states

A Monte–Carlo method is used to study phase transitions and critical phenomena in a three-dimensional Potts model with 4 spin states and nonmagnetic impurities. Systems with linear sizes L=20–32 and spin concentrations p=1.00, 0.90, and 0.65 are studied. A fourth order Binder cumulant method is used to show that this model yields a second order phase transition for strong dilution with a spin concentration p=0.65, while the pure model (p=1.00) and the model with weak dilution (p=0.90) yield a first order phase transition. The static critical indices for the heat capacity α, susceptibility γ, magnetization β, and correlation radius ν are calculated using a finite-dimensional scaling theory.

Frustrations and phase transitions in the three-vertex Potts model with next-nearest-neighbor interactions on a triangular lattice

Phase transitions in two-dimensional arrays described by the three-vertex Potts model involving the interactions between magnetic moments located at the nearest-neighbor and next-nearest-neighbor sites of a triangular lattice are studied using the Monte Carlo method. The ratio of the next-nearest-neighbor and nearest-neighbor exchange constants r = J2/J1 is chosen within the 0–1 range. The analysis of the low temperature behavior of the entropy and density of states in the system, as well as of the fourth-order Binder cumulants, shows that in the range 0 ≤ r < 0.2, this model exhibits a first-order phase transition, whereas at r ≥ 0.2, frustrations arise in such a system.

Histogram data analysis for a three-dimensional diluted ferromagnetic 3- and 4-state potts models

On the basis of a histogram data analysis, phase transitions (PTs) in a three-dimensional diluted ferromagnetic 3- and 4-state Potts models are investigated. Systems with linear dimensions of L = 20–60 and spin concentrations of p = 1.00, 0.95, and 0.65 are studied. It is shown that the introduction of weak disorder (p ∼ 0.95) into the system in the three-dimensional Potts model with q = 3 is sufficient to change a first-order phase transition to a second-order one, whereas, in the three-dimensional Potts model with q = 4, the change of a first-order phase transition to a second-order one occurs only in the presence of strong disorder (p ∼ 0.65).

Phase transitions in a two-dimensional antiferromagnetic Potts model on a triangular lattice with next-nearest neighbor interactions

Phase transitions (PTs) and frustrations in two-dimensional structures described by a three-vertex antiferromagnetic Potts model on a triangular lattice are investigated by the Monte Carlo method with regard to nearest and next-nearest neighbors with interaction constants J1 and J2, respectively. PTs in these models are analyzed for the ratio r = J2/J1 of next-nearest to nearest exchange interaction constants in the interval |r| = 0–1.0. On the basis of the analysis of the low-temperature entropy, the density of states function of the system, and the fourth-order Binder cumulants, it is shown that a Potts model with interaction constants J1 < 0 and J2 < 0 exhibits a first-order PT in the range of 0 ⩽ r < 0.2, whereas, in the interval 0.2 ⩽ r ⩽ 1.0, frustrations arise in the system. At the same time, for J1 > 0 and J2 < 0, frustrations arise in the range 0.5 < |r| < 1.0, while, in the interval 0 ⩽ |r| ⩽ 1/3, the model exhibits a second-order PT.

Investigation of the Critical Properties in the 3D Site-Diluted Potts Model

The effect of quenched nonmagnetic impurities on phase transitions in the three-dimensional Potts model with the number of spin states q=3 is studied using the Wolff single-cluster algorithm of the Monte Carlo method. By the method of fourth-order Binder cumulants, it is demonstrated that the second-order phase transition occurs in the model under study at spin concentrations p=0.9, 0.8, 0.7, and 0.65, while the first-order phase transition is observed in the pure model (p=1.0). The static critical exponents (CEs) α (heat capacity), γ (susceptibility), β (magnetization), and ν (correlation length) are calculated based on the finite-size scaling theory.

Frustrations and Phase Transitions in Low-Dimensional Magnetic Systems

We studied magnetic orderings and frustrations on 1D chain and 2D lattices: square, triangular, kagome, and hexagonal in the Ising, 3-state Potts and standard 4-state Potts models. The spins interrelate with one another via the nearest-neighbor, the next-nearest-neighbor or the third-neighbor exchange interactions and by an external magnetic field. For problem solving we mainly calculated the entropy and specific heat using the rigorous analytical solutions for maximum eigenvalue of Kramers-Wannier transfer-matrix and exploiting computer simulation, par excellence, by Wang-Landau algorithm. Whether a system is ordered or frustrated is related to the signs and values of exchange interactions. An external magnetic field may both favor the ordering of a system and create frustrations. With the help of calculations of the entropy, the specific heat and magnetic parameters, we obtained the points and ranges of frustrations, the frustration fields and the phase transition points. The results obtained also show that the same exchange interactions my either be competing or noncompeting which depends on the topology of a lattice.