Imagine you're a factory inspector checking that every widget rolling off the line matches the blueprint. You could measure every dimension of every widget — expensive. Or you could develop a clever sampling strategy that checks just enough features, in just the right order, that you can certify with confidence the widget is correct. This paper does that for fermionic Gaussian states: the quantum analog of checking a blueprint, where the 'widgets' are quantum states and the 'features' are single-qubit measurements chosen adaptively. The committed claim: for any d-mode pure fermionic Gaussian state, you can certify it using O(d²/ε) copies, where ε is your desired precision. The ε⁻¹ dependence is provably optimal — no measurement strategy, not even one using entangled measurements across copies, can do better. Previous protocols either scaled as ε⁻² (from fidelity estimation), worked only for random states, or required entangled measurements across many qubits. This paper gets optimal scaling with single-qubit measurements, measuring only one qubit adaptively per copy. The engine under the hood is beautiful. The sample complexity is controlled by the spectral gap of a d↔d−2 down-up walk — a Markov chain studied in high-dimensional expander theory, equivalent to two-site Glauber dynamics in statistical physics. The connection is not cosmetic: the worst case for certification corresponds to the slowest-mixing instance of this walk. The authors show the worst-case bound O(d²/ε) is tight for their protocol and is attained by physically relevant states — ground states of the fully dimerized SSH chain and BCS pair states — not contrived adversarial examples. On the ladder: the main comparison is to prior certification protocols. Direct fidelity estimation gives O(d²/ε²) — this paper shaves an entire factor of ε⁻¹. O'Donnell and Wright's 2015 tomography lower bound of Ω(d²/ε) for the related state certification problem confirms the optimality. The Haar-typical case is even better: numerics up to d=14 suggest O(d/ε) copies suffice, a factor-of-d improvement, though this remains a conjecture. Integrity is strong for a theory paper. The ε⁻¹ optimality is a mathematical proof, not a simulation. The worst-case tightness is demonstrated by exhibiting explicit state families (SSH, BCS) that saturate the bound. The Haar-typical improvement is supported by numerics to d=14 but flagged honestly as unproven. Classical processing per copy is O(d³), which is practical. No code release is mentioned, but the protocol is fully specified and reproducible from the math. The paper extends beyond Gaussian states to 'phase-dressed' states — fermionic Gaussian states with an efficiently computable diagonal phase injected. This class includes non-Gaussian magic states for matchgate computation, which is a meaningful extension: it says the certification protocol works for states that are computationally more powerful than the Gaussian class, without additional sample cost. The same O(d²/ε) bound holds because it depends only on the computational-basis distribution, which the phase dressing leaves unchanged. The successor question is clear: proving the O(d/ε) Haar-typical conjecture, scaling the numerics past d=14, and extending to mixed states or noisy settings. The authors test numerics only to d=14, which is likely a compute limitation for the random sampling involved. The mixed-state extension is harder and may be the next paper. For practitioners building quantum simulation pipelines, this result means certification is no longer the bottleneck — the sample cost scales linearly in precision, matching the information-theoretic floor.
