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:2505.04298

Magnetization-resolved density of states and quasi-first order transition in the two-dimensional random bond Ising model: an entropic sampling study

Yi Liu, Ding Wang, Xin Wang, Dao‐Xin Yao, Lei‐Han Tang

Published May 7, 2025• 0 citations
arXivHTMLPDF

Abstract

Systems with quenched disorder possess complex energy landscapes that are challenging to explore under the conventional Monte Carlo method. In this work, we implement an efficient entropy sampling scheme for accurate computation of the entropy function in low-energy regions. The method is applied to the two-dimensional $\pm J$ random-bond Ising model, where frustration is controlled by the fraction $p$ of ferromagnetic bonds. We investigate the low-temperature paramagnetic--ferromagnetic phase boundary below the multicritical point at $T_N = 0.9530(4)$, $P_N = 0.89078(8)$, as well as the zero-temperature ferromagnetic--spin-glass transition. Finite-size scaling analysis reveals that the phase boundary for $T < T_N$ exhibits reentrant behavior. By analyzing the evolution of the magnetization-resolved density of states $g(E, M)$ and ground-state spin configurations against increasing frustration, we provide strong evidence that the zero-temperature transition is quasi-first order. Finite-size scaling conducted on the spin-glass side supports the validity of $β= 0$, with a correlation length exponent $ν= 1.50(8)$. Our results provide new insights into the nature of the ferromagnetic-to-spin-glass phase transition in an extensively degenerate ground state.

DOI
10.48550/arxiv.2505.04298

Related papers

Magnetization-resolved density of states and mixed-order transition in the two-dimensional random bond Ising model: an entropic sampling study

Yi Liu, Ding Wang, Xin Wang +2 more

Abstract Systems with quenched disorder possess complex energy landscapes that are challenging to explore under conventional Monte Carlo methods. In this work, we implement an efficient entropy sampling scheme for accurate computation of the entropy function in low-energy regions. The method is applied to the two-dimensional ± J random-bond Ising model, where frustration is controlled by the fraction p of ferromagnetic bonds. We investigate the low-temperature paramagnetic–ferromagnetic phase boundary below the multicritical point at T N = 0.9530(4), P N = 0.89078(8), as well as the zero-temperature ferromagnetic–spin-glass transition. Finite-size scaling analysis reveals that the phase boundary for T < T N exhibits reentrant behavior. By analyzing the evolution of the magnetization-resolved density of states g ( E , M ) and ground-state spin configurations against increasing frustration, we provide strong evidence that the zero-temperature transition is a mixed-order. Finite-size scaling conducted on the spin-glass side supports the validity of β = 0, where β is the magnetization exponent, with a correlation length exponent ν = 1.50(8). Our results provide new insights into the nature of the ferromagnetic-to-spin-glass phase transition in an extensively degenerate ground state.

Communications in Theoretical Physics•Jun 14, 2025•1 citation

Low-temperature transition of 2d random-bond Ising model and quantum infinite randomness

Akshat Pandey, Aditya Mahadevan, A. Alan Middleton +1 more

At low temperatures, the classical two-dimensional random bond Ising model undergoes a frustration-driven ferromagnet-to-paramagnet transition controlled by a zero-temperature fixed point separating ferromagnet and spin glass phases. We show that this critical point can be understood through a renormalization group transformation that constructs the ground state of the Ising model through a sequence of Hamiltonians that, starting with an unfrustrated model, iteratively adds in frustration until the target Hamiltonian is reached. Via a mapping of the thermodynamics of the 2d Ising model to the spectral properties of a related Hermitian matrix -- the Hamiltonian of a noninteracting quantum problem -- this RG procedure corresponds to an iterative diagonalization of the quantum Hamiltonian. The flow toward zero temperature in the Ising picture manifests as a flow toward infinite randomness in the spectrum of the quantum Hamiltonian, with the log gap of the Hamiltonian scaling as a power of the system size: $\log \varepsilon_{\it min}^{-1} \sim L^ψ$. The tunneling exponent $ψ$ is equal to the spin stiffness exponent $θ_c$ characterizing the zero-temperature fixed point.

cond-mat.stat-mechCross-listedarXiv:2603.02308•Mar 4, 2026•0 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

Critical scaling of the mutual information in two-dimensional disordered Ising models

Sriluckshmy, P., Mandal, I.

Renyi mutual information, computed from second Renyi entropies, can identify classical phase transitions from their finite-size scaling at critical points. We apply this technique to examine the presence or absence of finite temperature phase transitions in various two-dimensional models on a square lattice, which are extensions of the conventional Ising model by adding a quenched disorder. When the quenched disorder causes the nearest neighbor bonds to be both ferromagnetic and antiferromagnetic, (a) a spin glass phase exists only at zero temperature, and (b) a ferromagnetic phase exists at a finite temperature when the antiferromagnetic bond distributions are sufficiently dilute. Furthermore, finite temperature paramagnetic-ferromagnetic transitions can also occur when the disordered bonds involve only ferromagnetic couplings of random strengths. In our numerical simulations, the 'zero temperature only' phase transitions are identified when there is no consistent finite-size scaling of the Renyi mutual information curves, while for finite temperature critical points, the curves can identify the critical temperature T-c by their crossings at T-c and 2 T-c.

•Apr 3, 2018•8 citations

Domain-wall energies and magnetization of the two-dimensional random-bond Ising model

Carlo Amoruso, Alexander K. Hartmann

We study ground-state properties of the two-dimensional random-bond Ising model with couplings having a concentration $p∊[0,1]$ of antiferromagnetic and $(1\ensuremath{-}p)$ of ferromagnetic bonds. We apply an exact matching algorithm which enables us the study of systems with linear dimension $L$ up to $700$. We study the behavior of the domain-wall energies and of the magnetization. We find that the paramagnet-ferromagnet transition occurs at ${p}_{c}\ensuremath{\sim}0.103$ compared to the concentration ${p}_{n}\ensuremath{\sim}0.109$ at the Nishimori point, which means that the phase diagram of the model exhibits a re-entrance. Furthermore, we find no indications for an (intermediate) spin-glass ordering at finite temperature.

arXiv:cond-mat/0401464•Physical Review B•Oct 29, 2004•48 citations

Phase transitions in the two-dimensional XY model with random phases: A Monte Carlo study

J. Maucourt, Daniel R. Grempel

We study the two-dimensional $\mathrm{XY}$ model with quenched random phases by Monte Carlo simulation and finite-size scaling analysis. We determine the phase diagram of the model and study its critical behavior as a function of disorder and temperature. If the strength of the randomness is less than a critical value, ${\ensuremath{\sigma}}_{c},$ the system has a Kosterlitz-Thouless (KT) phase transition from the paramagnetic phase to a state with quasi-long-range order. Our data suggest that the latter exists down to $T=0$ in contradiction with theories that predict the appearance of a low-temperature reentrant phase. At the critical disorder ${T}_{\mathrm{KT}}\ensuremath{\rightarrow}0$ and for $\ensuremath{\sigma}>{\ensuremath{\sigma}}_{c}$ there is no quasiordered phase. At zero temperature there is a phase transition between two different glassy states at ${\ensuremath{\sigma}}_{c}.$ The functional dependence of the correlation length on $\ensuremath{\sigma}$ suggests that this transition corresponds to the disorder-driven unbinding of vortex pairs.

arXiv:cond-mat/9703109•Physical review. B, Condensed matter•Aug 1, 1997•30 citations
Systems with quenched disorder possess complex energy landscapes that are challenging to explore under the conventional Monte Carlo method. In this work, we implement an efficient entropy sampling scheme for accurate computation of the entropy function in low-energy regions. The method is applied to the two-dimensional ±J\pm J±J random-bond Ising model, where frustration is controlled by the fraction ppp of ferromagnetic bonds. We investigate the low-temperature paramagnetic--ferromagnetic phase boundary below the multicritical point at TN=0.9530(4)T_N = 0.9530(4)TN​=0.9530(4), PN=0.89078(8)P_N = 0.89078(8)PN​=0.89078(8), as well as the zero-temperature ferromagnetic--spin-glass transition. Finite-size scaling analysis reveals that the phase boundary for T<TNT < T_NT<TN​ exhibits reentrant behavior. By analyzing the evolution of the magnetization-resolved density of states g(E,M)g(E, M)g(E,M) and ground-state spin configurations against increasing frustration, we provide strong evidence that the zero-temperature transition is quasi-first order. Finite-size scaling conducted on the spin-glass side supports the validity of β=0β= 0β=0, with a correlation length exponent ν=1.50(8)ν= 1.50(8)ν=1.50(8). Our results provide new insights into the nature of the ferromagnetic-to-spin-glass phase transition in an extensively degenerate ground state.
Abstract Systems with quenched disorder possess complex energy landscapes that are challenging to explore under conventional Monte Carlo methods. In this work, we implement an efficient entropy sampling scheme for accurate computation of the entropy function in low-energy regions. The method is applied to the two-dimensional ± J random-bond Ising model, where frustration is controlled by the fraction p of ferromagnetic bonds. We investigate the low-temperature paramagnetic–ferromagnetic phase boundary below the multicritical point at T N = 0.9530(4), P N = 0.89078(8), as well as the zero-temperature ferromagnetic–spin-glass transition. Finite-size scaling analysis reveals that the phase boundary for T < T N exhibits reentrant behavior. By analyzing the evolution of the magnetization-resolved density of states g ( E , M ) and ground-state spin configurations against increasing frustration, we provide strong evidence that the zero-temperature transition is a mixed-order. Finite-size scaling conducted on the spin-glass side supports the validity of β = 0, where β is the magnetization exponent, with a correlation length exponent ν = 1.50(8). Our results provide new insights into the nature of the ferromagnetic-to-spin-glass phase transition in an extensively degenerate ground state.
At low temperatures, the classical two-dimensional random bond Ising model undergoes a frustration-driven ferromagnet-to-paramagnet transition controlled by a zero-temperature fixed point separating ferromagnet and spin glass phases. We show that this critical point can be understood through a renormalization group transformation that constructs the ground state of the Ising model through a sequence of Hamiltonians that, starting with an unfrustrated model, iteratively adds in frustration until the target Hamiltonian is reached. Via a mapping of the thermodynamics of the 2d Ising model to the spectral properties of a related Hermitian matrix -- the Hamiltonian of a noninteracting quantum problem -- this RG procedure corresponds to an iterative diagonalization of the quantum Hamiltonian. The flow toward zero temperature in the Ising picture manifests as a flow toward infinite randomness in the spectrum of the quantum Hamiltonian, with the log gap of the Hamiltonian scaling as a power of the system size: log⁡εmin−1∼Lψ\log \varepsilon_{\it min}^{-1} \sim L^ψlogεmin−1​∼Lψ. The tunneling exponent ψψψ is equal to the spin stiffness exponent θcθ_cθc​ characterizing the zero-temperature fixed point.
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.
Renyi mutual information, computed from second Renyi entropies, can identify classical phase transitions from their finite-size scaling at critical points. We apply this technique to examine the presence or absence of finite temperature phase transitions in various two-dimensional models on a square lattice, which are extensions of the conventional Ising model by adding a quenched disorder. When the quenched disorder causes the nearest neighbor bonds to be both ferromagnetic and antiferromagnetic, (a) a spin glass phase exists only at zero temperature, and (b) a ferromagnetic phase exists at a finite temperature when the antiferromagnetic bond distributions are sufficiently dilute. Furthermore, finite temperature paramagnetic-ferromagnetic transitions can also occur when the disordered bonds involve only ferromagnetic couplings of random strengths. In our numerical simulations, the 'zero temperature only' phase transitions are identified when there is no consistent finite-size scaling of the Renyi mutual information curves, while for finite temperature critical points, the curves can identify the critical temperature T-c by their crossings at T-c and 2 T-c.
We study ground-state properties of the two-dimensional random-bond Ising model with couplings having a concentration p∊[0,1]p∊[0,1]p∊[0,1] of antiferromagnetic and $(1\ensuremath{-}p)$ of ferromagnetic bonds. We apply an exact matching algorithm which enables us the study of systems with linear dimension LLL up to 700700700. We study the behavior of the domain-wall energies and of the magnetization. We find that the paramagnet-ferromagnet transition occurs at ${p}_{c}\ensuremath{\sim}0.103$ compared to the concentration ${p}_{n}\ensuremath{\sim}0.109$ at the Nishimori point, which means that the phase diagram of the model exhibits a re-entrance. Furthermore, we find no indications for an (intermediate) spin-glass ordering at finite temperature.
We study the two-dimensional XY\mathrm{XY}XY model with quenched random phases by Monte Carlo simulation and finite-size scaling analysis. We determine the phase diagram of the model and study its critical behavior as a function of disorder and temperature. If the strength of the randomness is less than a critical value, ${\ensuremath{\sigma}}_{c},$ the system has a Kosterlitz-Thouless (KT) phase transition from the paramagnetic phase to a state with quasi-long-range order. Our data suggest that the latter exists down to T=0T=0T=0 in contradiction with theories that predict the appearance of a low-temperature reentrant phase. At the critical disorder ${T}_{\mathrm{KT}}\ensuremath{\rightarrow}0$ and for $\ensuremath{\sigma}>{\ensuremath{\sigma}}_{c}$ there is no quasiordered phase. At zero temperature there is a phase transition between two different glassy states at ${\ensuremath{\sigma}}_{c}.$ The functional dependence of the correlation length on $\ensuremath{\sigma}$ suggests that this transition corresponds to the disorder-driven unbinding of vortex pairs.