© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2605.07867

Cluster Dynamics Stay Fast-Until Tricriticality

Minjun Jeon, Alexandros Vasilopoulos, Dong-Hee Kim, Víctor Martín-Mayor, Nikolaos G. Fytas

Published May 8, 2026• 0 citations

Abstract

Cluster Monte Carlo algorithms are widely regarded as the most effective route to overcoming critical slowing down in lattice spin systems. Whether this acceleration persists in the presence of vacancies and multicritical fluctuations, however, remains unresolved. We address this question through a systematic dynamic-scaling study of hybrid cluster-local update schemes in the two-dimensional Blume-Capel model, which exhibits a line of continuous Ising-like transitions terminating at a tricritical point. Along the entire critical line, hybrid dynamics retain the near-optimal efficiency of pure cluster updates despite the presence of annealed vacancies. Strikingly, this acceleration collapses precisely at tricriticality, where the dynamic critical exponent reverts to the local-update value. We trace this breakdown to the correlated percolation of vacancies, whose emergent system-spanning geometry obstructs nonlocal relaxation in the spin sector. Our results identify a fundamental geometric limitation of cluster acceleration at tricriticality and establish vacancy percolation as the mechanism controlling dynamic universality in hybrid Monte Carlo dynamics.

DOI: 10.48550/arxiv.2605.07867

Related papers

Percolation between vacancies in the two-dimensional Blume-Capel model

Youjin Deng, Wenan Guo, Henk W. J. Blöte

Using suitable Monte Carlo methods and finite-size scaling, we investigate the Blume-Capel model on the square lattice. We construct percolation clusters by placing nearest-neighbor bonds between vacancies with a variable bond probability p(b) . At the tricritical point, we locate the percolation threshold of these vacancy clusters at p(bc) =0.706 33 (6) . At this point, we determine the fractal dimension of the vacancy clusters as Xf =0.1308 (5) approximately equal to 21/160, and the exponent governing the renormalization flow in the p(b) direction as y(p) =0.426 (2) approximately equal to 17/40 . For bond probability p(b) > p(bc) , the vacancy clusters maintain strong critical correlations; the fractal dimension is Xf =0.0750 (2) approximately equal to 3/40 and the leading correction exponent is y(p) =-0.45 (2) approximately equal to -19/40 . The above values fit well in the Kac table for the tricritical Ising model. These vacancy clusters have much analogy with those consisting of Ising spins of the same sign, although the associated quantities rho and magnetization m are energylike and magnetic quantities, respectively. However, along the critical line of the Blume-Capel model, the vacancies are more or less uniformly distributed over the whole lattice. In this case, no critical percolation correlations are observed in the vacancy clusters, at least in the physical region p(b) < or = 1 .

Physical Review E•Jul 1, 2005•17 citations

Constrained tricritical phenomena in two dimensions

Youjin Deng, Jouke R. Heringa, Henk W. J. Blöte

We investigate several tricritical models on the square lattice by means of Monte Carlo simulations. These include the Blume-Capel model, Baxter's hard-square model, and the q=1 , 3, and 4 Potts models with vacancies. We use a combination of the Wolff and geometric cluster methods, which conserves the total number of vacancies or lattice-gas particles and suppresses critical slowing down. Several quantities are sampled, such as the specific heat C and the structure factor C(s) , which accounts for the large-scale spatial inhomogeneity of the energy fluctuations. We find that the constraint strongly modifies some of the critical singularities. For instance, the specific heat C reaches a finite value at tricriticality, while C(s) remains divergent as in the unconstrained system. We are able to explain these observed constrained phenomena on the basis of the Fisher renormalization mechanism generalized to include a subleading relevant thermal scaling field. In this context, we find that, under the constraint, the leading thermal exponent y(t1) is renormalized to 2- y(t1) , while the subleading exponent y(t2) remains unchanged.

Physical Review E•Mar 15, 2005•10 citations

Dynamic critical exponent z of the three-dimensional Ising universality class: Monte Carlo simulations of the improved Blume-Capel model

Martin Hasenbusch B

We study purely dissipative relaxational dynamics in the three-dimensional Ising universality class. To this end, we simulate the improved Blume-Capel model on the simple cubic lattice by using local algorithms. We perform a finite size scaling analysis of the integrated autocorrelation time of the magnetic susceptibility in equilibrium at the critical point. We obtain z=2.0245(15) for the dynamic critical exponent. As a complement, fully magnetized configurations are suddenly quenched to the critical temperature, giving consistent results for the dynamic critical exponent. Furthermore, our estimate of z is fully consistent with recent field theoretic results.

arXiv:1908.01702•Physical review. E•Feb 24, 2020•39 citations

Critical behavior of the spin-3/2 Blume-Capel model on a random two-dimensional lattice

F. W. S. Lima, J. A. Plascak

We investigate the critical properties of the spin-3/2 Blume-Capel model in two dimensions on a random lattice with quenched connectivity disorder. The disordered system is simulated by applying the cluster hybrid Monte Carlo update algorithm and re-weighting techniques. We calculate the critical temperature as well as the critical point exponents gamma/ n, b / n, a /nu, and nu. We find that, contrary of what happens to the spin-1/2 case, this random system does not belong to the same universality class as the regular two-dimensional ferromagnetic model.

arXiv:cond-mat/0604145•Brazilian Journal of Physics•Sep 1, 2006•7 citations

Critical Behavior of the Spin-3/2 Blume-Capel Model on a Random Two-Dimensional Lattice

F. W. S. Lima, J. A. Plascak

We investigate the critical properties of the spin-3/2 Blume-Capel model in two dimensions on a random lattice with quenched connectivity disorder. The disordered system is simulated by applying the cluster hybrid Monte Carlo update algorithm and re-weighting techniques. We calculate the critical temperature as well as the critical point exponents gamma/ n, b / n, a /nu, and nu. We find that, contrary of what happens to the spin-1/2 case, this random system does not belong to the same universality class as the regular two-dimensional ferromagnetic model.

•Jan 1, 2006•4 citations

Invaded cluster algorithm for a tricritical point in a diluted Potts model

Ivan Balog, Katarina Uzelac

The invaded cluster approach is extended to the two-dimensional Potts model with annealed vacancies by using the random-cluster representation. Geometrical arguments are used to propose the algorithm which converges to the tricritical point in the two-dimensional parameter space spanned by temperature and the chemical potential of the vacancies. The tricritical point is identified as a simultaneous onset of the percolation of a Fortuin-Kasteleyn cluster and of a percolation of the "geometrical disorder cluster." The location of the tricritical point and the concentration of vacancies for q=1,2,3 are found to be in good agreement with the best known results. Scaling properties of the percolating scaling cluster and related critical exponents are also presented.

arXiv:cond-mat/0703759•Physical Review E•Jul 3, 2007•8 citations

Cluster Monte Carlo algorithms are widely regarded as the most effective route to overcoming critical slowing down in lattice spin systems. Whether this acceleration persists in the presence of vacancies and multicritical fluctuations, however, remains unresolved. We address this question through a systematic dynamic-scaling study of hybrid cluster-local update schemes in the two-dimensional Blume-Capel model, which exhibits a line of continuous Ising-like transitions terminating at a tricritical point. Along the entire critical line, hybrid dynamics retain the near-optimal efficiency of pure cluster updates despite the presence of annealed vacancies. Strikingly, this acceleration collapses precisely at tricriticality, where the dynamic critical exponent reverts to the local-update value. We trace this breakdown to the correlated percolation of vacancies, whose emergent system-spanning geometry obstructs nonlocal relaxation in the spin sector. Our results identify a fundamental geometric limitation of cluster acceleration at tricriticality and establish vacancy percolation as the mechanism controlling dynamic universality in hybrid Monte Carlo dynamics.

Using suitable Monte Carlo methods and finite-size scaling, we investigate the Blume-Capel model on the square lattice. We construct percolation clusters by placing nearest-neighbor bonds between vacancies with a variable bond probability p(b) . At the tricritical point, we locate the percolation threshold of these vacancy clusters at p(bc) =0.706 33 (6) . At this point, we determine the fractal dimension of the vacancy clusters as Xf =0.1308 (5) approximately equal to 21/160, and the exponent governing the renormalization flow in the p(b) direction as y(p) =0.426 (2) approximately equal to 17/40 . For bond probability p(b) > p(bc) , the vacancy clusters maintain strong critical correlations; the fractal dimension is Xf =0.0750 (2) approximately equal to 3/40 and the leading correction exponent is y(p) =-0.45 (2) approximately equal to -19/40 . The above values fit well in the Kac table for the tricritical Ising model. These vacancy clusters have much analogy with those consisting of Ising spins of the same sign, although the associated quantities rho and magnetization m are energylike and magnetic quantities, respectively. However, along the critical line of the Blume-Capel model, the vacancies are more or less uniformly distributed over the whole lattice. In this case, no critical percolation correlations are observed in the vacancy clusters, at least in the physical region p(b) < or = 1 .

We investigate several tricritical models on the square lattice by means of Monte Carlo simulations. These include the Blume-Capel model, Baxter's hard-square model, and the q=1 , 3, and 4 Potts models with vacancies. We use a combination of the Wolff and geometric cluster methods, which conserves the total number of vacancies or lattice-gas particles and suppresses critical slowing down. Several quantities are sampled, such as the specific heat C and the structure factor C(s) , which accounts for the large-scale spatial inhomogeneity of the energy fluctuations. We find that the constraint strongly modifies some of the critical singularities. For instance, the specific heat C reaches a finite value at tricriticality, while C(s) remains divergent as in the unconstrained system. We are able to explain these observed constrained phenomena on the basis of the Fisher renormalization mechanism generalized to include a subleading relevant thermal scaling field. In this context, we find that, under the constraint, the leading thermal exponent y(t1) is renormalized to 2- y(t1) , while the subleading exponent y(t2) remains unchanged.

We study purely dissipative relaxational dynamics in the three-dimensional Ising universality class. To this end, we simulate the improved Blume-Capel model on the simple cubic lattice by using local algorithms. We perform a finite size scaling analysis of the integrated autocorrelation time of the magnetic susceptibility in equilibrium at the critical point. We obtain z=2.0245(15) for the dynamic critical exponent. As a complement, fully magnetized configurations are suddenly quenched to the critical temperature, giving consistent results for the dynamic critical exponent. Furthermore, our estimate of z is fully consistent with recent field theoretic results.

We investigate the critical properties of the spin-3/2 Blume-Capel model in two dimensions on a random lattice with quenched connectivity disorder. The disordered system is simulated by applying the cluster hybrid Monte Carlo update algorithm and re-weighting techniques. We calculate the critical temperature as well as the critical point exponents gamma/ n, b / n, a /nu, and nu. We find that, contrary of what happens to the spin-1/2 case, this random system does not belong to the same universality class as the regular two-dimensional ferromagnetic model.

We investigate the critical properties of the spin-3/2 Blume-Capel model in two dimensions on a random lattice with quenched connectivity disorder. The disordered system is simulated by applying the cluster hybrid Monte Carlo update algorithm and re-weighting techniques. We calculate the critical temperature as well as the critical point exponents gamma/ n, b / n, a /nu, and nu. We find that, contrary of what happens to the spin-1/2 case, this random system does not belong to the same universality class as the regular two-dimensional ferromagnetic model.

The invaded cluster approach is extended to the two-dimensional Potts model with annealed vacancies by using the random-cluster representation. Geometrical arguments are used to propose the algorithm which converges to the tricritical point in the two-dimensional parameter space spanned by temperature and the chemical potential of the vacancies. The tricritical point is identified as a simultaneous onset of the percolation of a Fortuin-Kasteleyn cluster and of a percolation of the "geometrical disorder cluster." The location of the tricritical point and the concentration of vacancies for q=1,2,3 are found to be in good agreement with the best known results. Scaling properties of the percolating scaling cluster and related critical exponents are also presented.