The Hamiltonian
H = −t ∑⟨ij⟩,σ (c†iσcjσ + h.c.) + U ∑i ni↑ni↓
t = 1, U = 2.7825594, N↑ = N↓ = 4 on 14 sites.
- lattice
- chain
- sites
- 14
- boundary
- periodic
- Nf
- 4
- U
- 2.7825594
- energies
- 2
Stored as a total energy; quoted here as
E/site, which is what the papers report
(how the conversion works).
Every published energy
Grouped by what the number is. The exact energy comes first where there is one; in every other group the lowest unflagged energy is its best value, and each row's distance is to that best (rules §6).
Numerically exact (1)
Exact diagonalization, an exact solution, or sign-problem-free QMC where that is established. An exact energy is shown first, before the best of each other class. QMC is exact only within its error bar, which it must state, and reads exact (stochastic); an exact diagonalization outranks it. For a diagonalization resolved by symmetry sector only the ground state's sector is listed; the lowest energies of the other sectors are the instance's spectrum, in its JSON, and are not shown as the exact energy.
Strict variational upper bounds (1)
Energies that a variational state reached (VMC, DMRG and MPS at finite bond dimension, PEPS, NQS, VQE). Each is an upper bound on the ground-state energy.
How these numbers were read
Where each number below was read, and how it was checked.
-0.8300032, ED
Checked 2026-09-29 · QMBL ED cost run checks/cost/ed/results/Hubbard--chain_14_P_4_2.7825594.json (checks/cost/ed/run_ed.py: NetKet 3.22.4 operator, sparse CSR built in blocks, SciPy eigsh/ARPACK k=1 SA; commit 716e043); results file read locally, no rerun
Reported as. energy -11.620044877345844 (stored_exact -11.620044877345785, relative_difference 5.04e-15, reproduces_stored true)
|reproduced - stored| = 5.86e-14 (relative 5.04e-15); within 1e-9 relative. Sector N_up = N_dn = 4, dimension 1002001, 181 Lanczos steps (the file states no residual). The stored value is the source's own ED (varbench@2024-10-22), so two codes agree; a third code, E1's node Lanczos (qmbl-runs/qmbl-verify/out/E1-exact-recompute/results/Hubbard__chain_14_P_4_2.7825594.json), gives -11.620044877344789 (1.06e-12 from the cost run). Also read by VB2-varbench-code-hubbard#4 (scope value: VarBench Hubbard/chain_14_P_4_2.7825594.md row 'Exact diagonalization'; checks/hubbard-u-labels/ed_results.jsonl job {lattice: chain, L: 14, n: 4, U: 2.7825594}).
-0.8300032, DMRG
Checked 2026-09-29 · run script read in full: programs/dmrg_itensors_hubbard/chain_14_P_4_2.7825594.jl (varbench/methods@ed31bb0); same physics in MaxSilvester/Hubbard-and-Heisenberg-codes@b5c28f3 hubbard_chain_14_P_4_2.7825594.jl; VarBench history (git clone of github.com/varbench/varbench, partial, read with git log --follow; old names normalised); energy compared with the QMBL full-space recompute
Reported as. -11.6200448772 | | 2.0e-10 | DMRG(MaxBondDim ~1500) (VarBench upload 053b62d/755133b (Silvester, 2022-09-09; variance first printed 1.964224338735221e-10, then rounded to 2.0e-10))
H = -t sum (Cdag C + h.c.) + U Nupdn with t = 1 on a 14-site ring with the [1,N] bond (ITensors OpSum, automatic fermion signs, so fermionic PBC), N_up = N_dn = 4: the instance. The energy is the reached variational energy after the last sweep at maxdim 1550 ('~1500' in the string), above the exact energy by 1.5e-10, consistent with its variance: Temple check (E - E0 <= Var/(E1 - E), E1 = -11.520287464180326 the next level, eigsh lowest): E - E0 = 1.5e-10 vs bound 2.0e-9, consistent. bound_type variational is right; the extrapolated value in the method string is a separate row (answers the M8/M6 flags; the tiny V-score is genuine 1D convergence). No output log exists in either repository, so the number itself is not located independently (scope config).