Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review E

Constrained tricritical phenomena in two dimensions

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

Published Mar 15, 2005• 10 citations
Publisher

Abstract

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.

DOI
10.1103/physreve.71.036115

Related papers

Constrained tricritical Blume-Capel model in three dimensions

Youjin Deng, Henk W. J. Blöte

Using the Wolff and geometric cluster Monte Carlo methods, we investigate the tricritical Blume-Capel model in three dimensions. Since these simulations conserve the number of vacancies and thus effectively introduce a constraint, we generalize the Fisher renormalization for constrained critical behavior to tricritical systems. We observe that, indeed, the tricritical behavior is significantly modified under this constraint. For instance, at tricriticality, the specific heat has only a finite cusp and the Binder ratio assumes a different value from that in unconstrained systems. Since 3 is the upper tricritical dimensionality of Ising systems, we expect that the mean-field theory correctly predicts a number of universal parameters in three dimensions. Therefore, we calculate the partition sum of the mean-field tricritical Blume-Capel model, and accordingly obtain the exact value of the Binder ratio. Under the constraint, we show that this mean-field tricritical system reduces to the mean-field critical Ising model. However, our three-dimensional data do not agree with this mean-field prediction. Instead, they are successfully explained by the generalized Fisher renormalization mechanism.

Physical Review E•Oct 21, 2004•37 citations

Finite-size scaling in the canonical ensemble

Youjin Deng, Henk W. J. Blöte

We investigate the critical scaling behavior of finite systems in the canonical ensemble. The essential difference with the grand canonical ensemble. i.e., the constraint on the number of particles, is already known to lead to the Fisher renormalization phenomenon that modifies the thermal critical singularities. We show that, in observables that are not Fisher renormalized, it also leads to a finite-size effect governed by an exponent $y_1$ that depends on the temperature exponent $y_t$ and the dimensionality $d$ as $y_1=-|2y_t-d|$. We verify this prediction by a Monte Carlo analysis of several two-dimensional lattice models in the percolation, the Ising and the 3-state Potts universality classes.

arXiv:cond-mat/0508348•Aug 15, 2005•1 citation

Edge phase transitions of the tricritical Potts model in two dimensions

Youjin Deng, Henk W. J. Blöte

Using Monte Carlo techniques and finite-size analysis, we investigate several two-dimensional lattice models with open edges, including the Blume-Capel model and the q=1 and 3 Potts models with vacancies. At bulk tricriticality, we find that the open edges are dominated by the vacancies when the surface coupling K(s) and the chemical potential D(s) of the vacancies assume the bulk values. When K(s) and/or D(s) is sufficiently enhanced, an edge phase transition takes place, beyond which spontaneous one-dimensional order occurs on the edges. Edge phase transitions can also be induced by a surface magnetic field H(s) . We numerically determine a number of edge critical exponents and derive phase diagrams in terms of K(s) , D(s) , and H(s) . In the low-temperature region, we observe first-order transitions when K(s) and D(s) are varied; the associated hysteresis loops of surface quantities are remarkably asymmetric. Some further insight into these edge transitions is provided by the exact equivalence of the tricritical q=1 Potts model and the Ising model.

Physical Review E•Feb 10, 2005•7 citations

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

Red-bond exponents of the critical and the tricritical Ising model in three dimensions

Youjin Deng, Henk W. J. Blöte

Using the Wolff and geometric cluster algorithms and finite-size scaling analysis, we investigate the critical Ising and the tricritical Blume-Capel models with nearest-neighbor interactions on the simple-cubic lattice. The sampling procedure involves the decomposition of the Ising configuration into geometric clusters, each of which consists of a set of nearest-neighboring spins of the same sign connected with bond probability p. These clusters include the well-known Kasteleyn-Fortuin clusters as a special case for p=1-exp(-2K) , where K is the Ising spin-spin coupling. Along the critical line K=Kc , the size distribution of geometric clusters is investigated as a function of p . We observe that, unlike in the case of two-dimensional tricriticality, the percolation threshold in both models lies at pc =1-exp(-2Kc) . Further, we determine the corresponding red-bond exponents as yr =0.757(2) and 0.501(5) for the critical Ising and the tricritical Blume-Capel models, respectively. On this basis, we conjecture yr =1/2 for the latter model.

Physical Review E•Nov 24, 2004•9 citations

Critical behavior of the geometrical spin clusters and interfaces in the two-dimensional thermalized bond Ising model

S Davatolhagh, M Moshfeghian, A A Saberi

The fractal dimensions and the percolation exponents of geometrical spin clusters of like sign at criticality are obtained numerically for an Ising model with temperature dependent annealed bond dilution, also known as the thermalized bond Ising model (TBIM), in two dimensions. For this purpose, a modified Wolff single-cluster Monte Carlo simulation is used to generate equilibrium spin configurations on square lattices in the critical region. A tie-breaking rule is employed to identify nonintersecting spin cluster boundaries along the edges of the dual lattice. The values obtained for the fractal dimensions of the spanning geometrical clusters D c , and their interfaces D I , are in perfect agreement with those reported for the standard two-dimensional ferromagnetic Ising model. Furthermore, the variance of the winding angles results in a diffusivity κ = 3 for the two-dimensional thermalized bond Ising model, thus placing it in the universality class of the regular Ising model. A finite-size scaling analysis of the largest geometrical clusters results in a reliable estimation of the critical percolation exponents for the geometrical clusters in the limit of an infinite lattice size. The percolation exponents thus obtained are also found to be consistent with those reported for the regular Ising model. These consistencies are explained in terms of the Fisher renormalization relations, which express the thermodynamic critical exponents of systems with annealed bond dilution in terms of those of the regular model system.

arXiv:1202.6683•Journal of Statistical Mechanics Theory and Experiment•Feb 21, 2012•2 citations
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.
Using the Wolff and geometric cluster Monte Carlo methods, we investigate the tricritical Blume-Capel model in three dimensions. Since these simulations conserve the number of vacancies and thus effectively introduce a constraint, we generalize the Fisher renormalization for constrained critical behavior to tricritical systems. We observe that, indeed, the tricritical behavior is significantly modified under this constraint. For instance, at tricriticality, the specific heat has only a finite cusp and the Binder ratio assumes a different value from that in unconstrained systems. Since 3 is the upper tricritical dimensionality of Ising systems, we expect that the mean-field theory correctly predicts a number of universal parameters in three dimensions. Therefore, we calculate the partition sum of the mean-field tricritical Blume-Capel model, and accordingly obtain the exact value of the Binder ratio. Under the constraint, we show that this mean-field tricritical system reduces to the mean-field critical Ising model. However, our three-dimensional data do not agree with this mean-field prediction. Instead, they are successfully explained by the generalized Fisher renormalization mechanism.
We investigate the critical scaling behavior of finite systems in the canonical ensemble. The essential difference with the grand canonical ensemble. i.e., the constraint on the number of particles, is already known to lead to the Fisher renormalization phenomenon that modifies the thermal critical singularities. We show that, in observables that are not Fisher renormalized, it also leads to a finite-size effect governed by an exponent y1y_1y1​ that depends on the temperature exponent yty_tyt​ and the dimensionality ddd as y1=−∣2yt−d∣y_1=-|2y_t-d|y1​=−∣2yt​−d∣. We verify this prediction by a Monte Carlo analysis of several two-dimensional lattice models in the percolation, the Ising and the 3-state Potts universality classes.
Using Monte Carlo techniques and finite-size analysis, we investigate several two-dimensional lattice models with open edges, including the Blume-Capel model and the q=1 and 3 Potts models with vacancies. At bulk tricriticality, we find that the open edges are dominated by the vacancies when the surface coupling K(s) and the chemical potential D(s) of the vacancies assume the bulk values. When K(s) and/or D(s) is sufficiently enhanced, an edge phase transition takes place, beyond which spontaneous one-dimensional order occurs on the edges. Edge phase transitions can also be induced by a surface magnetic field H(s) . We numerically determine a number of edge critical exponents and derive phase diagrams in terms of K(s) , D(s) , and H(s) . In the low-temperature region, we observe first-order transitions when K(s) and D(s) are varied; the associated hysteresis loops of surface quantities are remarkably asymmetric. Some further insight into these edge transitions is provided by the exact equivalence of the tricritical q=1 Potts model and the Ising model.
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 .
Using the Wolff and geometric cluster algorithms and finite-size scaling analysis, we investigate the critical Ising and the tricritical Blume-Capel models with nearest-neighbor interactions on the simple-cubic lattice. The sampling procedure involves the decomposition of the Ising configuration into geometric clusters, each of which consists of a set of nearest-neighboring spins of the same sign connected with bond probability p. These clusters include the well-known Kasteleyn-Fortuin clusters as a special case for p=1-exp(-2K) , where K is the Ising spin-spin coupling. Along the critical line K=Kc , the size distribution of geometric clusters is investigated as a function of p . We observe that, unlike in the case of two-dimensional tricriticality, the percolation threshold in both models lies at pc =1-exp(-2Kc) . Further, we determine the corresponding red-bond exponents as yr =0.757(2) and 0.501(5) for the critical Ising and the tricritical Blume-Capel models, respectively. On this basis, we conjecture yr =1/2 for the latter model.
The fractal dimensions and the percolation exponents of geometrical spin clusters of like sign at criticality are obtained numerically for an Ising model with temperature dependent annealed bond dilution, also known as the thermalized bond Ising model (TBIM), in two dimensions. For this purpose, a modified Wolff single-cluster Monte Carlo simulation is used to generate equilibrium spin configurations on square lattices in the critical region. A tie-breaking rule is employed to identify nonintersecting spin cluster boundaries along the edges of the dual lattice. The values obtained for the fractal dimensions of the spanning geometrical clusters D c , and their interfaces D I , are in perfect agreement with those reported for the standard two-dimensional ferromagnetic Ising model. Furthermore, the variance of the winding angles results in a diffusivity κ = 3 for the two-dimensional thermalized bond Ising model, thus placing it in the universality class of the regular Ising model. A finite-size scaling analysis of the largest geometrical clusters results in a reliable estimation of the critical percolation exponents for the geometrical clusters in the limit of an infinite lattice size. The percolation exponents thus obtained are also found to be consistent with those reported for the regular Ising model. These consistencies are explained in terms of the Fisher renormalization relations, which express the thermodynamic critical exponents of systems with annealed bond dilution in terms of those of the regular model system.