Quantum Eigenvalue Transformations for Arbitrary Matrices
Pith reviewed 2026-05-10 03:02 UTC · model grok-4.3
The pith
Special block encodings let quantum signal processing apply polynomials to eigenvalues of any matrix.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that an n-regular block encoding of an arbitrary square matrix A turns the application of quantum signal processing into the direct application of a polynomial of degree at most n to A, with the transformation acting on the eigenvalues associated with the Jordan normal form of A and independent of any further details of A's structure.
What carries the argument
The n-regular block encoding: a unitary whose successive powers up to order n reproduce the corresponding powers of the block-encoded matrix.
If this is right
- Any polynomial of degree at most n can be applied to the eigenvalues of an arbitrary square matrix through quantum signal processing.
- The transformation works without requiring the matrix to be diagonalizable or any explicit knowledge of its Jordan structure.
- Standard block encodings can be upgraded to n-regular form using only O(log n) ancillary qubits and gates.
- The resulting operator equals the polynomial applied to the original matrix regardless of the matrix's internal details.
Where Pith is reading between the lines
- This approach could be combined with existing quantum linear-system solvers to handle non-Hermitian operators directly.
- The same regularity idea might extend to rectangular matrices or to functions beyond polynomials.
- Small quantum simulations on two-by-two Jordan blocks would provide an immediate experimental test of the eigenvalue mapping.
Load-bearing premise
That an n-regular block encoding can be built from a standard one without adding errors that would stop the polynomial from correctly transforming the eigenvalues of non-diagonalizable matrices.
What would settle it
A direct computation or small-scale circuit run on a non-diagonalizable Jordan block matrix showing that the output state after the quantum signal processing step fails to match the eigenvalues obtained by applying the target polynomial to the Jordan form would falsify the claim.
Figures
read the original abstract
Quantum Signal Processing (QSP) and Quantum Singular Value Transformation (QSVT) provide an efficient framework for implementing polynomials of block-encoded matrices, and thus offer a systematic approach to quantum algorithm design. However, despite a number of recent advances, important limitations remain. In particular, QSP can only transform unitary matrices, by applying a polynomial to their eigenvalues, while QSVT is a singular-value transformation and thus one can only obtain the polynomial of Hermitian matrices. As a consequence, these techniques do not directly apply to an arbitrary non-Hermitian matrix that is not diagonalizable. In this work, we propose a simple yet powerful method to extend these ideas to arbitrary square matrices by acting on their eigenvalues. To this end, we introduce the notion of an $n$-regular block encoding, namely, a block encoding whose $k$-th power reproduces the $k$-th power of the encoded matrix for every $0 < k < n$. We show that applying QSP to any unitary with this property is equivalent to applying a polynomial of degree at most $n$ to the block-encoded matrix, independently of its internal structure. Moreover, we provide a simple construction that transforms any block encoding into an $n$-regular one using only $O(\log n)$ ancillary qubits and operations. Finally, we show that this construction induces the desired transformation on the eigenvalues associated with the Jordan normal form of the matrix.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces n-regular block encodings and shows that QSP applied to any unitary with this property implements a polynomial of degree at most n on the block-encoded matrix A, independently of its internal structure. It provides an explicit construction converting any block encoding into an n-regular one with O(log n) ancillary qubits and operations, and claims this induces the desired eigenvalue transformation on the Jordan normal form of arbitrary (including defective) matrices.
Significance. If the central construction and equivalence hold exactly, the result would meaningfully extend QSP/QSVT beyond unitary and Hermitian cases, enabling systematic polynomial transformations on the eigenvalues of general square matrices via Jordan calculus. The O(log n) overhead and structure-independent claim are efficient and broadly applicable strengths; the work supplies a concrete algorithmic primitive that could support new quantum linear-algebra routines.
major comments (2)
- [§3] §3 (Construction of the n-regular block encoding): the argument that the O(log n)-ancilla unitary U satisfies [U^k]_{top-left} = A^k exactly for all k ≤ n must be verified for defective matrices. For a Jordan block J = D + N with nilpotent N ≠ 0, any ancillary leakage into the generalized eigenspace would alter the action of higher powers on N and therefore change p(A) away from the intended Jordan-functional-calculus result. An explicit invariance lemma or direct computation on the generalized eigenvectors is required.
- [Theorem 4.1] Theorem 4.1 (Equivalence of QSP on n-regular unitary to polynomial action on A): the proof that the top-left block of the QSP circuit equals p(A) relies on the n-regular property holding without error. If the construction in §3 fails to preserve exact powers for non-diagonalizable A, the claimed independence from internal structure does not follow. The manuscript should supply a short error-bound or commutator calculation showing that the ancillary registers do not mix the Jordan chains.
minor comments (2)
- [Definition 2.3] The definition of n-regular block encoding (Eq. (3) or equivalent) would be clearer if written as an explicit block-matrix condition on U^k rather than a verbal statement.
- [Figure 2] Figure 2 (circuit for the n-regular construction) would benefit from explicit labeling of the O(log n) ancilla registers and a brief caption explaining how the controlled operations enforce the power-reproduction property.
Simulated Author's Rebuttal
We thank the referee for their careful reading of our manuscript and for providing constructive feedback. We are pleased that the referee recognizes the potential significance of extending QSP to arbitrary matrices. We address each of the major comments below and outline the revisions we plan to make.
read point-by-point responses
-
Referee: [§3] §3 (Construction of the n-regular block encoding): the argument that the O(log n)-ancilla unitary U satisfies [U^k]_{top-left} = A^k exactly for all k ≤ n must be verified for defective matrices. For a Jordan block J = D + N with nilpotent N ≠ 0, any ancillary leakage into the generalized eigenspace would alter the action of higher powers on N and therefore change p(A) away from the intended Jordan-functional-calculus result. An explicit invariance lemma or direct computation on the generalized eigenvectors is required.
Authors: We appreciate the referee's emphasis on rigorously verifying the construction for defective matrices. Our construction of the n-regular block encoding is designed such that the ancillary qubits are entangled in a controlled manner that ensures the top-left block exactly reproduces A^k for k ≤ n, without introducing errors from generalized eigenspaces. This holds because the additional operations are applied in a way that commutes with the Jordan structure in the relevant subspace. To make this explicit, we will add a new lemma in §3 that provides a direct computation on generalized eigenvectors, demonstrating invariance of the Jordan chains under the powered unitary. This will confirm that no ancillary leakage alters the nilpotent part. revision: yes
-
Referee: [Theorem 4.1] Theorem 4.1 (Equivalence of QSP on n-regular unitary to polynomial action on A): the proof that the top-left block of the QSP circuit equals p(A) relies on the n-regular property holding without error. If the construction in §3 fails to preserve exact powers for non-diagonalizable A, the claimed independence from internal structure does not follow. The manuscript should supply a short error-bound or commutator calculation showing that the ancillary registers do not mix the Jordan chains.
Authors: We agree that strengthening the proof of Theorem 4.1 with an explicit calculation would enhance clarity. The proof proceeds by showing that each step in the QSP sequence preserves the n-regular property, leading to the top-left block being p(A) by the definition of polynomial application via the Jordan functional calculus. We will incorporate a short commutator calculation in the revised manuscript to explicitly show that the ancillary registers do not mix the Jordan chains, thereby confirming the structure-independent nature of the transformation. revision: yes
Circularity Check
No circularity; derivation relies on explicit new definition and construction
full rationale
The paper defines n-regular block encoding directly via the power-reproduction property on blocks, supplies an O(log n)-ancilla construction to realize it from any standard block encoding, and derives the QSP-polynomial equivalence from that definition plus the algebraic action of QSP on the unitary. The Jordan-functional-calculus claim follows from the preserved block powers rather than from any self-citation chain, fitted parameter, or renamed prior result. No load-bearing step reduces to its own input by construction; the central claims are supported by the provided construction and are therefore self-contained.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Standard block-encoding properties and QSP polynomial implementation hold for the unitary part of the n-regular encoding.
invented entities (1)
-
n-regular block encoding
no independent evidence
Reference graph
Works this paper leans on
-
[1]
G. H. Low and I. L. Chuang, Optimal Hamiltonian Sim- ulation by Quantum Signal Processing, Physical Review Letters118, 010501 (2017)
2017
-
[2]
J. M. Martyn, Z. M. Rossi, A. K. Tan, and I. L. Chuang, A Grand Unification of Quantum Algorithms, PRX Quantum2, 40203 (2021)
2021
-
[3]
G. H. Low and I. L. Chuang, Hamiltonian Simulation by Qubitization, Quantum3, 163 (2019)
2019
-
[4]
A. M. Childs, R. Kothari, and R. D. Somma, Quan- tum Algorithm for Systems of Linear Equations with Exponentially Improved Dependence on Precision, SIAM Journal on Computing46, 1920 (2017)
1920
-
[5]
Lin and Y
L. Lin and Y. Tong, Near-optimal ground state prepara- tion, Quantum4, 372 (2020)
2020
-
[6]
Y. Dong, L. Lin, and Y. Tong, Ground-State Preparation and Energy Estimation on Early Fault-Tolerant Quan- tum Computers via Quantum Eigenvalue Transformation of Unitary Matrices, PRX Quantum3, 40305 (2022)
2022
-
[7]
S. McArdle, A. Gily´ en, and M. Berta, Quantum state preparation without coherent arithmetic (2022), arXiv:2210.14892 [quant-ph]
-
[8]
Robust black-box quantum-state preparation via quantum signal processing
L. Laneve, Robust black-box quantum-state preparation via quantum signal processing (2023), arXiv:2305.04705 [quant-ph]
-
[9]
G. H. Low, T. J. Yoder, and I. L. Chuang, Methodology of Resonant Equiangular Composite Quantum Gates, Phys- ical Review X6, 41067 (2016)
2016
-
[10]
Gily´ en, Y
A. Gily´ en, Y. Su, G. H. Low, and N. Wiebe, Quantum singular value transformation and beyond: Exponential improvements for quantum matrix arithmetics, inPro- ceedings of the 51st Annual ACM SIGACT Symposium on Theory of Computing(ACM, 2019) pp. 193–204
2019
-
[11]
Motlagh and N
D. Motlagh and N. Wiebe, Generalized Quantum Signal Processing, PRX Quantum5, 020368 (2024)
2024
-
[12]
Alexis, G
M. Alexis, G. Mnatsakanyan, and C. Thiele, Quantum signal processing and nonlinear Fourier analysis, Revista Matem´ atica Complutense37, 655 (2024)
2024
-
[13]
Alexis, L
M. Alexis, L. Lin, G. Mnatsakanyan, C. Thiele, and J. Wang, Infinite quantum signal processing for arbitrary Szeg¨ o functions, Communications on Pure and Applied Mathematics79, 123 (2026)
2026
-
[14]
Generalized quantum signal processing and non- linear fourier transform are equivalent,
L. Laneve, Generalized Quantum Signal Processing and Non-Linear Fourier Transform are equivalent (2025), arXiv:2503.03026 [quant-ph]
- [15]
-
[16]
arXiv:2312.00723 [quant-ph] doi:10
C. S¨ underhauf, Generalized Quantum Singular Value Transformation (2023), arXiv:2312.00723 [quant-ph]
-
[17]
G. H. Low and Y. Su, Quantum Eigenvalue Processing, in2024 IEEE 65th Annual Symposium on Foundations of Computer Science (FOCS)(2024) pp. 1051–1062
2024
-
[18]
An, J.-P
D. An, J.-P. Liu, and L. Lin, Linear Combination of Hamiltonian Simulation for Nonunitary Dynamics with Optimal State Preparation Cost, Physical Review Letters 131, 150603 (2023)
2023
-
[19]
https://arxiv.org/abs/2508.19596
X. Huang and D. An, Fourier transform-based lin- ear combination of Hamiltonian simulation (2025), arXiv:2508.19596 [quant-ph]
- [20]
-
[21]
D. An, A. M. Childs, L. Lin, and L. Ying, Laplace Transform–Based Quantum Eigenvalue Transformation via Linear Combination of Hamiltonian Simulation, SIAM Journal on Computing55, 376 (2026)
2026
-
[22]
A. M. Childs and N. Wiebe, Hamiltonian simulation us- ing linear combinations of unitary operations, Quantum Information and Computation12, 901 (2012)
2012
-
[23]
J.-P. Liu, H. O. Kolden, H. K. Krovi, N. F. Loureiro, K. Trivisa, and A. M. Childs, Efficient quantum al- gorithm for dissipative nonlinear differential equations, Proceedings of the National Academy of Sciences118, e2026805118 (2021)
2021
-
[24]
Krovi, Improved quantum algorithms for linear and nonlinear differential equations, Quantum7, 913 (2023)
H. Krovi, Improved quantum algorithms for linear and nonlinear differential equations, Quantum7, 913 (2023)
2023
-
[25]
A. M. Childs, J.-P. Liu, and A. Ostrander, High-precision quantum algorithms for partial differential equations, Quantum5, 574 (2021)
2021
-
[26]
Weimer, A
H. Weimer, A. Kshetrimayum, and R. Or´ us, Simulation methods for open quantum many-body systems, Reviews of Modern Physics93, 015008 (2021)
2021
- [27]
-
[28]
A. M. Dalzell, S. McArdle, M. Berta, P. Bienias, C.-F. Chen, A. Gily´ en, C. T. Hann, M. J. Kastoryano, E. T. Khabiboulline, A. Kubica, G. Salton, S. Wang, and F. G. S. L. Brand˜ ao,Quantum Algorithms: A Survey of Appli- cations and End-to-end Complexities(Cambridge Uni- versity Press, Cambridge, 2025)
2025
-
[29]
Takahira, A
S. Takahira, A. Ohashi, T. Sogabe, and T. S. Usuda, Quantum Algorithms based on the Block-Encoding Framework for Matrix Functions by Contour Integrals, Quantum Information and Computation22, 965 (2022)
2022
-
[30]
S. Jiang and D. An, Contour-integral based quantum eigenvalue transformation: Analysis and applications (2026), arXiv:2601.11959 [quant-ph]
-
[31]
Jetzek,Galois Fields, Linear Feedback Shift Registers and Their Applications(Hanser, 2019)
U. Jetzek,Galois Fields, Linear Feedback Shift Registers and Their Applications(Hanser, 2019). End Matter Proof of Lemma 4.For a given polynomial of degreen P(z) = Pn k=0 αkzk, ifUis block-encoding ˜Aexactly, then by linearity and the definition ofn-regular encoding it block-encodesP( ˜A) exactly. The error in the block- encoding is bounded by P( ˜A)−P(A)...
2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.