Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1909.03898 · 中国科学通报：英文版

Variational algorithms for linear algebra

Xiaosi Xu, Jinzhao Sun, Suguru Endo, Ying Li, Simon C. Benjamin, Xiao Yuan

Published Jun 26, 2021• 164 citations
arXivHTMLPDF

Abstract

Quantum algorithms have been developed for efficiently solving linear algebra tasks. However, they generally require deep circuits and hence universal fault-tolerant quantum computers. In this work, we propose variational algorithms for linear algebra tasks that are compatible with noisy intermediate-scale quantum devices. We show that the solutions of linear systems of equations and matrix-vector multiplications can be translated as the ground states of the constructed Hamiltonians. Based on the variational quantum algorithms, we introduce Hamiltonian morphing together with an adaptive ansätz for efficiently finding the ground state, and show the solution verification. Our algorithms are especially suitable for linear algebra problems with sparse matrices, and have wide applications in machine learning and optimisation problems. The algorithm for matrix multiplications can be also used for Hamiltonian simulation and open system simulation. We evaluate the cost and effectiveness of our algorithm through numerical simulations for solving linear systems of equations. We implement the algorithm on the IBM quantum cloud device with a high solution fidelity of 99.95%.

DOI
10.1016/j.scib.2021.06.023

Related papers

Variational Quantum Simulation of General Processes

Suguru Endo, Jinzhao Sun, Ying Li +2 more

Variational quantum algorithms have been proposed to solve static and dynamic problems of closed many-body quantum systems. Here we investigate variational quantum simulation of three general types of tasks-generalized time evolution with a non-Hermitian Hamiltonian, linear algebra problems, and open quantum system dynamics. The algorithm for generalized time evolution provides a unified framework for variational quantum simulation. In particular, we show its application in solving linear systems of equations and matrix-vector multiplications by converting these algebraic problems into generalized time evolution. Meanwhile, assuming a tensor product structure of the matrices, we also propose another variational approach for these two tasks by combining variational real and imaginary time evolution. Finally, we introduce variational quantum simulation for open system dynamics. We variationally implement the stochastic Schrödinger equation, which consists of dissipative evolution and stochastic jump processes. We numerically test the algorithm with a 6-qubit 2D transverse field Ising model under dissipation.

arXiv:1812.08778•Physical Review Letters•Jun 29, 2020•237 citations

Solving a Linear System of Equations on a Quantum Computer by Measurement

Alain Giresse Tene, Thomas Konrad

We present a variational algorithm for fault tolerant quantum computing to solve a system of linear equations which directly maximises the parameters of the target fidelity. This so-called measurement test algorithm can be applied to any computational task with a solution that is represented as eigenvector of a self-adjoint matrix. The solution is prepared as state of a register in the quantum computer by a von Neumann measurement of a corresponding observable, which is implemented using the phase estimation algorithm. The probability to project the system thus into the unknown target state, which equals the target fidelity, is measured in terms of relative frequencies and iteratively optimised to read out the target state. The new algorithm overcomes three issues of previous variational quantum algorithms: i) It does not rely on a decomposition in terms of Pauli strings and therefore can compute eigenvectors of dense matrices. ii) The accuracy is not limited by the condition number $κ$ of the matrix, provided a logarithmic number ($O(\logκ)$) of qubits is used to encode the eigenvalues and iii) the target fidelity $F_T = 1-ε$ can be reached with an accuracy $ε$ that scales with $1/N$ for $N$ measurements per iteration. We demonstrate this by numerical simulations for dense random real-valued $16\times 16$ matrices with non-vanishing determinant.

quant-pharXiv:2604.26098•Apr 30, 2026•0 citations

A Variational Approach to Solving Differential Equations on a Quantum Computer

James Nelson

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.

•Jan 1, 2021•0 citations

Near term algorithms for linear systems of equations

Aidan Pellow-Jarman, Ilya Sinayskiy, Anban W. Pillay +1 more

Abstract Finding solutions to systems of linear equations is a common problem in many areas of science and engineering, with much potential for a speed up on quantum devices. While the Harrow–Hassidim–Lloyd (HHL) quantum algorithm yields up to an exponential speed up over classical algorithms in some cases, it requires a fault-tolerant quantum computer, which is unlikely to be available in the near-term. Thus, attention has turned to the investigation of quantum algorithms for noisy intermediate-scale quantum (NISQ) devices where several near-term approaches to solving systems of linear equations have been proposed. This paper focuses on the Variational Quantum Linear Solvers (VQLS), and other closely related methods and adaptions. Several contributions are made in this paper, which include: the first application of the Evolutionary Ansatz to the VQLS (EAVQLS), the first implementation of the Logical Ansatz VQLS (LAVQLS), based on the Classical Combination of Quantum States (CQS) method, a proof of principle demonstration of the CQS method on real quantum hardware and a method for the implementation of the Adiabatic Ansatz on the VQLS (AAVQLS). These approaches are implemented and contrasted. The CQS method is run with moderate success on a real quantum device. The EAVQLS and AAVQLS show promise as possible improvements to the standard VQLS algorithm once refined.

arXiv:2108.11362•Quantum Information Processing•Jun 24, 2023•7 citations

Variational quantum algorithm for the Poisson equation

Hailing Liu, Yusen Wu, Lin‐Chun Wan +4 more

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 $(2{log}_{2}n+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.

arXiv:2012.07014•Physical Review A•Aug 18, 2021•127 citations

A quantum algorithm for linear differential equations with layerwise parameterized quantum circuits

Junxiang Xiao, Jingwei Wen, Zengrong Zhou +4 more

Abstract Solving linear differential equations is a common problem in almost all fields of science and engineering. Here, we present a variational algorithm with shallow circuits for solving such a problem: given an $$N \times N$$ N × N matrix $${\varvec{A}}$$ A , an N-dimensional vector $$\varvec{b}$$ b , and an initial vector $$\varvec{x}(0)$$ x ( 0 ) , how to obtain the solution vector $$\varvec{x}(T)$$ x ( T ) at time T according to the constraint $$\textrm{d}\varvec{x}(t)/\textrm{d} t = {\varvec{A}}\varvec{x}(t) + \varvec{b}$$ d x ( t ) / d t = A x ( t ) + b . The core idea of the algorithm is to encode the equations into a ground state problem of the Hamiltonian, which is solved via hybrid quantum-classical methods with high fidelities. Compared with the previous works, our algorithm requires the least qubit resources and can restore the entire evolutionary process. In particular, we show its application in simulating the evolution of harmonic oscillators and dynamics of non-Hermitian systems with $$\mathcal{P}\mathcal{T}$$ P T -symmetry. Our algorithm framework provides a key technique for solving so many important problems whose essence is the solution of linear differential equations.

AAPPS bulletin•Mar 4, 2024•3 citations
Quantum algorithms have been developed for efficiently solving linear algebra tasks. However, they generally require deep circuits and hence universal fault-tolerant quantum computers. In this work, we propose variational algorithms for linear algebra tasks that are compatible with noisy intermediate-scale quantum devices. We show that the solutions of linear systems of equations and matrix-vector multiplications can be translated as the ground states of the constructed Hamiltonians. Based on the variational quantum algorithms, we introduce Hamiltonian morphing together with an adaptive ansätz for efficiently finding the ground state, and show the solution verification. Our algorithms are especially suitable for linear algebra problems with sparse matrices, and have wide applications in machine learning and optimisation problems. The algorithm for matrix multiplications can be also used for Hamiltonian simulation and open system simulation. We evaluate the cost and effectiveness of our algorithm through numerical simulations for solving linear systems of equations. We implement the algorithm on the IBM quantum cloud device with a high solution fidelity of 99.95%.
Variational quantum algorithms have been proposed to solve static and dynamic problems of closed many-body quantum systems. Here we investigate variational quantum simulation of three general types of tasks-generalized time evolution with a non-Hermitian Hamiltonian, linear algebra problems, and open quantum system dynamics. The algorithm for generalized time evolution provides a unified framework for variational quantum simulation. In particular, we show its application in solving linear systems of equations and matrix-vector multiplications by converting these algebraic problems into generalized time evolution. Meanwhile, assuming a tensor product structure of the matrices, we also propose another variational approach for these two tasks by combining variational real and imaginary time evolution. Finally, we introduce variational quantum simulation for open system dynamics. We variationally implement the stochastic Schrödinger equation, which consists of dissipative evolution and stochastic jump processes. We numerically test the algorithm with a 6-qubit 2D transverse field Ising model under dissipation.
We present a variational algorithm for fault tolerant quantum computing to solve a system of linear equations which directly maximises the parameters of the target fidelity. This so-called measurement test algorithm can be applied to any computational task with a solution that is represented as eigenvector of a self-adjoint matrix. The solution is prepared as state of a register in the quantum computer by a von Neumann measurement of a corresponding observable, which is implemented using the phase estimation algorithm. The probability to project the system thus into the unknown target state, which equals the target fidelity, is measured in terms of relative frequencies and iteratively optimised to read out the target state. The new algorithm overcomes three issues of previous variational quantum algorithms: i) It does not rely on a decomposition in terms of Pauli strings and therefore can compute eigenvectors of dense matrices. ii) The accuracy is not limited by the condition number κκκ of the matrix, provided a logarithmic number (O(log⁡κ)O(\logκ)O(logκ)) of qubits is used to encode the eigenvalues and iii) the target fidelity FT=1−εF_T = 1-εFT​=1−ε can be reached with an accuracy εεε that scales with 1/N1/N1/N for NNN measurements per iteration. We demonstrate this by numerical simulations for dense random real-valued 16×1616\times 1616×16 matrices with non-vanishing determinant.
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.
Abstract Finding solutions to systems of linear equations is a common problem in many areas of science and engineering, with much potential for a speed up on quantum devices. While the Harrow–Hassidim–Lloyd (HHL) quantum algorithm yields up to an exponential speed up over classical algorithms in some cases, it requires a fault-tolerant quantum computer, which is unlikely to be available in the near-term. Thus, attention has turned to the investigation of quantum algorithms for noisy intermediate-scale quantum (NISQ) devices where several near-term approaches to solving systems of linear equations have been proposed. This paper focuses on the Variational Quantum Linear Solvers (VQLS), and other closely related methods and adaptions. Several contributions are made in this paper, which include: the first application of the Evolutionary Ansatz to the VQLS (EAVQLS), the first implementation of the Logical Ansatz VQLS (LAVQLS), based on the Classical Combination of Quantum States (CQS) method, a proof of principle demonstration of the CQS method on real quantum hardware and a method for the implementation of the Adiabatic Ansatz on the VQLS (AAVQLS). These approaches are implemented and contrasted. The CQS method is run with moderate success on a real quantum device. The EAVQLS and AAVQLS show promise as possible improvements to the standard VQLS algorithm once refined.
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)(2{log}_{2}n+1)(2log2​n+1) items, of its coefficient matrix under a specific set of simple operators, where nnn 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.
Abstract Solving linear differential equations is a common problem in almost all fields of science and engineering. Here, we present a variational algorithm with shallow circuits for solving such a problem: given an N×NN \times NN×N N × N matrix $${\varvec{A}}$$ A , an N-dimensional vector $$\varvec{b}$$ b , and an initial vector $$\varvec{x}(0)$$ x ( 0 ) , how to obtain the solution vector $$\varvec{x}(T)$$ x ( T ) at time T according to the constraint $$\textrm{d}\varvec{x}(t)/\textrm{d} t = {\varvec{A}}\varvec{x}(t) + \varvec{b}$$ d x ( t ) / d t = A x ( t ) + b . The core idea of the algorithm is to encode the equations into a ground state problem of the Hamiltonian, which is solved via hybrid quantum-classical methods with high fidelities. Compared with the previous works, our algorithm requires the least qubit resources and can restore the entire evolutionary process. In particular, we show its application in simulating the evolution of harmonic oscillators and dynamics of non-Hermitian systems with PT\mathcal{P}\mathcal{T}PT P T -symmetry. Our algorithm framework provides a key technique for solving so many important problems whose essence is the solution of linear differential equations.