We propose a distinct approach to solving linear and nonlinear differential equations (DEs) on quantum computers by encoding the problem into ground states of effective Hamiltonian operators. Our algorithm relies on constructing such operators in the Chebyshev space, where an effective Hamiltonian is a sum of global differential and data constraints. Once the effective Hamiltonian is formed, solutions of differential equations can be obtained using the ground state preparation techniques (e.g. imaginary-time evolution and quantum singular value transformation), bypassing variational search. Unlike approaches based on discrete grids, the algorithm enables evaluation of solutions beyond fixed grid points and implements constraints in the physics-informed way. Our proposal inherits the best traits from quantum machine learning-based DE solving (compact basis representation, automatic differentiation, nonlinearity) and quantum linear algebra-based approaches (fine-grid encoding, provable speed-up for state preparation), offering a robust strategy for quantum scientific computing in the early fault-tolerant era.
We propose a physics-informed quantum algorithm to solve nonlinear and multidimensional differential equations (DEs) in a quantum latent space. We suggest a strategy for building quantum models as state overlaps, where exponentially large sets of independent basis functions are used for implicitly representing solutions. By measuring the overlaps between states which are representations of DE terms, we construct a loss that does not require independent sequential function evaluations on grid points. In this sense, the solver evaluates the loss in an intrinsically parallel way, utilizing a global type of the model. When the loss is trained variationally, our approach can be related to the differentiable quantum circuit protocol, which does not scale with the training grid size. Specifically, using the proposed model definition and feature map encoding, we represent function- and derivative-based terms of a differential equation as corresponding quantum states. Importantly, we propose an efficient way for encoding nonlinearity, for some bases requiring only an additive linear increase of the system size
O(N+p) in the degree of nonlinearity
p. By utilizing basis mapping, we show how the proposed model can be evaluated explicitly. This allows to implement arbitrary functions of independent variables, treat problems with various initial and boundary conditions, and include data and regularization terms in the physics-informed machine learning setting. On the technical side, we present toolboxes for exponential Chebyshev and Fourier basis sets, developing tools for automatic differentiation and multiplication, implementing nonlinearity, and describing multivariate extensions. The approach is compatible with, and tested on, a range of problems including linear, nonlinear and multidimensional differential equations.
Differential equations can be encoded as ground-state problems by constructing a positive-semidefinite effective Hamiltonian whose minimum-energy state represents the solution. We extend this framework in several directions. First, we show how multiple nonzero initial, boundary, and data conditions can be incorporated into the homogeneous formulation
A∣f⟩=0 using a known nonzero reference condition. We then analyse an alternative formulation based on
A∣f⟩=∣b⟩, with Hamiltonian
Hb=A†(I−∣b⟩⟨b∣)A, which incorporates source terms and nonzero constraints through an augmented system. For nonlinear differential equations, we examine the unphysical ground-state degeneracy introduced by tensor-product representations and consider both Hamiltonian constraints and ansatz-level restrictions for targeting physically valid product states. Finally, we develop grid-value amplitude-encoded versions of both Hamiltonian constructions and compare them with spectral coefficient encoding. Through linear and nonlinear examples, we assess how the choice of Hamiltonian formulation and function representation affects solution recovery, spectral gap, degeneracy, and readout. These results broaden the applicability of ground-state-based quantum differential-equation solvers while clarifying their principal practical trade-offs.
Quantum algorithms based on linear-system approaches for solving differential equations demand qubit and precision resources beyond near-term capabilities. To address these challenges, this work proposes a physics-informed quantum machine learning (PIQML) framework with hard constraint embedding, specifically designed for NISQ era. Within this framework, parameterized quantum circuits serve as machine learning models, where the input variable is encoded into a high-dimensional feature space via a Fourier feature map. Subsequently, to eliminate approximation errors in critical physical conditions, the solution is constructed through a rigorously designed function mapper that analytically enforces initial conditions as hard constraints. Crucially, we compute derivatives with respect to the input variable using the parameter-shift rule---a quantum native gradient evaluation technique that avoids classical discretization. Unlike generic loss functions that target abstract data patterns, our loss function focuses on the differential equation residual and reference data. This design ensures that the trained model not only approximates the data but also intrinsically satisfies the physical constraint expressed by the DE itself. Our method is validated on several differential equations, including highly oscillatory ones, demonstrating its capability to tackle challenging nonlinear dynamics. Results demonstrate that our quantum model successfully learns the solution, showing close agreement with a high-precision classical numerical benchmark.
As quantum computing moves toward scientific computing applications, nonlinear differential equations remain a central challenge since quantum evolution is intrinsically linear. In this work, we introduce a reduced basis algorithm (RBA) for polynomial nonlinear ordinary differential equations (ODEs) and spatially discretized partial differential equations (PDEs). After time discretization, the method composes the resulting polynomial update map over
m timesteps, identifies the reduced monomial basis appearing in this composed map, and constructs a linear RBA operator whose action recovers the exact
m-timestep nonlinear dynamics. Thus, at the level of the chosen discrete update rule, the method introduces no additional approximation error beyond the time discretization error. The qubit number requirement is governed by the size of the reduced monomial basis. For an
n-dimensional polynomial ODE system of degree
p>1, the lifted register requires at most
qmODE=O(nmlogp) qubits in the full basis scenario. For PDEs discretized on
ND grid points, a locality-based construction requires at most
qmPDE=O(DlogN+nmD+1logp) qubits. Hence, the dependence on the grid size remains logarithmic, while the nonlinear overhead is controlled by local reduced basis size. The main computational burden is moved from the quantum computer to a classical preprocessing step, where the reduced monomial basis and RBA operator are constructed for the chosen timestep window. Through numerical tests on the Lorenz system and the one-dimensional Burgers equation, we verify that the RBA reproduces the corresponding discrete time nonlinear dynamics exactly, while exposing the trade-off between timestep composition, reduced basis growth, and locality.
From weather to neural networks, modeling is not only useful for understanding various phenomena but also has a wide range of potential applications. Although nonlinear differential equations are extremely useful tools in modeling, their solutions are difficult to obtain. Based on the expectation of quantum transcendence, quantum algorithms for efficiently solving nonlinear differential equations continue to be developed. However, even the latest promising algorithm has been pointed out to have an evolution time limit. This limit is the theoretically predestined divergence of solutions. We propose algorithms of divergence-free simulation for nonlinear differential equations in quantum computers. For Hamiltonian simulations, a pivot state s is introduced in the neighborhood of state x to be solved. Divergence of the solutions is prevented by moving s to a neighborhood of x whenever x leaves the neighborhood of s . Since updating s is directly related to computational cost, to minimize the number of updates, the nonlinear differential equations are approximated by nonlinear polynomials around s , which are then Carleman linearized. Hamiltonian simulations of nonlinear differential equations based on several representative models are performed to show that the proposed methods break through the theoretical evolution time limit. Furthermore, we demonstrate that the proposed methods work well on an actual quantum computer. The solution of nonlinear differential equations free from evolution time constraints opens the door to practical applications of quantum computers.
Quantum computers have long been expected to efficiently solve complex classical differential equations. Most digital, fault-tolerant approaches use Carleman linearization to map nonlinear systems to linear ones and then apply quantum linear-system solvers. However, provable speedups typically require digital truncation and full fault tolerance, rendering such linearization approaches challenging to implement on current hardware. Here we present an analog, continuous-variable algorithm based on coupled bosonic modes with qubit-based adaptive measurements that avoids Hilbert-space digitization. This method encodes classical fields as coherent states and, via Kraus-channel composition derived from the Koopman-von Neumann (KvN) formalism, maps nonlinear evolution to linear dynamics. Unlike many analog schemes, the algorithm is provably efficient: advancing a first-order,
L-grid point,
d-dimensional, order-
K spatial-derivative, degree-
r polynomial-nonlinearity, strongly dissipative partial differential equations (PDEs) for
T time steps costs
O(T(logL+drlogK)). The capability of the scheme is demonstrated by using it to simulate the one-dimensional Burgers' equation and two-dimensional Fisher-KPP equation. The resilience of the method to photon loss is shown under strong-dissipation conditions and an analytic counterterm is derived that systematically cancels the dominant, experimentally calibrated noise. This work establishes a continuous-variable framework for simulating nonlinear systems and identifies a viable pathway toward practical quantum speedup on near-term analog hardware.