Adding a fast-forwardable penalty projection to the generator of a quantum-simulated ODE enforces boundary conditions up to error ε, with gate complexity overhead O(log λ).
Compact quantum algorithms for time-dependent differential equations
1 Pith paper cite this work. Polarity classification is still indexing.
abstract
Many claims of computational advantages have been made for quantum computing over classical, but they have not been demonstrated for practical problems. Here, we present algorithms for solving time-dependent PDEs, with particular reference to fluid equations. We build on an idea based on linear combination of unitaries to simulate non-unitary, non-Hermitian quantum systems, and generate hybrid quantum-classical algorithms that efficiently perform iterative matrix-vector multiplication and matrix inversion operations. These algorithms are end-to-end, with relatively low-depth quantum circuits that demonstrate quantum advantage, with the best-case asymptotic complexities, which we show are near-optimal. We demonstrate the performance of the algorithms by conducting: (a) fully gate level, state-vector simulations using an in-house, high performance, quantum simulator called QFlowS; (b) experiments on a real quantum device; and (c) noisy simulations using Qiskit Aer. We also provide device specifications such as error-rates (noise) and state sampling (measurement) to accurately perform convergent flow simulations on noisy devices. The results offer evidence that the proposed algorithm is amenable for use on near-term quantum devices.
citation-role summary
citation-polarity summary
fields
quant-ph 1years
2025 1verdicts
CONDITIONAL 1roles
background 1polarities
unclear 1representative citing papers
citing papers explorer
-
Arbitrary Boundary Conditions and Constraints in Quantum Algorithms for Differential Equations via Penalty Projections
Adding a fast-forwardable penalty projection to the generator of a quantum-simulated ODE enforces boundary conditions up to error ε, with gate complexity overhead O(log λ).