The harness's model zoo: 29 canonical quantum lattice models, each with its Hamiltonian, a structured A1–D16 property table, phases, canonical observables, and the methods that work on it. Every row is generated from the model's knowledge-base card.
nearest-neighbor quantum magnets — critical chains, Néel order, Haldane physics
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain / quasi-1D ladder / 2D square (Z=4), triangular, kagome / 3D | Geometry is the master axis: it sets entanglement scaling and frustration. |
| A2 boundary conditions | OBC (DMRG default) · PBC (ED/QMC) · cylinder (2D DMRG) | Cylinder width caps the 2D-DMRG bond-dim budget. |
| A3 statistics & local dim | spin-S; d = 2S+1 (d=2 for S=1/2) | Maps to hard-core bosons via Matsubara–Matsuda; no statistics sign by itself. |
| A4 interaction range | short-range (nearest-neighbor) | Local — area-law-compatible. |
| B5 entanglement scaling | 1D gapless: area+log, S=(c/3)log ℓ, c=1 · 2D ordered: area law (∝ L) | 1D AFM chain is critical (c=1); 2D Néel state obeys area law. |
| B6 spectral gap | 1D S=1/2 AFM: gapless (power-law correlations) · 2D square AFM: gapless Goldstone magnons over ordered GS | Half-integer 1D chain gapless (LSM); 2D has long-range order + gapless spin waves. |
| B7 ground-state order | 1D S=1/2: quasi-long-range (no SSB, Mermin–Wagner) · 2D square/3D: Néel SSB · frustrated 2D: candidate spin liquid / VBS | Bipartite unfrustrated lattices order (2D/3D); 1D and frustrated cases do not. |
| B8 frustration | none on bipartite (chain, square, cubic) · geometric on triangular/kagome/pyrochlore | Frustration is what turns on the QMC sign problem for spins. |
| C9 global symmetry | SU(2) (total S), U(1) (S^z_tot) | Full SU(2) at the isotropic point; a field breaks SU(2)→U(1). |
| C10 spatial symmetry | translation (k), point group (D_4 square, D_6 triangular) | Block-diagonalizes ED sectors. |
| C11 integrability | 1D S=1/2 chain: Bethe-ansatz integrable · 2D / S≥1 / frustrated: non-integrable | Exact 1D thermodynamics via TBA; everything else numerical. |
| C12 sign problem | sign-free on bipartite (Marshall rule) · sign-ful on frustrated lattices | Unfrustrated → exact QMC at scale; frustrated → QMC blocked. |
| D13 regime | ground state (T=0) default; finite-T and dynamics out of card scope | GS energy / order parameter are the canonical targets. |
| D14 filling / doping | N/A (spin model, no charge) | Doping appears only after fermionization (t-J / Hubbard). |
| D15 disorder | clean (translation-invariant) by default | Bond/site disorder → random-singlet physics (out of scope). |
| D16 hermiticity | Hermitian / closed | — |
⟨S_0·S_r⟩ ∼ (-1)^r / r (log corrections).m_s and structure-factor peak S(π,π).spin-liquid.E/N.m_s / sublattice magnetization (order parameter).S(q), peaked at the ordering wavevector.⟨S_i·S_j⟩; central charge c (1D, from entanglement scaling).method-property-map.md §MPS).N): sign-free QMC (SSE) — bipartite Marshall sign rule makes it exact at scale (§QMC, C12).E/N = 1/4 − ln 2 ≈ −0.443147 — exact Bethe ansatz (Hulthén 1938; convention H = J Σ S_i·S_j, J=1).E/N ≈ −0.6694 (high-precision −0.669441857(7)), staggered magnetization m_s ≈ 0.3074 — QMC/SSE (Sandvik, Phys. Rev. B 56, 11678 (1997)).Key reference: manousakis_1991_spin — the authoritative review of the spin-½ square-lattice Heisenberg antiferromagnet (spin-wave, Schwinger boson, series, QMC, ED) and its connection to the cuprate parents.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain (Z=2) | The defining 1D integrable spin chain. |
| A2 boundary conditions | PBC (Bethe ansatz / ED) · OBC (DMRG default) | PBC needed for the clean Bethe-ansatz spectrum and momentum sectors. |
| A3 statistics & local dim | spin-1/2; d = 2 | Jordan–Wigner maps to spinless fermions with NN hopping + NN interaction (the t-V chain). |
| A4 interaction range | short-range (nearest-neighbor) | Local — area-law compatible. |
| B5 entanglement scaling | critical phase (−1<Δ≤1): area+log, S=(c/3)log ℓ, c=1 · gapped phases (Δ<−1, Δ>1): area law (constant S) | Luttinger-liquid central charge c=1 in the whole XY/critical window. |
| B6 spectral gap | gapless (−1<Δ≤1, power-law correlations) · gapped FM (Δ<−1) · gapped Ising-Néel AFM (Δ>1) | Gap opens exponentially as Δ→1⁺ (BKT); FM transition at Δ=−1 is first-order. |
| B7 ground-state order | gapless XY/Luttinger liquid (quasi-long-range, no SSB) · ferromagnet (Δ<−1, fully polarized SSB) · Ising-Néel AFM (Δ>1, staggered SSB) | Quasi-LRO in the critical phase (Mermin–Wagner forbids true LRO in 1D). |
| B8 frustration | none (bipartite chain) | NN-only on a bipartite lattice; the AFM case is sign-free (Marshall). |
| C9 global symmetry | U(1) (S^z_tot); full SU(2) only at Δ=1; Z_2 spin-flip | Anisotropy Δ≠1 breaks SU(2)→U(1). |
| C10 spatial symmetry | translation (k), inversion/parity | Block-diagonalizes ED sectors; CDW/Néel image breaks translation by one site. |
| C11 integrability | Bethe-ansatz integrable for all Δ | Exact spectrum, ground-state energy, gap, and thermodynamics (TBA) at every anisotropy. |
| C12 sign problem | sign-free for the bipartite AFM (Marshall sign rule) | Unfrustrated bipartite chain → QMC-exact; DMRG is the practical workhorse. |
| D13 regime | ground state (T=0) default; finite-T accessible via TBA / LTRG (out of card scope) | E/N, gap, c, correlators are the canonical targets. |
| D14 filling / doping | N/A (spin model; S^z_tot plays the role of magnetization/filling) | A longitudinal field tunes magnetization, staying integrable. |
| D15 disorder | clean (translation-invariant) by default | Bond disorder → random-singlet phase (out of scope). |
| D16 hermiticity | Hermitian / closed | — |
−1 < Δ ≤ 1) : quasi-long-range order, no SSB; spin–spin correlations decay as a power law with Δ-dependent exponent; central charge c = 1.Δ > 1) : staggered magnetization m_s = (1/N)Σ_i (−1)^i ⟨S^z_i⟩, structure-factor peak S(π), finite gap. Onset is a BKT transition at Δ = 1.Δ < −1) : fully polarized ⟨S^z_i⟩ = ±1/2. The transition at Δ = −1 is first-order (ground-state level crossing).E/N.Δ>1); central charge c and Luttinger parameter K (critical phase, from entanglement / correlation scaling).S(q) (peaks at q=π in the Néel phase); spin–spin correlators ⟨S^z_0 S^z_r⟩, ⟨S^+_0 S^-_r⟩.S^z conservation; converges fast in the gapped phases, χ-hungry but reliable in the critical phase (per method-property-map.md §MPS, B5).Δ (C11); sign-free QMC/SSE for the AFM side at scale.Δ = 1 (Heisenberg, PBC, thermodynamic limit): E/N = 1/4 − ln 2 ≈ −0.443147 — exact Bethe ansatz (convention H = J Σ S_i·S_j, J = 1).Δ = 0 (free fermions via Jordan–Wigner, half-filling, thermodynamic limit): E/N = −1/π ≈ −0.318310 (convention J = 1, in-plane terms only).Key reference: franchini_2016_introduction — pedagogical all-details monograph on integrable techniques; works the XXZ chain explicitly from the coordinate Bethe ansatz to the algebraic Bethe ansatz, with the ground-state energy, the Δ-tuned phase structure, and finite-temperature thermodynamics.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain (Z=2) | Quasi-1D ladders extend the family but soften the topological order. |
| A2 boundary conditions | OBC (DMRG default; exposes edge spins) · PBC (clean entanglement-spectrum cut) | OBC reveals the S=1/2 Haldane edge modes. |
| A3 statistics & local dim | spin-1; d = 3 | Larger d than S=1/2 → higher per-site MPS/ED cost. |
| A4 interaction range | short-range (nearest-neighbor) + on-site D | Local. |
| B5 entanglement scaling | Haldane phase: area law (constant S), entanglement spectrum doubly degenerate | SPT signature lives in the entanglement spectrum, not a local order parameter. |
| B6 spectral gap | Haldane phase: gapped (Haldane gap) · transitions (Ising/Gaussian): gap closes | Integer-spin chain is gapped — Haldane's conjecture. |
| B7 ground-state order | Haldane = symmetry-protected topological (SPT) · large-D trivial · Néel (large Δ) SSB | Protected by Z_2×Z_2 / inversion / time reversal; diagnose via string order. |
| B8 frustration | none (default) | NNN coupling could add competition (out of default scope). |
| C9 global symmetry | U(1) (S^z_tot); full SU(2) only at Δ=1, D=0; Z_2×Z_2 (π-rotations) protects the Haldane SPT | D breaks SU(2)→U(1) even at Δ=1. |
| C10 spatial symmetry | translation, inversion (protects SPT), reflection | Inversion is one of the protecting symmetries. |
| C11 integrability | non-integrable (S=1 Heisenberg chain not Bethe-solvable; cf. exactly-solvable AKLT point with added biquadratic term) | The pure Heisenberg S=1 chain has no exact solution; AKLT is a nearby solvable model. |
| C12 sign problem | sign-free (unfrustrated bipartite chain → QMC applicable; DMRG is the workhorse) | No frustration → no spin sign problem. |
| D13 regime | ground state (T=0) default | Gap, string order, entanglement spectrum are the targets. |
| D14 filling / doping | N/A (spin model) | — |
| D15 disorder | clean by default | — |
| D16 hermiticity | Hermitian / closed | — |
O_string^z = −lim_{|i−j|→∞} ⟨S_i^z exp(iπ Σ_{i<k<j} S_k^z) S_j^z⟩; doubly-degenerate entanglement spectrum; fractional S=1/2 edge spins (OBC).D trivial : product-like Π|S^z=0⟩; string order vanishes; entanglement spectrum non-degenerate.Δ) : staggered magnetization, conventional SSB.E/N; Haldane gap Δ_H.method-property-map.md B7-SPT row).L ≲ 14 for exact gap/spectrum; TEBD imaginary-time route; sign-free QMC also available (unfrustrated).Δ=1, D=0): Δ_H/J ≈ 0.41048(6) — DMRG/QMC (White & Huse, Phys. Rev. B 48, 3844 (1993); Todo & Kato 2001).E/N ≈ −1.401484039 per spin (White & Huse DMRG, same H with in-plane coupling 1).Key reference: wierschem_2014_characterizing — concise review of the Haldane phase in spin-1 Heisenberg antiferromagnets: string order, SPT classification, entanglement spectrum, and quasi-1D phase diagram.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain (Z=2) · 2D square (Z=4) · higher-D | The canonical 1D quantum-critical model; 2D is Wilson–Fisher 3D-Ising universality. |
| A2 boundary conditions | OBC (DMRG) · PBC (ED, free-fermion) · cylinder (2D) | PBC matters for the exact Jordan–Wigner mapping. |
| A3 statistics & local dim | spin-1/2; d = 2 | Maps to free fermions in 1D via Jordan–Wigner. |
| A4 interaction range | short-range (nearest-neighbor); long-range 1/r^α variant | Long-range version inflates bond dimension (use TDVP). |
| B5 entanglement scaling | gapped phases: area law (const in 1D) · 1D critical Γ=J: area+log, c=1/2 | c=1/2 is the Ising-CFT central charge (one Majorana). |
| B6 spectral gap | gapped FM (Γ<J) and PM (Γ>J) · gapless at the QCP Γ=J (1D) | Quantum critical point separates ordered and disordered phases. |
| B7 ground-state order | FM (Γ<J): Z_2 SSB · PM (Γ>J): trivial paramagnet | Order parameter ⟨σ^z⟩ onsets below the critical field. |
| B8 frustration | none on bipartite FM · geometric if AFM on triangular | Default FM is unfrustrated. |
| C9 global symmetry | Z_2 spin-flip (P = Π_i σ^x_i, parity) | The symmetry whose breaking defines the FM phase. |
| C10 spatial symmetry | translation (k), inversion/parity | Conserved momentum in PBC. |
| C11 integrability | free-fermion / quadratic (1D, exact via Jordan–Wigner) · 2D non-integrable | 1D diagonalizable in O(N)/O(N³); the textbook exactly-solvable QPT. |
| C12 sign problem | sign-free (ferromagnetic / bipartite → QMC applicable) | SSE/QMC works at scale; 1D is exact anyway. |
| D13 regime | ground state (T=0) + gap; dynamics/finite-T out of card scope | E/N and gap are canonical targets. |
| D14 filling / doping | N/A (spin model) | After Jordan–Wigner: free fermions at fixed filling. |
| D15 disorder | clean by default; random-bond/field → infinite-randomness fixed point | Disordered 1D TFIM is the canonical strong-disorder RG example. |
| D16 hermiticity | Hermitian / closed | — |
Γ < J) : Z_2-broken; order parameter ⟨σ^z⟩ ≠ 0 (magnetization).Γ > J) : trivial, field-polarized along x, ⟨σ^z⟩ = 0.Γ = J) : Ising CFT, c = 1/2, exponents ν = 1, β = 1/8, z = 1.E/N; spectral gap Δ (closes at the QCP).⟨σ^z⟩ (order parameter); longitudinal correlations ⟨σ^z_i σ^z_j⟩.c from entanglement scaling at criticality.Z_2 parity sector reduces cost (per method-property-map.md §MPS).Γ_c/J = 1 exactly (self-dual / Jordan–Wigner); at criticality c = 1/2, ν = 1, β = 1/8. Ground-state energy density at Γ = J = 1: E/N = −4/π ≈ −1.2732 (Pauli convention H = −J Σ σ^z σ^z − Γ Σ σ^x; from the free-fermion dispersion, consistent with this card's Verification note).(Γ/J)_c = 3.04438(2), 3D-Ising universality (Blöte & Deng, Phys. Rev. E 66, 066110 (2002)).Key reference: dutta_2010_quantum — comprehensive downloadable review of quantum phase transitions in transverse-field spin models (1D exact solution, scaling, higher-D, dynamics, quantum information), preferred over the Sachdev textbook for an all-details source.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain (Z=2) | 2D AKLT (e.g. honeycomb) variants exist but the chain is the canonical case. |
| A2 boundary conditions | OBC (exposes the edge spins) · PBC (clean entanglement-spectrum cut) | OBC gives the 4-fold-degenerate edge-state manifold. |
| A3 statistics & local dim | spin-1; d = 3 | The ground state is an exact bond-dimension-2 MPS. |
| A4 interaction range | short-range (nearest-neighbor, bilinear + biquadratic) | Local. |
| B5 entanglement scaling | area law (constant S); entanglement spectrum exactly 2-fold degenerate; S = 2 ln 2 (open chain, two cut ends) | Each cut contributes ln 2; the exact-degeneracy is the SPT fingerprint. |
| B6 spectral gap | gapped (Haldane gap above the unique bulk ground state) | Bulk gap ≈ 0.35 (numerical); the existence of a gap is rigorously established for AKLT. |
| B7 ground-state order | symmetry-protected topological (SPT) — valence-bond solid in the Haldane phase | Protected by SO(3)/Z_2×Z_2 / time-reversal / inversion; diagnosed by hidden string order, not a local order parameter. |
| B8 frustration | none | Unfrustrated bilinear-biquadratic chain. |
| C9 global symmetry | SU(2) (total spin) — exact at the AKLT point; U(1) (S^z_tot); Z_2×Z_2 (π-rotations) protects the SPT | The VBS is an SU(2) singlet on a closed chain. |
| C10 spatial symmetry | translation, inversion (protects the SPT), reflection | Inversion is one of the protecting symmetries. |
| C11 integrability | not Bethe-integrable, but the ground state is exactly constructed (VBS / MPS) | Exact GS, energy, and correlators; the full spectrum is not solvable. |
| C12 sign problem | sign-free (unfrustrated bipartite chain → QMC applicable; DMRG is the workhorse) | No frustration. |
| D13 regime | ground state (T=0) default | Energy, string order, gap, entanglement spectrum, edge states are the targets. |
| D14 filling / doping | N/A (spin model) | — |
| D15 disorder | clean by default | — |
| D16 hermiticity | Hermitian / closed | — |
O_string^z = −lim_{|i−j|→∞} ⟨S^z_i exp(iπ Σ_{i<k<j} S^z_k) S^z_j⟩; doubly-degenerate entanglement spectrum; on an open chain, 4-fold ground-state degeneracy from emergent free spin-1/2 edge modes (the two ends carry S=1/2 each).E/N (exactly −2/3).4/9 for the AKLT state.χ; MPS natively measures string order and the entanglement spectrum (per method-property-map.md B7-SPT row).L ≲ 14 for the exact spectrum/gap and edge-state degeneracy; the analytic VBS construction provides closed-form benchmarks (E/N, string order).E/N = −2/3 — exact (VBS construction; convention H = Σ [S_i·S_{i+1} + (1/3)(S_i·S_{i+1})²]).O_string = 4/9 — exact for the AKLT VBS state.Δ ≈ 0.35 (numerical, units of J); two-point correlations decay with correlation length ξ = 1/ln 3 ≈ 0.91 sites (exact).Key reference: affleck_1987_rigorous — the defining paper: rigorous construction of the valence-bond ground state, proof of the spectral gap, exponentially-decaying correlations, and the hidden topological order of the Haldane phase.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain (Z=2) · 2D square | 1D q=3 is the canonical parafermion-CFT critical chain. |
| A2 boundary conditions | PBC ring (criticality work) · OBC · cylinder (2D) | PBC preferred for finite-size scaling of the QPT. |
| A3 statistics & local dim | qudit; d = q (3, 4, 5, …) | Local dimension grows with q, raising ED/MPS cost. |
| A4 interaction range | short-range (nearest-neighbor) | Local. |
| B5 entanglement scaling | gapped phases: area law · 1D critical (q≤4): area+log (c=4/5 at q=3, c=1 at q=4) | Central charge set by the parafermion/Potts CFT; q=3 → c=4/5. |
| B6 spectral gap | gapped FM (h<h_c) and PM (h>h_c) · gapless at QCP (q≤4) · q≥5: first-order (no gapless point) | Order of transition changes at q=5 in 1D. |
| B7 ground-state order | FM (h<h_c): Z_q SSB, q-fold degenerate · PM (h>h_c): trivial | Z_q clock order parameter ⟨Z⟩. |
| B8 frustration | none on bipartite FM · possible on non-square (e.g. triangular Potts) | Default square/chain FM unfrustrated. |
| C9 global symmetry | Z_q clock symmetry (generator Π_i Z_i) | The symmetry whose breaking defines the FM phase. |
| C10 spatial symmetry | translation (k), inversion/reflection | Conserved momentum in PBC rings. |
| C11 integrability | non-integrable in general; q=2 free-fermion (TFIM); 1D critical point described by Z_q parafermion / Potts CFT | Self-dual critical point (Kramers–Wannier-type) located analytically. |
| C12 sign problem | sign-free (ferromagnetic / bipartite) | QMC/SSE applicable; DMRG-qudit is the default workhorse. |
| D13 regime | ground state (T=0) + order parameter scan | Criticality via finite-size scaling. |
| D14 filling / doping | N/A (spin/qudit model) | — |
| D15 disorder | clean by default | — |
| D16 hermiticity | Hermitian / closed | — |
h < h_c) : Z_q-broken, q-fold degenerate; clock order parameter ⟨Z⟩ ≠ 0.h > h_c) : trivial, field-polarized; ⟨Z⟩ = 0.q ≤ 4) : continuous transition, parafermion/Potts CFT (q=3: c=4/5, ν=5/6). For q ≥ 5 (1D) the transition is first-order.E/N; Z_q order parameter ⟨Z⟩ across an h-scan.⟨Z_i Z_j^†⟩.c / critical exponent ν from finite-size collapse (q=3: ν=5/6).Z_q quantum number (per method-property-map.md §MPS, C9).N ≲ 16 for q=3) for exact spectrum and the order of the transition; TEBD imaginary-time; sign-free QMC also applicable.q=3 Clock (= 3-state Potts) QCP: h_c = 1 (self-dual), Z_3 parafermion CFT c = 4/5, ν = 5/6 ≈ 0.833 (literature range ν ∈ [0.83, 0.85]).q=4: critical c = 1; q ≥ 5 (1D): first-order transition (no continuous critical point) — Wu RMP 1982.h = 0 → q-fold-degenerate FM, E/N = −2; h → ∞ → PM, E/N = −2h; q = 2 → recovers TFIM.Key reference: wu_1982_potts — the authoritative review of the Potts model (classical and quantum, all q, duality, order of transition, exact results, CFT connections).
competing exchanges and Z₂ topological order — where QMC meets the sign problem
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 2D square lattice (Z=4 NN + 4 NNN); also 1D zigzag chain | NNN couplings span both sublattices, breaking bipartiteness for QMC. |
| A2 boundary conditions | cylinder (2D DMRG default) · torus (ED) · OBC | Cylinder wrapping/width strongly affects the intermediate regime. |
| A3 statistics & local dim | spin-1/2; d = 2 | Default S=1/2; the contested regime is specific to S=1/2. |
| A4 interaction range | short-range (NN + NNN) | Still local, but two competing couplings. |
| B5 entanglement scaling | ordered phases: 2D area law · intermediate g≈0.5: enhanced/area-violating (candidate gapless QSL → power-law) | Entanglement growth in the window is what makes it method-limited. |
| B6 spectral gap | Néel/stripe: gapless (Goldstone) over ordered GS · intermediate: gapped Z2 vs gapless U(1) QSL — unresolved | Order-of-gap in the window is part of the open question. |
| B7 ground-state order | g≲0.4 Néel SSB · g≳0.6 stripe (collinear) SSB · g≈0.5 candidate spin liquid / VBS (contested) | Two ordered phases flank a debated nonmagnetic window. |
| B8 frustration | interaction-driven (competing J_1 vs J_2) | The canonical interaction-frustrated benchmark. |
| C9 global symmetry | SU(2) (total S), U(1) (S^z_tot) | Isotropic Heisenberg couplings → full SU(2). |
| C10 spatial symmetry | translation, C_4v point group; stripe phase breaks C_4 → C_2 | Lattice-rotation breaking is a stripe-order diagnostic. |
| C11 integrability | non-integrable | No exact solution at any g≠0. |
| C12 sign problem | severe (frustration breaks the Marshall sign rule) | QMC blocked → DMRG/PEPS/VMC + PolyOpt bounds. |
| D13 regime | ground state (T=0) | Phase identification is the goal. |
| D14 filling / doping | N/A (spin model) | — |
| D15 disorder | clean by default | — |
| D16 hermiticity | Hermitian / closed | — |
g ≲ 0.4) : staggered magnetization m_s, structure-factor peak at (π,π).g ≳ 0.6) : peak at (π,0)/(0,π), breaks C_4 lattice rotation.g ∈ [0.45,0.55] (frontier) : candidate gapless U(1) spin liquid vs gapped Z2 vs valence-bond solid — diagnose with spin-liquid; dimer order parameter / VBS structure factor and entanglement spectrum.E/N; spin structure factor S(q) (locating (π,π) vs (π,0) peaks).m_s (Néel), stripe magnetization, VBS/dimer order parameter.method-property-map.md B8/C12 gate).N ≤ 40 clusters; PolyOpt for a certified lower bound. Cross-method disagreement in the window is a known, reportable phenomenon — do not average it away.g ≲ 0.4, nonmagnetic window g ∈ ~[0.4, 0.6], stripe for g ≳ 0.6 (Morita, Kaneko, Imada, J. Phys. Soc. Jpn. 84, 024720 (2015)).g = 0: reduces to NN square Heisenberg, E/N ≈ −0.6694 (Sandvik QMC benchmark) — a built-in limit check.Key reference: morita_2014_quantum — many-variable VMC with quantum-number projection mapping the full square-lattice J_1–J_2 phase diagram; a canonical, downloadable entry into the contested intermediate-regime debate.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 2D square lattice with orthogonal diagonal-dimer bonds (Z=4 from the square + dimer bonds) | The orthogonality of neighboring dimers is the defining structural feature. |
| A2 boundary conditions | torus / PBC (ED) · cylinder (DMRG) | Torus for clean momentum sectors; cylinders for 2D DMRG. |
| A3 statistics & local dim | spin-1/2; d = 2 | — |
| A4 interaction range | short-range (NN square exchange J + dimer exchange J') | Local. |
| B5 entanglement scaling | dimer phase: area law (product of singlets, S → 0) · plaquette/Néel: 2D area law (∝ L) | The exact dimer phase is essentially unentangled across dimer-respecting cuts. |
| B6 spectral gap | gapped (dimer-singlet and plaquette phases) · gapless (Néel phase, Goldstone magnons) | A spin gap protects the dimer/plaquette states; closes in the Néel phase. |
| B7 ground-state order | exact dimer-singlet product (small J/J') · plaquette-singlet (intermediate) · Néel AFM (large J/J') | The dimer phase is a short-range-entangled product; the others are conventional. |
| B8 frustration | strong geometric frustration (orthogonal dimers) | Frustration stabilizes the dimer product and the magnetization plateaus. |
| C9 global symmetry | SU(2) (total spin), U(1) (S^z_tot) | A field breaks SU(2)→U(1), tuning the magnetization plateaus. |
| C10 spatial symmetry | translation, point group of the Shastry–Sutherland lattice | Plateau superstructures break translation. |
| C11 integrability | not integrable, but the dimer-singlet product is an exact ground state for small J/J' (Shastry–Sutherland 1981) | Exact GS in a regime; the rest of the phase diagram is numerical. |
| C12 sign problem | sign-ful (geometric frustration) → QMC blocked for the frustrated regime; the exact dimer state and DMRG/ED sidestep it | Frustration turns on the spin sign problem; sign-free QMC is restricted. |
| D13 regime | ground state (T=0) default; magnetization process in a field | E/spin, spin gap, and the m(H) plateau structure are the targets. |
| D14 filling / doping | N/A (spin model; magnetization m/m_sat is the field-tuned analog) | Plateaus appear at commensurate m/m_sat. |
| D15 disorder | clean by default | — |
| D16 hermiticity | Hermitian / closed | — |
J/J', J/J' ≲ 0.675) : product of NN singlets on the dimer bonds — an exact eigenstate; full spin gap, no magnetic order.0.675 ≲ J/J' ≲ 0.76) : resonating singlets on empty plaquettes; gapped, breaks lattice symmetry.J/J', J/J' ≳ 0.76) : staggered magnetization, gapless Goldstone modes.m/m_sat (e.g. 1/8, 1/4, 1/3, …) from crystallization of triplet excitations.E/spin; spin gap.J/J'; plaquette / dimer order parameters.m(H) and plateau values m/m_sat; triplet dispersion (nearly flat → localized triplets).method-property-map.md B8/C12 frustrated-2D row).J/J' regime; VMC/NQS for variational comparison; PolyOpt for a certified energy bound in the frustrated regime.J/J'): E/spin = −3J'/8 — exact dimer-singlet energy (each dimer singlet contributes −3J'/4 per bond = −3J'/8 per spin; convention H = J Σ_NN + J' Σ_dimer).m/m_sat = 1/8, 1/3, … (and others such as 1/4, 1/2) observed in SrCu2(BO3)2 and reproduced by ED/DMRG.J/J' ≈ 0.675 (dimer→plaquette) and ≈ 0.76 (plaquette→Néel). The intermediate plaquette phase and these two boundaries are established by later studies (e.g. Corboz & Mila, PRB 87, 115144 (2013)); the cited 1998 Miyahara–Ueda paper reports a single dimer→Néel transition near J/J' ≈ 0.7.Key reference: miyahara_1998_exact — Miyahara & Ueda, the foundational analysis identifying the Shastry–Sutherland (orthogonal-dimer) model as the description of SrCu2(BO3)2: the exact dimer ground state, the spin gap, and the critical J/J' for the dimer→Néel transition.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 2D honeycomb lattice (Z=3, two-site unit cell) | Three inequivalent bond directions — the source of the compass exchange. |
| A2 boundary conditions | torus / PBC (flux sectors, exact solution) · cylinder (DMRG) | Torus is needed for the Z2 ground-state degeneracy and flux sectors. |
| A3 statistics & local dim | spin-1/2; d = 2; solved via Majorana fermions (4 Majoranas/site, projected) | Fermionization is the engine of exact solvability. |
| A4 interaction range | short-range (nearest-neighbor, bond-dependent) | Local. |
| B5 entanglement scaling | area law + topological term; gapped phase has TEE γ = ln 𝒟 = ln 2 (total quantum dimension 𝒟 = 2) | The toric-code (gapped) phase is a Z2 topological state. |
| B6 spectral gap | gapless Z2 spin liquid (Dirac cones, isotropic regime) · gapped Z2 (anisotropic regime) · field-induced gap in the gapless phase | A magnetic field gaps the gapless phase into a chiral non-Abelian (Ising-anyon) phase. |
| B7 ground-state order | Z2 quantum spin liquid — no local order; fractionalized Majorana + static Z2 flux excitations | Gapless (B-phase, Dirac) or gapped (A-phase, Abelian toric-code anyons); field → non-Abelian Ising anyons. |
| B8 frustration | bond-dependent (compass) exchange frustration | Competing x/y/z Ising axes cannot be simultaneously satisfied → exchange frustration. |
| C9 global symmetry | Z2 gauge structure (static flux W_p per plaquette, conserved); time-reversal; lattice-encoded | Each plaquette flux W_p = ±1 is a constant of motion; ground state is flux-free (Lieb's theorem). |
| C10 spatial symmetry | translation, C_3 rotation (permutes x/y/z bonds), inversion | C_3 exchanges the three coupling channels. |
| C11 integrability | exactly solvable (free Majorana fermions in a static Z2 gauge field) | Quadratic after fermionization → exact spectrum, phase diagram, and anyon content. |
| C12 sign problem | the spin model has a sign problem; the exact Majorana solution sidesteps it entirely | Direct spin QMC is sign-blocked (frustration); the free-fermion mapping makes the model exactly tractable instead. |
| D13 regime | ground state (T=0) default; the exact solution also gives finite-T and dynamics | Flux-sector energetics, gap, anyon statistics are the targets. |
| D14 filling / doping | N/A (spin model; the Majorana sector is at its natural filling) | — |
| D15 disorder | clean by default | Bond disorder studied as a perturbation (out of scope). |
| D16 hermiticity | Hermitian / closed | — |
|J_x| ≤ |J_y|+|J_z| (and cyclic), Majorana Dirac cones; no local order parameter — diagnosed by the flux structure and gaplessness.|J_z| > |J_x|+|J_y|; Abelian anyons (e, m), TEE γ = ln 2.ν = ±1; excitations are Ising (non-Abelian) anyons.γ; spectral Chern number ν (field-induced phase); anyon braiding / fusion data.O(N³)), giving the spectrum, phase diagram, and anyon content exactly (per method-property-map.md C11 free-fermion row).|J_x| ≤ |J_y|+|J_z| and cyclic permutations; the gapped A phases lie outside, e.g. the boundary |J_z| = |J_x|+|J_y| (convention H = −Σ J_a σ^a σ^a).ν = ±1, Ising (non-Abelian) anyons; the gapped A phase has ν = 0, Abelian toric-code anyons with TEE γ = ln 2.Key reference: kitaev_2005_anyons — the foundational paper: exact Majorana solution, the full phase diagram (gapless vs gapped Z2), Lieb's flux-free ground state, the field-induced non-Abelian phase, and the Chern-number classification of anyons.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 3D pyrochlore lattice (corner-sharing tetrahedra, 4-site unit cell, Z=6) | The canonical 3D frustrated geometry; FCC array of tetrahedra. |
| A2 boundary conditions | PBC (Monte Carlo on L×L×L cubic cells) · open (real crystals) | Long-range dipolar sums use Ewald summation under PBC. |
| A3 statistics & local dim | classical Ising spin (d=2, S=±1 along local ⟨111⟩); quantum spin ice → effective spin-1/2 | Classical version is a statistical-mechanics model; quantum version is a genuine QMB problem. |
| A4 interaction range | NN effective exchange (minimal model) · long-range dipolar 1/r^3 (realistic dipolar spin ice) | Remarkably, the dipolar interaction "projects" onto essentially the same ice manifold (near self-screening). |
| B5 entanglement scaling | classical: n/a (thermal ensemble) · quantum spin ice: area law with emergent-gauge structure | The classical Coulomb phase is characterized by power-law (dipolar) *correlations*, not entanglement. |
| B6 spectral gap | classical: gapless ice manifold (macroscopic degeneracy) · monopole excitations cost a finite energy Δ above the ice rules · quantum spin ice: gapless emergent photon | The 2-in-2-out manifold is the degenerate ground space; flipping to 3-in-1-out creates a monopole pair. |
| B7 ground-state order | classical spin liquid / Coulomb phase with an emergent U(1) gauge field (dipolar correlations, pinch points); excitations are fractionalized emergent magnetic monopoles. Quantum spin ice → U(1) quantum spin liquid with an emergent photon | No conventional symmetry-broken order; the "order" is the emergent gauge constraint (∇·B = 0 ice rule). |
| B8 frustration | strong geometric frustration (corner-sharing tetrahedra) | The ice rule cannot pick a unique state → extensive ground-state degeneracy (the defining feature). |
| C9 global symmetry | global Ising Z_2 (classical); U(1) emergent gauge symmetry in the Coulomb phase; spin-rotation broken to the local ⟨111⟩ axes | Local axes are fixed by crystal field; the emergent U(1) is a *low-energy* gauge structure, not microscopic. |
| C10 spatial symmetry | cubic point group Fd-3m; pyrochlore translations; large degenerate manifold per the symmetry | The structure factor's pinch points sit at high-symmetry zone-boundary points. |
| C11 integrability | not integrable | Solved by Monte Carlo (classical) / numerics (quantum), not by exact ansatz; Pauling's count is an approximation. |
| C12 sign problem | classical → sign-free Monte Carlo (positive Boltzmann weights, loop/worm updates) · quantum spin ice has a sign problem (frustrated transverse terms) | Classical dipolar spin ice is a Monte Carlo workhorse; the U(1) QSL needs sign-problem-aware methods. |
| D13 regime | mostly finite-temperature thermodynamics (specific heat, residual entropy) + dynamics (monopole transport, neutron scattering); quantum spin ice → ground state + low-T | The residual-entropy plateau and pinch points are finite-T equilibrium signatures. |
| D14 filling / doping | n/a (localized moments); "doping" = nonmagnetic dilution (e.g. Y substitution) which modifies the residual entropy | Dilution studies probe the robustness of the Pauling count. |
| D15 disorder | clean (ideal crystal) by default; real materials have stuffing/dilution disorder | Dilution and the slow-equilibration of Dy₂Ti₂O₇ are active experimental subtleties. |
| D16 hermiticity | classical (Boltzmann) / Hermitian quantum | — |
T) : thermally disordered, no ice correlations.T, above any ordering) : extensive 2-in-2-out manifold; residual Pauling entropy; emergent U(1) gauge field with dipolar spin correlations producing pinch points in the (neutron) structure factor; fractionalized magnetic monopole excitations (deconfined defects of the ice rule). No local order parameter — the "order" is the divergence-free emergent field ∇·B = 0.T orderings, e.g. all-in-all-out, can pre-empt the ideal Coulomb phase.)C(T) and the integrated residual entropy S(T→0).S(q) from (polarized) neutron scattering — the pinch-point singularities are the smoking gun of the Coulomb phase.method-property-map.md §QMC/MCRG, C12 sign-free).S ≈ (1/2) ln(3/2) ≈ 0.202 k_B per spin (≈ 0.202 R per mole of spins) — Pauling's ice-rule count, consistent with the residual entropy measured in Dy₂Ti₂O₇ by Ramirez et al. (Nature 399, 333 (1999)). Convention: nearest-neighbor / dipolar spin-ice. (Caveat: thermally well-equilibrated Dy₂Ti₂O₇ shows the residual entropy can be released at very low T — Pomaranski et al., Nat. Phys. 9, 353 (2013) — so the plateau is a quasi-equilibrium feature.)S(q) at high-symmetry zone-boundary points, the real-space dipolar-correlation signature of the emergent U(1) Coulomb phase (neutron scattering on Ho₂Ti₂O₇/Dy₂Ti₂O₇) castelnovo_2011_spin.1/r potential; their density and transport set the low-T dynamics castelnovo_2011_spin.Key reference: castelnovo_2011_spin — Annual Review of Condensed Matter Physics review "Spin Ice, Fractionalization, and Topological Order"; the all-details source for the ice rule and pyrochlore geometry, the Coulomb-phase emergent gauge field and pinch points, the residual Pauling entropy, and (centrally) the fractionalized magnetic-monopole excitations and their deconfinement. (The monopole-defining paper is Castelnovo-Moessner-Sondhi, Nature 2008, arXiv:0710.5515; this review subsumes it with the full Coulomb-phase / topological context.)
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 2D square lattice, qubits on edges (Z = 4 edges per star/plaquette) | Generalizes to any 2D surface; genus controls the degeneracy. |
| A2 boundary conditions | torus (PBC×PBC) for the topological degeneracy · planar with boundaries for codes | Ground-state degeneracy is 2^{2g} on a genus-g surface. |
| A3 statistics & local dim | spin-1/2 qubit; d = 2 | The qubit-on-edge layout is the defining feature. |
| A4 interaction range | short-range (local 4-body stabilizers) | Strictly local. |
| B5 entanglement scaling | area law with a topological correction: S = α L − γ, γ = ln 2 (𝒟 = 2) | TEE γ = ln 2 is the entanglement signature of the Z2 order. |
| B6 spectral gap | gapped (constant gap above the ground space) | Excitations cost a fixed energy per violated stabilizer. |
| B7 ground-state order | intrinsic topological order (Z2 / long-range entangled) | No local order parameter; characterized by GSD, anyons, TEE — not symmetry breaking. |
| B8 frustration | none in the usual sense (commuting stabilizers) | The model is exactly frustration-free (all terms minimized simultaneously). |
| C9 global symmetry | Z2 × Z2 (1-form) symmetries; the conserved stabilizers A_v, B_p | The Wilson/'t Hooft loop operators generate the topological sectors. |
| C10 spatial symmetry | translation, square point group D_4 | Not needed for solvability; the topological structure is symmetry-independent. |
| C11 integrability | exactly solvable (commuting stabilizer Hamiltonian) | Entire spectrum known; ground states are common +1 eigenstates of all stabilizers. |
| C12 sign problem | N/A — solved exactly; commuting projectors mean no Monte Carlo is needed | The exact stabilizer structure bypasses any sampling. |
| D13 regime | ground state (T=0); finite-T destroys the order in 2D (no thermal stability) | The 2D toric code has no finite-T topological order; 4D variant does. |
| D14 filling / doping | N/A (spin/qubit model) | — |
| D15 disorder | clean by default; disorder/perturbations studied for code thresholds | Stability under perturbation underlies its use as a quantum memory. |
| D16 hermiticity | Hermitian / closed | — |
γ.e (electric charge, violated A_v), m (magnetic flux, violated B_p), and the composite fermion ε = e × m; e and m are mutual semions (braiding phase −1).4 on the torus (2^{2g} on genus g).γ = ln 2 (total quantum dimension 𝒟 = 2).A_v, B_p (no numerics required).γ; DMRG on cylinders to extract γ from the entanglement entropy as a benchmark of the topological-order toolkit.GSD = 4 on the torus (exact; 2^{2g} on genus g).γ = ln 2 (exact; 𝒟 = 2).Δ = 2 (in units of the stabilizer coupling — one e plus one m excitation, or 2 per single anyon pair depending on convention).Key reference: kitaev_1997_fault — the founding paper: introduces the toric code, the stabilizer formalism, topological degeneracy on surfaces, Abelian anyons, and fault-tolerant quantum computation by anyon braiding.
itinerant fermions with on-site repulsion — Mott, pairing, Kondo physics
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain / quasi-1D ladder / 2D square (Z=4), triangular / 3D cubic / Z→∞ (DMFT) | Geometry sets the difficulty; 2D square is the cuprate-parent target. |
| A2 boundary conditions | OBC (DMRG) · PBC (ED/QMC) · cylinder (2D DMRG) · infinite (iDMRG/DMFT) | Cylinder width caps the 2D-DMRG bond-dim budget. |
| A3 statistics & local dim | fermion; d = 4 per site (∅, ↑, ↓, ↑↓) | The four-state local space sets the ED 4^N wall and MPS per-site cost. |
| A4 interaction range | short-range: on-site U + NN hopping (extended Hubbard adds NN V, t') | Local — area-law-compatible. |
| B5 entanglement scaling | 1D: area+log near criticality (c=1 charge + c=1 spin off half-filling; Mott-gapped charge at half-filling) · 2D: area law | Spin-charge separation gives two gapless modes in the 1D metal. |
| B6 spectral gap | half-filling: Mott charge gap for any U>0 in 1D, for U≳U_c in 2D (gapless spin) · doped: gapless (metal) | 1D Mott gap opens at infinitesimal U (Lieb–Wu); spin sector stays gapless. |
| B7 ground-state order | half-filled bipartite: AF Mott insulator (Néel SSB in 2D/3D) · doped 2D: stripes / d-wave SC candidate / pseudogap | Doped 2D ground state is the central open problem; stripes vs uniform SC nearly degenerate. |
| B8 frustration | none on bipartite (chain, square, cubic) · geometric on triangular/t'≠0 · fermionic sign always present | Fermionic antisymmetry is the intrinsic frustration. |
| C9 global symmetry | U(1)_charge × SU(2)_spin (N↑, N↓ separately conserved; S^z) · half-filled bipartite: SO(4) = spin SU(2) × η-pairing SU(2) | SO(4) (Yang–Zhang) adds pseudospin/η-pairing at half-filling on bipartite lattices. |
| C10 spatial symmetry | translation (k), point group (D_4 square), inversion | Block-diagonalizes ED sectors. |
| C11 integrability | 1D: Bethe-ansatz integrable (Lieb–Wu) · 2D / 3D: non-integrable | 1D has exact spectrum and thermodynamics (TBA); higher D fully numerical. |
| C12 sign problem | half-filled bipartite: sign-free in DQMC (particle-hole) · attractive U<0: sign-free at any filling · doped repulsive: severe sign problem | Sign-free at half-filling is what makes that regime numerically exact at scale. |
| D13 regime | ground state (T=0) default; finite-T (DQMC/DMFT) and dynamics out of card scope | E/N + double occupancy are canonical GS targets. |
| D14 filling / doping | half-filling (Mott) is the symmetric reference; doping is the key axis (turns on sign problem, opens competing orders) | Mott→doped is the decisive control parameter. |
| D15 disorder | clean by default; disorder → Anderson–Hubbard / MIT (out of scope) | — |
| D16 hermiticity | Hermitian / closed | — |
m_s, spin structure-factor peak S(π,π), charge (Mott) gap Δ_c.d_{x²−y²} superconducting pairing, pseudogap — diagnose via pair-field and charge/spin structure factors.E/N; double occupancy ⟨n↑ n↓⟩.m_s, spin structure factor S(q); charge gap Δ_c.d-wave pair correlations and charge/spin stripe order parameters (doped 2D).n(k); central charge c (1D, entanglement scaling).method-property-map.md §MPS).N): sign-free DQMC/AFQMC — particle-hole symmetry makes it numerically exact at scale (§QMC, C12).U/t=4): exact ground-state energy per site E/N ≈ −0.5737 t from the Lieb–Wu integral equations (Lieb & Wu, Phys. Rev. Lett. 20, 1445 (1968); convention H = -tΣc†c + UΣn↑n↓). At U→∞ half-filling, E/N → 4 ln 2 · (−t²/U) (Heisenberg AFM, J=4t²/U).U/t=8, δ=1/8 doping (t'=0): ground state is a period-8 stripe with E/N = −0.767 ± 0.004 t — four-method agreement (DMRG, AFQMC, iPEPS, DMET), reported in the Qin review qin_2021_hubbard from the Simons-collaboration benchmark (Zheng et al., Science 358, 1155 (2017)).Key reference: qin_2021_hubbard — the authoritative multi-method computational review of the Hubbard model (DMRG, AFQMC, DMFT/DCA, tensor networks) covering half-filling and the doped 2D problem, with cross-method consensus benchmarks.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain / 2D square (Z = 4) / 3D cubic; bipartite lattices are the canonical setting | 2D square is the standard testbed for the superconducting T_c dome. |
| A2 boundary conditions | PBC (ED / QMC) · OBC / cylinder (DMRG) · infinite (iDMRG) | Bipartite PBC clusters used for sign-free DQMC. |
| A3 statistics & local dim | fermion; d = 4 per site (∅, ↑, ↓, ↑↓) | Same four-state local space as the repulsive Hubbard; the doubly-occupied state is now energetically favored. |
| A4 interaction range | short-range: on-site U < 0 + NN hopping | Local — area-law-compatible. |
| B5 entanglement scaling | 1D: area + log near criticality (gapless pairing/charge mode) · 2D: area law | 1D Luther–Emery liquid (gapped spin, gapless charge) off half-filling. |
| B6 spectral gap | spin gap (pairs are spin singlets) for any U < 0; charge sector gapless off half-filling (superconductor), gapped at half-filling (pair CDW) | The pairing gap is the order; the spin gap is a hallmark of local-pair formation in the BEC regime. |
| B7 ground-state order | s-wave superconductor (off half-filling); at half-filling on a bipartite lattice the SC and charge-density-wave orders are degenerate (pseudospin SU(2)) | Off half-filling the pseudospin field tilts toward SC; at half-filling SC and CDW are exactly degenerate. |
| B8 frustration | none on bipartite lattices; fermionic statistics present but does not produce a DQMC sign here (see C12) | The attraction pairs up/down spins symmetrically, removing the usual fermionic sign. |
| C9 global symmetry | U(1)_charge × SU(2)_spin; PLUS an exact pseudospin (η-pairing) SU(2) at half-filling on bipartite lattices | The pseudospin SU(2) rotates between pairing and charge density — its three components are the SC order parameter (two) and the density (one); it forces the SC–CDW degeneracy. |
| C10 spatial symmetry | translation (k), point group (D_4 square), inversion | Bipartite sublattice structure underlies the sign-free property. |
| C11 integrability | 1D: Bethe-ansatz integrable (Lieb–Wu with U < 0) · 2D / 3D: non-integrable | 1D has the exact spectrum / thermodynamics; higher D is fully numerical. |
| C12 sign problem | SIGN-FREE in DQMC at ANY filling (the attraction gives identical up- and down-spin determinants, \det M_↑ \det M_↓ = (\det M)^2 ≥ 0) | This is the headline computational feature — contrast the doped repulsive Hubbard, which has a severe sign problem. |
| D13 regime | ground state (T = 0) and finite temperature (the T_c of the SC transition is the central target) | Finite-T DQMC for T_c; ground-state DMRG/ED for pairing correlations. |
| D14 filling / doping | filling tunes the order: off half-filling → s-wave SC; half-filling → the SC/CDW-degenerate (pseudospin-symmetric) point | Filling is the key control axis between the SC and SC+CDW regimes (and along the T_c dome). |
| D15 disorder | clean by default; disorder studied for the SC–insulator (localization) transition of pairs | — |
| D16 hermiticity | Hermitian / closed | — |
U < 0) : pairing order parameter Δ = ⟨c_{i↑} c_{i↓}⟩; pair-field (s-wave) correlation function and superfluid density; spin gap. Diagnose by the q = 0 pair structure factor and superfluid stiffness.|U| \ll t) to tightly-bound on-site bosonic pairs (|U| \gg t), tracked by the pair size, double occupancy, and the two temperature scales T_p (pairing) and T_c (condensation).P_s(r) = ⟨Δ_i^\dagger Δ_j⟩ and s-wave pair structure factor; superfluid density ρ_s.T_c vs filling ⟨n⟩ and |U|/t (the T_c dome).⟨n_↑ n_↓⟩ (grows toward 1 in the BEC limit); momentum distribution n(k); quasiparticle weight.T_p; at half-filling the CDW structure factor S(π,π) (degenerate with SC).N, finite-T and ground state): sign-free DQMC — the attractive interaction makes the determinant a perfect square, so DQMC is numerically exact at scale at all fillings (per method-property-map.md §QMC / C12) — the decisive advantage over the doped repulsive model.\det M_↑ \det M_↓ = (\det M)^2 ≥ 0), so the average sign is exactly 1 for any ⟨n⟩ (convention H = -tΣc†c − |U|Σn↑n↓) — the contrast with the sign-ful doped repulsive Hubbard.T_c dome (DQMC): a broad maximum T_c ≈ 0.16\,t near |U|/t ≈ 5 ± 1 and band filling ⟨n⟩ ≈ 0.79 ± 0.09 — Fontenele et al. fontenele_2022_attractive.⟨n⟩ = 1 the s-wave-SC and CDW order parameters are exactly degenerate (pseudospin / η-pairing SU(2)).|U| \ll t, large overlapping pairs, T_c rising with |U|) to BEC (|U| \gg t, tightly-bound local pairs, T_c ∝ t^2/|U| falling) with the maximum in between — Fontenele et al. fontenele_2022_attractive.Key reference: fontenele_2022_attractive — Fontenele, Costa, dos Santos & Paiva, "The 2D attractive Hubbard model and the BCS-BEC crossover" (Phys. Rev. B 105, 184502, 2022): a downloadable all-details study using sign-free DQMC to map the superconducting T_c across band filling and |U|/t, the BCS–BEC crossover (pairing T_p vs degeneracy T_d scales, double occupancy, n(k), quasiparticle weight), chosen as the key reference because it directly delivers the headline benchmarks. Broader context: Micnas, Ranninger & Robaszkiewicz, "Superconductivity in narrow-band systems with local nonretarded attractive interactions", Rev. Mod. Phys. 62, 113 (1990) (doi:10.1103/RevModPhys.62.113; no arXiv preprint).
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | single-site/impurity + bath · lattice (square/cubic) · Z→∞ (DMFT/CDMFT) | Runtime scope is the impurity/DMFT-embedded problem; lattice multiorbital is DMFT territory. |
| A2 boundary conditions | impurity: none (0D + bath) · lattice: PBC/cylinder/infinite (DMFT) | Bath discretization, not boundaries, dominates the impurity problem. |
| A3 statistics & local dim | fermion; local dim 4^M per site (64 for M=3) | The 4^M local wall is the master cost axis — cost it out before any method. |
| A4 interaction range | short-range: on-site Kanamori (U, U', J_H) + NN hopping | Local interactions; area-law compatible. |
| B5 entanglement scaling | impurity: area-law bath chain · lattice: area law (2D) / DMFT local | Impurity-as-chain is MPS-friendly; 4^M inflates per-site MPS cost. |
| B6 spectral gap | metal (Hund's metal, gapless) · Mott insulator (U > U_c(J_H)) · orbital-selective Mott (some bands gapped, others metallic) | J_H suppresses the Fermi-liquid coherence scale → bad metal above it. |
| B7 ground-state order | Hund's metal (incoherent correlated metal) · Mott / orbital-selective Mott insulator · magnetic order at low T | "Spin-freezing" non-Fermi-liquid regime is the hallmark of Hund's-metal physics. |
| B8 frustration | fermionic sign always; orbital + spin degeneracy enlarges the low-energy manifold | Multiorbital low-energy degeneracy is the source of Hund's-metal correlations. |
| C9 global symmetry | U(1)_charge × SU(2)_spin × orbital symmetry; J_H breaks full orbital rotation (keeps SO(3) only with full rotationally-invariant Kanamori) | J_H is what lowers orbital symmetry; density-density-only further breaks it. |
| C10 spatial symmetry | impurity: orbital point group (t_{2g}/e_g crystal field) · lattice: translation + point group | Crystal-field splitting labels the orbital sectors. |
| C11 integrability | non-integrable (multiorbital interactions) | No exact solution; numerical throughout. |
| C12 sign problem | generically severe in multiorbital DQMC; CT-HYB (hybridization-expansion CTQMC) sign-free for density-density, sign-ful with spin-flip/pair-hopping & off-diagonal hybridization | The severe multiorbital sign problem is why CTQMC/ED-bath dominate over lattice DQMC. |
| D13 regime | ground state + finite-T (CTQMC/DMFT); dynamics out of card scope | Hund's-metal coherence scale is a finite-T phenomenon. |
| D14 filling / doping | shell filling N (per-orbital occupancy) is the control parameter; strongest correlations away from N=M or N=1 ("Janus" fillings N=2,4 for M=3) | Average shell occupancy, not just U, sets the correlation strength. |
| D15 disorder | clean by default | — |
| D16 hermiticity | Hermitian / closed | — |
Z, large effective mass m*/m, reduced coherence temperature T_coh; non-Fermi-liquid self-energy in the spin-freezing regime (Im Σ(iω_n) ∝ ω_n^α, α ≈ 0.5 at the boundary).U > U_c(J_H); charge gap Δ_c ~ U_{eff} = U+(M-1)J_H at N=M.Z_m and spectral weight.⟨n_m⟩, double occupancy ⟨n_{m↑} n_{m↓}⟩.⟨S^2⟩, instantaneous ⟨S_z^2⟩; total spin (high-spin Hund's-rule multiplet at large J_H).Z_m = (1 - ∂Σ/∂ω)^{-1}, effective mass m*/m; coherence scale T_coh.A(ω), self-energy Σ(iω_n) (spin-freezing diagnostic).M, finite bath): ED — exact within the discretized bath, U(1)×SU(2)×orbital sectors cut the 4^M·d_bath space (per method-property-map.md §ED).M-orbital shell (N=M): the atomic/Mott gap is enhanced by Hund's coupling to U_{eff} = U + (M-1)J_H (Georges-de'Medici-Mravlje, Annu. Rev. Condens. Matter Phys. 4, 137 (2012), Eq. 11; Kanamori convention U' = U - 2J_H). For the half-filled 3-orbital Hubbard-Kanamori model, DMFT finds U_c strongly *reduced* by J_H (their Fig. 2), whereas at N=1 U_c *increases* quasi-linearly with J_H.Im Σ(iω_n) ∝ ω_n^α with α ≈ 1/2 at the frozen-moment phase boundary (Werner-Gull-Troyer-Millis, PRL 101, 166405 (2008); reviewed in georges_2012_strong).Key reference: georges_2012_strong — the authoritative review of Hund's-coupling physics (Hund's metals, spin-freezing, orbital-selective Mott, the U_c(J_H) and U_{eff}=U+(M-1)J_H relations) for multiorbital correlated metals, with DMFT benchmarks across 3d/4d TMOs and iron pnictides.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain / quasi-1D ladder / 2D square (Z=4); triangular when frustrated | 2D square doped is the cuprate-physics target. |
| A2 boundary conditions | OBC (DMRG) · PBC (ED) · cylinder (2D DMRG) | Cylinder width caps the 2D-DMRG budget; stripes need long cylinders. |
| A3 statistics & local dim | fermion; d = 3 per site (∅, ↑, ↓) — no double occupancy | The no-double-occupancy constraint is a Hilbert-space restriction, not a symmetry; d=3 (not 4) is the projected local space. |
| A4 interaction range | short-range: NN hopping + NN exchange J (extended t-J adds t', J') | Local — area-law compatible. |
| B5 entanglement scaling | 1D: area+log (gapless Luttinger liquid; c=1 charge + c=1 spin) · 2D: area law | Doped 1D t-J is a Luttinger liquid with spin-charge separation. |
| B6 spectral gap | doped: gapless · zero doping: reduces to Heisenberg (gapless 1D / 2D AFM) | At half-filling (n=1) kinetic term is fully projected → pure Heisenberg. |
| B7 ground-state order | 2D large-J/t: d_{x²-y²} superconductivity / stripes / phase separation · doped: Luttinger liquid (1D) · half-filled: AFM | Phase separation at large J/t; stripe vs uniform d-SC near degenerate in the cuprate window. |
| B8 frustration | none on bipartite lattices · geometric on triangular · fermionic sign on doping | Fermionic antisymmetry under doping is the intrinsic frustration. |
| C9 global symmetry | U(1)_charge × SU(2)_spin (N↑, N↓ conserved; S^z) | At the supersymmetric point J=2t, an enlarged su(2|1) superalgebra appears (1D). |
| C10 spatial symmetry | translation (k), point group (D_4 square), inversion | Block-diagonalizes ED sectors. |
| C11 integrability | 1D supersymmetric integrable point at J = 2t (Bethe ansatz, Sutherland/Schlottmann); otherwise non-integrable | J=2t gives an exact spectrum/thermodynamics benchmark; generic J/t is fully numerical. |
| C12 sign problem | half-filling (n=1, → Heisenberg, bipartite): sign-free · doped: severe sign problem | Doping turns on the fermion sign — QMC is blocked away from the Heisenberg limit. |
| D13 regime | ground state (T=0) default; finite-T / dynamics out of card scope | E/N + spin/charge correlations are the canonical targets. |
| D14 filling / doping | doping δ = 1-n is the central axis; n=1 (half-filling) → Heisenberg | Doping breaks particle-hole symmetry and drives SC/stripe competition. |
| D15 disorder | clean by default | — |
| D16 hermiticity | Hermitian / closed | — |
n=1) : antiferromagnet (the kinetic term is projected away → Heisenberg) — staggered magnetization, S(π,π) (2D) or power-law spin correlations (1D).K_ρ-controlled correlations, spin-charge separation; check exponents.J/t : d_{x²-y²} pairing (pair-field correlations), stripe order (charge/spin density modulation), and phase separation at large J/t (J/t ≳ 3 per Emery–Kivelson–Lin) — diagnose via density profile and pair/charge structure factors.E/N; hole density δ.⟨S_i·S_j⟩, spin structure factor S(q); charge/density correlations N(q).d-wave pair-field correlations (2D doped); single-hole dispersion / quasiparticle weight.K_ρ, central charge c (1D entanglement scaling).d=3 projected local space + U(1)×SU(2) conservation keep cost low (per method-property-map.md §MPS).J = 2t at half-filling: exactly solvable (Bethe ansatz) with ground-state energy per site E/N = -2t·ln 2 + ... reducing to the Heisenberg value at n=1 — the canonical integrable benchmark (Sutherland, Phys. Rev. B 12, 3795 (1975); Schlottmann, Phys. Rev. B 36, 5177 (1987)).J/t ≳ 3 (and at all J/t near half-filling), established by ED/numerics (Emery, Kivelson & Lin, PRL 64, 475 (1990)); the small-J/t cuprate window (J/t ≈ 0.3-0.4) hosts the contested stripe-vs-uniform-d-SC competition reviewed in dagotto_1993_correlated.Key reference: dagotto_1993_correlated — the authoritative review of computational studies of the t-J and Hubbard models for the cuprates (d-wave SC, phase separation, spin correlations vs doping, photoemission), the standard reference for what numerics established about correlated electrons in high-T_c materials.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain / quasi-1D ladder / 2D square (Z=4) / triangular (frustrated) | 1D chain is the exactly-tractable reference (→ XXZ). |
| A2 boundary conditions | OBC (DMRG) · PBC (ED/QMC) · cylinder (2D) | OBC gives Friedel oscillations; PBC needed for clean momentum sectors. |
| A3 statistics & local dim | spinless fermion; d = 2 per site (empty / occupied) | The smallest fermionic local space — cheap per-site cost. |
| A4 interaction range | short-range: NN hopping + NN repulsion V (extended adds V') | Local — area-law compatible. |
| B5 entanglement scaling | 1D: area+log (Luttinger liquid, c=1) for V < V_c; area law (gapped CDW) for V > V_c · 2D: area law | At half-filling the 1D transition (V_c = 2t, i.e. XXZ Δ=1) is BKT-type. |
| B6 spectral gap | 1D half-filling: gapless metal (|V| < 2t) / gapped CDW (V > 2t) · away from half-filling: gapless Luttinger liquid | Maps to the XXZ gap: gapless XY phase |Δ|<1, gapped Ising-AFM Δ>1. |
| B7 ground-state order | Luttinger liquid (metal) · charge-density wave at half-filling for V > V_c=2t · phase separation for V < -V_c (attractive) | CDW = staggered density ⟨n_i⟩ = n̄ ± δ(-1)^i; the XXZ Ising-AFM Néel state. |
| B8 frustration | none on bipartite (chain, square) · geometric on triangular · fermionic sign on doping/2D | Bipartite half-filling is unfrustrated and sign-free. |
| C9 global symmetry | U(1)_charge (N_f conserved); particle-hole at half-filling on bipartite lattices | No spin SU(2) (spinless); JW image has only U(1) S^z (XXZ has no SU(2) except Δ=1). |
| C10 spatial symmetry | translation (k), inversion/parity, sublattice exchange (bipartite) | CDW spontaneously breaks the Z_2 sublattice (translation by one site). |
| C11 integrability | 1D: Bethe-ansatz integrable (exact map to XXZ via Jordan–Wigner) · 2D/3D: non-integrable | 1D gives exact ground-state energy, gap, and V_c=2t from the XXZ solution. |
| C12 sign problem | 1D and 2D bipartite at half-filling: sign-free (particle-hole / bipartite) · doped / frustrated / non-bipartite: sign problem can appear | Sign-free half-filled bipartite case is QMC-exact at scale. |
| D13 regime | ground state (T=0) default; finite-T / dynamics out of card scope | E/N + charge structure factor N(q) are canonical targets. |
| D14 filling / doping | half-filling (CDW reference) is the symmetric point; doping → incommensurate Luttinger liquid | Commensurate half-filling is where the CDW can lock in. |
| D15 disorder | clean by default | — |
| D16 hermiticity | Hermitian / closed | — |
|V| < V_c = 2t at half-filling (and generically off half-filling) — power-law density correlations, no gap; characterized by K_ρ.V > V_c = 2t at half-filling — staggered density order parameter m_{CDW} = (1/N)Σ_i (-1)^i ⟨n_i⟩, peak in N(q=π), finite charge gap. Maps to the XXZ Néel state for Δ > 1.V < -2t.E/N; density profile ⟨n_i⟩.N(q) = (1/N)Σ_{ij} e^{iq(i-j)}⟨n_i n_j⟩_c; CDW order parameter m_{CDW}.Δ_c = E(N_f+1)+E(N_f-1)-2E(N_f); Luttinger parameter K_ρ, central charge c (1D).N_f conservation, small d=2 (per method-property-map.md §MPS).N): sign-free QMC — bipartite particle-hole symmetry → numerically exact at scale (§QMC, C12).Δ = V/2t). The metal→CDW (BKT) transition sits at Δ = 1, i.e. V_c = 2t; for V > 2t the CDW gap opens exponentially (XXZ Ising-AFM gap), the Yang–Yang exact result (C. N. Yang & C. P. Yang, Phys. Rev. 150, 321 (1966)). Convention H = -t Σ(c†c+h.c.) + V Σ n_i n_j.V = 0 (free fermions): exact tight-binding band, E/N = -2t/π·sin(πn) → -2t/π ≈ -0.6366 t at half-filling (n=1/2), the free-fermion limit benchmark (Voit review voit_1995_one, Luttinger-liquid V→0 endpoint).Key reference: voit_1995_one — the authoritative review of 1D Fermi (Luttinger) liquids, covering the Luttinger-liquid universality class, the exact Bethe-ansatz solution, spin-charge separation, and the CDW/Mott alternatives — the all-details source for the gapless and CDW physics of 1D spinless-fermion (t-V) and related models.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain / 2D square (Z=4) / 3D / Z→∞ (DMFT — exact) | The infinite-dimensional limit is exactly solvable (DMFT); finite-D studied by Monte Carlo over f-configs. |
| A2 boundary conditions | PBC / OBC (finite clusters) · infinite (DMFT, Bethe / hypercubic DOS) | DMFT works directly in the thermodynamic limit. |
| A3 statistics & local dim | itinerant c (spinless d_c = 2; d_c = 4 with spin) plus a classical static f-occupation {0,1} per site | The f-electrons are not a quantum degree of freedom in the dynamics — a classical Ising-like variable. |
| A4 interaction range | short-range: on-site c–f repulsion U + NN c-hopping | Local. |
| B5 entanglement scaling | the c-subsystem (for a fixed f-config) is a free-fermion area-law state; the full model = classical average over f-configurations | No genuine c–f entanglement growth — the f-config is a conserved classical field. |
| B6 spectral gap | metal at small U · c-spectral (Mott-like) gap opens at large U (band splits into lower/upper Hubbard-like sub-bands) | The metal–insulator feature is driven by U; at half-filling the CDW also gaps the spectrum at low T. |
| B7 ground-state order | half-filling, bipartite, low T: checkerboard charge-density wave (staggered f-occupation) · large U: Mott-like insulator | The simplest correlated model with both a CDW transition and a metal–insulator transition. |
| B8 frustration | none on bipartite lattices (checkerboard CDW); geometric frustration of the f-arrangement on non-bipartite lattices | The f-config ordering can be frustrated by lattice geometry. |
| C9 global symmetry | c-charge U(1) plus each n^f_i separately conserved → extensively many conserved quantities (f-electrons are static) | The macroscopic set of conserved n^f_i is the model's defining structural feature ([H, n^f_i] = 0 for all i). |
| C10 spatial symmetry | translation (k), point group (D_4 square), inversion; sublattice (bipartite) for the checkerboard CDW | The CDW spontaneously breaks the sublattice (translation) symmetry. |
| C11 integrability | exactly solvable in infinite dimensions (DMFT) — Brandt–Mielsch / Freericks–Zlatić | The Z→∞ self-energy is local and the impurity problem closes exactly; finite-D is not integrable but is sign-free. |
| C12 sign problem | sign-free: the static f-configuration reduces the quantum problem to a classical Monte Carlo over f-configs (each config = a free-fermion c-determinant, manifestly positive) | No fermion sign — the f-config sum is a classical statistical-mechanics problem. |
| D13 regime | ground state and finite-T (the CDW ordering temperature T_c is a central target); real-time / spectral via DMFT | Finite-T phase diagram (CDW T_c(U)) and the T=0 MIT are the canonical targets. |
| D14 filling / doping | half-filling (ρ_c = ρ_f = ½) → checkerboard CDW, symmetric reference; doping the c or f density changes the ordered pattern / melts the CDW | Both c-filling and f-filling are control axes. |
| D15 disorder | clean by default; an annealed/quenched random f-config maps the model onto a binary-alloy / Anderson-disorder problem | The c-electrons effectively see the f-config as a (self-consistent) disorder potential. |
| D16 hermiticity | Hermitian / closed | — |
T < T_c) : staggered f-occupation (f-electrons localize on one sublattice) — order parameter is the staggered f-density / CDW structure-factor peak S_f(π,…); also gaps the c-spectrum.U, T > T_c) : disordered f-config, metallic c-band.U) : the c-density of states splits into lower/upper sub-bands separated by a U-driven gap — metal–insulator transition with U.c-electron spectral function A(ω) / density of states (DMFT) — band splitting / MIT.f-density) and structure factor S_f(q); ordering temperature T_c(U).f-occupation pattern; c-charge gap.Z→∞; gives the c-spectral function, the metal–insulator transition, and the CDW transition analytically/numerically (per method-property-map.md A1 Z→∞ row).f-configurations — for each f-config diagonalize the free-fermion c-Hamiltonian (a positive weight); sign-free, numerically exact at scale (§C12).f-configs exactly); free-fermion c-diagonalization for any fixed f-config.T=0: ground state is the checkerboard CDW (f-electrons occupy one sublattice) for any U>0 — exact for the infinite-D / bipartite case (Brandt–Mielsch; reviewed in freericks_2003_exact).T_c(U) separates the ordered checkerboard phase from the disordered phase; T_c is non-monotonic in U (rises then falls, peaking at intermediate U) — DMFT result freericks_2003_exact.U: the c-density of states splits into lower and upper sub-bands separated by a gap ≈ U (metal–insulator transition) — exact DMFT spectral function freericks_2003_exact (convention H = −tΣc†c + U Σ n^c n^f).Key reference: freericks_2003_exact — Freericks & Zlatić, "Exact dynamical mean-field theory of the Falicov–Kimball model" (the RMP review): the all-details source for the DMFT solution, the formalism, the CDW transition, and the metal–insulator physics. (arXiv preprint cond-mat/0301188; published RMP 75, 1333 (2003).)
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain / quasi-1D ladder / 2D square (Z=4) / 3D / Z→∞ (DMFT/DMFT-DCA) | 1D KLM is the cleanest DMRG target; higher D is the heavy-fermion / quantum-criticality arena. |
| A2 boundary conditions | OBC (DMRG) · PBC (ED) · cylinder (2D DMRG) · infinite (DMFT) | Cylinder width caps the 2D-DMRG bond-dim budget. |
| A3 statistics & local dim | fermion conduction band (d_c = 4: ∅, ↑, ↓, ↑↓) plus a localized spin-½ at each site → effective d = 8 per site | The extra two-dimensional local spin space doubles the per-site cost over plain Hubbard. |
| A4 interaction range | short-range: on-site Kondo exchange J_K + NN hopping | Local — area-law-compatible. |
| B5 entanglement scaling | 1D: area law (gapped Kondo insulator at half-filling) · area+log near gapless metallic / critical regimes · 2D: area law | Heavy-Fermi-liquid metal is gapless (Fermi-surface log corrections). |
| B6 spectral gap | half-filling: spin gap and charge gap for all J_K>0 (Kondo insulator, 1D) · doped / metallic: gapless (heavy Fermi liquid) | 1D half-filled KLM is gapped for any J_K>0 (Tsunetsugu–Sigrist–Ueda). |
| B7 ground-state order | heavy Fermi liquid (paramagnetic, large Fermi surface) · RKKY antiferromagnet (small J_K) · ferromagnet (low conduction filling) · Kondo insulator (half-filling) | Doniach picture: AF order at small J_K, screened paramagnet (heavy FL) at large J_K. |
| B8 frustration | none on bipartite lattices · interaction-driven competition (Kondo vs RKKY) · fermionic sign always present | The Kondo–RKKY competition is the model's defining tension, not geometric frustration. |
| C9 global symmetry | U(1)_charge (conduction N) × SU(2)_spin (total spin of conduction + local moments); S^z conserved | Particle-hole symmetry at half-filling on bipartite lattices (sign-free DQMC). |
| C10 spatial symmetry | translation (k), point group (D_4 square), inversion | Block-diagonalizes ED sectors. |
| C11 integrability | non-integrable (the lattice model) | Contrast the single Kondo impurity, which IS Bethe-ansatz solvable; the lattice destroys integrability → full numerics required. |
| C12 sign problem | half-filled particle-hole-symmetric KLM: sign-free in DQMC · doped: severe sign problem in general | Half-filling on a bipartite lattice is the sign-free reference point. |
| D13 regime | ground state (T=0) default; finite-T (DQMC/DMFT) for the Kondo crossover T_K and T_coh | E/N, spin/charge gaps, and the Doniach phase boundary are the targets. |
| D14 filling / doping | half-filling → Kondo insulator (symmetric reference); doping turns on the sign problem and the heavy-FL / magnetic competition | Conduction filling is the key control axis (along with J_K/t). |
| D15 disorder | clean by default; disorder → Kondo-disorder / non-Fermi-liquid physics (out of scope) | — |
| D16 hermiticity | Hermitian / closed | — |
J_K ≳ J_K*); large Fermi surface enclosing conduction + local-moment count, strongly enhanced effective mass m*. Diagnostics: m* / coherence temperature T_coh, large-Fermi-surface volume.J_K, indirect RKKY exchange (∝ J_K²) orders the local moments — staggered magnetization m_s, spin structure-factor peak S(π,…).J_K>0 — measure Δ_spin, Δ_charge.E/N; spin gap Δ_spin and charge gap Δ_charge (half-filling).m_s, spin structure factor S(q) (RKKY phase).m* / quasiparticle weight Z; coherence / Kondo temperature T_K ∝ exp(−1/J_K ρ).J_K* where T_K ~ T_RKKY (RKKY-AF ↔ heavy-FL boundary).method-property-map.md §MPS); the standard tool for the 1D-KLM phase diagram.Z limit (§A1 Z→∞).J_K > 0 (no magnetic order, no Doniach transition at half-filling) — Tsunetsugu, Sigrist & Ueda, Rev. Mod. Phys. 69, 809 (1997) (convention H = −tΣc†c + J_K Σ S·s, J_K>0).T_RKKY ∝ J_K² vs Kondo scale T_K ∝ exp(−1/J_K ρ) cross at a coupling J_K*; for J_K < J_K* the ground state is an RKKY antiferromagnet, for J_K > J_K* a paramagnetic heavy Fermi liquid (Doniach 1977; reviewed in coleman_2006_heavy).Key reference: coleman_2006_heavy — Coleman's "Heavy Fermions: electrons at the edge of magnetism" review: the all-details pedagogical source covering the Kondo lattice, the Doniach phase diagram, Kondo screening vs RKKY, heavy Fermi liquids and Kondo insulators. Chosen over Tsunetsugu–Sigrist–Ueda RMP 69, 809 (1997) because it is freely downloadable (arXiv) and broader in scope, while still covering the 1D-KLM benchmarks cited below.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 0D impurity + bath (not a lattice) — bath represented as a star (from Γ(ω)) or a 1D Wilson/Lanczos chain | The bath geometry (star vs chain) is a representation choice, not a lattice. |
| A2 boundary conditions | none in the lattice sense; finite-bath ED truncates the bath, NRG uses a semi-infinite chain | Bath discretization quality is the dominant concern. |
| A3 statistics & local dim | fermion; impurity local dim d = 4 (∅, ↑, ↓, ↑↓) + bath sites (d=4 each) | Cost is impurity 4 × bath 4^{L_bath}. |
| A4 interaction range | interaction is purely local (on the impurity); hybridization is local impurity↔bath | Only the impurity is interacting — this is what makes NRG/CTQMC tractable. |
| B5 entanglement scaling | area-law along the Wilson/Lanczos bath chain (single cut) → MPS/NRG friendly | Logarithmic discretization gives energy-scale separation NRG exploits. |
| B6 spectral gap | gapless (metallic bath); below T_K a Kondo resonance pins at the Fermi level | The Kondo screening scale T_K, not a gap, is the relevant low-energy scale. |
| B7 ground-state order | unique non-degenerate Kondo singlet (local moment screened by the bath) | Local-moment regime above T_K, screened singlet below — a crossover, not SSB. |
| B8 frustration | none (single impurity); fermionic, but sign-controlled by NRG/ED | No geometric frustration; multichannel variants can give non-Fermi-liquid fixed points. |
| C9 global symmetry | total N and total S^z conserved (impurity + bath); symmetric point ε_d=-U/2 adds particle-hole symmetry; SU(2) spin with no field | PH symmetry at the symmetric point fixes ⟨n_d⟩=1 and protects the NRG iteration. |
| C10 spatial symmetry | n/a (0D); channel/orbital symmetry in multichannel/multiorbital variants | Single-channel symmetric SIAM is the canonical case. |
| C11 integrability | Bethe-ansatz solvable (Wiegmann; Tsvelick–Wiegmann) — exact thermodynamics | Provides exact T_K, susceptibility, and Wilson-ratio benchmarks. |
| C12 sign problem | NRG/ED: none · CT-HYB CTQMC: sign-free for the single-orbital SIAM | Single-orbital is benign; multiorbital with spin-flip/pair-hopping is sign-ful (→ multiorbital-hubbard). |
| D13 regime | ground state + finite-T (NRG gives full T-dependence) and dynamics (spectral functions) | NRG natively produces thermodynamics and A(ω) down to T→0. |
| D14 filling / doping | impurity occupancy ⟨n_d⟩ set by ε_d/U; ⟨n_d⟩=1 at the symmetric point (Kondo regime); mixed-valence near ε_d≈0 or ε_d≈-U | Tuning ε_d moves between Kondo, empty/full, and mixed-valence regimes. |
| D15 disorder | clean single impurity by default (disorder enters only via lattice/DMFT embedding) | — |
| D16 hermiticity | Hermitian / closed | — |
T ≫ T_K) : unscreened impurity spin — Curie-like susceptibility χ_imp ~ 1/T, impurity entropy S_imp → ln 2.T ≪ T_K) : moment quenched — Pauli-like (χ_imp finite), S_imp → 0, Kondo resonance in A(ω) at the Fermi level. Crossover, not a phase transition; characterized by the Kondo scale T_K.⟨n_d⟩ ≠ 1).⟨n_d⟩, double occupancy ⟨n_{d↑} n_{d↓}⟩; local moment ⟨S_z^2⟩.T_K; impurity contributions to entropy S_imp(T), specific heat, magnetic susceptibility χ_imp(T).R_W = 4π²χ_imp/(3γ_imp); impurity spectral function A(ω) (Kondo resonance, Hubbard satellites).T_K and the full crossover; the method of record for quantum-impurity problems (per method-property-map.md §ED/RG reasoning).L_bath ≲ 8) as an exact oracle (§ED); CT-HYB CTQMC for a continuous bath at finite-T (§QMC, sign-free single-orbital); DMRG/MPS impurity solver for a long bath chain (§MPS).T impurity entropy S_imp(T→∞) = ln 2 and Wilson ratio R_W = 2 are reproduced by NRG with high precision (Wilson, Rev. Mod. Phys. 47, 773 (1975); Bulla–Costi–Pruschke review bulla_2007_numerical, Fig. 5). NRG with discretization Λ ≈ 2 already gives static properties to within a few percent.T_K ∝ √{UΓ}·exp(-πU/8Γ + πΓ/2U) (Haldane, Phys. Rev. Lett. 40, 416 (1978)); the leading exponential T_K ~ exp(-πU/8Γ) sets the exponentially small low-energy scale that NRG is built to resolve. The Bethe-ansatz solution (Tsvelick–Wiegmann, Adv. Phys. 32, 453 (1983)) gives the exact universal thermodynamics against which NRG χ_imp(T) is checked.Key reference: bulla_2007_numerical — the authoritative review of the numerical renormalization group for quantum-impurity systems (logarithmic discretization, iterative diagonalization, fixed points, thermodynamics, spectral functions, Wilson ratio), the all-details methods source for the Anderson/Kondo impurity.
quadratic hopping with topology — edge zero modes and Chern numbers
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain, two-site (A/B sublattice) unit cell | The dimerization is the entire story. |
| A2 boundary conditions | OBC (edge states, bulk–boundary correspondence) · PBC (Bloch bands, winding number) | OBC is essential to expose the zero-energy edge modes; PBC for the bulk invariant. |
| A3 statistics & local dim | spinless fermion; d = 2 per site (single orbital) | Single-particle problem per k; the many-body state is a Slater determinant. |
| A4 interaction range | short-range: nearest-neighbor hopping only (alternating t_1, t_2) | Local; no interactions. |
| B5 entanglement scaling | area law (constant, gapped) — entanglement spectrum carries the topological signature (degeneracy in the topological phase) | Free-fermion ground state; entanglement-spectrum degeneracy ↔ edge modes. |
| B6 spectral gap | gapped (bulk gap = 2|t_1 − t_2|) everywhere except the critical point t_1 = t_2 (Dirac point, gapless) | Gap closing at t_1=t_2 is the topological phase transition. |
| B7 ground-state order | 1D symmetry-protected topological (SPT) phase (t_2>t_1) vs trivial (t_1>t_2) — class BDI / AIII (chiral/sublattice symmetry) | Topological phase: winding number ν=1, Zak phase π, protected zero-energy edge states. Trivial: ν=0, Zak phase 0. |
| B8 frustration | none (free fermions, bipartite) | — |
| C9 global symmetry | chiral / sublattice (A/B) symmetry Γ = σ_z (the protecting symmetry; quantizes the winding/Zak phase) + U(1) charge (particle number); time-reversal + particle-hole → class BDI | Chiral symmetry is what protects the SPT; breaking it (e.g. an on-site staggered potential) trivializes the topology. |
| C10 spatial symmetry | translation by one unit cell (k); inversion | Inversion also quantizes the Zak phase (inversion-symmetric SPT). |
| C11 integrability | free-fermion / quadratic → exactly solvable in O(N³) despite the lattice | Diagonalize the 2N×2N (or 2×2 Bloch) single-particle Hamiltonian; no many-body numerics needed. |
| C12 sign problem | N/A — free fermions, no Monte Carlo required | — |
| D13 regime | ground state (T=0) default; quench / dynamics also exactly tractable (quadratic) | Half-filled (one fermion per cell) lower band is the canonical ground state. |
| D14 filling / doping | half-filling (lower band filled) is the insulating, topologically nontrivial reference; the gap sits at zero energy | — |
| D15 disorder | clean by default; bond disorder preserving chiral symmetry keeps the topology (studied as a perturbation) | — |
| D16 hermiticity | Hermitian / closed by default | Non-Hermitian SSH is a heavily studied extension (asymmetric hopping → non-Hermitian skin effect, breakdown of conventional bulk–boundary correspondence). |
t_2 > t_1) : winding number ν = 1, Zak phase π; under OBC, two exponentially-localized zero-energy edge modes (one per end) — diagnose by the bulk winding number, the Zak/Berry phase, or the OBC edge-state count.t_1 > t_2) : winding number ν = 0, Zak phase 0; no edge states.t_1 = t_2) : bulk gap closes (Dirac point), topological transition.= 2|t_1 − t_2|.ν (chiral invariant) / Zak (Berry) phase (0 or π).O(N³) (Bloch bands for the invariant under PBC; real-space 2N×2N matrix for the OBC edge modes), per method-property-map.md C11 free-fermion row.N as a sanity check; DMRG/MPS to confirm entanglement-spectrum degeneracy and to extend to the interacting / extended-SSH case (where the free-fermion solution no longer applies).Δ = 2|t_1 − t_2| (PBC, convention H = Σ t_1 c†_A c_B + t_2 c†_B c_{A,n+1} + h.c.); gap closes at t_1 = t_2.ν = 1, Zak = π for t_2 > t_1 (topological); ν = 0, Zak = 0 for t_1 > t_2 (trivial) — Asbóth et al. asboth_2015_short.t_2 > t_1): exactly two zero-energy edge modes (one localized at each end), exponentially decaying with length ∝ 1/ln(t_2/t_1); none in the trivial phase.Key reference: asboth_2015_short — Asbóth, Oroszlány & Pályi, "A Short Course on Topological Insulators": the canonical pedagogical all-details source; its opening SSH chapter develops the dimerized chain, chiral symmetry, the winding number, the Zak phase, bulk–boundary correspondence, and the edge states in full.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 2D honeycomb lattice, two-site (A/B) unit cell, coordination Z = 3 | Two Dirac points (K, K') in the Brillouin zone are the seat of the topology. |
| A2 boundary conditions | PBC / torus (Bloch bands, Chern number) · ribbon / cylinder (chiral edge modes, bulk–boundary correspondence) | The ribbon geometry exposes the chiral edge states; the torus defines C. |
| A3 statistics & local dim | spinless fermion; d = 2 per site; quadratic (single-particle) | d = 2 (single orbital per site); the many-body state is a Slater determinant filling the lower band. |
| A4 interaction range | short-range: real NN hopping + complex NNN hopping + on-site mass | Local; no interactions in the canonical model. |
| B5 entanglement scaling | 2D area law (gapped); the half-filled lower-band ground state is short-range entangled but carries a nonzero Chern number | Free-fermion ground state; the Chern number is the bulk topological signature (no local order parameter). |
| B6 spectral gap | gapped (band gap at K, K') except on the phase boundary M = ±3\sqrt{3}\, t_2 \sin\phi where one Dirac point closes | Gap closing at a single Dirac cone = topological transition (C jumps by ±1). |
| B7 ground-state order | Chern insulator / quantum anomalous Hall (topological band, class A, C = ±1) for |M| < |3\sqrt{3}\, t_2 \sin\phi|; trivial band insulator otherwise | Topological phase: one chiral edge mode per edge, quantized σ_xy = ±e²/h. No symmetry breaking — the order is purely band-topological. |
| B8 frustration | none (free fermions) | — |
| C9 global symmetry | charge U(1) (particle number); time-reversal is broken by the complex t_2 (this is what allows C ≠ 0) | Broken time-reversal → class A; the staggered flux carries no net flux but breaks T. The mass term M instead breaks inversion. |
| C10 spatial symmetry | translation (k); C_3 rotation; inversion broken by M, T broken by φ | At M = 0, \sin\phi \ne 0 inversion is preserved but T is broken (topological); at \phi = 0, M \ne 0 T is preserved but inversion broken (trivial). |
| C11 integrability | free-fermion / quadratic → exactly solvable (diagonalize the 2×2 Bloch Hamiltonian per k, O(N³) real-space) | Berry curvature integrated over the BZ gives the integer Chern number. |
| C12 sign problem | N/A — free fermions, no Monte Carlo required (interacting/Haldane–Hubbard extensions are sign-ful) | — |
| D13 regime | ground state (T = 0) default; quench / Hall-response dynamics also exactly tractable (quadratic) | The Chern number and edge spectrum are the targets. |
| D14 filling / doping | half-filling (lower band filled, gap at zero energy) is the topological insulating reference | At half-filling the Fermi level sits in the gap; the quantized Hall response requires the gap to be at the chemical potential. |
| D15 disorder | clean by default; weak disorder preserves the quantized σ_xy (topological protection) until it closes the mobility gap | Topological robustness against disorder is a defining feature of the QAH plateau. |
| D16 hermiticity | Hermitian / closed | — |
|M| < |3\sqrt{3}\, t_2 \sin\phi|) : Chern number C = ±1 (sign set by \mathrm{sgn}(\phi)); one chiral edge mode per edge; quantized Hall conductance σ_xy = C e²/h. Diagnose by the BZ-integrated Berry curvature (Chern number), the ribbon edge spectrum, or the Hall conductance.|M| > |3\sqrt{3}\, t_2 \sin\phi|) : C = 0; no chiral edge modes; σ_xy = 0.M = ±3\sqrt{3}\, t_2 \sin\phi) : one Dirac point closes, C jumps by ±1.K, K'.C (BZ-integrated Berry curvature of the filled band).σ_xy = C e²/h (TKNN / Kubo).2×2 Bloch Hamiltonian on a k-grid for the Chern number / bands, and the real-space ribbon for edge modes, O(N³), per method-property-map.md C11 free-fermion row.C = ±1 (QAH) for |M| < |3\sqrt{3}\, t_2 \sin\phi|, C = 0 (trivial) otherwise (convention H = t_1 Σ_⟨ij⟩ c†c + t_2 Σ_⟨⟨ij⟩⟩ e^{iφ} c†c + M Σ ξ_i n_i) — Haldane 1988; reviewed in Hasan & Kane hasan_2010_topological.σ_xy = C e²/h = ±e²/h in the topological phase, 0 in the trivial phase (TKNN integer).|C| = 1 chiral edge mode per edge in the topological phase (none in the trivial phase) — the central self-consistency check.Key reference: hasan_2010_topological — Hasan & Kane, "Colloquium: Topological insulators" (RMP 82, 3045, 2010): the authoritative all-details downloadable review. Section II.B.2 ("Graphene, Dirac electrons, Haldane model") develops the honeycomb Dirac structure, the complex-NNN staggered-flux term, the time-reversal-breaking gap, the Chern number, the chiral edge states, and the quantized Hall response in full — chosen over the (no-arXiv) original because it is downloadable and covers the model end-to-end. The defining paper is Haldane, Phys. Rev. Lett. 61, 2015 (1988) (doi:10.1103/PhysRevLett.61.2015; predates arXiv).
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 2D square lattice; magnetic unit cell of q sites at flux p/q (Z = 4) | The enlarged magnetic cell (→ q subbands) is the defining structure; reduces to the 1D Harper equation in Landau gauge. |
| A2 boundary conditions | torus / PBC (magnetic Brillouin zone, Chern/TKNN labels) · cylinder / strip (chiral edge modes in the gaps) | The magnetic BZ requires magnetic translation operators; a strip exposes the IQHE edge modes. |
| A3 statistics & local dim | spinless fermion; d = 2 per site; quadratic (single-particle) | Single orbital per site; the spectrum is a single-particle band problem (Slater-determinant many-body state). |
| A4 interaction range | short-range: nearest-neighbor hopping with Peierls phases | Local; no interactions (Hofstadter–Hubbard is the interacting extension). |
| B5 entanglement scaling | 2D area law (gapped, when the Fermi level sits in a gap); filled-subband ground states carry nonzero Chern numbers | Free-fermion ground state; topology lives in the gap labels, not a local order parameter. |
| B6 spectral gap | self-similar fractal of gaps and q subbands (the butterfly); gapped whenever E_F sits in a butterfly gap, gap closings on subband touchings | The fractal gap structure as a function of α is the signature; each gap is robust where it is open. |
| B7 ground-state order | topological bands labeled by Chern numbers (integer quantum Hall on a lattice, class A) | Each filled set of subbands has a total Chern number C; the Hall conductance is σ_xy = C e²/h. No symmetry breaking — purely band-topological. |
| B8 frustration | none (free fermions) | — |
| C9 global symmetry | charge U(1); magnetic translation symmetry (a projective / ray representation — the two magnetic translations commute only up to a phase) | Ordinary translation is broken by the gauge field; the magnetic translation group is the residual symmetry and forces the q-site magnetic cell. Time-reversal is broken by the field. |
| C10 spatial symmetry | magnetic translations (Z_q-enlarged unit cell → magnetic BZ); point-group symmetry reduced by the field | Magnetic translations replace ordinary translations; the magnetic BZ is 1/q of the original. |
| C11 integrability | free-fermion / quadratic → exactly solvable via the Harper equation (a q×q Bloch matrix per magnetic crystal momentum), O(N³) | Diagonalizing the q×q Harper matrix over the magnetic BZ gives the full butterfly and the Chern labels. |
| C12 sign problem | N/A — free fermions, no Monte Carlo required (Hofstadter–Hubbard / interacting extensions become sign-ful, magnetic-flux sign problem) | The magnetic field is itself a generic source of the QMC sign problem in interacting variants. |
| D13 regime | ground state (T = 0) default; spectrum/transport are the targets; dynamics exactly tractable (quadratic) | The butterfly and the gap Chern labels are computed once and for all. |
| D14 filling / doping | gaps open at fillings set by the Diophantine equation r = q\, s + p\, C (TKNN; C = gap Chern number, s integer); the Hall plateau sits at each such commensurate filling | Filling selects which butterfly gap (and hence which Chern number C) determines σ_xy. |
| D15 disorder | clean by default; disorder broadens subbands into Landau-like bands with localized states between the (topologically protected) extended levels — the mechanism of IQHE plateaus | Disorder is what makes the experimental quantum Hall plateaus flat; the Chern number is robust. |
| D16 hermiticity | Hermitian / closed | — |
C of the filled subbands; quantized Hall conductance σ_xy = C e²/h; chiral edge modes in a strip geometry. Diagnose by the TKNN integers from the Diophantine equation and the strip edge spectrum.α ∈ [0,1], a self-similar fractal; at α = p/q the band splits into exactly q subbands.α — the Hofstadter butterfly.= q at α = p/q) and the butterfly gap structure.C of each gap (TKNN / Diophantine equation); Hall conductance σ_xy = C e²/h.q×q magnetic-Bloch matrix over the magnetic BZ gives the spectrum, subbands, and Chern labels in O(q³) per k; a strip diagonalization gives the edge modes. Free-fermion / quadratic, per method-property-map.md C11 row.α = p/q (coprime) the spectrum splits into exactly q subbands separated by q − 1 gaps (Hofstadter 1976) hofstadter_1976_energy.α ∈ [0,1] is a fractal with the same structure recurring at all scales (Hofstadter 1976).C solving the Diophantine equation r = q\, s + C\, p with |C| ≤ q/2 (TKNN 1982); the Hall conductance at that filling is σ_xy = C e²/h.Key reference: hofstadter_1976_energy — Hofstadter, "Energy levels and wave functions of Bloch electrons in rational and irrational magnetic fields", Phys. Rev. B 14, 2239 (1976): the defining paper that derives the Harper equation, the q-subband structure at rational flux, and the self-similar butterfly spectrum. The TKNN topological labeling (Thouless, Kohmoto, Nightingale, den Nijs, Phys. Rev. Lett. 49, 405 (1982)) supplies the Chern-number / Hall-conductance interpretation. _bib stub — original PRB predates arXiv, no PDF reachable._
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain | The defining model of a 1D topological superconductor. |
| A2 boundary conditions | OBC (Majorana end modes, ground-state parity degeneracy) · PBC (BdG Bloch bands, Z2 invariant) | OBC is essential to expose the Majorana zero modes; PBC for the bulk invariant. |
| A3 statistics & local dim | spinless fermion; d = 2 per site; BdG-quadratic | Solved by Bogoliubov–de Gennes diagonalization (pairing → particle-hole-doubled single-particle problem). |
| A4 interaction range | short-range: NN hopping + NN p-wave pairing | Local; no interactions. |
| B5 entanglement scaling | area law (constant, gapped); ground-state fermion-parity (Majorana) degeneracy under OBC | Free-fermion ground state; the topological degeneracy is the entanglement/edge signature. |
| B6 spectral gap | gapped (bulk BdG gap) except at the transition μ = ±2t where the gap closes | Gap closing at |μ| = 2t is the topological phase transition. |
| B7 ground-state order | 1D topological superconductor (class D, Z2 invariant) for |μ| < 2t vs trivial for |μ| > 2t | Topological phase: Majorana zero modes localized at the two ends; near-degenerate even/odd-parity ground states. |
| B8 frustration | none (free fermions) | — |
| C9 global symmetry | fermion parity Z2 (P = ∏(1−2n_i)) — pairing breaks charge U(1); particle-hole (BdG) symmetry; time-reversal → class D | No charge conservation: superconducting pairing breaks U(1) down to Z2 parity. The two Majorana-degenerate ground states differ by total fermion parity. |
| C10 spatial symmetry | translation (k, PBC); inversion | Bulk Z2 invariant defined from the BdG Bloch Hamiltonian. |
| C11 integrability | free-fermion / quadratic → exactly solvable (BdG / Bogoliubov diagonalization, O(N³)) | Quadratic in the fermions; exact spectrum, edge modes, and phase diagram. |
| C12 sign problem | N/A — free fermions, no Monte Carlo required | — |
| D13 regime | ground state (T=0) default; quench / braiding dynamics also exactly tractable (quadratic) | Parity-resolved ground states and the Majorana spectrum are the targets. |
| D14 filling / doping | controlled by μ (no fixed filling — pairing does not conserve particle number); μ = 0 is the particle-hole-symmetric sweet spot | μ is the tuning parameter across the topological transition, not a filling constraint. |
| D15 disorder | clean by default; on-site disorder (preserving class D) shifts but does not immediately destroy the topological phase | — |
| D16 hermiticity | Hermitian / closed | — |
|μ| < 2t, Δ ≠ 0) : Z2 invariant nontrivial; under OBC, two Majorana zero modes (one at each end) combine into a single delocalized fermion → near-degenerate even/odd fermion-parity ground states. Diagnose by the bulk Z2 invariant, the OBC zero-mode spectrum, and the ground-state parity (near-)degeneracy.|μ| > 2t) : trivial Z2 invariant; no Majorana end modes, unique ground state.μ = ±2t) : bulk gap closes.N).μ=0, t=Δ).O(N³) (Bogoliubov transformation; Bloch BdG bands for the invariant, real-space Nambu matrix for the Majorana end modes), per method-property-map.md C11 free-fermion row.N (resolve the parity sectors and degeneracy); DMRG/MPS for the interacting Kitaev chain (NN repulsion breaks the free-fermion solution) and to confirm the topological degeneracy via entanglement. The Jordan–Wigner map to the transverse-field Ising chain (transverse-field-ising) gives an independent analytic benchmark.|μ| < 2t (with Δ ≠ 0); trivial for |μ| > 2t; bulk gap closes at μ = ±2t (convention H = Σ −t(c†c+h.c.) − μ(n−½) + Δ(cc+h.c.)) — Kitaev kitaev_2000_unpaired.∝ e^{−N/ξ} (exponentially small in chain length).μ = 0, t = Δ: the end Majoranas are perfectly localized on the two terminal sites (zero localization length), exact zero-energy modes for any N.|μ|<2t) phase corresponds to the Ising ordered (ferromagnetic) phase, the transition at |μ|=2t to the Ising critical point.Key reference: kitaev_2000_unpaired — Kitaev, "Unpaired Majorana fermions in quantum wires": the defining paper introducing the model, the topological vs trivial phases, the Majorana end modes, the bulk Z2 invariant, and the Jordan–Wigner connection to the Ising chain.
random potentials, many-body localization, SYK, and driven-dissipative steady states
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | dimension is THE master axis: 1D / 2D / 3D (hyper)cubic, Z=2d | 1D & 2D (orthogonal): all states localized for any W>0 (lower critical dimension d_c=2); 3D: extended states + mobility edge → metal–insulator transition. Nothing about this model matters more than d. |
| A2 boundary conditions | OBC · PBC · quasi-1D bars / cylinders (transfer-matrix MacKinnon–Kramer) | Transfer-matrix finite-size scaling on long quasi-1D bars of cross-section M^{d-1} is the workhorse for the 3D transition. |
| A3 statistics & local dim | a single non-interacting fermion (d=2 per site formally, but no many-body space) | The "Hilbert space" is just the L^d lattice sites; cost is polynomial (O(L^{3d}) dense, far less with sparse/transfer methods), not exponential. |
| A4 interaction range | short-range (nearest-neighbor hopping) · long-range hopping 1/r^a is a studied variant (power-law random banded matrices) | NN hopping with on-site disorder is the standard Anderson model; long-range hopping can change the critical dimension. |
| B5 entanglement scaling | n/a in the many-body sense (single particle). Eigenfunctions: exponentially localized (localized phase) with localization length ξ; extended/critical at the mobility edge | The relevant "size" is the localization length ξ and the multifractal wavefunction structure at criticality, not entanglement entropy. |
| B6 spectral gap | gapless single-particle band; the key feature is the mobility edge E_c separating localized (band tails) from extended (band center) states (3D) | E_c moves with W; when E_c sweeps through the Fermi energy the system undergoes the Anderson metal–insulator transition. |
| B7 ground-state order | no symmetry-breaking order — the transition is between an Anderson insulator (localized, exponentially small DC conductivity) and a diffusive metal (extended) | Order parameter ≈ the typical (geometric-mean) local density of states / inverse participation ratio, not a Landau order parameter. |
| B8 frustration | none (single particle, bipartite hypercubic lattice) | Disorder, not frustration, is the only complication. |
| C9 global symmetry | U(1) (single-particle number) — but the decisive structure is the symmetry class: orthogonal (time-reversal, no SOC), unitary (broken T, e.g. magnetic field), symplectic (T + spin–orbit) | The Wigner-Dyson class controls the critical behavior and whether 2D can delocalize (unitary/symplectic classes can; orthogonal cannot). |
| C10 spatial symmetry | none per realization (disorder breaks translation); statistical translation/isotropy after disorder averaging | Each sample is inhomogeneous; observables are averaged or, better, characterized by their full distribution (broad at criticality). |
| C11 integrability | non-interacting → directly diagonalizable (O(L^{3d}) dense; sparse/transfer-matrix/kernel-polynomial much cheaper) | There is no many-body correlation to capture — the entire physics is the localization of single-particle eigenfunctions, accessible from the one-body matrix. |
| C12 sign problem | n/a (single particle, no Monte Carlo over a many-body amplitude) | Studied by exact one-body diagonalization, transfer matrices, kernel-polynomial DOS, and the self-consistent / supersymmetric field theories. |
| D13 regime | eigenstates and single-particle spectrum (T=0 localization properties); conductance/transport via Landauer / Kubo | One asks where states localize and how the conductance scales with system size, not a finite-T many-body quantity. |
| D14 filling / doping | the Fermi energy / mobility edge position sets metal vs insulator (3D); tuning E_F across E_c is the transition | Filling matters only through which single-particle states are occupied relative to E_c. |
| D15 disorder | quenched on-site disorder is the defining axis — the localization–delocalization transition is tuned by W (and d, and symmetry class) | Requires averaging over realizations; distributions are broad and non-self-averaging at criticality (typical ≠ mean) — use geometric means / finite-size scaling. |
| D16 hermiticity | Hermitian / closed | Non-Hermitian Anderson models (e.g. Hatano–Nelson, imaginary gauge field) are a distinct, actively-studied variant with a different (skin-effect) localization transition. |
W; 3D strong W or band tails) : eigenfunctions decay as e^{-r/ξ}; vanishing DC conductivity; Poisson level statistics.E_c (3D) : separates the two within a single spectrum; at the critical point wavefunctions are multifractal (a continuous set of exponents τ_q/f(α)), level statistics are scale-invariant ("critical statistics").ξ from transfer matrices, typical DOS (geometric mean), Thouless/dimensionless conductance g, and the multifractal spectrum at the transition.ξ(W, E) from transfer-matrix Lyapunov exponents on quasi-1D bars; finite-size scaling of Λ = ξ_M/M (the MacKinnon–Kramer ratio whose M-flow direction distinguishes metal from insulator).P_2 = Σ_i |ψ_i|^4 and its disorder average / multifractal exponents τ_q at the mobility edge.g and its scaling β-function (the basis of the scaling theory of localization).M^{d-1}), extracting Λ=ξ_M/M and fitting the crossing / scaling to get W_c and ν. The classic high-precision route (not a many-body harness solver; a one-body linear-algebra calculation).O(L^{3d}) for IPR/multifractality on modest L; sparse / kernel-polynomial / shift-invert for the DOS and interior states on large L). This is "ED" applied to the L^d one-body matrix — polynomial, not d^N.E=0), orthogonal class: critical disorder W_c ≈ 16.5 t — the standard transfer-matrix value (Slevin–Ohtsuki, Phys. Rev. Lett. 82, 382 (1999); MacKinnon–Kramer) evers_2007_anderson.ν ≈ 1.57 ± 0.02 — the high-precision transfer-matrix finite-size-scaling result quoted in the review evers_2007_anderson (contrast the ε-expansion ν≈1, noted there as asymptotically inaccurate).W>0 (no transition), from the one-parameter scaling theory (Abrahams, Anderson, Licciardello & Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979); d_c=2); the unitary and symplectic classes *can* support a 2D transition evers_2007_anderson.Key reference: evers_2007_anderson — Evers & Mirlin, "Anderson transitions" (Rev. Mod. Phys. 80, 1355, 2008): the authoritative downloadable all-details review — the scaling theory of localization, the full symmetry-class classification (Wigner-Dyson, chiral, Bogoliubov–de Gennes), critical exponents, multifractality of wavefunctions at the mobility edge, and the σ-model field theory.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain (Z=2) | 1D is where MBL is best established; higher-d MBL stability is debated (avalanche instability). |
| A2 boundary conditions | OBC (DMRG-X / shift-invert ED) · PBC (clean ED level statistics) | PBC restores translation for level-statistics ensembles; OBC is convenient for entanglement cuts and DMRG-X. |
| A3 statistics & local dim | spin-½; d = 2 per site (Jordan-Wigner → interacting spinless fermions) | Same 2^N Hilbert space as the clean chain; disorder lifts the need for momentum sectors. |
| A4 interaction range | short-range (nearest-neighbor exchange + on-site random field) | Locality is essential: it is what allows quasi-local integrals of motion (l-bits) to form in the MBL phase. |
| B5 entanglement scaling | excited eigenstates: volume-law (thermal, small W) → area-law (MBL, large W) | The eigenstate-entanglement transition is the defining MBL diagnostic; area-law excited states are why MPS/DMRG-X can target them. |
| B6 spectral gap | no protecting gap — physics is at finite energy density (mid-spectrum), many-body level spacing ∼ 2^{-N} | The relevant scale is the level statistics of bulk eigenstates, not a ground-state gap. |
| B7 ground-state order | thermal phase: no order (ergodic) · MBL phase: emergent integrability; possible eigenstate (localization-protected) order / spin-glass order in eigenstates | MBL can stabilize order and "forbidden" eigenstate order at high energy density that equilibrium would wash out. |
| B8 frustration | none (unfrustrated NN chain); the complication is quenched randomness, not geometric frustration | Disorder, not frustration, drives the physics. |
| C9 global symmetry | U(1) (S^z_tot conserved; random field breaks SU(2)→U(1) and breaks translation) | S^z-resolution is the cheap ED/DMRG-X reduction; the random field removes momentum as a good quantum number. |
| C10 spatial symmetry | none (disorder breaks translation and point group); restored only after disorder averaging | Each realization is inhomogeneous; observables are averaged (median/typical) over realizations. |
| C11 integrability | clean chain: Bethe-ansatz integrable · thermal phase (W small): non-integrable, chaotic (ETH) · MBL phase: emergent integrability — extensive set of quasi-local integrals of motion (l-bits) | The l-bit picture is the central theoretical structure: H = Σ h̃_i τ^z_i + Σ J̃_{ij} τ^z_i τ^z_j + … in dressed pseudospins τ^z (exponentially localized), contrasting the clean chain's Bethe-ansatz integrability. |
| C12 sign problem | n/a — real fields, but studied by ED / shift-invert / Krylov / DMRG-X, not QMC (mid-spectrum excited states + real-time at finite energy density are off-limits to QMC) | The target (highly-excited eigenstates, long-time dynamics) is exactly the regime where Monte Carlo has no foothold. |
| D13 regime | excited eigenstates / quench dynamics at finite energy density (mid-spectrum, infinite-T ensemble) — the headline regime | Not a ground-state problem: one targets the middle of the many-body spectrum and post-quench dynamics, not E_0. |
| D14 filling / doping | fixed S^z sector (≈ half-filling in fermion language); usually the S^z_tot=0 / largest sector | Filling fixes the sector size for shift-invert ED. |
| D15 disorder | disordered (quenched random field) — the defining axis; ergodic ↔ MBL transition tuned by W | Requires averaging over many realizations (sample multiplier); the ergodic→MBL crossover is the whole point of the model. |
| D16 hermiticity | Hermitian / closed (unitary dynamics) | Coupling to a bath / dissipation destabilizes MBL — a separate open-system question. |
W) : obeys ETH; excited eigenstates have volume-law entanglement; GOE/Wigner-Dyson level statistics (level-spacing ratio ⟨r⟩ ≈ 0.53); DC transport and thermalization after a quench.W) : ETH-violating; excited eigenstates have area-law entanglement; emergent l-bits (quasi-local integrals of motion); Poisson level statistics (⟨r⟩ ≈ 0.386); no DC transport; only logarithmic-in-time entanglement growth S(t) ∝ ln t after a quench; memory of initial-state imbalance.⟨r⟩, mid-spectrum entanglement entropy and its variance, post-quench imbalance/return probability, the logarithmic entanglement-growth slope, l-bit localization length.⟨r_n⟩ = ⟨min(δ_n,δ_{n+1})/max(δ_n,δ_{n+1})⟩ for mid-spectrum eigenvalues (GOE ≈ 0.5307, Poisson ≈ 0.3863).I(t) / staggered-magnetization memory, and entanglement entropy S(t) (linear in thermal phase, ∝ ln t in MBL).L ≲ 22–24 after S^z resolution) — gives ⟨r⟩, eigenstate entanglement, and level statistics; the standard MBL workhorse (per method-property-map.md §ED, B5 volume↔area, D15).L): DMRG-X / MPS — exploits the area-law of MBL excited eigenstates to target individual highly-excited states variationally (§MPS, D15 "excited-state area law").L (exact spectrum / oracle); always average over realizations and report the distribution, not a single sample.W_c ≈ 3.5 J (box disorder h_i∈[-W,W]) — the much-cited Pal–Huse / Luitz–Laflorencie–Alet value, with a strong caveat that finite-size drifts and avalanche arguments cast doubt on whether a sharp transition survives L→∞ abanin_2018_colloquium (Luitz, Laflorencie & Alet, Phys. Rev. B 91, 081103 (2015)).⟨r⟩ ≈ 0.5307 (GOE, thermal) → ⟨r⟩ ≈ 0.3863 (Poisson, MBL) — the standard ergodicity diagnostic (Oganesyan–Huse, Phys. Rev. B 75, 155111 (2007); Atas et al. 2013) abanin_2018_colloquium.S(t) ∝ ln t in the MBL phase (vs linear in the thermal phase) — the slow-dynamics signature explained by the l-bit dephasing (Žnidarič et al. 2008; Bardarson–Pollmann–Moore 2012) abanin_2018_colloquium.Key reference: abanin_2018_colloquium — Abanin, Altman, Bloch & Serbyn, "Colloquium: Many-body localization, thermalization, and entanglement" (Rev. Mod. Phys. 91, 021001, 2019): the authoritative downloadable all-details review — ETH and its MBL breakdown, the random-field XXZ/Heisenberg chain, the l-bit (quasi-local integrals of motion) picture, area-law excited-state entanglement, logarithmic entanglement growth, level statistics, eigenstate order, and the experimental platforms.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 0+1-dimensional (quantum mechanics); no lattice, no geometry — all-to-all coupling among N Majoranas | There is no notion of position, distance, or dimension; "geometry" is the complete graph on N sites. |
| A2 boundary conditions | N/A — no lattice, hence no boundary conditions | A single quantum-mechanical degree-of-freedom cluster of N Majoranas. |
| A3 statistics & local dim | Majorana fermions (χ_i); N Majoranas → Hilbert dimension 2^{N/2} (complex-fermion variant: Dirac fermions, d = 2 per mode) | Pairs of Majoranas form a qubit, so N Majoranas span 2^{N/2} states — the ED wall. |
| A4 interaction range | all-to-all (infinite-range), 4-body random — the opposite extreme of short-range | Every quartet (i,j,k,l) interacts; this infinite-range randomness is what enables the melonic large-N solution. |
| B5 entanglement scaling | volume-law (eigenstates are highly entangled; the thermofield-double / ground state has near-maximal entanglement) | No spatial area law (no geometry); tensor-network methods do not apply. |
| B6 spectral gap | gapless / no quasiparticles; conformal (scale-invariant) IR with a continuum of excitations | The IR is governed by a conformal Green's function, not by particle-like poles — a non-Fermi liquid. |
| B7 ground-state order | non-Fermi liquid / strange metal: conformal IR, no quasiparticles, extensive T→0 entropy; holographically dual to AdS₂ / near-extremal black holes | No symmetry-breaking order parameter; characterized by the conformal exponents, the chaos rate, and the residual entropy. |
| B8 frustration | disorder-induced (the random all-to-all couplings produce a frustrated, glass-free but highly entangled landscape) | The randomness, not geometry, drives the complexity; the model is a non-glassy quantum liquid. |
| C9 global symmetry | Majorana version: fermion parity Z_2 only (no charge) · complex-fermion version: U(1) charge | Disorder average restores statistical homogeneity; individual realizations have only parity (or U(1)). |
| C10 spatial symmetry | N/A — no lattice; large-N replica / O(N) structure after disorder averaging | No translation or point group; the relevant structure is the replica-diagonal G–Σ bilocal field. |
| C11 integrability | solvable in the large-N limit (melonic dominance → Schwinger–Dyson equations, conformal IR); finite-N is non-integrable / chaotic → ED | Large N: closed G–Σ equations. Finite N: fully chaotic (random-matrix level statistics), requires exact diagonalization. |
| C12 sign problem | N/A in the large-N (analytic) solution; finite-N studied by ED (no Monte Carlo). Real-time / many-replica numerics can be sign-ful | The standard route is analytic large-N + finite-N ED, not QMC. |
| D13 regime | finite temperature (conformal IR, free energy, entropy) and real-time chaos / dynamics (OTOC, spectral form factor) are the focus; ground state for the residual entropy | The chaos, scrambling, and thermodynamics — not a ground-state energy — are the headline targets. |
| D14 filling / doping | Majorana version: half-filling fixed by particle–hole structure · complex-fermion version: charge Q is a tuning parameter (compressibility) | Filling matters only in the complex-fermion (charged) variant. |
| D15 disorder | quenched disorder is intrinsic — the random couplings J_{ijkl}; physical quantities are disorder-averaged (self-averaging at large N) | The disorder average is part of the model definition, not an optional complication; it is what produces the O(N)/replica structure. |
| D16 hermiticity | Hermitian / closed | Non-Hermitian and Lindbladian SYK are studied extensions. |
Δ = 1/q (= 1/4 for q = 4), the maximal Lyapunov exponent, and the extensive residual entropy.G(τ) = \overline{⟨T\chi_i(τ)\chi_i(0)⟩} and its conformal IR form (fermion dimension Δ = 1/q).S_0 = N s_0.λ_L (chaos / scrambling).N.N Schwinger–Dyson (melonic) equations — the disorder-averaged G–Σ self-consistent equations solved numerically (iteration on the imaginary-time grid) give the Green's function, free energy, entropy, and conformal exponents. Analytic / numerical-saddle, not a harness lattice solver.N (build the 2^{N/2} Majorana Hamiltonian for one or many disorder realizations, diagonalize, average) — gives the spectral form factor, level statistics, OTOC, and entropy, and tests the approach to the large-N predictions, per method-property-map.md ED row (volume-law, small N).λ_L = 2πk_BT/ℏ (Eq. 12.52 of the review; the bound λ_L ≤ 2πk_BT/ℏ is conjectured for all strongly-interacting systems) — chowdhury_2021_sachdev.Δ = 1/q = 1/4 for the q = 4 model (the IR Green's function G(τ) ∝ \mathrm{sgn}(τ)/|τ|^{2Δ}) — chowdhury_2021_sachdev.S_0/N = s_0 ≈ 0.2324\,k_B per fermion for the q = 4 Majorana model (Kitaev; Maldacena–Stanford) — a non-Fermi-liquid hallmark, reviewed in chowdhury_2021_sachdev.Key reference: chowdhury_2021_sachdev — Chowdhury, Georges, Parcollet & Sachdev, "Sachdev-Ye-Kitaev models and beyond: Window into non-Fermi liquids" (RMP 94, 035004, 2022): the authoritative downloadable all-details review covering the model definition, the large-N melonic / Schwinger–Dyson solution, the conformal IR, maximal chaos, the residual entropy, finite-N numerics, the holographic dual, and the connections to strange metals.
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain (boundary-driven transport) · 2D lattice (driven-dissipative arrays) | Boundary-driven XXZ chain is the canonical transport setup; 2D driven-dissipative lattices (cavity/Rydberg arrays) host steady-state phase transitions. |
| A2 boundary conditions | OBC (boundary driving — jumps act only at the ends) · PBC (bulk dissipation) · semi-infinite leads | Boundary driving requires OBC by definition; bulk dissipation can use PBC. |
| A3 statistics & local dim | spin-½; the density matrix doubles the space: d²=4 per site, total 4ᴺ (vectorized) | The d² blow-up (vs d for a closed system) is what makes vectorized-Liouvillian ED hit a wall around N≈8–10 spins. |
| A4 interaction range | short-range coherent coupling (NN) + local (single-site) jump operators | Locality of both H and L_k is what makes MPDO / dissipative-TEBD viable. |
| B5 entanglement scaling | steady state generally has operator-space entanglement (mixed-state correlations); area-law-like for gapped dissipative steady states, can grow near criticality | Tensor-network cost is set by the operator-space (MPDO) bond dimension, the open-system analog of χ. |
| B6 spectral gap | Liouvillian gap Δ_𝓛 = -max Re λ_{≠0} (the smallest nonzero |Re λ|) sets the asymptotic relaxation rate; closes at a dissipative phase transition | The relevant gap is of the non-Hermitian 𝓛, not of H; a closing Liouvillian gap (in the thermodynamic limit) signals a dissipative phase transition or slow/algebraic relaxation. |
| B7 ground-state order | steady-state order ρ_ss instead of a ground state: dissipative phases (e.g. ordered vs disordered NESS) separated by dissipative phase transitions; bistability/limit cycles possible | No variational ground state — the analog of "the phase" is the structure of the nonequilibrium steady state and any transition between such states. |
| B8 frustration | none required; competition is between coherent drive and incoherent dissipation, not geometric frustration | The driving-vs-dissipation balance plays the role frustration plays in closed models. |
| C9 global symmetry | weak vs strong symmetries of the Liouvillian — a *strong* symmetry (U commutes with H and every L_k) gives conserved quantities and multiple steady states; a *weak* symmetry (only the full 𝓛 is invariant) just block-diagonalizes 𝓛 | This weak/strong distinction is special to open systems and governs steady-state degeneracy / decoherence-free subspaces. |
| C10 spatial symmetry | translation (bulk dissipation) · only end-to-end / inversion for boundary-driven chains | Translation block-diagonalizes 𝓛 for bulk-dissipative lattices; boundary driving breaks it. |
| C11 integrability | generic driven-dissipative 𝓛: non-integrable · special cases (boundary-driven XXZ) admit exact matrix-product steady states | Prosen's matrix-product-operator ansatz gives exact NESS and current scaling for certain boundary-driven integrable chains — rare analytic benchmarks. |
| C12 sign problem | n/a in the standard routes (vectorized-𝓛 ED, MPDO/dissipative-TEBD, quantum trajectories); real-time/open Monte Carlo can have sign issues | The workhorses are ED of 𝓛, tensor networks for ρ, and trajectory averaging — not equilibrium QMC. |
| D13 regime | nonequilibrium dynamics / steady state (NESS) — relaxation under 𝓛 and the t→∞ ρ_ss; the defining regime | Not equilibrium and not a ground state: the targets are the steady state, the relaxation (Liouvillian gap), and transport. |
| D14 filling / doping | S^z (magnetization) is generally not conserved once decay/injection jumps act; a boundary chemical-potential bias drives a spin current | Drive/dissipation pump magnetization; the analog of "doping" is the imposed boundary bias Γ_L≠Γ_R. |
| D15 disorder | clean by default; disordered jump rates / fields are a studied extension (dissipation-assisted or -hindered localization) | Disorder + dissipation is its own research direction (stability of MBL to a bath). |
| D16 hermiticity | non-Hermitian / open (Lindblad) — the defining axis: complex Liouvillian spectrum, biorthogonal eigenmodes, ρ_ss = zero-mode right eigenvector, no standard variational principle | This is the whole point of the card: closed-system ground-state machinery does not apply; one works with the density matrix / superoperator or unravels into quantum trajectories. |
ρ_ss : the unique (generically) zero-eigenvalue right eigenvector of 𝓛; its observables (steady-state magnetization, correlations, current) characterize the "phase".ρ_ss (and a closing Liouvillian gap in the thermodynamic limit) as drive/dissipation are tuned — e.g. ordered ↔ disordered steady-state magnetization, bistability, or a limit cycle.j and its system-size scaling j ∼ N^{-α} — ballistic (α=0), diffusive (α=1), or anomalous/superdiffusive depending on the XXZ anisotropy Δ (KPZ superdiffusion at the isotropic point).Δ_𝓛 (relaxation rate), steady-state order parameters, two-point steady-state correlators, the full spectrum of 𝓛 in the complex plane.⟨O⟩_ss = Tr(O ρ_ss) — magnetization, density, correlations.{λ} (complex) and the Liouvillian gap Δ_𝓛 → asymptotic relaxation time τ ∼ 1/Δ_𝓛.j and its scaling with N (transport class).Tr(ρ²), operator-space (mixed-state) entanglement, and full counting statistics of the current.N): vectorized-Liouvillian ED — build 𝓛 as a 4ᴺ×4ᴺ matrix, get ρ_ss as its zero (right) eigenvector and Δ_𝓛 from the next eigenvalue; exact but capped near N≈8–10 by the d²ᴺ wall (per method-property-map.md §ED with the open-system d² doubling, D16).N): matrix-product density operators (MPDO) + dissipative TEBD/TDVP — evolve ρ (or its purification) as a tensor network in operator space to reach the steady state; the open-system analog of DMRG/TEBD (§MPS, A1 1D, A4 local jumps).ρ; trades the d² density-matrix cost for d-dimensional states × samples, and integrates naturally with t-DMRG.ρ_ss is set by Δ_𝓛 = -max_{λ≠0} Re λ, with relaxation time τ ∼ 1/Δ_𝓛; a closing gap (in N→∞) marks a dissipative phase transition (Kessler et al., Phys. Rev. A 86, 012116 (2012)).j∼N^{-α} is ballistic (Δ<1), diffusive (Δ>1), and superdiffusive / KPZ (α=2/3) at the isotropic point Δ=1 — exact matrix-product NESS and large-scale simulation (Prosen, Phys. Rev. Lett. 106, 217206 (2011); Žnidarič 2011; reviewed in the open-many-body literature, cf. daley_2014_quantum).Key reference: daley_2014_quantum — Daley, "Quantum trajectories and open many-body quantum systems" (Adv. Phys. 63, 77, 2014): the best downloadable lattice/methods-focused all-details source — the Lindblad master equation, the quantum-trajectories (Monte Carlo wavefunction) method and its physical (continuous-measurement) interpretation, and crucially its integration with time-dependent DMRG for open many-body lattices, with worked AMO examples. Chosen over Sieberer–Buchhold–Diehl (Keldysh field theory, arXiv:1512.00637) because this card's operational focus is the lattice numerical toolkit (trajectories + tensor networks), not the field-theoretic formulation.
lattice bosons and Rydberg blockade — superfluid–Mott transition and quantum scars
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain · 2D square (Z=4) · 3D cubic (Z=6); any optical-lattice geometry | Each dimension has a well-characterized SF–Mott transition. |
| A2 boundary conditions | OBC (DMRG) · PBC (ED/QMC) · cylinder (2D DMRG); trap potential in experiments | A harmonic trap produces the "wedding-cake" shells of alternating SF/Mott regions. |
| A3 statistics & local dim | soft-core boson, local dim d = n_max + 1 (truncated) | Hard-core limit d = 2 maps to XY/XXZ; soft-core needs n_max convergence (cost ∝ d·χ³ in MPS). |
| A4 interaction range | short-range: NN hopping + on-site U (extended Bose-Hubbard adds NN V, giving supersolid/density-wave) | Local — area-law compatible. |
| B5 entanglement scaling | Mott (gapped): area law · superfluid (gapless): area-law-ish in 2D/3D, area+log in 1D (c=1 Luttinger liquid) | 1D SF is a critical Luttinger liquid; gapped Mott is the cheap regime for tensor networks. |
| B6 spectral gap | Mott insulator: gapped (incompressible, charge gap ∝ U) · superfluid: gapless (sound mode, compressible) | The Mott gap and compressibility κ = ∂n/∂μ are the order diagnostics. |
| B7 ground-state order | superfluid (off-diagonal long-range order, ⟨b^\dagger_i b_j⟩ → const, condensate) ↔ Mott insulator (no ODLRO, integer filling, gapped) | The SF–Mott quantum phase transition; SF spontaneously breaks the U(1) phase. |
| B8 frustration | none on bipartite lattices; geometric frustration possible (triangular → frustrated/supersolid) | Bipartite Bose-Hubbard is unfrustrated and sign-free. |
| C9 global symmetry | U(1) particle number (N = Σ n_i conserved); the SF phase spontaneously breaks U(1) phase rotation | The conserved N blocks the Hamiltonian and is exploited by ED/DMRG/QMC. |
| C10 spatial symmetry | translation (k), point group, inversion/parity | Block-diagonalizes ED; Mott lobes respect translation, SF is uniform. |
| C11 integrability | 1D non-integrable (generic) — except the hard-core limit, which is free-fermion (Jordan-Wigner) | The hard-core point gives exact 1D benchmarks; the soft-core model needs full numerics. |
| C12 sign problem | sign-free (bosonic, positive worldline weights) on any lattice/filling | The defining advantage: worm-algorithm / SSE QMC is numerically exact at scale — the workhorse method. |
| D13 regime | ground state (T=0, SF–Mott phase diagram) default; finite-T (thermal SF, BKT in 2D) and real-time quench dynamics all standard | Cold-atom quench experiments probe SF→Mott dynamics directly. |
| D14 filling / doping | commensurate (integer) filling → Mott insulator at small t/U; incommensurate filling → superfluid | Mott lobes occur only at integer n; off-integer the system is always (compressible) SF. |
| D15 disorder | clean by default; on-site disorder + interactions → Bose glass (compressible, gapless, insulating) | The disordered Bose-Hubbard model is the canonical Bose-glass platform. |
| D16 hermiticity | Hermitian / closed | — |
⟨b^\dagger_i b_j⟩ (condensate fraction / superfluid stiffness ρ_s), gapless, compressible (κ > 0); spontaneously breaks U(1). Occurs at large t/U and at any incommensurate filling.n ∈ ℤ, gapped (charge gap Δ ∝ U), incompressible (κ = 0), no ODLRO. Occurs in lobes at small t/U and commensurate filling.(d+1)D XY universality class with dynamical exponent z=1; the generic lobe boundary (density-changing) is mean-field-like with z=2.E/N; density n and compressibility κ = ∂n/∂μ (zero in Mott, finite in SF).⟨b^\dagger_i b_j⟩ / condensate fraction; superfluid stiffness ρ_s (winding number in QMC); momentum distribution n(k) (the experimentally imaged quantity).Δ = E(N+1) + E(N-1) - 2E(N); phase-boundary (t/U)_c of each lobe.N and finite T, the standard benchmark for the SF–Mott phase diagram (per method-property-map.md §QMC, C12 sign-free).N conservation; converge n_max and χ (§MPS).N sector); mean-field / Gutzwiller for a quick phase-boundary sketch (becomes exact as Z → ∞).n = 1 SF–Mott (Mott-lobe tip): (t/U)_c ≈ 0.2974(3) — DMRG (Kühner-White-Monien, PRB 61, 12474 (2000); consistent with Läuchli-Kollath). Convention H = -t Σ(b†b+h.c.) + (U/2)Σ n(n-1).n = 1: (t/U)_c = (J/U)_c = 0.05974(3) — worm-algorithm QMC (Capogrosso-Sansone et al., PRA 77, 015602 (2008)).U_c/(z t) ≈ 5.83 at the n=1 tip (z = 6), the Gutzwiller/strong-coupling result — exact-in-Z→∞ baseline, corrected downward by fluctuations bloch_2007_many.Key reference: bloch_2007_many — Reviews of Modern Physics "Many-Body Physics with Ultracold Gases"; the all-details downloadable source for the Bose-Hubbard model in optical lattices — derives H from the continuum, lays out the SF–Mott transition and phase diagram, the Mott-lobe structure, mean-field/Gutzwiller treatment, and the experimental signatures (momentum distribution, the Greiner SF→Mott observation). (The defining theory paper is Fisher-Weichman-Grinstein-Fisher, PRB 40, 546 (1989); no arXiv preprint exists, so it would be a stub — this RMP review is preferred as the full-text source covering the same physics.)
| Axis | Value | Note |
|---|---|---|
| A1 dimension & geometry | 1D chain (Z=2); also 2D arrays (square / other tweezer geometries) | The 1D chain is the canonical scar platform; 2D PXP also shows scarring. |
| A2 boundary conditions | OBC (tweezer chain / DMRG) · PBC (ED, clean momentum sectors) | Revival fidelity and the scar tower are cleanest with translation symmetry (PBC). |
| A3 statistics & local dim | spin-1/2 / two-level atom; d = 2 per site, but a constrained Hilbert space (no two adjacent excitations) | Constraint → Fibonacci-dimensional space, dim ∼ φ^N, φ = (1+√5)/2 (golden ratio), not 2^N. |
| A4 interaction range | short-range (NN blockade + on-site drive in the PXP limit) | Physical van-der-Waals tail is 1/r^6; the PXP limit keeps only the NN hard constraint. |
| B5 entanglement scaling | volume law for generic eigenstates (ETH bulk) · sub-thermal / anomalously low entanglement for the scar-tower states | The scar subspace is the exception that makes MPS capture the revival dynamics to long times. |
| B6 spectral gap | no protecting gap (non-integrable chaotic spectrum); the scar tower sits at ≈ equal energy spacing inside the bulk | Scars are special excited states embedded in a thermal continuum, not a low-energy gap structure. |
| B7 ground-state order | disordered (paramagnetic) at small/negative Δ · Z₂ (period-2) ordered "antiferromagnetic" phase at suitable Δ > 0 | The detuned ground state breaks translation by one site (⟨Z₂⟩ order); the scar physics is a dynamical, not ground-state, phenomenon. |
| B8 frustration | none (constraint, not competing couplings) | The blockade is a kinematic constraint rather than frustrated exchange. |
| C9 global symmetry | Z_2 (spatial inversion / reflection); particle-hole-like spectral reflection of the PXP spectrum | No U(1): the drive does not conserve excitation number (Δ=0 PXP). |
| C10 spatial symmetry | translation (k), inversion/parity | Scar states carry definite momentum (k=0 and k=π); used to resolve the tower in ED. |
| C11 integrability | non-integrable (level repulsion, ETH for the bulk) — but hosts quantum many-body scars (weak ergodicity breaking, ETH violated only in the scar subspace) | The scarred subspace is approximately decoupled from the thermal bulk (an approximate su(2) "spectrum-generating" algebra). |
| C12 sign problem | n/a — this is a real-time-dynamics / ED-MPS target, not a QMC target | Quench dynamics from |Z₂⟩ is the workhorse calculation; no Monte Carlo sampling. |
| D13 regime | real-time quench dynamics (revivals from |Z₂⟩) is the defining regime; also full-spectrum (ED) and ground state (detuned phase) | |Z₂⟩ is an infinite-temperature state for the constrained ensemble, yet shows periodic revivals. |
| D14 filling / doping | n/a (spin/qubit model; detuning Δ plays the role of a chemical potential for excitations) | Tuning Δ drives the disordered ↔ Z₂-ordered transition. |
| D15 disorder | clean (translation-invariant) by default | Site/detuning disorder can be added (tweezer arrays), tuning toward localization. |
| D16 hermiticity | Hermitian / closed (unitary quench) | Dissipation (atom loss, spontaneous emission) is a separate open-system extension. |
Δ) : no broken symmetry; short-range correlations.Δ > 0) : staggered Rydberg density ⟨n_i⟩ = n̄ ± δ(-1)^i, order parameter ⟨Z₂⟩ (structure-factor peak at q=π); the ground state is close to the |•◦•◦…⟩ pattern. The disordered→Z₂ transition is in the (1+1)D Ising universality class.Δ≈0) : a tower of ≈ equally-spaced special eigenstates with anomalously large overlap on |Z₂⟩; diagnosed by long-lived revivals of the |Z₂⟩ autocorrelation/fidelity and by the sub-thermal entanglement of the tower states.|Z₂⟩: return probability / fidelity |⟨Z₂|ψ(t)⟩|², domain-wall or staggered-magnetization dynamics, entanglement-entropy growth.|⟨Z₂|E_n⟩|² vs energy (the scar tower), entanglement entropy vs energy (scars are outliers), level-spacing statistics (Wigner-Dyson for the bulk → non-integrable).⟨Z₂⟩ order parameter; Hilbert-space dimension growth ∼ φ^N.N ≲ 32, Fibonacci-reduced dimension after symmetry) give the full spectrum, the scar tower, level statistics, and exact quench dynamics; the universal oracle here (per method-property-map.md §ED, A3/C10).N): MPS / TDVP — the scar subspace's sub-thermal entanglement keeps the bond dimension manageable out to many revival periods (§MPS, D13); standard TEBD/TDVP for the constrained chain.dim ∼ φ^N with φ = (1+√5)/2 ≈ 1.618 (PBC: dim = L_N Lucas number; OBC: Fibonacci) — the kinematic signature of the infinite blockade serbyn_2020_quantum.|Z₂⟩ quench shows periodic revivals of the return probability with period T_rev ≈ 2π/(1.3–1.5 Ω) (i.e. an effective frequency slightly above the bare Ω) — the original observation in Rydberg arrays serbyn_2020_quantum.≈ ΔE_scar, the inverse revival period), with anomalously large |⟨Z₂|E_n⟩|² standing out against the thermal continuum serbyn_2020_quantum.Key reference: serbyn_2020_quantum — pedagogical Nature Physics review of quantum many-body scars and weak ergodicity breaking; works the PXP model explicitly (Eq. 1, the projector-dressed drive), the |Z₂⟩-quench revivals, the scar tower of special eigenstates with enhanced |Z₂⟩ overlap, the approximate spectrum-generating algebra, and the MPS/TDVP analysis — the best single downloadable all-details source. (The defining theory paper is Turner-Michailidis-Abanin-Serbyn-Papić, arXiv:1711.03528; this review subsumes and contextualizes it.)