Recognition: unknown
Numerical approximation of the Koopman-von Neumann equation: Operator learning and quantum computing
Pith reviewed 2026-05-10 17:18 UTC · model grok-4.3
The pith
The Koopman-von Neumann operator projects to a unitary matrix that quantum circuits can implement and data can approximate.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The Koopman-von Neumann operator is unitary even if the dynamics are non-Hamiltonian. Projecting this operator onto a finite-dimensional subspace allows us to represent it by a unitary matrix, which in turn can be expressed as a quantum circuit. We exploit relationships between the Koopman-von Neumann framework and classical transfer operators in order to derive numerical methods to approximate the Koopman-von Neumann operator and its eigenvalues and eigenfunctions from data. Furthermore, the choice of basis functions and domain are crucial to ensure that the operator is well-defined.
What carries the argument
Finite-dimensional projection of the Koopman-von Neumann operator onto a subspace of basis functions, producing a unitary matrix that admits a quantum-circuit representation.
If this is right
- Data-driven approximations recover the operator, eigenvalues, and eigenfunctions for both undamped and damped systems.
- The same projection technique applies to nonlinear models such as Lotka-Volterra.
- The resulting unitary matrix supplies an explicit quantum-circuit implementation of the evolution.
- Relationships with Koopman and Perron-Frobenius operators supply concrete numerical schemes for operator learning.
Where Pith is reading between the lines
- The approach may scale to high-dimensional or partially observed systems where classical matrix methods become prohibitive.
- Hybrid quantum-classical pipelines could use the learned unitary to accelerate long-time predictions or control tasks.
- Similar projections might be tested on stochastic differential equations to check whether unitarity is preserved under noise.
Load-bearing premise
Suitable basis functions and domains exist such that the projected operator stays unitary and the data-driven approximations converge.
What would settle it
For the damped oscillator or Lotka-Volterra examples, the finite matrix obtained from standard polynomial or Fourier bases either fails to remain unitary or its computed eigenvalues deviate from the known spectral properties of the true operator.
read the original abstract
The Koopman-von Neumann equation describes the evolution of wavefunctions associated with autonomous ordinary differential equations and can be regarded as a quantum physics-inspired formulation of classical mechanics. The main advantage compared to conventional transfer operators such as Koopman and Perron-Frobenius operators is that the Koopman-von Neumann operator is unitary even if the dynamics are non-Hamiltonian. Projecting this operator onto a finite-dimensional subspace allows us to represent it by a unitary matrix, which in turn can be expressed as a quantum circuit. We will exploit relationships between the Koopman-von Neumann framework and classical transfer operators in order to derive numerical methods to approximate the Koopman-von Neumann operator and its eigenvalues and eigenfunctions from data. Furthermore, we will show that the choice of basis functions and domain are crucial to ensure that the operator is well-defined. We will illustrate the results with the aid of guiding examples, including simple undamped and damped oscillators and the Lotka-Volterra model.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops numerical methods for approximating the Koopman-von Neumann (KvN) operator for autonomous ODEs. It argues that the KvN operator remains unitary even for non-Hamiltonian dynamics, so that its projection onto a finite-dimensional subspace yields a unitary matrix representable as a quantum circuit. Connections to classical transfer operators are used to obtain data-driven approximations of the operator, its eigenvalues, and eigenfunctions. Basis and domain selection are identified as critical for well-definedness, with concrete illustrations on undamped and damped oscillators plus the Lotka-Volterra system.
Significance. If the unitarity and convergence claims are placed on a rigorous footing, the work would supply a concrete bridge between classical dynamical systems and quantum operator learning, exploiting unitarity where standard Koopman or Perron-Frobenius operators lose it. The data-driven component aligns with existing operator-learning literature and could enable spectral analysis of non-Hamiltonian flows on quantum hardware once suitable bases are constructible.
major comments (2)
- [§3] §3 (Finite-dimensional projection): the statement that projection onto any finite subspace produces a unitary matrix is not accompanied by a general theorem. Unitarity holds only when the subspace is invariant under the KvN operator and the basis is orthonormal with respect to the appropriate measure; the manuscript demonstrates this on the oscillator and Lotka-Volterra examples but supplies neither a constructive selection algorithm nor a proof that such bases exist and can be learned from data for arbitrary autonomous ODEs without prior knowledge of the flow.
- [§4] §4 (Data-driven approximation via transfer operators): the claimed convergence of the data-driven KvN approximation to the true eigenvalues and eigenfunctions is asserted but not supported by error bounds, consistency proofs, or quantitative numerical verification beyond the guiding examples. Without these, it is unclear whether the method improves upon or merely reformulates existing DMD-type estimators.
minor comments (2)
- Notation for the inner product and measure with respect to which unitarity is claimed should be introduced explicitly at first use rather than left implicit.
- The abstract promises to 'show that the choice of basis functions and domain are crucial'; the corresponding discussion in the main text should be expanded with a counter-example where an ill-chosen basis destroys unitarity.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments correctly identify areas where additional rigor is needed to support the claims on unitarity of projections and convergence of the data-driven approximations. We have revised the manuscript to include clarifications on the conditions for unitarity, a consistency result for the estimator, and enhanced numerical validation. Our responses to the major comments are as follows.
read point-by-point responses
-
Referee: [§3] the statement that projection onto any finite subspace produces a unitary matrix is not accompanied by a general theorem. Unitarity holds only when the subspace is invariant under the KvN operator and the basis is orthonormal with respect to the appropriate measure; the manuscript demonstrates this on the oscillator and Lotka-Volterra examples but supplies neither a constructive selection algorithm nor a proof that such bases exist and can be learned from data for arbitrary autonomous ODEs without prior knowledge of the flow.
Authors: We agree with the referee that unitarity of the finite-dimensional projection requires the subspace to be invariant under the KvN operator and the basis to be orthonormal with respect to the KvN inner product. The manuscript implicitly relies on this by choosing appropriate bases for the examples (e.g., trigonometric polynomials for the oscillators and suitable functions for Lotka-Volterra). We have added a new Proposition 3.1 in Section 3 stating the precise conditions under which the projected matrix is unitary, along with a proof sketch based on the definition of the KvN operator as a unitary group on L2. However, we do not claim a general constructive algorithm for selecting such bases for arbitrary ODEs without some knowledge of the flow or domain; this is acknowledged as a limitation similar to other basis-dependent methods like EDMD. The data-driven part allows approximation of the matrix elements once the basis is fixed. We have expanded the discussion on basis and domain selection in the revised manuscript to make this explicit. revision: partial
-
Referee: [§4] the claimed convergence of the data-driven KvN approximation to the true eigenvalues and eigenfunctions is asserted but not supported by error bounds, consistency proofs, or quantitative numerical verification beyond the guiding examples. Without these, it is unclear whether the method improves upon or merely reformulates existing DMD-type estimators.
Authors: We acknowledge that the original manuscript lacked explicit error bounds and detailed proofs. In the revision, we have added Theorem 4.1 providing a consistency result: under the assumption that the chosen basis spans an invariant subspace and the data is sampled from the invariant measure, the data-driven estimator converges in operator norm to the true projected KvN operator as the number of data points goes to infinity, with a rate depending on the sampling. This is derived by relating the KvN approximation to the transfer operator via the connection mentioned in the paper, and leveraging existing convergence results for EDMD. For quantitative verification, we have included additional figures in Section 5 showing eigenvalue approximation errors versus data size for the damped oscillator, demonstrating convergence. Regarding improvement over DMD: while the estimator shares similarities with DMD/EDMD, the preservation of unitarity (when the basis conditions are met) is a distinguishing feature that enables the quantum circuit representation, which standard DMD does not guarantee for dissipative systems. We have added a paragraph comparing the two approaches. revision: yes
Circularity Check
No load-bearing circularity; minor self-citation only
full rationale
The derivation exploits known relationships between the Koopman-von Neumann operator and classical transfer operators to obtain data-driven approximations, then projects onto finite subspaces to obtain unitary matrices. No equations or fitted quantities are shown that reduce the central claims (unitary representation, quantum-circuit encoding, or eigenvalue approximation) to tautologies or self-referential fits. Basis and domain selection is explicitly flagged as crucial and is illustrated on concrete examples rather than asserted via a general existence theorem that would collapse into the method itself. Self-citations appear but are not invoked as the sole justification for the core approximation steps or uniqueness claims.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The Koopman-von Neumann operator is unitary for any autonomous ODE, even non-Hamiltonian ones.
- domain assumption Finite-dimensional projections of the operator remain unitary when suitable bases and domains are chosen.
Reference graph
Works this paper leans on
-
[1]
P. W. Shor. Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum computer.SIAM Journal on Computing, 26(5):1484–1509, 1997.doi:10.1137/ S0097539795293172. 20
1997
-
[2]
Quantum algorithms: an overview,
A. Montanaro. Quantum algorithms: an overview.npj Quantum Information, 2:15023, 2016. doi:10.1038/npjqi.2015.23
-
[3]
L. K. Grover. A fast quantum mechanical algorithm for database search. InProceedings of the Twenty-Eighth Annual ACM Symposium on Theory of Computing, STOC ’96, pages 212–219, New York, NY, USA, 1996. Association for Computing Machinery.doi:10.1145/237814. 237866
-
[4]
G. Brassard, P. Høyer, M. Mosca, and A. Tapp. Quantum amplitude amplification and estima- tion. InQuantum computation and information, volume 305 ofContemp. Math., pages 53–74. American Mathematical Society, Providence, RI, USA, 2002.doi:10.1090/conm/305/05215
-
[5]
A. W. Harrow, A. Hassidim, and S. Lloyd. Quantum algorithm for linear systems of equations. Physical Review Letters, 103:150502, 2009.doi:10.1103/PhysRevLett.103.150502
-
[6]
A. M. Dalzell, S. McArdle, M. Berta, P. Bienias, C. Chen, A. Gily´ en, C. T. Hann, M. J. Kastoryano, E. T. Khabiboulline, A. Kubica, G. Salton, S. Wang, and F. G. S. L. Brand˜ ao. Quantum Algorithms: A Survey of Applications and End-to-end Complexities. Cambridge University Press, 2025
2025
-
[7]
D. Mauro. On Koopman–von Neumann waves.International Journal of Modern Physics A, 17(09):1301–1325, 2002.doi:10.1142/S0217751X02009680
-
[8]
Y. I. Bogdanov and N. A. Bogdanova. The study of classical dynamical systems using quantum theory. InInternational Conference on Micro-and Nano-Electronics 2014, volume 9440, pages 476–484, 2014
2014
-
[9]
U. Klein. From Koopman–von Neumann theory to quantum theory.Quantum Studies: Math- ematics and Foundations, 5(2):219–227, 2018.doi:10.1007/s40509-017-0113-2
-
[10]
I. Joseph. Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics.Physical Review Research, 2:043102, 2020.doi:10.1103/PhysRevResearch.2. 043102
-
[11]
Y. T. Lin, R. B. Lowrie, D. Aslangil, Y. Suba¸ sı, and A. T. Sornborger. Challenges for quantum computation of nonlinear dynamical systems using linear representations, 2022.arXiv:2202. 02188
2022
-
[12]
Lasota and M
A. Lasota and M. C. Mackey.Chaos, fractals, and noise: Stochastic aspects of dynamics, volume 97 ofApplied Mathematical Sciences. Springer, New York, 2nd edition, 1994
1994
-
[13]
Dellnitz and O
M. Dellnitz and O. Junge. On the approximation of complicated dynamical behavior.SIAM Journal on Numerical Analysis, 36(2):491–515, 1999
1999
-
[14]
I. Mezi´ c. Spectral properties of dynamical systems, model reduction and decompositions. Nonlinear Dynamics, 41(1):309–325, 2005.doi:10.1007/s11071-005-2824-x
-
[15]
M. Budiˇ si´ c, R. Mohr, and I. Mezi´ c. Applied Koopmanism.Chaos: An Interdisciplinary Journal of Nonlinear Science, 22(4), 2012.doi:10.1063/1.4772195
-
[16]
Sch¨ utte and M
C. Sch¨ utte and M. Sarich.Metastability and Markov State Models in Molecular Dynamics: Modeling, Analysis, Algorithmic Approaches. Number 24 in Courant Lecture Notes. American Mathematical Society, 2013
2013
-
[17]
I. Mezi´ c. Analysis of fluid flows via spectral properties of the Koopman operator.Annual Review of Fluid Mechanics, 45:357–378, 2013.doi:10.1146/annurev-fluid-011212-140652. 21
-
[18]
S. Klus, P. Koltai, and C. Sch¨ utte. On the numerical approximation of the Perron–Frobenius and Koopman operator.Journal of Computational Dynamics, 3(1):51–79, 2016.doi:10. 3934/jcd.2016003
2016
-
[20]
I. Mezi´ c. A transfer operator approach to relativistic quantum wavefunction.Journal of Physics A: Mathematical and Theoretical, 56(9):094001, 2023.doi:10.1088/1751-8121/acb675
-
[21]
M. Stengl, P. Gelß, S. Klus, and S. Pokutta. Existence and uniqueness of solutions of the Koopman–von Neumann equation on bounded domains.Journal of Physics A: Mathematical and Theoretical, 57(39):395302, 2024.doi:10.1088/1751-8121/ad6f7d
-
[22]
McCaul, A
G. McCaul, A. Pechen, and D. I. Bondar. Entropy nonconservation and boundary conditions for Hamiltonian dynamical systems.Physical Review E, 99:062121, 2019.doi:10.1103/ PhysRevE.99.062121
2019
-
[23]
F. Grotto. Essential self-adjointness of Liouville operator for 2D Euler point vortices.Journal of Functional Analysis, 279(6):108635, 2020.doi:10.1016/j.jfa.2020.108635
-
[24]
S. Simon, R. Santagati, M. Degroote, N. Moll, M. Streif, and N. Wiebe. Improved precision scaling for simulating coupled quantum-classical dynamics.PRX Quantum, 5:010343, 2024. doi:10.1103/PRXQuantum.5.010343
-
[25]
S. Cochran, J. Stokes, P. Jayakumar, and S. Veerapaneni. An application of continuous- variable gate synthesis to quantum simulation of classical dynamics.AVS Quantum Science, 7(2):023801, 2025.doi:10.1116/5.0234007
-
[26]
I. Novikau and I. Joseph. Quantum algorithm for the advection-diffusion equation and the Koopman–von Neumann approach to nonlinear dynamical systems.Computer Physics Com- munications, 309:109498, 2025.doi:10.1016/j.cpc.2025.109498
-
[27]
No´ e and F
F. No´ e and F. N¨ uske. A variational approach to modeling slow processes in stochastic dynam- ical systems.Multiscale Modeling & Simulation, 11(2):635–655, 2013
2013
-
[28]
M. O. Williams, I. G. Kevrekidis, and C. W. Rowley. A data-driven approximation of the Koopman operator: Extending dynamic mode decomposition.Journal of Nonlinear Science, 25(6):1307–1346, 2015.doi:10.1007/s00332-015-9258-5
-
[29]
M. O. Williams, C. W. Rowley, and I. G. Kevrekidis. A kernel-based method for data-driven Koopman spectral analysis.Journal of Computational Dynamics, 2(2):247–265, 2015.doi: 10.3934/jcd.2015005
-
[30]
A. Mauroy and J. Goncalves. Linear identification of nonlinear systems: A lifting technique based on the Koopman operator. In2016 IEEE 55th Conference on Decision and Control (CDC), pages 6500–6505, 2016.doi:10.1109/CDC.2016.7799269
-
[31]
S. Klus, I. Schuster, and K. Muandet. Eigendecompositions of transfer operators in re- producing kernel Hilbert spaces.Journal of Nonlinear Science, 2019.doi:10.1007/ s00332-019-09574-z
2019
-
[32]
S. Klus, F. N¨ uske, S. Peitz, J.-H. Niemann, C. Clementi, and C. Sch¨ utte. Data-driven ap- proximation of the Koopman generator: Model reduction, system identification, and control. Physica D: Nonlinear Phenomena, 406:132416, 2020.doi:10.1016/j.physd.2020.132416. 22
-
[33]
S. Klus, F. N¨ uske, and B. Hamzi. Kernel-based approximation of the Koopman generator and Schr¨ odinger operator.Entropy, 22(7), 2020.doi:10.3390/e22070722
-
[34]
S. Klus, F. N¨ uske, and S. Peitz. Koopman analysis of quantum systems.Journal of Physics A: Mathematical and Theoretical, 55(31):314002, 2022.doi:10.1088/1751-8121/ac7d22
-
[35]
Meyer, G
K. Meyer, G. Hall, and D. Offin.Introduction to Hamiltonian Dynamical Systems and the N-Body Problem. Number 90 in Applied Mathematical Sciences. Springer, New York, 2008
2008
-
[36]
Kaiser, J
E. Kaiser, J. N. Kutz, and S. L. Brunton. Data-driven discovery of Koopman eigenfunctions for control.Machine Learning: Science and Technology, 2(3):035023, 2021.doi:10.1088/ 2632-2153/abf0f5
2021
-
[37]
S. Jin, N. Liu, and Y. Yu. Quantum simulation of partial differential equations via Schr¨ odingerization.Physical Review Letters, 133:230602, 2024.doi:10.1103/PhysRevLett. 133.230602
-
[38]
Giannakis, A
D. Giannakis, A. Ourmazd, P. Pfeffer, J. Schumacher, and J. Slawinska. Embedding classical dynamics in a quantum computer.Physical Review A, 105:052404, 2022.doi:10.1103/ PhysRevA.105.052404
2022
-
[39]
A. Mauroy and I. Mezi´ c. Global stability analysis using the eigenfunctions of the Koopman operator.IEEE Transactions on Automatic Control, 61(11):3356–3369, 2016.doi:10.1109/ TAC.2016.2518918
-
[40]
Korda and I
M. Korda and I. Mezi´ c. On convergence of Extended Dynamic Mode Decomposition to the Koopman operator.Journal of Nonlinear Science, 28(2):687–710, 2018
2018
-
[41]
J. J. Bramburger and G. Fantuzzi. Auxiliary functions as Koopman observables: Data-driven analysis of dynamical systems via polynomial optimization.Journal of Nonlinear Science, 34(1):8, 2023.doi:10.1007/s00332-023-09990-2
-
[42]
Waldvogel
J. Waldvogel. The period in the Lotka–Volterra system is monotonic.Journal of Mathematical Analysis and Applications, 114(1):178–184, 1986
1986
-
[43]
Y.-A. Ma and H. Qian. A thermodynamic theory of ecology: Helmholtz theorem for Lotka–Volterra equation, extended conservation law, and stochastic predator–prey dynam- ics.Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 471(2183):20150456, 2015.doi:10.1098/rspa.2015.0456
-
[44]
M. J. Colbrook, L. J. Ayton, and M. Sz˝ oke. Residual dynamic mode decomposition: robust and verified Koopmanism.Journal of Fluid Mechanics, 955:A21, 2023.doi:10.1017/jfm. 2022.1052
work page doi:10.1017/jfm 2023
-
[45]
C. Geuzaine and J.-F. Remacle. Gmsh: A 3-d finite element mesh generator with built-in pre- and post-processing facilities.International Journal for Numerical Methods in Engineering, 79(11):1309–1331, 2009.doi:10.1002/nme.2579
-
[46]
I. A. Baratta, J. P. Dean, J. S. Dokken, M. Habera, J. Hale, C. N. Richardson, M. E. Rognes, M. W. Scroggs, N. Sime, and G. N. Wells. DOLFINx: the next generation FEniCS problem solving environment, 2023.doi:10.5281/zenodo.10447666
-
[47]
S. Das and D. Giannakis. Koopman spectra in reproducing kernel Hilbert spaces.Applied and Computational Harmonic Analysis, 49(2):573–607, 2020.doi:https://doi.org/10.1016/j. acha.2020.05.008. 23
work page doi:10.1016/j 2020
-
[48]
F. N¨ uske and S. Klus. Efficient approximation of molecular kinetics using random Fourier features.The Journal of Chemical Physics, 159(7):074105, 2023.doi:10.1063/5.0162619
-
[49]
Q. Li, F. Dietrich, E. M. Bollt, and I. G. Kevrekidis. Extended dynamic mode decomposition with dictionary learning: A data-driven adaptive spectral decomposition of the Koopman operator.Chaos: An Interdisciplinary Journal of Nonlinear Science, 27(10):103111, 2017. doi:10.1063/1.4993854
-
[50]
A. Mardt, L. Pasquali, H. Wu, and F. No´ e. VAMPnets for deep learning of molecular kinetics. Nature Communications, 9, 2018.doi:10.1038/s41467-017-02388-1
- [51]
-
[52]
P´ erez-Hern´ andez, F
G. P´ erez-Hern´ andez, F. Paul, T. Giorgino, G. De Fabritiis, and F. No´ e. Identification of slow molecular order parameters for Markov model construction.The Journal of Chemical Physics, 139(1), 2013
2013
-
[53]
F. N¨ uske, B. G. Keller, G. P´ erez-Hern´ andez, A. S. J. S. Mey, and F. No´ e. Variational approach to molecular kinetics.Journal of Chemical Theory and Computation, 10(4):1739–1752, 2014. doi:10.1021/ct4009156
-
[54]
C. Sch¨ utte, S. Klus, and C. Hartmann. Overcoming the timescale barrier in molecular dy- namics: Transfer operators, variational principles and machine learning.Acta Numerica, 32:517–673, 2023.doi:10.1017/S0962492923000016
-
[55]
N¨ uske, S
F. N¨ uske, S. Peitz, F. Philipp, M. Schaller, and K. Worthmann. Finite-data error bounds for Koopman-based prediction and control.Journal of Nonlinear Science, 33(1):14, 2023
2023
-
[56]
Zhang and E
C. Zhang and E. Zuazua. A quantitative analysis of Koopman operator methods for system identification and predictions.Comptes Rendus M´ ecanique, 351(S1):1–31, 2023.doi:10.5802/ crmeca.138
2023
-
[57]
L. Llamazares-Elias, S. Llamazares-Elias, J. Latz, and S. Klus. Data-driven approximation of Koopman operators and generators: Convergence rates and error bounds, 2024.arXiv: 2405.00539
-
[58]
D. C. Freeman, D. Giannakis, and J. Slawinska. Quantum Mechanics for Closure of Dynamical Systems.Multiscale Modeling & Simulation, 22(1):283–333, 2024.doi:10.1137/22M1514246
-
[59]
Giannakis and M
D. Giannakis and M. Montgomery.Koopman and Transfer Operator Techniques from the Perspective of Quantum Theory, pages 1–35. Springer, Basel, 2026.doi:10.1007/ 978-3-0348-0692-3_125-1. A Comparison with the QECD approach Thequantum embedding of classical dynamics(QECD) approach described in [38, 58, 59] allows for an embedding of classical states and obser...
2026
-
[60]
SinceA 3 =−r 2A, the exponential series for exp(t A) reduces to a combination ofI,A, andA 2, i.e., exp(tA) =I+ sin(rt) r A+ 1−cos(rt) r2 A2 = cos(rt)− sin(rt) r z1 − sin(rt) r z2 − sin(rt) r z3 sin(rt) r z1 1− 1−cos(rt) r2 z2 1 − 1−cos(rt) r2 z1z2 − 1−cos(rt) r2 z1z3 sin(rt) r z2 − 1−cos(rt) r2 z1z2 1− 1−cos(rt) r2 z2 2 − 1−cos(rt) r2 z2z3 sin(r...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.