Posted

Hongkang Ni, Lexing Ying (Jun 30 2026).
Abstract: Linear combination of Hamiltonian simulation (LCHS) provides an efficient method for implementing matrix exponentials etAe^{-tA} on quantum computers. In this paper, we develop LCHS formulas for computing general matrix functions f(A)f(A) when ff is analytic on the numerical range of AA, with AA possibly non-normal. The essential technical tool is Weyl calculus, which reduces the construction of LCHS formulas for noncommuting operators to scalar Fourier approximation problems. Our construction yields a quantum eigenvalue transformation algorithm with optimal O(log1ϵ)\mathcal{O}(\log\frac{1}{\epsilon}) query complexity scaling. Furthermore, our Weyl-calculus-based theory gives rise to an ansatz-free convex optimization framework that directly produces discrete LCHS formulas. This circumvents the inefficiencies of traditional quadrature rules and yields formulas highly optimized for coherent implementation on quantum computers. In addition, both our theory and optimization framework apply to the simulation of time-dependent dissipative ODE ddtψ(t)=A(t)ψ(t)\frac{\mathrm{d}}{\mathrm{d} t} \psi(t) = -A(t)\psi(t), for which we achieve a 2.1×2.1\times cost reduction over prior art.

Order by:

Want to join this discussion?

Join our community today and start discussing with our members by participating in exciting events, competitions, and challenges. Sign up now to engage with quantum experts!