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

Communications in Theoretical Physics

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, Dao‐Xin Yao, Lei‐Han Tang

Published Jun 14, 2025• 1 citation
Publisher

Abstract

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.

DOI
10.1088/1572-9494/ade49d

Related papers

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 +2 more

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.

arXiv:2505.04298•May 7, 2025•0 citations

Ferromagnetic-glassy transitions in three-dimensional Ising spin glasses

Giacomo Ceccarelli, Andrea Pelissetto C, Ettore Vicari

We investigate the ferromagnetic-glassy transitions that separate the low-temperature ferromagnetic and spin-glass phases in the temperature-disorder phase diagram of three-dimensional Ising spin-glass models. For this purpose, we consider the cubic-lattice $\ifmmode\pm\else\textpm\fi{}J$ (Edwards-Anderson) Ising model with bond distribution $P(J)=p\ensuremath{\delta}(J\ensuremath{-}1)+(1\ensuremath{-}p)\ensuremath{\delta}(J+1)$, and present a numerical Monte Carlo study of the critical behavior along the line that marks the onset of ferromagnetism. The finite-size scaling analysis of the Monte Carlo data shows that the ferromagnetic-glassy transition line is slightly reentrant. As a consequence, for an interval of the disorder parameter $p$, around $p\ensuremath{\approx}0.77$, the system presents a low-temperature glassy phase, an intermediate ferromagnetic phase, and a high-temperature paramagnetic phase. Along the ferromagnetic-glassy transition line, magnetic correlations show a universal critical behavior with critical exponents $\ensuremath{\nu}=0.96(2)$ and $\ensuremath{\eta}=\ensuremath{-}0.39(2)$. The hyperscaling relation $\ensuremath{\beta}/\ensuremath{\nu}=(1+\ensuremath{\eta})/2$ is satisfied at the transitions, so that $\ensuremath{\beta}/\ensuremath{\nu}=0.305(10)$. This magnetic critical behavior represents a new universality class for ferromagnetic transitions in Ising-type disordered systems. Overlap correlations are apparently not critical and show a smooth behavior across the transition.

arXiv:1107.3005•Physical Review B•Oct 7, 2011•29 citations

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

P V Sriluckshmy, Ipsita Mandal

Abstract Rényi mutual information, computed from second Rényi 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 Rényi 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 .

arXiv:1711.02352•Journal of Statistical Mechanics Theory and Experiment•Apr 1, 2018•9 citations

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

In and beyond the Griffiths phase: A large-deviation study of the magnetic susceptibility of the two-dimensional bond-diluted Ising model

Lambert Münster, Alexander K. Hartmann, Martin Weigel

The Griffiths phase in systems with quenched disorder occurs below the ordering transition of the pure system down to the ordering transition of the actual disordered system. While it does not exhibit long-range order, large fluctuations in the disorder degrees of freedom result in exponentially rare, long-range ordered states and hence the occurrence of broad distributions in response functions. Inside the Griffiths phase of the two-dimensional bond-diluted Ising model the distribution of the magnetic susceptibility is expected to have such a broad, exponential tail. A large-deviation Monte Carlo algorithm is used to sample this distribution and the exponential tail is extracted over a wide range of the support down to very small probabilities of the order of 10^{-300}. We study the behavior of the susceptibility distribution across the full phase diagram, from the paramagnetic state through the Griffiths phase to the ferromagnetically ordered system and down to the zero-temperature point. We extract the rate function of large-deviation theory as well as its finite-size scaling behavior and we reveal interesting differences and similarities between the cases. A connection between the fraction of ferromagnetic bonds in a given disorder sample and the size of the magnetic susceptibility is demonstrated numerically.

arXiv:2405.02889•Physical review. E•Nov 13, 2024•5 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
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.
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.
We investigate the ferromagnetic-glassy transitions that separate the low-temperature ferromagnetic and spin-glass phases in the temperature-disorder phase diagram of three-dimensional Ising spin-glass models. For this purpose, we consider the cubic-lattice $\ifmmode\pm\else\textpm\fi{}J$ (Edwards-Anderson) Ising model with bond distribution $P(J)=p\ensuremath{\delta}(J\ensuremath{-}1)+(1\ensuremath{-}p)\ensuremath{\delta}(J+1)$, and present a numerical Monte Carlo study of the critical behavior along the line that marks the onset of ferromagnetism. The finite-size scaling analysis of the Monte Carlo data shows that the ferromagnetic-glassy transition line is slightly reentrant. As a consequence, for an interval of the disorder parameter ppp, around $p\ensuremath{\approx}0.77$, the system presents a low-temperature glassy phase, an intermediate ferromagnetic phase, and a high-temperature paramagnetic phase. Along the ferromagnetic-glassy transition line, magnetic correlations show a universal critical behavior with critical exponents $\ensuremath{\nu}=0.96(2)$ and $\ensuremath{\eta}=\ensuremath{-}0.39(2)$. The hyperscaling relation $\ensuremath{\beta}/\ensuremath{\nu}=(1+\ensuremath{\eta})/2$ is satisfied at the transitions, so that $\ensuremath{\beta}/\ensuremath{\nu}=0.305(10)$. This magnetic critical behavior represents a new universality class for ferromagnetic transitions in Ising-type disordered systems. Overlap correlations are apparently not critical and show a smooth behavior across the transition.
Abstract Rényi mutual information, computed from second Rényi 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 Rényi 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 .
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.
The Griffiths phase in systems with quenched disorder occurs below the ordering transition of the pure system down to the ordering transition of the actual disordered system. While it does not exhibit long-range order, large fluctuations in the disorder degrees of freedom result in exponentially rare, long-range ordered states and hence the occurrence of broad distributions in response functions. Inside the Griffiths phase of the two-dimensional bond-diluted Ising model the distribution of the magnetic susceptibility is expected to have such a broad, exponential tail. A large-deviation Monte Carlo algorithm is used to sample this distribution and the exponential tail is extracted over a wide range of the support down to very small probabilities of the order of 10^{-300}. We study the behavior of the susceptibility distribution across the full phase diagram, from the paramagnetic state through the Griffiths phase to the ferromagnetically ordered system and down to the zero-temperature point. We extract the rate function of large-deviation theory as well as its finite-size scaling behavior and we reveal interesting differences and similarities between the cases. A connection between the fraction of ferromagnetic bonds in a given disorder sample and the size of the magnetic susceptibility is demonstrated numerically.
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.