{
  "model": "Hubbard",
  "lattice": "square",
  "n_sites": 64,
  "boundary": "PA",
  "params": {
    "Nf": 32,
    "U": 6
  },
  "instance_id": "Hubbard/square_64_PA_32_6",
  "rows": [
    {
      "energy": -42.676,
      "sigma": 0.007,
      "energy_variance": 0.82,
      "dof": 64,
      "einf": 96,
      "v_score": 0.0027289225735036777,
      "method": "VMC Hidden Fermion Determinant State Ansatz (N_hidden = 16. Single hidden layer fully connected net with alpha = 1). Soft mean-field constraint for Neel order.",
      "bound_type": "variational",
      "bound_type_reason": "vmc|\\brbm\\b|\\brnn\\b|jast",
      "reference": "[paper](https://www.pnas.org/doi/full/10.1073/pnas.2122059119)",
      "source": "varbench@2024-10-22",
      "provenance": "imported",
      "defect": {
        "flag": "below-exact",
        "finding": "VarBench stores -42.676 (E/N = -0.666813) and cites the HFDS paper, Moreno et al., PNAS 119, e2122059119 (2022). The paper gives -0.6574(2) for this lattice (8x8, half filling, U = 6, periodic along one side and antiperiodic along the other) in both the PNAS supplement (Table 5, read via Europe PMC PMC9371695) and the arXiv version (SI Table V). The stored number also sits 0.5 below the sign-free AFQMC total for this lattice, -42.17(2) (Qin, Shi & Zhang, PRB 94, 085103 (2016), Table IV) - 1.2%, about 25 sigma - which no variational energy can do. It breaks the paper's own size trend at U = 6 as well (-0.68135, -0.6609, -0.6574 for L = 4, 6, 8).",
        "diagnosis": "The published value is a minimum of the optimization trace, not a converged measurement. Three seeds each, netket 3.22.4, alpha=1 complex RBM, 2000 SR steps. Evaluating the trained parameters by FULL SUMMATION over all 1024 basis states (zero Monte-Carlo error) puts every converged energy ABOVE the exact value, as the variational principle requires, while every training trace dips 1.5e-3 to 3.5e-3 BELOW it. Each published value lies between the two.",
        "ruled_out": "Not an autocorrelation-underestimated error bar: tau_corr <= 0.05 and R_hat = 1.0000 across all six runs.",
        "evidence": "checks/results-tfising-rbm.json, checks/tfising_rbm_check.py"
      }
    },
    {
      "energy": -42.0736,
      "sigma": 0.0128,
      "energy_variance": null,
      "dof": 64,
      "einf": 96,
      "v_score": null,
      "method": "VMC Hidden Fermion Determinant State Ansatz (N_hidden = 16, fully parametrized hidden sub-matrix, hidden-unit density alpha = 1)",
      "bound_type": "variational",
      "bound_type_reason": "variational ansatz; energy is a strict upper bound (assigned during source verification)",
      "reference": "Robledo Moreno, Carleo, Georges & Stokes, Fermionic wave functions from neural-network constrained hidden states, Proc. Natl. Acad. Sci. U.S.A. 119, e2122059119 (2022), arXiv:2111.10420, doi:10.1073/pnas.2122059119",
      "peer_reviewed": true,
      "source": "sweep-worklist-2026-09-15",
      "provenance": "primary",
      "verified": {
        "checked_on": "2026-09-15",
        "method": "PNAS supplementary PDF (via Europe PMC, PMC9371695) extracted locally, cross-read against the arXiv version; no LLM transcription",
        "reported_as": "-0.6574(2) per site, no variance",
        "note": "PNAS SI Table 5 (identical to arXiv SI Table V), \"Variational energy per site in the L x L Hubbard model at half filling with periodic boundary conditions along one of the sides of the square and anti-periodic boundary conditions along the other side\", row L = 8, column U = 6. SI Sec. 6: N_hidden = 16, alpha = 1 for 8x8, the settings in the VarBench method string. Column alignment checked as for square_36_PA_18_2. As a total, -42.0736 sits 2.3e-3 (relative) above the sign-free AFQMC value -42.17(2) for this lattice (Qin, Shi & Zhang, PRB 94, 085103 (2016), Table IV, PBC-APBC), as a variational energy must, and fits the paper's size trend at U = 6 (-0.68135, -0.6609, -0.6574 for L = 4, 6, 8).",
        "secondary_of": null
      }
    },
    {
      "energy": -42.17,
      "sigma": 0.02,
      "energy_variance": null,
      "dof": 64,
      "einf": 96,
      "v_score": null,
      "method": "AFQMC, sign-problem-free at half filling (tau = 0.01, Trotter error below the statistical error)",
      "bound_type": "exact",
      "bound_type_reason": "sign-problem-free AFQMC at half filling on a bipartite lattice; the paper states the results are numerically exact (RULES.md 4)",
      "reference": "Qin, Shi & Zhang, Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method, Phys. Rev. B 94, 085103 (2016), arXiv:1605.09421",
      "peer_reviewed": true,
      "source": "exact-2026-09-15",
      "provenance": "primary",
      "verified": {
        "checked_on": "2026-09-15",
        "method": "arXiv PDF extracted locally with pypdf in layout mode, E rows of Table IV parsed by column position; no LLM transcription",
        "reported_as": "E = -42.17(2) (total), 8 x 8, U = 6, PBC-APBC",
        "note": "Table IV, \"Total ground state energy (E) ... in the Hubbard model at half-filling, for U = 2, 4, 6, 8 ... for both PBC and PBC-APBC. Statistical errors are on the last digit.\" Row 8 x 8, column U = 6 PBC-APBC. Half filling: N_up = N_dn = 32. Cross-checks: the 8 x 8 PBC values at U = 4 and 8 agree with VarBench's own AFQMC exact rows to 1.3 sigma, and the 8 x 8, 10 x 10 and 12 x 12 PBC U = 8 values with the AFQMC references quoted by arXiv:2507.02644 (Table S1) to the printed digits.",
        "secondary_of": null
      }
    }
  ],
  "url": "https://qmbl.org/i/Hubbard/square_64_PA_32_6/",
  "per_site_divisor": 64,
  "per_site_label": "E/site",
  "record": {
    "energy": -42.0736,
    "sigma": 0.0128,
    "method": "VMC Hidden Fermion Determinant State Ansatz (N_hidden = 16, fully parametrized hidden sub-matrix, hidden-unit density alpha = 1)",
    "reference": "Robledo Moreno, Carleo, Georges & Stokes, Fermionic wave functions from neural-network constrained hidden states, Proc. Natl. Acad. Sci. U.S.A. 119, e2122059119 (2022), arXiv:2111.10420, doi:10.1073/pnas.2122059119",
    "energy_per_site": -0.6574
  },
  "no_record_reason": null
}
