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:0811.1626 · Journal of the Physical Society of Japan

Ground-State Properties of a Heisenberg Spin Glass Model with a Hybrid Genetic Algorithm

Yuh-ichi Iyama, Fumitaka Matsubara

Published Jan 13, 2009• 1 citation
arXivHTMLPDF

Abstract

We developed a genetic algorithm (GA) in the Heisenberg model that combines a triadic crossover and a parameter-free genetic algorithm. Using the algorithm, we examined the ground-state stiffness of the $\pm J$ Heisenberg model in three dimensions up to a moderate size range. Results showed the stiffness constant of $θ= 0$ in the periodic-antiperiodic boundary condition method and that of $θ\sim 0.62$ in the open-boundary-twist method. We considered the origin of the difference in $θ$ between the two methods and suggested that both results show the same thing: the ground state of the open system is stable against a weak perturbation.

DOI
10.1143/jpsj.78.014703

Related papers

Ground state and low-lying excitations in a Heisenberg spin-glass model in three dimensions

Fumitaka Matsubara, Takayuki Shirakura, S. Takahashi +1 more

We developed a hybrid genetic algorithm and studied the ground state and low-lying excitations of the $\ifmmode\pm\else\textpm\fi{}J$ Heisenberg model on a three-dimensional $L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L(L\ensuremath{\leqslant}11)$ lattice. Our results suggest occurrence of a Parisi state. That is, there exists a metastable state with a finite energy of $\ensuremath{\Delta}E\ensuremath{\sim}O(J)$ in the thermodynamic limit and the energy barrier $\ensuremath{\Delta}W$ between the metastable state and the ground state is $\ensuremath{\Delta}W\ensuremath{\propto}J{L}^{\ensuremath{\theta}}$ with $\ensuremath{\theta}\ensuremath{\sim}0.53$. Our results also suggest that the energy $\ensuremath{\Delta}E$ of small scale excitations are related to a scale of the excitation $V$ in a similar manner to those of the ferromagnetic model, i.e., $\ensuremath{\Delta}E\ensuremath{\propto}J{L}^{{\ensuremath{\theta}}_{L}}{V}^{{\ensuremath{\theta}}_{v}}$ with ${\ensuremath{\theta}}_{L}\ensuremath{\sim}\ensuremath{-}4.7$ and ${\ensuremath{\theta}}_{v}\ensuremath{\sim}2.0$. This fact suggests that similarity exists in the ground state structure between the $\ifmmode\pm\else\textpm\fi{}J$ and the ferromagnetic Heisenberg models, though the spin directions of those two models differ greatly from one another.

Physical Review B•Nov 11, 2004•7 citations

Scaling of stiffness energy for three-dimensional ±J Ising spin glasses

Alexander K. Hartmann

Large numbers of ground states of three-dimensional (3D) Edwards-Anderson Ising spin glasses are calculated for sizes up to ${10}^{3}$ using a combination of a genetic algorithm and cluster-exact approximation. A detailed analysis shows that true ground states are obtained. The ground-state stiffness (or domain wall) energy \ensuremath{\Delta} is calculated. A $|\ensuremath{\Delta}|\ensuremath{\sim}{L}^{{\ensuremath{\Theta}}_{S}}$ behavior with ${\ensuremath{\Theta}}_{S}=0.19(2)$ is found that strongly indicates that the 3D model has an equilibrium spin-glass-paramagnet transition for nonzero ${T}_{c}.$

arXiv:cond-mat/9806114•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 1, 1999•93 citations

Ground States and Defect Energies of the Two-Dimensional XY Spin Glass from a Quasiexact Algorithm

Martin Weigel, Michel J. P. Gingras

We employ a novel algorithm using a quasiexact embedded-cluster matching technique as minimization method within a genetic algorithm to reliably obtain numerically exact ground states of the Edwards-Anderson XY spin-glass model with bimodal coupling distribution for square lattices of up to 28 x 28 spins. Contrary to previous conjectures, the ground state of each disorder replica is nondegenerate up to a global O(2) rotation. The scaling of spin and chiral defect energies induced by applying several different sets of boundary conditions exhibits strong crossover effects. This suggests that previous calculations have yielded results far from the asymptotic regime. The novel algorithm and the aspect-ratio scaling technique consistently give theta(s)=-0.308(30) and theta(c)=-0.114(16) for the spin and chiral stiffness exponents, respectively.

arXiv:cond-mat/0510614•Physical Review Letters•Mar 9, 2006•16 citations

Sampling the ground-state magnetization of d -dimensional p -body Ising models

Creighton K. Thomas, Helmut G. Katzgraber

We demonstrate that a recently introduced heuristic optimization algorithm [Phys. Rev. E 83, 046709 (2011)] that combines a local search with triadic crossover genetic updates is capable of sampling nearly uniformly among ground-state configurations in spin-glass-like Hamiltonians with $p$-spin interactions in $d$ space dimensions that have highly degenerate ground states. Using this algorithm we probe the zero-temperature ferromagnet to spin-glass transition point ${q}_{c}$ of two example models, the disordered version of the two-dimensional three-spin Baxter-Wu model [${q}_{c}=0.1072(1)$] and the three-dimensional Edwards-Anderson model [${q}_{c}=0.2253(7)$], by computing the Binder ratio of the ground-state magnetization.

arXiv:1108.1953•Physical Review B•Nov 4, 2011•6 citations

Zero-temperature phase of the XY spin glass in two dimensions: Genetic embedded matching heuristic

Martin Weigel, Michel J. P. Gingras

For many real spin-glass materials, the Edwards-Anderson model with continuous-symmetry spins is more realistic than the rather better understood Ising variant. In principle, the nature of an occurring spin-glass phase in such systems might be inferred from an analysis of the zero-temperature properties. Unfortunately, with few exceptions, the problem of finding ground-state configurations is a nonpolynomial problem computationally, such that efficient approximation algorithms are called for. Here, we employ the recently developed genetic embedded matching (GEM) heuristic to investigate the nature of the zero-temperature phase of the bimodal $XY$ spin glass in two dimensions. We analyze bulk properties such as the asymptotic ground-state energy and the phase diagram of disorder strength vs disorder concentration. For the case of a symmetric distribution of ferromagnetic and antiferromagnetic bonds, we find that the ground state of the model is unique up to a global $\mathrm{O}(2)$ rotation of the spins. In particular, there are no extensive degeneracies in this model. The main focus of this work is on an investigation of the excitation spectrum as probed by changing the boundary conditions. Using appropriate finite-size scaling techniques, we consistently determine the stiffness of spin and chiral domain walls and the corresponding fractal dimensions. Most noteworthy, we find that the spin and chiral channels are characterized by two distinct stiffness exponents and, consequently, the system displays spin-chirality decoupling at large length scales. Results for the overlap distribution do not support the possibility of a multitude of thermodynamic pure states.

arXiv:0706.0227•Physical Review B•Mar 27, 2008•16 citations

Low-energy excitations in spin glasses from exact ground states

Matteo Palassini, Frauke Liers, Michael Juenger +1 more

We investigate the nature of the low-energy, large-scale excitations in the three-dimensional Edwards-Anderson Ising spin glass with Gaussian couplings and free boundary conditions, by studying the response of the ground state to a coupling-dependent perturbation introduced previously. The ground states are determined exactly for system sizes up to ${12}^{3}$ spins using a branch-and-cut algorithm. The data are consistent with a picture where the surface of the excitations is not space filling, such as the droplet or the trivial-non-trivial (TNT) picture, with only minimal corrections to scaling. When allowing for very large corrections to scaling, the data are also consistent with a picture with space-filling surfaces, such as replica symmetry breaking. The energy of the excitations scales with their size with a small exponent ${\ensuremath{\theta}}^{\ensuremath{'}},$ which is compatible with zero if we allow moderate corrections to scaling. We compare the results with data for periodic boundary conditions obtained with a genetic algorithm, and discuss the effects of different boundary conditions on corrections to scaling. Finally, we analyze the performance of our branch-and-cut algorithm, finding that it is correlated with the existence of large-scale, low-energy excitations.

arXiv:cond-mat/0212551•Physical review. B, Condensed matter•Aug 13, 2003•21 citations
We developed a genetic algorithm (GA) in the Heisenberg model that combines a triadic crossover and a parameter-free genetic algorithm. Using the algorithm, we examined the ground-state stiffness of the ±J\pm J±J Heisenberg model in three dimensions up to a moderate size range. Results showed the stiffness constant of θ=0θ= 0θ=0 in the periodic-antiperiodic boundary condition method and that of θ∼0.62θ\sim 0.62θ∼0.62 in the open-boundary-twist method. We considered the origin of the difference in θθθ between the two methods and suggested that both results show the same thing: the ground state of the open system is stable against a weak perturbation.
We developed a hybrid genetic algorithm and studied the ground state and low-lying excitations of the $\ifmmode\pm\else\textpm\fi{}J$ Heisenberg model on a three-dimensional $L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L(L\ensuremath{\leqslant}11)$ lattice. Our results suggest occurrence of a Parisi state. That is, there exists a metastable state with a finite energy of $\ensuremath{\Delta}E\ensuremath{\sim}O(J)$ in the thermodynamic limit and the energy barrier $\ensuremath{\Delta}W$ between the metastable state and the ground state is $\ensuremath{\Delta}W\ensuremath{\propto}J{L}^{\ensuremath{\theta}}$ with $\ensuremath{\theta}\ensuremath{\sim}0.53$. Our results also suggest that the energy $\ensuremath{\Delta}E$ of small scale excitations are related to a scale of the excitation VVV in a similar manner to those of the ferromagnetic model, i.e., $\ensuremath{\Delta}E\ensuremath{\propto}J{L}^{{\ensuremath{\theta}}_{L}}{V}^{{\ensuremath{\theta}}_{v}}$ with ${\ensuremath{\theta}}_{L}\ensuremath{\sim}\ensuremath{-}4.7$ and ${\ensuremath{\theta}}_{v}\ensuremath{\sim}2.0$. This fact suggests that similarity exists in the ground state structure between the $\ifmmode\pm\else\textpm\fi{}J$ and the ferromagnetic Heisenberg models, though the spin directions of those two models differ greatly from one another.
Large numbers of ground states of three-dimensional (3D) Edwards-Anderson Ising spin glasses are calculated for sizes up to 103{10}^{3}103 using a combination of a genetic algorithm and cluster-exact approximation. A detailed analysis shows that true ground states are obtained. The ground-state stiffness (or domain wall) energy \ensuremath{\Delta} is calculated. A $|\ensuremath{\Delta}|\ensuremath{\sim}{L}^{{\ensuremath{\Theta}}_{S}}$ behavior with ${\ensuremath{\Theta}}_{S}=0.19(2)$ is found that strongly indicates that the 3D model has an equilibrium spin-glass-paramagnet transition for nonzero Tc.{T}_{c}.Tc​.
We employ a novel algorithm using a quasiexact embedded-cluster matching technique as minimization method within a genetic algorithm to reliably obtain numerically exact ground states of the Edwards-Anderson XY spin-glass model with bimodal coupling distribution for square lattices of up to 28 x 28 spins. Contrary to previous conjectures, the ground state of each disorder replica is nondegenerate up to a global O(2) rotation. The scaling of spin and chiral defect energies induced by applying several different sets of boundary conditions exhibits strong crossover effects. This suggests that previous calculations have yielded results far from the asymptotic regime. The novel algorithm and the aspect-ratio scaling technique consistently give theta(s)=-0.308(30) and theta(c)=-0.114(16) for the spin and chiral stiffness exponents, respectively.
We demonstrate that a recently introduced heuristic optimization algorithm [Phys. Rev. E 83, 046709 (2011)] that combines a local search with triadic crossover genetic updates is capable of sampling nearly uniformly among ground-state configurations in spin-glass-like Hamiltonians with ppp-spin interactions in ddd space dimensions that have highly degenerate ground states. Using this algorithm we probe the zero-temperature ferromagnet to spin-glass transition point qc{q}_{c}qc​ of two example models, the disordered version of the two-dimensional three-spin Baxter-Wu model [qc=0.1072(1){q}_{c}=0.1072(1)qc​=0.1072(1)] and the three-dimensional Edwards-Anderson model [qc=0.2253(7){q}_{c}=0.2253(7)qc​=0.2253(7)], by computing the Binder ratio of the ground-state magnetization.
For many real spin-glass materials, the Edwards-Anderson model with continuous-symmetry spins is more realistic than the rather better understood Ising variant. In principle, the nature of an occurring spin-glass phase in such systems might be inferred from an analysis of the zero-temperature properties. Unfortunately, with few exceptions, the problem of finding ground-state configurations is a nonpolynomial problem computationally, such that efficient approximation algorithms are called for. Here, we employ the recently developed genetic embedded matching (GEM) heuristic to investigate the nature of the zero-temperature phase of the bimodal XYXYXY spin glass in two dimensions. We analyze bulk properties such as the asymptotic ground-state energy and the phase diagram of disorder strength vs disorder concentration. For the case of a symmetric distribution of ferromagnetic and antiferromagnetic bonds, we find that the ground state of the model is unique up to a global O(2)\mathrm{O}(2)O(2) rotation of the spins. In particular, there are no extensive degeneracies in this model. The main focus of this work is on an investigation of the excitation spectrum as probed by changing the boundary conditions. Using appropriate finite-size scaling techniques, we consistently determine the stiffness of spin and chiral domain walls and the corresponding fractal dimensions. Most noteworthy, we find that the spin and chiral channels are characterized by two distinct stiffness exponents and, consequently, the system displays spin-chirality decoupling at large length scales. Results for the overlap distribution do not support the possibility of a multitude of thermodynamic pure states.
We investigate the nature of the low-energy, large-scale excitations in the three-dimensional Edwards-Anderson Ising spin glass with Gaussian couplings and free boundary conditions, by studying the response of the ground state to a coupling-dependent perturbation introduced previously. The ground states are determined exactly for system sizes up to 123{12}^{3}123 spins using a branch-and-cut algorithm. The data are consistent with a picture where the surface of the excitations is not space filling, such as the droplet or the trivial-non-trivial (TNT) picture, with only minimal corrections to scaling. When allowing for very large corrections to scaling, the data are also consistent with a picture with space-filling surfaces, such as replica symmetry breaking. The energy of the excitations scales with their size with a small exponent ${\ensuremath{\theta}}^{\ensuremath{'}},$ which is compatible with zero if we allow moderate corrections to scaling. We compare the results with data for periodic boundary conditions obtained with a genetic algorithm, and discuss the effects of different boundary conditions on corrections to scaling. Finally, we analyze the performance of our branch-and-cut algorithm, finding that it is correlated with the existence of large-scale, low-energy excitations.