Quantum algorithm for anisotropic diffusion and convection equations with vector norm scaling

Abstract

In this work, we tackle the resolution of partial differential equations (PDEs) on digital quantum computers. Two fundamental PDEs are addressed: the anisotropic diffusion equation and the anisotropic convection equation. We present a quantum numerical scheme consisting of three steps: quantum state preparation, evolution with diagonal operators, and measurement of observables of interest. The evolution step relies on a high-order centered finite difference and a product formula approximation, also known as Trotterization. We provide novel vector-norm analysis to bound the different sources of error. We prove that the number of time-steps required in the evolution can be reduced by a factor $\Theta (16^n)$ for the diffusion equation, and $\Theta (4^n)$ for the convection equation, where $n$ is the number of qubits per dimension, an exponential reduction compared to the previously established operator-norm analysis.

This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this contribution is published in Quantum Engineering Sciences and Technologies for Industry and Services (QUEST-IS 2025), and is available online at https://doi.org/10.1007/978-3-032-13855-2_23

References in corpus (15)

• Hamiltonian Simulation by Qubitization
• Optimal Hamiltonian Simulation by Quantum Signal Processing
• Polynomial-time quantum algorithm for the simulation of chemical dynamics
• Iterative Quantum Amplitude Estimation
• Quantum State Preparation with Optimal Circuit Depth: Implementations and Applications
• Efficient Quantum Circuits for Diagonal Unitaries Without Ancillas
• Quantum Simulation of Chemistry with Sublinear Scaling in Basis Size
• Exponentially More Precise Quantum Simulation of Fermions in the Configuration Interaction Representation
• Low depth algorithms for quantum amplitude estimation
• Time-dependent unbounded Hamiltonian simulation with vector norm scaling
• Bounding the costs of quantum simulation of many-body physics in real space
• Optimal (controlled) quantum state preparation and improved unitary synthesis by quantum circuits with any number of ancillary qubits
• Amplitude Estimation from Quantum Signal Processing
• Linear-depth quantum circuits for loading Fourier approximations of arbitrary functions
• Quantum state preparation for multivariate functions