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

High-temperature series expansion of the dynamic Matsubara spin correlator

Ruben Burkard, Benedikt Schneider, Björn Sbierski

Published Oct 1, 2025• 3 citations
arXivHTMLPDF

Abstract

The high-temperature series expansion for quantum spin models is a well-established tool to compute thermodynamic quantities and equal-time spin correlations, in particular for frustrated interactions. We extend the scope of this expansion to the dynamic Matsubara spin-spin correlator and develop an algorithm that yields exact expansion coefficients in the form of rational numbers. We focus on Heisenberg models with a single coupling constant J and spin lengths S=1/2,1. The expansion coefficients up to 12th order in J/T are precomputed on all possible $\sim 10^6$ graphs embeddable in arbitrary lattices and are provided in a repository. This enables calculation of static momentum-resolved susceptibilities for arbitrary site-pairs or wavevectors. We test our results for the antiferromagnetic S=1/2 chain and triangular lattice model. An important application that we discuss in a companion letter is the calculation of real-frequency dynamic structure factors. This is achieved by identifying the high-frequency expansion coefficients of the Matsubara correlator with frequency moments of the spectral function.

DOI
10.1103/1l92-z6qd

Related papers

Dynamic Correlations of Frustrated Quantum Spins from High-Temperature Expansion

Ruben Burkard, Benedikt Schneider, Björn Sbierski

For quantum spin systems in equilibrium, the dynamic structure factor (DSF) is among the most feature-packed experimental observables. However, from a theory perspective it is often hard to simulate in an unbiased and accurate way, especially for frustrated and high-dimensional models at intermediate temperature. To address this challenge, we compute the DSF from a dynamic extension of the high-temperature expansion to frequency moments. We focus on nearest-neighbor Heisenberg models with spin lengths S=1/2 and 1. We provide comprehensive benchmarks and consider a variety of frustrated two- and three-dimensional antiferromagnets as applications. In particular we shed new light on the anomalous intermediate temperature regime of the S=1/2 triangular lattice model and reproduce the DSF measured recently for the S=1 pyrochlore material NaCaNi_{2}F_{7}. An open-source numerical implementation for arbitrary lattice geometries is also provided.

arXiv:2505.14571•Physical Review Letters•Jan 7, 2026•3 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

Diagrammatic approach to frustrated spin systems

Benedikt Schneider

This thesis is devoted to diagrammatic frameworks for computing the dynamic structure factor (DSF) and the static susceptibility of spin systems at finite temperature in the thermodynamic limit. We consider two complementary diagrammatic frameworks. First, we study functional renormalization group (FRG) approaches based on pseudo-fermion (pf-FRG) and pseudo-Majorana (pm-FRG) representations. Second, we develop high-order strong-coupling and high-temperature expansions (HTE) using free and restricted-graph techniques. To enable finite-temperature calculations in pf-FRG, we eliminate unphysical states by incorporating the Popov–Fedotov trick. Within pm-FRG, we introduce a temperature-flow scheme that significantly improves the efficiency and accuracy of finite-size scaling analyses of thermal phase transitions. We further extend pm-FRG to include external magnetic fields, allowing treatment of XXZ models with a field along the z direction. While the pseudo-Majorana representation maps local spin-1/2 operators to free Majorana fermions, higher spin values S>1/2 can instead be represented as strongly interacting real bosons, an idea underlying the recently proposed spin-FRG. Rather than implementing the spin-FRG itself, we analyze its perturbation series using a strong-coupling expansion. To evaluate the high-dimensional imaginary-time integrals arising in this expansion analytically with efficient bookkeeping, we introduce the kernel trick, which rewrites any Matsubara diagram in terms of the partial spectral functions of its constituent Green’s functions. The resulting fourth-order strong-coupling expansion for general spin S provides systematic beyond-mean-field corrections to critical temperatures and fields in both 1/d and large-interaction-range expansions. More importantly, it provides an analytic explanation of the long-standing quantum-to-classical correspondence observed in diagrammatic Monte Carlo simulations: the normalized static susceptibilities of quantum and classical Heisenberg models coincide after a simple temperature rescaling. Our analysis shows that this arises from the fact that the static susceptibility is well approximated by a simple functional form akin to a mean-field susceptibility, which we call the renormalized-mean-field form. To reach higher orders, we generalize the HTE to Matsubara-frequency-resolved Green’s functions (Dyn-HTE). This includes the static susceptibility (at zero Matsubara frequency) as a special case and validates the renormalized-mean-field form up to 12th order in perturbation theory. By analytically continuing the frequency dependence obtained from Dyn-HTE using a continued-fraction expansion, we compute DSFs for one-parameter Heisenberg models on arbitrary lattices, which can be directly compared with experimental inelastic neutron-scattering spectra. All Dyn-HTE calculations are implemented in a publicly available, well-documented open-source code. Finally, we clarify the relation between pseudo-Majorana and spin-based diagrammatic approaches by deriving the exact pseudo-Majorana action directly from the spin Hamiltonian, without invoking any parton representation. This derivation paves the way for extending Majorana-based methods to higher-spin systems.

•Jan 1, 2025•0 citations

High Temperature Expansion for the SU(n) Heisenberg Model in One Dimension

Noboru Fukushima, Yoshio Kuramoto

Thermodynamic properties of the SU($n$) Heisenberg model in one dimension is studied by means of high-temperature expansion for arbitrary $n$. The specific heat up to $O[(βJ)^{23}]$ and the correlation function up to $O[(βJ)^{18}]$ are derived with $βJ$ being the antiferromagnetic exchange in units of temperature. It is found for $n>2$ that the specific heat shows a shoulder in the high-temperature side of a peak. The origin of this structure is clarified by deriving the temperature dependence of the correlation function. With decreasing temperature, the short-range correlation with two-site periodicity develops first, and then another correlation with $n$-site periodicity at lower temperature. This behavior is in contrast to that of the inverse square interaction model, where the specific heat shows a single peak according to the exact solution. Our algorithm has an advantage that neither computational time nor memory depends on the multiplicity $n$ per site; the series coefficients are obtained as explicit functions of $n$.

arXiv:cond-mat/0110550•Journal of the Physical Society of Japan•May 15, 2002•10 citations

Finite-temperature phase transitions in S=12 three-dimensional Heisenberg magnets from high-temperature series expansions

M. G. Gonzalez, B. Bernu, Laurent Pierre +1 more

Many frustrated spin models on three-dimensional (3D) lattices are currently being investigated, both experimentally and theoretically, and develop new types of long-range orders in their respective phase diagrams. They present finite-temperature phase transitions, most likely in the Heisenberg 3D universality class. However, the combination between the 3D character and frustration makes them hard to study. We present here several methods derived from high-temperature series expansions (HTSEs), which give exact coefficients directly in the thermodynamic limit up to a certain order; for several 3D lattices, supplementary orders than in previous literature are reported for the HTSEs. We introduce an interpolation method able to describe thermodynamic quantities at $T>{T}_{c}$, which we use here to reconstruct the magnetic susceptibility and the specific heat and to extract universal and nonuniversal quantities (for example, critical exponents, temperature, energy, entropy, and other parameters related to the phase transition). While the susceptibility associated with the order parameter is not usually known for more exotic long-range orders, the specific heat is indicative of a phase transition for any kind of symmetry breaking. We present examples of applications on ferromagnetic and antiferromagnetic models on various 3D lattices and benchmark our results whenever possible.

arXiv:2303.03135•Physical review. B./Physical review. B•Jun 29, 2023•7 citations

Dynamics of strong-coupling models for cuprate superconductors: Exact results on finite lattices

H. B. Schuettler, C.-X. Chen, A. J. Fedro

We discuss recent applications of exact numerical continued fraction expansion (CFE) techniques to calculate dynamical correlation functions of various strong-coupling models related to the high-temperature cuprate superconductors. For the two-dimensional square-lattice spin-1/2 Heisenberg antiferromagnet, we present exact results for the zero-temperature dynamical structure factor on finite-sized lattices and compare them to approximate results from a Schwinger boson mean-field theory, recently proposed by Arovas and Auerbach. We find that the mean-field theory represents a very good approximation to the exact spin excitation spectra and to the static spin correlations. We then investigation the dynamical spin-spin structure factor and the single-particle spectral function for finite model clusters with dopant induced hole-type charge carriers, in the strong-coupling limit of both the single-band Hubbard model (t-J-model) and the three-band Hubbard model (Kondo-Heisenberg model). Our results are consistent with the physical picture recently proposed by Zhang and Rice which implies an approximate mapping of the low-energy states in the three-band model onto an effective single-band theory. 33 refs., 4 figs.

•Jan 1, 1989•0 citations
The high-temperature series expansion for quantum spin models is a well-established tool to compute thermodynamic quantities and equal-time spin correlations, in particular for frustrated interactions. We extend the scope of this expansion to the dynamic Matsubara spin-spin correlator and develop an algorithm that yields exact expansion coefficients in the form of rational numbers. We focus on Heisenberg models with a single coupling constant J and spin lengths S=1/2,1. The expansion coefficients up to 12th order in J/T are precomputed on all possible ∼106\sim 10^6∼106 graphs embeddable in arbitrary lattices and are provided in a repository. This enables calculation of static momentum-resolved susceptibilities for arbitrary site-pairs or wavevectors. We test our results for the antiferromagnetic S=1/2 chain and triangular lattice model. An important application that we discuss in a companion letter is the calculation of real-frequency dynamic structure factors. This is achieved by identifying the high-frequency expansion coefficients of the Matsubara correlator with frequency moments of the spectral function.
For quantum spin systems in equilibrium, the dynamic structure factor (DSF) is among the most feature-packed experimental observables. However, from a theory perspective it is often hard to simulate in an unbiased and accurate way, especially for frustrated and high-dimensional models at intermediate temperature. To address this challenge, we compute the DSF from a dynamic extension of the high-temperature expansion to frequency moments. We focus on nearest-neighbor Heisenberg models with spin lengths S=1/2 and 1. We provide comprehensive benchmarks and consider a variety of frustrated two- and three-dimensional antiferromagnets as applications. In particular we shed new light on the anomalous intermediate temperature regime of the S=1/2 triangular lattice model and reproduce the DSF measured recently for the S=1 pyrochlore material NaCaNi_{2}F_{7}. An open-source numerical implementation for arbitrary lattice geometries is also provided.
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.
This thesis is devoted to diagrammatic frameworks for computing the dynamic structure factor (DSF) and the static susceptibility of spin systems at finite temperature in the thermodynamic limit. We consider two complementary diagrammatic frameworks. First, we study functional renormalization group (FRG) approaches based on pseudo-fermion (pf-FRG) and pseudo-Majorana (pm-FRG) representations. Second, we develop high-order strong-coupling and high-temperature expansions (HTE) using free and restricted-graph techniques. To enable finite-temperature calculations in pf-FRG, we eliminate unphysical states by incorporating the Popov–Fedotov trick. Within pm-FRG, we introduce a temperature-flow scheme that significantly improves the efficiency and accuracy of finite-size scaling analyses of thermal phase transitions. We further extend pm-FRG to include external magnetic fields, allowing treatment of XXZ models with a field along the z direction. While the pseudo-Majorana representation maps local spin-1/2 operators to free Majorana fermions, higher spin values S>1/2 can instead be represented as strongly interacting real bosons, an idea underlying the recently proposed spin-FRG. Rather than implementing the spin-FRG itself, we analyze its perturbation series using a strong-coupling expansion. To evaluate the high-dimensional imaginary-time integrals arising in this expansion analytically with efficient bookkeeping, we introduce the kernel trick, which rewrites any Matsubara diagram in terms of the partial spectral functions of its constituent Green’s functions. The resulting fourth-order strong-coupling expansion for general spin S provides systematic beyond-mean-field corrections to critical temperatures and fields in both 1/d and large-interaction-range expansions. More importantly, it provides an analytic explanation of the long-standing quantum-to-classical correspondence observed in diagrammatic Monte Carlo simulations: the normalized static susceptibilities of quantum and classical Heisenberg models coincide after a simple temperature rescaling. Our analysis shows that this arises from the fact that the static susceptibility is well approximated by a simple functional form akin to a mean-field susceptibility, which we call the renormalized-mean-field form. To reach higher orders, we generalize the HTE to Matsubara-frequency-resolved Green’s functions (Dyn-HTE). This includes the static susceptibility (at zero Matsubara frequency) as a special case and validates the renormalized-mean-field form up to 12th order in perturbation theory. By analytically continuing the frequency dependence obtained from Dyn-HTE using a continued-fraction expansion, we compute DSFs for one-parameter Heisenberg models on arbitrary lattices, which can be directly compared with experimental inelastic neutron-scattering spectra. All Dyn-HTE calculations are implemented in a publicly available, well-documented open-source code. Finally, we clarify the relation between pseudo-Majorana and spin-based diagrammatic approaches by deriving the exact pseudo-Majorana action directly from the spin Hamiltonian, without invoking any parton representation. This derivation paves the way for extending Majorana-based methods to higher-spin systems.
Thermodynamic properties of the SU(nnn) Heisenberg model in one dimension is studied by means of high-temperature expansion for arbitrary nnn. The specific heat up to O[(βJ)23]O[(βJ)^{23}]O[(βJ)23] and the correlation function up to O[(βJ)18]O[(βJ)^{18}]O[(βJ)18] are derived with βJβJβJ being the antiferromagnetic exchange in units of temperature. It is found for n>2n>2n>2 that the specific heat shows a shoulder in the high-temperature side of a peak. The origin of this structure is clarified by deriving the temperature dependence of the correlation function. With decreasing temperature, the short-range correlation with two-site periodicity develops first, and then another correlation with nnn-site periodicity at lower temperature. This behavior is in contrast to that of the inverse square interaction model, where the specific heat shows a single peak according to the exact solution. Our algorithm has an advantage that neither computational time nor memory depends on the multiplicity nnn per site; the series coefficients are obtained as explicit functions of nnn.
Many frustrated spin models on three-dimensional (3D) lattices are currently being investigated, both experimentally and theoretically, and develop new types of long-range orders in their respective phase diagrams. They present finite-temperature phase transitions, most likely in the Heisenberg 3D universality class. However, the combination between the 3D character and frustration makes them hard to study. We present here several methods derived from high-temperature series expansions (HTSEs), which give exact coefficients directly in the thermodynamic limit up to a certain order; for several 3D lattices, supplementary orders than in previous literature are reported for the HTSEs. We introduce an interpolation method able to describe thermodynamic quantities at T>TcT>{T}_{c}T>Tc​, which we use here to reconstruct the magnetic susceptibility and the specific heat and to extract universal and nonuniversal quantities (for example, critical exponents, temperature, energy, entropy, and other parameters related to the phase transition). While the susceptibility associated with the order parameter is not usually known for more exotic long-range orders, the specific heat is indicative of a phase transition for any kind of symmetry breaking. We present examples of applications on ferromagnetic and antiferromagnetic models on various 3D lattices and benchmark our results whenever possible.
We discuss recent applications of exact numerical continued fraction expansion (CFE) techniques to calculate dynamical correlation functions of various strong-coupling models related to the high-temperature cuprate superconductors. For the two-dimensional square-lattice spin-1/2 Heisenberg antiferromagnet, we present exact results for the zero-temperature dynamical structure factor on finite-sized lattices and compare them to approximate results from a Schwinger boson mean-field theory, recently proposed by Arovas and Auerbach. We find that the mean-field theory represents a very good approximation to the exact spin excitation spectra and to the static spin correlations. We then investigation the dynamical spin-spin structure factor and the single-particle spectral function for finite model clusters with dopant induced hole-type charge carriers, in the strong-coupling limit of both the single-band Hubbard model (t-J-model) and the three-band Hubbard model (Kondo-Heisenberg model). Our results are consistent with the physical picture recently proposed by Zhang and Rice which implies an approximate mapping of the low-energy states in the three-band model onto an effective single-band theory. 33 refs., 4 figs.