Pith. sign in

REVIEW 3 major objections 5 minor 46 references

Finite-Temperature Kinetic Ferromagnetism in the Square Lattice Hubbard Model

T0 review · 3 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The square-lattice Hubbard model hosts finite-temperature kinetic ferromagnetism in a doping window that widens with interaction strength, reaching about 33% in the infinite-U limit.

desk verdict A solid finite-temperature NLCE study of short-range kinetic ferromagnetism on the square lattice; the qualitative picture is convincing, but the quantitative boundaries are less secure than the abstract suggests. read the letter →

arxiv 2502.07252 v2 pith:3PQZLUCK submitted 2025-02-11 cond-mat.str-el cond-mat.quant-gas

Robin C. Newby , Ehsan Khatami This is my paper · ORCID
classification cond-mat.str-elcond-mat.quant-gas
keywords NagaokaferromagnetismHubbardmodelt-Jnumericallinked-clusterexpansionfinitetemperaturespincorrelationspolaronsquarelattice
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper is trying to establish that Nagaoka's kinetic ferromagnetism, usually proved for a single hole in an infinite-$U$ Hubbard band at zero temperature, survives as a finite-temperature, finite-doping phenomenon on the square lattice. Using numerical linked-cluster expansions (NLCE), the authors find that for $U/t$ between 80 and 240 and temperatures down to $T/t \sim 0.3$, the nearest-neighbor spin correlations turn ferromagnetic in a doping interval $\delta_{c1} < \delta < \delta_{c2}$, with $\delta_{c2}$ growing with $U/t$ toward about 33% and $\delta_{c1}$ shrinking toward zero. They also identify Nagaoka polarons, ferromagnetic bubbles bound to holes, persisting to roughly 20% doping. A reader should care because this maps where cold-atom square-lattice quantum simulators could look for kinetic magnetism, and it complements recent triangular-lattice observations by showing that a bipartite geometry still hosts the effect, just over a finite doping window.

What carries the argument

The numerical linked-cluster expansion (NLCE) is the central tool: it expresses an extensive property of the infinite lattice as a sum over all topologically distinct connected clusters, with inclusion-exclusion weights, and computes each cluster exactly by diagonalization, up to 9 sites for the Hubbard model and 11 sites for the $t$-$J$ model. Wynn's epsilon algorithm resums the site series to reach lower temperatures, and results are accepted when successive resummation cycles differ by less than 0.003. The order parameters are normalized correlation functions: $C_{ss}$, the two-point nearest-neighbor spin correlation normalized by the squared fraction of singly occupied sites, and $C_{hss}$, the three-point hole-spin-spin correlation normalized by the hole fraction times the squared singly occupied fraction. The zero crossings of $C_{ss}$ as a function of doping define $\delta_{c1}$ and $\delta_{c2}$, the boundaries of the ferromagnetic window.

What would settle it

A determinant quantum Monte Carlo calculation or a square-lattice ultracold-atom experiment at $U/t \approx 120$, $T/t \approx 0.4$, and hole doping near 18% should measure a positive nearest-neighbor spin correlation $C_{ss}$ if the claim is right; a negative value would falsify the ferromagnetic window. The same measurement at 35% doping in the $U\to\infty$ limit would test the predicted $\delta_{c2} \sim 33\%$ boundary, since a positive $C_{ss}$ there would contradict it.

Watch

Extended reading notes

Core claim

The paper argues that, at interaction strengths $U/t \gtrsim 80$ and temperatures as low as $T/t \sim 0.3$, the square-lattice Hubbard model and its strong-coupling $t$-$J$ limit develop a window of hole doping, $\delta_{c1} < \delta < \delta_{c2}$, in which the nearest-neighbor spin correlation $C_{ss}$ is ferromagnetic. The window widens with increasing $U/t$; in the $U \to \infty$ ($J \to 0$) limit $C_{ss}$ is ferromagnetic for any infinitesimal doping and the upper boundary approaches $\delta_{c2} \sim 33\%$. The three-point hole-spin-spin correlation $C_{hss}$ shows ferromagnetic bubbles around dopants, Nagaoka polarons, that persist to about 20% doping, and the polarons grow as temperature drops. The central claim is that kinetic ferromagnetism on a bipartite lattice is not confined to the single-hole, zero-temperature Nagaoka limit, but appears as a finite-temperature, finite-doping short-range correlated phase.

Load-bearing premise

The load-bearing assumption is that the truncated cluster expansion, clusters of up to 9 sites for the Hubbard model and 11 for the t-J model resummed with Wynn's algorithm, correctly represents the infinite lattice at the reported temperatures, with the only evidence being agreement between successive resummation cycles to within 0.003.

Editorial extensions

If this is right

  • A square optical lattice filled with ultracold fermions at $U/t \sim 100$-$240$ and $T/t \sim 0.4$ should show positive nearest-neighbor $C_{ss}$ for hole dopings roughly between a few percent and 20-30%, with the window shrinking as temperature rises.
  • The $t$-$J$ model with $J=0$ predicts $\delta_{c1}=0$ and $\delta_{c2}\sim 1/3$ in the $U\to\infty$ ground state, giving a concrete target for ground-state methods to confirm or refute.
  • Extrapolating the finite-temperature boundary $\delta_{c2}$ to zero yields a lower bound $U/t \sim 20$ for nearest-neighbor ferromagnetic correlations, so experiments below that interaction should see none.
  • Because the ferromagnetic correlations remain short-ranged at these temperatures, the predicted phase is a short-range correlated state rather than a fully polarized ferromagnet; detecting it requires nearest-neighbor correlation measurements, not bulk magnetization alone.
  • The close agreement between Hubbard and $t$-$J$ results at small doping validates using the $t$-$J$ model as a computational shortcut for strong-coupling kinetic-magnetism physics.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The 0.003 threshold between Wynn cycles is an internal consistency check; if it were tightened, the accessible temperature range would shrink and the reported $\delta_{c1}$ values at low temperature might shift upward, so the extrapolation $\delta_{c1}\to 0$ at zero temperature is the most fragile part of the picture.
  • A direct extension would be to compute the same correlators on a tilted or weakly anisotropic square lattice to test whether the finite ferromagnetic window is specific to the perfect bipartite square geometry or survives mild asymmetry.
  • The $J=0$ $t$-$J$ results can be read as a parameter-free prediction: at $U\to\infty$, $C_{ss}$ should be positive for any infinitesimal doping at the lowest temperatures, a statement an optical-lattice experiment with tunable interaction could test directly by measuring the sign of nearest-neighbor correlations at very low doping.
  • The apparent lower bound $U/t \sim 20$ suggests the kinetic-ferromagnetism window is an intermediate-to-strong-coupling phenomenon; a natural follow-up is to search for how the window depends on temperature just above that threshold.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This manuscript uses numerical linked-cluster expansions (NLCE) up to 9 sites for the Hubbard model and 11 sites for the t-J model to compute finite-temperature nearest-neighbor spin correlations Css and hole-spin-spin correlations Chss on the square lattice. The authors identify a doping window, bounded by critical dopings delta_c1 and delta_c2, within which Css is ferromagnetic, and they study its temperature and interaction dependence for U/t between 60 and 240 (and the J=0 t-J limit). They report evidence for finite-temperature Nagaoka polarons and conclude that ferromagnetic nearest-neighbor correlations survive in the thermodynamic limit up to dopings of about 33% in the U to infinity limit. The paper includes a comparison with recent cold-atom data on the triangular lattice and provides a public data repository.

Significance. If the quantitative boundaries are correct, this is a useful finite-temperature extension of the Nagaoka problem and provides concrete doping and temperature targets for cold-atom experiments in square optical lattices. The study is largely parameter-free: the correlation functions are computed directly from the Hubbard and t-J Hamiltonians, and no fitted parameters enter the central claim. The use of an established cluster method, the comparison with t-J results, the connection to existing variational and DMRG ground-state literature, and the deposition of data in Zenodo are all strengths. However, the quantitative deliverables—delta_c1, delta_c2, and the U to infinity extrapolation—rest on an internal Wynn-resummation consistency threshold and an incompletely specified outlier-removal step, which need to be addressed before the central finite-temperature claim is fully secured.

major comments (3)
  1. [Section III.A, Figs. 1-4] The convergence criterion that successive Wynn resummation cycles differ by less than 0.003 is an internal consistency check, not a systematic error estimate. Near the delta_c1 and delta_c2 crossings in Figs. 1-4, Css has magnitude roughly 0.01-0.02, so the 0.003 threshold is a substantial fraction of the signal; a modest resummation bias could move the zero crossings by several percent in doping and shift the phase boundaries plotted in Fig. 4. The manuscript should provide a quantitative convergence study, for example by comparing site-expansion orders, varying the acceptance threshold, or cross-checking against an independent method such as determinant quantum Monte Carlo at T/t around 0.4-0.5 for U/t = 80-120.
  2. [Fig. 1 caption] The caption states that 'A small number of outliers at various temperatures are removed' without specifying how many outliers, at which temperatures, or whether the removal was sign-blind. If outlier removal is correlated with the sign or magnitude of Css, the extracted delta_c1 and delta_c2 values in Fig. 4 could be systematically biased. Please report the exact removal criteria or show the raw data alongside the cleaned data.
  3. [Fig. 4, Conclusions] The quantitative claims that delta_c2 approaches about 33% in the U to infinity limit and that there is a lower bound of U/t around 20 are read off from the gray crosses and the dashed linear-regression line in Fig. 4, but no uncertainties or regression parameters are given for this extrapolation. The only error bars mentioned are from the mu-T interpolation and are smaller than the symbols; they do not include NLCE truncation or resummation uncertainty. Please propagate the convergence uncertainty into delta_c1 and delta_c2 and into the extrapolated quantities, or state explicitly that the gray crosses and linear regression are only qualitative guides.
minor comments (5)
  1. [Abstract vs. Conclusions] The abstract says temperatures as low as 0.2 of the hopping energy, but Section V states T/t ~ 0.3 and the data in Fig. 1 appear to begin near 0.3-0.4; please reconcile this inconsistency.
  2. [Fig. 6 caption] The sentence 'Here, we do not show data for which the difference between the two resummations is smaller than 0.003' appears to reverse the criterion used elsewhere in the paper; it should presumably read 'larger than 0.003'.
  3. [Fig. 7 caption] There is a typo: 'neaerest-neighbor' should be 'nearest-neighbor'.
  4. [Fig. 10 caption] The caption refers to yellow bonds placed where 'the agreement between Wynn orders exceeds our error thresholds,' but the threshold value is not defined in the text; please specify it or refer explicitly to Section III.A.
  5. [Reference [37]] Reference [37] is listed as '(In preparation)'; if the cluster enumeration algorithm is important for reproducibility, please provide a preprint identifier or DOI when available.

Circularity Check

0 steps flagged · score 0.0 of 10

No material circularity: the central correlation functions are computed directly from the Hubbard and t-J Hamiltonians via NLCE; self-citations serve only as method benchmarks.

full rationale

The paper's central quantities, Css and Chss, are defined in Eqs. (9)-(10) as expectation values of products of spin and hole operators, evaluated by exact diagonalization of finite clusters and combined through the inclusion-exclusion weights of Eqs. (3)-(4). No parameter is fitted to the FM window, δc1 or δc2; those boundaries are read off from the sign of the computed Css. The only regression in the paper, the dashed gray line in Fig. 4, is explicitly labeled a guide for the eye and is not used to derive the U→∞ δc2∼1/3 statement, which instead comes from the J=0 t-J calculation. Self-citations (Refs. 32-33 for the NLCE formalism and Refs. 39-42 for prior benchmarks) validate the method against independent DQMC and experimental results rather than supplying the target physics. The Wynn-resummation agreement threshold of 0.003 is an internal consistency criterion, not a fitted input, and while the unspecified outlier removal and the abstract/conclusion temperature mismatch (0.2 vs 0.3) are legitimate correctness concerns, they are not instances of circular reasoning. No step in the derivation chain reduces to its own input, so the circularity score is 0.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The central claim rests on numerical convergence assumptions and on the use of short-range correlation signs as a proxy for ferromagnetic correlations, not on fitted physical parameters. No new particles, forces, or conserved quantities are introduced; the Nagaoka polaron is an existing concept used as an interpretive label. The Wynn threshold is an operational numerical tolerance, not a physical free parameter.

assumptions (5)
  • domain assumption The NLCE site expansion at order 9 (Hubbard) and 11 (t-J), together with Wynn resummation, converges to the thermodynamic limit for the reported temperatures and densities.
    Invoked in Section III.A; the convergence check is internal consistency between Wynn cycles rather than an external exact result.
  • ad hoc to paper The difference between successive Wynn resummation cycles being below 0.003 is a sufficient error criterion for all reported quantities.
    Used throughout to decide which data to show; the threshold is chosen by the authors and not independently justified.
  • domain assumption The sign of the normalized nearest-neighbor spin correlation Css is a valid proxy for the existence of Nagaoka ferromagnetic correlations.
    Central to defining delta_c1 and delta_c2 in Section IV; a short-range correlation sign is not a phase transition marker.
  • standard math The t-J model with J=4t^2/U accurately represents the strong-coupling limit of the Hubbard model for the U/t values considered.
    Standard strong-coupling expansion, used in Section II.B and in comparisons in Figs. 2-4.
  • standard math Particle-hole symmetry on the square lattice maps doublon-spin-spin correlations to hole-spin-spin correlations.
    Used in Section III.B to relate Cdss and Chss.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Finite-Temperature Kinetic Ferromagnetism in the Square Lattice Hubbard Model." pith.science (2026). https://pith.science/paper/3PQZLUCK

@misc{pith2026250207252,
  author       = {Pith},
  title        = {Pith review of: Finite-Temperature Kinetic Ferromagnetism in the Square Lattice Hubbard Model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3PQZLUCK}},
  note         = {Machine review of arXiv:2502.07252}
}
abstract

While the exact phase diagram of the Fermi-Hubbard model remains poorly understood despite decades of progress, nearly 60 years ago, Nagaoka proved that a single dopant in an otherwise half-filled Hubbard system can bring about ferromagnetism through kinetic means. The phenomenon was recently observed with ultracold atoms in triangular optical lattices. Here, we explore the kinetic ferromagnetism within the square lattice Hubbard model and its strong-coupling counterpart, the $t-J$ model, at finite temperatures in the thermodynamic limit via numerical linked-cluster expansions. We find evidence of ferromagnetic Nagaoka polarons at dopings up to $\sim 30\%$ away from half filling for a variety of interaction strengths and at temperatures as low as $0.2$ of the hopping energy. We map out the boundaries of this phase through analyzing various correlation functions.

Figures

Figures reproduced from arXiv: 2502.07252 by the authors.

Figure 1
Figure 1. FIG. 1. The normalized nearest-neighbor spin correlations [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 4
Figure 4. FIG. 4. The critical dopings [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 3
Figure 3. FIG. 3. The normalized nearest-neighbor spin correlations [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (6 more)
Figure 5
Figure 5. Figure 5: FIG. 5. The normalized hole-spin-spin correlation function [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. The normalized hole-spin-spin correlation function of [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. The normalized neaerest-neighbor spin correlations [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Spatial map of [PITH_FULL_IMAGE:figures/full_fig_p008_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. The normalized hole-spin-spin correlation function of [PITH_FULL_IMAGE:figures/full_fig_p008_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10. Spatial map of [PITH_FULL_IMAGE:figures/full_fig_p009_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

46 extracted references · 34 canonical work pages

  1. [1]

    Keimer, S

    B. Keimer, S. A. Kivelson, M. R. Norman, S. Uchida, and J. Zaanen, From quantum matter to high-temperature superconductivity in copper oxides, Nature 518, 179 (2015)

    2015

  2. [2]

    D. P. Arovas, E. Berg, S. A. Kivelson, and S. Raghu, The Hubbard Model, Annual Review of Condensed Matter Physics 13, 239 (2022)

    2022

  3. [3]

    M. Qin, T. Sch¨ afer, S. Andergassen, P. Corboz, and E. Gull, The Hubbard Model: A Computational Per- spective, Annual Review of Condensed Matter Physics 13, 275 (2022)

    2022

  4. [4]

    Bohrdt, L

    A. Bohrdt, L. Homeier, C. Reinmoser, E. Demler, and F. Grusdt, Exploration of doped quantum magnets with ultracold atoms, Annals of Physics Special Issue on Philip W. Anderson, 435, 168651 (2021)

    work page 2021

  5. [5]

    Xu, C.-M

    H. Xu, C.-M. Chung, M. Qin, U. Schollw¨ ock, S. R. White, and S. Zhang, Coexistence of superconductivity with par- tially filled stripes in the Hubbard model, Science 384, eadh7691 (2024)

    work page 2024

  6. [6]

    Nagaoka, Ferromagnetism in a narrow, almost half- filled s band, Physical Review 147, 392 (1966)

    Y. Nagaoka, Ferromagnetism in a narrow, almost half- filled s band, Physical Review 147, 392 (1966)

    work page 1966

  7. [7]

    J. P. Dehollain, U. Mukhopadhyay, V. P. Michal, Y. Wang, B. Wunsch, C. Reichl, W. Wegscheider, M. S. Rudner, E. Demler, and L. M. K. Vandersypen, Nagaoka ferromagnetism observed in a quantum dot plaquette, Nature 579, 528 (2020)

    2020

  8. [8]

    Y. Tang, L. Li, T. Li, Y. Xu, S. Liu, K. Barmak, K. Watanabe, T. Taniguchi, A. H. MacDonald, J. Shan, and K. F. Mak, Simulation of Hubbard model physics in WSe2/WS2 moir´ e superlattices, Nature579, 353 (2020)

    2020

Show all 46 references
  1. [9]

    Ciorciaro, T

    L. Ciorciaro, T. Smole´ nski, I. Morera, N. Kiper, S. Hies- tand, M. Kroner, Y. Zhang, K. Watanabe, T. Taniguchi, E. Demler, and A. ˙Imamo˘ glu, Kinetic magnetism in tri- angular moir´ e materials, Nature623, 509 (2023)

    2023

  2. [10]

    Z. Tao, W. Zhao, B. Shen, T. Li, P. Kn¨ uppel, K. Watan- abe, T. Taniguchi, J. Shan, and K. F. Mak, Observation of spin polarons in a frustrated moir´ e Hubbard system, Nature Physics 20, 783 (2024)

    2024

  3. [11]

    Lebrat, M

    M. Lebrat, M. Xu, L. H. Kendrick, A. Kale, Y. Gang, P. Seetharaman, I. Morera, E. Khatami, E. Demler, and M. Greiner, Observation of Nagaoka polarons in a Fermi– Hubbard quantum simulator, Nature 629, 317 (2024)

    2024

  4. [12]

    M. L. Prichard, B. M. Spar, I. Morera, E. Demler, Z. Z. Yan, and W. S. Bakr, Directly imaging spin polarons in a kinetically frustrated Hubbard system, Nature 629, 323 (2024)

    2024

  5. [13]

    A. F. Andreev, Structure of vacancies in solid He3, JETP Lett. 24, 608 (1976)

    1976

  6. [14]

    S. R. White and I. Affleck, Density matrix renormaliza- tion group analysis of the Nagaoka polaron in the two- dimensional t − J model, Physical Review B 64, 024411 (2001)

    2001

  7. [15]

    Hanisch, B

    T. Hanisch, B. Kleine, A. Ritzl, and E. M¨ uller-Hartmann, Ferromagnetism in the Hubbard model: Instability of the Nagaoka state on the triangular, honeycomb and kagome lattices, Annalen der Physik 507, 303 (1995)

    1995

  8. [16]

    J. O. Haerter and B. S. Shastry, Kinetic Antiferromag- netism in the Triangular Lattice, Physical Review Letters 95, 087202 (2005)

    2005

  9. [17]

    Carlstr¨ om, In situ controllable magnetic phases in doped twisted bilayer transition metal dichalcogenides, Physical Review Research 4, 043126 (2022)

    J. Carlstr¨ om, In situ controllable magnetic phases in doped twisted bilayer transition metal dichalcogenides, Physical Review Research 4, 043126 (2022)

    2022

  10. [18]

    Doucot and X

    B. Doucot and X. G. Wen, Instability of the Nagaoka state with more than one hole, Physical Review B 40, 2719 (1989)

    1989

  11. [19]

    B. S. Shastry, H. R. Krishnamurthy, and P. W. Anderson, Instability of the nagaoka ferromagnetic state of the U = ∞ Hubbard model, Physical Review B 41, 2375 (1990)

    1990

  12. [20]

    W. O. Putikka, M. U. Luchini, and M. Ogata, Ferro- magnetism in the two-dimensional t - J model, Physical Review Letters 69, 2288 (1992)

    1992

  13. [21]

    Tian, The Nagaoka state in the one-hand Hubbard model with two and more holes, Journal of Physics A: Mathematical and General 24, 513 (1991)

    G.-S. Tian, The Nagaoka state in the one-hand Hubbard model with two and more holes, Journal of Physics A: Mathematical and General 24, 513 (1991)

    1991

  14. [22]

    X. Y. Zhang, E. Abrahams, and G. Kotliar, Quantum Monte Carlo algorithm for constrained fermions: Appli- cation to the infinite-U Hubbard model, Physical Review Letters 66, 1236 (1991)

    1991

  15. [23]

    W. V. D. Linden and D. M. Edwards, Ferromagnetism in the Hubbard model, Journal of Physics: Condensed Matter 3, 4917 (1991)

    1991

  16. [24]

    Becca and S

    F. Becca and S. Sorella, Nagaoka ferromagnetism in the two-dimensional infinite-U Hubbard model, Physical Re- view Letters 86, 3396 (2001)

    2001

  17. [25]

    Wurth, G

    P. Wurth, G. Uhrig, and E. M¨ uller-Hartmann, Ferromag- netism in the Hubbard model on the square lattice: Im- proved instability criterion for the Nagaoka state, An- nalen der Physik 508, 148 (1996)

    1996

  18. [26]

    S. Yun, W. Dobrautz, H. Luo, V. Katukuri, N. Lieber- mann, and A. Alavi, Ferromagnetic domains in the large- U Hubbard model with a few holes: A full configuration interaction quantum Monte Carlo study, Physical Review B 107, 064405 (2023)

    2023

  19. [27]

    Samajdar and R

    R. Samajdar and R. N. Bhatt, Polaronic mechanism of Nagaoka ferromagnetism in Hubbard models, Physical Review B 109, 235128 (2024)

    2024

  20. [28]

    Samajdar and R

    R. Samajdar and R. N. Bhatt, Nagaoka ferromagnetism in doped Hubbard models in optical lattices, Physical Review A 110, L021303 (2024)

    2024

  21. [29]

    T. J. Harris, U. Schollw¨ ock, A. Bohrdt, and F. Grusdt, Kinetic magnetism and stripe order in the antiferro- magnetic bosonic t − J model (2024), arXiv:2410.00904 [cond-mat.quant-gas]

    2024 arXiv

  22. [30]

    Morera and E

    I. Morera and E. Demler, Itinerant magnetism and mag- netic polarons in the triangular lattice Hubbard model (2024), arXiv:2402.14074 [cond-mat.str-el]

    2024 arXiv

  23. [31]

    Q. Chen, S. A. Chen, and Z. Zhu, Ferromagnetism Mech- anism in a Geometrically Frustrated Triangular Lattice (2024), arXiv:2408.05971 [cond-mat]

    2024

  24. [32]

    Rigol, T

    M. Rigol, T. Bryant, and R. R. P. Singh, Numerical linked-cluster approach to quantum lattice models, Phys. Rev. Lett. 97, 187202 (2006)

    2006

  25. [33]

    B. Tang, E. Khatami, and M. Rigol, A short introduc- tion to numerical linked-cluster expansions, Computer Physics Communications 184, 557 (2013)

    2013

  26. [34]

    Hubbard, Electron Correlations in Narrow Energy Bands, Proceedings of the Royal Society of London

    J. Hubbard, Electron Correlations in Narrow Energy Bands, Proceedings of the Royal Society of London. Se- ries A, Mathematical and Physical Sciences 276, 238 (1963)

    1963

  27. [35]

    M. F. Sykes, J. W. Essam, B. R. Heap, and B. J. Hiley, Lattice constant systems and graph theory, Journal of Mathematical Physics 7, 1557 (1966). 11

    1966

  28. [36]

    B. D. McKay and A. Piperno, Practical graph isomor- phism, II, Journal of Symbolic Computation 60, 94 (2014)

    2014

  29. [37]

    Seetharaman and E

    P. Seetharaman and E. Khatami, A versatile algorithm for numerical linked-cluster expansions (2025), (In prepa- ration)

    2025

  30. [38]

    Wynn, On the Convergence and Stability of the Ep- silon Algorithm, SIAM Journal on Numerical Analysis 3, 91 (1966)

    P. Wynn, On the Convergence and Stability of the Ep- silon Algorithm, SIAM Journal on Numerical Analysis 3, 91 (1966)

    1966

  31. [39]

    Khatami and M

    E. Khatami and M. Rigol, Thermodynamics of strongly interacting fermions in two-dimensional optical lattices, Phys. Rev. A 84, 053611 (2011)

    2011

  32. [40]

    B. Tang, T. Paiva, E. Khatami, and Rigol, Finite- temperature properties of strongly correlated fermions in the honeycomb lattice, Phys. Rev. B 88, 125127 (2013)

    2013

  33. [41]

    Khatami, R

    E. Khatami, R. T. Scalettar, and R. R. P. Singh, Finite- temperature superconducting correlations of the Hub- bard model, Phys. Rev. B(R) 91, 241107 (2015)

    2015

  34. [42]

    L. W. Cheuk, M. A. Nichols, K. R. Lawrence, M. Okan, H. Zhang, E. Khatami, N. Trivedi, T. Paiva, M. Rigol, and M. W. Zwierlein, Observa- tion of spatial charge and spin correlations in the 2d Fermi-Hubbard model, Science 353, 1260 (2016), http://science.sciencemag.org/content/...

    2016

  35. [43]

    C. N. Varney, C.-R. Lee, Z. J. Bai, S. Chiesa, M. Jarrell, and R. T. Scalettar, Quantum Monte Carlo study of the two-dimensional fermion Hubbard model, Phys. Rev. B 80, 075116 (2009)

    2009

  36. [44]

    T. A. Hilker, G. Salomon, F. Grusdt, A. Omran, M. Boll, E. Demler, I. Bloch, and C. Gross, Revealing hidden an- tiferromagnetic correlations in doped Hubbard chains via string correlators, Science 357, 484 (2017)

    2017

  37. [45]

    Omran, M

    A. Omran, M. Boll, T. A. Hilker, K. Kleinlein, G. Sa- lomon, I. Bloch, and C. Gross, Microscopic Observation of Pauli Blocking in Degenerate Fermionic Lattice Gases, Physical Review Letters 115, 263001 (2015)

    2015

  38. [46]

    Data plotted in figures of this manuscript can be found at https://doi.org/10.5281/zenodo.15495644

    doi

Pith tools

Reviewed August 8, 2026 · model on record in the stance chip above.