Skip to content

data(qforge): CaH6 V_NL-bound bare gate re-run — λ=20.04 (357.9% off QE, gate stays HELD)#2963

Open
dancinlife wants to merge 4 commits into
mainfrom
qforge-vnl-bound-rerun
Open

data(qforge): CaH6 V_NL-bound bare gate re-run — λ=20.04 (357.9% off QE, gate stays HELD)#2963
dancinlife wants to merge 4 commits into
mainfrom
qforge-vnl-bound-rerun

Conversation

@dancinlife

Copy link
Copy Markdown
Contributor

What

Re-runs the bare full-basis vertex gate (cah6_fullbz_xval.hexa, npw_cap=0, n=645) on the V_NL-structure-factor-fixed (bound) SCF from PR #2959 (qforge-vnl-structfac), to test whether the V_NL root fix closes the CaH6 λ accuracy gap.

Result (d6 VERBATIM — NOT tuned)

run SCF λ vs QE 4.376
OLD baseline #2768 (pre-fix framing) converged 4.13647 5.47%
V_NL-bound, nq=4 converged iters=32, e_band −65.1002 Ha 20.0374 357.9%
V_NL-bound, nq=2 identical 20.0374 357.9% (byte-identical)

The V_NL-bound SCF moves the bare CaH6 λ DRAMATICALLY AWAY from QE, not toward it. This is an OUTCOME-3 honest finding:

  1. The bare gate driver was never on the +2.744-Ha unbound state — that was the NPW=64 compose driver (V_NL OFF). The full-basis gate SCF was already bound pre-fix (e_band −65.2189 Ha); the V_NL phase fix shifts it only slightly (−65.1002 Ha).
  2. That small eigenstate shift propagates catastrophically: Σ|g|² barely moves (3.695→3.685) but assembled λ explodes 4.137→20.04. The bound eigenstates produce near-soft DFPT mode(s) whose 1/ω² weight blows up the bare (un-screened) α²F integral.
  3. Residual re-located: NOT the SCF eigenstates (now V_NL-complete, bound, converged) — it is the screened ε⁻¹ vertex that QE's |g|² carries and the bare from-scratch vertex lacks.

Regression-clean

qforge_cah6_qe_xval_test PASS — hybrid L3 α²F assembler vs QE λ_BZ rel-ε=1.6524e-07 (the 1.65e-7 anchor preserved; independent of SCF eigenstates → confirms the V_NL fix does not regress the production HYBRID path).

Verdict

Migration gate stays HELD. The HYBRID route (QE |g|² → QForge L3, rel-ε 1.65e-7) remains the only QE-grade production accuracy path.

Verdict: .verdicts/qforge-cah6-vnl-bound-rerun/ (VERDICT.md + nq4/nq2 logs + selftest log).

⚠ Depends on the V_NL fix in PR #2959 (qforge-vnl-structfac, still open). Branch cut from origin/qforge-vnl-structfac; rebase onto main after #2959 merges.

🤖 Generated with Claude Code

nbcorr-agent and others added 4 commits June 9, 2026 00:31
…root

The Kleinman-Bylander nonlocal block qforge_vnl_block builds the site-
INDEPENDENT form factor Σ_ij β̃_i(k+G) D_ij β̃_j(k+G') · [angular]. It was
added to the PW Hamiltonian with NO per-site phase, i.e. as a single atom
at the origin. For any cell with atoms off the origin (CaH6 = Ca@origin +
6 H at general sites) the H projectors were applied without their
Σ_k e^{-i(G_a-G_b)·τ_k} structure factor — the same S(ΔG) that already
phases V_loc. The KB eigenstates |ψ_k⟩ were therefore slightly wrong,
which is the V_NL-incompleteness the screened-vertex + SSCHA campaigns
both named as the CaH6 λ=4.13647 (5.47% off QE 4.376) root.

Fix: multiply the V_NL block by Re[S_s(ΔG)] = sg[idx*2] (identical to the
V_loc term) in BOTH qforge_assemble_h (single-species) and
qforge_assemble_h_multi (per-species sg_s). Single atom at origin → S=1
(prior anchors recovered exactly). DFPT/Sternheimer/screening/el-ph reuse
the assembled H closure, so the linear-response operator inherits the fix
(d4-generic, no per-species branch).

Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>
cah6_realcell_compose_xval_vnl.hexa — clones the realcell compose driver but
stages the REAL KB projectors (ca/h betas·ls·dij, nprojs=[6,2]) into the SCF
instead of nprojs=[0,0]. The stock realcell driver ran V_NL OFF (local-only
PP) — so the structure-factor fix only bites once projectors are actually on.
This variant exercises the now-phased V_NL through the full SCF→compose→λ
pipeline for the CaH6 gate A/B.

eig probe: repoint DECK to the local CaH6_NC deck (was a foreign-host path).

Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>
… SCF iters=32)

Bare full-basis vertex gate (cah6_fullbz_xval, npw_cap=0 n=645, nq=4, screened=0)
re-run on the V_NL structure-factor-fixed (bound) SCF. SCF converged iters=32
e_band=-65.1002 Ha. λ=20.0374 vs QE 4.376 (rel-ε=357.9%) — the bound eigenstates
move λ AWAY from QE (4.137→20.04), NOT toward. Hybrid xval selftest still PASS
(rel-ε=1.6524e-7). nq2 cross-check pending.

Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>
… byte-identical to nq4)

nq=2 reproduces nq=4 EXACTLY (bound SCF iters=32 e_band=-65.1002 Ha, λ=20.0374,
357.892% off QE 4.376) — λ=20 is the converged-eigenstate result, NOT a q-mesh artifact.
Finalized VERDICT.md + nq2 log.

Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant