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:1709.03609 · Physical review. B./Physical review. B

Efficient generation of series expansions for ±J Ising spin glasses in a classical or a quantum field

Rajiv R. P. Singh, A. P. Young

Published Dec 29, 2017• 2 citations
arXivHTMLPDF

Abstract

We discuss generation of series expansions for Ising spin glasses with a symmetric $\ifmmode\pm\else\textpm\fi{}J$ (i.e., bimodal) distribution on $d$-dimensional hypercubic lattices using linked-cluster methods. Simplifications for the bimodal distribution allow us to go to higher order than for a general distribution. We discuss two types of problems, one classical and one quantum. The classical problem is that of the Ising spin glass in a longitudinal magnetic field $h$, for which we obtain high temperature series expansions in variables $tanh(J/T)$ and $tanh(h/T)$. The quantum problem is a $T=0$ study of the Ising spin glass in a transverse magnetic field ${h}_{T}$ for which we obtain a perturbation theory in powers of $J/{h}_{T}$. These methods require (i) enumeration and counting of all connected clusters that can be embedded in the lattice up to some order $n$, and (ii) an evaluation of the contribution of each cluster for the quantity being calculated, known as the weight. We discuss a general method that takes the much smaller list (and count) of all no free-end (NFE) clusters on a lattice up to some order $n$ and automatically generates all other clusters and their counts up to the same order. The weights for finite clusters in both cases have a simple graphical interpretation that allows us to proceed efficiently for a general configuration of the $\ifmmode\pm\else\textpm\fi{}J$ bonds and at the end perform suitable disorder averaging. The order of our computations is limited by the weight calculations for the high-temperature expansions of the classical model, while they are limited by graph counting for the $T=0$ quantum system. Details of the calculational methods are presented.

DOI
10.1103/physrevb.96.214210

Related papers

Linked cluster expansion and 1/d expansions for classical spin systems

Einar Halvorsen, M. Bartkowiak

We consider $1/d$ expansions for classical spin systems based on the vertex renormalized linked cluster expansion (LCE). The free multiplicities of the LCE graphs on a hypercubic lattice in an arbitrary dimension d are calculated using generating functions. The technique is applied to the Ising model and to a two-component classical lattice gas corresponding to an extended Hubbard model at half filling in the zero-bandwidth limit. We use a resummation of the LCE to generate $1/d$ expansions for the equation of state and for the critical temperature. The method, which is rather general and applicable to a wide range of models, proves convenient for calculating asymptotic power series expansions in $1/d.$ The vertex renormalized expansion is shown to break down at low temperature in higher order approximations, barring attempts to construct simple approximations that are both self-consistent and exact to some finite order in $1/d.$

Physical review. B, Condensed matter•Dec 11, 2000•8 citations

High temperature series expansions of S = 1/2 Heisenberg spin models: Algorithm to include the magnetic field with optimized complexity

Laurent Pierre, B. Bernu, Laura Messio

This work presents an algorithm for calculating high temperature series expansions (HTSE) of Heisenberg spin models with spin S=1/2 S = 1 / 2 in the thermodynamic limit. This algorithm accounts for the presence of a magnetic field. The paper begins with a comprehensive introduction to HTSE and then focuses on identifying the bottlenecks that limit the computation of higher order coefficients. HTSE calculations involve two key steps: graph enumeration on the lattice and trace calculations for each graph. The introduction of a non-zero magnetic field adds complexity to the expansion because previously irrelevant graphs must now be considered: bridged graphs. We present an efficient method to deduce the contribution of these graphs from the contribution of sub-graphs, that drastically reduces the time of calculation for the last order coefficient (in practice increasing by one the order of the series at almost no cost). Previous articles of the authors have utilized HTSE calculations based on this algorithm, but without providing detailed explanations. The complete algorithm is publicly available, as well as the series on many lattice and for various interactions.

arXiv:2404.02271•SciPost Physics•Oct 7, 2024•9 citations

Thermodynamics of two-dimensional spin models with bimodal random-bond disorder

Baoming Tang, Deepak Iyer, Marcos Rigol

We use numerical linked cluster expansions to study the thermodynamic properties of the two-dimensional classical Ising, quantum $XY$, and quantum Heisenberg models with bimodal random-bond disorder on the square and honeycomb lattice. In all cases, the nearest-neighbor coupling between the spins takes values $\ifmmode\pm\else\textpm\fi{}J$ with equal probability. We obtain the disorder-averaged (over all disorder configurations) energy, entropy, specific heat, and uniform magnetic susceptibility in each case. These results are compared with the corresponding ones in the clean models. Analytic expressions are obtained for low orders in the expansion of these thermodynamic quantities in inverse temperature.

arXiv:1501.00990•Physical Review B•May 14, 2015•18 citations

High-temperature series expansion for spin glasses. I. Derivation of the series

Rajiv R. P. Singh, Sudip Chakravarty

It is shown that the free energy and the Edwards-Anderson susceptibility for Ising spin glasses have a star-graph expansion. With the use of this, a method for generating the high-temperature series for these quantities is described. The scope and the applicability of this method to other properties involving quenched disorder is discussed.

Physical review. B, Condensed matter•Jul 1, 1987•35 citations

Test of Universality in the Ising Spin Glass Using High Temperature Graph Expansion.

D. Daboul, Iksoo Chang, Amnon Aharony A

We calculate high-temperature graph expansions for the Ising spin glass model with 4 symmetric random distribution functions for its nearest neighbor interaction constants Jij. Series for the Edwards-Anderson susceptibility χ are obtained to order 13 in the expansion variable (J/(kBT)) EA 2 for the general d-dimensional hyper-cubic lattice, where the parameter J determines the width of the distributions. We explain in detail how the expansions are calculated. The analysis, using the Dlog-Padé approximation and the techniques known as M1 and M2, leads to estimates for the critical threshold (J/(kBTc)) 2 and for the critical exponent γ in dimensions 4, 5, 7 and 8 for all the distribution functions. In each dimension the values for γ agree, within their uncertainty margins, with a common value for the different distributions, thus confirming universality. PACS: 05.70.Jk Critical point phenomena – 75.10.Nr Spin-glass and other random models 1

•Jan 1, 2004•15 citations

Linked-cluster expansions for quantum magnets on the hypercubic lattice

K. Coester, Darshan G. Joshi, Matthias Vojta C +1 more

For arbitrary space dimension $d$, we investigate the quantum phase transitions of two paradigmatic spin models defined on a hypercubic lattice, the coupled-dimer Heisenberg model and the transverse-field Ising model. To this end, high-order linked-cluster expansions for the ground-state energy and the one-particle gap are performed up to order 9 about the decoupled-dimer and high-field limits, respectively. Extrapolations of the high-order series yield the location of the quantum phase transition and the correlation-length exponent $\ensuremath{\nu}$ as a function of space dimension $d$. The results are complemented by $1/d$ expansions to next-to-leading order of observables across the phase diagrams. Remarkably, our analysis of the extrapolated linked-cluster expansion allows to extract the coefficients of the full $1/d$ expansion for the phase-boundary location in both models exactly in leading order and quantitatively for subleading corrections.

arXiv:1605.05848•Physical review. B./Physical review. B•Sep 6, 2016•17 citations
We discuss generation of series expansions for Ising spin glasses with a symmetric \ifmmode±\else\textpm\fiJ\ifmmode\pm\else\textpm\fi{}J\ifmmode±\else\textpm\fiJ (i.e., bimodal) distribution on ddd-dimensional hypercubic lattices using linked-cluster methods. Simplifications for the bimodal distribution allow us to go to higher order than for a general distribution. We discuss two types of problems, one classical and one quantum. The classical problem is that of the Ising spin glass in a longitudinal magnetic field hhh, for which we obtain high temperature series expansions in variables tanh(J/T)tanh(J/T)tanh(J/T) and tanh(h/T)tanh(h/T)tanh(h/T). The quantum problem is a T=0T=0T=0 study of the Ising spin glass in a transverse magnetic field hT{h}_{T}hT​ for which we obtain a perturbation theory in powers of J/hTJ/{h}_{T}J/hT​. These methods require (i) enumeration and counting of all connected clusters that can be embedded in the lattice up to some order nnn, and (ii) an evaluation of the contribution of each cluster for the quantity being calculated, known as the weight. We discuss a general method that takes the much smaller list (and count) of all no free-end (NFE) clusters on a lattice up to some order nnn and automatically generates all other clusters and their counts up to the same order. The weights for finite clusters in both cases have a simple graphical interpretation that allows us to proceed efficiently for a general configuration of the \ifmmode±\else\textpm\fiJ\ifmmode\pm\else\textpm\fi{}J\ifmmode±\else\textpm\fiJ bonds and at the end perform suitable disorder averaging. The order of our computations is limited by the weight calculations for the high-temperature expansions of the classical model, while they are limited by graph counting for the T=0T=0T=0 quantum system. Details of the calculational methods are presented.
We consider 1/d1/d1/d expansions for classical spin systems based on the vertex renormalized linked cluster expansion (LCE). The free multiplicities of the LCE graphs on a hypercubic lattice in an arbitrary dimension d are calculated using generating functions. The technique is applied to the Ising model and to a two-component classical lattice gas corresponding to an extended Hubbard model at half filling in the zero-bandwidth limit. We use a resummation of the LCE to generate 1/d1/d1/d expansions for the equation of state and for the critical temperature. The method, which is rather general and applicable to a wide range of models, proves convenient for calculating asymptotic power series expansions in 1/d.1/d.1/d. The vertex renormalized expansion is shown to break down at low temperature in higher order approximations, barring attempts to construct simple approximations that are both self-consistent and exact to some finite order in 1/d.1/d.1/d.
This work presents an algorithm for calculating high temperature series expansions (HTSE) of Heisenberg spin models with spin S=1/2 S = 1 / 2 in the thermodynamic limit. This algorithm accounts for the presence of a magnetic field. The paper begins with a comprehensive introduction to HTSE and then focuses on identifying the bottlenecks that limit the computation of higher order coefficients. HTSE calculations involve two key steps: graph enumeration on the lattice and trace calculations for each graph. The introduction of a non-zero magnetic field adds complexity to the expansion because previously irrelevant graphs must now be considered: bridged graphs. We present an efficient method to deduce the contribution of these graphs from the contribution of sub-graphs, that drastically reduces the time of calculation for the last order coefficient (in practice increasing by one the order of the series at almost no cost). Previous articles of the authors have utilized HTSE calculations based on this algorithm, but without providing detailed explanations. The complete algorithm is publicly available, as well as the series on many lattice and for various interactions.
We use numerical linked cluster expansions to study the thermodynamic properties of the two-dimensional classical Ising, quantum XYXYXY, and quantum Heisenberg models with bimodal random-bond disorder on the square and honeycomb lattice. In all cases, the nearest-neighbor coupling between the spins takes values $\ifmmode\pm\else\textpm\fi{}J$ with equal probability. We obtain the disorder-averaged (over all disorder configurations) energy, entropy, specific heat, and uniform magnetic susceptibility in each case. These results are compared with the corresponding ones in the clean models. Analytic expressions are obtained for low orders in the expansion of these thermodynamic quantities in inverse temperature.
We calculate high-temperature graph expansions for the Ising spin glass model with 4 symmetric random distribution functions for its nearest neighbor interaction constants Jij. Series for the Edwards-Anderson susceptibility χ are obtained to order 13 in the expansion variable (J/(kBT)) EA 2 for the general d-dimensional hyper-cubic lattice, where the parameter J determines the width of the distributions. We explain in detail how the expansions are calculated. The analysis, using the Dlog-Padé approximation and the techniques known as M1 and M2, leads to estimates for the critical threshold (J/(kBTc)) 2 and for the critical exponent γ in dimensions 4, 5, 7 and 8 for all the distribution functions. In each dimension the values for γ agree, within their uncertainty margins, with a common value for the different distributions, thus confirming universality. PACS: 05.70.Jk Critical point phenomena – 75.10.Nr Spin-glass and other random models 1
For arbitrary space dimension ddd, we investigate the quantum phase transitions of two paradigmatic spin models defined on a hypercubic lattice, the coupled-dimer Heisenberg model and the transverse-field Ising model. To this end, high-order linked-cluster expansions for the ground-state energy and the one-particle gap are performed up to order 9 about the decoupled-dimer and high-field limits, respectively. Extrapolations of the high-order series yield the location of the quantum phase transition and the correlation-length exponent $\ensuremath{\nu}$ as a function of space dimension ddd. The results are complemented by 1/d1/d1/d expansions to next-to-leading order of observables across the phase diagrams. Remarkably, our analysis of the extrapolated linked-cluster expansion allows to extract the coefficients of the full 1/d1/d1/d expansion for the phase-boundary location in both models exactly in leading order and quantitatively for subleading corrections.