pith. sign in

arxiv: 2605.19187 · v1 · pith:MDGVPNQAnew · submitted 2026-05-18 · ⚛️ physics.flu-dyn · quant-ph

Unitary discretization of the Koopman-von Neumann equation for quantum simulation of fluid and plasma dynamics

Aleksandar Jemcov , Scott C. Morris This is my paper

Pith reviewed 2026-05-20 06:55 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn quant-ph
keywords Koopman-von Neumann equationquantum simulationfluid dynamicsplasma dynamicsunitary discretizationsummation-by-partsWeyl ordering
0
0 comments X

The pith

A Weyl-ordered discretization of the Koopman-von Neumann equation produces exactly unitary operators on any grid for quantum fluid and plasma simulations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that symmetrizing the Koopman-von Neumann generator in Weyl order and pairing it with summation-by-parts finite differences yields discrete evolution operators whose exponentials remain unitary to machine precision. This unitarity is an algebraic identity that survives arbitrary grid spacing and arbitrary stencil order. The formulation evolves directly in physical amplitude space, avoiding the phase-space doubling that would otherwise impose an uncertainty-principle limit on resolution during time stepping. A single-ancilla Kraus absorbing layer handles boundaries while preserving the overall unitary character. Numerical tests on viscous Navier-Stokes, incompressible Euler, and Hasegawa-Mima triads recover the expected convergence rates with unitarity errors at round-off levels.

Core claim

For real velocity fields the Weyl-ordered Koopman-von Neumann generator admits a unique anti-Hermitian symmetrization; when this generator is discretized by summation-by-parts, the resulting operators satisfy an algebraic identity that guarantees exact discrete unitarity independent of grid resolution and stencil order.

What carries the argument

Weyl-ordered KvN generator symmetrized to anti-Hermitian form for real velocities, discretized via summation-by-parts to enforce discrete anti-Hermiticity.

Load-bearing premise

The velocity fields must be real so that a unique anti-Hermitian symmetrization of the generator exists.

What would settle it

Implementation on a coarse grid using an odd-order stencil that produces an operator whose norm deviates from unity by more than machine epsilon would falsify the algebraic unitarity claim.

Figures

Figures reproduced from arXiv: 2605.19187 by Aleksandar Jemcov, Scott C. Morris.

Figure 1
Figure 1. Figure 1: FIG. 1. Worst-case norm ratio [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Grid convergence for the NS triad ( [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Domain sizing study at fixed grid count for the NS triad ( [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. NS triad ( [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. PDF evolution for the NS triad ( [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Euler triad ( [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Hasegawa–Mima triad ( [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
read the original abstract

The Koopman--von Neumann (KvN) formulation of spectrally truncated fluid and plasma dynamics is considered as a potential approach for quantum computation. The KvN framework embeds the Liouville equation into a Hilbert space with norm-preserving, unitary evolution. Here, we propose a Weyl-ordered KvN generator along with a summation-by-parts discretization, which ensures that the resulting operators are exactly unitary as required for quantum computers. The Weyl-ordered KvN generator is derived as the unique anti-Hermitian operator symmetrization for real velocity fields. The formulation operates directly in the physical amplitude space without phase-space doubling, so the Heisenberg uncertainty principle does not constrain the grid resolution during evolution. This limitation re-enters only at the measurement stage on a quantum computer. Exact discrete unitarity is proved as a purely algebraic identity that holds regardless of grid resolution or stencil order. To manage boundaries, a split-step Kraus absorbing layer is introduced via a Stinespring dilation requiring only one ancilla qubit. Validation on three test cases spanning dissipative and Hamiltonian regimes (a viscous Navier--Stokes triad, an incompressible Euler triad, and a Hasegawa--Mima drift-wave triad) confirms fourth-order convergence and machine-precision unitarity.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proposes a Weyl-ordered Koopman-von Neumann (KvN) generator discretized via summation-by-parts (SBP) finite differences to achieve exact discrete unitarity for quantum simulation of fluid and plasma dynamics. It derives the generator as the unique anti-Hermitian symmetrization for real velocity fields, proves that the resulting operator is exactly anti-Hermitian (hence unitary) as a purely algebraic identity independent of grid resolution and stencil order, introduces a split-step Kraus absorbing layer for open boundaries using a single ancilla qubit via Stinespring dilation, and validates the scheme on three test cases (viscous Navier-Stokes triad, incompressible Euler triad, Hasegawa-Mima drift-wave triad) reporting fourth-order convergence and machine-precision unitarity.

Significance. If the central algebraic unitarity result holds for spatially varying real velocity fields, the work would be significant for quantum algorithms in fluids and plasmas by operating directly in physical amplitude space without phase-space doubling or uncertainty-principle constraints on resolution during evolution. The purely algebraic character of the unitarity claim (if fully demonstrated) and the numerical confirmation across dissipative and Hamiltonian regimes are strengths; the boundary treatment via one-ancilla Kraus layer is a practical contribution.

major comments (2)
  1. [Abstract and derivation of Weyl-ordered KvN generator] Abstract and the section deriving the discrete operator: the claim that exact discrete unitarity holds as a purely algebraic identity 'regardless of grid resolution or stencil order' for variable real velocity fields u(x) requires an explicit demonstration that the SBP discretization of the Weyl-ordered symmetrized product exactly cancels the Hermitian part in the chosen inner product. Standard SBP telescoping works for constant coefficients or periodic domains, but the manuscript must show that interpolation of u onto the difference operator and boundary closures introduce no residual Hermitian component at arbitrary orders.
  2. [Numerical validation on test cases] Numerical validation section: the reported fourth-order convergence and machine-precision unitarity on the three triads must be accompanied by the precise definition of the discrete inner product used to compute the norm and by confirmation that no data points or time steps were excluded from the error tables or plots; without this, it is difficult to assess whether the algebraic identity is verified to machine precision for the variable-coefficient cases.
minor comments (2)
  1. [Discretization section] The manuscript should include explicit low-order matrix representations of the SBP operators and the Weyl symmetrizer for a small grid to illustrate the cancellation mechanism.
  2. [Test case descriptions] Clarify whether the velocity fields in the test cases are exactly divergence-free at the discrete level or if any projection is applied, as this affects the interpretation of the Hamiltonian regime results.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address each major comment point by point below and indicate the revisions planned for the manuscript.

read point-by-point responses

  1. Referee: [Abstract and derivation of Weyl-ordered KvN generator] Abstract and the section deriving the discrete operator: the claim that exact discrete unitarity holds as a purely algebraic identity 'regardless of grid resolution or stencil order' for variable real velocity fields u(x) requires an explicit demonstration that the SBP discretization of the Weyl-ordered symmetrized product exactly cancels the Hermitian part in the chosen inner product. Standard SBP telescoping works for constant coefficients or periodic domains, but the manuscript must show that interpolation of u onto the difference operator and boundary closures introduce no residual Hermitian component at arbitrary orders.

    Authors: We agree that an explicit algebraic demonstration for variable u(x) would strengthen the presentation. The manuscript derives the anti-Hermitian character from the unique Weyl symmetrization combined with the SBP property as an algebraic identity, but we acknowledge that the term-by-term cancellation for interpolated velocity fields and boundary closures was not expanded in full detail. In the revised manuscript we will add an appendix containing this explicit expansion, confirming that the Hermitian part cancels exactly for any grid resolution and stencil order. revision: yes

  2. Referee: [Numerical validation on test cases] Numerical validation section: the reported fourth-order convergence and machine-precision unitarity on the three triads must be accompanied by the precise definition of the discrete inner product used to compute the norm and by confirmation that no data points or time steps were excluded from the error tables or plots; without this, it is difficult to assess whether the algebraic identity is verified to machine precision for the variable-coefficient cases.

    Authors: We thank the referee for this clarification request. The discrete inner product is the standard SBP inner product defined by the diagonal norm matrix with quadrature weights that enforce the summation-by-parts identity. In the revised numerical validation section we will state this definition explicitly. We also confirm that all grid points and all time steps were included in the convergence and unitarity computations; no data were excluded. A statement to this effect will be added to the manuscript. revision: yes

Circularity Check

0 steps flagged

No circularity: unitarity is an algebraic identity from SBP discretization, independent of inputs

full rationale

The paper derives the Weyl-ordered KvN generator as the unique anti-Hermitian symmetrization for real velocity fields and proves exact discrete unitarity as a purely algebraic identity that holds for any grid resolution or stencil order via summation-by-parts discretization. No load-bearing step reduces to a fitted parameter, self-citation chain, or input by construction; the result is presented as following directly from the discrete inner-product properties and symmetrization without reintroducing the target quantity. The assumption of real velocity fields is stated explicitly rather than smuggled in. This is the most common honest finding for a self-contained algebraic derivation against external SBP benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard operator symmetrization techniques and numerical analysis identities from prior literature, with the novel element being their specific combination for quantum compatibility.

axioms (2)
  • domain assumption The KvN formulation embeds the Liouville equation into a Hilbert space with norm-preserving, unitary evolution
    This is the foundational embedding used throughout the work.
  • domain assumption Velocity fields are real
    Invoked to obtain the unique anti-Hermitian symmetrization of the generator.

pith-pipeline@v0.9.0 · 5750 in / 1310 out tokens · 54862 ms · 2026-05-20T06:55:31.000042+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

51 extracted references · 51 canonical work pages · 2 internal anchors

  1. [1]

    For triadic truncations satisfying this resonance condition, the nonzero components of Cjmn reduce to cyclic entries Cj, so T (a) = a1a2a3 P j Cj

    The bound is uniform on [0, T∗) and precludes finite-time blow- up, soT ∗ =∞.□ Fourier modes couple through the quadratic advection nonlinearity only when their wavevectors satisfy a clo- sure condition: the product eikm·x eikn·x = ei(km+kn)·x projects onto a third mode kj only when kj = −(km+kn), so that k1+k2+k3 = 0. For triadic truncations satisfying t...

    work page

  2. [2]

    The cyclic structure of the resulting coupling coeffi- cients ensuresP j Cj = 0 (Appendices B 1 and B 2), so all three cases satisfy the structural condition T≡ 0 of Propo- sition 1. The cases differ in the constraint imposed on the parent PDE (incompressibility for Navier–Stokes and Eu- ler, modified vorticity inversion for Hasegawa–Mima) and in whether ...

    work page

  3. [3]

    Statistical moments were obtained from ⟨g⟩ =P i g(ai)|ψi|2∆V

    ≈ 10−9 initial probability from the nearest boundary, making boundary truncation negligible relative to all other error sources. Statistical moments were obtained from ⟨g⟩ =P i g(ai)|ψi|2∆V . Monte Carlo references were computed from RK4 trajectories sampled from the same Gaussian initial condition. All three triad cases ( N = 3) used M = 40 grid points p...

    work page

  4. [4]

    Phase-space norm bound Proposition 1 predicts two distinct confinement regimes for the three test cases. Figure 1 shows a time history of the worst-case norm ratio maxcorners ∥a(t)∥2/∥a(0)∥2 over all 2N corners of the support box ¯a± 6σ, where ¯aj is the initial condition mean in coordinate j, the most stringent test of confinement. For the NS and HM tria...

    work page

  5. [5]

    Grid convergence The NS triad was considered for convergence analysis for three reasons: it has a compressible phase-space flow (∇·f̸ = 0) unlike the Euler triad, it exercises the full three- dimensional tensor-product grid, and its antisymmetric coupling ensures Plost is negligible at the 6 σ margin, iso- lating resolution error as the sole source of dis...

    work page

  6. [6]

    Enlarging the box reduces the probability truncated at t = 0 but coarsens the grid spacing ∆ a at fixed M, degrading the resolution of the PDF

    Domain sizing Choosing the domain extent Ω a involves two compet- ing effects. Enlarging the box reduces the probability truncated at t = 0 but coarsens the grid spacing ∆ a at fixed M, degrading the resolution of the PDF. In practice the dissipative dynamics contract the PDF rapidly, leav- ing resolution error to dominate at any adequately-sized domain. ...

    work page

  7. [7]

    Contributions from the remaining boundaries are smaller by several orders of magnitude

    ≈ 10−9. Contributions from the remaining boundaries are smaller by several orders of magnitude. The total initial tail probability is therefore Ptail(6)≈ 1 2erfc(6/ √ 2)≈10 −9,(39) where erf(x) = (2/√π) R x 0 e−s2 ds. Since the dissipative dynamics contract the PDF over time, this sets a strict upper bound on the probability that can ever reach the bounda...

    work page

  8. [8]

    The truncation error is, for any smooth functiong: (D(2)g)i =g ′(ai) + (∆a)2 6 g′′′(ai) +O (∆a)4 .(A5)

    Second-order operator The second-order periodic SBP operator has two nonzero weights per row: w±1 =± 1 2∆a , w k = 0 for|k| ̸= 1.(A3) ForM= 6: D(2) = 1 2∆a   0 1 0 0 0−1 −1 0 1 0 0 0 0−1 0 1 0 0 0 0−1 0 1 0 0 0 0−1 0 1 1 0 0 0−1 0   .(A4) Skew-symmetryD (2) + (D(2))T = 0 is manifest. The truncation error is, for any smooth functiong: (D(...

    work page

  9. [9]

    The truncation error is, for any smooth functiong: (D(4)g)i =g ′(ai)− (∆a)4 30 g(5)(ai) +O (∆a)6 .(A8)

    Fourth-order operator The fourth-order periodic SBP operator has four nonzero weights per row: w±1 =± 2 3∆a , w ±2 =∓ 1 12∆a , wk = 0 for|k|>2.(A6) ForM= 6: D(4) = 1 12∆a   0 8−1 0 1−8 −8 0 8−1 0 1 1−8 0 8−1 0 0 1−8 0 8−1 −1 0 1−8 0 8 8−1 0 1−8 0   .(A7) Skew-symmetryD (4) +(D (4))T = 0 is again manifest. The truncation error is, for any...

    work page

  10. [10]

    The general circulant structure of Eq

    Verification of exact skew-symmetry This subsection verifies the conditionD † j = −Dj re- quired by Theorem 1 for the specific operators defined in Appendices A 1 and A 2. The general circulant structure of Eq. (A2) reduces exact skew-symmetry to the anti- symmetry of the weight sequence w−k = −wk. Setting Q(p) = ∆a D(p) for order p∈ { 2, 4} and reading o...

    work page

  11. [11]

    On the full N- dimensional tensor-product grid of size K = M N, the Weyl-ordered generator (20) inherits this sparsity via the Kronecker structure (17)

    Sparsity and qubit cost Both operators are sparse:D (2) has 2 nonzeros per row andD (4) has 4 nonzeros per row. On the full N- dimensional tensor-product grid of size K = M N, the Weyl-ordered generator (20) inherits this sparsity via the Kronecker structure (17). Each termF jDj +D jFj has the same sparsity pattern asD j itself, contributing 2 s nonze- ro...

    work page

  12. [12]

    Incompressibility 17 guarantees a vector potential; the relation ω = −∇2Ψ is the kinematic definition of vorticity

    Incompressible Navier–Stokes and Euler The two-dimensional incompressible Navier–Stokes equations in vorticity–streamfunction form are ∂ω ∂t +J(Ψ, ω) =ν∇ 2ω, ω=−∇ 2Ψ,(B1) where ω is the vorticity, Ψ is the streamfunction, and J(Ψ, ω) = Ψxωy −Ψyωx is the Jacobian. Incompressibility 17 guarantees a vector potential; the relation ω = −∇2Ψ is the kinematic de...

    work page

  13. [13]

    The term −ϕ arises from the ion polarization drift

    Hasegawa–Mima equation The Hasegawa–Mima equation [24] governs electrostatic drift-wave turbulence in magnetized plasmas: ∂ ∂t (∇2ϕ−ϕ) + [ϕ,∇ 2ϕ] =µ∇ 2(∇2ϕ−ϕ),(B6) where ϕ is the electrostatic potential, [ ·,· ] denotes the 2D Poisson bracket, and µ is a collisional damping rate. The term −ϕ arises from the ion polarization drift. The diamagnetic drift te...

    work page

  14. [14]

    Lasota and M

    A. Lasota and M. C. Mackey,Chaos, Fractals, and Noise: Stochastic Aspects of Dynamics, 2nd ed. (Springer, New York, 1994)

    work page 1994

  15. [15]

    Hopf, Statistical hydromechanics and functional calcu- lus, J

    E. Hopf, Statistical hydromechanics and functional calcu- lus, J. Rat. Mech. Anal.1, 87 (1952)

    work page 1952

  16. [16]

    S. B. Pope,Turbulent Flows(Cambridge University Press, Cambridge, UK, 2000)

    work page 2000

  17. [17]

    G. S. Fishman,Monte Carlo: Concepts, Algorithms, and Applications(Springer, New York, 1996)

    work page 1996

  18. [18]

    S. B. Pope, PDF methods for turbulent reactive flows, Prog. Energy Combust. Sci.11, 119 (1985)

    work page 1985

  19. [19]

    D. C. Haworth, Progress in probability density function methods for turbulent reacting flows, Prog. Energy Com- bust. Sci.36, 168 (2010). 18

    work page 2010

  20. [20]

    B. O. Koopman, Hamiltonian systems and transformation in Hilbert space, Proc. Natl. Acad. Sci. U.S.A.17, 315 (1931)

    work page 1931

  21. [21]

    von Neumann, Zur Operatorenmethode in der klassis- chen Mechanik, Ann

    J. von Neumann, Zur Operatorenmethode in der klassis- chen Mechanik, Ann. Math.33, 587 (1932)

    work page 1932

  22. [22]

    Mauro, On Koopman–von Neumann waves, Int

    D. Mauro, On Koopman–von Neumann waves, Int. J. Mod. Phys. A17, 1301 (2002)

    work page 2002

  23. [23]

    Joseph, Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics, Phys

    I. Joseph, Koopman–von Neumann approach to quantum simulation of nonlinear classical dynamics, Phys. Rev. Res.2, 043102 (2020)

    work page 2020

  24. [24]

    Morinishi, T

    Y. Morinishi, T. S. Lund, O. V. Vasilyev, and P. Moin, Fully conservative higher order finite difference schemes for incompressible flow, J. Comput. Phys.143, 90 (1998)

    work page 1998

  25. [25]

    Morinishi, Skew-symmetric form of convective terms and fully conservative finite difference schemes for variable density low-Mach number flows, J

    Y. Morinishi, Skew-symmetric form of convective terms and fully conservative finite difference schemes for variable density low-Mach number flows, J. Comput. Phys.229, 276 (2010)

    work page 2010

  26. [26]

    Brassard, P

    G. Brassard, P. Høyer, M. Mosca, and A. Tapp, Quan- tum amplitude amplification and estimation, inQuantum Computation and Quantum Information, AMS Contem- porary Mathematics, Vol. 305 (American Mathematical Society, Providence, RI, 2002) pp. 53–74

    work page 2002

  27. [27]

    Novikau and I

    I. Novikau and I. Joseph, Quantum algorithm for the advection-diffusion equation and the Koopman–von Neu- mann approach to nonlinear dynamical systems, Comput. Phys. Commun.309, 109498 (2025)

    work page 2025

  28. [28]

    An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

    I. Novikau and I. Joseph, An efficient explicit implementa- tion of a near-optimal quantum algorithm for simulating linear dissipative differential equations, arXiv preprint (2025), arXiv:2501.11146

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  29. [29]

    S. Jin, N. Liu, and Y. Yu, Quantum simulation of partial differential equations: Applications and detailed analysis, Phys. Rev. A108, 032603 (2023)

    work page 2023

  30. [30]

    S. Jin, N. Liu, and Y. Yu, Quantum simulation of partial differential equations via Schr¨ odingerization, Phys. Rev. Lett.133, 230602 (2024)

    work page 2024

  31. [31]

    J.-P. Liu, H. Ø. Kolden, H. K. Krovi, N. F. Loureiro, K. Trivisa, and A. M. Childs, Efficient quantum algorithm for dissipative nonlinear differential equations, Proc. Natl. Acad. Sci. U.S.A.118, e2026805118 (2021)

    work page 2021

  32. [32]

    Engel, G

    A. Engel, G. Smith, and S. E. Parker, Quantum algorithm for the Vlasov equation, Phys. Rev. A100, 062315 (2019)

    work page 2019

  33. [33]

    Giannakis, A

    D. Giannakis, A. Ourmazd, P. Pfeffer, J. Schumacher, and J. Slawinska, Embedding classical dynamics in a quantum computer, Phys. Rev. A105, 052404 (2022)

    work page 2022

  34. [34]

    Waleffe, The nature of triad interactions in homoge- neous turbulence, Phys

    F. Waleffe, The nature of triad interactions in homoge- neous turbulence, Phys. Fluids A4, 350 (1992)

    work page 1992

  35. [35]

    Holmes, J

    P. Holmes, J. L. Lumley, G. Berkooz, and C. W. Rowley, Turbulence, Coherent Structures, Dynamical Systems and Symmetry, 2nd ed. (Cambridge University Press, Cam- bridge, UK, 2012)

    work page 2012

  36. [36]

    A. D. D. Craik,Wave Interactions and Fluid Flows(Cam- bridge University Press, Cambridge, UK, 1985)

    work page 1985

  37. [37]

    Hasegawa and K

    A. Hasegawa and K. Mima, Pseudo-three-dimensional tur- bulence in magnetized nonuniform plasma, Phys. Fluids 21, 87 (1978)

    work page 1978

  38. [38]

    Weyl, Quantenmechanik und Gruppentheorie, Z

    H. Weyl, Quantenmechanik und Gruppentheorie, Z. Phys. 46, 1 (1927)

    work page 1927

  39. [39]

    Kreiss and G

    H.-O. Kreiss and G. Scherer, Finite element and finite difference methods for hyperbolic partial differential equa- tions, inMathematical Aspects of Finite Elements in Par- tial Differential Equations(Academic Press, 1974) pp. 195–212

    work page 1974

  40. [40]

    Strand, Summation by parts for finite difference ap- proximations for d/dx, J

    B. Strand, Summation by parts for finite difference ap- proximations for d/dx, J. Comput. Phys.110, 47 (1994)

    work page 1994

  41. [41]

    Sv¨ ard and J

    M. Sv¨ ard and J. Nordstr¨ om, Review of summation-by- parts schemes for initial-boundary-value problems, J. Comput. Phys.268, 17 (2014)

    work page 2014

  42. [42]

    M. A. Nielsen and I. L. Chuang,Quantum Computation and Quantum Information, 10th ed. (Cambridge Univer- sity Press, Cambridge, UK, 2010)

    work page 2010

  43. [43]

    T. D. Lee, On some statistical properties of hydrodynam- ical and magneto-hydrodynamical fields, Q. Appl. Math. 10, 69 (1952)

    work page 1952

  44. [44]

    R. H. Kraichnan, The structure of isotropic turbulence at very high Reynolds numbers, J. Fluid Mech.5, 497 (1959)

    work page 1959

  45. [45]

    H. K. Khalil,Nonlinear Systems, 3rd ed. (Prentice Hall, Upper Saddle River, NJ, 2002)

    work page 2002

  46. [46]

    A. H. Al-Mohy and N. J. Higham, Computing the action of the matrix exponential, with an application to exponential integrators, SIAM J. Sci. Comput.33, 488 (2011)

    work page 2011

  47. [47]

    W. F. Stinespring, Positive functions on C*-algebras, Proc. Am. Math. Soc.6, 211 (1955)

    work page 1955

  48. [48]

    D. W. Berry, A. M. Childs, R. Cleve, R. Kothari, and R. D. Somma, Simulating Hamiltonian dynamics with a truncated Taylor series, Phys. Rev. Lett.114, 090502 (2015)

    work page 2015

  49. [49]

    G. H. Low and I. L. Chuang, Optimal Hamiltonian sim- ulation by quantum signal processing, Phys. Rev. Lett. 118, 010501 (2017)

    work page 2017

  50. [50]

    H¨ aner, M

    T. H¨ aner, M. Roetteler, and K. M. Svore, Optimizing quantum circuits for arithmetic, Quantum Inf. Comput. 18, 673 (2018)

    work page 2018

  51. [51]

    L. K. Grover and T. Rudolph, Creating superpositions that correspond to efficiently integrable probability distri- butions, arXiv preprint (2002), arXiv:quant-ph/0208112

    work page internal anchor Pith review Pith/arXiv arXiv 2002