Where each number below was read, and how it was checked.
-1.2199062, DMRG
Checked 2026-09-29 · programs/dmrg_itensors_hubbard/rectangular-4x16_64_PO_30_2.jl (varbench/methods@ed31bb0) read in full; the author's own copy MaxSilvester/Hubbard-and-Heisenberg-codes@b5c28f3 rectangular-4x16_64_PO_30_2.jl has the same Hamiltonian, filling and schedule; VarBench history (git clone of github.com/varbench/varbench, partial, read with git log --follow; old names normalised)
Reported as. -78.074 | | 0.250 | DMRG(MaxLinkDim = 4000, MaxTruncErr~1.2E-04) (VarBench upload by J. Maxwell Silvester)
4x16 cylinder: periodic around the 4-site circumference, open along the 16-site length (PO), t = 1, U n_up n_dn with no shift, N_up = N_dn = 30, conserve_qns: the instance. Maximum bond dimension in the script equals the MaxLinkDim printed. Energy is the reached DMRG energy: variational is right. History: uploaded 2022-12-15 (6fb0167), unchanged. No output log in either repository (scope config). Provenance note: in VarBench's copies of the ten rectangular scripts the variance line is evaluated on the initial product state ('E,psi = dmrg(H,ψ,sweeps)' then 'var = inner(H,ψ,H,ψ) - E^2'; ITensors v0.6.16 dmrg optimises psi = copy(psi0), src/lib/ITensorMPS/src/dmrg.jl line 189, so ψ stays the product state), so as committed they would print <ψ0|H^2|ψ0> - E^2 (for the half-filled 32_8 start, a Néel product state, 2 x 124 allowed hops give 248 - 32.923^2 = -836), not the stored positive variance. The author's repository fixed it on 2024-08-05 ('E,ψ = dmrg(...)' plus the var line, e.g. commits e78c997, 4d3e522); the stored variances were uploaded in 2022, before either script existed. Suggest pointing the code link at MaxSilvester/Hubbard-and-Heisenberg-codes or fixing varbench/methods.