{"id":"z1-barrier-campaign-round4","title":"Z1 Barrier-Accurate MLIP Campaign — Round-4 Results (Two Precision Chains)","subtitle":"Refuted twice — all four models miss the 40 meV gate (135–243 meV) identically under float32 and float64; systematic under-prediction.","category":"validation","tags":["z1","barriers","neb","batteries","round-4","featured"],"source":"articles/docs/validation/z1-barrier-campaign-round4-results.md","lang":"en","words":608,"readMinutes":3,"toc":[{"depth":2,"text":"Preregistered question","id":"preregistered-question"},{"depth":2,"text":"Panel and protocol","id":"panel-and-protocol"},{"depth":2,"text":"Results: both precision chains","id":"results-both-precision-chains"},{"depth":2,"text":"The shape of the failure","id":"the-shape-of-the-failure"},{"depth":2,"text":"Receipts","id":"receipts"}],"html":"<h1 id=\"z1-barrier-accurate-mlip-campaign-round-4-results-two-precision-chains\">Z1 Barrier-Accurate MLIP Campaign — Round-4 Results (Two Precision Chains)</h1><p><strong>Status:</strong> completed campaign · verdict <strong>refuted</strong> — independently on two precision chains\n<strong>Campaign:</strong> <code>discovery.round-4.z1-barriers.v1</code> · executed 2026-07-19 (float32 chain and float64 chain, same locked panel)\n<strong>Gate:</strong> migration-barrier MAE ≤ 40 meV against published DFT-NEB references across 30 chemistry-held-out paths\n<strong>Verdict:</strong> all four available models FAIL, 3.4–6× over threshold, with a systematic under-prediction bias — unchanged under float64</p>\n<h2 id=\"preregistered-question\">Preregistered question</h2><p>Can the declared available foundation MLIPs (chgnet 0.4.2; mace-torch 0.3.16 small / medium / mpa-0-medium) predict solid-state ion migration barriers accurately enough for battery-electrolyte screening — a mean absolute error of 40 meV or better against DFT-NEB references on 30 held-out chemical systems?</p>\n<h2 id=\"panel-and-protocol\">Panel and protocol</h2><ul>\n<li><strong>Reference panel:</strong> <code>data/candidates/z1_nebdft2k_barriers.lock.json</code> (SHA-256 <code>192fe54a…</code>) — 30 chemistry-held-out DFT-NEB paths, one deterministic path per official LiTraj nebDFT2k test chemistry (<em>npj Computational Materials</em>, 2025, DOI <a href=\"https://doi.org/10.1038/s41524-025-01571-z\">10.1038/s41524-025-01571-z</a>; source revision and archive hash pinned; deterministic rebuild verified byte-for-byte). DFT-relaxed endpoints/saddles/full energy profiles, reference barriers 0.068–3.251 eV, frozen CI-NEB protocol (FIRE, climb, improved tangent, k = 5 eV/Å², fmax 0.1 eV/Å).</li>\n<li><strong>Basis honesty:</strong> references are published DFT-NEB, not experiment. Barrier convention is max(image energy) − min(image energy) on both the reference and prediction sides; the panel builder validates the dataset&#39;s own barrier against the profile max−min within 0.5 meV.</li>\n<li><strong>Execution:</strong> per model, one Cloud Run cell on an isolated Round-4 job (NVIDIA L4), CI-NEB per path, per-path predicted barriers/signed errors/failures recorded <strong>without imputation</strong>; manifest and panel SHA validation fail closed. Images <code>z1-barrier-20260719</code> (float32 chain) and <code>z1-barrier-f64r2-20260719</code> (float64 chain).</li>\n</ul>\n<h2 id=\"results-both-precision-chains\">Results: both precision chains</h2><div class=\"table-wrap\"><table><thead><tr>\n<th>Model</th>\n<th>float32 MAE</th>\n<th>float64 MAE</th>\n<th>Paths (both chains)</th>\n<th>Gate</th>\n</tr>\n</thead><tbody><tr>\n<td data-label=\"Model\">mace-mpa-0-medium</td>\n<td data-label=\"float32 MAE\">135.0 meV</td>\n<td data-label=\"float64 MAE\">135.0 meV</td>\n<td data-label=\"Paths (both chains)\">28/30</td>\n<td data-label=\"Gate\">✗</td>\n</tr>\n<tr>\n<td data-label=\"Model\">mace-mp-small</td>\n<td data-label=\"float32 MAE\">152.0 meV</td>\n<td data-label=\"float64 MAE\">151.9 meV</td>\n<td data-label=\"Paths (both chains)\">26/30</td>\n<td data-label=\"Gate\">✗</td>\n</tr>\n<tr>\n<td data-label=\"Model\">mace-mp-medium</td>\n<td data-label=\"float32 MAE\">174.7 meV</td>\n<td data-label=\"float64 MAE\">174.7 meV</td>\n<td data-label=\"Paths (both chains)\">29/30</td>\n<td data-label=\"Gate\">✗</td>\n</tr>\n<tr>\n<td data-label=\"Model\">chgnet</td>\n<td data-label=\"float32 MAE\">242.5 meV</td>\n<td data-label=\"float64 MAE\">242.5 meV</td>\n<td data-label=\"Paths (both chains)\">28/30</td>\n<td data-label=\"Gate\">✗</td>\n</tr>\n</tbody></table></div><p>The first chain ran at float32 — against MACE vendor guidance for geometry optimization — a protocol defect found in post-execution review, fixed once (barrier row only; PR #32), and re-measured end-to-end at float64. <strong>The verdict is precision-independent</strong>: per-model MAEs are identical to within 0.1 meV. The miss is model error, not numerical artifact. Both chains stand as executed evidence with separate artifact prefixes (<code>z1/campaign/</code> and <code>z1/campaign-float64/</code>); checkpoint contexts bind calculator dtype (PR #36), so no float32 prediction can leak into a float64 record.</p>\n<h2 id=\"the-shape-of-the-failure\">The shape of the failure</h2><ul>\n<li><strong>Systematic under-prediction.</strong> For mace-mp-small (float32 chain), all 26 completed paths have <em>negative</em> signed error (−13 to −467 meV): the models under-predict migration barriers — consistent with training distributions dominated by near-equilibrium structures and with this program&#39;s Z3 finding of underbound interfaces. One systematic direction, two observables.</li>\n<li><strong>Convergence failures are honest, not hidden.</strong> 1–4 paths per model (the largest systems, 87–191 atoms) failed CI-NEB convergence under the frozen protocol and are recorded as failures — MAEs are computed on completed paths only, and <code>measurement_complete</code> is false everywhere because the protocol demands all 30.</li>\n<li><strong>Precision note.</strong> The f64 signed-error distribution partially relaxes (17/26 negative, mean −23.7 meV for mace-mp-small) while the MAE is unmoved — per-path pairing analysis is registered follow-up.</li>\n<li><strong>Precedent replicated.</strong> This confirms at 30-chemistry scale the Round-3 five-compound result (77.1 meV) whose claim was withdrawn: foundation MLIPs are not barrier-accurate, and the error is structural.</li>\n</ul>\n<h2 id=\"receipts\">Receipts</h2><ul>\n<li>Manifest: <code>campaigns/v1/z1.campaign-manifest.v1.json</code> (content hash <code>sha256:0a85044c…</code>, pins the panel lock); claim <code>registry/claims/discovery.z1.barrier-accuracy.v1.json</code> (unsupported pending ingestion; baseline bundle seeded).</li>\n<li>Panel + builder: <code>data/candidates/z1_nebdft2k_barriers.lock.json</code> (+ <code>.sha256</code>), <code>tools/build_z1_barrier_panel.py</code>.</li>\n<li>Raw artifacts: <code>gs://shed-489901-atlas-outputs/z1/campaign/&lt;model&gt;/</code> and <code>gs://shed-489901-atlas-outputs/z1/campaign-float64/&lt;model&gt;/</code> (<code>cell_result.json</code>, <code>cell_checkpoint.json</code>; execution metadata records <code>calculator_dtype</code>).</li>\n<li>Measurement rows: <code>data/candidates/z1/measurements.jsonl</code> (four RFC 8785 hash-chained aggregate fail rows, <code>tools/build_z1_measurement_rows.py</code>); the float64 chain has its own parallel set at <code>data/candidates/z1/f64/measurements.jsonl</code> (+ <code>f64/artifact-manifest.json</code>, built with <code>--source-root data/candidates/z1/f64 --result-uri-root gs://shed-489901-atlas-outputs/z1/campaign-float64</code>).</li>\n<li>Runner hardening from review (PRs #32, #36): float64 for barrier geometry optimization; non-finite energies fail the path; checkpoint dtype binding.</li>\n</ul>\n"}