Critical Ising on the Square Lattice Mixes in Polynomial Time

The Ising model is widely regarded as the most studied model of spin-systems in statistical physics. The focus of this paper is its dynamic (stochastic) version, the Glauber dynamics, introduced in 1963 and by now the most popular means of sampling the Ising measure. Intensive study throughout the last three decades has yielded a rigorous understanding of the spectral-gap of the dynamics on everywhere except at criticality. While the critical behavior of the Ising model has long been the focus for physicists, mathematicians have only recently developed an understanding of its critical geometry with the advent of SLE, CLE and new tools to study conformally invariant systems.     

Ising Critical Exponents on Random Trees and Graphs

We study the critical behavior of the ferromagnetic Ising model on random trees as well as so-called locally tree-like random graphs. We pay special attention to trees and graphs with a power-law offspring or degree distribution whose tail behavior is characterized by its power-law exponent τ > 2. We show that the critical inverse temperature of the Ising model equals the hyperbolic arctangent of the reciprocal of the mean offspring or mean forward degree distribution. In particular, the critical inverse temperature equals zero when ${\tau \in (2,3]}$ where this mean equals infinity. We further study the critical exponents δ, β and γ, describing how the (root) magnetization behaves close to criticality. We rigorously identify these critical exponents and show that they take the values as predicted by Dorogovstev et al. (Phys Rev E 66:016104, 2002) and Leone et al. (Eur Phys J B 28:191–197, 2002). These values depend on the power-law exponent τ, taking the same values as the mean-field Curie-Weiss model (Exactly solved models in statistical mechanics, Academic Press, London, 1982) for τ > 5, but different values for ${\tau \in (3,5)}$ .     

The Mixing Time Evolution of Glauber Dynamics for the Mean-Field Ising Model

We consider Glauber dynamics for the Ising model on the complete graph on n vertices, known as the Curie-Weiss model. It is well-known that the mixing-time in the high temperature regime (β < 1) has order n log n, whereas the mixing-time in the case β > 1 is exponential in n. Recently, Levin, Luczak and Peres proved that for any fixed β < 1 there is cutoff at time with a window of order n, whereas the mixing-time at the critical temperature β = 1 is Θ(n ^3/2). It is natural to ask how the mixing-time transitions from Θ(n log n) to Θ(n ^3/2) and finally to exp (Θ(n)). That is, how does the mixing-time behave when β = β(n) is allowed to tend to 1 as n → ∞. In this work, we obtain a complete characterization of the mixing-time of the dynamics as a function of the temperature, as it approaches its critical point β _c  = 1. In particular, we find a scaling window of order around the critical temperature. In the high temperature regime, β = 1 − δ for some 0 < δ < 1 so that δ ² n → ∞ with n, the mixing-time has order (n/δ) log(δ ² n), and exhibits cutoff with constant and window size n/δ. In the critical window, β = 1± δ, where δ ² n is O(1), there is no cutoff, and the mixing-time has order n ^3/2. At low temperature, β = 1 + δ for δ > 0 with δ ² n → ∞ and δ = o(1), there is no cutoff, and the mixing time has order . Research of J. Ding and Y. Peres was supported in part by NSF grant DMS-0605166.     

The High Temperature Ising Model on the Triangular Lattice is a Critical Bernoulli Percolation Model

We define a new percolation model by generalising the FK representation of the Ising model, and show that on the triangular lattice and at high temperatures, the critical point in the new model corresponds to the Ising model. Since the new model can be viewed as Bernoulli percolation on a random graph, our result makes an explicit connection between Ising percolation and critical Bernoulli percolation, and gives a new justification of the conjecture that the high temperature Ising model on the triangular lattice is in the same universality class as Bernoulli percolation.     

New Phase Transitions of the Ising Model on Cayley Trees

We show that the nearest neighbor Ising model on the Cayley tree exhibits new temperature–driven phase transitions. These transitions occur at various inverse temperatures different from the critical one. They are characterised by a change in the number of Gibbs states as well as by a drastic change of the behavior of free energies at these new transition points. We also consider the model in the presence of an external field and compute the free energies of translation invariant and some alternating boundary conditions.     

Holomorphic Spinor Observables in the Critical Ising Model

We introduce a new version of discrete holomorphic observables for the critical planar Ising model. These observables are holomorphic spinors defined on double covers of the original multiply connected domain. We compute their scaling limits, and show their relation to the ratios of spin correlations, thus providing a rigorous proof to a number of formulae for those ratios predicted by CFT arguments.     

Crossover Critical Behavior in the Three-Dimensional Ising Model

The character of critical behavior in physical systems depends on the range of interactions. In the limit of infinite range of the interactions, systems will exhibit mean-field critical behavior, i.e., critical behavior not affected by fluctuations of the order parameter. If the interaction range is finite, the critical behavior asymptotically close to the critical point is determined by fluctuations and the actual critical behavior depends on the particular universality class. A variety of systems, including fluids and anisotropic ferromagnets, belongs to the three-dimensional Ising universality class. Recent numerical studies of Ising models with different interaction ranges have revealed a spectacular crossover between the asymptotic fluctuation-induced critical behavior and mean-field-type critical behavior. In this work, we compare these numerical results with a crossover Landau model based on renormalization-group matching. For this purpose we consider an application of the crossover Landau model to the three-dimensional Ising model without fitting to any adjustable parameters. The crossover behavior of the critical susceptibility and of the order parameter is analyzed over a broad range (ten orders) of the scaled distance to the critical temperature. The dependence of the coupling constant on the interaction range, governing the crossover critical behavior, is discussed.     

Vanishing Critical Magnetization in the Quantum Ising Model

Adapting the recent argument of Aizenman, Duminil-Copin and Sidoravicius for the classical Ising model, it is shown here that the magnetization in the transverse-field Ising model vanishes at the critical point. The proof applies to the ground state in dimension d ≥ 2 and to positive-temperature states in dimension d ≥ 3, and relies on graphical representations as well as an infrared bound.     

Critical Value of the Quantum Ising Model on Star-Like Graphs

We present a rigorous determination of the critical value of the ground-state quantum Ising model in a transverse field, on a class of planar graphs which we call star-like. These include the junction of several copies of ℤ at a single point. Our approach is to use the graphical, or fk-, representation of the model, and the probabilistic and geometric tools associated with it. This research was carried out during the author’s Ph.D. studentship at the University of Cambridge, UK, and the Royal Institute of Technology (KTH), Sweden. The author gratefully acknowledges funding from KTH during this period. The author would also like to thank Riddarhuset, Stockholm, for generous support during his studies.     

Glauber Dynamics on Trees: Boundary Conditions and Mixing Time

We give the first comprehensive analysis of the effect of boundary conditions on the mixing time of the Glauber dynamics in the so-called Bethe approximation. Specifically, we show that the spectral gap and the log-Sobolev constant of the Glauber dynamics for the Ising model on an n-vertex regular tree with (+)-boundary are bounded below by a constant independent of n at all temperatures and all external fields. This implies that the mixing time is O(logn) (in contrast to the free boundary case, where it is not bounded by any fixed polynomial at low temperatures). In addition, our methods yield simpler proofs and stronger results for the spectral gap and log-Sobolev constant in the regime where the mixing time is insensitive to the boundary condition. Our techniques also apply to a much wider class of models, including those with hard-core constraints like the antiferromagnetic Potts model at zero temperature (proper colorings) and the hard–core lattice gas (independent sets).An extended abstract of this paper appeared under the title The Ising model on trees: Boundary conditions and mixing time in Proceedings of the 44th Annual IEEE Symposium on Foundations of Computer Science, October 2003, pp. 628–639.This work was done while this author was visiting the Departments of EECS and Statistics, University of California, Berkeley, supported in part by a Miller Visiting Professorship.Supported in part by NSF Grant CCR-0121555 and DARPA cooperative agreement F30602-00-2-0601.Supported in part by NSF Grant CCR-0121555.     

The Critical Ising Model via Kac-Ward Matrices

The Kac-Ward formula allows to compute the Ising partition function on any finite graph G from the determinant of 2^2g matrices, where g is the genus of a surface in which G embeds. We show that in the case of isoradially embedded graphs with critical weights, these determinants have quite remarkable properties. First of all, they satisfy some generalized Kramers-Wannier duality: there is an explicit equality relating the determinants associated to a graph and to its dual graph. Also, they are proportional to the determinants of the discrete critical Laplacians on the graph G, exactly when the genus g is zero or one. Finally, they share several formal properties with the Ray-Singer ${\overline{\partial}}$ -torsions of the Riemann surface in which G embeds.     

Critical Behavior of the Annealed Ising Model on Random Regular Graphs

In Giardinà et al. (ALEA Lat Am J Probab Math Stat 13(1):121–161, 2016), the authors have defined an annealed Ising model on random graphs and proved limit theorems for the magnetization of this model on some random graphs including random 2-regular graphs. Then in Can (Annealed limit theorems for the Ising model on random regular graphs, arXiv:1701.08639, 2017), we generalized their results to the class of all random regular graphs. In this paper, we study the critical behavior of this model. In particular, we determine the critical exponents and prove a non standard limit theorem stating that the magnetization scaled by \(n^{3/4}\) converges to a specific random variable, with n the number of vertices of random regular graphs.     

Phase Transition for the Ising Model on the Critical Lorentzian Triangulation

The ferromagnetic Ising model without external field on an infinite Lorentzian triangulation sampled from the uniform distribution is considered. We prove uniqueness of the Gibbs measure in the high temperature region and coexistence of at least two Gibbs measures at low temperature. The proofs are based on the disagreement percolation method and on a variant of the Peierls contour method. The critical temperature is shown to be constant a.s.     

Critical Region for Droplet Formation in the Two-Dimensional Ising Model

We study the formation/dissolution of equilibrium droplets in finite systems at parameters corresponding to phase coexistence. Specifically, we consider the 2D Ising model in volumes of size L ² , inverse temperature > _c and overall magnetization conditioned to take the value m L ² –2m v _L , where _c ^–1 is the critical temperature, m =m () is the spontaneous magnetization and v _L is a sequence of positive numbers. We find that the critical scaling for droplet formation/dissolution is when v _L ^3/2 L ^–2 tends to a definite limit. Specifically, we identify a dimensionless parameter , proportional to this limit, a non-trivial critical value _c and a function such that the following holds: For < _c , there are no droplets beyond log L scale, while for > _c , there is a single, Wulff-shaped droplet containing a fraction _c =2/3 of the magnetization deficit and there are no other droplets beyond the scale of log L. Moreover, and are related via a universal equation that apparently is independent of the details of the system.     

Universal Amplitude Ratios in the Critical Two-Dimensional Ising Model on a Torus

Using results from conformal field theory, we compute several universal amplitude ratios for the two-dimensional Ising model at criticality on a symmetric torus. These include the correlation-length ratio x =lim _L (L)/L and the first four magnetization moment ratios V _2n = ²ⁿ / ^{2 n} . As a corollary we get the first four renormalized 2n-point coupling constants for the massless theory on a symmetric torus, G*_2n . We confirm these predictions by a high-precision Monte Carlo simulation.     

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.     

Marginal Effect of Structural Defects on the Nonequilibrium Critical Behavior of the Two-Dimensional Ising Model

The influence of various initial magnetizations m0 and structural defects on the nonequilibrium critical behavior of the two-dimensional Ising model is numerically simulated by Monte Carlo methods. Based on analysis of the time dependence of magnetization and the two-time dependences of autocorrelation function and dynamic susceptibility, we revealed the influence of logarithmic corrections and the crossover phenomena of percolation behavior on the nonequilibrium characteristics and the critical exponents. Violation of the fluctuation–dissipation theorem is studied, and the limiting fluctuation–dissipation ratio is calculated for the case of high-temperature initial state. The influence of various initial states on the limiting fluctuation–dissipation ratio is investigated. The nonequilibrium critical dynamics of weakly disordered systems with spin concentrations p ≥ 0.9 is shown to belong to the universality class of the nonequilibrium critical behavior of the pure model and to be characterized by the same critical exponents and the same limiting fluctuation–dissipation ratios. The nonequilibrium critical behavior of systems with p ≤ 0.85 demonstrates that the universal characteristics of the nonequilibrium critical behavior depend on the defect concentration and the dynamic scaling is violated, which is related to the influence of the crossover effects of percolation behavior.     

The Critical Behavior of the Antiferromagnetic Ising Model with Long-Range Interaction Effects

In this paper, the antiferromagnetic Ising model with ferromagnetic long-range interaction is modeled by the Monte Carlo method. The case of ferromagnetic long-range forces decreasing by a power law is considered. The dependence of the phase-transition temperature on long-range interaction parameters is obtained. The phase diagram was constructed at different values of long-range interaction parameters. Conditions for the existence of the frustrated state of the system were revealed.     

Critical Exponents of the Diluted Ising Model between Dimensions 2 and 4

Within the massive field-theoretic renormalization-group approach the expressions for the and functions of the anisotropic mn-vector model are obtained for general space dimension d in three-loop approximation. Resumming corresponding asymptotic series, critical exponents for the case of the weakly diluted quenched Ising model (m = 1, n = 0), as well as estimates for the marginal order parameter component number m _c of the weakly diluted quenched m-vector model, are calculated as functions of d in the region 2 d < 4. Conclusions concerning the effectiveness of different resummation techniques are drawn.