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:2303.07754 · Physical review. E

Linked cluster expansion on trees

Deepak Iyer, Yuyi Wan

Published Feb 9, 2023• 1 citation
arXivHTMLPDF

Abstract

The linked cluster expansion has been shown to be highly efficient in calculating equilibrium and nonequilibrium properties of a variety of 1D and 2D classical and quantum lattice models. In this article, we extend the linked cluster method to the Cayley tree and its boundaryless cousin the Bethe lattice. We aim to (a) develop the linked cluster expansion for these lattices, a novel application, and (b) to further understand the surprising convergence efficiency of the linked cluster method, as well as its limitations. We obtain several key results. First, we show that for nearest-neighbor Hamiltonians of a specific form, all finite treelike clusters can be mapped to one dimensional finite chains. We then show that the qualitative distinction between the Cayley tree and Bethe lattice appears due to differing lattice constants that is a result of the Bethe lattice being boundaryless. We use these results to obtain the explicit closed-form formula for the zero-field susceptibility for the entire disordered phase up to the critical point for Bethe lattices of arbitrary degree; remarkably, only 1D chainlike clusters contribute. We also obtain the exact zero field partition function for the Ising model on both trees with only the two smallest clusters, similar to the 1D chain. Finally, these results achieve a direct comparison between an infinite lattice with a nonnegligible boundary and one without any boundary, allowing us to show that the linked cluster expansion eliminates boundary terms at each order of the expansion, answering the question about its surprising convergence efficiency. We conclude with some ramifications of these results, and possible generalizations and applications.

DOI
10.1103/physreve.107.024108

Related papers

Linked-cluster Expansions for Lattice Spin Models

Yuyi Wan

Similar to various series expansions that are used to approximate mathematical func- tions, the linked-cluster expansion is an approximation method that allows us to approach the actual values of a very large physical system’s different physical quan- tities by systematically studying smaller systems embedded in this larger system. The main concept in linked-cluster expansion, weight, represents the additional con- tribution to a certain physical quantity by increasing the system size by one unit. These weights are used to eventually build up the result on a larger system. In our case, we focus on the partition function, a quantity that can be used to calculate several essential thermodynamic aspects of the system such as average energy, spe- cific heat, magnetic susceptibility, etc. We study these weights for the Ising and Potts models on a one dimensional lattice and Bethe lattice, with and without an external magnetic field. Previous studies have shown that in the one dimensional Ising model, adding more than two lattice sites to the system does not create any ad- ditional contribution to the partition function (i.e. their weights are all zero), giving us the result for an infinitely large system with just two lattice sites, a remarkable simplification. In our study, we prove that this property holds not only for the Ising model, but also its generalization, the Potts model, and show that it is a result of a spin flip symmetry inherent in the system. In order to test this, we break this symmetry using an externally applied magnetic field and show that in this case, for magnetic dipole energies comparable with the exchange energy between neighboring spins, this special property vanishes.

•Jan 1, 2020•0 citations

Quantum Ising transition on the Cayley tree: boundary criticality, universality, and an effective Bethe line

Pietro Richelli, Natalia Chepiga

We investigate quantum criticality in the transverse-field Ising model on finite Cayley trees and demonstrate that their extensive boundary fundamentally changes the nature of the quantum phase transition compared with the infinite, boundary-free Bethe lattice. Using tree tensor-network simulations for systems with up to $N=12286$ sites, we identify a boundary-driven continuous transition at $h_x^c\approx1.76$, distinct from the mean-field Bethe-lattice transition at $h_x^c\approx2.23$. We characterize this transition through finite-size scaling of the magnetization and magnetic susceptibility and extract a consistent set of critical exponents. By introducing next-nearest-neighbor interactions, we demonstrate the robustness of the extracted critical exponents, providing numerical argument for the Ising boundary universality class on the Cayley tree. Furthermore, we show that the physics of the boundary-free Bethe lattice can be efficiently recovered on a finite Cayley tree by fine-tuning a local boundary field. The resulting \textit{effective Bethe line}, whose existence is tied to the mean-field nature of the bulk transition, effectively screens the extensive boundary and restores both the critical point and the mean-field critical scaling of the infinite Bethe lattice.

arXiv:2610.16774•Oct 12, 2026•0 citations

Yang-Lee zeros of two-dimensional nearest-neighbor antiferromagnetic Ising models: A numerical linked cluster expansion study

Mahmoud Abdelshafy, Muhammad Sedik

We study Yang-Lee zeros as they approach the thermodynamic limit of the two-dimensional (2D) nearest-neighbor antiferromagnetic Ising model on square and triangular lattices. We employ the numerical linked cluster expansion (NLCE) equipped with exact enumeration (EE) of the partition function to estimate the Laplacian of the free energy, which is proportional to the zeros density. Using a modified NLCE, where the expansion can be carried directly on the Yang-Lee zeros of the involved clusters, we estimate the density of Yang-Lee zeros in the thermodynamic limit. NLCE gives significantly more zeros than EE in the complex field plane providing more insights on how the root curves look in the thermodynamic limit. For the square lattice at $T\ensuremath{\ll}{T}_{c}$, the results suggest that two vertical lines at $\ifmmode\pm\else\textpm\fi{}{h}_{c}(T)$ in the complex field plane (i.e., two concentric circles in the complex fugacity plane) are the thermodynamic root curves. A similar picture is found for the triangular lattice for phase transitions at large values of magnetic field while further study is needed for phase transitions at smaller values of magnetic field. The convergence of the NLCE and (EE) calculations of the partition function to the thermodynamic limit is studied in both lattices and the temperature-field phase diagram is obtained from Yang-Lee zeros using both methods. This NLCE-based approach will facilitate the study of different types of phase transitions using Yang-Lee zeros in future research.

arXiv:2503.14602•Physical review. B./Physical review. B•Aug 8, 2025•2 citations

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

Random transverse field Ising model on the Cayley tree: analysis via boundary strong disorder renormalization

Cécile Monthus, Thomas Garel

Strong disorder renormalization for the random transverse field Ising model leads to a complicated topology of surviving clusters as soon as d > 1. Even if one starts from a Cayley tree, the network of surviving renormalized clusters will contain loops, so that no analytical solution can be obtained. Here we introduce a modified procedure called 'boundary strong disorder renormalization' that preserves the tree structure, so that one can write simple recursions with respect to the number of generations. We first show that this modified procedure allows one to recover exactly most of the critical exponents for the one-dimensional chain. After this important check, we study the RG equations for the quantum Ising model on a Cayley tree with a uniform ferromagnetic coupling J and random transverse fields with support [ h min , h max ]. We find the following picture: (i) for J > h max , only bonds are decimated, so that the whole tree is a quantum ferromagnetic cluster; (ii) for J < h min , only sites are decimated, so that no quantum ferromagnetic cluster is formed, and the ferromagnetic coupling to the boundary coincides with the partition function of a directed polymer model in a random medium; (iii) for h min < J < h max , both sites and bonds can be decimated, and the quantum ferromagnetic clusters can either remain finite (the physics is then similar to (ii), with a quantitative mapping to a modified directed polymer model) or an infinite quantum ferromagnetic cluster appears. We find that the quantum transition can be of two types: (a) either the quantum transition takes place in the region where quantum ferromagnetic clusters remain finite, and the singularity of the ferromagnetic coupling to the boundary involves the typical correlation length exponent ν typ = 1; (b) or the quantum transition takes place at the point where an extensive quantum ferromagnetic cluster appears, with a correlation length exponent ν ≃ 0.75.

arXiv:1205.4512•Journal of Statistical Mechanics Theory and Experiment•Oct 17, 2012•6 citations

Efficient simulation of infinite tree tensor network states on the Bethe lattice

Wei Li, Jan von Delft, Tao Xiang

We show that the simple update approach proposed by Jiang et al. [H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)] is an efficient and accurate method for determining the infinite tree tensor network states on the Bethe lattice. Ground-state properties of the quantum transverse Ising model and the Heisenberg XXZ model on the Bethe lattice are studied. The transverse Ising model is found to undergo a second-order quantum phase transition with a diverging magnetic susceptibility but a finite correlation length which is upper bounded by $1/\mathrm{ln}(q\ensuremath{-}1)$ even at the transition point ($q$ is the coordinate number of the Bethe lattice). An intuitive explanation on this peculiar ``critical'' phenomenon is given. The XXZ model on the Bethe lattice undergoes a first-order quantum phase transition at the isotropic point. Furthermore, the simple update scheme is found to be related with the Bethe approximation. Finally, by applying the simple update to various tree tensor clusters, we can obtain rather nice and scalable approximations for two-dimensional lattices.

arXiv:1209.2387•Physical Review B•Nov 26, 2012•84 citations
The linked cluster expansion has been shown to be highly efficient in calculating equilibrium and nonequilibrium properties of a variety of 1D and 2D classical and quantum lattice models. In this article, we extend the linked cluster method to the Cayley tree and its boundaryless cousin the Bethe lattice. We aim to (a) develop the linked cluster expansion for these lattices, a novel application, and (b) to further understand the surprising convergence efficiency of the linked cluster method, as well as its limitations. We obtain several key results. First, we show that for nearest-neighbor Hamiltonians of a specific form, all finite treelike clusters can be mapped to one dimensional finite chains. We then show that the qualitative distinction between the Cayley tree and Bethe lattice appears due to differing lattice constants that is a result of the Bethe lattice being boundaryless. We use these results to obtain the explicit closed-form formula for the zero-field susceptibility for the entire disordered phase up to the critical point for Bethe lattices of arbitrary degree; remarkably, only 1D chainlike clusters contribute. We also obtain the exact zero field partition function for the Ising model on both trees with only the two smallest clusters, similar to the 1D chain. Finally, these results achieve a direct comparison between an infinite lattice with a nonnegligible boundary and one without any boundary, allowing us to show that the linked cluster expansion eliminates boundary terms at each order of the expansion, answering the question about its surprising convergence efficiency. We conclude with some ramifications of these results, and possible generalizations and applications.
Similar to various series expansions that are used to approximate mathematical func- tions, the linked-cluster expansion is an approximation method that allows us to approach the actual values of a very large physical system’s different physical quan- tities by systematically studying smaller systems embedded in this larger system. The main concept in linked-cluster expansion, weight, represents the additional con- tribution to a certain physical quantity by increasing the system size by one unit. These weights are used to eventually build up the result on a larger system. In our case, we focus on the partition function, a quantity that can be used to calculate several essential thermodynamic aspects of the system such as average energy, spe- cific heat, magnetic susceptibility, etc. We study these weights for the Ising and Potts models on a one dimensional lattice and Bethe lattice, with and without an external magnetic field. Previous studies have shown that in the one dimensional Ising model, adding more than two lattice sites to the system does not create any ad- ditional contribution to the partition function (i.e. their weights are all zero), giving us the result for an infinitely large system with just two lattice sites, a remarkable simplification. In our study, we prove that this property holds not only for the Ising model, but also its generalization, the Potts model, and show that it is a result of a spin flip symmetry inherent in the system. In order to test this, we break this symmetry using an externally applied magnetic field and show that in this case, for magnetic dipole energies comparable with the exchange energy between neighboring spins, this special property vanishes.
We investigate quantum criticality in the transverse-field Ising model on finite Cayley trees and demonstrate that their extensive boundary fundamentally changes the nature of the quantum phase transition compared with the infinite, boundary-free Bethe lattice. Using tree tensor-network simulations for systems with up to N=12286N=12286N=12286 sites, we identify a boundary-driven continuous transition at hxc≈1.76h_x^c\approx1.76hxc​≈1.76, distinct from the mean-field Bethe-lattice transition at hxc≈2.23h_x^c\approx2.23hxc​≈2.23. We characterize this transition through finite-size scaling of the magnetization and magnetic susceptibility and extract a consistent set of critical exponents. By introducing next-nearest-neighbor interactions, we demonstrate the robustness of the extracted critical exponents, providing numerical argument for the Ising boundary universality class on the Cayley tree. Furthermore, we show that the physics of the boundary-free Bethe lattice can be efficiently recovered on a finite Cayley tree by fine-tuning a local boundary field. The resulting \textit{effective Bethe line}, whose existence is tied to the mean-field nature of the bulk transition, effectively screens the extensive boundary and restores both the critical point and the mean-field critical scaling of the infinite Bethe lattice.
We study Yang-Lee zeros as they approach the thermodynamic limit of the two-dimensional (2D) nearest-neighbor antiferromagnetic Ising model on square and triangular lattices. We employ the numerical linked cluster expansion (NLCE) equipped with exact enumeration (EE) of the partition function to estimate the Laplacian of the free energy, which is proportional to the zeros density. Using a modified NLCE, where the expansion can be carried directly on the Yang-Lee zeros of the involved clusters, we estimate the density of Yang-Lee zeros in the thermodynamic limit. NLCE gives significantly more zeros than EE in the complex field plane providing more insights on how the root curves look in the thermodynamic limit. For the square lattice at $T\ensuremath{\ll}{T}_{c}$, the results suggest that two vertical lines at $\ifmmode\pm\else\textpm\fi{}{h}_{c}(T)$ in the complex field plane (i.e., two concentric circles in the complex fugacity plane) are the thermodynamic root curves. A similar picture is found for the triangular lattice for phase transitions at large values of magnetic field while further study is needed for phase transitions at smaller values of magnetic field. The convergence of the NLCE and (EE) calculations of the partition function to the thermodynamic limit is studied in both lattices and the temperature-field phase diagram is obtained from Yang-Lee zeros using both methods. This NLCE-based approach will facilitate the study of different types of phase transitions using Yang-Lee zeros in future research.
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.
Strong disorder renormalization for the random transverse field Ising model leads to a complicated topology of surviving clusters as soon as d > 1. Even if one starts from a Cayley tree, the network of surviving renormalized clusters will contain loops, so that no analytical solution can be obtained. Here we introduce a modified procedure called 'boundary strong disorder renormalization' that preserves the tree structure, so that one can write simple recursions with respect to the number of generations. We first show that this modified procedure allows one to recover exactly most of the critical exponents for the one-dimensional chain. After this important check, we study the RG equations for the quantum Ising model on a Cayley tree with a uniform ferromagnetic coupling J and random transverse fields with support [ h min , h max ]. We find the following picture: (i) for J > h max , only bonds are decimated, so that the whole tree is a quantum ferromagnetic cluster; (ii) for J < h min , only sites are decimated, so that no quantum ferromagnetic cluster is formed, and the ferromagnetic coupling to the boundary coincides with the partition function of a directed polymer model in a random medium; (iii) for h min < J < h max , both sites and bonds can be decimated, and the quantum ferromagnetic clusters can either remain finite (the physics is then similar to (ii), with a quantitative mapping to a modified directed polymer model) or an infinite quantum ferromagnetic cluster appears. We find that the quantum transition can be of two types: (a) either the quantum transition takes place in the region where quantum ferromagnetic clusters remain finite, and the singularity of the ferromagnetic coupling to the boundary involves the typical correlation length exponent ν typ = 1; (b) or the quantum transition takes place at the point where an extensive quantum ferromagnetic cluster appears, with a correlation length exponent ν ≃ 0.75.
We show that the simple update approach proposed by Jiang et al. [H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)] is an efficient and accurate method for determining the infinite tree tensor network states on the Bethe lattice. Ground-state properties of the quantum transverse Ising model and the Heisenberg XXZ model on the Bethe lattice are studied. The transverse Ising model is found to undergo a second-order quantum phase transition with a diverging magnetic susceptibility but a finite correlation length which is upper bounded by $1/\mathrm{ln}(q\ensuremath{-}1)$ even at the transition point (qqq is the coordinate number of the Bethe lattice). An intuitive explanation on this peculiar ``critical'' phenomenon is given. The XXZ model on the Bethe lattice undergoes a first-order quantum phase transition at the isotropic point. Furthermore, the simple update scheme is found to be related with the Bethe approximation. Finally, by applying the simple update to various tree tensor clusters, we can obtain rather nice and scalable approximations for two-dimensional lattices.