Variable-moment fluid closures with Hamiltonian structure

Based on ideas due to Scovel–Weinstein, I present a general framework for constructing fluid moment closures of the Vlasov–Poisson system that exactly preserve that system’s Hamiltonian structure. Notably, the technique applies in any space dimension and produces closures involving arbitrarily-large finite collections of moments. After selecting a desired collection of moments, the Poisson bracket for the closure is uniquely determined. Therefore data-driven fluid closures can be constructed by adjusting the closure Hamiltonian for compatibility with kinetic simulations.


Introduction
The simplest non-dissipative kinetic model for a single-component plasma in N space dimensions is the Vlasov-Poisson system, Here f = f (q, p) denotes the Vlasov distribution function and V = V (q) is the electrostatic potential.Particles that comprise the sole plasma species have charge e, mass m, position q = (q 1 , . . ., q N ), and momentum p = (p 1 , . . ., p N ).A uniform neutralizing background n 0 ensures overall charge neutrality.The explicit time dependence of f and V has been suppressed.This single-component model generalizes simply to allow for multiple plasma species.Kinetic plasma models, like the Vlasov-Poisson system, are formally equivalent to infinite hierarchies of hydrodynamic equations for the fluid moments.The degree-k fluid moment M M M k = M i 1 ...i k dq i 1 ⊗ • • • ⊗ dq i k is a symmetric tensor whose components are given by the momentum space integrals M i 1 ...i k (q) = p i 1 . . .p i k f (q, p) d p, i ℓ = 1, . . ., N, k = 0, 1, 2, . ... In general, the evolution equation for the degree-k fluid moment couples to the degree-(k + 1) fluid moment, a phenomenon that epitomizes the so-called moment closure problem.For example, evolution equations for the first three fluid moments (M 0 , M 1 , M 2 ) = (n, P P P, S) are ∂ t n + m −1 ∇ • P P P = 0 (1) ∂ t S + m −1 ∇ • T = −e(∇V P P P + P P P∇V ), where T = M 3 denotes the degree-3 fluid moment.Any prescription for expressing the degree-(k + 1) fluid moment as a function of lower-degree fluid moments is known as a fluid closure.Upon introducing a fluid closure, the rate of change of the first k fluid moments becomes a function of only those first k fluid moments.Thus, it becomes feasible to solve a finite collection of hydrodynamic equations to find time evolution of the first k fluid moments.Of course, there is no guarantee that moment evolution determined by a fluid closure bears any relation to moment evolution determined by the parent kinetic model, unless the fluid closure is chosen wisely.Traditional methods for developing fluid closures involve some combination of asymptotic expansions [1][2][3][4] and intuitive guesswork.More recently, advances in machine learning 5,6 have accelerated the development of data-driven fluid closures, as well as other data-driven reduced-order models 7 .But fluid closures built using black-box neural networks are known to deliver poor results when employed over large simulation time intervals.There is therefore a growing effort to incorporate inductive bias into data-driven fluid closures, with the hope of improving long-term predictions.For example, Huang et al. 8,9 and Porteous et al. 10 developed variable-moment neural-network-based fluid closures that automatically enforce hyperbolicity of the ensuing moment evolution equations.
This Article presents a framework for building variable-moment fluid closures that automatically enjoy an important structural property inherent to non-dissipative kinetic plasma models: Hamiltonian structure.The framework involves two essential components.(1) Specification of a Poisson bracket for the first m 0 fluid moments, where m 0 is any integer greater than 1. (2) Specification of the plasma system energy as a function of the first m 0 fluid moments.The form of the Poisson bracket only depends on the the number of moments m 0 included in a desired closure.On the other hand, the form of the energy function, or Hamiltonian H , is free to specify.Parameterizing H using a neural network as in Hamiltonian neural networks 11 , or a library of candidate terms 12,13 , offers a novel path toward data-driven closure modeling for non-dissipative kinetic plasmas.As presented here, the framework applies specifically to fluid closures of the Vlasov-Poisson system in N space dimensions.However, similar ideas apply to other non-dissipative kinetic plasma models, such as the Vlasov-Maxwell system [14][15][16] .
The framework still applies in the presence of dissipation caused by particle collisions since a dissipative fluid closure should respect the metriplectic 17 structure of collisional kinetic theory.Notably, metriplectic dynamical systems always decompose as the sum of a Hamiltonian part and a dissipative part; the framework developed here applies to the Hamiltonian part.
The Hamiltonian structure underlying the infinite hierarchy of fluid moments for the Vlasov-Poisson system was first identified by Gibbons 18 .However, finding true Poisson brackets for finite, but arbitrarily large collections of fluid moments for the Vlasov-Poisson system represents an outstanding theoretical challenge.The Poisson bracket for the closure framework presented here provides the first complete solution.In the context of slab drift kinetic models, Tassi [19][20][21] previously found Hamiltonian fluid moment closures with any number of fluid moments.But the methods underlying this result fail to extend to other dissipation-free kinetic models in any obvious manner.In the Vlasov-Poisson context, Holm and Tronci 22 found Hamiltonian fluid moment closures that include arbitrarily-large collections of fluid moments that do not include the density moment.Otherwise, the best available partial solutions for Vlasov-Poisson work in one space dimension (N = 1), where the tensorial nature of fluid moments does not appear, and fall into two categories.The first category contains Hamiltonian moment closures derived by way of direct imposition of the Jacobi identity -Jacobi closures.The second category contains closures based on exact solutions of the parent kinetic model that may be parameterized by fluid moments -exact reductions.
At a schematic level, derivations of Jacobi closures work by first substituting an unknown closure relation into the Poisson bracket for the parent kinetic model, and then finding conditions the closure relation needs to satisfy in order for the closure bracket to satisfy the Jacobi identity.These conditions take the form of nonlinear partial differential equations (PDEs) that, in general, do not admit solutions or cannot be solved in closed form.However, there are notable exceptions.Tassi 23 studied Jacobi closures for drift kinetics that include two moments.Chandre and Shadwick 24 found a Jacobi closure that accommodates all fluid moments with degree k < 5 and ensures stability of uniform solutions.This extends the earlier result 25 27 .Jacobi closures involving moments with degree k < 2 are automatic; the algebra of functionals that depend only on density and momentum density closes under the Vlasov-Maxwell Poisson bracket.The same algebraic closure property is not automatic for models that include moments with degrees k ≥ 2. de Guillebon and Chandre 28 dramatically highlighted this last point by revealing an earlier fluid moment model 29 as a false-positive; a claimed Hamiltonian fluid closure including moments with degree k = 2 that was provably not Hamiltonian.Jacobi closures that include fluid moments of arbitrarily large degree have not yet been discovered.
Derivations of Hamiltonian moment closures based on exact reductions work like the classical variation of parameters method.After postulating a parametric form for the distribution function f , evolution equations for the parameters are derived that ensure a time-evolving parametric f solves the kinetic equation exactly.This procedure only works for very special parametric f .For example, a distribution function of the form f (q, p) = n(q) δ (p − P P P(q)/n(q)), where n is a positive function and P P P is a vector field, solves the Vlasov-Poisson system if and only if (n, P P P) satisfies the Euler equations with a particular barotropic equation of state -the exact reduction.The parameters for this δ -closure are just the fluid moments with degrees k < 2; the Poisson bracket among these parameters induced by the Vlasov-Poisson bracket recovers the k < 2 closure bracket mentioned above.Perin et al. 30 initiated a program aimed at finding moment closures with arbitrarily high-degree moments by re-parameterizing a multiple waterbag exact reduction.(Perin et al. 31 also studied waterbags in the drift kinetic context.)Simple explicit formulas express the fluid moments of multiple waterbags as polynomials in waterbag parameters.Whenever these polynomial relationships can be inverted to give waterbag parameters in terms of fluid moments a Hamiltonian moment closure results.However, due to the general insolubility of large-degree polynomial equations, the inversion has only been demonstrated for a single waterbag.Like Jacobi closures, Hamiltonian closures based on exact reductions have yet to 2/18 accommodate fluid moments of arbitrarily large degree.
The Hamiltonian moment closures presented here fall into the exact reduction category.The relevant parametric form for f is given by Here, n denotes the density, the gradient of the scalar function ψ determines the center of the distribution function, and µ i 1 ...i k denotes the degree-k centered fluid moment.Section 5 contains a detailed exposition.A simple upper-triangular transformation, Eq. ( 20), relates the first m 0 centered fluid moments to the first m 0 standard fluid moments.Hamiltonian closures based on this form of f therefore accommodate arbitrarily high-degree moments by increasing the integer m 0 .This feature of the closures, together with their applicability in every space dimension N, sets them apart from all previous Hamiltonian fluid closures that currently appear in the literature.
The ansatz (4) can reproduce the first m 0 moments of any distribution.The Poisson bracket associated with this exact reduction may therefore be used to model the first m 0 moments of any distribution.In particular, it is not necessary to assume that the distributions modeled by the new closures are exactly of the form (4).This observation forms the basis for the proposed closure modeling approach, described with greater detail in Section 6: fix the Poisson bracket and incorporate all modeling choices into the form of the Hamiltonian H .However, it is important to keep in mind that the theory developed here only guarantees the existence of closures that apply to f close to the form (4). Applications of the new closures to, e.g., bump-on-tail distributions should be considered exploratory efforts.
Demonstrating that (4) comprises an exact reduction and uncovering the corresponding Poisson bracket structure for general m 0 and N requires considerable bookkeeping technology.Fortunately, Scovel and Weinstein 32 developed the necessary general machinery in the context of distributions f localized near a point z = (q, p) in single-particle phase space.(Channell 33 provides additional insights into elements of the Scovel-Weinstein construction.)Such distributions, which are relevant to self-consistent beam dynamics in particle accelerators, also give rise to a moment closure problem, where the moments are redefined as the phase space integrals In order to find a finite-dimensional Poisson bracket that describes dynamics of the first m 0 phase space moments, Scovel-Weinstein developed a reduction strategy that can be applied to any Lie-Poisson Hamiltonian system whose associated Lie algebra g decomposes as a sum of subalgebras g = b + s.(The Lie bracket on g need not factor as the sum of brackets from b and s.)For the beam problem, g comprises the phase space functions with multiplication provided by the usual single-particle Poisson bracket, b contains the affine phase space functions, and s is the collection of phase space functions that vanish at (q, p) = (0, 0), together with their first partial derivatives.For fluid moments, the Lie algebra g remains the same, but the subalgebras b, s differ substantially.Section 2 summarizes the Lie theory necessary to apply Scovel-Weinstein theory to fluid moments, and in particular defines the subalgebras b and s.Section 3 precisely states the elements of Scovel-Weinstein theory that enable subsequent development of the new closure method in Sections 4, 5, and 6.
The following notation is used in what follows without further comment.If V = R d is any finite-dimensional vector space then C(V ) contains smooth functions on V with mild growth (e.g. at most polynomial grown) at infinity.The dual space C * (V ) contains rapidly decaying distributional functions on V .The duality pairing between κ ∈ C(V ) and ρ ∈ C * (V ) is ρ, κ = V κ(v) ρ(v) dv, where dv denotes the Lebesgue measure on R d .

Lie theory for the contactomorphism group
A complete understanding of the Hamiltonian moment closures described in this Article requires familiarity with some of the rich structure of the contactomorphism group G.This Section aims to provide a self-contained introduction to G and its properties with greatest relevance to fluid moments.For this Section and Section 3 only, readers should be familiar with the basic mathematical theory of Lie groups at the level of Chapter 4 in the book 34 by Abraham and Marsden.Marsden and Weinstein 15 were first to recognize that G is the Lie group that integrates the Lie algebra underlying the Vlasov equation.Subsequently, Gay-Balmaz and Tronci 35 used structural properties of contactomorphisms to show that the integrable Bloch-Iserles system is geodesic flow on a subgroup of G. Gay-Balmaz and Tronci 36, 37 also provide additional exposition concerning the contactomorphism group as part of their ongoing efforts in quantum-classical hybrid modeling.
Let Q ∋ q be an N-dimensional vector space representing the configuration space for an N-degree of freedom Hamiltonian system.The corresponding momentum phase space is P = Q × Q * ∋ (q, p), where Q * denotes the dual vector space to Q.The cotangent projection map is π : P → Q : (q, p) → q.Each p ∈ Q * is a linear mapping p : Q → R. Application of the linear map p ∈ Q * to a vector q ∈ Q will be denoted p, q ∈ R. The canonical 1-form ϑ on P is defined by ϑ = p, dq .The corresponding canonical symplectic form is the 2-form ω = −dϑ given by exterior differentiation of ϑ .If h : P → R is any smooth function on P, its partial derivatives δ h/δ q and δ h/δ p are defined by the directional derivative formula Note that the value of δ h/δ q at a particular (q, p) ∈ P lies in Q * , while δ h/δ p lies in Q.The Poisson bracket associated with the canonical symplectic form ω is given on pairs of functions f , g : P → R by A (strict) contactomorphism g = (Φ, ϕ) is a pair comprising a diffeomorphism Φ : P → P and a smooth function ϕ : P → R that together satisfy the identity The product of a pair of contactomorphisms g 1 , g 2 , denoted g 1 g 2 , is When equipped with this product, the set of all contactomorphisms comprise a group known as the (strict) contactomorphism group.
The contactomorphisms b = (Φ, ϕ) such that form a subgroup B, known as the group of fiber translations.Note that if τ ψ : P → P denotes the mapping (q, p) → (q, p − δ ψ/δ q), then every fiber translation may be written as b = (τ ψ , π * ψ) for some unique ψ.The mapping I B : B → C(Q) : (τ ψ , π * ψ) → ψ, from B to the Abelian group of smooth functions on configuration space is a group isomoprhism.Let Z = {(q, p) ∈ P | p = 0} denote the zero section in phase space.The submanifold Z contains the static system states.The contactomorphisms s = (Φ, ϕ) such that Φ(Z) = Z and ϕ(Z) = 0 form a Lie subgroup S ⊂ G.In this work, S will be referred to as the group of isostatic contactomorphisms.The isostatic contactomorphisms complement the fiber translations in the following sense.
Theorem 1.The product mapping m : B × S → G : (b, s) → bs, restricts to a diffemorphism from a neighborhood of (e, e) ∈ B × S to a neighborhood of e ∈ G.In particular, for each g ∈ G sufficiently close to the identity there is a unique fiber translation b and a unique isostatic contactomorphism s such that g = bs.
Proof.The proof proceeds by constructing the inverse of m 0 on an appropriate neighborhood U G of the identity in G.In particular, given g ∈ G the key step is finding a cooresponding Observe that the zero section Z is a graph over In other words, ψ is (minus) a generating function for the Lagrangian submanifold Φ(Z).A generating function ψ for Φ(Z) determines a fiber translation b = (τ ψ , π * ψ), which very nearly furnishes the b = b(g) mentioned above.However, there is more to do because (τ ψ , π * ψ) −1 g is not necessarily in S.Moreover, the fiber translation (τ ψ , π * ψ) is not uniquely determined by g.In fact, ψ ′ = ψ + c, c ∈ R, also generates Φ(Z), but determines a different fiber translation The following argument shows that for each g ∈ U G there is a unique generating function ψ = ψ(g) for Φ(Z) such that s = (τ ψ , π * ψ) −1 g ∈ S. Let ι Q : Q → P : q → (q, 0) denote the map that embeds Q in P as Z, and set As mentioned already, there exists a (non-unique) ψ that generates Φ(Z).Since Φ(Z) is the graph of −δ ψ/δ q, the mapping τ −1 ψ maps the Lagrangian submanifold This means ϕ s (q, 0) = ϕ(q, 0) − ψ(η s (q)) is constant as a function of q.Eliminate the ambiguity in the definition of ψ by requiring that ϕ s (q, 0) = 0. Once this is done, b = b(g) = (τ ψ(g) , π * ψ(g)) is uniquely determined by g and b −1 g ∈ S, as desired.
The inverse of The Lie algebra g = T e G of infinitesmial contactomorphisms comprises pairs (X, Σ), where X is a vector field on P, Σ : P → R is a function on P, and dΣ + L X ϑ = 0.The adjoint action of G and Lie bracket on g are given by The linear map I g : g → C(P) : (X, Σ) → Σ + ϑ (X), from infinitesimal contactomorphisms to phase space functions under Poisson bracket is a Lie algebra isomorphism.This observation recovers the well-known fact that the contactomorphism group integrates the Lie algebra of phase space functions equipped with the Poisson bracket.
The Lie subalgebra b = T e B ⊂ g of infinitesimal fiber translations comprises pairs (X, Σ) ∈ g, where Σ = π * σ for some σ ∈ C(Q) and X = X π * σ .The Lie bracket on b vanishes because B is Abelian.When identifying g with phase space functions via I g , b corresponds to p-independent functions.Denote the image Lie subalgebra of b in (C(P), {•, •}) as The Lie subalgebra s = T e S ⊂ g of isostatic infinitesimal contactomorphisms comprises pairs (X, Σ) ∈ g such that X is tangent to Z and Σ(Z) = 0. (In fact, the latter condition implies the former, and so can be taken as the sole property characterizing s as a subspace of g.)The image Lie subalgebra Theorem 1 implies that, as a vector space, the Lie algebra of infinitesimal contactomorphisms is the direct sum of b and s, g = b + s.The corresponding statement in C(P), namely C(P) = C b (P) + C s (P), may be understood more simply as a byproduct of Taylor's theorem with remainder.
The dual space g * comprises linear functions µ : g → R. Since g = b + s, the duals to b and s are subspaces of g * given by Similarly, the dual space to C(P), C * (P), contains linear functions f : C(P) → R, and the duals C * b (P), C * s (P), are the subspaces of C * (P) given by The coadjoint action of G on g * is given by Ad * g : Provided the space of phase space functions is restricted to smooth functions with slow growth at infinity, the dual C * (P) contains all of the typical smooth Vlasov densities.It includes more exotic functionals as well, such as Dirac δ -functions and their derivatives.

5/18 3 Scovel-Weinstein theory
The (infinite-dimensional) Hamiltonian structure underlying the Vlasov-Poisson system is dual to the Lie algebra C(P) of phase space functions under the standard (finite-dimensional) Poisson bracket.The problem of finding reduced models, in particular reduced fluid models, for Vlasov-Poisson may therefore be approached at the level of the algebra C(P).Adopting this point of view requires addressing how the Lie algebra C(P) might be "reduced."For the purposes of this Article, any reduction should somehow encode the objects dual to the first m 0 fluid moments, namely the phase space functions h(q, p) that are polynomial in momentum p with degree at most m 0 .Plausible reductions of C(P) include proper Lie subalgebras: subspaces of C(P) that close under Poisson bracket.For example, the phase space functions that only depend on the configuration variable q (i.e.degree-0 momentum polynomials) comprise an Abelian Lie subalgebra dual to the space of density moments.More generally, the momentum polynomials with degree at most 1 close under Poisson bracket because However, by itself, subalgebra reduction of C(P) provides insufficient flexibility to accommodate fluid moments with degree between 0 and m 0 when m 0 > 1.This follows from the failure of momentum polynomials with degree at most 2 to close under Poisson bracket.Quotient algebra reduction offers a more nuanced alternative to subalgebra reduction.Instead of attempting to encapsulate moments in a small subalgebra, quotient algebra reduction aims to find large ideals I in C(P) such that the quotient space C I (P) = C(P)/I encodes fluid moments.Here, ideals are subalgebras with an additional absorbing property: the bracket of functions in I with arbitrary functions in C(P) must be contained in I .The subspace I being an ideal ensures that the quotient C I (P) inherits a Lie algebra structure from C(P).Unfortunately, as for subalgebra reduction, quotient algebra reduction cannot accommodate the fluid moments with degree between 0 and m 0 by itself.As described in detail within the proof of Theorem 2, it can successfully encode moments with degree between 1 and m 0 .
The paper 32 by Scovel and Weinstein devises an alternative Lie-theoretic reduction strategy that hybridizes subalgebra reduction and quotient algebra reduction.Developed originally for application to phase space moments, Scovel-Weinstein theory adapts to the fluid moment context in straightforward fashion.In fact, their methods apply any time a Lie group G admits a decomposition by Lie subgroups B, S ⊂ G. Here, decomposition means that for each g ∈ G there are unique b ∈ B and s ∈ S such that g = bs.In the context of fluid moments, the group G is the contactomorphism group, B is the group of fiber translations, and S is the group of isostatic contactomorphisms.Each of these groups is described with some detail in Section 2. The technical heart of Scovel-Weinstein theory may be described as follows.
Suppose that G is a Lie group that admits a decomposition by Lie subgroups B, S ⊂ G, G = BS.The Lie algebra g then admits a corresponding direct sum decomposition g = b + s by Lie subalgebras.Likewise, the dual space to g admits the direct sum decomposition g * = b * + s * .The Scovel-Weinstein map Γ : Here Ad * g = (Ad Note that the space (B × b * ) × s * is larger that g * by the factor B. This extra variable b is the price paid for "untangling" the Lie-Poisson structure on g * in terms of standard Poisson structures associated with b * and s * .Scovel-Weinstein advocate the following interpretation of Γ and its source and target spaces.In the space (B × b * ) × s * the factor g * 0 = b * × s * represents a "reference" copy of g * sitting at the identity e in the group B. The mapping Γ advects the reference space g * 0 from e to b.This advection interpretation warrants further elaboration.Along a solution of the Vlasov-Poisson system, the distribution function f (q, p) is advected by canonical transformations that prescribe the motion of individual particles in phase space.These canonical transformations naturally reside with the contactomorphism group G.The Scovel-Weinstein map Γ may therefore be interpreted as a general way of identifying special system motions advected by a subgroup B ⊂ G. Now suppose that the Lie subalgebra s contains an ideal I ⊂ s.Note that ideals in s are not necessarily ideals in g!The quotient algebra s I = s/I then inherits a Lie algebra structure from s by requiring that the quotient map Π I : s → s I is a

Poisson bracket on the space of truncated fluid moment vectors
This Section uses Poisson-geometric and Lie-theoretic methods to identify a Poisson bracket on a space that encodes density, fluid momentum, pressure, and higher fluid moments up to some truncation order m 0 .The precise relationship between the Poisson bracket constructed here and the Lie Poisson bracket on the space of distribution functions f ∈ C * (P) will be described in Section 5.
A degree-k fluid moment M k , k ∈ {0, 1, 2, . ..}, is a symmetric, covariant, degree-k tensor field on configuration space Q. (Geometrically, fluid moments are naturally tensor densities 38 instead of ordinary tensors, but here these two kinds of objects are identified using the standard flat metric on configuration space.)If {q i } i=1,...,N denotes the standard coordinate system on Q = R N , M k may be written as Any phase space distribution f = f (q, p) that decays sufficiently fast as |p| → ∞ determines fluid moments of each degree by the formula Conversely, under certain conditions a sequence of fluid moments M k , k = 0, . . ., ∞, determines a distribution f .This is true, for example, when the moment generating function associated with the sequence M k exists.A fluid moment vector with degree between n ≥ 0 and m ≥ n is a tuple M M M = (M n , M n+1 , . . ., M m ) whose i th entry is a fluid moment of degree i. Denote the set of all such fluid moment vectors M m 0 n .Note that M k k is the vector space of degree-k fluid moments, and that Also note that elements of M 0 0 are just fluid densities n ∈ C * (Q).Although specifying moments of each degree k = 0, . . ., ∞ formally specifies a unique phase space distribution, specifying only the moments with degree k ≤ m 0 does not.Additional data must be supplied to recover a distribution function from a fluid moment vector in M m 0 0 .In this work that additional data will comprise a scalar configuration space function ψ ∈ C(Q) whose derivative dψ(q) specifies the "center" of the distribution f (q, p) at q.The precise goal of this Section is constructing a Poisson bracket on the space for any non-negative integer m 0 .Here the symbol µ µ µ is used in place of M M M to emphasize that the moment vector µ µ µ contains centered moments.
The Poisson bracket on P m 0 originates from a Lie algebra structure on the objects dual to fluid moment vectors: fluid comoment vectors.A fluid comoment C k of degree k ∈ {0, 1, 2, . ..} is a symmetric, contravariant, degree-k tensor field on 7/18 configuration space Q.If {e i } i∈{1,...,N} denotes the standard basis for Q, evey degree-k fluid comoment may be expressed using index notation as where the array of component functions C i 1 i 2 ... i k is symmetric under index permutation.A fluid comoment vector with degree between n ≥ 0 and m ≥ n is a tuple C C C = (C n ,C n+1 , . . .,C m ) whose i th entry is a fluid comoment of degree i.The vector space of all fluid comoment vectors with degree between n and m 0 is C m 0 n .Fluid comoments with degrees k = 0, 1, 2 will be denoted using the special symbols σ , w w w, Q, respectively.The duality pairing between fluid moment vectors M m 0 n and fluid comoment vectors C m 0 n is given by By way of this pairing, the fluid moment vectors arise as the dual space to fluid comoment vectors, M m 0 n = (C m 0 n ) * .Fluid comoment vectors with degree between 1 and m 0 ≥ 1 comprise a Lie algebra.This result is striking because the comoment vectors with degree between 0 and m 0 do not comprise a Lie algebra!The following Theorem exhibits the relevant Lie bracket and proves that it satisfies the Jacobi identity.
Theorem 2. Let m 0 ≥ 1 be an integer.The space C m 0 1 of fluid comoment vectors with degree between 1 and m 0 is a Lie algebra.The Lie bracket between pairs of comoment vectors C C C 1 ,C C C 2 is given by where the sum is over all integer pairs (k Remark 1.This formula for the Lie bracket on C m 0 1 reproduces the indicial form of the so-called Schouten concomitant 38,39 between symmetric contravariant tensors, modulo degree m 0 + 1.This observation is consistent with the method of proof given below since the Schouten concomitant is the unique Lie algebra structure on symmetric tensors that makes the map sending symmetric tensors to momentum polynomials a Lie algebra homomorphism.However, note that the symmetric contravariant tensors with degree between 1 and m 0 are not closed under the Schouten concomitant.Instead, these tensors possess a quotient Lie algebra structure.A similar quotient algebra structure was first noticed by Channell 33 in the context of phase space moments.Partially motivated by Channell's observations, Holm and Tronci 22 found the Lie-Poisson bracket associated with (5) without identifying its corresponding Lie-algebraic origin.Therefore, Theorem 2 may be formulated in the following alternative manner: the Poisson bracket for fluid moments with degree between 1 and m 0 presented in Holm and Tronci's Section 2.3 is a Lie Poisson bracket for the Lie algebra whose bracket is given by (2).
Proof.First observe that for any non-negative integers n 1 ≤ n 2 there is a linear isomorphism I n 2 n 1 : n 1 and the space C n 2 n 1 (P) of momentum polynomials on phase space P with degree between n 1 and n Note the the sum over comoment degree is denoted explicitly, while the sums over tensor indices are implicit.The proof now proceeds by examining momentum polynomials in place of comoment vectors.
Let C s (P) denote the smooth phase space functions that vanish when p = 0. Note that C ℓ 1 (P) ⊂ C s (P) for each integer ℓ ≥ 1.The space C s (P) a Lie subalgebra of the Lie algebra of phase space functions under Poisson bracket.This follows immediately from the discussion of Lie subgroups of contactomorphisms in Section 2. Alternatively it can be verified directly as follows.Suppose s 1 , s 2 ∈ C s (P).The Poisson bracket of these functions is The sums on the right-hand-side each vanish when p = 0 because (∂ q i s ℓ )(q, 0) = ∂ q i (s ℓ (q, 0)) = ∂ q i (0) = 0, ℓ = 1, 2. Therefore {s 1 , s 2 } ∈ C s (P), as claimed.

8/18
For each integer ℓ ≥ 0 let I ℓ ⊂ C s (P) denote the phase space functions that vanish together with their first ℓ momentum derivatives when p = 0.If r ∈ I ℓ and s ∈ C s (P) the Poisson bracket {r, s} is {r, s} = ∂ q i r ∂ p i s − ∂ q i s ∂ p i r.
Differentiating the first sum with respect to momentum ℓ times or less produces a function that vanishes when p = 0 because ∂ q i r vanishes together with its first ℓ momentum derivatives when p = 0. Differentiating the second sum with respect to momentum fewer than ℓ times likewise produces a function that vanishes when p = 0 because ∂ p i r vanishes when p = 0 together with its first ℓ − 1 momentum derivatives.The only contribution to the ℓ th -order momentum derivative of the second sum that does not vanish at p = 0 due to vanishing of ∂ p i r is a product of ℓ momentum derivatives of ∂ p i r with ∂ q i s; ℓ th -order derivatives of ∂ p i r need not vanish when p = 0.But this product vanishes when p = 0 because ∂ q i s vanishes when p = 0. Therefore {r, s} ∈ I ℓ .In other words I ℓ is an ideal inside the Lie algebra C s (P) for each integer ℓ ≥ 0.
The quotient space C s m 0 (P) = C s (P)/I m 0 inherits the structure of a Lie algebra because I m 0 is an ideal.The Lie bracket is given by Poisson bracket mod There is a canonical isomorphism that sends each equivalence class in C s m (P) to its unique representative that is a momentum space polynomial with degree between 1 and m 0 .Equivalently, C s m 0 (P) ≈ C m 0 1 (P).Therefore C m 0 1 (P) is a Lie algebra with Lie bracket given by Here mod p m+1 denotes setting all polynomial coefficients with degree greater than m 0 equal to 0. In particular, if Here round braces around a grouping of indices denotes averaging over all permutations of that grouping.Finally, since C 1 that makes I m 0 1 an isomorphism of Lie algebras.In symbols, the bracket is Directly computing the right-hand-side of this formula with the aid of Eq. ( 6) reproduces Eq. ( 5), as claimed.
The Lie algebra structure on the comoment vectors C m 0 1 with degree between 1 and m 0 implies that the dual space (C m 0 1 ) * is a Poisson manifold with a Lie-Poisson bracket.But this dual space coincides with the space of moment vectors with degree between 1 and m 0 , (C m 0 1 ) * = M m 0 1 .The corresponding Lie-Poisson bracket {•, •} m 0 1 on M m 0 1 comes very close to providing the sought-after Poisson structure on P m 0 .However, the degree-0 moments, those that correspond to density, are conspicuously missing from M m 0 1 .The scalar field ψ ∈ C(Q) that specifies the distribution center is also not contained in M m 0 1 .This shortcoming of the bracket on M m 0 1 , as well as the analogous bracket on phase space moments with degree between 2 and m 0 , has been recognized previously 22,33 .
Further progress comes upon noticing that the space of fluid densities n ∈ C * (Q) is dual to the space of scalars ψ ∈ C(Q).This means there is a natural Poisson bracket {•, •} 0 on the space C(Q) × C * (Q) given on pairs of functionals F, G by It follows that there is a natural Poisson bracket on each of the factors in 1 .The simplest Poisson bracket on P m 0 that incorporates the Lie Poisson bracket on M

9/18
The Jacobi identity for this bracket is automatic since the direct sum of Poisson manifolds is again a Poisson manifold.Equation ( 7) provides the Poisson bracket for the Hamiltonian fluid moment closures presented in this Article.Its relevance to moment dynamics will be established in Section 5, where it will become clear that this bracket naturally addresses the shortcoming of the Lie-Poisson bracket on M m 0 1 in isolation.Note that, for each m 0 , this bracket accommodates all fluid moments with degree less than m 0 , including the density moment.Also note that there is no restriction on the dimension of configuration space N.
When m 0 = 0, the bracket {•, •} m 0 is given by which provides the appropriate Poisson bracket for Hamiltonian fluid closures involving ψ and n only.This is also the bracket for the Madelung transform of the Schrodinger equation.
When m 0 = 1, with µ 1 = P P P 0 , the bracket {•, •} m 0 is given by which provides the appropriate Poisson bracket for Hamiltonian fluid closures involving ψ and fluid moments with degree k < 2.
When m 0 = 2, with µ 2 = S 0 , the bracket {•, •} m 0 is given by where A s = A + A T for a dyad A.

Hamiltonian dynamics of localized momentum distributions
Extending the observations of Channell 33 , Scovel-Weinstein 32 provide a technical device that can be used to relate the product Poisson bracket on 1 described in Section 4 to the Lie-Poisson bracket on C * (P) that underlies the Vlasov-Poisson system.The following Theorem represents an application of that device, which was originally developed for treating phase space moments, to fluid moments.
is a Poisson map.Proof.Let G be a Lie group with Lie subgroups B, S ⊂ G such that the multiplication map B × S → G : (b, s) → bs is a diffeomorphism.As summarized in Section 3, Scovel-Weinstein proved that defines a Poisson map from the product Poisson manifold (B × b * ) × (s * ) to the Lie-Poisson space g * .Here Ad * g denotes the (left) coadjoint action of G on g * .Let G be the contactomorphism group, B ⊂ G the fiber translations, and S ⊂ G the isostatic contactomorphisms, as defined in Section 2. Identify the fiber translations B with the Abelian group of scalar functions ψ ∈ C(Q).Identify the Lie algebras b, s, g, with configuration space scalars σ ∈ C(Q) ≈ b, phase space functions s ∈ C s (P) ≈ s that vanish when p = 0, and general phase space functions h ∈ C(P) ≈ g.The dual spaces are then (distributional) configuration space functions n ∈ C * (Q) ≈ b * , (distributional) phase space functions f0 ∈ C * s (P) with zero fluid density, and general (distributional) phase space functions f ∈ C * (P) ≈ g * .The corresponding duality pairings are the usual L 2 integration pairings.Under these identifications, the Poisson map Γ : Note that C(Q)×C * (Q)×C * s (P) is endowed with the product Poisson structure corresponding to the canonical poisson bracket on C(Q) × C * (Q) and the Lie Poisson bracket on C * s (P).As shown in the proof of Theorem 2, the phase space functions r ∈ I m 0 that vanish along with their first m 0 momentum derivatives when p = 0 form an ideal in the Lie algebra C s (P).Therefore the quotient space C s (P)/I m 0 = C s m 0 (P) has a natural Lie algebra structure and the quotient map Π m 0 : C s (P) → C s m 0 (P) is a Lie algebra homomorphism.Also from the proof of Theorem 2, the quotient space C s m 0 (P) is naturally isomorphic as a Lie algebra to the space of comoment vectors C m 0 1 with degree between 1 and m 0 .(The Lie algebra structure on C m 0 1 is the one provided by Theorem 2.) Under this identification, Π m 0 : C s (P) → C m 0 1 is the map that assigns to each s ∈ C s (P) its weighted momentum Taylor series coefficients with degrees between 1 and m 0 : Since Π m 0 is a Lie algebra homomorphism, the adjoint Π * m 0 : M m 0 1 → C * s (P) is a Poisson mapping between Lie-Poisson spaces.A simple direct calculation leads to the explicit formula Since Π * m 0 and Γ are Poisson maps, the composition is a Poisson map as well.Evaluating the composition using the above formulas shows that Γ m 0 is the map given in the Theorem statement.
Theorem 3 is quite powerful due to a construction 41 known as "collectivisation."Given a Poisson mapping Γ : Λ → Z between Poisson manifolds Λ, Z, every Hamiltonian function h on Z gives rise to a corresponding Hamiltonian function H = Γ * h on Λ known as the "collective" Hamiltonian.Remarkably, if λ (t) is any integral curve for the collective Hamiltonian vector field X H then the image curve z(t) = Γ(λ (t)) is an integral curve for the Hamiltonian vector field X h .This follows from the simple chain of identities for arbitrary q ∈ C(Z): In the context of fluid moments, the collectivisation construction and Theorem 3 have the following implication.Let 1 .If (ψ(t), n(t), µ µ µ(t)) is any solution of the Hamiltonian system on P m 0 defined by H m 0 then f (t) = Γ m 0 (ψ(t), n(t), µ µ µ(t)) is an exact solution of the Vlasov-Poisson system.In other words, Theorem 3 produces exact solutions of Vlasov-Poisson with any finite number m 0 of fluid moments!

11/18
When m 0 = 0, the space P m 0 is just P 0 = C(Q) × C * (Q) ∋ (ψ, n), with the standard canonical Poisson bracket (8).Pulling back the physical Hamiltonian H ( f ) along Γ 0 leads to the collective Hamiltonian The functional derivatives of the collective Hamiltonian are given by The corresponding equations of motion are therefore Upon introducing the velocity variable v v v = −m −1 ∇ψ, these equations may be re-written as which reveals equivalence with Euler's equations for a cold, irrotational fluid.It is a well-known fact 42 that such irrotational flows correspond to exact solutions of the Vlasov-Poisson system.When m 0 = 1, with µ 1 = P P P 0 , the space P m 0 becomes 1 equipped with the Poisson bracket in Eq. ( 9).Pulling back the physical Hamiltonian along Γ 1 leads to the collective Hamiltonian H 1 (ψ, n, P P P 0 ) = 1 2mn The functional derivatives of the collective Hamiltonian are The corresponding equations of motion are therefore In terms of these variables, the equations of motion are instead ∂ t (mn) + ∇ • P P P = 0, ∂ t P P P + ∇ These equations reveal that the exact reduction with m 0 = 1 again recovers familiar hydrodynamic equations, but now with a non-trivial closure relation for the momentum flux tensor S.This Hamiltonian moment closure appears to be new, but it would be interesting to compare it with previously-derived Hamiltonian moment closures that include moments with degree k < 2.
When m 0 = 2, with µ 2 = S 0 , the space 12/18 The equations of motion are therefore given by Note that these equations imply Upon introducing the conventional (non-centered) moments, and performing a tedious explicit calculation, the evolution equations for (ψ, n, P P P 0 , S 0 ) may be re-written as ∂ t n = −m −1 ∇ • P P P (14) ∂ t S = −m −1 ∇ • T − e (∇ V P P P + P P P∇ V ).
These evolution equations for (n, P P P, S) exactly recover the first three equations in the moment hierarchy (3) with a novel closure for the third-degree fluid moment T = p p p f d p given by Here ⊙ denotes the symmetrized tensor product, e.g.
The system (13)-( 16) admits stationary homogeneous equilibria of the form ψ = 0, n = n 0 = const.,P P P = 0, S = S 0 = const.When N = 1 and S 0 ≥ 0, these equilibria are stable if and only if S 0 = 0.When S 0 > 0 the linearized spectrum for Fourier mode k contains a pair of imaginary eigenvalues corresponding to oscillations at the Bohm-Gross frequency, and a pair of real eigenvalues corresponding to exponential growth and decay with rates (Curiously, equilibria with S 0 < 0 are linearly stable!)So the exact closure with m 0 = 2 already begins to reflect the wellknown fact that only certain homogeneous equilibria for the Vlasov-Poisson system, e.g.monotone-decreasing distributions 43 , are linearly stable.Note that the condition S 0 > 0 does not imply a stable distribution; there are many linearly-unstable homogeneous Vlasov-Poisson equilibria with positive second moment.Also note that the growth rate γ = O( |k|) for large k, implying ill-posedness of the initial value problem.Physical solutions must therefore be found within the center manifolds 2,4,44 attached to such equilibria, much as physical solutions of the Abraham-Lorentz-Dirac equation for individual electrons are contained in a center manifold 45 .Further work is needed to determine linear stability properties of homogeneous equilibria with N > 1, m 0 > 2. It is unclear whether the data-driven Hamiltonians described in Section 6 will also require reduction to a center manifold to describe physical dynamics.

A closure scheme for general distributions
The theory developed so far leads to a framework for constructing data-driven fluid closures of the Vlasov-Poisson system with Hamiltonian structure.This framework assumes availability of a legacy simulation code for the Vlasov-Poisson system; either particle-in-cell (PIC) codes or continuum codes will suffice.The framework presented here should be contrasted with the structure-preserving reduced-order modeling techniques developed by Hesthaven, Pagliantini, and Ripamonti [46][47][48] .Those reduced-order models assume the reduced, data-driven dependent variables admit constant symplectic or Poisson structures, or else rely on transformations to such variables.In contrast, the method presented here assumes the reduced dependent variables admit the non-constant Poisson structure provided by Scovel-Weinstein theory.In addition, the reduced dependent variables used here are always simply related to physical quantities of interest, namely the fluid moments.Developing a closure that includes fluid moments up to degree-m 0 requires a parametric energy function and a model that assigns a centered distribution f 0 = F 0 (µ µ µ) to each centered moment vector µ µ µ ∈ M m 0 0 .This terminology has been borrowed from Scovel-Weinstein 32 , who introduced the notion of a model for phase space moments.For consistency, the model distribution function F 0 (µ µ µ) must have moments compatible with µ µ µ: Given a field ψ ∈ C(Q), the ordinary distribtion function f associated with F 0 (µ µ µ) is f = τ ψ * F 0 (µ µ µ).To fully specify the closure, the variable W that parameterizes the energy must be adjusted to an optimal value W * .The closure itself will then comprise a Hamiltonian system on C(Q) × M m 0 0 with Hamiltonian H W * and Poisson bracket given by Eq. ( 7).Finding the optimal parameter value W * proceeds as follows.Start by selecting a finite collection D of samples (ψ, µ µ µ) ∈ D from the space C(Q) × M m 0 0 , together with a timestep parameter ∆t ∈ R. For each sample (ψ, µ µ µ) ∈ D perform two computations.(1) Approximately evolve (ψ, µ µ µ) along the Hamiltonian vector field X H W = ( ψH W , μ µ µ H W ) for ∆t seconds using the forward-Euler formula Then compute the corresponding candidate evolved moment vector M M M(∆t) = (M 0 (∆t), . . ., M m 0 (∆t)), where the degree-k candidate moment is given by Note that M k (∆t) depends explicitly on W . (2) Reconstruct the distribution function f = τ ψ * F 0 (µ µ µ).Use the legacy simulation code to find the time evolution f (∆t) of f .Then compute the corresponding "ground-truth" evolved moment vector M M M(∆t), where the degree-k ground-truth evolved moment is given by With each sample (ψ, µ µ µ) ∈ D now processed, construct the mean-squared loss and find W * by solving the optimization problem In this way, the energy parameters W are chosen in order to maximize agreement between moment evolution predicted by the Vlasov-Poisson system and moment evolution predicted by the Hamiltonian hydrodynamic model for the field ψ and the centered moments µ µ µ.
Through application of the binomial formula for the degree-k candidate evolved moment the W -dependence of the loss S(W ) can be made more explicit.For example, when m 0 = 1, with µ 0 = n, µ 1 ≡ P P P 0 , the

Discussion
This Article has described a new method for developing fluid moment closures of the Vlasov-Poisson system with Hamiltonian structure.These closures apply in any number of space dimensions, and allow for collections of fluid moments with degree between 0 and m 0 , where m 0 is any non-negative integer.The closures work by introducing an irrotational velocity field v v v that specifies the "center" of the distribution function in moment space.The theory presented here leads to both novel exact solutions of the Vlasov-Poisson system parameterized by arbitrarily-large collections of fluid moments, as well as a general framework for finding data-driven fluid closures with Hamiltonian structure.Variable-moment data-driven fluid closures garner increasing attention in the machine learning community.Recent results [8][9][10] indicate that it is possible to construct such closures while strictly enforcing hyperbolicity.The results presented in the present Article therefore suggest investigation of the following intriguing question.Is it possible to develop expressive variable-moment data-driven fluid closures that strictly enforce both hyperbolicity and Hamiltonian structure?Within the closure framework presented here, such a constraint would limit the allowable forms of the closure Hamiltonian H W (ψ, n, µ µ µ).
The fluid closures presented here always include all fluid moments up to and including some maximum degree m 0 .It would be interesting to determine whether these closures can be specialized to only include a subset of these moments.For example, are there subclosures with scalar pressure?
This Article has treated the electrostatic limit of collisionless plasma kinetic theory.Variable-moment fluid closures with Hamiltonian structure would also be interesting to develop for the electromagnetic problem, as embodied by the Vlasov-Maxwell system.Forthcoming publications will reveal that the Lie-theoretic methods used here transfer to the electromagnetic problem in a mathematically and physically appealing fashion.
Each of the fluid closures developed here is associated with a certain direct-sum Poisson bracket {•, •} m 0 .(See Eq. (7).)It would be useful to identify Casimir functions for this bracket and use them to study stability of solutions using the energy-Casimir method 49 .Finding the Casimirs would also be a necessary step when including dissipation in the fluid closures presented here using a metriplectic 17 bracket.
As mentioned in the introduction, Tassi [19][20][21] has found variable-moment fluid closures for slab drift-kinetic systems.But the methods used by Tassi to find closures are seemingly unrelated to the Lie-theoretic methods of Scovel-Weinstein.Is it possible to find Tassi's drift kinetic closures using Scovel-Weinstein theory 32 as was done for the closures presented here?An affirmative result might suggest a path toward variable-moment gyrofluid closures for more intricate gyroaveraged kinetic theories such as fullf gyrokinetics [50][51][52] .In a similar vein, it would be interesting to determine whether the closures presented here can be understood as limits of the multiple waterbag reductions described by Perin et al. 30 .
More broadly, this Article suggests that the methods developed by Scovel-Weinstein 32 in the context of beam dynamics apply to a wider class of dissipation-free statistical closure problems than originally envisioned.There are many possible extensions.For example, consider the problem of solving the functional Vlasov equation for a variational weakly-nonlinear wave equation with random initial data.The algebra of functionals on the field-theoretic phase space enjoys a canonical Poisson bracket.By decomposing this algebra into affine functionals and functionals that vanish together with their first functional derivatives for the vacuum state, Scovel-Weinstein theory may be applied.The result will be Hamiltonian evolution equations for a finite collection of field correlation functions.These non-dissipative correlation dynamics do not seem to have been studied previously.

− 1 )
* denotes the (left) coadjoint action of G on g * .Let {•, •} T * B denote the Poisson bracket on B × b * given by pulling back the canonical Poisson bracket on T * B along the left trivialization map.Let {•, •} s * denote the Lie-Poisson bracket on s * .Finally, let {•, •} g * denote the Lie-Poisson bracket on g * .The central result due to Scovel-Weinstein is the following.When (B × b * ) × s * is equipped with the evident product Poisson manifold structure, the Scovel-Weinstein map Γ is a Poisson map.In other words, the direct product Poisson bracket is compatible with the Lie-Poisson bracket on g according to

6 / 18
Lie algebra homomorphism.If s 1 + I , s 2 + I ∈ s I are elements of the quotient space their Lie bracket is[s 1 + I , s 2 + I ] = [s 1 , s 2 ] + I .Scovel-Weinstein highlight that since Π I is a Lie algebra homomorphism, the dual map Π * I : s * I → s * is a Poisson map.Therefore the map id × Π * I : (B × b * ) × s * I → (B × b * ) × s * , is Poisson as well, relative to the product Poisson manifold structures on the source and target spaces.It follows that there is a second Poisson map related to Γ via the quotient algebra construction.The reduced Scovel-Weinstein map Γ I is the Poisson map Γ I : (B × b * ) × s * I → g * given as the composition Γ I = Γ • [id × Π * I ].Theorem 3 below represents a restatement of the Poisson map property for Γ I within the fluid moment context.
), there is a unique Lie algebra structure on comoment vectors C m 0

Remark 2 .
The right-hand-side of Eq.(11) recovers the ansatz (4) for the novel exact reduction of Vlasov-Poisson dynamics mentioned in the introduction.The fact that the map(11) is Poisson justifies the claim from Section 4 that the bracket {•, •} m 0 repairs the shortcoming of the Hamiltonian fluid moment closures due to Holm and Tronci22 that do not include the density moment.It is interesting to note that Chong40 previously constructed a Poisson map that sends phase space distributions into fluid-like quantities; that Poisson map appears to be unrelated to this one.
) denote the Hamiltonian for the Vlasov-Poisson system.Here V : C * (Q) → C(Q) is the self-adjoint linear operator that sends a density n to the corresponding solution V of Poisson's equation, −∆V = 4πe n, where n denotes the number density less its average.Pulling back H along the Poisson map Γ m 0 produces a collective Hamiltonian H m 0 on P m 0 = C(Q) ×C * (Q) × M m 0
of Perin et al. that also found Hamiltonian closures with k < 5, but without stability guarantees.The feasibility of finding Jacobi closures was first established 26 by Perin et al. through a Jacobi-closure derivation of the earliest Hamiltonian fluid model due to Morrison and Greene general Hamiltonian vector field X H = ( ψH , ṅH , Ṗ P P 0 H ) is given by 2dq.