We give detailed analysis and circuit design of structure-preserving quantum algorithms for second-order linear evolutionary PDEs, including parabolic equations and hyperbolic equations with mixed Dirichlet, Neumann, and periodic boundary conditions and source terms. While prior quantum algorithms usually neglect the stability problem from the PDE-to-ODE reduction, our method-of-lines approach investigates the boundary lifting via Coons interpolation and boundary-aware discretization, so that the resulting semi-discrete systems are stable and compatible with efficient quantum ODE primitives. For the parabolic problem, we use a diagonal similarity transform to ensure the semi-discrete generator must have a positive semi-definite Hermitian part, and then solve the resulting ODE system by the optimal linear combination of Hamiltonian simulation (LCHS). For the hyperbolic problem, we rewrite the semi-discrete equation as an equivalent first-order system and solve it by Hamiltonian simulation. We implement our quantum algorithms with explicit block-encoding constructions and circuit implementations, as well as demonstrating the end-to-end complexity bounds together with spatial and quadrature error estimates. We conduct classical numerical experiments on the convection-diffusion equation, inhomogeneous heat equation, and Klein-Gordon equation to validate our structure-preserving analysis and algorithmic constructions.
We give detailed analysis and circuit design of structure-preserving quantum algorithms for second-order linear evolutionary PDEs, including parabolic equations and hyperbolic equations with mixed Dirichlet, Neumann, and periodic boundary conditions and source terms. While prior quantum algorithms usually neglect the stability problem from the PDE-to-ODE reduction, our method-of-lines approach investigates the boundary lifting via Coons interpolation and boundary-aware discretization, so that the resulting semi-discrete systems are stable and compatible with efficient quantum ODE primitives. For the parabolic problem, we use a diagonal similarity transform to ensure the semi-discrete generator must have a positive semi-definite Hermitian part, and then solve the resulting ODE system by the optimal linear combination of Hamiltonian simulation (LCHS). For the hyperbolic problem, we rewrite the semi-discrete equation as an equivalent first-order system and solve it by Hamiltonian simulation. We implement our quantum algorithms with explicit block-encoding constructions and circuit implementations, as well as demonstrating the end-to-end complexity bounds together with spatial and quadrature error estimates. We conduct classical numerical experiments on the convection-diffusion equation, inhomogeneous heat equation, and Klein-Gordon equation to validate our structure-preserving analysis and algorithmic constructions.
This paper explores the feasibility of quantum simulation for partial differential equations (PDEs) with physical boundary or interface conditions. Semi-discretisation of such problems does not necessarily yield Hamiltonian dynamics and even alters the Hamiltonian structure of the dynamics when boundary and interface conditions are included. This seemingly intractable issue can be resolved by using a recently introduced Schrödingerisation method (Jin et al. 2022) -- it converts any linear PDEs and ODEs with non-Hermitian dynamics to a system of Schrödinger equations, via the so-called warped phase transformation that maps the equation into one higher dimension. We implement this method for several typical problems, including the linear convection equation with inflow boundary conditions and the heat equation with Dirichlet and Neumann boundary conditions. For interface problems, we study the (parabolic) Stefan problem, linear convection, and linear Liouville equations with discontinuous and even measure-valued coefficients. We perform numerical experiments to demonstrate the validity of this approach, which helps to bridge the gap between available quantum algorithms and computational models for classical and quantum dynamics with boundary and interface conditions.
Abstract We propose an explicit quantum framework for numerically simulating general linear partial differential equations (PDEs), extending previous work (Guseynov et al 2025 Phys. Rev. Res. 7 033100) to incorporate (a) Robin boundary conditions—which include Neumann and Dirichlet conditions as special cases–(b) inhomogeneous terms, and (c) variable coefficients in space and time. Our approach begins with a general finite-difference discretization and applies the Schrödingerisation technique to transform the resulting system into one that admits unitary quantum evolution, enabling quantum simulation. For the Schrödinger equation corresponding to the discretized PDE, we construct an efficient block-encoding of the Hamiltonian H that scales polylogarithmically with the number of grid points N . This encoding is compatible with quantum signal processing and allows for the implementation of the evolution operator e − i H t . The explicit circuit construction in our method permits complexity to be measured in fundamental gate units–namely, CNOT gates and single-qubit rotations–bypassing the inefficiencies of oracle queries. Consequently, the overall algorithm scales polynomially with N and linearly with the spatial dimension d . Under certain input/output assumptions our method achieves a polynomial speedup in N and an exponential advantage in d for a wide class of PDEs, thereby mitigating the classical curse of dimensionality. The validity and efficiency of the proposed approach are further substantiated by numerical simulations. By explicitly defining the quantum operations and quantifying their resource requirements, our approach offers a practical alternative for numerically solving PDEs, distinct from others that rely on oracle queries and purely asymptotic scaling methods.
Quantum algorithms offer new avenues for solving partial differential equations (PDEs). While the potential for end-to-end quantum advantage is at present not well understood, recent literature presents explicit circuit constructions for solving certain classes of linear PDEs in the frequency domain and thus offers concrete examples to study. In this work, we develop end-to-end implementations of these quantum circuits compiled to machine-level instructions and benchmark them in both numerical simulations and IBMQ hardware experiments. We focus on the advection, wave, and Poisson equations and study quantum circuits that propagate the dynamics in frequency space via the quantum Fourier transform using approximate methods based on a first-order approximation which offer compact representations with uncontrollable approximation error, and polynomial approximation methods based on quantum signal processing (QSP) leading to deeper circuits with tunable algorithmic error. In addition, we experimentally demonstrate that the QSP-augmented algorithm can provide accurate solutions under realistic hardware constraints. Finally, we extend our method to address non-homogeneous Dirichlet boundary conditions and verify it numerically for a Poisson equation with source term obtained from high-fidelity physics simulations of a capacitively coupled plasma.
Partial differential equations (PDEs) are crucial for modeling various physical phenomena such as heat transfer, fluid flow, and electromagnetic waves. In computer-aided engineering (CAE), the ability to handle fine resolutions and large computational models is essential for improving product performance and reducing development costs. However, solving large-scale PDEs, particularly for systems with spatially varying material properties, poses significant computational challenges. In this paper, we propose a quantum algorithm for solving second-order linear PDEs of nonconservative systems with spatially varying parameters, using the linear combination of the Hamiltonian simulation (LCHS) method. Our approach transforms those PDEs into ordinary differential equations represented by qubit operators, through spatial discretization using the finite-difference method. Then, we provide an algorithm that efficiently constructs the operator corresponding to the spatially varying parameters of PDEs via a logic minimization technique, which reduces the number of terms and subsequently the circuit depth. We also develop a scalable method for realizing a quantum circuit for LCHS, using a tensor-network-based technique, specifically a matrix product state (MPS). We validate our method with applications to the acoustic equation with spatially varying parameters and the dissipative heat equation. Our approach includes a detailed recipe for constructing quantum circuits for PDEs, leveraging efficient encoding of spatially varying parameters of PDEs and scalable implementation of LCHS, which we believe marks a significant step towards advancing quantum computing's role in solving practical engineering problems.
This paper proposes a systematic and explicit quantum circuit framework for solving advection-diffusion equations with boundary conditions, based on the Linear Combination of Hamiltonian Simulations (LCHS) method.By employing the Finite Volume Method (FVM) combined with various flux construction schemes, we elaborate the design of quantum circuits tailored explicitly for Robin boundary conditions (including Dirichlet and Neumann boundary conditions as special cases) and periodic boundary conditions. The Robin boundary elimination is encoded by newly constructed endpoint-projector circuits, and all physical and boundary-induced inhomogeneous terms are incorporated through an outer linear combination of unitaries (LCU). The error analysis adopts an LCU-weighted form, preserving the decay of the LCHS kernel in the product-formula bounds. Complementing this analysis, a comprehensive gate complexity analysis demonstrates the quantum advantages over classical computing in high-dimensional scenarios. We simulate the proposed circuits on a fault-tolerant emulator, and numerical results validate the effectiveness of the proposed framework across homogeneous, inhomogeneous, and high-dimensional cases.The proposed framework is compatible with numerous spatial discretization methods and numerical schemes, extends naturally to other linear PDEs, and establishes a practical foundation for solving large-scale PDE problems on future fault-tolerant quantum computers.
The biharmonic equation is a prototypical fourth-order partial differential equation whose high-dimensional discretization suffers from rapidly growing degrees of freedom and severe ill-conditioning. We develop QSVT--VTAA quantum linear-system algorithms by constructing explicit block-encodings tailored to periodic, simply supported, and Dirichlet--Neumann boundary conditions. For periodic and simply supported problems, Fourier and sine-transform diagonalizations yield augmented Poisson systems with the condition-number scaling of a second-order operator. For Dirichlet--Neumann problems, we introduce a second-order boundary-corrected finite-difference discretization, establish mesh-independent stability, and construct an explicit block-encoding of the resulting nonsymmetric matrix. We also formulate a coupled-Laplace system with additional boundary unknowns and characterize its complexity in terms of the condition number of the complete augmented matrix. The analysis covers discretization error, block-encoding normalization, gate complexity, and solution extraction under an amplitude-input and quantum-state-output model. Numerical experiments validate the proposed discretizations and the corresponding linear solves.