Fourier extensions for matrix-function block encodings with error-independent subnormalization bounds

Abstract:Block encodings of non-unitary matrix functions are central to quantum numerical linear algebra. Hamiltonian simulation is a natural input model for Hermitian matrices, but accurate block encodings often incur large subnormalization. We decouple the accuracy from the subnormalization by formulating the matrix-function block encoding as a Fourier-extension approximation problem, yielding a linear combination of unitaries for Hermitian matrix inputs. Fourier extensions approximate non-periodic functions by a Fourier series on a larger periodic domain, creating redundant coefficients that can be optimized for their absolute sum, and hence subnormalization. The coefficients may be chosen for optimal subnormalization with algebraic convergence, by tuning the subnormalization bound for increasing rates of exponential convergence, or by Sobolev-regularized fitting to accommodate more general spectral sets. Fourier-extension block encodings apply to eigenvalue transforms of Hermitian matrices, or to odd singular-value transforms of general matrices, including as a quantum linear systems algorithm.

arXiv:2601.18024 [quant-ph]

(or arXiv:2601.18024v2 [quant-ph] for this version)

https://doi.org/10.48550/arXiv.2601.18024 Focus to learn more arXiv-issued DOI via DataCite

arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.

Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.

Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.