Long-range connections, quantum magnets and dilute contact processes

In this article, we briefly review the critical behaviour of a long-range percolation model in which any two sites are connected with a probability that falls off algebraically with the distance. The results of this percolation transition are used to describe the quantum phase transitions in a dilute transverse Ising model at the percolation threshold p_c of the long-range connected lattice. In the similar spirit, we propose a new model of a contact process defined on the same long-range diluted lattice and explore the transitions at p_c. The long-range nature of the percolation transition allows us to evaluate some critical exponents exactly in both the above models. Moreover, mean field theory is valid for a wide region of parameter space. In either case, the strength of Griffiths McCoy singularities are tunable as the range parameter is varied.     

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.     

Kauffman Cellular Automata on Quasicrystal Topology

In this paper, we perform numerical simulations to study Kauffman cellular automata (KCA) on quasiperiod lattices. In particular, we investigate phase transition, magnetic entropy, and propagation speed of the damage on these lattices. Both the critical threshold parameter \(p_{c}\) and the critical exponents are estimated with good precision. In order to investigate the increase of statistical fluctuations and the onset of chaos in the critical region of the model, we have also defined a magnetic entropy to these systems. It is seen that the magnetic entropy behaves in a different way when one passes from the frozen regime (p < p_c) to the chaotic regime (p > p_c). For a further analysis, the robustness of the propagation of failures is checked by introducing a quenched site dilution probability q on the lattices. It is seen that the damage spreading is quite sensitive when a small fraction of the lattice sites are disconnected. A finite-size scaling analysis is employed to estimate the critical exponents. From these numerical estimates, we claim that on both pure (q =?0) and diluted (q =?0.05) quasiperiodic lattices, the KCA model belongs to the same universality class than on square lattices. Furthermore, with the aim of comparing the dynamical behavior between periodic and quasiperiodic systems, the propagation speed of the damage is also calculated for the square lattice assuming the same conditions. It is found that on square lattices the propagation speed of the damage obeys a power law as \(v\sim (p-p_{c})^{\alpha }\), whereas on quasiperiod lattices, it follows a logarithmic law as \(v \sim \ln (p-p_{c})^{\alpha }\).     

Strong-Disorder Paramagnetic-Ferromagnetic Fixed Point in the Square-Lattice ±J Ising Model

We consider the random-bond ±J Ising model on a square lattice as a function of the temperature T and of the disorder parameter p (p=1 corresponds to the pure Ising model). We investigate the critical behavior along the paramagnetic-ferromagnetic transition line at low temperatures, below the temperature of the multicritical Nishimori point at T ^*=0.9527(1), p ^*=0.89083(3). We present finite-size scaling analyses of Monte Carlo results at two temperature values, T≈0.645 and T=0.5. The results show that the paramagnetic-ferromagnetic transition line is reentrant for T<T ^*, that the transitions are continuous and controlled by a strong-disorder fixed point with critical exponents ν=1.50(4), η=0.128(8), and β=0.095(5). This fixed point is definitely different from the Ising fixed point controlling the paramagnetic-ferromagnetic transitions for T>T ^*. Our results for the critical exponents are consistent with the hyperscaling relation 2β/ν?η=d?2=0.     

Fat Fisher zeroes

We show that it is possible to determine the locus of Fisher zeroes in the thermodynamic limit for the Ising model on planar (“fat”) φ⁴ random graphs and their dual quadrangulations by matching up the real part of the high and low temperature branches of the expression for the free energy. The form of this expression for the free energy also means that series expansion results for the zeroes may be obtained with rather less effort than might appear necessary at first sight by simply reverting the series expansion of a function g(z) which appears in the solution and taking a logarithm.Unlike regular 2D lattices where numerous unphysical critical points exist with non-standard exponents, the Ising model on planar φ⁴ graphs displays only the physical transition at c=exp(−2β)=1/4 and a mirror transition at c=−1/4 both with KPZ/DDK exponents (α=−1, β=1/2, γ=2). The relation between the φ⁴ locus and that of the dual quadrangulations is akin to that between the (regular) triangular and honeycomb lattices since there is no self-duality.     

Disorder driven roughening transitions of elastic manifolds and periodic elastic media

The simultaneous effect of both disorder and crystal-lattice pinning on the equilibrium behavior of oriented elastic objects is studied using scaling arguments and a functional renormalization group technique. Our analysis applies to elastic manifolds, e.g., interfaces, as well as to periodic elastic media, e.g., charge-density waves or flux-line lattices. The competition between both pinning mechanisms leads to a continuous, disorder driven roughening transition between a flat state where the mean relative displacement saturates on large scales and a rough state with diverging relative displacement. The transition can be approached by changing the impurity concentration or, indirectly, by tuning the temperature since the pinning strengths of the random and crystal potential have in general a different temperature dependence. For D dimensional elastic manifolds interacting with either random-field or random-bond disorder a transition exists for 2<D<4, and the critical exponents are obtained to lowest order in . At the transition, the manifolds show a superuniversal logarithmic roughness. Dipolar interactions render lattice effects relevant also in the physical case of D=2. For periodic elastic media, a roughening transition exists only if the ratio p of the periodicities of the medium and the crystal lattice exceeds the critical value . For p<p _c the medium is always flat. Critical exponents are calculated in a double expansion in and and fulfill the scaling relations of random field models. Received 28 August 1998     

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 properties of the three-dimensional Ising model with quenched disorder

The static critical properties of the three-dimensional Ising model with quenched disorder are studied by the Monte-Carlo (MC) method on a simple cubic lattice, in which the quenched disorder is distributed as nonmagnetic impurities by the canonical manner. The calculations are carried out for systems with periodic boundary conditions and spin concentrations p=1.0; 0.95; 0.9; 0.8; 0.7; 0.6. The systems of non-linear sizes L×L×L, L=20-60 are researched. On the basis of the finite-size scaling (FSS) theory, the static critical exponents of specific heat α, susceptibility γ, magnetization β, and an exponent of the correlation radius in a studied interval of concentrations p are calculated. It is shown that the three-dimensional Ising model with quenched disorder has two regimes of the critical behavior universality in a dependence on nonmagnetic impurities.     

Critical finite-size scaling of magnetization distribution function for Baxter–Wu model

The distribution function P_L(m) of the order parameter for the Baxter–Wu model is studied using blocks of linear dimension L of a larger triangular lattice. At a given temperature, we use the Metropolis algorithm for the generation of a sample of configurations on the triangular lattice. The similarities and differences of this distribution with the usual cases of Ising lattices are investigated. We conclude that the present model obeys, at the critical temperature, a finite-scaling law with the known critical exponents as expected. However, our numerical data strongly indicate that the analytic form of the scaling function does not conform to the corresponding function for the usual Ising model. An analytic expression that gives a good fit is presented.     

Absence of re-entrant phase transition of the antiferromagnetic Ising model on the simple cubic lattice: Monte Carlo study of the hard-sphere lattice gas

We perform Monte Carlo simulations of the hard-sphere lattice gas on the simple cubic lattice with nearest neighbour exclusion. The critical activity is estimated, z_c = 1.0588 ± 0.0003. Using a relation between the hard-sphere lattice gas and the antiferromagnetic Ising model in an external magnetic field, we conclude that there is no re-entrant phase transition of the latter on the simple cubic lattice.     

Monte Carlo analysis of critical properties of the two-dimensional randomly site-diluted Ising model via Wang-Landau algorithm

The influence of random site dilution on the critical properties of the two-dimensional Ising model on a square lattice was explored by Monte Carlo simulations with the Wang-Landau sampling. The lattice linear size was L=20-120 and the concentration of diluted sites q=0.1,0.2,0.3. Its pure version displays a second-order phase transition with a vanishing specific heat critical exponent α, thus, the Harris criterion is inconclusive, in that disorder is a relevant or irrelevant perturbation for the critical behaviour of the pure system. The main effort was focused on the specific heat and magnetic susceptibility. We have also looked at the probability distribution of susceptibility, pseudocritical temperatures and specific heat for assessing self-averaging. The study was carried out in appropriate restricted but dominant energy subspaces. By applying the finite-size scaling analysis, the correlation length exponent ν was found to be greater than one, whereas the ratio of the critical exponents (α/ν) is negative and (γ/ν) retains its pure Ising model value supporting weak universality.     

Majority-vote on directed Erd?s-Rényi random graphs

Through Monte Carlo Simulation, the well-known majority-vote model has been studied with noise on directed random graphs. In order to characterize completely the observed order-disorder phase transition, the critical noise parameter q_c, as well as the critical exponents β/ν, γ/ν and 1/ν have been calculated as a function of the connectivity z of the random graph.     

Critical behavior of a three-dimensional dimer model

The phase transition behavior of a dimer model on a three-dimensional lattice is studied. This model is of biological interest because of its relevance to the lipid bilayer main phase transition. The model has the same kind of inactive low-temperature behavior as the exactly solvable Kasteleyn dimer model on a two-dimensional honeycomb lattice. Because of low-temperature inactivity, determination of the lowest-lying excited states allows one to locate the critical temperature. In this paper the second-lowest-lying excited states are studied and exact asymptotic results are obtained in the limit of large lattices. These results together with a finite-size scaling ansatz suggest a logarithmic divergence of the specific heat aboveT _c for the three-dimensional model. Use of the same ansatz recovers the exact divergence (α=1/2) for the two-dimensional model.     

Critical wetting in the square Ising model with a boundary field

The Ising square lattice with nearest-neighbor exchangeJ>0 and a free surface at which a boundary magnetic fieldH ₁ acts has a second-order wetting transition. We study the surface excess magnetization and the susceptibility ofL×M lattices by Monte Carlo simulation and probe the critical behavior of this wetting transition, applying finite-size scaling methods. For the cases studied, the results are not consistent with the presumably exactly known values of the critical exponents, because the asymptotic critical region has not yet been reached. Implication of our results for critical wetting in three dimensions and for the application of the present model to adsorbed wetting layers at surface steps are briefly discussed.Alexander von Humboldt-Fellow     

Nonequilibrium model on Archimedean lattices

On (4, 6, 12) and (4, 8²) Archimedean lattices, the critical properties of the majority-vote model are considered and studied using the Glauber transition rate proposed by Kwak et al. [Kwak et al., Phys. Rev. E, 75, 061110 (2007)] rather than the traditional majority-vote with noise [Oliveira, J. Stat. Phys. 66, 273 (1992)]. We obtain T _c and the critical exponents for this Glauber rate from extensive Monte Carlo studies and finite size scaling. The calculated values of the critical temperatures and Binder cumulant are T _c = 0.651(3) and U ₄ ^* = 0.612(5), and T _c = 0.667(2) and U ₄ ^* = 0.613(5), for (4, 6, 12) and (4, 8²) lattices respectively, while the exponent (ratios) β/ν, γ/ν and 1/ν are respectively: 0.105(8), 1.48(11) and 1.16(5) for (4, 6, 12); and 0.113(2), 1.60(4) and 0.84(6) for (4, 8²) lattices. The usual Ising model and the majority-vote model on previously studied regular lattices or complex networks differ from our new results.     

Dynamic critical phenomena in two-dimensional fully frustrated Coulomb gas model with disorder

The dynamic critical phenomena near depinning transition in two-dimensional fully frustrated square lattice Coulomb gas model with disorders was studied using Monte Carlo technique. The ground state of the model system with disorder σ=0.3 is a disordered state. The dependence of charge current density J on electric field E was investigated at low temperatures. The nonlinear J-E behavior near critical depinning field can be described by a scaling function proposed for three-dimensional flux line system [M.B. Luo, X. Hu, Phys. Rev. Lett. 98 (2007) 267002]. We evaluated critical exponents and found an Arrhenius creep motion for field region Ec/2<E<Ec. The scaling law of the depinning transition is also obtained from the scaling function.     

The simplified Hubbard model in one and two dimensions

Thermodynamic and dynamic properties of the one and two-dimensional simplified Hubbard model are studied. At zero temperature and half filling, no metal-insulator transition occurs for nonzero couplingU and the system is an antiferromagnetic insulator. The behavior of the gap in the single-particle density of states is investigated as a function ofU, temperature and band fillingp. For weak to intermediate coupling the gap at half filling closes for increasing temperatures. The ground state of doped lattices exhibits a metal-insulator transition at ?4d<U _c (p)≦?2d (d is the lattice dimensionality) and displays ferromagnetism without long-range order forU>U _c . The co-existence for variable temperatures and electron densities of metallic behavior and magnetic and charge-density long-range order is demonstrated. The critical temperature for long-range order is calculated for the half-filled two-dimensional case. Results for the optical conductivity and several thermodynamic properties are presented.     

Universality of the ratio of the critical amplitudes of the magnetic susceptibility in a two-dimensional ising model with nonmagnetic impurities

The behavior of the magnetic susceptibility of a two-dimensional Ising model with nonmagnetic impurities is investigated numerically. A new method for determining the critical amplitudes and critical temperature is developed. The results of a numerical investigation of the ratio of the critical amplitudes of the magnetic susceptibility are presented. It is shown that the ratio of the critical amplitudes is universal right up to impurity concentrations q ≤ 0.25 (the percolation point of a square lattice is q _c = 0.407254). The behavior of the effective critical exponent γ(q) of the magnetic susceptibility is discussed. Apparently, a transition from Ising-type universal behavior to percolation behavior should occur in a quite narrow concentration range near the percolation point of the lattice.     

The deconfinement transition in SU(2) gauge theory: A finite-size analysis

The finite-temperature deconfinement transition in pure SU(2) gauge theory is studied using finite-size analysis. The theory is formulated on a symmetric lattice with asymmetric couplings so the physical size of the spatial box can be changed continuously at a given physical temperature. Finite-size analysis based on the three-dimensional Ising model then permits the infinite-volume limit to be taken. The resulting critical temperature T_c scales according to asymptotic freedom and T_c/Λ_L = 44.0 ± 2.0.     

New results for the susceptibility of the two-dimensional Ising model at criticality

We report several new results for the susceptibility of the two-dimensional Ising model in the critical region including the exact evaluation of the constant (first regular term) in the asymptotic expansion near T_c for the ferromagnetic case.