Abstract Variational quantum algorithms (VQAs) have attracted attention as quantum algorithms employing on noisy intermediate-scale quantum devices. VQAs for solving the Poisson equation were recently proposed, as this method transforms the linear equation with a symmetric positive definite matrix obtained by discretizing the Poisson equation into a minimization problem. In this study, we propose a VQA for second-order linear differential equations including the Poisson equation on the basis of the referenced study, where the coefficient matrix obtained through discretization is a non-symmetric matrix. Furthermore, in our VQAs, we succeeded in decomposing the coefficient matrices arising from the periodic and the Dirichlet boundary condition into linear combinations of measurement operators and unitary matrices.
Abstract Variational quantum algorithms (VQAs) have attracted attention as quantum algorithms employing on noisy intermediate-scale quantum devices. VQAs for solving the Poisson equation were recently proposed, as this method transforms the linear equation with a symmetric positive definite matrix obtained by discretizing the Poisson equation into a minimization problem. In this study, we propose a VQA for second-order linear differential equations including the Poisson equation on the basis of the referenced study, where the coefficient matrix obtained through discretization is a non-symmetric matrix. Furthermore, in our VQAs, we succeeded in decomposing the coefficient matrices arising from the periodic and the Dirichlet boundary condition into linear combinations of measurement operators and unitary matrices.
Abstract To solve the Poisson equation it is usually possible to discretize it into solving the corresponding linear system Ax = b. Variational quantum algorithms (VQAs) for the discretized Poisson equation have been studied before. We present a VQA based on the banded Toeplitz systems for solving the Poisson equation with respect to the structural features of matrix A. In detail, we decompose the matrices A and A 2 into a linear combination of the corresponding banded Toeplitz matrix and sparse matrices with only a few non-zero elements. For the one-dimensional Poisson equation with different boundary conditions and the d-dimensional Poisson equation with Dirichlet boundary conditions, the number of decomposition terms is less than that reported in [Phys. Rev. A 2023 108, 032 418 ]. Based on the decomposition of the matrix, we design quantum circuits that efficiently evaluate the cost function. Additionally, numerical simulation verifies the feasibility of the proposed algorithm. Finally, the VQAs for linear systems of equations and matrix-vector multiplications with the K-banded Toeplitz matrix T n K are given, where T n K ∈ R n × n and K ∈ O ( ploylog n ) .
The Poisson equation has wide applications in many areas of science and engineering. Although there are some quantum algorithms that can efficiently solve the Poisson equation, they generally require a fault-tolerant quantum computer, which is beyond the current technology. We propose a variational quantum algorithm (VQA) to solve the Poisson equation, which can be executed on noisy intermediate-scale quantum devices. In detail, we first adopt the finite-difference method to transform the Poisson equation into a linear system. Then, according to the special structure of the linear system, we find an explicit tensor product decomposition, with only (2log2n+1) items, of its coefficient matrix under a specific set of simple operators, where n is the dimension of the coefficient matrix. This implies that the proposed VQA needs fewer quantum measurements, which dramatically reduces the required quantum resources. Additionally, we design observables to efficiently evaluate the expectation values of the simple operators on a quantum computer. Numerical experiments demonstrate that our algorithm can solve the Poisson equation.
Solving a Poisson equation is generally reduced to solving a linear system with a coefficient matrix A of entries aij,i,j=1,2,...,n, from the discretized Poisson equation. Although the variational quantum algorithms are promising algorithms to solve the discretized Poisson equation, they generally require that A be decomposed into a sum of O[poly(log2n)] simple operators in order to evaluate efficiently the loss function. A tensor product decomposition of A with 2log2n+1 terms has been explored in previous works. In this paper, based on the decomposition of sparse Hamiltonians we greatly reduce the number of terms. We first write the loss function in terms of the operator ${\ensuremath{\sigma}}_{x}\ensuremath{\bigotimes}A$ with ${\ensuremath{\sigma}}_{x}$ denoting the standard Pauli operator. Then for the one-dimensional Poisson equations with different boundary conditions and for the d-dimensional Poisson equations with Dirichlet boundary conditions, we decompose ${\ensuremath{\sigma}}_{x}\ensuremath{\bigotimes}A$ into a sum of at most 7 and (4d+1) Hermitian, one-sparse, and self-inverse operators, respectively. We design explicitly the quantum circuits to evaluate efficiently the loss function. The decomposition method and the design of quantum circuits can also be easily extended to linear systems with Hermitian and sparse coefficient matrices satisfying ai,i+c=ac for $c=0,1,...,n\ensuremath{-}1$ and $i=0,\ensuremath{\cdots},n\ensuremath{-}1\ensuremath{-}c$.
As quantum computing technology continues to develop, quantum algorithms have the potential to outperform classical algorithms in future for certain computational tasks. One of these potential applications of quantum computing is numerically solving differential equations and integral equations. In this thesis the variational quantum linear solver algorithm (VQLS) proposed by Bravo-Prieto et al. is adapted to solve the Laplace and Poisson equations, ordinary differential equation initial value problems (IVPs) with both constant and function coefficients, and a Fredholm integral equation of the second type. As this thesis demonstrates, variational quantum algorithms for solving differential equations and integral equations such as VQLS are highly promising for use on Noisy Intermediate-Scale Quantum (NISQ) devices. The quantum circuit components of the VQLS algorithm were simulated using IBM Quantum's statevector simulator with the optimized solution circuits implemented on the IBM Quantum Guadalupe device (in most cases). A discussion of errors and convergence is also included. The application of VQLS to these differential equations and integral equations is novel and, to my knowledge, work using VQLS to solve the Poisson equation has only recently been published once by authors Liu et al. which was conducted concurrently and independently of this research.
Computer-aided engineering techniques are indispensable in modern engineering developments. In particular, partial differential equations are commonly used to simulate the dynamics of physical phenomena, but very large systems are often intractable within a reasonable computation time, even when using supercomputers. To overcome the inherent limit of classical computing, we present a variational quantum algorithm for solving the Poisson equation that can be implemented in noisy intermediate-scale quantum devices. The proposed method defines the total potential energy of the Poisson equation as an expectation of certain observables, which are decomposed into a linear combination of Pauli operators and simple observables. The expectation value of these observables is then minimized with respect to a parameterized quantum state. Because the number of decomposed terms is independent of the size of the problem, this method requires relatively few quantum measurements. Numerical experiments demonstrate the faster computing speed of this method compared with classical computing methods and a previous variational quantum approach. We believe that our approach brings quantum computer-aided techniques closer to future applications in engineering developments.
Differential equations underpin a vast range of scientific and engineering domains, including structural mechanics, fluid dynamics, and financial modelling. Their intricate nature often renders classical numerical methods computationally demanding, requiring fine discretization or large basis expansions, which can lead to prohibitive costs. Quantum computing offers a fundamentally different paradigm, leveraging superposition and entanglement with the promise of computational speedup. However, in the current Noisy Intermediate-Scale Quantum (NISQ) era, fully fault-tolerant quantum computers remain decades away, motivating the development of hybrid strategies such as Variational Quantum Algorithms (VQAs), which train parametrised quantum circuits using classical optimizers. This work investigates spectral VQAs for solving differential equations, extending the Chebyshev-encoded VQA of Kyriienko et al. and introducing a novel quantum spectral approach: the Lagrange-encoded VQA. This approach shifts the paradigm by incorporating entanglement directly within the encoding block while maintaining a shallow variational ansatz. Inspired by the Hadamard test structure, the design enables an efficient representation of the complete set of Lagrange polynomial basis functions using a limited number of gates and facilitates the use of the Hadamard test differentiation method, thereby reducing the computational load. VQA performance is analysed through the accuracy of the result in combination with the gate count and circuit count metrics, which serve as proxies for runtime and noise accumulation on NISQ devices. The proposed VQA is benchmarked on multiple problems: (i) a first-order linear ODEs and (ii) the second-order ODEs damped mass-spring system, both comparing Lagrange- and Chebyshev-encoded VQAs, (iii) the Poisson equation, comparing the Lagrange-encoded VQA with the discretized approach of Sato et al., and (iv) a first-order ODE with nonlinear coefficients addressed with the Lagrange-encoded VQA. Across these cases, the Lagrange-encoded VQA consistently achieves comparable or improved accuracy while significantly reducing circuit depth and computational overload, making it a promising candidate for NISQ devices.