Quantum Signal Processing Phase Factorization via Logarithmic Derivatives
Abstract
Quantum signal processing (QSP) requires the efficient computation of phase factors associated with a polynomial transformation. One existing stable factorization method avoids direct root finding but relies on Prony's method to recover a polynomial from weighted power sums of its roots (arXiv:2202.02671). We present a modification of this approach that eliminates Prony's method by using the logarithmic derivative of an associated polynomial. Because the logarithmic derivative has the same residue at every root after accounting for multiplicity, its Fourier coefficients directly yield the unweighted power sums of the roots inside the unit disk. The Newton--Girard identities then recover the coefficients of the corresponding polynomial without computing its roots. The resulting algorithm has O(dlog2 d) computational complexity and simplifies the factorization procedure of the original method. In numerical experiments on Hamiltonian simulation, eigenstate filtering, matrix inversion, and the Fermi--Dirac operator, our implementation computes phase factors in seconds with maximum errors on the order of 10−12.
Related
- Particle-Number-Preserving Quantum Signal Processing for Thermal Properties: Canonical-Ensemble Estimation and Grand-Canonical Reconstruction
- Representation Learning with Quantum Signal Processing
- Constrained minimax approximation for quantum signal processing
- Automated Verification of Hamiltonian Simulation by Quantum Signal Processing
- Robustness of Quantum Signal Processing and Resource-aware Polynomial Design
BibTeX
@misc{peng2026quantum,
title = {{Quantum Signal Processing Phase Factorization via Logarithmic Derivatives}},
author = {Christopher Peng},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.14259},
eprint = {2610.14259},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.14259},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Quantum Signal Processing Phase Factorization via Logarithmic Derivatives
AU - Christopher Peng
PY - 2026
DA - 2026/10/09
PB - arXiv
DO - 10.48550/arxiv.2610.14259
UR - https://arxiv.org/abs/2610.14259
AB - Quantum signal processing (QSP) requires the efficient computation of phase factors associated with a polynomial transformation. One existing stable factorization method avoids direct root finding but relies on Prony's method to recover a polynomial from weighted power sums of its roots (arXiv:2202.02671). We present a modification of this approach that eliminates Prony's method by using the logarithmic derivative of an associated polynomial. Because the logarithmic derivative has the same residue at every root after accounting for multiplicity, its Fourier coefficients directly yield the unweighted power sums of the roots inside the unit disk. The Newton--Girard identities then recover the coefficients of the corresponding polynomial without computing its roots. The resulting algorithm has O(dlog^2 d) computational complexity and simplifies the factorization procedure of the original method. In numerical experiments on Hamiltonian simulation, eigenstate filtering, matrix inversion, and the Fermi--Dirac operator, our implementation computes phase factors in seconds with maximum errors on the order of 10^−12.
ER -