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.     

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.     

Investigation of the influence of quenched nonmagnetic impurities on phase transitions in the three-dimensional Potts model

The influence of quenched nonmagnetic impurities on phase transitions in the three-dimensional Potts model with the number of spin states q = 3 is investigated using the Wolff single-cluster algorithm of the Monte Carlo method. The systems with linear sizes L = 20–44 at the spin concentrations p = 1.0, 0.9, 0.8, and 0.7 are analyzed. It is demonstrated with the use of the method of fourth-order Binder cumulants that the second-order phase transition occurs in the model under consideration at the spin concentrations p = 0.9, 0.8, and 0.7 and that the first-order phase transition is observed in the pure model (p = 1.0). The static critical exponents α (heat capacity), γ (susceptibility), β (magnetization), and ν (correlation length) are calculated in the framework of the finite-size scaling theory. The problem regarding the universality classes of the critical behavior of weakly diluted systems is discussed.     

Influence of quenched nonmagnetic impurities on phase transitions in a three-dimensional diluted Potts model with spin state number q = 4

The influence of quenched nonmagnetic impurities on phase transitions and critical phenomena in the 3D Potts model with the spin state number q = 4 is studied using the Monte Carlo method. Systems with the linear size L = 20–32 and spin concentrations p = 1.00, 0.90, 0.65 are considered. The fourth order Binder cumulant method is used to demonstrate that in the strongly diluted regime, a phase transition of the second kind is observed in this model for the spin concentration p = 0.65, and a phase transition of the first kind is observed for the pure (p = 1.00) and weakly diluted (p = 0.90) models. The theory of finite-dimensional scaling is used to calculate the static critical parameters of heat capacity α, susceptibility γ, magnetization β, and correlation radius ν.     

Phase transitions and critical phenomena in a three-dimensional site-diluted Potts model

The phase transitions and critical phenomena in the three-dimensional (3D) site-diluted q-state Potts models on a simple cubic lattice are explored. We systematically study the phase transitions of the models for q=3 and q=4 on the basis of Wolff high-effective algorithm by the Monte–Carlo (MC) method. The calculations are carried out for systems with periodic boundary conditions and spin concentrations p=1.00–0.65. It is shown that introducing of weak disorder (p∼0.95) into the system is sufficient to change the first order phase transition into a second order one for the 3D 3-state Potts model, while for the 3D 4-state Potts model, such a phase transformation occurs when introducing strong disorder (p∼0.65). Results for 3D pure 3-state and 4-state Potts models (p=1.00) agree with conclusions of mean field theory. The static critical exponents of the specific heat α, susceptibility γ, magnetization β, and the exponent of the correlation radius ν are calculated for the samples on the basis of finite-size scaling theory.     

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 diagram of the two-dimensional Ising model with random competing interactions

An Ising model with ferromagnetic nearest-neighbor interactions J₁ (J₁>0) and random next-nearest-neighbor interactions [+J₂ with probability p and −J₂ with probability (1−p); J₂>0] is studied within the framework of an effective-field theory based on the differential-operator technique. The order parameters are calculated, considering finite clusters with n=1,2, and 4 spins, using the standard approximation of neglecting correlations. A phase diagram is obtained in the plane temperature versus p, for the particular case J₁=J₂, showing both superantiferromagnetic (low p) and ferromagnetic (higher values of p) orderings at low temperatures.     

Variational calculations on the energy levels of graphene quantum antidots

The critical properties of the two-dimensional Ising and Blume-Capel model on directedsmall-world lattices with quenched connectivity disorder are investigated. The disordered system is simulated by applying the Monte Carlo method with heat bath update algorithm and histogram re-weighting techniques. The critical temperature, as well as the critical exponents are obtained. For both models the critical parameters have been obtained for several values of the rewiring probability p. It is found that these disorder systems do not belong to the same universality class as two-dimensional ferromagnetic model on regular lattices. In particular, the Blume-Capel model, with zero crystal field interaction, on a directedsmall-world lattice presents a second-order phase transition for p < p _c , and a first-order phase transition for p > p _c , where p _c  ≈ 0.25. The critical exponents for p < p _c are different from those of the same model on a regular lattice, but are identical to the exponents of the Ising model on directedsmall-world lattice.     

Wang-Landau study of the 3D Ising model with bond disorder

We implement a two-stage approach of the Wang-Landau algorithm to investigate the critical properties of the 3D Ising model with quenched bond randomness. In particular, we consider the case where disorder couples to the nearest-neighbor ferromagnetic interaction, in terms of a bimodal distribution of strong versus weak bonds. Our simulations are carried out for large ensembles of disorder realizations and lattices with linear sizes L in the range L=8-64L=8{-}64. We apply well-established finite-size scaling techniques and concepts from the scaling theory of disordered systems to describe the nature of the phase transition of the disordered model, departing gradually from the fixed point of the pure system. Our analysis (based on the determination of the critical exponents) shows that the 3D random-bond Ising model belongs to the same universality class with the site- and bond-dilution models, providing a single universality class for the 3D Ising model with these three types of quenched uncorrelated disorder.     

Weak Universality in the Disordered Two-Dimensional Antiferromagnetic Potts Model on a Triangular Lattice

The critical behavior of the disordered two-dimensional antiferromagnetic Potts model with the number of spin states q= 3 on a triangular lattice with disorder in the form of nonmagnetic impurities is studied by the Monte Carlo method. The critical exponents for the susceptibility γ, magnetization β, specific heat α, and correlation radius ν are calculated in the framework of the finite-size scaling theory at spin concentrations p = 0.90, 0.80, 0.70, and 0.65. It is found that the critical exponents increase with the degree of disorder, whereas the ratios and do not change, thus holding the scaling equality \(\frac{{2\beta }}{\nu } + \frac{\gamma }{\nu } = d\). Such behavior of the critical exponents is related to the weak universality of the critical behavior characteristic of disordered systems. All results are obtained using independent Monte Carlo algorithms, such as the Metropolis and Wolff algorithms.     

Monte Carlo simulation of Ising model on decagonal covering structure

We study the Ising model on a two-dimensional quasilattice developed from the decagonal covering structure. The periodic boundary conditions are applied to a patch of rhombus-like covering pattern. By means of the Monte Carlo simulation and the finite-size scaling analysis the critical temperature is estimated as 2.317±0.002. Two critical exponents are obtained being 1/v=0.992±0.003 and η=0.247±0.002, which are close to the values of the two-dimensional regular lattices as well as the Penrose tilings.     

Features of the phase transitions in the three-dimensional diluted potts model with number of spin states <Emphasis Type="Italic">q</Emphasis> = 3

The effect of frozen disorder, implemented in the form of nonmagnetic impurities, on phase transitions in the three-dimensional Potts model with number of states q = 3 has was investigated by the Monte Carlo method, using the Wolf single-cluster algorithm. Systems with linear sizes L = 20–44 were considered at spin concentrations p = 1.00–0.65. The method of fourth-order Binder cumulants was used to demonstrate that a first-order phase transition is observed for the pure model (p = 1.00) and a second-order phase transition occurs at concentrations p = 0.90, 0.80, 0.70, and 0.65.     

Critical behavior of the 3D Ising model on a poissonian random lattice

The single-cluster Monte Carlo algorithm and the reweighting technique are used to simulate the 3D ferromagnetic Ising model on 3D Voronoi-Delauney lattices. It is assumed that the coupling factor J varies with the distance r between the first neighbors as J(r)∝e^−ar, with a≥0. The critical exponents γ/ν, β/ν, and ν are calculated, and according to the present estimates for the critical exponents, we argue that this random system belongs to the same universality class of the pure 3D ferromagnetic Ising model.     

Nonmagnetic heavy-fermion superconductor, CePt3Si

We have studied the effect of heat treatment and the Ge substitution in place of Si in the recently discovered heavy-fermion superconductor CePt₃Si. The annealed CePt₃Si exhibited nonmagnetic heavy-fermion behavior instead of the antiferromagnetism (AF) found in quenched samples. The AF state was destroyed by only about 1 at.% of Ge-substitution and may not be a stable phase. Specific-heat measurements on the annealed CePt₃Si and the Ge-substituted samples revealed a large hump around 2.2 K, originally claimed as Néel temperature. Its true nature is not clarified yet but conjectured at present as a sort of quadrupolar transition rather than AF long-range order. The superconducting transition around 0.75 K was equally sharp with ΔC_p/γT_c=0.7 for clean quenched and annealed samples. The interplay between the 2.2 K-anomaly and the superconductivity is discussed.     

GPU-based single-cluster algorithm for the simulation of the Ising model

We present the GPU calculation with the common unified device architecture (CUDA) for the Wolff single-cluster algorithm of the Ising model. Proposing an algorithm for a quasi-block synchronization, we realize the Wolff single-cluster Monte Carlo simulation with CUDA. We perform parallel computations for the newly added spins in the growing cluster. As a result, the GPU calculation speed for the two-dimensional Ising model at the critical temperature with the linear size L = 4096 is 5.60 times as fast as the calculation speed on a current CPU core. For the three-dimensional Ising model with the linear size L = 256, the GPU calculation speed is 7.90 times as fast as the CPU calculation speed. The idea of quasi-block synchronization can be used not only in the cluster algorithm but also in many fields where the synchronization of all threads is required.     

Critical Relaxation of a Three-Dimensional Fully Frustrated Ising Model

Critical relaxation from the low-temperature ordered state of the three-dimensional fully frustrated Ising model on a simple cubic lattice is studied by the short-time dynamics method. Cubic systems with periodic boundary conditions and linear sizes of L = 32, 64, 96, and 128 in each crystallographic direction are studied. Calculations were carried out by the Monte Carlo method using the standard Metropolis algorithm. The static critical exponents for the magnetization and correlation radius and the dynamic critical exponents are calculated.     

Relaxation times in a finite Ising system with random impurities

Finite-size scaling effects of the Ising model with quenched random impurities are studied, focusing on critical dynamics. In contrast to the pure Ising model, disordered systems are characterized by continuous relaxation time spectra. Dynamic field theory is applied to compute the spectral densities of the magnetizationM(t) and ofM ²(t). In addition, universal cumulant ratios are calculated to second order in ^1/4, where =4–d andd<4 denotes the spatial dimension.     

Critical parameters near the ferromagnetic-paramagnetic phase transition in La0.7A0.3(Mn1−xbx)O3 (A=Sr; B=Ti and Al; x=0.0 and 0.05) compounds

The critical parameters provide important information concerning the interaction mechanisms near the paramagnetic-to-ferromagnetic transition. In this paper, we present a thorough study for the critical behavior of La_0.7A_0.3(Mn_1−xB_x)O₃ (A=Sr; B=Ti and Al; x=0.0 and 0.05) polycrystalline samples near ferromagnetic-paramagnetic phase transition temperature by analyzing isothermal magnetization data. We have analyzed our dc-magnetization data near the transition temperature with the help of the modified Arrot plot, Kouvel-Fisher method. We have determined the critical temperature T_C and the critical parameters β, γ and δ. With the values of T_C, β and γ, we plot M×(1−T/T_C)^−β vs. H×(1−T/T_C)^−γ. All the data collapse on one of the two curves. This suggests that the data below and above T_C obey scaling, following a single equation of state. Critical parameters for x=0 and x_Ti=0.05 samples are between those predicted for a 3D-Heisenberg model and mean-field theory and for x_Al=0.05 samples the values obtained for the critical parameters are close to those predicted by the mean-field theory.     

Confinement-Higgs transition in a disordered gauge theory and the accuracy threshold for quantum memory

We study the ±J random-plaquette Z₂ gauge model (RPGM) in three spatial dimensions, a three-dimensional analog of the two-dimensional ±J random-bond Ising model (RBIM). The model is a pure Z₂ gauge theory in which randomly chosen plaquettes (occurring with concentration p) have couplings with the “wrong sign” so that magnetic flux is energetically favored on these plaquettes. Excitations of the model are one-dimensional “flux tubes” that terminate at “magnetic monopoles” located inside lattice cubes that contain an odd number of wrong-sign plaquettes. Electric confinement can be driven by thermal fluctuations of the flux tubes, by the quenched background of magnetic monopoles, or by a combination of the two. Like the RBIM, the RPGM has enhanced symmetry along a “Nishimori line” in the p-T plane (where T is the temperature). The critical concentration p_c of wrong-sign plaquettes at the confinement-Higgs phase transition along the Nishimori line can be identified with the accuracy threshold for robust storage of quantum information using topological error-correcting codes: if qubit phase errors, qubit bit-flip errors, and errors in the measurement of local check operators all occur at rates below p_c, then encoded quantum information can be protected perfectly from damage in the limit of a large code block. Through Monte-Carlo simulations, we measure p_c0, the critical concentration along the T=0 axis (a lower bound on p_c), finding p_c0=.0293±.0002. We also measure the critical concentration of antiferromagnetic bonds in the two-dimensional RBIM on the T=0 axis, finding p_c0=.1031±.0001. Our value of p_c0 is incompatible with the value of p_c=.1093±.0002 found in earlier numerical studies of the RBIM, in disagreement with the conjecture that the phase boundary of the RBIM is vertical (parallel to the T axis) below the Nishimori line. The model can be generalized to a rank-r antisymmetric tensor field in d dimensions, in the presence of quenched disorder.     

Critical exponents of the Ising model on low-dimensional fractal media

The critical behavior of the Ising model on fractal substrates with noninteger Hausdorff dimension d_H<2 and infinite ramification order is studied by means of the short-time critical dynamic scaling approach. Our determinations of the critical temperatures and critical exponents β, γ, and ν are compared to the predictions of the Wilson-Fisher expansion, the Wallace-Zia expansion, the transfer matrix method, and more recent Monte Carlo simulations using finite-size scaling analysis. We also determined the effective dimension (d_ef), which plays the role of the Euclidean dimension in the formulation of the dynamic scaling and in the hyperscaling relationship d_ef=2β/ν+γ/ν. Furthermore, we obtained the dynamic exponent z of the nonequilibrium correlation length and the exponent θ that governs the initial increase of the magnetization. Our results are consistent with the convergence of the lower-critical dimension towards d=1 for fractal substrates and suggest that the Hausdorff dimension may be different from the effective dimension.