Higher Order Geometric Theory of Information and Heat Based on Poly-Symplectic Geometry of Souriau Lie Groups Thermodynamics and Their Contextures: The Bedrock for Lie Group Machine Learning

We introduce poly-symplectic extension of Souriau Lie groups thermodynamics based on higher-order model of statistical physics introduced by Ingarden. This extended model could be used for small data analytics and machine learning on Lie groups. Souriau geometric theory of heat is well adapted to describe density of probability (maximum entropy Gibbs density) of data living on groups or on homogeneous manifolds. For small data analytics (rarified gases, sparse statistical surveys, …), the density of maximum entropy should consider higher order moments constraints (Gibbs density is not only defined by first moment but fluctuations request 2nd order and higher moments) as introduced by Ingarden. We use a poly-sympletic model introduced by Christian Günther, replacing the symplectic form by a vector-valued form. The poly-symplectic approach generalizes the Noether theorem, the existence of moment mappings, the Lie algebra structure of the space of currents, the (non-)equivariant cohomology and the classification of G-homogeneous systems. The formalism is covariant, i.e., no special coordinates or coordinate systems on the parameter space are used to construct the Hamiltonian equations. We underline the contextures of these models, and the process to build these generic structures. We also introduce a more synthetic Koszul definition of Fisher Metric, based on the Souriau model, that we name Souriau-Fisher metric. This Lie groups thermodynamics is the bedrock for Lie group machine learning providing a full covariant maximum entropy Gibbs density based on representation theory (symplectic structure of coadjoint orbits for Souriau non-equivariant model associated to a class of co-homology).


Introduction
These two Pierre Duhem's citations (see [4] for English translations) make reference to Aristotle definition of "motion" (which can be found in The Physics) to designate not only a change of position in space, but also any change of state, even if not accompanied by any displacement. In this case, dynamics appears as a special case of General Thermodynamics [2,3,5], to describe in common principles all changes in the state of the body, both changes of place and changes in physical qualities. Making reference to Duhem's "Energetics", Stefano Bordini write in [6]: "This theoretical design led Duhem to rediscover and reinterpret the tradition of Aristotle's natural philosophy and Pascal's epistemology . . . This outcome was surprising and clearly echoed the Aristotelian language and concept of motion as change and transformation: within the framework of Aristotelian natural philosophy, motion in the modern physical sense was actually a special case of the general concept of motion. The mathematisation of thermodynamics coincided with a generalisation of mechanics, and this generalisation led to an unexpected connection between modern mathematical physics and ancient natural philosophy" (see [7,8] for more developments on the affiliation between Aristotle, Pascal and Duhem philosophies). This conceptual and epistemology point of view was enlightened 75 years later by Jean-Marie Souriau through the symplectic model of geometric mechanics applied to statistical mechanics and used to build a "Lie groups thermodynamics" of dynamical systems, where the Gibbs density is covariant with respect to the action of the Lie group on the system (dynamical groups as Galileo group). This Souriau theory is based on tools related to non-equivariant model associated to a class of co-homology and affine representation of Lie groups and Lie algebra (last approach was independently studied in mathematical domain by Koszul to characterize homogeneous convex cones geometry [9][10][11]). Duhem [12] and Souriau [13,14] also both studied how to extend Thermodynamics for a continuous media.
In this paper, we will explore and compare the joint geometric contextures shared in information theory (based on Koszul's information geometry) and heat theory (based on Souriau's Lie groups thermodynamics) to highlight their joint elementary structures. Classically, we address analogies between mathematical or physical models by comparing their "structures" defined as the arrangement of and relations between the parts or elements, or as the way in which the parts are arranged or organized. My personal concept of "contexture" is more general and phenomenological and could be defined as the act, process, or manner of weaving parts into a whole. We have then replaced the relations between objects by the act to build these relations. Based on Souriau's general definition of entropy as the Legendre transform of the logarithm of generalized Laplace transform and symplectic structure associated to Lie group coadjoint orbits, we will see how geometric structures of information and heat theories are generated by these Souriau's "generative processes". We will extend theses contextures in the vector-valued case based on poly-symplectic model of higher order Souriau's Lie groups thermodynamics.
In this paper, we identify the Riemanian metric introduced by Souriau based on co-homology, in the framework of "Lie groups thermodynamics" as an extension of classical Fisher metric introduced in information geometry. We have observed that Souriau metric preserves Fisher metric structure as the Hessian of the minus logarithm of a partition function, where the partition function is defined as a generalized Laplace transform on a convex cone. Souriau's definition of Fisher metric extends the classical one in case of Lie groups or homogeneous manifolds. Souriau has developed "Lie groups thermodynamics" in the framework of homogeneous symplectic manifolds in geometric statistical mechanics for dynamical systems, but as observed by Souriau, these model equations are no longer linked to the symplectic manifold but only depend on the Lie group and the associated co-cycle.
This analogy with Fisher metric opens potential applications in machine learning, where the Fisher metric is used by information geometry, to define the "natural gradient" tool to improve ordinary stochastic gradient descent sensitivity to rescaling or changes of variable in parameter space [15][16][17][18][19][20][21][22]. In machine learning revised by natural gradient of information geometry, the ordinary gradient is designed to integrate the Fisher matrix. Amari has theoretically proved the asymptotic optimality of the natural gradient compared to classical gradient. With the Souriau approach, the Fisher metric could be extended, by Souriau-Fisher metric, to design natural gradients for data on homogeneous manifolds. Information geometry has been derived from invariant geometrical structure involved in statistical inference. The Fisher metric defines a Riemannian metric as the Hessian of two dual potential functions, linked to dually coupled affine connections in a manifold of probability distributions. With the Souriau model, this structure is extended preserving the Legendre transform between two dual potential function parametrized in Lie algebra of the group acting transentively on the homogeneous manifold.
Classically, to optimize the parameter θ of a probabilistic model, based on a sequence of observations y t , is an online gradient descent: with learning rate η t , and the loss function l t = − log p(y t /ŷ t ). This simple gradient descent has a first drawback of using the same non-adaptive learning rate for all parameter components, and a second drawback of non invariance with respect to parameter re-encoding inducing different learning rates. Amari has introduced the natural gradient to preserve this invariance to be insensitive to the characteristic scale of each parameter direction. The gradient descent could be corrected by I(θ) −1 where I is the Fisher information matrix with respect to parameter θ, given by: with natural gradient: Amari has proved that the Riemannian metric in an exponential family is the Fisher information matrix defined by: and the dual potential, the Shannon entropy, is given by the Legendre transform: In geometric statistical mechanics, Souriau has developed a "Lie groups thermodynamics" of dynamical systems where the (maximum entropy) Gibbs density is covariant with respect to the action of the Lie group. In the Souriau model, previous structures of information geometry are preserved: S(Q) = β, Q − Φ(β) with Q = ∂Φ(β) ∂β ∈ g * and β = ∂S(Q) ∂Q ∈ g In the Souriau Lie groups thermodynamics model, β is a "geometric" (Planck) temperature, element of Lie algebra g of the group, and Q is a "geometric" heat, element of dual Lie algebra g * of Entropy 2018, 20, 840 4 of 60 the group. Souriau has proposed a Riemannian metric that we have identified as a generalization of the Fisher metric: The tensor Θ used to define this extended Fisher metric is defined by the moment map J(x), M (homogeneous symplectic manifold) to the dual Lie algebra g * , given by: This tensor Θ is also defined in tangent space of the cocycle θ(g) ∈ g * (this cocycle appears due to the non-equivariance of the coadjoint operator Ad * g , action of the group on the dual lie algebra): Q Ad g (β) = Ad * g (Q) + θ(g) (11) Θ(X, Y) : g × g → with Θ(X) = T e θ(X(e)) X, Y → Θ(X), Y In Souriau's Lie groups thermodynamics, the invariance by re-parameterization in information geometry has been replaced by invariance with respect to the action of the group. When an element of the group g acts on the element β ∈ g of the Lie algebra, given by adjoint operator Ad g . Under the action of the group Ad g (β), the entropy S(Q) and the Fisher metric I(β) are invariant: β ∈ g → Ad g (β) ⇒ S Q Ad g (β) = S(Q) I Ad g (β) = I(β) (13) In the case of small data analytics, we propose to parameterized the (maximum entropy) Gibbs density with higher order "geometric" temperature β k and higher order heat Q k , that parameterized higher order entropy S(Q 1 , . . . , Q n ) and dual potential function Φ(β 1 , . . . , β n ): We will develop in the paper that the geometric approach of statistical thermodynamics, introduced by Souriau, offers an advantage over traditional formulations. Classical thermodynamics has been developed for static systems taking into accound only the time evolution, but in case of dynamical systems (e.g., a centrifuge system), this statistical physics is no longer valid because the Gibbs density (the density of the maximum entropy) is not covariant. In case of only time translation, what is preserved is only the energy, but for dynamical systems where a group is acting, invariants are given by components of the "moment map" (which is a geometrization of the Noether theorem providing invariants if there are symmetries). The "moment map" has been introduced in parallel by Kostant in mathematics and by Souriau in physics. Souriau has developed the non-equivariant case, and has applied it to statistical mechanics. The main advantages of "Lie groups thermodynamics" of dynamical systems, is that this statistical physics is a coordinate-free model preserving invariances with respect to the action of the (dynamical) Lie group acting on the system. We give in appendix the development of the centrifuge thermodynamics with classical approach given by Roger Balian, and prove that with the Souriau approach, the problem is solved by only applying Lie groups thermodunamics equations through moment map computation, where in classical case, we The geometric definition and extension of Fisher metric has been recently studied in the framework of quantum information geometry, but this community seems unaware of Souriau's work on Lie groups thermodynamics for the study of statistical physics of dynamical systems based on symplectic geometry and c-homology tools in the 70s, and in particular the non-equivariant case developed by Souriau and Koszul. We can make reference to the following recent works on the symplectic formulation of the Fisher information theory [56][57][58][59].
The structure of the paper is the following: • In Section 1, we introduce seminal idea on Symplectic geometry used in mechanics and in statistical mechanics, as introduced by Jean-Marie-Souriau during the 60s. From previous work of François Gallissot extending Cartan's results on integral invariant (theorem on types of differential forms generating equations of movement of a material point invariant in the transformations of the Galilean group), we present the Lagrange 2-form and moment map elaborated by Souriau to build a geometric mechanics theory, where a dynamical system is then represented by a foliation of the evolution, determined by an antisymmetric covariant second order tensor. Souriau has applied this tool for mechanical statistics to build a thermodynamics of dynamical systems, where the classical notion of Gibbs canonical ensemble is extended for a homogeneous symplectic manifold on which a Lie group (dynamical group) has a symplectic action. In case of Galileo group, the symmetry is broken, and new "co-homological" Souriau relations should be verified in Lie algebra of the group.

•
In Section 2, we synthetize results on higher order thermodynamics based on higher order temperatures and heats, as introduced by Ingarden and Jaworski for mesoscopic systems. This model is based on higher order maximum entropy Gibbs density definition constraining solution with respects to higher order moments.

•
In Section 3, we develop "Lie groups thermodynamics" model, developed to describe Gibbs state for dynamical systems, where Souriau introduced the concept of co-adjoint action of a group on its momentum space that allows designing physical observables like energy, heat and momentum or moment as pure geometrical objects. The Souriau model then generalizes the Gibbs equilibrium state to all symplectic manifolds that have a dynamical group, with a "geometric" (Planck) temperature as an element of the Lie algebra and "geometric heat" as an element of the dual Lie algebra. We have observed that Souriau has introduced a symmetric tensor that is an extension of classical Fisher metric in information geometry. This new Fisher-Souriau metric is invariant with respect to the action of the group. These equations are universal, because they are not dependent on the symplectic manifold but only on the dynamical group and its associated two-cocycle. Souriau  and geometric developments of geometric mechanics. Koszul has defined a skew symmetric bilinear form by a closed expression depending only on the cocycle and related to the Souriau antisymmetric bilinear map introduced previously in Section 3. This Koszul study of the moment map non-equivariance, and the existence of an affine action of G on g* is at the cornerstone of Souriau theory of Lie groups thermodynamics.

•
In Section 5, at the step of the Souriau Lie groups thermodynamics presentation, we will introduce a generalized Souriau definition of entropy, as the Legendre transform of the logarithm of the Laplace transform, making the connection with information geometry. This definition is a general contexture that can be extended to highly abstract spaces preserving Legendre structure, if we are able to generalize the Laplace transform.

•
In Section 6, we illustrate Souriau's Lie groups thermodynamics for a centrifuge system. The main Souriau idea was to define the Gibbs states for one-parameter subgroups of the Galilean group, because he proved that the action of the full Galilean group on the space of motions of an isolated mechanical system is not related to any equilibrium Gibbs state (the open subset of the Lie algebra, associated to this Gibbs state, is empty).

•
In Section 7, we have defined an higher-order model of Lie groups thermodynamics based on a poly-symplectic vector valued approach. This multi-symplectic extension, is based on a multi-valued one that preserve the notion of (poly-)moment map built by Günther based on an n-symplectic model. We replace the symplectic form of the Souriau model by a vector valued form that is called poly-symplectic. We consider the non-equivariance of poly-moment map by introducing poly-cocycle. We finally conclude with poly-symplectic definition extension of the Fisher-Souriau metric.

•
In Section 8, we conclude with potential extension to Lie group machine learning.
To facilitate understanding of previous results, we add some additional complements: • In Appendix A, we recall a synthesis of Günther's poly-symplectic model with initial notation • In Appendix B, we develop computation of the Fisher metric for multivariate Gaussian density, to establish links with Souriau's Lie groups Gibbs density model.

•
In Appendix C, we give more details on the Legendre transform, the basic tool of information geometry and Souriau Lie groups thermodynamics. More especially, we give a definition of the Legendre transform with projective geometry definition by Chasles as reciprocal polar with respect to a paraboloid.

•
In Appendix D, we give solution of a centrifuge system thermodynamics, given by Roger Balian based on a classical approach, to make the link with the Souriau approach.

•
In Appendix E, we recall the main proofs of Souriau's Lie groups thermodynamics and its poly-symplectic extension.

•
In Appendix F, we present another Souriau statistical physics model, developed for relativistic thermodynamics of continua, which preserves the Legendre transform, where temperature is given by a killing vector.

Seminal Idea of Symplectic Geometry in Mechanics and in Statistical Mechanics by Gallissot and Souriau
The symplectic structure has been introduced in mathematics much earlier than the word symplectic, in works of the French physicist Joseph Louis Lagrange (see paper on the slow changes of the orbital elements of planets in the solar system), who showed that this geometry is a fundamental tool in the mathematical model of any problem in mechanics. Jean-Marie Souriau has shown that Lagrange's parentheses (nowdays called Lagranges bracket) are the components of the canonical symplectic 2-form on the manifold of motions of the mechanical system, in the chart of that manifold [60,61]. As early as 1966, influenced by François Gallissot's work, Souriau applied his theory of geometric mechanics to statistical mechanics, developed in Chapter IV of his book "Structure of Dynamical Systems" [43,64], what he called "Lie groups thermodynamics". We have discovered that Souriau and Galissot both attended the 1954 International Congress of Mathematicians (ICM'54) in Moscow. We could assume that they have discussed 1952 Gallisot's paper introducing three types of differential forms generating equations of movement of a material point invariant in the transformations of the Galilean group and their links with Poincaré-Cartan integral invariant. This seminal work of Gallisot helped Souriau to formulate his new geometric mechanics and its extenxion to geometric statistical physics. Using Lagrange's viewpoint, in Souriau statistical mechanics, a statistical state is a probability measure on the manifold of motions. As we can read in his book, Souriau was influenced by François Gallisot to introduce the Lagrange(-Souriau) 2-form.
In place of classical mechanical equations of a material point subjected to a force F, defined by its mass m and its position r at time t, the second order differential equations As early as 1966, influenced by François Gallissot's work, Souriau applied his theory of geometric mechanics to statistical mechanics, developed in Chapter IV of his book "Structure of Dynamical Systems" [43,64], what he called "Lie groups thermodynamics". We have discovered that Souriau and Gallissot both attended the 1954 International Congress of Mathematicians (ICM'54) in Moscow. We could assume that they have discussed 1952 Gallissot's paper introducing three types of differential forms generating equations of movement of a material point invariant in the transformations of the Galilean group and their links with Poincaré-Cartan integral invariant. This seminal work of Gallissot helped Souriau to formulate his new geometric mechanics and its extenxion to geometric statistical physics. Using Lagrange's viewpoint, in Souriau statistical mechanics, a statistical state is a probability Entropy 2018, 20, 840 8 of 60 measure on the manifold of motions. As we can read in his book, Souriau was influenced by François Gallissot to introduce the Lagrange(-Souriau) 2-form.
In place of classical mechanical equations of a material point subjected to a force F, defined by its mass m and its position r at time t, the second order differential equations m d 2 r dt 2 = F is rewritten by a system of first order differential equations in phase space r v : If the force F is derived from a potential w, we have classical equations: Ldt and Hamilton-Jacobi equations This idea of Lagrange, rediscovered by Souriau was to consider time t like the others variables. One should use then the 7-dimensional space V (evolution space) (see Figure 2): Classical system of first order differential equations in phase space can then be rewritten in evolution space V by the homogeneous form: At each point y of V, these equations define the tangent direction to the curve x described by the point y during the evolution of the system. These curves are the leaves (lines of force) of the field of directions defined by the equations of the homogeneous form, as defined for foliated manifolds. See [43], for more details on definition of the different derivatives used.
Entropy 2018, 20, x FOR PEER REVIEW 9 of 62 Classical system of first order differential equations in phase space can then be rewritten in evolution space V by the homogeneous form: At each point y of V, these equations define the tangent direction to the curve x described by the point y during the evolution of the system. These curves are the leaves (lines of force) of the field of directions defined by the equations of the homogeneous form, as defined for foliated manifolds. See [43], for more details on definition of the different derivatives used.  σ is considered as a bilinear operator on tangent vectors of V. If we choose two such vectors: σ associates to them an antisymmetric scalar product: In the Souriau-Lagrange model, σ is a 2-form on the evolution space V, and the differential equation of motion δy ∈ ε implies: σ(δy)(δ y) = 0 , ∀δ y which can be written as: For study of this Souriau-Lagrange 2-form, readers should see the papers of Obădeanu [65][66][67]. Souriau has observed that this 2-form was introduced by Lagrange in a different language in his study of celestial mechanics in 1808. Souriau was also influenced by François Gallissot that used this 2-form in [68,69]. We will see in the following the Souriau's "moment map µ" in dual Lie algebra of the group G, and the study of coadjoint orbits of G. For the definition of moment map, we make reference to [45]. Souriau has extended this model for thermodynamics. For this new phenomenological approach of mechanics, thermodynamics and information theory, we can give reference to Souriau introduction of his paper "Quantique? Alors c'est géométrique" [70] and a video of his talk [71]: "Plaçons-nous d'abord dans le cadre de la mécanique classique. Étudions un système mécanique isolé, non dissipatif-nous dirons brièvement une «chose». L'ensemble des mouvements de cette «chose» est une variété symplectique. Pourquoi? Il suffit de se reporter à la Mécanique Analytique de Lagrange (1811); l'espace des mouvements y est traité comme variété différentiable; les coordonnées covariantes et contravariantes de la forme symplectique y sont écrites (Ce sont les "parenthèses" et "crochets" de Lagrange). Évoquons maintenant la géométrie du 20 éme siècle. Soit G un groupe difféologique (par exemple un groupe de Lie); µ un moment de G (un moment, c'est une 1-forme invariante à gauche sur G); alors l'action du groupe sur µ engendre canoniquement un espace symplectique (ces groupes pourront avoir une dimension infinie). Présomption épistémologique: derrière chaque «chose» est caché un groupe G (sa "source"), et les mouvements de la «chose» sont simplement des moments de G (doublet latin mnémotechnique : momentum-movimentum). L'isolement de la «chose» indique alors que le groupe de Poincaré (respectivement de Galilée-Bargman) est inséré dans G; voilà l'origine des grandeurs conservées relativistes (respectivement classiques) associées à un mouvement x: elles constituent simplement le moment induit sur le groupe spacio-temporel par le moment-mouvement x." (In English: Let's put ourselves first in the framework of classical mechanics. Let's study an isolated, non-dissipative mechanical system-we will briefly say a "thing". The set of movements of this "thing" is a symplectic manifold. Why? It is enough to refer to the Analytical Mechanics of Lagrange (1811); the space of movements is treated as a differentiable manifold; the covariant and contravariant coordinates of the symplectic form are written there (these are the "parentheses" and "brackets" of Lagrange). Let's now talk about the geometry of the 20th century. Let G be a diffeological group (for example a Lie group); µ a moment of G (a moment is a left invariant 1-form on G); then the action of the group on µ canonically generates a symplectic space (these groups can have an infinite dimension). Epistemological presumption: behind each "thing" is hidden a group G (its "source"), and the movements of the "thing" are simply moments of G (mnemonic Latin doublet: momentum-movimentum). The isolation of the "thing" then indicates that the group of Poincaré (respectively Galileo-Bargman) is inserted in G; here is the origin of the relativistic (respectively classical) conserved magnitudes associated with a movement x: they simply constitute the moment induced on the spacio-temporal group by the moment-motion x.) "Il y a un théorème qui remonte au XXème siècle. Si on prend une orbite coadjointe d'un groupe de Lie, elle est pourvue d'une structure symplectique. Voici un algorithme pour produire des variétés symplectiques: prendre des orbites coadjointes d'un groupe. Donc cela laisse penser que derrière cette structure symplectique de Lagrange, il y avait un groupe caché. Prenons le mouvement classique d'un moment du groupe, alors ce groupe est très «gros» pour avoir tout le système solaire. Mais dans ce groupe est inclus le groupe de Galilée, et tout moment d'un groupe engendre des moments d'un sous-groupe. On va retrouver comme cela les moments du groupe de Galilée, et si on veut de la mécanique relativiste, cela va être du groupe de Poincaré. En fait avec le groupe de Galilée, il y a un petit problème, ce ne sont pas les moments du groupe de Galilée qu'on utilise, ce sont les moments d'une extension centrale du groupe de Galilée, qui s'appelle le groupe de Bargman, et qui est de dimension 11. C'est à cause de cette extension, qu'il y a cette fameuse constante arbitraire figurant dans l'énergie. Par contre quand on fait de la relativité restreinte, on prend le groupe de Poincaré et il n'y a plus de problèmes car parmi les moments il y a la masse et l'énergie c'est mc 2 . Donc le groupe de dimension 11 est un artéfact qui disparait, quand on fait de la relativité restreinte." (In Engish: There is a theorem dating back to the twentieth century. If we take a coadjoint orbit of a Lie group, it is provided with a symplectic structure. Here is an algorithm to produce symplectic manifolds: take coadjoint orbits from a group. So it suggests that behind this symplectic structure of Lagrange, there was a hidden group. Take the classic movement of a moment of the group, so this group is very "big" to have the whole solar system. But in this group is included the Galileo group, and any moment of a group generates moments of a subgroup. We will find like that the moments of the group of Galileo, and if we want relativistic mechanics, it will be Poincaré group. In fact with Galileo group, there is a small problem, it is not the moments of the Galileo group that are used, it is the moments of a central extension of the Galileo group, which is called the Bargman group, and that is of dimension 11. It is because of this extension, that there is this famous arbitrary constant appearing in the energy. On the other hand, when we do special relativity, we take Poincaré group and there are no more problems because among the moments there is the mass and the energy is mc 2 . So the 11-dimensional group is an artifact that disappears, when we do special relativity.) François Gallissot has observed that in his famous lessons on integral invariants, Elie Cartan has shown that all the properties of the differential equations of the dynamics of holonomic systems result from the existence of the integral invariant: Thus every holonomic system whose forces derive from a force function is associated to a form ω, the equations of motion being the characteristics of the exterior form dω. Around 1950, the theory of exterior forms on differentiable manifolds has been established on new foundations under the influence of topologists. The question was then to wonder: • if classical mechanics cannot benefit from these models by placing an exterior form of degree two at its base • if thanks to the notion of manifold, the notion of connection cannot be introduced in a more natural way • if the paradoxal indeterminations/impossibilities in the Lagrangian framework could be explained more clearly • if the problem of integration of equations of motion could be enlightened, generated by a form Ω of degree two.
To reach these various objectives, Gallissot has resumed first the study of the logical bases on which the Galilean mechanics is built. He thus shown that when it is proposed to find generating forms of the equations of motion of a material invariant point in the transformations of the Galilean group, the most interesting form is an exterior form of degree two defined on a variety E 3 × E × T (E 3 Euclidean space, T temporal). Gallissot had shown that any holonomic parametric system with n degrees of freedom is associated with a form Ω of degree 2n defined on a differentiable manifold whose characteristics are the equations of the movement. This form is expressed by means of 2n Pfaff forms and by dt, the Hamiltonian form being a simple special case. He gave a summary of how we can get rid of the servitude of coordinates in the study of dynamical systems and the important role played by the operator i() antiderivative introduced by Cartan, the characteristic field E of the form Ω being defined by the relation i(E)Ω = 0. Gallissot has then introduced the following theorem: There are three types of differential forms generating equations of movement of a material point invariant in the transformations of the Galilean group: A : If we consider the last form "C": dω = 0 constraints Pfaff form δ ij F i dx j to be closed, and to reduce the differential of function U: with H = T − U and T = 1 2 It proves that the exterior derivative of ω is: The form ω * = dω generates Elie Cartan integral invariant. In Chapter IV of his book, Souriau applied this model based on symplectic geometry for statistical mechanics. Souriau observed that Gibbs equilibrium is not covariant with respect to dynamic groups of physics. To solve this breaking of symmetry, Souriau introduced a new "geometric theory of heat" where the equilibrium states are indexed by a parameter β with values in the Lie algebra of the group, generalizing the Gibbs equilibrium states, where β plays the role of a geometric (Planck) temperature. Souriau observed that the group of time translations of the classical thermodynamics is not a normal subgroup of the Galileo group, proving that if a dynamical system is conservative in an inertial reference frame, it need not be conservative in another. Based on this fact, Souriau generalized the formulation of the Gibbs principle to become compatible with Galileo's relativity in classical mechanics and with Poincaré relativity in relativistic mechanics. The maximum entropy principle is preserved, and the Gibbs density is given by the density of maximum entropy (among the equilibrium states for which the average value of the energy takes a prescribed value, the Gibbs measures are those which have the largest entropy), but with a new principle "If a dynamical system is invariant under a Lie subgroup G' of the Galileo group, then the natural equilibria of the system forms the Gibbs ensemble of the dynamical group G'". The classical notion of Gibbs canonical ensemble is extended for a homogneous symplectic manifold on which a Lie group (dynamic group) has a symplectic action. In case of a Galileo group, the symmetry is broken, and new "cohomological" relations should be verified in Lie algebra of the group. A natural equilibrium state will thus be characterized by an element of the Lie algebra of the Lie group, determining the equilibrium temperature β. The entropy s(Q), parametrized by Q the geometric heat (mean of energy U, element of the dual Lie algebra) is defined by the Legendre transform of the Massieu potential given by Φ(β), parametrized by β (Φ(β) is the minus logarithm of the partition function ψ Ω (β)): Souriau has proposed to study the statistical mechanics from the new point of view of symplectic geometry, completing the work of Poincaré and Cartan on integral invariant, reinventing the Lagrangian symplectic form in place of classical variational formulation and geometrizing the Noether theorem with a moment map as new conserved quantities. Firstly, Souriau Lie groups thermodynamics gives geometrical status to the (Planck) temperature and the entropy with a new general definition of the Fisher Metric. Secondly, Souriau's relativistic thermodynamics of continua provides a geometrization of the smecond principle by the permanence of the entropy current, whose flux has positive divergence [13,14,[72][73][74]. This 2nd model of Souriau's thermodynamics is described in the Appendix. Other authors have studied this relativistic thermodynamics of continua [75][76][77][78][79][80][81][82].

Higher Order Thermodynamics Based on Higher Order Temperatures
We will generalize Souriau's theory [43,64], reconsidered in [23] and with links to information geometry in [91], in the framework of higher order thermodynamics as introduced by Ingarden [29][30][31] and Jaworski [32][33][34][35] for mesoscopic systems. We can make also reference to other publications of Ingarden [36][37][38][39][40], Jaworsky [92][93][94] and Nakagomi [95] on higher order thermodynamics. The Gibbs canonical state results from the maximum entropy principle when the statistical mean value of the energy is supposed to be known. A Polish school has studied the maximum entropy inference with higher-order moments of energy (when not only mean values but also statistical moments of higher order of some physical quantities are taken into account). Ingarden in 1963 and Jaworski in 1981 have introduced the concept of second and higher-order temperatures, by assuming a distribution function which includes information not only on the average of the energy but also on higher-order moments, in particular 2nd moment related to fluctuations. This case should be considered in situations where fluctuations are not negligible, such as near phase transitions or critical points, in metastable states in systems with a small number of degrees of freedom. Ingarden's idea is that if we can measure more details, such as the first n cumulants of the energy, we can then introduce n high-order temperature, as the Lagrange multipliers when we maximize the entropy with respect to these values: Ingarden proposed that if we can measure the second cumulant of the energy (the fluctuation of the energy), the equilibrium state is not the canonical state, but would need two temperatures.
Ingarden argues that for a macroscopic system there is very little difference between the two states, and that we would need a mesoscopic or microscopic system to be able to detect the higher temperature. Jaworski [27,28] has shown that the contribution to the total entropy, arising from the extra information corresponding to the higher-order moments, is o(N), when N tends to infinity and N/V ratio is constant, with N the number of particles and V the volume. The main result of Jaworski is that from a purely thermodynamic point of view, the information corresponding to the higher-order moments of extensive physical quantities is not essential and can be neglected in the maximum entropy procedure. Jaworski showed that the maximum entropy inference has a certain stability property with respect to information corresponding to higher order moments of extensive quantities. It can serve as an argument in favor of the maximum entropy method in statistical physics and to understand better why these methods are successful. Streater [96] has prefered to say that the states with generalized temperatures are not in equilibrium, assuming that the final state, at large times, will be the canonical or grand canonical state depending on mixing properties. Streater [96] intends that this occur even for a mesoscopic system, such as a few atoms, adding that his approach is equivalent to Ingarden model if the relaxation time from the state with generalized temperatures to the final equilibrium is very long.
Some examples of higher order maximum Entropy are given by Ingarden: • 1st Example of Higher Oder Maximum Entropy Density: Density of maximum Entropy under the constraints: is given by: with the following parameters: where: We illustrate this higher order maximum entropy density in Figure 3. where: (35) We illustrate this higher order maximum entropy density in Figure 3. • 2nd Example of Higher Oder Maximum Entropy Density: under the constraints: is given by: with the following parameters:  (38) where: is given by: with the following parameters: where: We illustrate this higher order maximum entropy density in Figure 4. We illustrate this higher order maximum entropy density in Figure 4. As soon as 1963, Ingarden has introduced this concept of higher order temperatures for statistical systems such as thermodynamics. In physics, the concept of temperature is connected with the mean value of kinetic energy of molecules in an ideal gas. For a general physical system with interactions among particles (the case of non-ideal gas: liquid or solid), an equilibrium probability distribution depends on temperature T as the only statistical parameter of the Gibbs state:  As soon as 1963, Ingarden has introduced this concept of higher order temperatures for statistical systems such as thermodynamics. In physics, the concept of temperature is connected with the mean value of kinetic energy of molecules in an ideal gas. For a general physical system with interactions among particles (the case of non-ideal gas: liquid or solid), an equilibrium probability distribution depends on temperature T as the only statistical parameter of the Gibbs state: where p is position, q the mechanical momentum and k β the Boltzmann constant (a factor to insure that β.H is dimensionless). If there are no stochastic interactions between particles (ideal gas), the partition function Z has the property to be integrable and we can obtain Gauss distribution in the momentum space deduced from the result of the limit theorem for large N. The ideal gas model of Boltzmann can fail if the number of particles is not large enough in the case of mesoscopic systems, and also if the interactions between particles are not weak enough. Gibbs hypothesis can also fail in other cases when stochastic interactions with the environment are not sufficiently weak. As remarked by Ingarden, nobody has ever observed thermal Gibbs equilibrium in large and complex systems (cosmic systems, Earth's atmosphere, biological organisms), but only in cases of turbulence, flows or pumping, by replacing classical approach by local temperature and concept of thermodynamic flows (non-equilibrium thermodynamics and thermo-hydrodynamics), that is non-coherent with the classical concept of temperature which is, by definition, global/intensive and does not depend on position. R.S. Ingarden proposed to consider the stationary case using of the concept of higher order temperatures given by the Gibbs density: where U = E(H) is the mean energy. This mean energy has been introduced to preserve the the total energy invariance with respect to an arbitrary additive constant, and β 0 = − log Z(β 1 , . . . , β n ) the constant of normalization. The new constants β k are said to be β-temperatures of order k. H(x) is usually defined as a quadratic function of x. The probability distribution is uniquely defined from statistical moments which should be measured experimentally. But if values number is too high to make this method practical, we are only able to measure the lowest moments up to some order (if we can neglect the higher orders that do not change the result to a given accuracy), and to fix β-temperatures defined as Lagrange multipliers by maximization of entropy of distribution S = − P (β 1 ,...,β n ) (x) log P (β 1 ,...,β n ) (x)dx, with the given moments as constraints. R.S. Ingarden observed that the entropy maximization randomizes higher moments in a symmetric way, and it cancel any possible bias with respect to their special values, and it gives the best estimate to a given accuracy. The values of β can be found by: β k x k dx and the relation : Ingarden has applied this model for linguistic statistics, assuming the appearance of higher order temperatures since there occur rather strong statistical correlations between phonemes and words as elements of these statistics. He argued his choice observing that in the case of word statistics, the existence of strong correlations is given by grammatical or semantical studies [9]. Ingarden made the conjecture that his high order thermodynamics is the model of statistically interacting, biological living systems, and small systems although the calculation/observation are more difficult.
Ingarden higher order temperatures could be defined in the case when no variation is considered, but when a probability distribution depending on more than one parameter. It has been observed by Ingarden, that Gibbs assumption can fail if the number of components of the sum goes to infinity and the components of the sum are stochastically independent, and if stochastic interactions with the environment are not sufficiently weak. In all these cases, we never observe absolute thermal equilibrium of Gibbs type but only flows or turbulence. Non-equilibrium thermodynamics could be indirectly addressed by means of high order temperatures.

Model of Souriau Lie Groups Thermodynamics
For introduction to symplectic geometry, we make reference to Marle's book [45] and Koszul's book [44]. In 1969, Souriau [43,64] introduced the concept of co-adjoint action of a group on its momentum space, based on the orbit method works, that allows to define physical observables like energy, heat and momentum or moment as pure geometrical objects. The moment map is a constant of the motion and is associated to symplectic cohomology. In a first step to establish new foundations of thermodynamics, Souriau has defined a Gibbs canonical ensemble on a symplectic manifold M for a Lie group action on M. In classical statistical mechanics, a state is given by the solution of Liouville equation on the phase space, the partition function. As symplectic manifolds have a completely continuous measure, invariant by diffeomorphisms, the Liouville measure λ, all statistical states will be the product of the Liouville measure by the scalar function given by the generalized partition function e Φ(β)− β,U(ξ) defined by the energy U (defined in the dual of the Lie algebra of this dynamical group) and the geometric temperature β, where Φ is a normalizing constant such the mass of probability is equal to 1, Souriau then generalizes the Gibbs equilibrium state to all symplectic manifolds that have a dynamical group. Souriau has observed that if we apply this theory for a Galileo group, the symmetry has been broken. For each temperature β, element of the Lie algebra g, Souriau has introduced a tensor Θ β , equal to the sum of the cocycle Θ and the heat coboundary (with [.,.] Lie bracket): This tensor Θ β has the following properties: Θ(X, Y) = Θ(X), Y where the map Θ is the symplectic one-cocycle of the Lie algebra g with values in g * , with Θ(X) = T e θ(X(e)) where θ the one-cocycle of the Lie group G. Θ(X, Y) is constant on M and the map Θ(X, Y) : g × g → is a skew-symmetric bilinear form, and is called the symplectic two-cocycle of Lie algebra g associated to the moment map J, with the following properties: where J X linear application from g to differential function on M:g → C ∞ (M, R), X → J X and the associated differentiable application J, called moment(um) map: The geometric temperature, element of the algebra g, is in the the kernel of the tensor Θ β : The following symmetric tensor g β ([β, ] is positive definite, and defines extension of the classical Fisher metric in information geometry (as the Hessian of the logarithm of partition function): with: These equations are universal, because they are not dependent on the symplectic manifold but only on the dynamical group G, the symplectic two-cocycle Θ, the temperature β and the heat Q. Souriau called it "Lie groups thermodynamics" (see Figures 5 and 6). • Action of Lie group on Lie algebra: • Characteristic function after Lie group action: • Invariance of entropy with respect to action of Lie group: • Action of Lie group on geometric heat: Souriau's equations of Lie groups thermodynamics are summarized in the following figures.     In the framework of Lie group action on a symplectic manifold, equivariance of moment could be studied to prove that there is a unique action a(.,.) of the Lie group G on the dual * g of its Lie algebra for which the moment map J is equivariant, that means for each When the group is not abelian (non-commutative group), the symmetry is broken, and new "cohomological" relations should be verified in Lie algebra of the group. A natural equilibrium state will thus be characterized by an element of the Lie algebra of the Lie group, determining the equilibrium temperature β . The entropy In the framework of Lie group action on a symplectic manifold, equivariance of moment could be studied to prove that there is a unique action a(.,.) of the Lie group G on the dual g * of its Lie algebra for which the moment map J is equivariant, that means for each x ∈ M: When the group is not abelian (non-commutative group), the symmetry is broken, and new "cohomological" relations should be verified in Lie algebra of the group. A natural equilibrium state will thus be characterized by an element of the Lie algebra of the Lie group, determining the equilibrium temperature β. The entropy s(Q), parametrized by Q the geometric heat (mean of energy U, element of the dual Lie algebra) is defined by the Legendre transform [97][98][99][100][101][102][103] of the Massieu potential Φ(β) parametrized by β (Φ(β) is the minus logarithm of the partition function ψ Ω (β)): Souriau completed his "geometric heat theory" by introducing a 2-form in the Lie algebra, that is a Riemannian metric tensor in the values of adjoint orbit of β, [β, Z] with Z an element of the Lie algebra. This metric is given for (β, Q): where Θ is a cocycle of the Lie algebra, defined by Θ = T e θ with θ a cocycle of the Lie group defined by θ(M) = Q(Ad M (β)) − Ad * M Q. We observe that Souriau Riemannian metric, introduced with symplectic cocycle, is a generalization of the Fisher metric, that we call the Souriau-Fisher metric, that preserves the property to be defined as a Hessian of the partition function logarithm g β = − ∂ 2 Φ ∂β 2 = ∂ 2 log ψ Ω ∂β 2 as in classical information geometry. We will establish the equality of two terms, between Souriau definition based on Lie group cocycle Θ and parameterized by "geometric heat" Q (element of dual Lie algebra) and "geometric temperature" β (element of Lie algebra) and hessian of characteristic function Φ(β) = − log ψ Ω (β) with respect to the variable β: If we differentiate this relation of Souriau theorem Q Ad g (β) = Ad * g (Q) + θ(g), this relation occurs: As the entropy is defined by the Legendre transform of the characteristic function, this Souriau-Fisher metric is also equal to the inverse of the hessian of "geometric entropy" s(Q) with respect to the variable Q: For the maximum entropy density (Gibbs density), the following three terms coincide: that describes the convexity of the log-likelihood function, the Fisher metric that describes the covariance of the log-likelihood gradient, whereas I(β) = E (ξ − Q)(ξ − Q) T = Var(ξ) that describes the covariance of the observables.
We can also observe that the Fisher metric I(β) = − ∂Q ∂β is exactly the Souriau metric defined through symplectic cocycle: The Fisher metric has been considered by Souriau as a generalization of "heat capacity". Souriau called it K the "geometric capacity". We could observe that Souriau Lie groups thermodynamics is compatible with Balian and Valentin's theory of thermodynamics [24], that is obtained by symplectization in dimension 2n + 2 of contact manifold in dimension 2n + 1. All elements of the Souriau geometric temperature vector are multiplied by the same gauge parameter. The Balian and Valentin model was first explored in [104] and has been recently developed by der Schaft and Maschke in [26,105].

Extended Koszul Study of Souriau Non-Equivariant Model Associated to a Class of Cohomology
Koszul has deepened Souriau's model in his book "Introduction to symplectic geometry" [44] as explained in [10]. In the historical foreword of this book, Koszul write "The development of analytical mechanics provided the basic concepts of symplectic structures. The term symplectic structure is due largely to analytical mechanics. But in this book, the applications of symplectic structure theory to mechanics is not discussed in any detail". Koszul considers in this book purely algebraic and geometric developments of geometric/analytic mechanics developed during the 60s, more especially in Jean-Marie Souriau's works detailed in chapters 4 and 5. The originality of this book lies in the fact that Koszul develops new points of view, and demonstrations not considered initially by Souriau and after by the geometrical mechanics community.
To highlight the importance of this Koszul book, we will illustrate the links of the detailed tools, including demonstrations or original Koszul extensions, with Souriau's Lie groups thermodynamics. Koszul originally developed Souriau's model, in the case of non-equivariance, of the action of the group G on the moment map. As explained in [106] by Thomas Delzant at the 2010 CIRM conference "Action Hamiltoniennes: invariants et classification", organized with Michel Brion: "The definition of the moment map is due to Jean-Marie Souriau . . . . In the book of Souriau, we find a proof of the proposition: the map J is equivariant for an affine action of G on g* whose linear part is Ad* . . . . In Souriau's book, we can also find a study of the non-equivariant case and its applications to classical and quantum mechanics. In the case of the Galileo group operating in the phase space of space-time, obstruction to equivariance (a class of cohomology) is interpreted as the inert mass of the object under study". We can uniquely define the moment map up to an additive constant of integration, that can always be chosen to make the moment map equivariant (a moment map is G-equivariant, when G acts on g* via the coadjoint action) if the group is compact or semi-simple. In 1969, Souriau has considered the non-equivariant case where the coadjoint action must be modified to make the map equivariant by a 1-cocycle on the group with values in dual Lie algebra g*.
In this book in Chapter 4, Koszul calls symplectic G-space a symplectic manifold (M; ω) on which a Lie group G acts by a symplectic action (an action which leaves unchanged the symplectic form ω). Koszul then introduces and develop properties of the moment map µ (Souriau's invention) of a Hamiltonian action of the Lie algebra g. Koszul also defines the Souriau 2-cocycle, considering that the difference of two moments of the same Hamiltonian action is a locally constant application on M ,showing that when µ is a moment map, for every pair (a;b) of elements of g, the function c µ (a, b) = { µ, a , µ, b } − µ, {a, b} is locally constant on M, defining an antisymmetric bilinear application of gxg in H 0 (M; R) which verifies Jacobi's identity. This is the 2-cocycle introduced by Jean-Marie Souriau in Geometric Mechanics, that will play a fundamental role in Souriau Lie Groups Thermodynamics to define an extension of the Fisher Metric from Information Geometry: "Fisher-Souriau metric".
The antisymmetric bilinear map (31) and (32), with definition (27) and (28), introduced by Souriau is exactly equal to the mathematical object extensively studied in Chapter 4 of Koszul's book: In this book, Koszul has studied this antisymmetric bilinear map considering the following developments. For any moment map µ, Koszul defines the skew symmetric bilinear form c µ (a, b) on Lie algebra by: Koszul observes that if he uses: by developing dµ(ax) = t ad a µ(x) + dθ µ (a) , x ∈ M, a ∈ g, he obtains: He has then: and the property: Koszul concludes by observing that if the moment map is transform as µ = µ + φ then we have: koszul highlights the property that: In Chapter 4, Koszul introduces the equivariance of the moment map µ. Based on the definitions of the adjoint and coadjoint representations of a Lie group or a Lie algebra, Koszul proves that when (M; ω) is a connected Hamiltonian G-space and µ : M → g * a moment of the action of G, there exists an affine action of G on g*, whose linear part is the coadjoint action, for which the moment µ is equivariant. This affine action is obtained by modifying the coadjoint action by means of a cocycle. This notion is also developed in Chapter 5 for studying Poisson manifolds.
Defining classical operation Ad s a = sas −1 , s ∈ G, a ∈ g, ad a b = [a, b] , a ∈ g, b ∈ g and Ad * s = t Ad s −1 , s ∈ G with classical properties: Koszul considers: From which, he obtains: Koszul then study µ • s M − Ad * s • µ : M → g * , and develops: d Ad * s • µ, a = d µ • s M , a and then proves that Koszul considers the cocycle given by θ µ (s) = µ(sx) − Ad * s µ(x) , s ∈ G, and observes that: From this action of the group on dual Lie algebra: Koszul introduces the following properties: This Koszul study of the moment map µ equivariance, and the existence of an affine action of G on g*, whose linear part is the coadjoint action, for which the moment µ is equivariant, is at the cornerstone of Souriau theory of geometric mechanics and Lie groups thermodynamics.
We compare Souriau and Koszul notations in Figure 7.
Koszul considers the cocycle given by , and observes that: From this action of the group on dual Lie algebra: Koszul introduces the following properties: This Koszul study of the moment map μ equivariance, and the existence of an affine action of G on g*, whose linear part is the coadjoint action, for which the moment μ is equivariant, is at the cornerstone of Souriau theory of geometric mechanics and Lie groups thermodynamics.
We compare Souriau and Koszul notations in Figure 7. We have also to make reference to Muriel Casalis' papers [41,42] on this topic.

Souriau Model of Generalized Entropy Based on Legendre and Laplace Transforms
At the step of the development of Souriau Lie groups thermodynamics, we will introduce generalized Souriau definition of entropy. Souriau first start to define "Laplace transform": Let E a vector space of finite size, µ a measure of its dual E * , then the function given by: for all α ∈ E such that the integral is convergent. This function is called (generalized) Laplace transform. This transform F of the measure µ is differentiable inside is definition set de f (F). Its p-th derivative is given by the following convergent integral for all point inside de f (F): -F is semi-definite convex function, F(α) > 0, ∀α ∈ de f (F) (84) f = log F is convex and semi-continuous -Let α an interior point of de f (F) then: See [107], for links between dual convex functions and optimization.
Before introducing Entropy, Souriau introduced the following lemma: Lemma 1. Let X be a locally compact space, Let λ a positive measure of X, having X as support, then the following function Φ is convex: such that the integral is converging.
The integral is strictly positive when it converges, and then insures existence of its logarithm. The epigraph of Φ is the set of h y such that X e h(x)−y λ(x)dx ≤ 1. Convexity of exponential shows that this epigraph is convex. Finally, Souriau introduced the "negentropy" as Legendre transform of the function Φ:

Definition 1. [Souriau Entropy Definition]
We call "Boltzmann Law" (relative to λ) all measure µ of X such that the set of real values: This definition of entropy by Souriau is a general scheme that can be extended to highly abstract spaces preserving Legendre structure [108], if we can define generalized Laplace transform. These operations of Laplace and Legendre transforms are the core contextures of theory of Information and Heat, generating the well-defined structures, from which we can preserve the definition of "average value". Jean-Marie Souriau explained this contexture property in the following sentence: "Il est évident que l'on ne peut définir de valeurs moyennes que sur des objets appartenant à un espace vectoriel (ou affine); donc-si bourbakiste que puisse sembler cette affirmation-que l'on n'observera et ne mesurera de valeurs moyennes que sur des grandeurs appartenant à un ensemble possédant physiquement une structure affine. Il est clair que cette structure est nécessairement unique-sinon les valeurs moyennes ne seraient pas bien définies." (In English: It is obvious that one can only define average values on objects belonging to a vector (or affine) space; Therefore-so this assertion may seem Bourbakist-that we will observe and measure average values only as quantity belonging to a set having physically an affine structure. It is clear that this structure is necessarily unique-if not the average values would not be well defined.) See also papers of Kostant [109] and Leray [100] for generalized Laplace transforms.

Illustration of Souriau Thermodynamics of a Centrifuge System
Duhem [110][111][112][113] and Poincaré [114] have studied statistical mechanics model of centrifuges. We will illustrate Souriau's Lie groups thermodynamics for Souriau Gibbs states for Hamiltonian actions of subgroups of the Galilean group, as illustrated in Souriau's book [43] and more recentltly by Charles-Michel Marle [23].
Consider a Galilean Lie group: Galilean Lie algebra: Action of Lie group: Galilean transformation on position and speed is given by: Souriau has proved that this action is Hamiltonian, with the map J, defined on the evolution space of the particle, with value in the dual g* of the Lie algebra G, as momentum map: where the coupling formula is given by: Entropy 2018, 20, 840 25 of 60 with: Souriau gave the demonstration for the Galilean moment map for a free particle, considering the definition of moment map: σ(dp)(δp) = −d J, Z , ∀dp (97) and the definition of tangent vector field: Then, as General Lagrange 2 form for a force F is: If F is equal to zero, we obtain: and the co-cycle is given by: The main Souriau idea was to define the Gibbs states for one-parameter subgroups of the Galilean group. Souriau has proved that action of the full Galilean group on the space of motions of an isolated mechanical system is not related to any equilibrium Gibbs state (the open subset of the Lie algebra, associated to this Gibbs state, is empty). Then, if we consider the 1-parameter subgroup of the Galilean group generated by b element of Lie algebra, is the set of matrices: and: Then, Gibbs state defined for a gas enclosed in a moving box could be computed by Souriau formula. If we fix the affine Euclidean reference frame 0, speed → δ /ε and an acceleration → α /ε. If we consider a gas made of N point particles, indexed by i ∈ {1,2, . . . , N}, enclosed in a box with rigid and undeformable walls, whose motion is described by the action of the 1-parameter subgroup of the Galilean group, A(t/ε) where t ∈ R. If we consider m i , r i (t), v i (t), respectively the mass, position vector and velocity vector of the ith particle at time t. If we assume free particle and we neglect contributions given by the collisions of the particles between themselves collisions with the walls, then we can write: The important idea is to observe that J i , β is invariant by the action of 1-parameter subgroup. The proof of J i , β invariance is based on Souriau equation for default of equivariance with cocyle. If the action of the 1-parameter subgroup is exp t ε β , according to Souriau equation: We obtain for: that can be reduded by using the properties: and: To obtain Souriau's Gibbs maximum entropy density, we have to use the following change of variables: We can then write: and finally, the Souriau Gibbs density is given by: If we consider the case of the centrifuge (as for a butter churn, device used to convert cream into butter), the parameter of Galilean group Lie algebra are reduced to: with variables: We obtain the closed form for maximum entropy Souriau-Gibbs density: This equation describes the behaviour of a gas made of point particles of various masses in a centrifuge rotating at a constant angular velocity and explains the observation that the heavier particles concentrate farther from the rotation axis than the lighter ones. Souriau made reference to thermodynamics of butter churn (see Figure 8). Souriau Lie groups thermodynamics provides right results if we apply it to subgroups of Galileo group, as previous example of a cylindrical box with fluid with an invariance sub-group of size 2 (rotation along the axis, time translation) providing a 2-dimensional Souriau (Planck) temperature-vector. Souriau has observed that the process, by which a refrigerated centrifuge transmits its own temperature-vector to its content, has two names: thermal conduction and viscosity, depending on the temperature-vector component that is considered. Conduction and viscosity should therefore be unified in a fundamental theory of irreversible processes (theory that remains to be constructed).
In the Appendix, we develop a solution given by Roger Balian [25] for the previous case of centrifuge thermodynamics based on classical methods. Balian recover the same Gibbs density but by introducing an additional Lagrange hyper-parameter associated to total angular momentum. Balian has computed the Boltzmann-Gibbs distribution without knowing the Souriau equations (exercice 7b of). Balian started by considering the constants of motion that are the energy and the component z J of the total angular momentum ( ) Balian observed that he must add Souriau Lie groups thermodynamics provides right results if we apply it to subgroups of Galileo group, as previous example of a cylindrical box with fluid with an invariance sub-group of size 2 (rotation along the axis, time translation) providing a 2-dimensional Souriau (Planck) temperature-vector. Souriau has observed that the process, by which a refrigerated centrifuge transmits its own temperature-vector to its content, has two names: thermal conduction and viscosity, depending on the temperature-vector component that is considered. Conduction and viscosity should therefore be unified in a fundamental theory of irreversible processes (theory that remains to be constructed).
In the Appendix, we develop a solution given by Roger Balian [25] for the previous case of centrifuge thermodynamics based on classical methods. Balian recover the same Gibbs density but by introducing an additional Lagrange hyper-parameter associated to total angular momentum.
Balian has computed the Boltzmann-Gibbs distribution without knowing the Souriau equations (exercice 7b of). Balian started by considering the constants of motion that are the energy and the component J z of the total angular momentum J = ∑ i (r i × p i ). Balian observed that he must add to the Lagrangian parameter, given by (Planck) temperature β for energy, an additional one associated with J z . He identifies this additional multiplier with −βω by evaluating the mean velocity at each point. He then introduced the same results also by changing the frame of reference, the Lagrangian and the Hamiltonian in the rotating frame and by writing down the canonical equilibrium in that frame. He uses the resulting distribution to find, through integration, over the momenta, an expression for the particles density as the function of the distance from the cylinder axis. The main Souriau model advantage is that we can define covariant Gibbs density for dynamical systems, only by applying formulas without any considerations [64].

Higher-Order Model of Lie Groups Thermodynamics Based on Poly-Symplectic Vector Valued Model
As observed by Souriau in Chapter IV of [43], the Gausian density is a maximum entropy density of 1st order. Considering multivariate Gaussian density, this remark is clear if we replace classical parameterization z and (m, R) by the new parameterization, linked to information geometry coordinates, ξ and β: We can observe in previous equations that classical multivariate Gaussian density, classically expressed by p (m,R) (z) = To understand the meaning of these tensors, we can consider them as homeomorph to the following respective matrices ξ = zz T z 0 1×n 0 ,ξ = R + mm T m 0 1×n 0 and β = 1 2 R −1 −R −1 m 0 1×n 0 with β, ξ = Tr βξ T (see [91] for more details). Z is the classical normalization constant that is equal to log(Z) = n 2 log(2π) + 1 2 log det(R) + where we recover classical definition of entropy by Clausius dS = dQ T when β = 1 T andξ = Q heat). We can also defined Fisher metric of information geometry by can observe that classical multivariate Gaussian Density pξ (ξ) = 1 Z e − β,ξ is a maximum entropy Gibbs density of 1st order with respect to the tensorial variableξ = E[ξ] = m R + mm T . Classically Gaussian density is considered as a maximum entropy Gibbs density of 2nd order where p (m,R) (z) = 1 The question is then, could we define a Gaussian density of higher order?
We have seen that Souriau has replaced classical maximum entropy approach by replacing Lagrange parameters by only one geometric "temperature vector" as element of Lie algebra. In parallel, Ingarden has introduced second and higher order temperature of the Gibbs state that could be extended to Souriau's theory of thermodynamics. The question is then, how to extend the Souriau model to define an higher order Lie groups thermodynamics. For this purpose, we propose to consider multi-symplectic geometry and more particularly poly-symplectic geometry [115]. The variational problems generalization with several variables was developed by Volterra in two papers [116,117] where two different generalizations of the Hamilton system of equations are introduced. In parallel, De Donder [53] has also studied this approach in a geometrical framework based on Elie Cartan's idea of invariant structure with no dependence to local coordinates and based on affine multisymplectic manifold. We can also formalize the multisymplectic geometry with an extension of the Poincaré-Cartan invariant integrals. Frédéric Hélein has observed the fact that different theories could cohabitate was considered jointly by Lepage [54], Dedecker [118,119] and Kijowski [92][93][94]. The Lepage-Dedecker theory was developed by Hélein [120], and the modern formulation using the multisymplectic (n + 1)-form as the fundamental structure of the theory starts with Kijowski's papers. The geometrical multisymplectic approach uses the generalized Legendre correspondence introduced by Lepage and Dedecker and Hamiltonian formalism developed by Hélein [55]. We can also make references to poly-symplectic formulation of physical systems by Carathéodory [121] and Weyl [122].
Among all multi-symplectic models, the more natural multi-valued one that preserve the notion of (poly-)moment map has been initiated by Günther based on n-symplectic model. Günther has shown that the symplectic structure on the phase space remains true, if we replace the symplectic form by a vector valued form, that is called poly-symplectic. The Günther formalism is based on the notion of a poly-symplectic form, which is a vector valued generalization of symplectic forms. Hamiltonian formalism for multiple integral variational problems and field theory is presented in a global geometric setting. Günther has introduced in this poly-symplectic formalism: Hamiltonian equations, canonical transformations, Lagrange systems, symmetries, Field theoretic moment mappings, a classification of G-homogeneous field theoretic systems on a generalization of coadjoint orbits.
Günther has defined six conditions for a multidimensional Hamiltonian formalism: • C0: For each field system, an evolution space can be constructed, which describes the states of the system completely. • C1: The evolution space carries a geometric structure, which assigns to each function (Hamiltonian density) its Hamiltonian equations. Günther has observed that Hamiltonian field theory by Marsden is not covariant, because C3 is not verified and causes problems in relativistic theories, and by the multisymplectic approach by Tulczyjew, based on the general theory by Dedecker, does not satisfy C1 and C2.
The key idea of Günther for this generalized Hamiltonian formalism is to replace the symplectic form in classical mechanics by a vector valued, so called poly-symplectic form with the property that: • the evolution space of a classical field will appear as the dual of a jet bundle, which carries naturally a polysymplectic structure. • canonical transformations are bundle isomorphisms leaving this poly-symplectic form invariant.
The polysymplectic approach recovers all classical results also generalize the Noether theorem based on canonical transformations and preserve the existence of momentum mappings. Christian Günther's work was inspired by the symplectic formulation of classical mechanics by Souriau and by the work of Edelen [52,123] and Rund [124] on a local Hamiltonian formulation of field theory. Edelen's work is a coordinate version of the local polysymplectic approach of Günther.
Initiated by Gunther [48,49] based on n-symplectic model [50,51], it has been shown that the symplectic structure on the phase space remains true, if we replace the symplectic form by a vector valued form, that is called polysymplectic.
In Günther's poly-symplectic model, we set: P : space of field values , φ : U → P and we consider the bundle of linear maps from R n into the tangent spaces of P: The base of R n is interpreted as n-tangent υectors of M, there is the isomophy: The natural projection is given by: The cojet space Hom(R n , TP) carries a natural R n -valued: • one-form: Θ 0 (canonical one-form): • two-form: Ω 0 = −dΘ 0 closed & non-degenerate (canonical polysymplectic form) • A closed nondegenerate R n -valued two-form Ω on a manifold M is called a polysymplectic form. The pair (M, Ω) is a polysymplectic manifold. • A polysymplectic form Ω on a manifold M is called a standard form iff M has an atlas of canonical charts for Ω, i.e., charts in which locally Ω is written as the canonical evaluation form on P x Lin (P,R n ). (M, Ω) is called a standard polysymplectic manifold.
The classification of symplectic homogeneous spaces by coadjoint orbits by Souriau belong to the major achievements in Hamiltonian mechanics. Günther has extended these results to polysymplectic manifolds. Let Ad : G × LG → LG be the adjoint action. We denote by Ad n induced action on Lin(R n , LG): Ad n g : G × Lin(R n , LG) → Lin(R n , LG) The dual of Ad n is denoted by Ad

Corollary 1. [Günther Corollary]
Let the moment map J (n) : M → Lin(LG, R n ) = LG * ⊗ R n , there is a smooth map θ (n) : with the following properties: θ (n) is a 1-cocyle for all g, h ∈ G then:

Theorem 4. [Günther Theorem (Vector-Valued Extension of Souriau Theorem)]
The map: is an affine operation of G on LG * ⊗ R n , and commutes for all g ∈ G.
is a poly-symplectic extension of Souriau-Fisher Metric.
Compared to the Souriau model, heat is replaced by previous polysymplectic model: Proposition 2. The characteristic function: exists.
Proof. We extrapolate Souriau's results, who proved in [1,2] that M U ⊗k (ξ).e − β k ,U ⊗k (ξ) dω is locally normally convergent using multi-linear norm U ⊗k = Sup U E, U k and where U ⊗k = U ⊗ (k) U . . . ⊗ U is defined as a tensorial product [43]. Entropy is defined by the Legendre transform of the Souriau-Massieu characteristic function: Definition 6. The poly-entropy is given by Legendre transform of the poly-symplectic characteristic function: The Gibbs density could be then extended with respect to high order temperatures.

Definition 7.
Gibbs density is defined as the maximum entropy density of poly-Entropy:

Conclusions and Possible Extensions
We have introduced contextures of geometric theory of information and heat based on Souriau's approach, but information geometry is at the interface between different geometries. First, information geometry is at the intersection between "Riemannian geometry", "complex geometry" and "symplectic geometry". Based on seminal work of Cartan on homogeneous domains and other works [125][126][127][128], information geometry is jointly founded by (see Figure 9): • Geometry of Jean-Marie Souriau: Study of homogeneous symplectic manifolds geometry with the action of dynamical groups. Introduction of the Lie groups thermodynamics in statistical mechanics [43,44]. • Geometry of Jean-Louis Koszul: Study of homogeneous bounded domains geometry, symmetric homogeneous spaces and sharp convex cones. Introduction of an invariant 2-form [9][10][11]97,98,129]. • Geometry of Erich Kähler: Study of differential manifolds geometry equipped with a unitary structure satisfying a condition of integrability. The homogeneous Kähler case studied by André Lichnerowicz [130]. The poly-entropy is given by Legendre transform of the poly-symplectic characteristic function: ( ) ) ,..., ( , ,..., The Gibbs density could be then extended with respect to high order temperatures.

Definition 7.
Gibbs density is defined as the maximum entropy density of poly-Entropy:

Conclusions and Possible Extensions
We have introduced contextures of geometric theory of information and heat based on Souriau's approach, but information geometry is at the interface between different geometries. First, information geometry is at the intersection between "Riemannian geometry", "complex geometry" and "symplectic geometry". Based on seminal work of Cartan on homogeneous domains and other works [125][126][127][128], information geometry is jointly founded by (see Figure 9): • Geometry of Jean-Marie Souriau: Study of homogeneous symplectic manifolds geometry with the action of dynamical groups. Introduction of the Lie groups thermodynamics in statistical mechanics [43,44].

•
Geometry of Jean-Louis Koszul: Study of homogeneous bounded domains geometry, symmetric homogeneous spaces and sharp convex cones. Introduction of an invariant 2-form [9][10][11]97,98,129]. • Geometry of Erich Kähler: Study of differential manifolds geometry equipped with a unitary structure satisfying a condition of integrability. The homogeneous Kähler case studied by André Lichnerowicz [130]. We have extended Souriau's Lie groups thermodynamics by a vector-valued model based on poly-symplectic geometry, introducing higher order Souriau-Gibbs density with higher order Souriau temperatures, and elements of Lie algebra. This model preserves all contextures of Souriau's We have extended Souriau's Lie groups thermodynamics by a vector-valued model based on poly-symplectic geometry, introducing higher order Souriau-Gibbs density with higher order Souriau temperatures, and elements of Lie algebra. This model preserves all contextures of Souriau's thermodynamics with covariance of Gibbs density with respect to dynamical groups in physics. Poly-moment maps are compliant with the Noether theorem generalization in vector-valued cases.
The Jean-Marie Souriau model and equations were extensively studied in the Koszul Lecture given in China in 1986 "Introduction to Symplectic Geometry", in Chinese (see Figure 10). This book should be translated in English in 2019. Chuan Yu Ma has written on the Koszul book: "This beautiful, modern book should not be absent from any institutional library. . . . . During the past eighteen years there has been considerable growth in the research on symplectic geometry. Recent research in this field has been extensive and varied. This work has coincided with developments in the field of analytic mechanics. Many new ideas have also been derived with the help of a great variety of notions from modern algebra, differential geometry, Lie groups, functional analysis, differentiable manifolds and representation theory. [Koszul's book] emphasizes the differential-geometric and topological properties of symplectic manifolds. It gives a modern treatment of the subject that is useful for beginners as well as for experts." thermodynamics with covariance of Gibbs density with respect to dynamical groups in physics.
Poly-moment maps are compliant with the Noether theorem generalization in vector-valued cases. The Jean-Marie Souriau model and equations were extensively studied in the Koszul Lecture given in China in 1986 "Introduction to Symplectic Geometry", in Chinese (see figure 10). This book should be translated in English in 2019. Chuan Yu Ma has written on the Koszul book: "This beautiful, modern book should not be absent from any institutional library. …. During the past eighteen years there has been considerable growth in the research on symplectic geometry. Recent research in this field has been extensive and varied. This work has coincided with developments in the field of analytic mechanics. Many new ideas have also been derived with the help of a great variety of notions from modern algebra, differential geometry, Lie groups, functional analysis, differentiable manifolds and representation theory. [Koszul's book] emphasizes the differential-geometric and topological properties of symplectic manifolds. It gives a modern treatment of the subject that is useful for beginners as well as for experts." We have seen that in geometrical mechanics, the Galileo group related to classical mechanics: and its central extension given by the Bargman group: and Poincaré group in relativity. We then observe, that affine group or its sub-groups are at cornerstone of different disciplines such as:

•
In robotics, the special Euclidean group SE(3) which is the homogeneous Galileo group (robotics also consider the group of similitudes SIM(3)): In information geometry, the general affine group is involved A(n,R) for exponential family: We have seen that in geometrical mechanics, the Galileo group related to classical mechanics: and its central extension given by the Bargman group: and Poincaré group in relativity. We then observe, that affine group or its sub-groups are at cornerstone of different disciplines such as: • In robotics, the special Euclidean group SE(3) which is the homogeneous Galileo group (robotics also consider the group of similitudes SIM(3)): • In information geometry, the general affine group is involved A(n,R) for exponential family: with particular case of Gaussian density, associated by Cholesky factorisation of covariance matrix, where covariance matrix square root is triangular matrix with positive elements on its diagonal (it is a group): • In the study of homogeneous bounded domains, as the simplest one given by Poincaré upper-half plane: As illustrated in Figure 11, Jean-Marie Souriau developed these models at Carthage in Tunisia and at Marseilles in France during 50's and 60's. Jean-Marie Souriau was motivated by group invariance, not only in physics but also in neuroscience. Souriau intuition was highly premonitory, because this neuroscience domain has been developed few decades after by Alain Berthoz at College de France (http://public.weconext.eu/academie-sciences/2017-10-03_5a7/video_id_002/index.html) and by Daniel Bennequin (https://www.youtube.com/watch?v=a-ctwxBpJxE) to study the brain sense of movment. We can read in Souriau's text the very interesting remarks on geometry and neuroscience: "Je me suis dit, à force de rencontrer des groupes, il y a quelque chose de caché là-dessous. La catégorie métaphysique des groupes qui plane dans l'empyrée des mathématiques, que nous découvrons et que nous adorons, elle doit se rattacher à quelque chose de plus proche de nous. En écoutant de nombreux exposés faits par des neurophysiologistes, j'ai fini par apprendre le rôle primitif du déplacement des objets. Nous savons manipuler ces déplacements mentalement avec une très grande virtuosité. Ce qui nous permet de nous manipuler nous-même, de marcher, de courir, de sauter, de nous rattraper quand nous tombons, etc. Ce n'est pas vrai seulement pour nous, c'est vrai aussi pour les singes ; ils sont beaucoup plus adroits que nous pour anticiper les résultats d'un déplacement. Pour certaines opérations élémentaires de «lecture», ils vont même dix fois plus vite que nous. Beaucoup de neurophysiologistes pensent qu'il y a une structure spéciale génétiquement inscrite dans le cerveau, le câblage d'un groupe . . . Lorsque il y un tremblement de terre, nous assistons à la mort de l'Espace. . . . Nous vivons avec nos habitudes que nous pensons universelles. . . . La neuroscience s'occupe rarement de la géométrie . . . Pour les singes qui vivent dans les arbres, certaines propriétés du groupe d'Euclide sont mieux câblées dans leurs cerveaux." (In Engish: "I said to myself, because of meeting groups everywhere, there is something hidden there. The metaphysical category of groups that hovers in the empyrean of mathematics, which we discover and adore, must be connected with something closer to us. Listening to many presentations by neurophysiologists, I ended up learning the primitive role of moving objects. We know how to manipulate these movements mentally with great virtuosity. That allows us to manipulate ourselves, to walk, run, jump, catch up when we fall, and so on. This is not true only for us, it is true also for monkeys; they are much more adroit than we are to anticipate the results of a trip. For some basic "reading" operations, they are even ten times faster than us. Many neurophysiologists think that there is a special structure genetically inscribed in the brain, the wiring of a group . . . When there is an earthquake, we witness the death of Space. . . . We live with our habits that we think universal. . . . Neuroscience rarely deals with geometry . . . For monkeys living in trees, some of Euclid's group properties are better wired in their brains.) Our new research directions will concern extension of "Le Hasard et la Courbure (Randomness and Curvature)" (title of Yann Ollivier HDR), that we have synthetized in Souriau-Fisher metric to "Le Hasard et la Torsion (Randomness and Torsion)" based on Elie Cartan works founded on Cosserats brothers model of elasticity [125][126][127]131].
Funding: This research received no external funding.

Conflicts of Interest:
The author declares no conflict of interest.

Appendix A. Günther's Polysymplectic Model
We recall in this appendix, a synthesis of Christian Günther Poly-symplectic model with his initial notation [48]. Our new research directions will concern extension of "Le Hasard et la Courbure (Randomness and Curvature)" (title of Yann Ollivier HDR), that we have synthetized in Souriau-Fisher metric to "Le Hasard et la Torsion (Randomness and Torsion)" based on Elie Cartan works founded on Cosserats brothers model of elasticity [125][126][127]131].
Funding: This research received no external funding.

Conflicts of Interest:
The author declares no conflict of interest.

Appendix A. Günther's Polysymplectic Model
We recall in this appendix, a synthesis of Christian Günther Poly-symplectic model with his initial notation [48]. The bundle of linear maps from R n into the tangent spaces of Q.
If a base of R n is chosen, can also be interpreted as n-tangent υectors of Q, there is the isomophy: The natural projection is given by: In analogy to the canonical forms on the cotangent bundle, the cojet space Hom(R n , TQ) carries a natural R n -valued: In the natural bundle coordinates the canonical forms on Hom(R n , TQ) have the local representation: Following diffeomorphism leaves invariant one and two forms: f : Q → Q and I n * f : Hom(TQ, R n ) → Hom(TQ, R n ) (I n * f ) * Θ 0 = Θ 0 and (I n * f ) * Ω 0 = Ω 0 (A7) Definition A1. A closed nondegenerate R n -valued two-form Ω on a manifold M is called a polysymplectic form. The pair (M, Ω) is a polysymplectic manifold.
The classification of linear polysymplectic forms is not trivial, because two polysymplectic forms are not necessarily locally equivariant.
Definition A2. A polysymplectic form Ω on a manifold M is called a standard form iff M has an atlas of canonical charts for Ω, i.e., charts in which locally Ω is written as the canonical evaluation form on Q x Lin (Q,R n ). (M, Ω) is called a standard polysymplectic manifold.
The polysymplectic structure provides the procedure which assigns to a function on M, the Hamiltonian, its associated Hamiltonian equations. Let (M, Ω) a polysymplectic manifold: An affine sub bundle of Hom(R n , TQ) is defined by: Definition A3. Ω #−1 (dH) is called the system of Hamiltonian partial differential equations associated with the Hamiltonian function H. A smooth map ψ : U → Mis a solution of Ω #−1 (dH) iff: Theorem A1. Let (M, Ω) be a standard polysymplectic manifold, (p,q) canonical coordinates for Ω on M, and H a Hamiltonian function. A smooth map ψ : U → M is a solution of Ω #−1 (dH) iff in canonical coordinates: If a base e 1 , . . . , e n of R n is chosen and p(u) = (p 1 (u), . . . , p n (u)) with respect to this base, then the equations take the form: Proof.

Proposition A2 [Günther Proposition].
Let Λ : G × M → M be a strongly polysymplectic group action with momentum map µ : M → Lin(LG, R n ) = LG * ⊗ R n . Assume M is connected. Then the map: is a constant on M for all g ∈ G.
Corollary A1. There is a smooth map χ: with the following properties: • is a 1-cocyle for all g, h ∈ G then • bilinear map ϕ on LG: ϕ := L χ : LG → LG * ⊗ R n , ϕ : LG × LG → R n is a 2 cocycle

Theorem A2. [Günther Theorem (Vector-valued extension of Souriau Theorem)]
Let Λ : G × M → M be a polysymplectic action with momentum map µ : M → LG * ⊗ R n . Then the map: is an affine operation of G on LG * ⊗ R n , and commutes for all g ∈ G and µ is G-equivariant. Proof.
Ξ is an action.
Christian Günther in a never found 1987 paper wrote that "The mathematical framework developed in this paper is used in a separate publication to provide a rigorous foundation for field theory". For a more recent study of Günther's poly-symplectic model, we make reference to [133].

Appendix B. Fisher Metric for Multivariate Gaussian Density
We will in the following illustrate information geometry for multivariate Gaussian density: If we develop: We can write the density as a Gibbs density: We can then rewrite density with canonical variables: The first potential function (free energy/logarithm of characteristic function) is given by: We verify the relation between the first potential function and moment: The second potential function (Shannon entropy) is given as a Legendre transform of the first one: components of ξ will provide moments of the first and second order of the density of probability ˆ( ) p ξ ξ . He used this change of variable

Appendix C. Geometric Definition of Legendre Transform by Chasles as Reciprocal Polar with Respect to a Paraboloid
The Legendre transform plays a central role related to duality and convexity. Adrien-Marie Legendre [102] has introduced the Legendre transform to solve a minimal surface problem given by Monge (Monge requested him to consolidate its proof), with a link to Poncelet duality [103]. Chasles and Darboux interpreted the Legendre transform as reciprocal polar with respect to a paraboloid (re-used by Hadamard and Fréchet in calculus of variations). Before Legendre, Alexis Clairaut introduced a Clairaut Equation that has been developed by Maurice Fréchet to characterize «distinguished densities» (densities with parameters that have covariance matrix reaching the Fréchet-Cramer-Rao Bound) [9].
Legendre Transform transformes one fonction defined by its value in one point in a fonction defined by its tangent, as illustrated in Figure A2.

Appendix C. Geometric Definition of Legendre Transform by Chasles as Reciprocal Polar with Respect to a Paraboloid
The Legendre transform plays a central role related to duality and convexity. Adrien-Marie Legendre [102] has introduced the Legendre transform to solve a minimal surface problem given by Monge (Monge requested him to consolidate its proof), with a link to Poncelet duality [103]. Chasles and Darboux interpreted the Legendre transform as reciprocal polar with respect to a paraboloid (re-used by Hadamard and Fréchet in calculus of variations). Before Legendre, Alexis Clairaut introduced a Clairaut Equation that has been developed by Maurice Fréchet to characterize «distinguished densities» (densities with parameters that have covariance matrix reaching the Fréchet-Cramer-Rao Bound) [9].
Legendre Transform transformes one fonction defined by its value in one point in a fonction defined by its tangent, as illustrated in Figure A2. Darboux gave in his book one interpretation of Chasles: "Ce qui revient suivant une remarque de M. Chasles, à substituer à la surface sa polaire réciproque par rapport à un paraboloïde». In the lecture «Leçons sur le calcul des variations", Hadamard, followed by Vessiot, used the reciprocal polar of figurative, and figuratrice. This has also been developed by Belgodère as presented by Cartan on «Extrémale d'une surface» [134,135]. Polarity on the plane is a transformation taking points to lines Darboux gave in his book one interpretation of Chasles: "Ce qui revient suivant une remarque de M. Chasles, à substituer à la surface sa polaire réciproque par rapport à un paraboloïde». In the lecture «Leçons sur le calcul des variations", Hadamard, followed by Vessiot, used the reciprocal polar of figurative, and figuratrice. This has also been developed by Belgodère as presented by Cartan on «Extrémale d'une surface» [134,135]. Polarity on the plane is a transformation taking points to lines and dually lines to points. A polarity preserves incidence and has degree 2. For a point P (that we name the pole) a conic polarity transforms it to its image which is a line p (that we name the polar) as follows: from P we draw the two tangents to the conic, which touch it in the points Q, R. If we now connect points Q, R with a line p we obtain the polar line of the pole P. A Self-conjugate point Q is incident with its polar q; that is Q lies on q.
Geometric interpretation of the Legendre transform by reciprocal polar with respect to a paraboloid is given by the following simple development. First, let's consider the surface: z = f (x, y) with p = ∂z ∂x and q = ∂z ∂y (A41) We consider the equation of the paraboloid: Reciprocal polar with respect to paraboloid has coordinates: X, Y, Z The polar plan with respect to paraboloid of this reciprocal polar Xx + Yy − z − Z = 0 should be equal to tangent plan of the surface at point (x 0 , y 0 , z 0 ): This equality provides: This is the Legendre transform. So in classical thermodynamics, the Legendre transform S(Q) = β, Q − Φ(β) is linked with polar reciprocal with respect to the paraboloid: We can develop other properties of Legengre transform. Let's z = f (x, y) with p = ∂z ∂x and q = ∂z ∂y and X = p , Y = q , Z = px + qy − z the Legendre transform. We compute the first derivative of Z: dZ = PdX + QdY with P = ∂Z ∂X and Q = ∂Z ∂Y (A46) Z = px + qy − z ⇒ dZ = pdx + qdy − dz + xdp + ydq ⇒ dz = pdx + qdy X = p, Y = q dZ = xdX + ydQ ⇒ P = x, Q = y (A47) We compute the 2nd derivative of Z: The link with with contact transformations is then the following. Considering new variables X,Y,Z and P,Q the derivatives of Z with respect to X and Y, problem of finding in which case this five quantities could be express of x,y,z,p and q est the same problem where we look for five functions X,Y,Z,P and Q of five independant variables x,y,z,p and q satisfying the differential equation: where ρ is a function of x,y,z,p and q.

Proof.
and the reciprocal: Links with Ampere transformation is given then by the following developments. Let's consider Ampere transformation: dz − pdx − qdy = d(z − qy) − pdx + ydq Set Z = z − qy, X = x, Y = q P = p, Q = −y ⇒ dZ − pdX − QdY = dz − pdx − qdy (A53) Then ρ = 1, and we have a contact transformation, also valid when Legendre transform is no longer valide (when rt − s 2 = 0, p and q are not independant) The link between Legendre transformation and Ampere transformation is then deduced. Legendre transform is obtained by same equality: We can set: X = p, Y = q, Z = z − px − qy P = x, Q = y (A55) For complementary studies on the Legendre transform, we can make reference to [99,101]. and its n-th derivatices given by: J . . . ⊗ J is defined as a tensorial product.

Appendix F. Relativistic Souriau Thermodynamics of Continua
We will summarize in this appendix the Souriau relativistic thermodynamics of fluids. This Souriau model about relativistic thermodynamics of continua will give a solution to Duhem's general thermodynanics: Nous avons fait de la dynamique un cas particulier de la thermodynamique, une Science qui embrasse dans des principes communs tous les changements d'état des corps, aussi bien les changements de lieu que les changements de qualités physiques ". (In English: We made dynamics a special case of thermodynamics, and science that embraces common principles in all changes of state bodies, changes of places as well as changes in physical qualities".) The objective is not to make a survey of all literature on this topic. We give this model to compare Souriau's approaches related to invariance and symmetries in thermodynamics. I think that this is the first time that Souriau relativistic model is presented in English. My objective is to underline that Souriau has replaced the geometric temperature of "Lie groups thermodynamics", where the temperature is an element of Lie algebra, by a temperature that is defined as a killing vector. I also underline, that in both models, Souriau was motivated to search solutions where the "Legendre transform" structure is preserved between Massieu thermodynamics potentials.
Kinematics is defined by the vector field Θ and the measurement of number of molecules: using two state functions, Souriau has built a (thermo-)dynamic according to the two principles: conservation of the Noetherian quantities attached to the Poincaré group, positive Entropy production. Such a dissipative fluid has movements in which the entropy production is nil; Θ is then a killing vector; the equations of motion fully integrate; Souriau found in particular the results of kinetic theory at equilibrium. This method can be used to study perfect fluids; Souriau recover the classic Lichnerowicz results; moreover, we can build, even in the non-isentropic case, an space-time 2-Form Ω which is Integral invariant (in the sense of Cartan-Poincaré) of the temperature vector Θ; this provides a generalization of Helmholtz's theorem. In weakly dissipative movements, naturally occur the two viscosity coefficients, as well as the thermal conductivity coefficient; they are accompanied by two other coefficients that may be measurable on actual fluids.
Jean-Marie Sourias has first considered the kinematics of a relativistic simple fluid, considering the following space-time vectors field by temperature vector X → Θ with: Θ generates a group with a parameter of diffeomorphisms of space-time E 4 ; the group's orbits (the current lines of the fluid) form an abstract space V 3 (has a manifold structure of dimension 3, characterized by the fact that the following projection is a restricted submersion: Let the metric tensor g Lie derivative (for the vector field X ∈ E 4 → Θ ): The Killing formula gives the symmetric tensor: Let consider positive density n of quotient manifold V 3 : Integral of n on V 3 gives the number of molecules. Its reciprocal image by projection is defined by: Particules conservation is given by: Direction of U or Θ defines a foliation of space-time E 4 . Leaves are current lines solutions of: We illustrate in Figure A3, the Souriau's midel of Thermodynamics of continua. We illustrate in Figure A3, the Souriau's midel of Thermodynamics of continua. The energy-momentum density tensor T λμ has been built by Souriau using the kinematic quantities, such as to verify the second principle.  Thermodynamic 1st principle in this model is given by: The energy-momentum density tensor T λµ has been built by Souriau using the kinematic quantities, such as to verify the second principle.