The role of differential equations (DEs) in science and engineering is of paramount importance, as they provide the mathematical framework for a multitude of natural phenomena. Since quantum computers promise significant advantages over classical computers, quantum algorithms for the solution of DEs have received a lot of attention. Particularly interesting are algorithms that offer advantages in the current noisy intermediate-scale quantum (NISQ) era, characterized by small and error-prone systems. We consider a framework of variational quantum algorithms, quantum circuit learning (QCL), in conjunction with derivation methods, in particular the parameter shift rule, to solve DEs. As these algorithms were specifically designed for NISQ computers, we analyze their applicability on NISQ devices by implementing QCL on an IBM quantum computer. Our analysis of QCL without the parameter shift rule shows that we can successfully learn different functions with three-qubit circuits. However, the hardware errors accumulate with increasing number of qubits, and thus, only a fraction of the qubits available on the current quantum systems can be effectively used. We further show that it is possible to determine derivatives of the learned functions using the parameter shift rule on the IBM hardware. The parameter shift rule results in higher errors, which limits its execution to low-order derivatives. Despite these limitations, we solve a first-order DE on the IBM quantum computer. We further explore the advantages of using multiple qubits in QCL by learning different functions simultaneously and demonstrate the solution of a coupled DE on a simulator.
The role of differential equations (DEs) in science and engineering is of paramount importance, as they provide the mathematical framework for a multitude of natural phenomena. Since quantum computers promise significant advantages over classical computers, quantum algorithms for the solution of DEs have received a lot of attention. Particularly interesting are algorithms that offer advantages in the current noisy intermediate-scale quantum (NISQ) era, characterized by small and error-prone systems. We consider a framework of variational quantum algorithms, quantum circuit learning (QCL), in conjunction with derivation methods, in particular the parameter shift rule, to solve DEs. As these algorithms were specifically designed for NISQ computers, we analyze their applicability on NISQ devices by implementing QCL on an IBM quantum computer. Our analysis of QCL without the parameter shift rule shows that we can successfully learn different functions with three-qubit circuits. However, the hardware errors accumulate with increasing number of qubits, and thus, only a fraction of the qubits available on the current quantum systems can be effectively used. We further show that it is possible to determine derivatives of the learned functions using the parameter shift rule on the IBM hardware. The parameter shift rule results in higher errors, which limits its execution to low-order derivatives. Despite these limitations, we solve a first-order DE on the IBM quantum computer. We further explore the advantages of using multiple qubits in QCL by learning different functions simultaneously and demonstrate the solution of a coupled DE on a simulator.
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.
Quantum computing has evolved into a promising computational paradigm capable of addressing problems beyond the reach of classical methods. Among current approaches, hybrid quantum–classical algorithms—and in particular Variational Quantum Algorithms (VQAs)—play a central role in the Noisy Intermediate-Scale Quantum (NISQ) era. This thesis analyzes the research conducted during the PhD within the variational framework along two main directions. First, it introduces the Quantum Non-Demolition Measurement (QNDM) method, a novel strategy for efficient gradient estimation in VQAs. By coupling an auxiliary detector qubit through non-demolition operators, QNDM enables derivative extraction from a single measurement, reducing quantum resource requirements by up to an order of magnitude compared to standard direct-measurement methods. The second part explores the application of quantum machine learning to the solution of differential equations, extending the Differential Quantum Circuit (DQC) framework to stochastic systems. This work was carried out during a one-year research internship at the French scale-up Pasqal.
The potential of quantum computing for scientific and industrial breakthroughs is immense, however, we are still in the Noisy Intermediate-Scale Quantum (NISQ) era, where the currently available quantum devices contain small numbers of qubits, are very sensitive to environmental conditions and prone to quantum decoherence. Even so, existing NISQ computers have already been shown to outperform conventional computers on specific problems and the key question is how to make use of today’s NISQ devices to achieve quantum advantage in the field of computational science and engineering (CSE). In this direction, this work proposes a hybrid computing formulation by combining quantum computing with machine learning for accelerating the solution of parameterized linear systems in NISQ devices. In particular, it focuses on the Variational Quantum Linear Solver (VQLS), which is a hybrid quantum- classical algorithm for solving linear systems of equations. VQLS employs a short-depth quantum circuit to efficiently evaluate a cost function related to the system solution. The algorithm’s circuit consists of a quantum gate sequence (unitary operators) that involves a set of tunable parameters. Then, well-established classical optimizers are utilized to tune the parameters of the sequence in order to minimize the cost function, which is equivalent to finding the system solution at an acceptable level of accuracy. In this work, it is demonstrated that machine learning tools such as feed-forward neural networks and nearest-neighbor interpolation techniques can be successfully employed to accelerate the convergence of the VQLS algorithm towards the optimal values for the circuit parameters when applied to parameterized linear systems that need to be solved for multiple parameter instances.
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.
Noisy Intermediate-Scale Quantum (NISQ) devices show great promise for quantum advantage via Parametrized Quantum Circuits (PQCs) and Variational Quantum Algorithms (VQAs). Due to hardware imperfections, however, especially barren plateaus and environmental Pauli noise, the cost landscape is highly distorted, resulting in inaccurate classical optimization during training. To overcome these challenges, a robust quantum gradient descent framework is presented with the integration of error mitigation directly within the optimization loop. A noise-resilient parameter-shift rule is defined and modeled analytically, as is a variant of stochastic gradient descent with real-time Pauli noise mitigation. The framework is based on a 6-qubit system, which is tested on the Max-Cut problem in the presence of different amounts of depolarizing noise. It is shown experimentally that the proposed algorithm is able to achieve the final cost error of 0.045 under significant noise, whereas the standard quantum gradient descent fails to achieve better than 0.320 under this noise. Moreover, the methodology is consistent, converging without a noise-induced flat landscape and without local minima. The method offers a scalable route towards future applications of quantum machine learning on non-fault-tolerant hardware that does not suffer from the large sampling overhead often associated with conventional error mitigation techniques.
Solving and optimizing differential equations (DEs) is ubiquitous in both engineering and fundamental science. The promise of quantum architectures to accelerate scientific computing thus naturally involved interest towards how efficiently quantum algorithms can solve DEs. Differentiable quantum circuits (DQC) offer a viable route to compute DE solutions using a variational approach amenable to existing quantum computers, by producing a machine-learnable surrogate of the solution. Quantum extremal learning (QEL) complements such approach by finding extreme points in the output of learnable models of unknown (implicit) functions, offering a powerful tool to bypass a full DE solution, in cases where the crux consists in retrieving solution extrema. In this work, we provide the results from the first experimental demonstration of both DQC and QEL, displaying their performance on a synthetic usecase. Whilst both DQC and QEL are expected to require digital quantum hardware, we successfully challenge this assumption by running a closed-loop instance on a commercial analog quantum computer, based upon neutral atom technology.