In this paper, we present two Hamiltonian simulation algorithms for multiscale linear transport equations, combining the Schrödingerization method [S. Jin, N. Liu and Y. Yu, Phys. Rev. Lett, 133(2024), 230602] [S. Jin, N. Liu and Y. Yu, Phys. Rev. A, 108 (2023), 032603] and exponential integrator while incorporating incoming boundary conditions. These two algorithms each have advantages in terms of design easiness and scalability, and the query complexity of both algorithms, outperforms classical algorithms for solving this equation. In terms of the theoretical framework, these are the first quantum Hamiltonian simulation algorithms for multiscale linear transport equation to combine the Schrödingerization method with effective asymptotic-preserving schemes, which are efficient for handling multiscale problems with stiff terms.
In this paper, we present two Hamiltonian simulation algorithms for multiscale linear transport equations, combining the Schrödingerization method [S. Jin, N. Liu and Y. Yu, Phys. Rev. Lett, 133(2024), 230602] [S. Jin, N. Liu and Y. Yu, Phys. Rev. A, 108 (2023), 032603] and exponential integrator while incorporating incoming boundary conditions. These two algorithms each have advantages in terms of design easiness and scalability, and the query complexity of both algorithms, outperforms classical algorithms for solving this equation. In terms of the theoretical framework, these are the first quantum Hamiltonian simulation algorithms for multiscale linear transport equation to combine the Schrödingerization method with effective asymptotic-preserving schemes, which are efficient for handling multiscale problems with stiff terms.
Abstract. Quantum dynamics, typically expressed in the form of a time-dependent Schrödinger equation with a Hermitian Hamiltonian, is a natural application for quantum computing. However, when simulating quantum dynamics that involves the emission of electrons, it is necessary to use artificial boundary conditions (ABCs) to confine the computation within a fixed domain. The introduction of ABCs alters the Hamiltonian structure of the dynamics, and existing quantum algorithms cannot be directly applied since the evolution is no longer unitary. The current paper utilizes a recently introduced Schrödingerization method that converts non-Hermitian dynamics into a Schrödinger form for the artificial boundary problems [S. Jin, N. Liu, and Y. Yu, Quantum Simulation of Partial Differential Equations via Schrödingerisation, preprint, arXiv:2212.13969 , 2022], [S. Jin, N. Liu, and Y. Yu, Phys. Rev. A, 108 (2023), 032603]. We implement this method for three types of ABCs, including the complex absorbing potential technique, perfectly matched layer methods, and Dirichlet-to-Neumann approach. We analyze the query complexity of these algorithms and perform numerical experiments to demonstrate the validity of this approach. This helps to bridge the gap between available quantum algorithms and computational models for quantum dynamics in unbounded domains.
Abstract. We analyze the Schrödingerization method for quantum simulation of a general class of nonunitary dynamics with inhomogeneous source terms. The Schrödingerization technique, introduced in [S. Jin, N. Liu, and Y. Yu, Phys. Rev. Lett., 133 (2024), 230602], transforms any linear ordinary and partial differential equations with nonunitary dynamics into a system under unitary dynamics via a warped phase transition that maps the equations into a higher dimension, making them suitable for quantum simulation. This technique can also be applied to these equations with inhomogeneous terms modeling source or forcing terms, or boundary and interface conditions, and discrete dynamical systems such as iterative methods in numerical linear algebra, through extra equations in the system. Difficulty arises with the presence of inhomogeneous terms since they can change the stability of the original system. In this paper, we systematically study—both theoretically and numerically—the important issue of recovering the original variables from the Schrödingerized equations, even when the evolution operator contains unstable modes. We show that even with unstable modes, one can still construct a stable scheme; however, to recover the original variable, one needs to use suitable data in the extended space. We analyze and compare both the discrete and continuous Fourier transforms used in the extended dimension and derive corresponding error estimates, which allow one to use the more appropriate transform for specific equations. We also provide a smoother initialization for the Schrödingerized system to gain higher-order accuracy in the extended space. We homogenize the inhomogeneous terms with a stretch transformation, making it easier to recover the original variable. Our recovery technique also provides a simple and generic framework to solve general ill-posed problems in a computationally stable way.
We investigate time complexities of finite difference methods for solving the multiscale transport equation with quantum algorithms. We find that the time complexities of both the classical treatment and quantum treatment for a standard explicit scheme scale as O(1/ε), where ε is the small scaling parameter, while the complexities for the even-odd parity based Asymptotic-Preserving (AP) scheme do not depend on ε. This indicates that it is still of great importance to use AP (and probably other efficient multiscale) schemes for multiscale problems in quantum computing when solving multiscale transport or kinetic equations.
Quantum linear systems algorithms (QLSAs) aim to solve linear systems $A\bb{x}=\bb{b}$ exponentially faster than classical methods under certain conditions. In this work, we develop quantum algorithms for solving linear algebraic equations from an ODE-based perspective. Inspired by the linear combination of Hamiltonian simulation (LCHS) representation in the Fourier approach \cite{Childs2017QLSA}, we express the solution $\bb{x}$ as a linear combination of solutions to a system of linear convection equations, which become Schrödinger-type equations with unitary evolutions in the Fourier domain. We refer to this representation as LC-Schrödingerization. Based on this result, we construct an LCHS-based quantum algorithm with two LCHS instances: one for time-marching and one for numerical integration. The key construction uses the derivative of a Gaussian-smoothed hat function and recovers the solution over a fixed auxiliary interval. This permits a truncation time independent of the target accuracy and avoids the loss in success probability from selecting a single grid point. Periodization and explicit Fourier coefficients provide the corresponding projection error bounds. Under the stated oracle assumptions and given a constant-factor estimate of the solution norm, direct simulation of the select operators and block preconditioning achieve the optimal matrix-query complexity O(κAlogε1) without using variable-time amplitude amplification (VTAA). The same upper bound holds for queries to the right-hand-side preparation oracle.
In this paper, we present a quantum implicit-explicit (IMEX) scheme for multiscale ordinary and partial differential equations whose discretization parameters are independent of the scaling parameter ε. A key ingredient of our approach is a continuous-time formulation of classical IMEX schemes, which decouples the evolution time of the quantum algorithm from the physical time of the differential equation and is therefore particularly useful in multiscale settings. Building on this idea, we employ the Schrödingerization framework [Phys. Rev. Lett. 133 (2024), 230602] to implement IMEX schemes on quantum computers. Compared to previous HHL type quantum AP scheme [J. Comput. Phys. 471 (2022), 111641], this new method requires narrower -- an extra logarithmic factor -- auxiliary register numerical examples on linear heat and multiscale telegraph equations demonstrate the independence in ε of the method.
Abstract. We introduce a simple and stable computational method for ill-posed partial differential equation (PDE) problems. The method is based on Schrödingerization, introduced in [S. Jin, N. Liu, and Y. Yu, Phys. Rev. Lett., 133 (2024), 230602; S. Jin, N. Liu, and Y. Yu, Phys. Rev. A, 108 (2023), 032603], which maps all linear PDEs into Schrödinger-type equations in one higher dimension, for quantum simulations of these PDEs. Although the original problem is ill-posed, the Schrödingerized equations are Hamiltonian systems and time-reversible, allowing stable computation both forward and backward in time. The original variable can be recovered by data from suitably chosen domain in the extended dimension. We will use the (constant and variable coefficient) backward heat equation and the linear convection equation with imaginary wave speed as examples. Error analysis of these algorithms are conducted and verified numerically. The methods are applicable to both classical and quantum computers, and we also lay out quantum algorithms for these methods. Moreover, we introduce a smooth initialization for the Schrödingerized equation, which will lead to essentially spectral accuracy for the approximation in the extended space, if a spectral method is used. Consequently, the extra qubits needed due to the extra dimension, if a qubit-based quantum algorithm is used, for both well-posed and ill-posed problems, becomes almost [Formula: see text] where [Formula: see text] is the desired precision. This optimizes the complexity of the Schrödingerization based quantum algorithms for any nonunitary dynamical system introduced in [S. Jin, N. Liu, and Y. Yu, Phys. Rev. Lett., 133 (2024), 230602; S. Jin, N. Liu, and Y. Yu, Phys. Rev. A, 108 (2023), 032603].