REVIEW 2 major objections 4 minor 43 references
Code-space recovery for sample-based quantum diagonalization beyond native symmetry constraints
T0 review · 2 major / 4 minor · reviewed 2026-07-14 · grok-4.5
Pith's one-line read Dual-rail code-space recovery lets sample-based quantum diagonalization repair noisy samples without native symmetries, beating raw unencoded supports on Ising models up to 36 spins.
desk verdict Clear methods result: dual-rail recovery beats a strong unencoded sample-support baseline on symmetry-free Ising models, with fair matched hardware design and ordinary near-term limits. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Dual-rail code-space recovery: map each logical qubit as |0⟩→|01⟩ and |1⟩→|10⟩, implement sampling operations that intertwine with this encoding, then repair invalid pairs (00 or 11) stochastically using cluster-specific, self-consistently updated rail-occupation reference vectors before decoding and classical projected diagonalization.
What would settle it
On the same device, SqDRIFT protocol, and shot budget, if the encoded workflow’s lowest Ritz energy at target basis size 2.5 million were higher than the unencoded full-support Ritz energy on the reported 25- or 36-site Ising benchmarks, the central claim would be false.
Extended reading notes
Core claim
For Ising Hamiltonians without a native measurement-visible symmetry, dual-rail code-space recovery produces lower projected Ritz energies than direct diagonalization on the full unencoded sample support, even when the recovered logical basis is smaller than that support and the encoded circuits use more qubits and more two-qubit gates.
Load-bearing premise
Noise from the larger, deeper encoded circuits still leaves pair-violation patterns that self-consistent repair can turn into a better logical basis than the raw unencoded samples provide.
Editorial extensions
If this is right
- SQD-style projected diagonalization can be aimed at general Hermitian eigenvalue problems whose low-lying states are sparse in a computational basis, not only chemistry Hamiltonians with particle-number symmetry.
- Encoding overhead need not erase sample quality if the engineered pair constraint supplies enough repair information to improve the projected subspace.
- Native constraints and engineered code-space constraints can be composed: repair dual-rail pairs first, then enforce a logical sector such as particle number.
- Any sampling circuit whose operations admit encoded lifts that preserve the code space in the noiseless limit can feed the same recovery loop.
Reading between the lines
- On noisier devices or encodings, pair statistics may randomize faster than recovery can correct, erasing the subspace-quality advantage.
- Other error-detecting codes with local, measurement-visible constraints could replace dual-rail if their encoded operations remain practical.
- The residual gap to DMRG on the 36-site tests may close with denser Krylov sampling or larger shot budgets rather than only better recovery.
- The work reframes simple encodings as classical post-processing aids for hybrid subspace methods, not only as quantum memories.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes code-space recovery for sample-based quantum diagonalization (SQD): dual-rail encoding maps each logical qubit to a physical pair (|0⟩→|01⟩, |1⟩→|10⟩) so that Hamming-weight-one violations in noisy samples become detectable recovery signals. Encoded counterparts of logical sampling operations are implemented so that the pair constraint is preserved in the noiseless limit (Eqs. 1–4, 20–23). After sampling, invalid pairs are repaired by a self-consistent, cluster-adaptive stochastic procedure guided by reference vectors, decoded logical bitstrings span a projected subspace, and the original logical Hamiltonian is diagonalized. Benchmarks on 1D/2D TFIM and MFIM instances (n=25 and 36) use matched SqDRIFT logical sequences on IBM Heron r2. Across all five problems, recovered subspaces at D_tar ≤ 2.5×10^6 yield lower Ritz energies and smaller |E−E_ref| than the full unencoded sample-support baseline (Table 1, Figs. 2–4), despite roughly doubled qubits and higher transpiled depth/N2q (Table 2).
Significance. If the empirical advantage holds more broadly, the work supplies a concrete route to apply SQD-style recovery to eigenvalue problems that lack a native, measurement-visible constraint such as particle number. The matched logical SqDRIFT design, identical shot budgets, and Rayleigh–Ritz projection of the original logical Hamiltonian make the comparison fair and falsifiable; energy-variance intercepts for the 25-site cases agree closely with exact/DMRG references, strengthening the claim that the recovered subspaces track the target low-energy branch. The explicit resource trade-off (Table 2) and the demonstration that recovery can still win when fully valid encoded strings are rare (0.02–4.8 %) are useful contributions to the quantum-centric diagonalization literature.
major comments (2)
- The central claim rests on a single sampling protocol (SqDRIFT with fixed K, Δt, M_seq) and one IBM Heron device. Section 2.2 and Methods 4.1 do not test whether the subspace-quality advantage survives under alternative samplers (e.g., different Krylov constructions or non-randomized circuits) or under substantially different noise. A short additional experiment or a clear statement that the advantage is protocol- and device-specific would make the scope of the claim precise.
- For the 36-site stress tests (Fig. 4), residual gaps to DMRG remain large (≈0.22 and ≈0.17 at D_tar=2.5×10^6) and the manuscript reports that a controlled zero-variance extrapolation was not possible. The Discussion attributes this to sampling budget and 72-qubit noise exposure, but does not quantify how much of the residual is due to incomplete support versus residual recovery error. Clarifying this distinction (e.g., by comparing recovered support overlap with a classical high-quality basis when available) would strengthen the interpretation of the stress-test results.
minor comments (4)
- Table 2 reports median (min–max) active qubits, depth, and N2q; a brief note on how routing/layout inflation of active qubits affects the encoded vs unencoded comparison would help readers interpret the overhead.
- The modified-ReLU parameters (ρ,δ)=(1/2,0.01) and N_C=2 are fixed without a sensitivity check (Methods 4.2). A short remark or supplementary scan would reassure that the advantage is not finely tuned to these choices.
- Fig. 1 is dense; labeling the six numbered stages more explicitly in the caption would improve readability for readers new to SQD recovery loops.
- The Discussion correctly notes that engineered and native constraints can be complementary, but does not cite or discuss any prior dual-rail or pair-code use outside photonic QC; a sentence situating the recovery-only use of dual-rail would help.
Circularity Check
No significant circularity: central claim is an empirical encoded-vs-unencoded comparison against independent hardware baselines and external exact/DMRG references.
full rationale
The paper's load-bearing claim is that dual-rail code-space recovery produces lower projected Ritz energies than matched unencoded sample-support diagonalization on Ising Hamiltonians lacking native U(1) constraints (Abstract; Table 1; Figs. 2–4). That comparison is not forced by definition: both workflows use the same logical SqDRIFT term sequences and shots; the unencoded baseline diagonalizes over the full observed logical support (plus known GHZ support); the encoded workflow recovers logical bases via the engineered pair constraint and then projects the original logical Hamiltonian. Rayleigh–Ritz is stated explicitly (Methods §4.2: λ_b^(t) ≥ λ_min(O) independently of how the basis was built), so reported energies remain variational upper bounds, not tautologies of the recovery map. Self-consistent reference updates from the current global-best Ritz state are iterative subspace refinement, not a self-definitional loop that equates input to output. Energy-variance intercepts are diagnostics only. Self-citations ([21] modified-ReLU scores; [22] cluster-adaptive partitioning) supply methodological choices and are not uniqueness theorems that forbid alternatives or force the Table 1 ranking. Dual-rail encoding and SqDRIFT are standard external constructions. The result is therefore self-contained against external classical references and an independent hardware baseline; no step reduces the central claim to its own inputs by construction.
Assumptions & free parameters
free parameters (5)
- modified-ReLU recovery parameters (ρ, δ) =
(1/2, 0.01)
- number of recovery clusters N_C =
2
- carry-over threshold ε_carry =
10^{-3}
- qDRIFT sequence length M_seq and Krylov parameters =
M_seq=75/100, K=4, Δt=1.0
- target projected dimensions D_tar and batch count N_B =
D_tar up to 2.5e6, N_B=10
assumptions (4)
- standard math Encoded operators A' satisfy A'V = V A so that ideal encoded circuits stay in the dual-rail code space and reproduce the logical sampling distribution after decoding.
- standard math Rayleigh–Ritz: the lowest eigenvalue of the projected Hamiltonian is a variational upper bound on the true ground-state energy, independent of how the basis was obtained.
- domain assumption Under ground-state concentration and sufficient Krylov sampling budget, important computational-basis states appear with non-negligible probability in at least one Krylov state (SqDRIFT/qDRIFT sampling assumption).
- ad hoc to paper Noise-induced dual-rail pair violations remain partially structured so that reference-vector-guided repair improves logical support quality relative to unencoded samples.
invented entities (1)
-
code-space recovery workflow (dual-rail + self-consistent pair repair for SQD)
Cite this review
Pith. "Pith review of Code-space recovery for sample-based quantum diagonalization beyond native symmetry constraints." pith.science (2026). https://pith.science/paper/MXXV2BPL
@misc{pith2026260710227,
author = {Pith},
title = {Pith review of: Code-space recovery for sample-based quantum diagonalization beyond native symmetry constraints},
year = {2026},
howpublished = {\url{https://pith.science/paper/MXXV2BPL}},
note = {Machine review of arXiv:2607.10227}
}
abstract
Sample-based quantum diagonalization (SQD) diagonalizes a Hamiltonian in a compact subspace built from quantum samples, and its performance often relies on recovery procedures that exploit native constraints such as particle-number symmetry. For a broad class of eigenvalue problems, however, no analogous constraint is guaranteed, limiting the applicability of SQD-type recovery. Here, we introduce code-space recovery, which engineers recoverable structure through encoding rather than assuming it in the target problem. Using a dual-rail representation, each logical qubit is mapped to a physical pair, $|0\rangle \to |01\rangle$ and $|1\rangle \to |10\rangle$, making code-space violations in noisy samples detectable and repairable. We combine this encoding with self-consistent recovery and benchmark it on transverse- and mixed-field Ising models with up to 36 spin sites. Despite increased circuit overhead, code-space recovery yields lower projected Ritz energies than unencoded sample-support diagonalization even at smaller projected-basis dimensions, suggesting that engineered recoverable structure can extend SQD beyond native constraints.
Reference graph
Works this paper leans on
-
[1]
& van der Vorst, H.Templates for the solution of algebraic eigenvalue problems: a practical guide(SIAM, Philadelphia, PA, 2000)
Bai, Z., Demmel, J., Dongarra, J., Ruhe, A. & van der Vorst, H.Templates for the solution of algebraic eigenvalue problems: a practical guide(SIAM, Philadelphia, PA, 2000)
2000
-
[2]
Saad, Y.Numerical methods for large eigenvalue problems: revised edition(SIAM, Philadelphia, PA, 2011)
2011
-
[3]
Nielsen, M. A. & Chuang, I. L.Quantum computation and quantum information: 10th anniversary edition(Cambridge University Press, Cambridge, 2010). 17
2010
-
[4]
Abrams, D. S. & Lloyd, S. Quantum algorithm providing exponential speed increase for finding eigenvalues and eigenvectors.Physical Review Letters83, 5162 (1999)
1999
-
[5]
Cao, Y.et al.Quantum chemistry in the age of quantum computing.Chemical Reviews 119, 10856–10915 (2019)
2019
-
[6]
& Rebentrost, P
Lloyd, S., Mohseni, M. & Rebentrost, P. Quantum principal component analysis.Nature Physics10, 631–633 (2014)
2014
-
[7]
Biamonte, J.et al.Quantum machine learning.Nature549, 195–202 (2017)
2017
-
[8]
Kitaev, A. Y. Quantum measurements and the Abelian stabilizer problem.arXiv preprint quant-ph/9511026(1995)
arXiv 1995
Show all 43 references
-
[9]
& Mosca, M
Cleve, R., Ekert, A., Macchiavello, C. & Mosca, M. Quantum algorithms revisited.Pro- ceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences454, 339–354 (1998)
1998
-
[10]
Peruzzo, A.et al.A variational eigenvalue solver on a photonic quantum processor.Nature Communications5, 4213 (2014)
2014
-
[11]
R., Romero, J., Babbush, R
McClean, J. R., Romero, J., Babbush, R. & Aspuru-Guzik, A. The theory of variational hybrid quantum-classical algorithms.New Journal of Physics18, 023023 (2016)
2016
-
[12]
Cerezo, M.et al.Variational quantum algorithms.Nature Reviews Physics3, 625–644 (2021)
2021
-
[13]
Yoshioka, N.et al.Krylov diagonalization of large many-body Hamiltonians on a quantum processor.Nature Communications16, 5014 (2025)
2025
-
[14]
M., Wecker, D
Reiher, M., Wiebe, N., Svore, K. M., Wecker, D. & Troyer, M. Elucidating reaction mechanisms on quantum computers.Proceedings of the National Academy of Sciences 114, 7555–7560 (2017)
2017
-
[15]
D.et al.Improved fault-tolerant quantum simulation of condensed-phase correlated electrons via trotterization.Quantum4, 296 (2020)
Kivlichan, I. D.et al.Improved fault-tolerant quantum simulation of condensed-phase correlated electrons via trotterization.Quantum4, 296 (2020)
2020
-
[16]
R., Boixo, S., Smelyanskiy, V
McClean, J. R., Boixo, S., Smelyanskiy, V. N., Babbush, R. & Neven, H. Barren plateaus in quantum neural network training landscapes.Nature Communications9, 4812 (2018)
2018
-
[17]
Wang, S.et al.Noise-induced barren plateaus in variational quantum algorithms.Nature Communications12, 6961 (2021)
2021
-
[18]
& Coles, P
Holmes, Z., Sharma, K., Cerezo, M. & Coles, P. J. Connecting ansatz expressibility to gradient magnitudes and barren plateaus.PRX Quantum3, 010313 (2022)
2022
-
[19]
F.et al.Measurements as a roadblock to near-term practical quantum advantage in chemistry: Resource analysis.Physical Review Research4, 033154 (2022)
Gonthier, J. F.et al.Measurements as a roadblock to near-term practical quantum advantage in chemistry: Resource analysis.Physical Review Research4, 033154 (2022)
2022
-
[20]
Kanno, K.et al.Quantum-selected configuration interaction: Classical diagonalization of Hamiltonians in subspaces selected by quantum computers.Physical Review Research8, 023268 (2026)
2026
-
[21]
Robledo-Moreno, J.et al.Chemistry beyond the scale of exact diagonalization on a quantum-centric supercomputer.Science Advances11, eadu9991 (2025)
2025
-
[22]
Park, B.et al.Cluster-adaptive sample-based quantum diagonalization for strongly correlated systems.arXiv preprint arXiv:2603.09346(2026). 18
2026 arXiv
-
[23]
Yu, J.et al.Quantum-centric algorithm for sample-based Krylov diagonalization.arXiv preprint arXiv:2501.09702(2025)
2025
-
[24]
Piccinelli, S.et al.Quantum chemistry with provable convergence via randomized sample- based Krylov quantum diagonalization.arXiv preprint arXiv:2508.02578(2025)
2025
-
[25]
& Motta, M
Barison, S., Robledo-Moreno, J. & Motta, M. Quantum-centric computation of molecular excited states with extended sample-based quantum diagonalization.Quantum Science and Technology10, 025034 (2025)
2025
-
[26]
J., Motta, M
Danilov, D., Robledo-Moreno, J., Sung, K. J., Motta, M. & Shee, J. Enhancing the accuracy and efficiency of sample-based quantum diagonalization with phaseless auxiliary- field quantum Monte Carlo.Journal of Chemical Theory and Computation21, 11585– 11594 (2025)
2025
-
[27]
Shajan, A.et al.Toward quantum-centric simulations of extended molecules: Sample- based quantum diagonalization enhanced with density matrix embedding theory.Journal of Chemical Theory and Computation21, 6801–6810 (2025)
2025
-
[28]
Patra, C.et al.Physics-informed generative machine learning for accelerated quantum- centric supercomputing.arXiv preprint arXiv:2512.06858(2025)
2025
-
[29]
Ronˇ cevi´ c, I.et al.A molecule with half-M¨ obius topology.Science392, eaea3321 (2026)
2026
-
[30]
M.et al.Crossing the 12,000-atom barrier with heterogeneous quantum- classical supercomputing: quantum chemistry of protein-ligand complexes.arXiv preprint arXiv:2605.01138(2026)
Merz Jr, K. M.et al.Crossing the 12,000-atom barrier with heterogeneous quantum- classical supercomputing: quantum chemistry of protein-ligand complexes.arXiv preprint arXiv:2605.01138(2026)
2026 arXiv
-
[31]
Chuang, I. L. & Yamamoto, Y. Simple quantum computer.Physical Review A52, 3489 (1995)
1995
-
[32]
& Milburn, G
Knill, E., Laflamme, R. & Milburn, G. J. A scheme for efficient quantum computation with linear optics.Nature409, 46–52 (2001)
2001
-
[33]
Kok, P.et al.Linear optical quantum computing with photonic qubits.Reviews of Modern Physics79, 135–174 (2007)
2007
-
[34]
Random compiler for fast Hamiltonian simulation.Physical Review Letters 123, 070503 (2019)
Campbell, E. Random compiler for fast Hamiltonian simulation.Physical Review Letters 123, 070503 (2019)
2019
-
[35]
White, S. R. Density matrix formulation for quantum renormalization groups.Physical Review Letters69, 2863 (1992)
1992
-
[36]
quimb: a python package for quantum information and many-body calculations
Gray, J. quimb: a python package for quantum information and many-body calculations. Journal of Open Source Software3, 819 (2018)
2018
-
[37]
& Imada, M
Kashima, T. & Imada, M. Path-integral renormalization group method for numerical study on ground states of strongly correlated electronic systems.Journal of the Physical Society of Japan70, 2287–2299 (2001)
2001
-
[38]
Javadi-Abhari, A.et al.Quantum computing with Qiskit.arXiv preprint arXiv:2405.08810(2024)
2024 arXiv
-
[39]
& Lloyd, S
Viola, L. & Lloyd, S. Dynamical suppression of decoherence in two-state quantum systems. Physical Review A58, 2733 (1998)
1998
-
[40]
D., Kang, H., Sundaresan, N
Nation, P. D., Kang, H., Sundaresan, N. & Gambetta, J. M. Scalable mitigation of measurement errors on quantum computers.PRX Quantum2, 040326 (2021). 19
2021
-
[41]
& Govaert, G
Nadif, M. & Govaert, G. Clustering for binary data and mixture models—choice of the model.Applied Stochastic Models and Data Analysis13, 269–278 (1997)
1997
-
[42]
& McCombs, J
Stathopoulos, A. & McCombs, J. R. PRIMME: preconditioned iterative multimethod eigensolver—methods and software description.ACM Transactions on Mathematical Software (TOMS)37, 21:1–21:30 (2010)
2010
-
[43]
Morin, S.et al.StepMix: A Python package for pseudo-likelihood estimation of generalized mixture models with external variables.Journal of Statistical Software113, 1–39 (2025). 20
2025
Reviewed July 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.