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:1501.00990 · Physical Review B

Thermodynamics of two-dimensional spin models with bimodal random-bond disorder

Baoming Tang, Deepak Iyer, Marcos Rigol

Published May 14, 2015• 18 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.91.174413

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 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 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

Thermodynamic properties of the two-dimensional quantum Heisenberg ferromagnet and the effects of bond dilution

E. F. Talbot

The energy, the specific heat, and the magnetic susceptibility of the bond-diluted S=(1/2) Heisenberg ferromagnet on a two-dimensional square lattice are studied using a Monte Carlo sampling of the terms contributing to the partition function. The full lattice is studied in detail for sizes up to 20\ifmmode\times\else\texttimes\fi{}20, and evidence is presented to suggest that the low-temperature correlation length behaves as \ensuremath{\xi}\ensuremath{\sim}${e}^{\ensuremath{\nu}\ensuremath{\beta}J}$ with \ensuremath{\nu}\ensuremath{\sim}1.91\ifmmode\pm\else\textpm\fi{}0.06. The temperature dependence of the susceptibility, however, is not quite a simple exponential in \ensuremath{\beta}J. Knowledge of the correlation length and finite-size scaling allows the determination of the prefactor. Our conclusion is that the low-temperature susceptibility behaves as ${T}^{\mathrm{\ensuremath{-}}1}$${e}^{2\ensuremath{\nu}\ensuremath{\beta}J}$. Square lattices of size ranging from 6\ifmmode\times\else\texttimes\fi{}6 to 12\ifmmode\times\else\texttimes\fi{}12 are studied for ${k}_{B}$T/J ranging from 0.1 to 3 and for several values of the bond probability p. Not surprisingly, the most drastic thermodynamic behavior near the critical bond-percolation probability is the sharp increase in the low-temperature susceptibility for p>${p}_{c}$. This corresponds to much larger correlation effects in the percolating lattice than below the percolation threshold. We also compare the thermodynamics at p=${p}_{c}$ to the scaling theories of Lubensky and Stanley et al. which suggest that the correlations at the percolation threshold are predominantly one dimensional.

Physical review. B, Condensed matter•Apr 1, 1986•5 citations

Properties of the two-dimensional spin-1/2 Heisenberg model on a honeycomb lattice with interlayer coupling

Löw

The magnetic properties of the two-dimensional S = 1 2 (quantum) antiferromagnetic Heisenberg model on a honeycomb lattice with and without interlayer coupling are studied by means of a continuous Euclidean time Quantum-Monte-Carlo algorithm.The internal energy, the magnetic susceptibility and the staggered magnetization are determined in the full temperature range.For the two-dimensional system the groundstate energy/bond is found to be E hc 0 = -0.36303(13),and the zero temperature staggered magnetization m st = 0.2681(8).For coupled planes of honeycomb systems a phase transition from an ordered phase to a disordered phase is found at T /J = 0.695(10).

Condensed Matter Physics•Jan 1, 2009•32 citations

Thermodynamic properties of disordered quantum spin ladders

Erol Vatansever, Georgi Georgiev Grahovski, Nikolaos G. Fytas

Abstract In this paper, we study the thermodynamic properties of spin-1/2 antiferromagnetic Heisenberg ladders by means of the stochastic series expansion quantum Monte Carlo technique. This includes the thermal properties of the specific heat, uniform and staggered susceptibilities, spin gap, and structure factor. Our numerical simulations are probed over a large ensemble of random realizations in a wide range of disorder strengths r, from the clean ( $$r=0$$ r = 0 ) case up to the diluted ( $$r \rightarrow 1$$ r → 1 ) limit, and for selected choices of number of legs $$L_y$$ L y per site. Our results show some interesting phenomena, like the presence of crossing points in the temperature plane for both the specific heat and uniform susceptibility curves which appear to be universal in r, as well as a variable dependence of the spin gap in the amount of disorder upon increasing $$L_y$$ L y . Graphical Abstract

arXiv:2312.12834•The European Physical Journal B•Mar 1, 2024•1 citation
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±\else\textpm\fiJ\ifmmode\pm\else\textpm\fi{}J\ifmmode±\else\textpm\fiJ 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.
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 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.
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.
The energy, the specific heat, and the magnetic susceptibility of the bond-diluted S=(1/2) Heisenberg ferromagnet on a two-dimensional square lattice are studied using a Monte Carlo sampling of the terms contributing to the partition function. The full lattice is studied in detail for sizes up to 20\ifmmode\times\else\texttimes\fi{}20, and evidence is presented to suggest that the low-temperature correlation length behaves as \ensuremath{\xi}\ensuremath{\sim}${e}^{\ensuremath{\nu}\ensuremath{\beta}J}$ with \ensuremath{\nu}\ensuremath{\sim}1.91\ifmmode\pm\else\textpm\fi{}0.06. The temperature dependence of the susceptibility, however, is not quite a simple exponential in \ensuremath{\beta}J. Knowledge of the correlation length and finite-size scaling allows the determination of the prefactor. Our conclusion is that the low-temperature susceptibility behaves as ${T}^{\mathrm{\ensuremath{-}}1}$${e}^{2\ensuremath{\nu}\ensuremath{\beta}J}$. Square lattices of size ranging from 6\ifmmode\times\else\texttimes\fi{}6 to 12\ifmmode\times\else\texttimes\fi{}12 are studied for kB{k}_{B}kB​T/J ranging from 0.1 to 3 and for several values of the bond probability p. Not surprisingly, the most drastic thermodynamic behavior near the critical bond-percolation probability is the sharp increase in the low-temperature susceptibility for p>pc{p}_{c}pc​. This corresponds to much larger correlation effects in the percolating lattice than below the percolation threshold. We also compare the thermodynamics at p=pc{p}_{c}pc​ to the scaling theories of Lubensky and Stanley et al. which suggest that the correlations at the percolation threshold are predominantly one dimensional.
The magnetic properties of the two-dimensional S = 1 2 (quantum) antiferromagnetic Heisenberg model on a honeycomb lattice with and without interlayer coupling are studied by means of a continuous Euclidean time Quantum-Monte-Carlo algorithm.The internal energy, the magnetic susceptibility and the staggered magnetization are determined in the full temperature range.For the two-dimensional system the groundstate energy/bond is found to be E hc 0 = -0.36303(13),and the zero temperature staggered magnetization m st = 0.2681(8).For coupled planes of honeycomb systems a phase transition from an ordered phase to a disordered phase is found at T /J = 0.695(10).
Abstract In this paper, we study the thermodynamic properties of spin-1/2 antiferromagnetic Heisenberg ladders by means of the stochastic series expansion quantum Monte Carlo technique. This includes the thermal properties of the specific heat, uniform and staggered susceptibilities, spin gap, and structure factor. Our numerical simulations are probed over a large ensemble of random realizations in a wide range of disorder strengths r, from the clean ( r=0r=0r=0 r = 0 ) case up to the diluted ( r→1r \rightarrow 1r→1 r → 1 ) limit, and for selected choices of number of legs LyL_yLy​ L y per site. Our results show some interesting phenomena, like the presence of crossing points in the temperature plane for both the specific heat and uniform susceptibility curves which appear to be universal in r, as well as a variable dependence of the spin gap in the amount of disorder upon increasing LyL_yLy​ L y . Graphical Abstract