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.
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.
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 .
In this work we present a thorough analysis of the phase transitions that occur in a ferromagnetic two-dimensional Ising model, with only nearest-neighbors interactions, in the framework of the Tsallis nonextensive statistics. We performed Monte Carlo simulations on square lattices with linear sizes L ranging from 32 up to 512. The statistical weight of the Metropolis algorithm was changed according to the nonextensive statistics. Discontinuities in the m(T) curve are observed for $q\ensuremath{\le}0.5$. However, we have verified only one peak on the energy histograms at the critical temperatures, indicating the occurrence of continuous phase transitions. For the $0.5<q\ensuremath{\le}1.0$ regime, we have found continuous phase transitions between the ordered and the disordered phases, and determined the critical exponents via finite-size scaling. We verified that the critical exponents $\ensuremath{\alpha}$, $\ensuremath{\beta}$, and $\ensuremath{\gamma}$ depend on the entropic index q in the range $0.5<q\ensuremath{\le}1.0$ in the form $\ensuremath{\alpha}(q)=(10{q}^{2}\ensuremath{-}33q+23)/20$, $\ensuremath{\beta}(q)=(2q\ensuremath{-}1)/8$, and $\ensuremath{\gamma}(q)=({q}^{2}\ensuremath{-}q+7)/4$. On the other hand, the critical exponent $\ensuremath{\nu}$ does not depend on q. This suggests a violation of the scaling relations $2\ensuremath{\beta}+\ensuremath{\gamma}=d\ensuremath{\nu}$ and $\ensuremath{\alpha}+2\ensuremath{\beta}+\ensuremath{\gamma}=2$ and a nonuniversality of the critical exponents along the ferro-paramagnetic frontier.
We investigate the zero-temperature glassy transitions in the square-lattice ±J Ising model, with bond distribution P(J{xy})=pδ(J{xy}-J)+(1-p)δ(J{xy}+J) ; p=1 and p=1/2 correspond to the pure Ising model and to the Ising spin glass with symmetric bimodal distribution, respectively. We present finite-temperature Monte Carlo simulations at p=4/5 , which is close to the low-temperature paramagnetic-ferromagnetic transition line located at p≈0.89 , and at p=1/2 . Their comparison provides a strong evidence that the glassy critical behavior that occurs for 1-p{0}<p<p{0} , p{0}≈0.897 , is universal, i.e., independent of p . Moreover, we show that glassy and magnetic modes are not coupled at the multicritical zero-temperature point where the paramagnetic-ferromagnetic transition line and the T=0 glassy transition line meet. On the theoretical side we discuss the validity of finite-size scaling in glassy systems with a zero-temperature transition and a discrete Hamiltonian spectrum. Because of a freezing phenomenon which occurs in a finite volume at sufficiently low temperatures, the standard finite-size scaling limit in terms of TL(1/ν) does not exist; the renormalization-group invariant quantity ξ/L should be used instead as basic variable.
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 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 TN=0.9530(4), 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<TN 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.
The random-bond Ising model on the square lattice has several disordered critical points, depending on the probability distribution of the bonds. There are a finite-temperature multicritical point, called the Nishimori point, and a zero-temperature fixed point, for both a binary distribution where the coupling constants take the values ± J and a Gaussian disorder distribution. Inclusion of dilution in the ± J distribution ( J = 0 for some bonds) gives rise to another zero-temperature fixed point which can be identified with percolation in the non-frustrated case ( J ≥ 0). We study these fixed points using numerical (transfer matrix) methods. We determine the location, critical exponents, and central charge of the different fixed points and study the spin–spin correlation functions. Our main findings are the following. (1) We confirm that the Nishimori point is universal with respect to the type of disorder, i.e. we obtain the same central charge and critical exponents for the ± J and Gaussian distributions of disorder. (2) The Nishimori point, the zero-temperature fixed point for the ± J and Gaussian distributions of disorder, and the percolation point in the diluted case all belong to mutually distinct universality classes. (3) The paramagnetic phase is re-entrant below the Nishimori point, i.e. the zero-temperature fixed points are not located exactly below the Nishimori point, either for the ± J distribution, or for the Gaussian distribution.