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

Numerical linked-cluster expansions for disordered lattice models

Michael Mulanix, Demetrius Almada, Ehsan Khatami

Published May 10, 2019• 8 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.99.205113

Related papers

Numerical Studies of Correlated Lattice Systems in One and Two Dimensions

Baoming Tang

We use numerical linked cluster expansions (NLCEs) to study correlated lattice systems in one and two dimensions.For fermions in the honeycomb lattice, we study its finite-temperature properties and short-range spin correlations using NLCEs and determinantal quantum Monte Carlo (DQMC) simulations.For the homogeneous system, we analyze a number of thermodynamic quantities, including the entropy, the specific heat, uniform and staggered spin susceptibilities, short-range spin correlations, and the double occupancy at and away from half filling.Employing a local density approximation (LDA), we examine the viability of adiabatic cooling by increasing the interaction strength for homogeneous as well as for trapped systems.We also use NLCEs to study thermodynamic properties of the two-dimensional spin-1/2 Ising, XY, and Heisenberg models with bimodal random-bond disorder on the square and honeycomb lattices.In all cases, the nearest-neighbor coupling between the spins takes values ±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.For many-body localization in disordered isolated systems, we study quantum quenches in the thermodynamic limit.By a quantum quench it is meant that the initial state is stationary with respect iii to an initial Hamiltonian, which is suddenly changed to a new (time-independent) Hamiltonian.The latter then drives the (unitary) dynamics of the system.We are interested in the time average of observables after relaxation following the quench, which can be computed under the so-called diagonal ensemble.

Georgetown University-Graduate School of Arts & Sciences•Jan 1, 2015•0 citations

Thermodynamics of the disordered Hubbard model studied via numerical linked-cluster expansions

Jacob Park, Ehsan Khatami

The interplay of disorder and strong correlations in quantum many-body systems remains an open question. That is despite much progress made in recent years with ultracold atoms in optical lattices to better understand phenomena such as many-body localization or the effect of disorder on Mott metal-insulator transitions. Here, we utilize the numerical linked-cluster expansion technique, extended to treat disordered quantum lattice models, and study exact thermodynamic properties of the disordered Fermi-Hubbard model on the square and cubic geometries. We consider box distributions for the disorder in the onsite energy, the interaction strength, as well as the hopping amplitude and explore how energy, double occupancy, entropy, heat capacity, and magnetic correlations of the system in the thermodynamic limit evolve as the strength of disorder changes. We compare our findings with those obtained from determinant quantum Monte Carlo simulations and discuss the relevance of our results to experiments with cold fermionic atoms in optical lattices.

arXiv:2101.12721•Physical review. B./Physical review. B•Oct 4, 2021•6 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 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

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

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
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 use numerical linked cluster expansions (NLCEs) to study correlated lattice systems in one and two dimensions.For fermions in the honeycomb lattice, we study its finite-temperature properties and short-range spin correlations using NLCEs and determinantal quantum Monte Carlo (DQMC) simulations.For the homogeneous system, we analyze a number of thermodynamic quantities, including the entropy, the specific heat, uniform and staggered spin susceptibilities, short-range spin correlations, and the double occupancy at and away from half filling.Employing a local density approximation (LDA), we examine the viability of adiabatic cooling by increasing the interaction strength for homogeneous as well as for trapped systems.We also use NLCEs to study thermodynamic properties of the two-dimensional spin-1/2 Ising, XY, and Heisenberg models with bimodal random-bond disorder on the square and honeycomb lattices.In all cases, the nearest-neighbor coupling between the spins takes values ±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.For many-body localization in disordered isolated systems, we study quantum quenches in the thermodynamic limit.By a quantum quench it is meant that the initial state is stationary with respect iii to an initial Hamiltonian, which is suddenly changed to a new (time-independent) Hamiltonian.The latter then drives the (unitary) dynamics of the system.We are interested in the time average of observables after relaxation following the quench, which can be computed under the so-called diagonal ensemble.
The interplay of disorder and strong correlations in quantum many-body systems remains an open question. That is despite much progress made in recent years with ultracold atoms in optical lattices to better understand phenomena such as many-body localization or the effect of disorder on Mott metal-insulator transitions. Here, we utilize the numerical linked-cluster expansion technique, extended to treat disordered quantum lattice models, and study exact thermodynamic properties of the disordered Fermi-Hubbard model on the square and cubic geometries. We consider box distributions for the disorder in the onsite energy, the interaction strength, as well as the hopping amplitude and explore how energy, double occupancy, entropy, heat capacity, and magnetic correlations of the system in the thermodynamic limit evolve as the strength of disorder changes. We compare our findings with those obtained from determinant quantum Monte Carlo simulations and discuss the relevance of our results to experiments with cold fermionic atoms in optical lattices.
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 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.
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.
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.