Koopman wavefunctions and Clebsch variables in Vlasov-Maxwell kinetic theory

Motivated by recent discussions on the possible role of quantum computation in plasma simulations, here we present different approaches to Koopman's Hilbert-space formulation of classical mechanics in the context of Vlasov-Maxwell kinetic theory. The celebrated Koopman-von Neumann construction is provided with two different Hamiltonian structures: one is canonical and recovers the usual Clebsch representation of the Vlasov density, the other is noncanonical and appears to overcome certain issues emerging in the canonical formalism. Furthermore, the canonical structure is restored for a variant of the Koopman-von Neumann construction that carries a different phase dynamics. Going back to van Hove's prequantum theory, the corresponding Koopman-van Hove equation provides an alternative Clebsch representation which is then coupled to the electromagnetic fields. Finally, the role of gauge transformations in the new context is discussed in detail.


Introduction
Ongoing discussions [43] on the potential of quantum information and computation in plasma physics have recently led to exploiting Hilbert-space approaches in the numerical simulation of magnetized plasmas [14,15,16,29], of the Navier-Stokes equations [20], and of arbitrary non-Hamiltonian systems of equations [29,33]. In particular, recent work [29] has emphasized the role of Koopman wavefunctions in classical dynamics while their usage in describing hybrid quantum-classical system was presented in [5,6,21,24,45,28]. In this paper we want to show how these enter in the variational and Hamiltonian formulation of Vlasov kinetic theory.
In the context of the classical Liouville equation ∂ t f = {H, f }, the idea of a Koopman wavefunction emerges naturally from the fact that the phase-space density is positive-definite, so that the relation f (x, p) = |Ψ(x, p)| 2 induces a Hilbert-space description of classical mechanics hinging on the Koopman-von Neumann (KvN) equation [31,38,49] where and {·, ·} denotes the canonical Poisson bracket. Here, the Liouvillian operator L H is selfadjoint so that the Koopman wavefunction undergoes unitary dynamics, thereby leading to a quantum analogy. While the Koopman propagator exp(−i −1 L H ) (Koopman operator) often appears in the dynamical systems literature [7,23], the role of Koopman wavefunctions is being recognized in physics only recently. Over the decades, prominent authors [3,10,12,46] discussed Koopman wavefunctions without referring to Koopman's original work. However, the interest in Koopman wavefunctions has been recently revived by their role in the quantumclassical divide [4]. One of the important motivations for the introduction of the Koopman formalism is to develop a rigorous Hilbert space theory for the approximation of the particle distribution function (PDF). For example, it is well known that the linearized Vlasov-Maxwell equation possesses Case-van Kampen eigenmodes and one might wish to understand how an expansion in eigenmodes, or some other complete set of eigenfunctions, converges to a solution of the nonlinear Vlasov-Maxwell equations. This point of view was taken by Koopman and von Neumann in their development of ergodic theory.
Yet another important application of the Koopman approach is the development of quantum algorithms for simulating nonlinear classical dynamical systems [29]. In principle, quantum computers can store and process an exponentially large amount of information efficiently and can perform certain calculations, such as the "quantum" Fourier transform, very efficiently [41]. While quantum computers hold great promise, they can only perform linear unitary operations on normalized wavefunctions. A number of authors [1,10,29] have reformulated nonlinear non-Hamiltonian dynamics as a linear unitary evolution of functions on phase space. Thus, as shown in Ref. [29], quantum Hamiltonian simulation algorithms for the Koopman evolution operator can achieve up to an exponential speedup over an Eulerian discretization of the PDF and using amplitude estimation of physical observables leads to up to a quadratic speedup over Lagrangian discretizations based on Monte Carlo and particle-in-cell (PIC).
Depending on the choice of phase, Koopman wavefunctions are related to the Clebsch representation of the Vlasov distribution [39], which may be of more common knowledge in the plasma physics community. However, as discussed in Section 2.1, the complete identification between Koopman wavefunctions and Clebsch variables may require the Koopman phase to be singular at some points or regions in phase space [29]. In turn, the phase is irrelevant in the Koopman-von Neumann theory and one should be able to set it to zero without affecting the general formulation of the theory. In order to overcome this possible issue, we develop an alternative noncanonical formulation that indeed allows for a zero phase. As shown in Section 2.2, this formulation follows from a suitable adaptation of the Euler-Poincaré variational principle for standard Vlasov dynamics [8,44,47]. In this case, the relation to standard Clebsch variables is replaced a momentum map structure that is discussed in Section 2.3. Then, Section 2.4 applies this construction to the Vlasov-Maxwell system.
In addition, a variant of the Koopman-von Neumann construction called the Koopmanvan Hove equation (KvH) will also be illustrated in Section 3. As shown in Section 3.1, in this variant the Koopman-von Neumann equation emerges as the amplitude equation, while the evolution of the phase is given by the classical action, i.e. the integral of the Lagrangian function, following Feynman's general prescription. The KvH formulation is based on a slight extension of standard canonical transformations first appeared in van Hove's thesis [48] and briefly discussed in Section 3.2. As the KvH construction is based on a canonical Hamiltonian structure, its relation to the Vlasov density comprises an extension of the standard Clebsch representation, as presented in Section 3.3. Then, Section 3.4 couples the KvH equation to the evolution of the electromagnetic fields and presents the resulting Hamiltonian structure. Finally, Section 3.5 discusses the role of gauge transformations in the KvH context.

Koopman-von Neumann theory
In this section, we present two alternative formulations of the Koopman-von Neumann (KvN) theory. The first formulation is canonical. In this case, the Hamiltonian functional depends on the phase and vanishes if the latter is initially set to zero. Indeed, this Hamiltonian generally differs from the expression of the total energy obtained by Koopman's prescription f = |Ψ| 2 . This ambiguity may be overcome by enforcing a specific constraint on the phase [29] thereby leading to branch-cut singularities. We show that these intricacies can be removed by resorting to an alternative noncanonical formulation that is entirely built upon the prescription f = |Ψ| 2 , so that the Hamiltonian functional is independent of the phase which can then be set to zero. In order to simplify the notation, most of our discussion will treat the case of a single degree of freedom, so that the phase space is two-dimensional. The generalization to multiple degrees of freedom is straightforward, and will be used when discussing the Vlasov-Maxwell system.

Clebsch variables and canonical structure
As first discussed in [13], the Koopman-von Neumann (KvN) equation is Hamiltonian with the canonical Hamiltonian structure arising from the following Dirac-Frenkel variational principle [19]: Here, Ψ is a square-integrable complex function, A, B = Re A|B is the real-valued pairing associated to the standard L 2 inner product A|B =´A(x, p) * B(x, p) dxdp, and we have Im A|B = iA, B . As discussed in [5,29], this Hamiltonian structure readily identifies a Hamiltonian functional given by This is accompanied by the usual Poisson bracket from standard quantum mechanics [9,18]: where we have introduced the double-bracket notation {{·, ·}} to distinguish from the canonical Poisson bracket. Then, the Madelung transform Ψ = √ De iS/ expresses the Koopman wavefunction in terms of the density D and the phase S, thereby leading to the celebrated Clebsch representation of the real-valued Vlasov density [26,35,39], that is Here, the classical Liouville equation ∂ t f = {H, f } follows by combining the Jacobi identity with the relations which arise from Eq. (1). For example, this formulation appeared in [42], where the Clebsch representation was used together with Koopman's original prescription f = |Ψ| 2 in the context of the Vlasov-Poisson system for electrostatic plasmas. In symplectic geometry, the map Ψ → Im{Ψ * , Ψ} identifies a momentum map for the action of canonical transformations on the Hilbert space H of Koopman wavefunctions, which is endowed with the canonical Poisson structure in Eq. (3). See Section 2.3 later on for the definition of momentum maps.
This picture leads us to two following important observations. First, the relation (4) differs from Koopman's prescription f = |Ψ| 2 for the phase-space density. Second, the Clebsch representation (4) does not generally identify a probability density, since´{D, S} dxdp = 0 whenever D and S are differentiable continuous functions. In [29], a solution to this apparent issue was provided by selecting a singular phase ensuring the following consistency condition: Due to the equations of motion (5), if this relation holds true as an initial condition, then it will hold for all time. In this case, the expression of the phase is necessarily singular and this is strictly necessary to produce the boundary terms generally required for the Clebsch representation (4) to hold true. Then, in this context, we realize that the KvN construction may be envisioned as a special case of Clebsch representation of the Vlasov density [29]. Notice that this representation requires a non-vanishing phase, which again contrasts with the phase-invariance of the original prescription f = |Ψ| 2 . Indeed, while one might wish to set the phase to zero and identify the Koopman wavefunction with a real-valued amplitude, this possibility is excluded by the Clebsch representation, which requires a nonzero phase by construction.
Whether or not the requirement of a nonzero, possibly singular phase is seen as a possible issue, it may still be desirable to develop a Hamiltonian structure of KvN theory that does not necessitate the presence of a nontrivial phase. In view of this, we will proceed by considering an alternative noncanonical Hamiltonian structure of Koopman-von Neumann theory, which appears to overcome this potential ambiguity.

Remark 2.1 (Guillemin-Sternberg collectivization)
We remark that, in both the canonical and the noncanonical case, Koopman wavefunctions provide a representation of the Liouville (Vlasov) density whose dynamics is then governed by the usual Lie-Poisson bracket on the Poisson algebra of phase-space functions [36,37,40]. The procedure taking the Koopman Hamiltonian structure into the Vlasov Lie-Poisson structure is an example of collectivization by a momentum map [25], which serves as the unifying geometric framework in this work; see Section 2.3. Specifically, the canonical action of a Lie group on a Poisson manifold induces a momentum map generalizing Noether's conserved quantity occurring in the particular case of a symmetry group. Then, when a Hamiltonian function(al) can be entirely written in terms of this momentum map, the Hamiltonian is called 'collective', borrowing the terminology from nuclear physics [25]. Here, we derive collective Hamiltonians for a series of models in Vlasov dynamics for which the Lie group is given by canonical transformations (or variants thereof, as in the case of Section 3.2). In this process, different Clebsch-type representations emerge from different actions on the Hilbert space of Koopman wavefunctions. Alternatively, depending on convenience, one may prefer to keep working with the Vlasov density itself.

Euler-Poincaré variational formulation
In this section, we shall start by presenting the Euler-Poincaré formulation of Koopman-von Neumann theory. We start from the Euler-Poincaré variational principle for the classical Liouville equation: As discussed extensively in [8,44,47], this variational principle arises from a symmetry reduction from the Lagrangian-path formulation in phase-space to the corresponding Eulerianvariable description. Indeed, if we introduce the notation J jℓ = {z j , z ℓ } and z 0 = (x 0 , p 0 ) is the phase-space label coordinate, then the equation of motionη = J∇H(η) of the Lagrangian trajectory η(z 0 , t) = (η q (z 0 , t), η p (z 0 , t)) follows from the variational principle δ´t 2 t 1´f 0 (η pηq − H(η q , η p )) dx 0 dp 0 dt = 0, where f 0 (z 0 ) is the reference Liouville density. In turn, the latter variational principle leads to the Eulerian counterpart (7) by a mere change of variables. Upon denoting by z = (x, p) the Eulerian phase-space coordinate, the Eulerian phase-space density evolves according to the Lagrange-to-Euler map: Here, det ∇η denotes the Jacobian determinant. In addition, we define the Eulerian phasespace vector field X = (u, σ) such that In this setting, we also define the infinitesimal displacement Ξ(z, t) = δη(z 0 , t) z 0 =η −1 (z,t) so that the Euler-Poincaré variations read which are obtained from the defining relations (8) and (9). Here, Ξ is an arbitrary displacement vanishing at the endpoints, i.e. Ξ(t 1 ) = Ξ(t 2 ) = 0. Then, the Euler-Poincaré variational principle (7) yields X = J∇H, which is accompanied by the continuity equation ∂ t f + div(f X) = 0 following from (8), thereby returning ∂ t f + {f, H} = 0.
The Euler-Poincaré approach for the Koopman-von Neumann theory arises immediately from the construction above. Indeed, the Koopman description f = |Ψ| 2 immediately leads to which is the standard Lagrange-to-Euler map for half-densities [2]. This evolution law leads to the relations One can introduce the operators Λ ℓ = −i ∂ ℓ to write the above as Schrödinger-like equations: where [·, ·] + denotes the anticommutator, X ℓ are simply the multiplicative operators identified by the components of X, and likewise for Ξ ℓ . Notice that here we are not asking for the vector field X to be Hamiltonian: as we shall see, this property will result as a consequence of the variational principle under consideration. In more generality, the evolution law (11) and its equation of motion in (12) also appear in other contexts involving non-Hamiltonian dynamical systems [1,29].
Since we are interested in the Hamiltonian structure, we shall write the Euler-Poincaré variational principle in terms of an arbitrary Hamiltonian functional h(Ψ) as follows: where we recall (9). Here, the arbitrary functional h(Ψ) generally depends on Ψ and its conjugate Ψ * . Upon combining (12) with the first equation in (10), the Euler-Poincaré variational principle yields Here, we have used the functional derivative notation so that δh := δh/δΨ, δΨ and, following Ref. [27], we have defined the diamond operator ⋄ for compactness of notation. Then, the wavefunction equation for an arbitrary Hamiltonian is obtained by substituting (14) in the first equation of (12), which leads to quite a cumbersome explicit form. However, in the case of Koopman-von Neumann theory, this simplifies substantially. Indeed, if the Hamiltonian functional depends only on f = |Ψ| 2 , then the chain rule relation δh/δΨ = 2(δh/δf )Ψ transforms (14) to the Hamiltonian vector field X = J∇ δh δf =: X δh/δf , thereby recovering the Koopman-von Neumann equation (1) in the general form In this case, the Lagrangian trajectory is a canonical transformation so that the propagator (11) reduces to Ψ(z, t) Notice that, while here we have used canonical coordinates associated to the canonical symplectic form dx∧dp, the possibility of a noncanonical structure does not pose any difficulty. Indeed, one would simply redefine canonical transformations as symplectomorphisms, that is smooth invertible transformations preserving the noncanonical symplectic structure. Then, the Liouvillian operator L H = {i H, } would naturally involve the corresponding noncanonical Poisson bracket.

Non-canonical Poisson bracket and momentum map structure
We realize that equation (15) arises naturally from the Euler-Poincaré formulation of the Vlasov equation, while its underlying structure is quite involved and this reflects in an intricate noncanonical Hamiltonian structure. Given an arbitrary functional k(Ψ), the usual relationk = {{k, h}} leads in this case to which reduces to the Lie-Poisson bracket {{k, h}}(f ) =ˆf δk δf , δh δf dx dp (17) in the case of functionals depending only on f = |Ψ| 2 . To obtain (16), we have used the relation Ψ ⋄ Υ, W = Υ, W · ∇Ψ + div(W)Ψ/2 for any Koopman wavefunctions Ψ and Υ, and any vector field W. Indeed, this relation leads tȯ where the last equality follows from (14). While expanding the terms in the noncanonical structure (16) does not lead to much insight, the present formulation has the advantage of restoring the use of the relation f = |Ψ| 2 everywhere in the problem. For example, the standard expression of the total energy now coincides unambiguously with the Hamiltonian functional, which no longer vanishes upon setting the phase to zero. We conclude this section by showing that the map Ψ → |Ψ| 2 is a momentum map [34] associated to the following (left) group action of canonical transformations: Ψ → Ψ • η −1 , where • denotes composition of functions. By the arguments in Remark 2.1, this ensures that the Poisson bracket (16) consistently reduces to (17) for any two functionals k and h of the type k(Ψ) =k(|Ψ| 2 ). Given the action of a Lie group G on a manifold M with Poisson structure {{·, ·}}, a momentum map J : M → g * is defined as for any function k on M and any element ξ of the Lie algebra g of G.
Here, g * denotes the dual space of g, ·, · is the duality pairing, d is the differential on M, and ξ M denotes the infinitesimal generator of the G−action on M. If G is a symmetry of the Hamiltonian, then J (x) is conserved in time by Noether's theorem. The Lie algebra of the group of canonical transformations is the space of Hamiltonian vector fields X = X ξ with ξ = ξ(x, p). Up to the addition of irrelevant constants, this space can be identified with the space of phase-space functions endowed with the canonical Poisson bracket. Then, the corresponding infinitesimal action on the Koopman Hilbert space H = L 2 (R 2 ) is given by ξ H (Ψ) = −X ξ · ∇Ψ = {ξ, Ψ}. Then, upon using the L 2 −pairing, we compute thereby proving (18). This picture extends the usual Clebsch representation treated in Section 2.1 to the noncanonical case. By proceeding analogously, one proves that the quantitieŝ comprise the plasma-to-fluid momentum map [37]. Specifically, these are momentum maps for the action of momentum translations (x, p) → (x, p − dϕ(x)) and configuration-space diffeomorphisms (x, p) → (η(x), −pdη(x)), respectively.

The KvN-Maxwell system
In this section, we want to apply the previous noncanonical formulation to the Maxwell-Vlasov system of magnetized plasmas. The canonical treatment is found in [42] for the case of the electrostatic limit. In the full electromagnetic case, gauge-invariance naturally acquires a prominent role. In particular, gauge freedom manifests itself in the Lagrangian formulation of Maxwell's equations through the fact that the Lagrangian is independent of ∂ t Φ, and, hence, variations with respect to Φ serve to enforce Gauss' law as a constraint. In the absence of sources, the Maxwell Lagrangian is written as is the Maxwell Hamiltonian expressed in terms of the standard L 2 −norm. If we express the Vlasov density in noncanonical coordinates, we can construct the variational principle by an immediate extension of Eq. (13). We write where is the Hamiltonian, while Gauss' law is enforced in (20) by the Lagrange multiplier Φ. Using (12) and the first equation in (10), arbitrary variations of the fields lead to Here, the subscript v denotes the velocity components of vector-valued quantities in phasespace. These equations reveal an intricate Hamiltonian structure whose Poisson bracket can be easily obtained by applying the usual relationk = {{k, h}} upon expandingk = δk/δΨ, ∂ t Ψ + δk/δE, ∂ t E + δk/δB, ∂ t B and using the equations of motion for an arbitrary gauge-invariant Hamiltonian h. Indeed, upon taking the curl of the first equation in (21) and using δh/δA = ∇ × δh/δB, one obtains the following structure: Then, one can simply evaluate which is accompanied by Faraday's and Ampére's law Here, the first equation is obtained by taking the curl of the first equation in (21).

Koopman-van Hove theory
As we saw in previous sections, the KvN phase remains constant along the phase-space Lagrangian trajectories; see the second equation in (5). Indeed, in the KvN construction the phase is entirely irrelevant and plays the role of a gauge freedom in the relation f = |Ψ| 2 . However, in classical mechanics one usually relates the classical phase to the Lagrangian function. Specifically, in Hamilton-Jacobi theory the phase is a function in configuration space and this function is given by the classical action integral. We observe that this is different from the Koopman phase, which instead is defined on phase-space. As we shall see, the Koopman-van Hove construction combines the theory of Koopman wavefunctions with Feynman's prescription of a phase expressed in terms of the Lagrangian. Again, to keep the notation as simple as possible, the phase space will be two-dimensional in most of our discussion. The extension to six dimensions is straightforward and will be used when coupling the KvH equation to the electromagnetic fields.

The Koopman-van Hove equation
An alternative theory of classical mechanics based on Koopman wavefunctions goes back to van Hove's thesis [48], where canonical transformations were extended to include phase factors. As later shown by Kostant [32], in this setting the phase function is again identified with the action integral, which is now defined in terms of the phase-space Lagrangian. In this setting, the KvN equation (1) becomes the Koopman-van Hove equation Here L H = ∂ x Hλ p − ∂ p Hλ x while L is the expression of the phase-space Lagrangian, which now identifies a phase term. The information contained in this equation may be unfolded by applying the Madelung transform Ψ = √ De iS/ , which leads to The second equation can be formally solved in terms of the Lagrangian trajectories as [21,22] S(z, t) =ˆt thereby showing how the phase evolution emerges from the integral of the Lagrangian, in analogy to Feynman's path-integral formulation of quantum mechanics. For a discussion of the relation between equation (24) and the Hamilton-Jacobi equation, see de Gosson's work in [11].
Remark 3.1 (Relation to hybrid quantum-classical dynamics) Notice that, for Hamiltonians of the type H = T + V , enforcing ∂ p Ψ = 0 and replacing p → −i ∂ x takes (22) into the quantum Schrödinger equation, thereby justifying the early name prequantum Schrödinger equation [32]. Indeed, equation (22) first appeared within the context of prequantization theory [30,48] and it has remained pretty unknown over the decades. Recently, it was recognized how this equation may actually lead to a consistent theory of quantum-classical coupling [5,21,24], where the phase plays a crucial role. Partly inspired by Kirillov [30], the authors of [5,21] called equation (

Geometric setting
While canonical transformations are enough to characterize the evolution of Koopman wavefunctions in KvN theory, the presence of the phase in the KvH formalism requires extending the KvN picture. A more detailed summary of the geometric setting of KvH theory is found in [17,21,22]. First, one introduces the U(1)−bundle R 2 × U(1), where R 2 is the Euclidean two-dimensional phase-space and U(1) is the group of complex phase factors. Gauge transformations are identified, as usual, with local phase factors so that the KvH wavefunction evolves according to compositions of gauge transformations and canonical transformations, that is Notice that, unlike (11), here we are restricting the Lagrangian trajectory to identify a canonical transformation at all times. The geometric characterization of the phase factor in (25) needs further discussion. Specifically, the relation between the phase factor and the phase-space Lagrangian emerges as follows. Upon defining the gauge connection A = pdx, it is well known [34] that d(η * A − A) = 0, where d is the exterior differential and η * A(z) := A ℓ (η(z))∇η ℓ (z) is the standard pullback of the connection one-form A by the canonical transformation η. Then, the phase factor in (25) is defined via Indeed, upon using the Lie derivative theorem d(η * A)/dt = η * £ X H A, we notice that Cartan's magic formula takes the time derivative of (26) into the form (z 0 , t)). Then, up to a time-dependent frequency ω(t), equation (24) follows from the defining relation S(z, t) : The evolution law (25) together with the definition (26) represents the KvH analogue of (11) from KvN theory. The propagator Ψ 0 (z) → Ψ(z, t) was first devised by van-Hove [48] and identifies a unitary transformation called a van Hove transformation in [21,22]. Given a canonical transformation η, the van Hove transformation reads Without going much into the details, here we shall simply point out that van Hove transformations possess a group structure whose Lie algebra is identified with the space of Hamiltonian functions endowed with the Lie bracket given by the canonical Poisson bracket. In turn, the self-adjoint operator in (22) identifies the infinitesimal generator −i −1 L H of van Hove transformations. Since KvH theory emerged historically in prequantization theory, here we shall keep the standard nomenclature by calling (28) the prequantum operator.

The phase-space density
So far, nothing has been said about the relation between KvH wavefunctions and the classical phase-space density. In principle, one could insist on following Koopman's original prescription f = |Ψ| 2 . However, as we will see shortly, this step poses questions similar to those arising in Section 2.1. At present, the only Hamiltonian structure available for the KvH equation (22) is given by the canonical bracket (3), which is accompanied by the Hamiltonian functional where we notice that ∂ p (p|Ψ| 2 ) = div(JA|Ψ| 2 ). If we follow Koopman's original prescription of using f = |Ψ| 2 for the phase-space density, then the Hamiltonian functional does not generally coincide with the total energy of the system. As shown in [5,21,22], the expression in parenthesis in (29) identifies the momentum map for the (left) unitary representation (27) of van Hove transformations on the Koopman Hilbert space H = L 2 (R 2 ). Consequently, by the arguments in Remark 2.1, the identification leads to the usual Liouville equation ∂ t f = {f, H}. While it may be objected that f is not positive definite, the flow of f preserves the sign of the initial condition thereby eliminating this apparent problem. Nevertheless, the identification (30) represents a change of perspective from the conventional Koopman prescription f = |Ψ| 2 in that the KvH phase S enters the expression of the phase-space density. However, while the f given in (30) comprises the entire physical information, we emphasize that in the present formalism the single terms in (30) do not possess any physical meaning despite the fact that both the first and the sum of the last two obey the Liouville equation. In particular, this representation of the phase-space density identifies an alternative Clebsch representation extending the usual case given by the last term in (30); see Section 2.1. Further discussions are found in [21,22]. A point of relevance for later purpose is that the momentum map f (Ψ) in (30) is covariant (or equivariant) with respect to canonical transformations. Indeed, using the notation of Eq. (27), we have The details of this property can be found in eg. [21].

Remark 3.2 (Constraints and phase singularities)
Notice that here we can follow the same arguments as in Section 2.1 in order to enforce the relation f = |Ψ| 2 as a specific constraint. Indeed, as hinted in [21], one may be tempted to write Ψ = √ De iS/ and choose the phase S so that ∂ p (p|Ψ| 2 ) + Im{Ψ * , Ψ} = div (DJ(∇S − A)) = 0.
As proven in Ref. [29], if this relation holds for the initial condition, then it will hold for all time. Before concluding this Section, we notice the expressions of the first two momentŝ f dp =ˆ|Ψ| 2 + Imˆ{Ψ * , Ψ} dp ,ˆpf dp = Imˆp{Ψ * , Ψ} dp , which will play a crucial role in the coupling to electromagnetic fields, as shown in the next section. In more generality, the relation p m f dp =ˆp m (1 − m)|Ψ| 2 + Im{Ψ * , Ψ} dp provides an alternative representation of Vlasov moments in terms of canonical variables.

KvH-Maxwell system and its Hamiltonian structure
In this section, we apply the KvH formalism to the Vlasov kinetic theory of magnetized plasmas.
In particular, we are interested in the Hamiltonian structure for the system comprising the electromagnetic component as well as the KvH wavefunction Ψ(x, v) expressed in terms of the velocity variable v = (p − qA)/m. As customary in the geometric approach to the Maxwell-Vlasov system, we start in terms of canonical variables and write the action principle for an arbitrary gauge where Here, h Max (E, A) is the standard Maxwell Hamiltonian (19) while h KvH (Ψ) = Ψ, L H Ψ involves the prequantum operator (28), so that Ψ, L H Ψ = f, H with f given in (30). Similarly Ψ, L Φ Ψ =´f Φ d 3 xd 3 p and Φ plays again the role of a Lagrange multiplier enforcing Gauss' law as in (20). Notice that here one may choose to include the Φ−terms in the Hamiltonian, which would then become a Routhian. However, in order to obtain an explicit Hamiltonian structure comprising a Poisson bracket, it is customary to fix the Hamiltonian gauge Φ = 0 in (32) so that the resulting Hamilton's equations yield We recall that the angle brackets denote the standard L 2 −pairing. Then, the KvH equation i ∂ t Ψ = L H Ψ is accompanied by ∂ t A = −E and Ampère's law in the form An immediate way of obtaining the Hamiltonian structure in terms of Ψ(x, v) is given by a direct change of coordinates. Here, we notice the convenient abuse in denoting the Koopman wavefunction of both canonical and noncanonical coordinates by the same symbol Ψ. The coordinate change (x, p) → (x, v) = (x, (p − qA)/m) corresponds to the replacement which in turn transforms (34) into Notice that, under the same change of variables (x, p) → (x, v), the prequantum operator (28) becomes Unless otherwise specified, in the remainder of this section we shall restrict to the case H = mv 2 /2, as in (33). Also, the expression (30) of the phase-space density becomes In conclusion, upon writing (λ which is accompanied by ∂ t A = −E and Ampère's law in the form More explicitly, one has In order to derive these relations, one makes repeated use of integration by parts in combination with (36), (37), and Here, the last equality follows from Im´{Ψ, and standard vector algebra. Notice that we have the relation which follows directly from the KvH equation i (∂ t Ψ + qm −1 E · ∇ v )Ψ = L H Ψ for an arbitrary function H.

Gauge invariance and charge conservation
So far, nothing has been said about the role of gauge transformations. In this section, we will show how Gauss's law arises as usual from gauge invariance. It is well known that Gauss law arises from the symmetry of the simultaneous action of gauge transformations over both vector potentials and phase-space quantities. In order to avoid unnecessary difficulties, as explained in [36], it is convenient to study the properties of gauge transformations in terms of the canonical coordinates (x, p).
In particular, besides the standard gauge transformation on the electromagnetic quantities (A, E) → (A + ∇χ, E), phase-space coordinates undergo momentum translations of the type (x, p) → (x, p + q∇χ). In turn, this produces an action of gauge transformations on the phasespace density that is f (x, p) = f (x, p − q∇χ). The latter is a type of canonical transformation, which will be the starting point of our discussion. In order to examine the role of gauge transformations, we have to construct an action of momentum translations on the space of KvH wavefunctions. If we were dealing with KvN theory, this would be simply given by Ψ KvN (x, p) → Ψ KvN (x, p − q∇χ). However, KvH wavefunctions also carry a phase factor which we now turn to. If η(x, p) = (x, p + q∇χ) identifies a momentum translation, equation (26) yields ∇ϕ = p − (p + q∇χ) = −q∇χ so that ϕ = −qχ and in this case the phase-space function ϕ depends only on the spatial coordinates x. Then, since χ(η −1 (z)) = χ(x), we are led to the following unitary action Ψ KvH → U η Ψ KvH of momentum translations on KvH wavefunctions: Upon dropping the subscript 'KvH', we will now show that the Hamiltonian (33) is invariant under the gauge transformation (Ψ(x, p), A, E) → e iqχ/ Ψ(x, p − q∇χ), A + ∇χ, E .
Evidently, h Max in (19) is manifestly gauge-invariant and thus here we consider the first term Ψ, L H Ψ . The gauge invariance of this functional is an immediate consequence of the fact that the phase-space density (30) identifies an equivariant momentum map. Indeed, since momentum translations are canonical transformations, we can use (31) to write Then, the transformation A → A + ∇χ leads to overall gauge invariance of Ψ, L H Ψ . At this point, we have characterized the gauge transformations that leave the Hamiltonian h(Ψ, A, E) = Ψ, L H Ψ + h Max (E, A) invariant and we are ready to present the associated conserved quantity. Here, we shall proceed once again by exploiting momentum maps: since the action of gauge transformations on (Ψ, A, E) leaves the Hamiltonian invariant, the momentum map associated to this action is conserved by the dynamics. This momentum map must satisfy the defining relation (18). In this case, the Poisson bracket is given by (34) and we have M = L 2 (T * Q) × T * Ω 1 (Q), where Ω 1 (Q) denotes the space of differential one-forms on Q and Q = R 3 . As discussed in [36], the Lie algebra of gauge transformations is identified with smooth scalar functions ξ(x) on Q, so that their infinitesimal generator reads Then, upon writing the momentum map and computing J (Ψ, A, E), ξ = −ǫ 0 E, ∇ξ − q Ψ, L ξ Ψ , one indeed verifies (18). Thus, Gauss Law (40) emerges as the zero-level set of a conserved momentum map associated to the action of gauge transformations.
We conclude our discussion by noticing that, while the Hamiltonian (33) is gauge invariant, its dependence on A cannot be generally expressed only in terms of the magnetic field B = ∇ × A. While this is precisely what happens also in the standard Hamiltonian treatment of the Maxwell-Vlasov system, here we observe that this feature persists after changing to noncanonical coordinates and this is due to the presence of the Lagrangian function in the prequantum operator (36). As proposed in [36], one can still write the Hamiltonian h(Ψ, A, E) in (39) in terms of (Ψ, B, E) at the expenses of fixing a convenient gauge such as the Coulomb gauge or the Poincaré gauge. For example, the Coulomb gauge yields A = ∇ × (∆ −1 B), which can then be replaced in the expression (36) of the prequantum operator. Then, the Poisson bracket (35) also changes according to the familiar chain rule relation δ/δA = ∇ × δ/δB.

Discussion
The KvN and KvH equations represent two valid approaches to developing a Hilbert space formulation of classical mechanics on phase space. Both approaches lead to a generalized Clebsch representation for the PDF f . A specific choice of the complex phase factor allows this Clebsch representation for f to become equal to the Koopman prescription |Ψ| 2 . However, this choice requires the phase factor to become singular.
In both formulations, the complex phase factor is generally involved in reproducing the classical dynamics. In fact, for the KvN formulation, the additional phase degree of freedom is formally required for obtaining a canonical variational formulation. This canonical KvN-Maxwell formulation parallels the development of the canonical KvH-Maxwell formulation, and can be obtained from the results of Sec. 3 by simply eliminating the Lagrangian from the definition of L H . In fact, this canonical KvN-Maxwell formulation has been treated in [42] in the electrostatic limit. Alternatively, in order to eliminate the need to include the phase in KvN dynamics, a noncanonical Poisson bracket was determined that reduces to the standard Vlasov bracket for functionals that depend only on |Ψ| 2 . In this case, the phase becomes completely irrelevant and can be set to be identically zero.
The "pre-quantum" KvH formulation begins to bridge the gap between the classical and quantum mechanical dynamics by providing a physically motivated prescription for the evolution of the phase factor that agrees with the semiclassical → 0 limit. Thus, the prequantum KvH equation can begin to describe some of the important physical consequences of coupling a classical system to a truly quantum system [5,21,24]. In contrast, the KvN formulation, with trivial phase dynamics, is perhaps better considered to correspond to the diagonal part of the density matrix.
In dealing with both canonical and noncanonical structures, some comments on their numerical aspects are also in place. Indeed, both symplectic numerical integrators and quantum simulation algorithms are well understood for canonical Hamiltonian systems, but not for noncanonical systems with an arbitrary Poisson bracket. While this makes the canonical KvN and KvH formulations more amenable to the development of numerical integration techniques that preserve conservation laws, it also motivates future research on developing numerical methods that target the new noncanonical KvN formulation derived here.
Both approaches can be used to develop a "quantum" representation of the classical Liouville equation and both can be simulated on a quantum computer. However, once the Koopman equation is coupled to Maxwell's equations, one obtains a coupled system of nonlinear partial differential equations. It is only possible to efficiently simulate these equations using a quantum computer if they are embedded within a unitary linear system of equations. This can be done by simulating the classical statistical probability density, F (Ψ, A, E; t, x), for the fields at every point in space-time. As described in Ref. [29], the Liouville equation for F can be simulated efficiently on a quantum computer. Understanding the complexity of quantum simulation for each of these Koopman-Maxwell formulations is an important topic for future research.