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

arXiv:2609.29323

Two temperature scales in the Ising model on the $\{5,4\}$ hyperbolic lattice with free boundaries: susceptibility peak and boundary-induced order

S. Jaroszewicz, Nahuel Daniel Mendez, Maria P. Beccar-Varela, Maria C. Mariani

Published Sep 24, 2026• 0 citations
arXivHTMLPDF

Abstract

On the $\{5,4\}$ hyperbolic lattice, the outermost generation holds 73\% of the sites at all system sizes, making macroscopic averages strictly boundary-dependent. We study the ferromagnetic Ising model on this geometry with free boundaries by Monte Carlo simulation, demonstrating that conventional observables cease to identify a single critical scale. Instead, the finite lattices are organized by two distinct temperature scales. The susceptibility maximum identifies the lower transition scale at $T_{c2}=1.4782(15)J/k_{B}$ without extrapolation. However, the order-parameter distribution lacks a conventional fixed point at this scale. A distinct pseudocritical scale emerges instead at $T^{*}\simeq1.67(3)$ inside the boundary-sensitive intermediate phase, where Binder cumulant curves cross and effective exponents are compatible with mean-field criticality. Between $T_{c2}$ and $T^{*}$, the interior amplifies boundary fluctuations into induced order. Furthermore, because the volume grows exponentially with depth, finite-size scaling must be formulated in the generation index rather than the total number of sites. The lack of asymptotic power-law convergence in the scaling stretch confirms through an independent observable that $T^{*}$ is a finite-size crossover rather than a thermodynamic fixed point. Finally, by applying a coherent field to the outermost generation alone, we recover the upper thermodynamic transition at $T_{pt}=2.81(1)J/k_{B}$, demonstrating that it remains accessible under an appropriate boundary perturbation.

DOI
10.48550/arxiv.2609.29323

Related papers

Boundary Condition dependent Universality Classes on a Hyperbolic Lattice

Anonymous, Debasish Banerjee, Arnab Kundu +1 more

We study the ferromagnetic Ising model on finite hyperbolic tessellations of Euclidean AdS$_2$ with open and wired boundary conditions. Because a finite fraction of spins remains at the boundary as the lattice grows, these conditions select distinct thermodynamic behaviours. In the $\{5,4\}$ tessellation, Monte-Carlo simulations using efficient worm algorithms on open boundaries (OBC) yield a transition near $β_c J=0.632(3)$, with susceptibility data collapsing under scaling by the total number of spins $N$. The fitted finite-size exponents, $1/\barν \simeq 0.162$ and $γ/\barν \simeq 0.743$, differ substantially from the mean-field volume scaling. Wired boundaries (WBC), which correlate boundary spins, instead show a transition around $β_c J=0.348(4)$ to an ordered phase consistent with the mean-field exponents. While the mean-field criticality with WBC is consistent with previous studies and with the suppression of independent boundary fluctuations, the OBC exponents hint at the presence of a new universality class. Results from other tessellations support the robustness of the observed scaling with OBC. We discuss a possible route to interpolate between the different boundary conditions. Our findings show that boundary dynamics must be specified when characterizing critical behaviour on hyperbolic lattices.

arXiv:2609.37325•Sep 29, 2026•0 citations

A nonequilibrium system on a restricted scale-free network

Rafael Alves Dumer, Marta Arantes Godoy

The nonequilibrium Ising model on a restricted scale-free network has been studied with one- and two-spin flip competing dynamics employing Monte Carlo simulations. The dynamics present in the system can be defined by the probability $q$ in which the one-spin flip process simulate the contact with a heat bath at a given temperature $T$, and with a probability ($1-q$) the two-spin flip process mimics the system subjected to an external flux of energy into it. The system network is described by a power-law degree distribution in the form $P(k)\sim k^{-α}$, and the restriction is made by fixing the maximum, $k_{m}$, and minimum, $k_{0}$, degree on distribution for the whole network size. This restriction keeps finite the second and fourth moment of degree distribution, allowing us to obtain a finite critical point for any value of $α$. For these critical points, we have calculated the thermodynamic quantities of the system, such as, the total ${m}_{N}^{F}$ and staggered ${m}_{N}^{AF}$ magnetizations per spin, susceptibility $χ_{N}$, and reduced fourth-order Binder cumulant ${U}_{N}$, for several values of lattice size $N$ and exponent $1\leα\le5$. Therefore, the phase diagram was built and a self-organization phenomena is observed from the transitions between antiferromagnetic AF to paramagnetic P, and P to ferromagnetic F phases. Using the finite-size scaling theory, we also obtained the critical exponents for the system, and a mean-field critical behavior is observed, exhibiting the same universality class of the system on the equilibrium and out of it.

arXiv:2306.04780•Jun 7, 2023•0 citations

The scaling window of the 5D Ising model with free boundary conditions

P. H. Lundow, K. Markström

The five-dimensional Ising model with free boundary conditions has recently received a renewed interest in a debate concerning the finite-size scaling of the susceptibility near the critical temperature. We provide evidence in favour of the conventional scaling picture, where the susceptibility scales as L2 inside a critical scaling window of width 1/L2. Our results are based on Monte Carlo data gathered on system sizes up to L=79 (ca. three billion spins) for a wide range of temperatures near the critical point. We analyse the magnetisation distribution, the susceptibility and also the scaling and distribution of the size of the Fortuin–Kasteleyn cluster containing the origin. The probability of this cluster reaching the boundary determines the correlation length, and its behaviour agrees with the mean field critical exponent δ=3, that the scaling window has width 1/L2.

arXiv:1601.04053•Nuclear Physics B•Aug 8, 2016•16 citations

Monte Carlo studies of finite-size scaling in the Gibbs ensemble at criticality for the two-dimensional Ising model

K. K. Mon

In this paper we study finite-size scaling at criticality in the Gibbs ensemble with Monte Carlo simulations. A recent lattice-gas version of Gibbs-ensemble sampling for the two-dimensional Ising-model magnet with nearest-neighbor ferromagnetic interactions on two subsystems of L\ifmmode\times\else\texttimes\fi{}L lattices is considered. Both square and triangular lattices with a range of lattice sizes are simulated at bulk criticality. We are interested in the question of universality for the fourth-order cumulant of the order-parameter distribution in the Gibbs ensemble. Estimates of -1.954\ifmmode\pm\else\textpm\fi{}0.005 and -1.9543\ifmmode\pm\else\textpm\fi{}0.005 are obtained for the square and triangular lattices, respectively. These values confirm universality within the Gibbs ensemble, but are statistically different from the canonical-ensemble values of \ensuremath{\sim}-1.83. We point out that one must interpret this apparent ensemble dependence with care. Unlike in the canonical ensemble, the order parameter in the Gibbs ensemble is the subsystem magnetization. Its fourth-order cumulant is not simply related to the usual fourth magnetic-field derivative of the total free energy. Standard finite-size scaling analysis of data for the total specific heat, subsystem magnetization, and susceptibility are also considered and found to be in good agreement with finite-size scaling.

Physical review. B, Condensed matter•Mar 1, 1993•8 citations

When correlations exceed system size: finite-size scaling in free boundary conditions above the upper critical dimension

Yulian Honchar, Bertrand Berche, Yurij V. Holovatch +1 more

We progress finite-size scaling in systems with free boundary conditions above their upper critical dimension, where in the thermodynamic limit critical scaling is described by mean-field theory. Recent works show that the correlation length is not bound by the system's physical size, a belief that long held sway. Instead, two scaling regimes can be observed — at the critical and pseudo-critical temperatures. We demonstrate that both are manifest for free boundaries. We use numerical simulations of the d = 5 Ising model to analyse the magnetization, susceptibility, magnetization Fourier modes and the partition function zeros. While some of the response functions hide the dual finite-size scaling, the precision enabled by the analysis of Lee–Yang zeros allows this be brought to the fore. In particular, finite-size scaling of leading zeros at the pseudo-critical point confirms recent predictions coming from correlations exceeding the system size. This paper is dedicated to Jaroslav Ilnytskyi on the occasion of his 60th birthday.

arXiv:2311.11721•Condensed Matter Physics•Mar 28, 2024•2 citations

Critical Behavior Simulation of the Two-Dimensional Ising Model—A Computational Physics Study Based on the Monte Carlo Method

佳钢 李

To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.

Modern Physics•Jan 1, 2026•0 citations
Two temperature scales in the Ising model on the {5,4}\{5,4\}{5,4} hyperbolic lattice with free boundaries: susceptibility peak and boundary-induced order
On the {5,4}\{5,4\}{5,4} hyperbolic lattice, the outermost generation holds 73% of the sites at all system sizes, making macroscopic averages strictly boundary-dependent. We study the ferromagnetic Ising model on this geometry with free boundaries by Monte Carlo simulation, demonstrating that conventional observables cease to identify a single critical scale. Instead, the finite lattices are organized by two distinct temperature scales. The susceptibility maximum identifies the lower transition scale at Tc2=1.4782(15)J/kBT_{c2}=1.4782(15)J/k_{B}Tc2​=1.4782(15)J/kB​ without extrapolation. However, the order-parameter distribution lacks a conventional fixed point at this scale. A distinct pseudocritical scale emerges instead at T∗≃1.67(3)T^{*}\simeq1.67(3)T∗≃1.67(3) inside the boundary-sensitive intermediate phase, where Binder cumulant curves cross and effective exponents are compatible with mean-field criticality. Between Tc2T_{c2}Tc2​ and T∗T^{*}T∗, the interior amplifies boundary fluctuations into induced order. Furthermore, because the volume grows exponentially with depth, finite-size scaling must be formulated in the generation index rather than the total number of sites. The lack of asymptotic power-law convergence in the scaling stretch confirms through an independent observable that T∗T^{*}T∗ is a finite-size crossover rather than a thermodynamic fixed point. Finally, by applying a coherent field to the outermost generation alone, we recover the upper thermodynamic transition at Tpt=2.81(1)J/kBT_{pt}=2.81(1)J/k_{B}Tpt​=2.81(1)J/kB​, demonstrating that it remains accessible under an appropriate boundary perturbation.
We study the ferromagnetic Ising model on finite hyperbolic tessellations of Euclidean AdS2_22​ with open and wired boundary conditions. Because a finite fraction of spins remains at the boundary as the lattice grows, these conditions select distinct thermodynamic behaviours. In the {5,4}\{5,4\}{5,4} tessellation, Monte-Carlo simulations using efficient worm algorithms on open boundaries (OBC) yield a transition near βcJ=0.632(3)β_c J=0.632(3)βc​J=0.632(3), with susceptibility data collapsing under scaling by the total number of spins NNN. The fitted finite-size exponents, 1/νˉ≃0.1621/\barν \simeq 0.1621/νˉ≃0.162 and γ/νˉ≃0.743γ/\barν \simeq 0.743γ/νˉ≃0.743, differ substantially from the mean-field volume scaling. Wired boundaries (WBC), which correlate boundary spins, instead show a transition around βcJ=0.348(4)β_c J=0.348(4)βc​J=0.348(4) to an ordered phase consistent with the mean-field exponents. While the mean-field criticality with WBC is consistent with previous studies and with the suppression of independent boundary fluctuations, the OBC exponents hint at the presence of a new universality class. Results from other tessellations support the robustness of the observed scaling with OBC. We discuss a possible route to interpolate between the different boundary conditions. Our findings show that boundary dynamics must be specified when characterizing critical behaviour on hyperbolic lattices.
The nonequilibrium Ising model on a restricted scale-free network has been studied with one- and two-spin flip competing dynamics employing Monte Carlo simulations. The dynamics present in the system can be defined by the probability qqq in which the one-spin flip process simulate the contact with a heat bath at a given temperature TTT, and with a probability (1−q1-q1−q) the two-spin flip process mimics the system subjected to an external flux of energy into it. The system network is described by a power-law degree distribution in the form P(k)∼k−αP(k)\sim k^{-α}P(k)∼k−α, and the restriction is made by fixing the maximum, kmk_{m}km​, and minimum, k0k_{0}k0​, degree on distribution for the whole network size. This restriction keeps finite the second and fourth moment of degree distribution, allowing us to obtain a finite critical point for any value of ααα. For these critical points, we have calculated the thermodynamic quantities of the system, such as, the total mNF{m}_{N}^{F}mNF​ and staggered mNAF{m}_{N}^{AF}mNAF​ magnetizations per spin, susceptibility χNχ_{N}χN​, and reduced fourth-order Binder cumulant UN{U}_{N}UN​, for several values of lattice size NNN and exponent 1≤α≤51\leα\le51≤α≤5. Therefore, the phase diagram was built and a self-organization phenomena is observed from the transitions between antiferromagnetic AF to paramagnetic P, and P to ferromagnetic F phases. Using the finite-size scaling theory, we also obtained the critical exponents for the system, and a mean-field critical behavior is observed, exhibiting the same universality class of the system on the equilibrium and out of it.
The five-dimensional Ising model with free boundary conditions has recently received a renewed interest in a debate concerning the finite-size scaling of the susceptibility near the critical temperature. We provide evidence in favour of the conventional scaling picture, where the susceptibility scales as L2 inside a critical scaling window of width 1/L2. Our results are based on Monte Carlo data gathered on system sizes up to L=79 (ca. three billion spins) for a wide range of temperatures near the critical point. We analyse the magnetisation distribution, the susceptibility and also the scaling and distribution of the size of the Fortuin–Kasteleyn cluster containing the origin. The probability of this cluster reaching the boundary determines the correlation length, and its behaviour agrees with the mean field critical exponent δ=3, that the scaling window has width 1/L2.
In this paper we study finite-size scaling at criticality in the Gibbs ensemble with Monte Carlo simulations. A recent lattice-gas version of Gibbs-ensemble sampling for the two-dimensional Ising-model magnet with nearest-neighbor ferromagnetic interactions on two subsystems of L\ifmmode\times\else\texttimes\fi{}L lattices is considered. Both square and triangular lattices with a range of lattice sizes are simulated at bulk criticality. We are interested in the question of universality for the fourth-order cumulant of the order-parameter distribution in the Gibbs ensemble. Estimates of -1.954\ifmmode\pm\else\textpm\fi{}0.005 and -1.9543\ifmmode\pm\else\textpm\fi{}0.005 are obtained for the square and triangular lattices, respectively. These values confirm universality within the Gibbs ensemble, but are statistically different from the canonical-ensemble values of \ensuremath{\sim}-1.83. We point out that one must interpret this apparent ensemble dependence with care. Unlike in the canonical ensemble, the order parameter in the Gibbs ensemble is the subsystem magnetization. Its fourth-order cumulant is not simply related to the usual fourth magnetic-field derivative of the total free energy. Standard finite-size scaling analysis of data for the total specific heat, subsystem magnetization, and susceptibility are also considered and found to be in good agreement with finite-size scaling.
We progress finite-size scaling in systems with free boundary conditions above their upper critical dimension, where in the thermodynamic limit critical scaling is described by mean-field theory. Recent works show that the correlation length is not bound by the system's physical size, a belief that long held sway. Instead, two scaling regimes can be observed — at the critical and pseudo-critical temperatures. We demonstrate that both are manifest for free boundaries. We use numerical simulations of the d = 5 Ising model to analyse the magnetization, susceptibility, magnetization Fourier modes and the partition function zeros. While some of the response functions hide the dual finite-size scaling, the precision enabled by the analysis of Lee–Yang zeros allows this be brought to the fore. In particular, finite-size scaling of leading zeros at the pseudo-critical point confirms recent predictions coming from correlations exceeding the system size. This paper is dedicated to Jaroslav Ilnytskyi on the occasion of his 60th birthday.
To demonstrate the application of computational physics in modern statistical physics, this paper investigates the critical behavior of a finite-size ferromagnetic system based on the two-dimensional square-lattice Ising model.A nearest-neighbor Ising Hamiltonian is constructed, and periodic boundary conditions are adopted to reduce finite-boundary effects.The Metropolis Monte Carlo algorithm is used to sample spin configurations at different temperatures, while energy, magnetization, specific heat, magnetic susceptibility and Binder cumulant are calculated to describe the transition from an ordered state to a disordered state.The simulation results show that the system exhibits clear long-range order at low temperatures, where the order parameter remains high.As temperature increases, thermal fluctuations become stronger, magnetization decreases rapidly, and the system gradually evolves into a disordered state.The specific heat curve shows a pronounced peak at approximately (T = 2.269), while the susceptibility peak shifts slightly toward the high-temperature side due to finite-size effects and the discrete temperature grid.Overall, the simulated critical region is consistent with the theoretical critical temperature of the two-dimensional square-lattice Ising model.Spin-configuration analysis further reveals enhanced fluctuations and domain reconstruction near the critical region.In addition, a thermalization criterion and an extended sampling test for L = 48 near the critical region are added to examine the influence of finite sampling and equilibration.These results indicate that Monte Carlo simulation of the two-dimensional Ising model can effectively reproduce continuous phase transition, critical fluctuation and finite-size effects, providing a basic computational case for understanding critical phenomena in complex systems and data-driven physical analysis.