Mathematical and numerical analysis of quantum signal processing
Lin Lin
TL;DR
This survey synthesizes the mathematical and numerical analysis of quantum signal processing (QSP), clarifying when phase factors exist and are unique, and unveiling deep connections to nonlinear Fourier analysis on SU(2). It covers infinite QSP (iQSP) and convergence (L1 and L2) via Szegő theory, complementary polynomials via the Weiss algorithm, and robust inverse NLFT algorithms (layer stripping, Riemann–Hilbert factorization, inverse nonlinear FFT) with stability analyses. It also surveys iterative phase-factor algorithms, the pivotal role of QSP in quantum singular value transformation (QSVT), and practical applications to Hamiltonian simulation, linear systems, and eigenvalue problems, as well as multiple generalizations and open directions. The work highlights how displacement-structure insights, outer function choices, and NLFT theory jointly enable efficient, stable representations and computational pipelines for polynomial and non-polynomial function mappings on quantum hardware.
Abstract
Quantum signal processing (QSP) provides a representation of scalar polynomials of degree $d$ as products of matrices in $\mathrm{SU}(2)$, parameterized by $(d+1)$ real numbers known as phase factors. QSP is the mathematical foundation of quantum singular value transformation (QSVT), which is often regarded as one of the most important quantum algorithms of the past decade, with a wide range of applications in scientific computing, from Hamiltonian simulation to solving linear systems of equations and eigenvalue problems. In this article we survey recent advances in the mathematical and numerical analysis of QSP. In particular, we focus on its generalization beyond polynomials, the computational complexity of algorithms for phase factor evaluation, and the numerical stability of such algorithms. The resolution to some of these problems relies on an unexpected interplay between QSP, nonlinear Fourier analysis on $\mathrm{SU}(2)$, fast polynomial multiplications, and Gaussian elimination for matrices with displacement structure.
