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

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

Mahmoud Abdelshafy, Marcos Rigol

Published May 24, 2024• 3 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreve.109.054127

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

Numerical linked cluster expansions for inhomogeneous systems

Johann Gan, Kaden R. A. Hazzard

We develop a numerical linked cluster expansion (NLCE) method that can be applied directly to inhomogeneous systems, for example, Hamiltonians with disorder and dynamics initiated from inhomogeneous initial states. We demonstrate the method by calculating dynamics for single-spin expectations and spin correlations in two-dimensional spin models on a square lattice, starting from a checkerboard state. We show that NLCE can give moderate to dramatic improvement over an exact diagonalization of comparable computational cost and that the advantage in computational resources grows exponentially as the size of the clusters included grows. Although the method applies to any type of NLCE, our explicit benchmarks use the rectangle expansion. Besides showing the capability to treat inhomogeneous systems, these benchmarks demonstrate the rectangle expansion's utility out of equilibrium.

arXiv:2005.03177•Physical Review A•Jul 17, 2020•16 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

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

Numerical linked cluster expansions for quantum quenches in one-dimensional lattices

Krishnanand Mallayya, Marcos Rigol

We discuss the application of numerical linked cluster expansions (NLCEs) to study one dimensional lattice systems in thermal equilibrium and after quantum quenches from thermal equilibrium states. For the former, we calculate observables in the grand canonical ensemble, and for the latter we calculate observables in the diagonal ensemble. When converged, NLCEs provide results in the thermodynamic limit. We use two different NLCEs: a maximally connected expansion introduced in previous works and a site-based expansion. We compare the effectiveness of both NLCEs. The site-based NLCE is found to work best for systems in thermal equilibrium. However, in thermal equilibrium and after quantum quenches, the site-based NLCE can diverge when the maximally connected one converges. We relate this divergence to the exponentially large number of clusters in the site-based NLCE and the behavior of the weights of observables in those clusters. We discuss the effectiveness of resummations to cure the divergence. Our NLCE calculations are compared to exact diagonalization ones in lattices with periodic boundary conditions. NLCEs are found to outperform exact diagonalization in periodic systems for all quantities studied.

arXiv:1701.02758•Physical review. E•Mar 7, 2017•20 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
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 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.
We develop a numerical linked cluster expansion (NLCE) method that can be applied directly to inhomogeneous systems, for example, Hamiltonians with disorder and dynamics initiated from inhomogeneous initial states. We demonstrate the method by calculating dynamics for single-spin expectations and spin correlations in two-dimensional spin models on a square lattice, starting from a checkerboard state. We show that NLCE can give moderate to dramatic improvement over an exact diagonalization of comparable computational cost and that the advantage in computational resources grows exponentially as the size of the clusters included grows. Although the method applies to any type of NLCE, our explicit benchmarks use the rectangle expansion. Besides showing the capability to treat inhomogeneous systems, these benchmarks demonstrate the rectangle expansion's utility out of equilibrium.
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.
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 discuss the application of numerical linked cluster expansions (NLCEs) to study one dimensional lattice systems in thermal equilibrium and after quantum quenches from thermal equilibrium states. For the former, we calculate observables in the grand canonical ensemble, and for the latter we calculate observables in the diagonal ensemble. When converged, NLCEs provide results in the thermodynamic limit. We use two different NLCEs: a maximally connected expansion introduced in previous works and a site-based expansion. We compare the effectiveness of both NLCEs. The site-based NLCE is found to work best for systems in thermal equilibrium. However, in thermal equilibrium and after quantum quenches, the site-based NLCE can diverge when the maximally connected one converges. We relate this divergence to the exponentially large number of clusters in the site-based NLCE and the behavior of the weights of observables in those clusters. We discuss the effectiveness of resummations to cure the divergence. Our NLCE calculations are compared to exact diagonalization ones in lattices with periodic boundary conditions. NLCEs are found to outperform exact diagonalization in periodic systems for all quantities studied.
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.