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:cond-mat/0611102 · Physical Review Letters

Numerical Linked-Cluster Approach to Quantum Lattice Models

Marcos Rigol, Tyler Bryant, Rajiv R. P. Singh

Published Nov 3, 2006• 171 citations
arXivHTMLPDF

Abstract

We present a novel algorithm that allows one to obtain temperature dependent properties of quantum lattice models in the thermodynamic limit from exact diagonalization of small clusters. Our numerical linked-cluster approach provides a systematic framework to assess finite-size effects and is valid for any quantum lattice model. Unlike high temperature expansions, which have a finite radius of convergence in inverse temperature, these calculations are accurate at all temperatures provided the range of correlations is finite. We illustrate the power of our approach studying spin models on kagomé, triangular, and square lattices.

DOI
10.1103/physrevlett.97.187202

Related papers

Numerical linked-cluster algorithms. I. Spin systems on square, triangular, and kagomé lattices

Marcos Rigol, Tyler Bryant, Rajiv R. P. Singh

We discuss recently introduced numerical linked-cluster (NLC) algorithms that allow one to obtain temperature-dependent properties of quantum lattice models, in the thermodynamic limit, from exact diagonalization of finite clusters. We present studies of thermodynamic observables for spin models on square, triangular, and kagomé lattices. Results for several choices of clusters and extrapolations methods, that accelerate the convergence of NLCs, are presented. We also include a comparison of NLC results with those obtained from exact analytical expressions (where available), high-temperature expansions (HTE), exact diagonalization (ED) of finite periodic systems, and quantum Monte Carlo simulations. For many models and properties NLC results are substantially more accurate than HTE and ED.

arXiv:0706.3254•Physical Review E•Jun 21, 2007•113 citations

Lanczos-boosted numerical linked-cluster expansion for quantum lattice models

Krishnakumar Bhattaram, Ehsan Khatami

Numerical linked-cluster expansions allow one to calculate finite-temperature properties of quantum lattice models directly in the thermodynamic limit through exact solutions of small clusters. However, full diagonalization is often the limiting factor for these calculations. Here we show that a partial diagonalization of the largest clusters in the expansion using the Lanczos algorithm can be as useful as full diagonalization for the method while mitigating some of the time and memory issues. As test cases, we consider the frustrated Heisenberg model on the checkerboard lattice and the Fermi-Hubbard model on the square lattice. We find that our approach can surpass state of the art in performance in a parallel environment.

arXiv:1810.06202•Physical review. E•Jul 11, 2019•10 citations

Numerical linked-cluster algorithms. II. t−J models on the square lattice

Marcos Rigol, Tyler Bryant, Rajiv R. P. Singh

We discuss the application of a recently introduced numerical linked-cluster (NLC) algorithm to strongly correlated itinerant models. In particular, we present a study of thermodynamic observables: chemical potential, entropy, specific heat, and uniform susceptibility for the t-J model on the square lattice, with Jt=0.5 and 0.3. Our NLC results are compared with those obtained from high-temperature expansions (HTE) and the finite-temperature Lanczos method (FTLM). We show that there is a sizeable window in temperature where NLC results converge without extrapolations whereas HTE diverges. Upon extrapolations, the overall agreement between NLC, HTE, and FTLM is excellent in some cases down to 0.25t . At intermediate temperatures NLC results are better controlled than other methods, making it easier to judge the convergence and numerical accuracy of the method.

arXiv:0706.3255•Physical Review E•Jun 21, 2007•61 citations

Numerical linked-cluster expansions for disordered lattice models

Michael Mulanix, Demetrius Almada, Ehsan Khatami

Imperfections in correlated materials can alter their ground state as well as finite-temperature properties in significant ways. Here, we develop a method based on numerical linked-cluster expansions for calculating exact finite-temperature properties of disordered lattice models directly in the thermodynamic limit. We show that a continuous distribution for disordered parameters can be achieved using a set of carefully chosen discrete modes in the distribution, which allows for the averaging of properties over all disorder realizations. We benchmark our results for thermodynamic properties of the square-lattice Ising and quantum Heisenberg models with bond disorder against Monte Carlo simulations and study them as the strength of disorder changes. We also apply the method to the disordered Heisenberg model on the frustrated checkerboard lattice, which is closely connected to ${\mathrm{Sr}}_{2}\mathrm{Cu}({\mathrm{Te}}_{0.5}{\mathrm{W}}_{0.5}){\mathrm{O}}_{6}$. Our method can be used to study finite-temperature properties of other disordered quantum lattice models, including those for interacting lattice fermions.

arXiv:1810.12901•Physical review. B./Physical review. B•May 10, 2019•8 citations

L -based numerical linked cluster expansion for square lattice models

Mahmoud Abdelshafy, Marcos Rigol

We introduce a numerical linked cluster expansion for square-lattice models whose building block is an L-shape cluster. For the spin-1/2 models studied in this work, we find that this expansion exhibits a similar or better convergence of the bare sums than that of the (larger) square-shaped clusters and can be used with resummation techniques (like the site- and bond-based expansions) to obtain results at even lower temperatures. We compare the performance of weak- and strong-embedding versions of this expansion in various spin-1/2 models and show that the strong-embedding version is preferable because of its convergence properties and lower computational cost. Finally, we show that the expansion based on the L-shape cluster can be naturally used to study properties of lattice models that smoothly connect the square and triangular lattice geometries.

arXiv:2303.02458•Physical review. E•Sep 29, 2023•3 citations

Numerical linked-cluster expansions for two-dimensional spin models with continuous disorder distributions

Mahmoud Abdelshafy, Marcos Rigol

We show that numerical linked cluster expansions (NLCEs) based on sufficiently large building blocks allow one to obtain accurate low-temperature results for the thermodynamic properties of spin lattice models with continuous disorder distributions. Specifically, we show that such results can be obtained computing the disorder averages in the NLCE clusters before calculating their weights. We provide a proof of concept using three different NLCEs based on L, square, and rectangle building blocks. We consider both classical (Ising) and quantum (Heisenberg) spin-1/2 models and show that convergence can be achieved down to temperatures that are up to two orders of magnitude lower than the relevant energy scale in the model. Additionally, we provide evidence that in one dimension one can obtain accurate results for observables such as the energy down to their ground-state values.

arXiv:2402.00931•Physical review. E•May 24, 2024•3 citations
We present a novel algorithm that allows one to obtain temperature dependent properties of quantum lattice models in the thermodynamic limit from exact diagonalization of small clusters. Our numerical linked-cluster approach provides a systematic framework to assess finite-size effects and is valid for any quantum lattice model. Unlike high temperature expansions, which have a finite radius of convergence in inverse temperature, these calculations are accurate at all temperatures provided the range of correlations is finite. We illustrate the power of our approach studying spin models on kagomé, triangular, and square lattices.
We discuss recently introduced numerical linked-cluster (NLC) algorithms that allow one to obtain temperature-dependent properties of quantum lattice models, in the thermodynamic limit, from exact diagonalization of finite clusters. We present studies of thermodynamic observables for spin models on square, triangular, and kagomé lattices. Results for several choices of clusters and extrapolations methods, that accelerate the convergence of NLCs, are presented. We also include a comparison of NLC results with those obtained from exact analytical expressions (where available), high-temperature expansions (HTE), exact diagonalization (ED) of finite periodic systems, and quantum Monte Carlo simulations. For many models and properties NLC results are substantially more accurate than HTE and ED.
Numerical linked-cluster expansions allow one to calculate finite-temperature properties of quantum lattice models directly in the thermodynamic limit through exact solutions of small clusters. However, full diagonalization is often the limiting factor for these calculations. Here we show that a partial diagonalization of the largest clusters in the expansion using the Lanczos algorithm can be as useful as full diagonalization for the method while mitigating some of the time and memory issues. As test cases, we consider the frustrated Heisenberg model on the checkerboard lattice and the Fermi-Hubbard model on the square lattice. We find that our approach can surpass state of the art in performance in a parallel environment.
We discuss the application of a recently introduced numerical linked-cluster (NLC) algorithm to strongly correlated itinerant models. In particular, we present a study of thermodynamic observables: chemical potential, entropy, specific heat, and uniform susceptibility for the t-J model on the square lattice, with Jt=0.5 and 0.3. Our NLC results are compared with those obtained from high-temperature expansions (HTE) and the finite-temperature Lanczos method (FTLM). We show that there is a sizeable window in temperature where NLC results converge without extrapolations whereas HTE diverges. Upon extrapolations, the overall agreement between NLC, HTE, and FTLM is excellent in some cases down to 0.25t . At intermediate temperatures NLC results are better controlled than other methods, making it easier to judge the convergence and numerical accuracy of the method.
Imperfections in correlated materials can alter their ground state as well as finite-temperature properties in significant ways. Here, we develop a method based on numerical linked-cluster expansions for calculating exact finite-temperature properties of disordered lattice models directly in the thermodynamic limit. We show that a continuous distribution for disordered parameters can be achieved using a set of carefully chosen discrete modes in the distribution, which allows for the averaging of properties over all disorder realizations. We benchmark our results for thermodynamic properties of the square-lattice Ising and quantum Heisenberg models with bond disorder against Monte Carlo simulations and study them as the strength of disorder changes. We also apply the method to the disordered Heisenberg model on the frustrated checkerboard lattice, which is closely connected to Sr2Cu(Te0.5W0.5)O6{\mathrm{Sr}}_{2}\mathrm{Cu}({\mathrm{Te}}_{0.5}{\mathrm{W}}_{0.5}){\mathrm{O}}_{6}Sr2​Cu(Te0.5​W0.5​)O6​. Our method can be used to study finite-temperature properties of other disordered quantum lattice models, including those for interacting lattice fermions.
We introduce a numerical linked cluster expansion for square-lattice models whose building block is an L-shape cluster. For the spin-1/2 models studied in this work, we find that this expansion exhibits a similar or better convergence of the bare sums than that of the (larger) square-shaped clusters and can be used with resummation techniques (like the site- and bond-based expansions) to obtain results at even lower temperatures. We compare the performance of weak- and strong-embedding versions of this expansion in various spin-1/2 models and show that the strong-embedding version is preferable because of its convergence properties and lower computational cost. Finally, we show that the expansion based on the L-shape cluster can be naturally used to study properties of lattice models that smoothly connect the square and triangular lattice geometries.
We show that numerical linked cluster expansions (NLCEs) based on sufficiently large building blocks allow one to obtain accurate low-temperature results for the thermodynamic properties of spin lattice models with continuous disorder distributions. Specifically, we show that such results can be obtained computing the disorder averages in the NLCE clusters before calculating their weights. We provide a proof of concept using three different NLCEs based on L, square, and rectangle building blocks. We consider both classical (Ising) and quantum (Heisenberg) spin-1/2 models and show that convergence can be achieved down to temperatures that are up to two orders of magnitude lower than the relevant energy scale in the model. Additionally, we provide evidence that in one dimension one can obtain accurate results for observables such as the energy down to their ground-state values.