pith. sign in

arxiv: 2606.31952 · v1 · pith:INOFUY3Unew · submitted 2026-06-30 · 🪐 quant-ph

An efficient Pauli decomposition algorithm for structured matrices

Pith reviewed 2026-07-01 04:57 UTC · model grok-4.3

classification 🪐 quant-ph
keywords Pauli decompositionsparse query accessrandomized algorithmquery complexityquantum algorithmsstructured matricesclassical preprocessingpolynomial time
0
0 comments X

The pith

A randomized classical algorithm recovers the exact Pauli decomposition of matrices with only polynomially many nonzero Pauli terms using sparse query access.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a randomized classical algorithm that, given a matrix promised to have support on only k = poly(n) Pauli strings and accessible through classical sparse queries, recovers the exact coefficients and strings with success probability at least 1 - δ. This matters because the generic Pauli decomposition problem requires time exponential in n for arbitrary matrices, yet many near-term quantum algorithms rely on structured classical inputs that satisfy the sparsity promise. The new algorithm exploits that promise to achieve query and runtime complexity polynomial in n, k, and log(1/δ).

Core claim

We present a randomized classical algorithm that recovers the exact Pauli decomposition with success probability at least 1 - δ, for any δ. Under the stated access model, the algorithm executes with query and runtime complexity that is polynomial in n, k, and log(1/δ). These results show that, even though finding the Pauli decomposition is exponentially hard for general matrices, it becomes efficiently solvable for matrices that are known to be sparse in the Pauli basis, a regime that is relevant to near-term quantum algorithms operating on structured classical input.

What carries the argument

Randomized sampling that queries matrix entries in the computational basis to identify the k supporting Pauli strings without exhaustive search over 4^n possibilities.

If this is right

  • Exact recovery succeeds with probability at least 1 - δ for any chosen δ.
  • Query and runtime scale polynomially in n, k, and log(1/δ) rather than exponentially in n.
  • The procedure directly applies to structured inputs that arise in near-term quantum algorithms.
  • Sparsity in the Pauli basis removes the exponential bottleneck that exists for dense matrices.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same sparsity promise could be leveraged for efficient classical preprocessing steps in hybrid quantum-classical workflows that prepare or verify input states.
  • Similar randomized sampling ideas might apply to decomposition tasks in other operator bases beyond Pauli strings.
  • Practical implementations on random sparse matrices of moderate size would allow direct measurement of observed scaling with n and k.

Load-bearing premise

The input matrix is guaranteed to have nonzero support on only polynomially many Pauli strings and is supplied through classical sparse query access.

What would settle it

Run the algorithm on an explicit matrix with k = n^2 nonzero Pauli coefficients; if it either fails to recover the exact decomposition with probability greater than δ or requires superpolynomial queries or time, the central claim is false.

Figures

Figures reproduced from arXiv: 2606.31952 by Alexey V. Gorshkov, Daniel J. Spencer, Kishor Bharti.

Figure 1
Figure 1. Figure 1: Schematic of the Pauli decomposition algorithm ( [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Runtime of our algorithm for 5 random instances each for values [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Query count of our algorithm for n = 2, . . . , 30. and so making a polynomial number of random checks would very likely miss this single nonzero entry. There are a few settings where our algorithm can be particularly powerful. The first is when the matrix is specified implicitly by an efficient classical procedure that can answer the required sparse row or column queries without needing to construct the f… view at source ↗
read the original abstract

Decomposing classical matrices into linear combinations of Pauli strings is a major bottleneck for end-to-end implementations of near-term quantum algorithms. In this work, we consider a promise version of this Pauli decomposition problem in which the matrix is guaranteed to have support on only $k = \mathsf{poly}(n)$ Pauli strings and is given through classical sparse query access. Existing Pauli decomposition algorithms are designed for the generic, dense problem and do not inherently take advantage of this promised sparsity, so these approaches take time that is exponential in $n$. We present a randomized classical algorithm that does take advantage of this sparsity and recovers the exact Pauli decomposition with success probability at least $1 - \delta$, for any $\delta$. Under the stated access model, the algorithm executes with query and runtime complexity that is polynomial in $n$, $k$, and $\log(1/\delta)$. These results show that, even though finding the Pauli decomposition is exponentially hard for general matrices, it becomes efficiently solvable for matrices that are known to be sparse in the Pauli basis, a regime that is relevant to near-term quantum algorithms operating on structured classical input.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The paper presents a randomized classical algorithm for the promise version of the Pauli decomposition problem: given an n-qubit matrix with exact support on k = poly(n) Pauli strings and classical sparse query access, the algorithm recovers the exact coefficients with success probability at least 1-δ in query and runtime complexity polynomial in n, k, and log(1/δ).

Significance. If the stated complexity and correctness claims hold, the result is significant for near-term quantum algorithms because it shows that Pauli decomposition, which is exponentially hard in the dense case, becomes tractable under an explicit sparsity promise that matches many structured classical inputs. The explicit access model and randomized guarantee with tunable δ are strengths; the work correctly distinguishes the promise setting from the generic dense case.

major comments (1)
  1. [Abstract / full text] The abstract and reader's summary state clear polynomial query/runtime bounds and success probability, but the provided manuscript text contains only the abstract and no algorithm pseudocode, proof of correctness, or analysis of edge cases (e.g., coefficient estimation under the sparse query model or handling of zero coefficients). Without these, the central claim cannot be verified as load-bearing for the polynomial bound.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive assessment of the work's significance and for identifying the issue with the provided manuscript content. We address the major comment below.

read point-by-point responses
  1. Referee: [Abstract / full text] The abstract and reader's summary state clear polynomial query/runtime bounds and success probability, but the provided manuscript text contains only the abstract and no algorithm pseudocode, proof of correctness, or analysis of edge cases (e.g., coefficient estimation under the sparse query model or handling of zero coefficients). Without these, the central claim cannot be verified as load-bearing for the polynomial bound.

    Authors: We acknowledge that the version of the manuscript provided for review contained only the abstract, which prevents verification of the central claims. This appears to be an error in the submission package. The complete manuscript includes a detailed description of the randomized algorithm (Section 3), pseudocode (Algorithm 1), a full proof of correctness and complexity (Theorem 1 and supporting lemmas in Section 4), and analysis of edge cases. In particular, the algorithm uses sparse query access to estimate coefficients via random sampling over the promised k-support and outputs only non-zero coefficients, naturally excluding zeros without additional handling. We will submit a revised version that incorporates all of these elements to make the polynomial bounds verifiable. revision: yes

Circularity Check

0 steps flagged

No significant circularity; direct algorithmic claim under explicit promise and access model

full rationale

The paper presents a randomized algorithm for exact Pauli decomposition under the promise that the input matrix has support on only k=poly(n) Pauli strings and is accessible via classical sparse queries. The central claim is a statement of query/runtime complexity polynomial in n, k, and log(1/δ) with success probability 1-δ. This is a standard promise-problem result in sparse recovery over the Pauli basis; no derivation step reduces by construction to fitted parameters, self-definitional equations, or load-bearing self-citations. The access model and sparsity promise are stated explicitly and independently of the algorithm's output. No self-citation chains or ansatzes are invoked to justify the core runtime bound. The result is self-contained against external benchmarks for sublinear sparse recovery.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

No free parameters, invented entities, or non-standard axioms are apparent from the abstract; the result rests on standard randomized algorithm techniques and query complexity models.

axioms (1)
  • standard math Standard assumptions of randomized algorithms (e.g., access to random bits) and classical query complexity models
    The algorithm is described as randomized and query-based, relying on these background models.

pith-pipeline@v0.9.1-grok · 5730 in / 1221 out tokens · 39241 ms · 2026-07-01T04:57:20.034374+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

49 extracted references · 12 canonical work pages · 3 internal anchors

  1. [1]

    silently

    0is the all-zero string. 4 Proof.Looking at the nonzero elementP(x, z) v,u, we re- quire each of (Z zj X xj)vj uj ̸= 0 for allj∈[n]. Then, looking at each term inP(x, z) individually, ifx j = 0, thenZ zj X xj =Z zj is diagonal and so the nonzero ele- ments are exactly specified byu j =v j. Sinceu j, vj ∈ {0,1}, this is equivalent tou j ⊕v j = 0. Likewise,...

  2. [2]

    residual

    The code was run with Python v. 3.14.0, NumPy v. 2.3.5, PennyLane v. 0.43.1, pandas v. 2.3.3, and matplotlib v. 3.10.7. The script used for these numerical experiments is publicly available athttps://github. com/dspencer2596/sparse-pauli-decomposition/ releases/tag/v1.0-paper. In Fig. 2, we compare the runtime of our algorithm with thePauliDecomposealgori...

  3. [3]

    Lanyon, J

    B. Lanyon, J. Whitfield, G. Gillett, M. Goggin, M. Almeida, I. Kassel, J. Biamonte, M. Mohseni, B. Pow- ell, M. Barbieri, A. Aspuru-Guzik, and A. White, To- wards quantum chemistry on a quantum computer, Na- ture Chemistry2, 106 (2010)

  4. [4]

    Y. Cao, J. Romero, J. P. Olson, M. Degroote, P. D. John- son, M. Kieferov´ a, I. D. Kivlichan, T. Menke, B. Per- opadre, N. P. Sawaya, S. Sim, L. Veis, and A. Aspuru- Guzik, Quantum chemistry in the age of quantum com- puting, Chemical Reviews119, 10856 (2019)

  5. [5]

    Bauer, S

    B. Bauer, S. Bravyi, M. Motta, and G. K.-L. Chan, Quantum algorithms for quantum chemistry and quan- tum materials science, Chemical Reviews120, 12685 (2020)

  6. [6]

    McArdle, S

    S. McArdle, S. Endo, A. Aspuru-Guzik, S. C. Benjamin, and X. Yuan, Quantum computational chemistry, Re- views of Modern Physics92, 015003 (2020)

  7. [7]

    Motta and J

    M. Motta and J. E. Rice, Emerging quantum comput- ing algorithms for quantum chemistry, Computational Molecular Science12, e1580 (2021)

  8. [8]

    Quantum algorithms for supervised and unsupervised machine learning

    S. Lloyd, M. Mohseni, and P. Rebentrost, Quantum algorithms for supervised and un- supervised machine learning, arXiv:1307.0411 https://doi.org/10.48550/arXiv.1307.0411 (2013)

  9. [9]

    Schuld, I

    M. Schuld, I. Sinayskiy, and F. Petruccione, An intro- duction to quantum machine learning, Contemporary Physics56, 172 (2014)

  10. [10]

    Biamonte, P

    J. Biamonte, P. Wittek, N. Pancotti, P. Rebentrost, N. Wiebe, and S. Lloyd, Quantum machine learning, Na- ture549, 195 (2017)

  11. [11]

    Dunjko and H

    V. Dunjko and H. J. Briegel, Machine learning & artifi- cial intelligence in the quantum domain: a review of re- cent progress, Reports on Progress in Physics81, 074001 (2018)

  12. [12]

    S. B. Ramezani, A. Sommers, H. K. Manchukonda, S. Rahimi, and A. Amirlatifi, Machine learning algo- rithms in quantum computing: a survey, 2020 Interna- tional Joint Conference on Neural Networks (ICNN) , 1 (2020)

  13. [13]

    Zaman, A

    K. Zaman, A. Marchisio, M. A. Hanif, and M. Shafique, A survey on quantum machine learning: current trends, challenges, opportu- nities, and the road ahead, arXiv:2310.10315 https://doi.org/10.48550/arXiv.2310.10315 (2023)

  14. [14]

    Wang and J

    Y. Wang and J. Liu, A comprehensive review of quantum machine learning: from NISQ to fault tolerance, Reports on Progress in Physics87, 116402 (2024)

  15. [15]

    Or´ us, S

    R. Or´ us, S. Mugel, and E. Lizaso, Quantum computing for finance: Overview and prospects, Reviews in Physics 4, 100028 (2019)

  16. [16]

    Bouland, W

    A. Bouland, W. van Dam, H. Joorati, I. Kereni- dis, and A. Prakash, Prospects and chal- lenges of quantum finance, arXiv2011.06492v1 https://doi.org/10.48550/arXiv.2011.06492 (2020)

  17. [17]

    Herman, C

    D. Herman, C. Googin, X. Liu, Y. Sun, A. Galda, I. Safro, M. Pistoia, and Y. Alexeev, Quantum computing for fi- nance, Nature Reviews Physics5, 450 (2023)

  18. [18]

    Weigold, J

    M. Weigold, J. Barzen, F. Leymann, and M. Salm, En- coding patterns for quantum algorithms, IET Quantum Communication2, 141 (2021)

  19. [19]

    Weigold, J

    M. Weigold, J. Barzen, F. Leymann, and M. Salm, Ex- panding data encoding patterns for quantum algorithms, 2021 IEEE 18th International Conference on Soft- ware Architecture Companion (ICSA-C) 10.1109/ICSA- C52384.2021.00025 (2021)

  20. [20]

    Nakaji, S

    K. Nakaji, S. Uno, Y. Suzuki, R. Raymond, T. Onodera, T. Tanaka, H. Tezuka, N. Mitsuda, and N. Yamamoto, Approximate amplitude encoding in shallow parameter- ized quantum circuits and its application to financial market indicators, Physical Review Research4, 023136 (2022)

  21. [21]

    Marin-Sanchez, J

    G. Marin-Sanchez, J. Gonzalez-Conde, and M. Sanz, Quantum algorithms for approximate function loading, Physical Review Research5, 033114 (2023)

  22. [22]

    Quantum computing in the NISQ era and beyond,

    J. Preskill, Quantum computing in the NISQ era and beyond, Quantum2, 10.22331/q-2018-08-06-79 (2018)

  23. [23]

    URL http://dx.doi.org/10.1103/RevModPhys.94.015004

    K. Bharti, A. Cervera-Lierta, T. H. Kyaw, T. Haug, S. Alperin-Lea, A. Anand, M. Degroote, H. Heimonen, J. S. Kottmann, T. Menke, W. K. Mok, S. Sim, L. C. Kwek, and A. Aspuru-Guzik, Noisy intermediate-scale quantum algorithms, Reviews of Modern Physics94, 10.1103/RevModPhys.94.015004 (2022)

  24. [24]

    Peruzzo, J

    A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O’Brien, A variational eigenvalue solver on a photonic quantum processor, Nature Communications5, 4213 (2014)

  25. [25]

    W. M. Kirby and P. J. Love, Variational quantum eigen- solvers for sparse Hamiltonians, Physical Review Letters 127, 110503 (2021)

  26. [26]

    Tilly, H

    J. Tilly, H. Chen, S. Cao, D. Picozzi, K. Setia, Y. Li, E. Grant, L. Wossnig, I. Rungger, G. H. Booth, and J. Tennyson, The variational quantum eigensolver: a re- view of methods and best practices, Physics Reports986, 1 (2022)

  27. [27]

    Patel, P

    D. Patel, P. J. Coles, and M. M. Wilde, Variational quan- tum algorithms for semidefinite programming, Quantum 8, 1374 (2024)

  28. [28]

    H. Y. Huang, K. Bharti, and P. Rebentrost, Near-term quantum algorithms for linear systems of equations with regression loss functions, New Journal of Physics23, 113021 (2021)

  29. [29]

    Bravo-Prieto, R

    C. Bravo-Prieto, R. LaRose, M. Cerezo, Y. Suba¸ si, L. Cincio, and P. J. Coles, Variational quantum linear solver, Quantum7, 1188 (2023). 15

  30. [30]

    Pellow-Jarman, I

    A. Pellow-Jarman, I. Sinayskiy, A. Pillay, and F. Petruc- cione, Near term algorithms for linear systems of equa- tions, Quantum Information Processing22, 258 (2023)

  31. [31]

    Bharti and T

    K. Bharti and T. Haug, Iterative quantum-assisted eigen- solver, Physical Review A104, L050401 (2021)

  32. [32]

    Bharti, T

    K. Bharti, T. Haug, V. Vedral, and L.-C. Kwek, Noisy in- termediate scale quantum algorithm for semidefinite pro- gramming, Physical Review A105, 052445 (2022)

  33. [33]

    Larocca, S

    M. Larocca, S. Thanasilp, S. Wang, K. Sharma, J. Bia- monte, P. J. Coles, L. Cincio, J. R. McClean, Z. Holmes, and M. Cerezo, Barren plateaus in variational quantum computing, Nature Reviews Physics7, 174 (2025)

  34. [34]

    PennyLane: Automatic differentiation of hybrid quantum-classical computations

    V. Bergholm, J. Izaac, M. Schuld, C. Gogolin, S. Ahmed, V. Ajith, M. S. Alam, G. Alonso-Linaje, B. Akash- Narayanan, A. Asadi, J. M. Arrazola, U. Azad, S. Ban- ning, C. Blank, T. R. Bromley, B. A. Cordier, J. Ceroni, A. Delgado, O. D. Matteo, A. Dusko, T. Garg, D. Guala, A. Hayes, R. Hill, A. Ijaz, T. Isacsson, D. Ittah, S. Ja- hangiri, P. Jain, E. Jiang,...

  35. [35]

    Gunlycke, M

    D. Gunlycke, M. C. Palenik, A. R. Emmert, and S. A. Fischer, Efficient algorithm for generating Pauli coordi- nates for an arbitrary linear operator, arXiv:2011.08942 https://doi.org/10.48550/arXiv.2011.08942 (2020)

  36. [36]

    S. V. Romero and J. Santos-Su´ arez, Paulicomposer: Compute tensor products of Pauli matrices efficiently, Quantum Information Processing22, 449 (2023)

  37. [37]

    Jones, Decomposing dense matrices into dense Pauli tensors, arXiv:2401.16378 https://doi.org/10.48550/arXiv.2401.16378 (2024)

    T. Jones, Decomposing dense matrices into dense Pauli tensors, arXiv:2401.16378 https://doi.org/10.48550/arXiv.2401.16378 (2024)

  38. [38]

    Koska, M

    O. Koska, M. Baboulin, and A. Gazda, A tree- approach Pauli decomposition algorithm with ap- plication to quantum computing, arXiv:2403.11644 https://doi.org/10.48550/arXiv.2403.11644 (2024)

  39. [39]

    T. N. Georges, B. K. Berntson, C. S¨ underhauf, and A. V. Ivanov, Pauli decomposition via the fast Walsh- Hadamard transform, New Journal of Physics27, 033004 (2025)

  40. [40]

    Aguilar, S

    G. Aguilar, S. Cichy, J. Eisert, and L. Bittel, Full classification of Pauli Lie algebras, arXiv:2408.00081 https://doi.org/10.48550/arXiv.2408.00081 (2024)

  41. [41]

    O’Donnell,Analysis of Boolean Functions(Cambridge University Press, 2014)

    R. O’Donnell,Analysis of Boolean Functions(Cambridge University Press, 2014)

  42. [42]

    Kushilevitz and Y

    E. Kushilevitz and Y. Mansour, Learning decision trees using the fourier spectrum, SIAM Journal of Computing 22, https://doi.org/10.1137/0222080 (1993)

  43. [43]

    Scheibler, S

    R. Scheibler, S. Haghighatshoar, and M. Vetterli, A fast Hadamard transform for signals with sublinear sparsity in the transform domain, IEEE Transactions on Informa- tion Theory61, 2115 (2015)

  44. [44]

    false positives

    R. Meshulam, An uncertainty inequality for finite abelian groups, European Journal of Combinatorics27, 63 (2006). 16 Appendix A: Correctness and error analysis In this Appendix, we show that our algorithm correctly recovers the Pauli decomposition with high probability. First, recall that for a fixedx∈ {0,1} n, we define bx(u) :=M x⊕u,u (A1) for a query o...

  45. [45]

    Ifp x = 1(i.e.,xis unique), the algorithm returns(Unique,ˆz, β)withˆz=z ⋆ andβ=β x(z⋆)deterministically, wherez ⋆ is the unique support element ofβ x

  46. [46]

    Proof.First, assumep x = 1

    Ifp x ≥2, the algorithm returnsDegeneratewith probability at least1−η x. Proof.First, assumep x = 1. Then, there exists a uniquez ⋆ ∈ {0,1} n such thatβ x(z⋆)̸= 0, or, in other words, bx(u) =β x(z⋆)(−1)z⋆·u (A14) for allu∈ {0,1} n. If we takeu= 0 n, then we have bx(0n) =β x(z⋆)̸= 0.(A15) Then, for each standard basis vectore j, we have bx(ej) bx(0n) = βx(...

  47. [47]

    Some true residual support element is never isolated in any of the folding rounds

  48. [48]

    Some collision bin is incorrectly accepted as a singleton

  49. [49]

    19 We allocate failure probability at mostδ x/3 to each of these events for simplicity

    After the peeling rounds, the residual is nonzero but the final randomized residual check fails to detect it. 19 We allocate failure probability at mostδ x/3 to each of these events for simplicity. If a binscontains exactly one support elementz ⋆, theng(s) =β x(z⋆) andg (w)(s) =β x(z⋆)(−1)z⋆·w for all shifts w. Choosingw=e j forj∈[n] and taking the ratio ...