Although quantum computing holds promise to accelerate a wide range of computational tasks, the quantum simulation of quantum dynamics as originally envisaged by Feynman remains the most promising candidate for achieving quantum advantage. A less explored possibility with comparably far-reaching technological applicability is the quantum simulation of classical nonlinear dynamics. Attempts to develop digital quantum algorithms based on the Koopman–von Neumann (KvN) formalism have met with challenges because of the necessary projection step from an infinite-dimensional Hilbert space to the finite-dimensional subspace described by a collection of qubits. This finitization produces numerical artifacts that limit solutions to very short time horizons. In this paper, we review continuous-variable quantum computing (CVQC), which naturally avoids such obstacles, and a CVQC algorithm for KvN simulation of classical nonlinear dynamics is advocated. In particular, we present explicit gate synthesis for product-formula Hamiltonian simulation of anharmonic vibrational dynamics.
Although quantum computing holds promise to accelerate a wide range of computational tasks, the quantum simulation of quantum dynamics as originally envisaged by Feynman remains the most promising candidate for achieving quantum advantage. A less explored possibility with comparably far-reaching technological applicability is the quantum simulation of classical nonlinear dynamics. Attempts to develop digital quantum algorithms based on the Koopman–von Neumann (KvN) formalism have met with challenges because of the necessary projection step from an infinite-dimensional Hilbert space to the finite-dimensional subspace described by a collection of qubits. This finitization produces numerical artifacts that limit solutions to very short time horizons. In this paper, we review continuous-variable quantum computing (CVQC), which naturally avoids such obstacles, and a CVQC algorithm for KvN simulation of classical nonlinear dynamics is advocated. In particular, we present explicit gate synthesis for product-formula Hamiltonian simulation of anharmonic vibrational dynamics.
The Koopman-von Neumann (KvN) formalism recasts classical mechanics in a Hilbert space framework using complex wavefunctions and linear operators, akin to quantum mechanics. Instead of evolving probability densities in phase space (as in Liouville's equation), KvN uses a Schrödinger-like equation for a classical wavefunction, with commuting position and momentum operators. Mapped to quantum computing, KvN offers a promising route to simulate classical dynamical systems using quantum algorithms by leveraging unitary evolution and quantum linear algebra tools, potentially enabling efficient classical-to-quantum mappings without invoking full quantum uncertainty. In this work, we specifically explore the implementation of the KvN approach on continuous-variable photonic quantum computing architectures, with the goals of leveraging quantum simulation for both sampling and computing intractable nonlinear dynamics. We will demonstrate its implementation and feasibility with two problems: the harmonic oscillator and a 1D partial differential equation governing nonlinear dynamics.
Solving Ordinary Differential Equations (ODE) is important for a broad range of domains, such as engineering, weather forecast, and finance. Most quantum algorithms proposed to solve them revolve around transforming ODEs into a system of linear equations, in order to benefit from the exponential advantage promised by the HHL algorithm. However, there are known limitations to this subroutine such as the linear scaling with condition number. In particular, for HHL-based ODE solvers, this dependency yields a complexity that generally grows exponentially with the integration time. In this paper, we present a scheme that does not rely on the HHL subroutine. We use the Koopman Von Neumann (KvN) formalism that maps arbitrary nonlinear dynamics to a Hilbert space where the dynamics are unitary. This effectively reduces the ODE problem to a Hamiltonian time evolution, which is a well-known problem in quantum computing. However, this comes at a cost as the KvN formalism is expressed in an infinite-dimensional Hilbert space, and involves nonphysical states with infinite energies. Previous works have tackled this by truncating the Hilbert space, and by approximating all the relevant states and operations on qubit-based systems. Instead of mapping to qubit-based computations, in this paper, we investigate the direct use of continuous-variable quantum computers for this problem. We provide an algorithm to compile a sequence of Gaussian and non-Gaussian continuous variables gates to solve an arbitrary one-dimensional polynomial ODE. We analyze the algorithm and propose that it is intrinsically better suited for solving so-called initial distribution problems, rather than initial condition problems. We propose the first steps towards a comprehensive complexity and specify which steps need to be developed for a complete analysis.
A number of recent studies have proposed that linear representations are appropriate for solving nonlinear dynamical systems with quantum computers, which fundamentally act linearly on a wave function in a Hilbert space. Linear representations, such as the Koopman representation and Koopman von Neumann mechanics, have regained attention from the dynamical-systems research community. Here, we aim to present a unified theoretical framework, currently missing in the literature, with which one can compare and relate existing methods, their conceptual basis, and their representations. We also aim to show that, despite the fact that quantum simulation of nonlinear classical systems may be possible with such linear representations, a necessary projection into a feasible finite-dimensional space will in practice eventually induce numerical artifacts which can be hard to eliminate or even control. As a result, a practical, reliable and accurate way to use quantum computation for solving general nonlinear dynamical systems is still an open problem.
Harmonic oscillators are promising continuous-variable (CV) quantum resources because their infinite-dimensional Hilbert spaces allow for resource-efficient quantum computing and simulation. To reach their full potential, CV platforms need to be able to efficiently implement non-Gaussian operations. Bosonic quantum-signal-processing schemes have emerged as attractive methods to construct arbitrary non-Gaussian operations; however, these schemes are restricted to single modes, i.e., the implementation of anharmonic potentials. Here, we introduce trigonometric-gate-implemented Fourier synthesis (TGIFS), a method for implementing arbitrary non-Gaussian operations applicable to both single- and multi-mode systems, allowing the generation of both anharmonicities and nonlinear multi-mode couplings. TGIFS synthesizes a target Hamiltonian by decomposing it into a Fourier series whose terms are implemented via bosonic quantum signal processing, which uses a discrete-variable (DV) system to induce a nonlinearity in the CV system. Our hybrid CV-DV protocol allows for the direct simulation of a broad range of CV phenomena (such as those in lattice gauge theory, chemical dynamics, and quantum chaos) and provides a richer toolbox for CV circuit compilation.
There is a powerful mathematical method for simulating quantum dynamics that goes beyond finite qubit representations called infinite-dimensional Hamiltonian systems. These systems make it possible to describe quantum fields, open quantum systems, and hybrid designs more accurately by simulating unbounded operators, spectral decompositions, and continuous-variable interactions. Quantum algorithms for partial differential equations, many-body physics, and error-tolerant models of infinite Hilbert spaces are all areas where they are useful. Adding Hamiltonians with infinite dimensions to computer systems makes it possible to use both classical approximation and quantum hardware. This method makes it easier to scale up, improves accuracy when dynamics aren’t simple, and sets the stage for more advanced robust quantum computing models in chemistry, information sciences, and physics.