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

Complex Replica Zeros of ±J Ising Spin Glass at Zero Temperature

Tomoyuki Obuchi, Yoshiyuki Kabashima

Published Aug 3, 2013• 7 citations
OpenAlex

Abstract

Abstract. Zeros of the nth moment of the partition function [Z n] are investigated in a vanishing temperature limit β → ∞, n → 0 keeping y = βn ∼ O(1). In this limit, the moment parameterized by y characterizes the distribution of the ground-state energy. We numerically investigate the zeros for ±J Ising spin glass models with several ladder and tree systems, which can be carried out with a feasible computational cost by a symbolic operation based on the Bethe–Peierls method. For several tree systems we find that the zeros tend to approach the real axis of y in the thermodynamic limit implying that the moment cannot be described by a single analytic function of y as the system size tends to infinity, which may be associated with breaking of the replica symmetry. However, examination of the analytical properties of the moment function and assessment of the spin-glass susceptibility indicate that the breaking of analyticity is relevant to neither one-step or full replica symmetry breaking. § obuchi@stat.phys.titech.ac.jpComplex Replica Zeros of ±J Ising Spin Glass at Zero Temperature 2 1.

Related papers

Replica symmetry breaking in the transverse-field Ising spin-glass model: Two fermionic representations

Eduardo Muccillo Monteiro Santos, Alba Theumann

We analyze the infinite range Ising spin glass in a transverse-field $\ensuremath{\Gamma}$ below the critical temperature by a one step replica symmetry breaking theory. The set of $n$ replicas is divided in $r$ blocks of $m$ replicas each. We present results for different values of the block-size parameter $m$. The spin operators are represented by bilinear combinations of fermionic fields and we compare the results of two models: In the four-states (4-S) model the diagonal ${S}_{i}^{z}$ operator has two unphysical vanishing eigenvalues, that are suppressed by a restraint in the two-states (2-S) model. In the static approximation we obtain qualitatively similar results for both models. They both exhibit a critical temperature ${T}_{c}(\ensuremath{\Gamma})$ that decreases when $\ensuremath{\Gamma}$ increases, until it reaches a quantum critical point at the same value of ${\ensuremath{\Gamma}}_{c}$ and they are both unstable under replica symmetry breaking in the whole spin glass phase. Below the critical temperature we present results for the order parameters and free energy.

arXiv:cond-mat/0505522•Physical Review B•Jul 19, 2005•4 citations

Computation of the Zero Temperature RSB parameter in Bethe Lattice Spin Glasses

F. De Santis, Giorgio Parisi A

Bethe Lattice Spin Glasses (BLSG) are models with finite connectivity which undergo a Replica Symmetry Breaking (RSB) phase transition in field, even at zero temperature. The order parameter describing this phase is a function related to the non-trivial organization of the states. We compute numerically the spin glass order parameter at zero temperature in BLSG with coordination number $z=3$ and couplings $J=\pm1$, near the transition. The method is based on a universal formula which relates the order parameter to the joint probability distribution of the energy gaps of the lowest lying states. Such distribution is derived numerically by introducing a convenient bulk perturbation to the Hamiltonian.

arXiv:1805.01228•May 3, 2018•0 citations

The Phase Transitions in a ppp spin Glass Model: A Numerical Study

Prerak Gupta, Auditya Sharma, Bharadwaj Vedula +2 more

We investigate the balanced $M=4$, $p=4$ spin-glass model for a one-dimensional long-range proxy for the finite dimensional short-range $p$-spin glass model to examine the nature of the glass transition beyond mean-field theory. We perform large-scale Monte Carlo equilibrated simulations for both fully connected and power-law diluted versions of the model. The critical temperatures extracted from the finite-size scaling (FSS) analysis of spin-glass susceptibility are in good agreement with theoretical predictions for $σ= 0, 0.25$, and 0.55. For these values of the long-range exponent $σ$ (which is the power of the decrease of the interactions between the spins with their separation), one might have expected that mean-field theory would provide a good description of the system. However, the spin-overlap distribution and the value of the $λ$-parameter do not provide numerical evidence for a one-step replica symmetry breaking (1RSB) phase transition. Instead, our results indicate a direct transition from the paramagnetic state to a full replica symmetry broken phase, with a renormalized value of $λ\equiv ω_2/ω_1 < 1$ suggesting a continuous FRSB transition, despite this ratio being equal to 2 at mean-field level. A value of $λ> 1$ is required for the discontinuous 1RSB transition. We argue that strong finite-size effects and closely spaced transition temperatures remove the expected 1RSB transition for the system sizes which we can study. For values of the exponent $σ= 0.85$, which roughly corresponds to a three dimensional system, we find that the renormalized value of $λ$ is again less than 1, with no signs of either the 1RSB transition or the continuous FRSB transition, suggesting that the Kauzmann temperature $T_K$ in three dimensions might be zero and the complete absence of phase transitions in structural glasses.

arXiv:2604.15433•Apr 16, 2026•0 citations

Low Temperature Phase of the m-component Spin Glass

A. Braun

We investigate the low temperature phase of a vector spin glass in the limit of large component numbers $m$ of the spins, where this model is expected to be replica symmetric. On a Bethe lattice, a random graph with fixed and finite connectivity, we employ a self-consistent effective-field approach to the cavity method to describe equilibrium states of the system in terms of the configuration of cavity fields, which are reduced mean fields onto each spin. We make use of an iterative algorithm to find such equilibrium states for finite system sizes $N$ and observe a phase transition, which is indicated by an order parameter describing the preferred direction each spin points into in terms of the configuration of cavity, resp. mean local fields. On single realizations of Bethe lattices the transition is sharp and accompanied by the Generalized-Bose-Einstein condensation, in which the cavity, resp. local fields of the spins condense into a subspace of the available $m$-dimensional space for large $m$. The dimension of this subspace, $m_{\rm eff}(T)$ depends on the temperature. Is either one or two right at the critical temperature and increases below. Furthermore, for fixed temperature it is expected to scale with the system size $N$ with an exponent $\mu_m$ which depends on the connectivity. Averaging over many realizations of disorder we find $\overline{m_{\rm eff}}(T_c)=2$ and $\overline{m_{\rm eff}}(T)\sim N^{\mu_m}$. Even more, the cavity fields are subject to rotations in the low temperature phase, i.e. there is no freezing in time as generally expected. In case of $m_{\rm eff}(T_c)=2$ the whole configuration rotates at the critical temperature about the same angle and the same axis of rotation corresponding to a single droplet. In addition, we find replica symmetry breaking below the large-$m$ limit, i.e. for $m_{\rm eff}(T) > m$. These results are supported by results of parallel tempering Monte Carlo simulations. Replica symmetry breaking is further investigated by an analytical method in the large connectivity limit, where we compare the free energy in terms of spin-spin-correlations of different replicas calculated with a replica symmetric and a one-step broken replica symmetry ansatz. In the large-$m$ limit, we find that both solutions coexist for finite connectivities and that replica symmetry holds for full connectivity up to corrections in second order around the saddle point for $m\to\infty$. In the last chapter we calculate the sample-to-sample fluctuations of the free energy of the SK model of the large-$m$ limit. We derive an exact connection between these fluctuations and bond chaos. Making use of a replica calculation for the {\it large deviations} and a scaling ansatz for the {\it small deviations}, we quantify bond chaos, i.e. calculate the probability of finding a given link overlap for two systems with possibly different sets of bonds. With this, we find the sample-to-sample fluctuations to scale with the number of spins, $N$, with an exponent $\mu$ in between $1/5$ and $3/10$.

•Jan 1, 2011•0 citations

Quantum Critical Behavior of the Infinite-range Transverse Ising Spin Glass : An Exact Numerical Diagonalization Study

Parongama Sen, Purusattam Ray, Bikas K. Chakrabarti

We report exact numerical diagonalization results of the infinite-range Ising spin glass in a transverse field $Γ$ at zero temperature. Eigenvalues and eigenvectors are determined for various strengths of $Γ$ and for system sizes $N \le 16$. We obtain the moments of the distribution of the spin-glass order parameter, the spin-glass susceptibility and the mass gap at different values of $Γ$. The disorder averaging is done typically over 1000 configurations. Our finite size scaling analysis indicates a spin glass transition at $Γ_c \simeq 1.5$. Our estimates for the exponents at the transition are in agreement with those known from other approaches. For the dynamic exponent, we get $z=2.1 \pm 0.1$ which is in contradiction with a recent estimate ($z=4$). Our cumulant analysis indicates the existence of a replica symmetric spin glass phase for $Γ< Γ_c$.

arXiv:cond-mat/9705297•May 29, 1997•0 citations

Analytic and computational studies of spin glasses

Shreyans Jain

The Sherrington-Kirkpatrick and the nearest-neighbour X-Y spin-glass models are studied.The effect of a small uniform symmetry breaking field on the infinite-ranged Sherrington-Kirkpatrick Ising spin-glass model is discussed.It is shown that the off-diagonal contributions to the fluctuation-dissipation formula vanish for a symmetric distribution of interactions and the limit of zero field.This contradicts the result of Young and Kirkpatrick based on numerical evidence for small samples.The dynamics of the Sherrington-Kirkpatrick spin-glass are investigated.Below the critical temperature this model exhibits non-ergodic behaviour.Consequently, in order to obtain timedependent correlation functions it is necessary to average over the initial states with a Boltzmann weight.This necessitates the use of "replicas".It is shown that in Parisi's scheme the "dynamical" Edwards-Anderson order parameter is given by q(x =1), whereas the statistical mechanics order parameter is given by r 1 q(x)dx.The result for the statistical mechanics order parameter j 0 implies that there is no violation of the fluctuation-dissipation theorem once an average over samples is performed.It is also shown that the system in general keeps its memory of the initial state for all time.However, this memory is "lost" if one makes a replica symmetric ansatz for the order parameters.This illus trates the connection between replica symmetry breaking and nonergodicity.The nearest neighbour X-Y spin-glass in both two-and threedimensions is investigated by Monte Carlo simulation.In two--11ln[ln Xcn + 1-42] versus In T SG } X2 (Tq ) for T (3D) 4.11 4.12a 4.12b 4.12c 4.13a ' 4.13b > Plots of the data for T (3D) 4.13c 4.14a 4.14b 4.15a 4.15b Plots of £ ^ for 3D SG A plot of the data for £ assuming that LCD = 3 SG A plot of the data for £ assuming that LCD = 3 SG and p__ = 1.2 SG A5.1a A5.

•Jan 1, 1985•0 citations
Abstract. Zeros of the nth moment of the partition function [Z n] are investigated in a vanishing temperature limit β → ∞, n → 0 keeping y = βn ∼ O(1). In this limit, the moment parameterized by y characterizes the distribution of the ground-state energy. We numerically investigate the zeros for ±J Ising spin glass models with several ladder and tree systems, which can be carried out with a feasible computational cost by a symbolic operation based on the Bethe–Peierls method. For several tree systems we find that the zeros tend to approach the real axis of y in the thermodynamic limit implying that the moment cannot be described by a single analytic function of y as the system size tends to infinity, which may be associated with breaking of the replica symmetry. However, examination of the analytical properties of the moment function and assessment of the spin-glass susceptibility indicate that the breaking of analyticity is relevant to neither one-step or full replica symmetry breaking. § obuchi@stat.phys.titech.ac.jpComplex Replica Zeros of ±J Ising Spin Glass at Zero Temperature 2 1.
We analyze the infinite range Ising spin glass in a transverse-field $\ensuremath{\Gamma}$ below the critical temperature by a one step replica symmetry breaking theory. The set of nnn replicas is divided in rrr blocks of mmm replicas each. We present results for different values of the block-size parameter mmm. The spin operators are represented by bilinear combinations of fermionic fields and we compare the results of two models: In the four-states (4-S) model the diagonal Siz{S}_{i}^{z}Siz​ operator has two unphysical vanishing eigenvalues, that are suppressed by a restraint in the two-states (2-S) model. In the static approximation we obtain qualitatively similar results for both models. They both exhibit a critical temperature ${T}_{c}(\ensuremath{\Gamma})$ that decreases when $\ensuremath{\Gamma}$ increases, until it reaches a quantum critical point at the same value of ${\ensuremath{\Gamma}}_{c}$ and they are both unstable under replica symmetry breaking in the whole spin glass phase. Below the critical temperature we present results for the order parameters and free energy.
Bethe Lattice Spin Glasses (BLSG) are models with finite connectivity which undergo a Replica Symmetry Breaking (RSB) phase transition in field, even at zero temperature. The order parameter describing this phase is a function related to the non-trivial organization of the states. We compute numerically the spin glass order parameter at zero temperature in BLSG with coordination number z=3z=3z=3 and couplings J=±1J=\pm1J=±1, near the transition. The method is based on a universal formula which relates the order parameter to the joint probability distribution of the energy gaps of the lowest lying states. Such distribution is derived numerically by introducing a convenient bulk perturbation to the Hamiltonian.
We investigate the balanced M=4M=4M=4, p=4p=4p=4 spin-glass model for a one-dimensional long-range proxy for the finite dimensional short-range ppp-spin glass model to examine the nature of the glass transition beyond mean-field theory. We perform large-scale Monte Carlo equilibrated simulations for both fully connected and power-law diluted versions of the model. The critical temperatures extracted from the finite-size scaling (FSS) analysis of spin-glass susceptibility are in good agreement with theoretical predictions for σ=0,0.25σ= 0, 0.25σ=0,0.25, and 0.55. For these values of the long-range exponent σσσ (which is the power of the decrease of the interactions between the spins with their separation), one might have expected that mean-field theory would provide a good description of the system. However, the spin-overlap distribution and the value of the λλλ-parameter do not provide numerical evidence for a one-step replica symmetry breaking (1RSB) phase transition. Instead, our results indicate a direct transition from the paramagnetic state to a full replica symmetry broken phase, with a renormalized value of λ≡ω2/ω1<1λ\equiv ω_2/ω_1 < 1λ≡ω2​/ω1​<1 suggesting a continuous FRSB transition, despite this ratio being equal to 2 at mean-field level. A value of λ>1λ> 1λ>1 is required for the discontinuous 1RSB transition. We argue that strong finite-size effects and closely spaced transition temperatures remove the expected 1RSB transition for the system sizes which we can study. For values of the exponent σ=0.85σ= 0.85σ=0.85, which roughly corresponds to a three dimensional system, we find that the renormalized value of λλλ is again less than 1, with no signs of either the 1RSB transition or the continuous FRSB transition, suggesting that the Kauzmann temperature TKT_KTK​ in three dimensions might be zero and the complete absence of phase transitions in structural glasses.
We investigate the low temperature phase of a vector spin glass in the limit of large component numbers mmm of the spins, where this model is expected to be replica symmetric. On a Bethe lattice, a random graph with fixed and finite connectivity, we employ a self-consistent effective-field approach to the cavity method to describe equilibrium states of the system in terms of the configuration of cavity fields, which are reduced mean fields onto each spin. We make use of an iterative algorithm to find such equilibrium states for finite system sizes NNN and observe a phase transition, which is indicated by an order parameter describing the preferred direction each spin points into in terms of the configuration of cavity, resp. mean local fields. On single realizations of Bethe lattices the transition is sharp and accompanied by the Generalized-Bose-Einstein condensation, in which the cavity, resp. local fields of the spins condense into a subspace of the available mmm-dimensional space for large mmm. The dimension of this subspace, meff(T)m_{\rm eff}(T)meff​(T) depends on the temperature. Is either one or two right at the critical temperature and increases below. Furthermore, for fixed temperature it is expected to scale with the system size NNN with an exponent μm\mu_mμm​ which depends on the connectivity. Averaging over many realizations of disorder we find meff‾(Tc)=2\overline{m_{\rm eff}}(T_c)=2meff​​(Tc​)=2 and meff‾(T)∼Nμm\overline{m_{\rm eff}}(T)\sim N^{\mu_m}meff​​(T)∼Nμm​. Even more, the cavity fields are subject to rotations in the low temperature phase, i.e. there is no freezing in time as generally expected. In case of meff(Tc)=2m_{\rm eff}(T_c)=2meff​(Tc​)=2 the whole configuration rotates at the critical temperature about the same angle and the same axis of rotation corresponding to a single droplet. In addition, we find replica symmetry breaking below the large-mmm limit, i.e. for meff(T)>mm_{\rm eff}(T) > mmeff​(T)>m. These results are supported by results of parallel tempering Monte Carlo simulations. Replica symmetry breaking is further investigated by an analytical method in the large connectivity limit, where we compare the free energy in terms of spin-spin-correlations of different replicas calculated with a replica symmetric and a one-step broken replica symmetry ansatz. In the large-mmm limit, we find that both solutions coexist for finite connectivities and that replica symmetry holds for full connectivity up to corrections in second order around the saddle point for m→∞m\to\inftym→∞. In the last chapter we calculate the sample-to-sample fluctuations of the free energy of the SK model of the large-mmm limit. We derive an exact connection between these fluctuations and bond chaos. Making use of a replica calculation for the {\it large deviations} and a scaling ansatz for the {\it small deviations}, we quantify bond chaos, i.e. calculate the probability of finding a given link overlap for two systems with possibly different sets of bonds. With this, we find the sample-to-sample fluctuations to scale with the number of spins, NNN, with an exponent μ\muμ in between 1/51/51/5 and 3/103/103/10.
We report exact numerical diagonalization results of the infinite-range Ising spin glass in a transverse field ΓΓΓ at zero temperature. Eigenvalues and eigenvectors are determined for various strengths of ΓΓΓ and for system sizes N≤16N \le 16N≤16. We obtain the moments of the distribution of the spin-glass order parameter, the spin-glass susceptibility and the mass gap at different values of ΓΓΓ. The disorder averaging is done typically over 1000 configurations. Our finite size scaling analysis indicates a spin glass transition at Γc≃1.5Γ_c \simeq 1.5Γc​≃1.5. Our estimates for the exponents at the transition are in agreement with those known from other approaches. For the dynamic exponent, we get z=2.1±0.1z=2.1 \pm 0.1z=2.1±0.1 which is in contradiction with a recent estimate (z=4z=4z=4). Our cumulant analysis indicates the existence of a replica symmetric spin glass phase for Γ<ΓcΓ< Γ_cΓ<Γc​.
The Sherrington-Kirkpatrick and the nearest-neighbour X-Y spin-glass models are studied.The effect of a small uniform symmetry breaking field on the infinite-ranged Sherrington-Kirkpatrick Ising spin-glass model is discussed.It is shown that the off-diagonal contributions to the fluctuation-dissipation formula vanish for a symmetric distribution of interactions and the limit of zero field.This contradicts the result of Young and Kirkpatrick based on numerical evidence for small samples.The dynamics of the Sherrington-Kirkpatrick spin-glass are investigated.Below the critical temperature this model exhibits non-ergodic behaviour.Consequently, in order to obtain timedependent correlation functions it is necessary to average over the initial states with a Boltzmann weight.This necessitates the use of "replicas".It is shown that in Parisi's scheme the "dynamical" Edwards-Anderson order parameter is given by q(x =1), whereas the statistical mechanics order parameter is given by r 1 q(x)dx.The result for the statistical mechanics order parameter j 0 implies that there is no violation of the fluctuation-dissipation theorem once an average over samples is performed.It is also shown that the system in general keeps its memory of the initial state for all time.However, this memory is "lost" if one makes a replica symmetric ansatz for the order parameters.This illus trates the connection between replica symmetry breaking and nonergodicity.The nearest neighbour X-Y spin-glass in both two-and threedimensions is investigated by Monte Carlo simulation.In two--11ln[ln Xcn + 1-42] versus In T SG } X2 (Tq ) for T (3D) 4.11 4.12a 4.12b 4.12c 4.13a ' 4.13b > Plots of the data for T (3D) 4.13c 4.14a 4.14b 4.15a 4.15b Plots of £ ^ for 3D SG A plot of the data for £ assuming that LCD = 3 SG A plot of the data for £ assuming that LCD = 3 SG and p__ = 1.2 SG A5.1a A5.